We study Riemann-Lebesgue integrability for interval-valued multifunctions relative to an interval-valued set multifunction. Some classic properties of the RL integral, such as monotonicity, order continuity, bounded variation, convergence are obtained. An application of interval-valued multifunctions to image processing is given for the purpose of illustration; an example is given in case of fractal image coding for image compression, and for edge detection algorithm. In these contexts, the image modelization as an interval valued multifunction is crucial since allows to take into account the presence of quantization errors (such as the so-called round-off error) in the discretization process of a real world analogue visual signal into a digital discrete one.