A method and processing device are presented for reconstructing an absorption and/or scattering image of a region of interest inside a scattering medium. A mathematical model is provided being representative of a relation between the distribution of the intensity and phase of electromagnetic radiation components scattered from a medium and a certain attenuation factor, which is function of spatial variations of scattering and absorption coefficients of the medium. The mathematical is used for processing a map of distribution of the intensity of electromagnetic radiation components scattered from known locations within the region of interest, thereby producing a halftone pattern of the region of interest.