A data processing method used for a holographic imaging device, comprising the steps of:
- acquiring data of a plurality of object light holograms (IjOR, j=1, . . . , N), which are off-axis holograms of interference fringes between object lights (Oj, j=1, . . . , N) emitted from an object sequentially illuminated with illumination lights (Qj, j=1, . . . , N) and an off-axis spherical wave reference light (R) being off-axis with respect to the object lights (Oj, j=1, . . . , N), without passing the object lights (Oj, j=1, . . . , N) through a lens, wherein the illumination lights (Qj, j=1, . . . , N) are composed of parallel lights, including one vertical incident light, having mutually different incident directions (θj, j=1, . . . , N) to the object and the data is acquired for each of the incident directions;
- acquiring data of a reference light hologram (ILR), which is an off-axis hologram of interference fringes between an in-line spherical wave reference light (L) being in-line with the object lights (Oj, j=1, . . . , N) and the off-axis spherical wave reference light (R);
- generating object light wave holograms (gj(x,y), j=1, . . . , N), which represent light waves of the object lights (Oj, j=1, . . . , N) at a hologram plane (z=0) by using the data of the reference light hologram (ILR) and the respective object light holograms (IjOR, j=1, . . . , N);
- generating spatial frequency spectra (Gj(u,v), j=1, . . . , N) at the hologram plane (z=0) by Fourier-transforming each of the object light wave holograms (gj(x,y), j=1, . . . , N) and generating object light spatial frequency spectra (Hj(u,v), j=1, . . . , N)) at an object position (z=zm) by propagating each of the spatial frequency spectra (Gj(u,v), j=1, . . . , N) at the hologram plane (z=0) to the object position (z=zm)
- generating a synthetic object light spatial frequency spectrum (HT(u,v)), at the object position (z=zm), based on a fact that the illumination lights used to acquire the data of the object lights are parallel lights, enlarged so as to occupy a wider frequency space, by moving and arranging each of the object light spatial frequency spectra (Hj(u,v), j=1, . . . , N) at the object position (z=zm) in a two-dimensional space of a spatial frequency space (u,v) using a shift amount (ujs, vjs) derived based on calculation of a cross correlation function so that each of those object light spatial frequency spectra overlaps with another over an area in which changes of amplitude and phase are common to the mutually overlapped spectra, and by making the object light spatial frequency spectra (Hj(u,v), j=1, . . . , N) match mutually in the overlap area using fitting coefficients (aαβ, α≠β, α, β=1, . . . , N) obtained for adjusting mutual amplitude and phase of the object light spatial frequency spectra (Hj(u,v), j=1, . . . , N) having the overlap area mutually;
- and
- generating a synthetic object light wave hologram (hT(x,y)) used for reconstruction of the image of the object, by inverse-Fourier transforming the synthetic object light spatial frequency spectrum (HT(u,v)).