Albert A. Michelson succeeded in determining the diameter of Betelgeuse (α Orionis) by using the Mount Wilson Observatory’s 100-inch telescope with the help of the movable outer mirrors in 1921 . The interference fringe formed in the focal plane of the telescope moved by air turbulence. But the moving fringe was visible and its contrast was perceptible. From the separation of the outer mirrors where the fringe contrast disappeared he determined the diameter of the fixed star. It is important to note that the images were not reconstructed due to the air turbulence. Here I will explain the principle of image reconstruction by interferometers and the problems in optical interferometers in the following and then its extensions to 3D and spectral imaging . I also introduce you the recent result of the nonlinear-optical spectral- imaging technique of the stimulated Raman scattering (SRS) microscopy [3-5].
The success of radio interferometers and problems in optical interferometric spectral image
The principle of interferometeric image reconstruction is based on the Van Cittert-Zernike (VCZ) Theorem . If you measure the spatial coherence function over the lateral plane with regard to the line of sight, the image intensity distribution is given by its 2D Fourier transform. Since the coherence function is an analytic function and the measurement of its imaginary part or the phase angle is crucial for the direct image reconstruction. Unfortunately, the air turbulence disturbs the phase of optical coherence function as is implied by the motion of optical fringe disturbed by the air turbulence. On the other hand, the radio interferometers are relatively free from the air turbulence and images are readily reconstructed. In the radio interferometers the radio signals are directrly detected by the antennas and the electric signals can be cabled to the desired locations and the coherence functions can be obtained by using the electronic correlators . High-resolution multispectral imaging can also be carried out by the radio interferometers, because if they are equipped with the electronic time-delay circuit, they can output the coherence function with the time delay from which we can retrieve the cross-spectral density function by the Fourier transform relationship. Since the cross-spectral density function is the spectral component of the coherence function, the high-resolution multispectral imaging is straightforward. The temporal correlation function or the coherence function for each base line is usually measured first by an analogue or digital correlator depending on the generations and the cross-power or cross-correlation functions is calculated by the Fourier transform relationship. Y. Chikada et al. constructed a modern FFT spectrum analyzer with the highly parallel pipeline architecture where the digitized radio signals are numerically Fourier transformed before the correlation process . The practical problems in optical interferometry are the weakness in air turbulence and the difficulty in high-resolution multi-spectral imaging.
THE FOURIER-TRANSFORM SPECTRAL IMAGING TECHNIQUE
If you properly combine the Michelson stellar interferometer and Michelson-Morley interferometer, you can obtain the coherence function with the long time delay and reconstruct a high-resolution spatio-spectral image, which is a two-dimensional (2D) array of high-resolution spectra at each location over the object surface. The principle is based on a Fourier-transform relationship derived from the basic propagation formula describing the field correlation function in terms of the source correlation function  and can be considered as a unification of the VCZ and Wiener-Khinchine theorems. The theory was published in 1986 with a simple experiment . Shortly after, Mariotti and Ridgway suggested a simple and quite effective combination of these interferometers for use in astronomy . The Long- baseline ground-based type of this spatio-spectral interferometers using optical single-mode fibers [12, 13] and a similar type of ground-based system without the use of single-mode fibers [14, 15] are under study. In the following, the author explains the principle of the Fourier-transform spectral imaging technique simply if at all possible and then introduce briefly the current status of the Fourier-transform spectral imaging techniques in space for astronomy [16-20].
The spectral imaging technique to be discussed here is, as already mentioned, based on the unified theory combining the principles of the Michelson’s two well-known interferometers. The unified theory states that if you can measure the 3D spatial coherence function of optical disturbances coming from the far zone you can retrieve at one sweep both of the spatial and spectral information of the light source. The principle will be explained in an intuitive and comprehensible way using a simple light source model in the following. Let us consider the 3D spatial coherence of optical field illuminated by a simple object in the far zone. Here we assume two point sources as shown in Fig. 1. The two point sources in the z = 0 plane emit mutually incoherent quasi-monochromatic lights. Here, “incoherent” means the lights from the two point sources are statistically independent and “quasi-monochromatic” means that the spectral width of the light is very small compared with the average frequency. For simplicity we assume that the observation area is located sufficiently far from the light sources and the optical waves coming from the two light sources can be considered to have flat or planar wavefronts. We also assume to focus on a certain component of polarization and deal with the optical field as scalar.
Let the two point sources be numbered (j = 1, 2) and the mean 3D wavenumbers and mean angular frequencies be and , respectively. Then, the optical wave field at the location R in the observation area at time t is given by the summation of two plane waves with a slowly-varying envelope functions that changes very slowly compared with the mean frequency or mean wavelength . Then the coherence function is given by 
where * denotes the complex conjugate, ⟨…⟩ means the statistical ensemble average, and the temporal stationarity of the optical field is assumed and the ensemble average can usually be considered to be replaceable with the time average.
From Eq. (1), we understand that the two mutually independent point sources produces in the far field the coherence function that depends only on the relative positions (R2 − R1) and produces a pair of uniform complex fringes. This is the typical coherence function produced by a couple of fixed stars in the night sky in front of our eyes. The period and direction of the complex fringe are determined by each 3D wave vector, k j (j = 1 or 2), the amplitude of the fringe being proportional to the brightness of the point source. It is now clear that if we Fourier transform Eq. (1), we get in the wave vector space the delta-like functions located at with the strength proportional to ⟨| V j | 2⟩. Every luminous objects or objects illuminated by substantially incoherent light can be considered as a collection of mutually independent light sources, the above discussion directly applies to general extended incoherent objects.
The relationship between the incoherent light source and the coherence function in the far zone has been explained by using the simple light source model. But for general discussions, the propagation of cross-spectral density is adequate. Itoh and Ohtsuka  derived a general formula resting on the work by Carter and Wolf  as
where Sω (rs⊥, 0) is the spectral intensity distribution of the incoherent light source in the 2D object plane at a point rs⊥. Equation (4) includes both of the spatial and spectral information of object and can be considered as the extension of the VCZ theorem to polychromatic objects.
They derived under appropriate conditions the following equation:
They showed that from the Fourier inversion relationship the function G(k) can be retrieved from the coherence function J (ρ) in the observation area, and named the function G(k), the spectral image. Eq. (3) is said to unify the VCZ and Wiener-Khinchine theorems.
Figure 2 shows the 2D spectral image reconstructed experimentally [10, 16]. A point object with two spectral components was produced by focusing a light beam composed of He-Cd (4416-A) and He-Ne (6328-A) laser spectra. You can see the two spectral components lie on the same radial line, which indicates the angular position of object. The authors also discussed briefly on the signal to the noise ratio for astronomical observation of stars.
Double-Fourier spatio-spectral interfereometry
The double-Fourier spatio-spectral technique was suggested in 1988 assuming aperture synthesis for the very high spatial and spectral resolution . The concept of the double Fourier interferometer is shown in Fig. 3 [11, 16]. As can be seen, the interferometer is a simple and effective combination of the two famous interferometers due to Michelson. The 3D spatial shear is introduced into the incoming optical field and spectral image is retrieved by the principle of the Fourier-transform spectral imaging.
Note that the path-difference introduced into the direction of the optical axis truly introduced to create 3D shear into the field not to create the time shear into the field as is expected in the simple Fourier spectrometry. This is the district difference between the Fourier-transform spectral imaging and the spatio-spectral radio interferometry or long-baseline spatio-spectral optical interferometry using the single-mode optical fibers. The latter ones generate the cross-spectral density by introducing the time delay into the incoming optical fields at two positions with the diffraction limited spatial resolution of the two telescopes or the siderostats.
Resolution limit and field of view
The resolution limit of the Fourier transform spectral imaging technique is briefly described here . Let the measurement range of ρ of J (ρ) be Δρ = (Lx, Ly, Lz). Then, the spatial lateral resolution about k is given by δkx = 1/ Lx, δky = 1/ Ly and the spectral resolution is approximately given by δkz = 1 / Lz. On the other hand, the field of view of is limited by the size of detection cell for J (ρ). If the detection cell size is given by δρ = (lx, ly, lz), the field of view is given by Δkx = 1/ lx, Δky = 1/ ly and the detectable spectral range is spproximately given by Δkz = 1 / lz. The resolution limit and the field of view of the double Fourier technique are determined in the same way.
Projects for high-resolution spectral imaging in space - WIIT, BETTII, SPIRIT, and SPECS -
Several projects of space-borne optical or infrared interferometers are also in progress [17-19]. The SPIRIT  will have the high angular resolution (0.3 arcsec) and the detection wavelength ranges from 25 to 400 µm. For details readers are referred to Ref. . The authors believe that the space-borne optical interferometers are the most promising application of the Fourier-transform spectral imaging or double Fourier technique.
The SPIRIT is a space-based spatio-spectral interferometer based on the Fourier-transform spectral imaging or the double-Fourier technique . The aperture synthesis technique  provides us much more spatial resolution than the single disc telescope. Figure 4 shows the concept of the Space Infrared Interferometric Telescope (SPIRIT). The SPIRIT is to provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view in the wavelength range of 25 - 400 µm. The three primary scientific objectives of SPIRIT is as follows: to (1) learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. The SPIRIT is to use the wide-field imaging technique employed in WIIT where a detector array was substituted for the usual single detector in the focal plane. Each detector element sees each area in the universe at the given 3D baseline and detects the 3D spatial coherence of incoming light in the parallel fashion . Then, you can get multiple sets of fringe data for wider field of view simultaneously. This simple but very effective technique is to be used in the future interferometers such as SPIRIT and the Balloon Experimental Twin Telescope for Infrared Interferometry (BETTII) . Please note that the SPIRIT is still a pathfinder to the Sub-millimeter Probe of the Evolution of Cosmic Structure (SPECS)  and the SPECS is conceived as a 1-kilometer maximum baseline far-IR interferometer requiring multiple spacecrafts flying in a tethered formation.
COHERENCE-BASED 3D AND SPECTRAL IMAGING TECHNIQUES
Interferometric methods to obtain both the 3D spatial and the spectral information of a usual polychromatic incoherent object are forming the new modality of imagers. These techniques can be considered to be a new extension of the Fourier transform spectral imaging technique. Those suggested recently are briefly introduced [22-24].
Inverse or Back Propagation of Cross-spectral Density for 3D imaging
The idea of back-propagation of cross-spectral density was first suggested by Arimoto, Yoshimori and Itoh in 1999 . The principle of the back-propagation of cross-spectral density is based on the formula derived in . They employed the computationally less intensive propagation formula in the Fourier domain. They succeeded in reconstructing a 1D light source simulated by the He-Ne laser light and rotating diffuser behind a slit. They measured the real and imaginary parts of the 1D coherence function or the cross-spectral density using the wave-front type of interferometer and then back-propagated the cross-spectrum into the source plane. The cross-spectra in the intermediate planes were also visualized.
Digital Holographic 3D Imaging Spectrometry
Teeranutraont and Yoshimori suggested the digital holographic technique of 3D and spectral imaging and demonstrated the reconstruct experiment . They first measured the asymmetric 3D coherence function of the optical field generated by the polychromatic objects illuminated by a natural light and they obtained the cross-spectral density by the Fourier transform relationship. They calculated the cross-spectral densities at the intermediate and the source planes by the back-propagation technique.
Mach-Zehnder Radial-shearing Interferometer for 3D & Spectral imaging
Just recently, Naik, Pendrini, Takada and Osten demonstrated the Radial-Shearing technique for the fast 3D and spectral imaging . Similarly to the previous work, they obtained first the shift-invariant cross-spectral density from the radially-sheared 3D coherence function via the Fourier transform relationship and then back-propagated the cross-spectral density function with the radial shear. They used the Mach-Zehnder configuration and created the radial shear by using reversed order of lenses. The shear along the optical axis was introduced by the shift of the piezo-electric mirror. The phase errors were monitored by the He-Ne laser source and images of the red and green LEDs were clearly reconstructed. The advantage of their system is the efficiency of coherence detection.
STIMULATED RAMAN SPECTRAL MICROSCOPY
Recently, stimulated Raman scattering (SRS) microscopy has been reported as a technique for overcoming the nonresonant background limitation of CARS microscopy [3-5]. In SRS microscopy, we use two-color pulses, which are called pump and Stokes, respectively. One of the pulses is modulated temporally and both of the pulses are precisely overlapped in time, and collinearly and tightly focused on a sample. When the difference in frequency corresponds to the intrinsic molecular vibrational frequency, the optical energy of the pump pulse is transferred to the Stokes pulse as a result of coherent excitation of molecular vibration through SRS process. The resultant intensity modulation is transferred to the other pulse, which is detected by a lock-in amplifier. SRS offers several advantages such as quantitative contrast and accessibility to vibrational spectrum. Video rate imaging is achievable because the sensitivity of SRS microscopy at the shot noise limit is comparable to that of CARS microscopy .
Principle of SRS microscopy
The SRS microscopy offers novel capabilities of the label-free and high-speed imaging with the chemical contrast. Fast imaging is especially important for biomedical applications where the degradation due to the sample movement is crucial. In SRS microscopy, two color laser pulses at pump frequency ωAS and ωS are used, and either the pump or the Stokes beam is intensity modulated and both of the pulses are focused into the sample. Then, the SRS-induced intensity-modulation of the other beam is detected by a lock-in amplifier. Figure 5 shows the resultant field due to the SRS process. The higher frequency component of the side bands is the CARS (Coherent Anti-Stokes Raman Scattering) signal.
Figure 2 shows the experimental setup for the fast SRS microsopy. A Ti:sapphire oscillator generated 76 MHz, 797 nm pulses with a 4.8 ps duration and a Yb fiber oscillator (YbF)  generated 38 MHz trains of broadband pulses at a central wavelength of 1030 nm, a spectral width of ∼30 nm. The YbF pulses were lead to TBPF (tunable band-pass filter), which consisted of a Galvanic Mirror, relay lenses, and a diffraction grating. The transmission bandwidth of the TBPF was as narrow as ∼0.3 nm. The corresponding Gaussian pulses had a duration of 5.2 ps at the Fourier-transform limit. Then, the filtered pulses were amplified by concatenated Yb-doped fiber amplifiers (YDFAs) to an average power of ∼300 mW.
A typical of the experimental results  is presented in Fig. 7. This is the colored composite image of a mouse pancreas tissue converted into a monochrome image. The three major components obtained by the ICA (Independent Component Analysis) technique were used for the composition of the color image. The dimensionality of the spectral image data obtained by the SRS microscope were first reduced by the PCA (Principal Component Analysis) and then processed by the ICA algorithm . The left side is tumor area and the intermediate layer is the boundary. The right hand side is healthy tissue. You can readily identify Nucleus, White blood cell, Healthy nucleus and so on. At the first glance, you can discriminate between the healthy and tumor areas even in the monochrome image.
The unprecedented progress in computer technology enabled the advanced stage of image processing and manipulation of big and multidimensional data. I have presented high-resolution multi-spectral imaging techniques that can supply us the multidimensional big data for astronomy and medical or life sciences. The former one is now going to be used as the tools for astronomy in space. The latter one is to be used for microscopic analysis of biological tissues. I hope these techniques be the useful eyes of human being for detection of very big and very small substances. It is noted that the optical field concerned with the Fourier transform spectral imaging is located in the far zone of the light source and is locally shift-invariant in the detection area. But the coherence function itself is globally shift-variant and depends on R of J (ρ) = J (R, R + ρ) in Eq. (4). Thus, 3D or tomographic imaging is passible by using the Fourier transform spectral imaging technique ingeniously .
A part of this work (Sec 4.3) was conducted in collaboration with Canon Inc. and Prof. Sakamoto of Keio Univ..