A method for 3D imaging of a scattering sample, in particular for determining the spatial scattering-power distribution S(x,y,z) of the sample, having the following steps a) splitting laser light having a wavelength that can be selected within a predetermined band width, into a sample and a reference light beam, b) illumining the sample at N different wavelengths λ n, n = 1,..., N, and returning the light scattered in the sample and also the reference light onto a 2D light detector array, c) superposing the scattered sample light with the reference light on the 2D light detector array and recording in each case one hologram for each of the N wavelengths, d) calculating the electric wave field of the sample light scattered in the sample for each of the N wavelengths in the plane of the 2D light detector array, e) propagating, using a propagation core to the sample light, all calculated wave fields of the sample light, in a predetermined reconstruction plane that is aligned perpendicular to the optical axis of the sample light beam, f) calculating the transverse space-frequency spectrum E n ( k x, k y, k n), with k n = 2 π / λ n of the electric wave field scattered in the sample and propagated in the reconstruction plane for each of the N wavelengths, characterized by the further steps: g) phase matching of the transverse space-frequency spectrum for each of the N wavelengths of the wave fields of the sample light for compensation of the propagation-time differences produced by steps d) and/or step e) between the calculated wave field of the sample light and wave field of the reference light, wherein a propagation time difference introduced in step d) is compensated by multiplication of the reference wave R ( x; y; k) of the reference light with a phase factor, wherein the phase factor is selected so that the resulting reference wave R' ( x; y; k) for a coordinate pair ( x,y) no longer depends on the wave number k and/or a propagation time difference introduced in step e) is compensated by multiplication of the propagation core with a phase factor, wherein the phase factor is selected so that the propagation core is independent from the wave number k for a fixed transverse spatial frequency ( kx,ky), h) calculating the depth-dependent scattering-power spectrum S ( kx,ky,z) by one-dimensional Fourier transformation of the space-frequency spectra for each transverse space frequency ( kx,ky) along the optical axis of the sample beam, i) calculating the scattering power S ( x,y,z) in at least one layer of the sample in the surroundings of the reconstruction plane. a) splitting laser light having a wavelength that can be selected within a predetermined band width, into a sample and a reference light beam, b) illumining the sample at N different wavelengths λ n, n = 1,..., N, and returning the light scattered in the sample and also the reference light onto a 2D light detector array, c) superposing the scattered sample light with the reference light on the 2D light detector array and recording in each case one hologram for each of the N wavelengths, d) calculating the electric wave field of the sample light scattered in the sample for each of the N wavelengths in the plane of the 2D light detector array, e) propagating, using a propagation core to the sample light, all calculated wave fields of the sample light, in a predetermined reconstruction plane that is aligned perpendicular to the optical axis of the sample light beam, f) calculating the transverse space-frequency spectrum E n ( k x, k y, k n), with k n = 2 π / λ n of the electric wave field scattered in the sample and propagated in the reconstruction plane for each of the N wavelengths, g) phase matching of the transverse space-frequency spectrum for each of the N wavelengths of the wave fields of the sample light for compensation of the propagation-time differences produced by steps d) and/or step e) between the calculated wave field of the sample light and wave field of the reference light, wherein a propagation time difference introduced in step d) is compensated by multiplication of the reference wave R ( x; y; k) of the reference light with a phase factor, wherein the phase factor is selected so that the resulting reference wave R' ( x; y; k) for a coordinate pair ( x,y) no longer depends on the wave number k and/or a propagation time difference introduced in step e) is compensated by multiplication of the propagation core with a phase factor, wherein the phase factor is selected so that the propagation core is independent from the wave number k for a fixed transverse spatial frequency ( kx,ky), h) calculating the depth-dependent scattering-power spectrum S ( kx,ky,z) by one-dimensional Fourier transformation of the space-frequency spectra for each transverse space frequency ( kx,ky) along the optical axis of the sample beam, i) calculating the scattering power S ( x,y,z) in at least one layer of the sample in the surroundings of the reconstruction plane.