Ординатура / Офтальмология / Английские материалы / Handbook of Optical Coherence Tomography_Bouma, Tearney_2002
.pdf
Spectral Radar |
341 |
Figure 2 Fiber-optic implementation of spectral radar.
regions of the object [11]. The waves scattered back from different depths of the object are coaxially observed.
Imaging with a finite aperture causes subjective speckle at the entrance of the fiber. The speckle contrast is the actual carrier of the information. With an accepting aperture NA ¼ 0:1 of the fiber core, there is only one speckle on the core. The confocal imaging of the backscattered photons onto the fiber core also has the advantage of spatially separating the photons. For photons that have been scattered several times, there exists hardly any correlation between run time and depth information [22]. Now just these photons will no longer hit the fiber.
The backscattered wave is superimposed with the reference wave. At the interferometer exit we use a grating spectrometer to locally separate the different wavelengths. The resolution of the spectrometer is 0.05 nm. The spectrum is imaged onto an array of 1024 photodiodes, each 2:5 mm 25 m. Each photodiode can collect 6 107 photons to give low photon shot noise and a large dynamic range. The signal is transferred to the host computer by a 14 bit A/D converter. Then the Fourier transformation is performed (which actually limits the dynamic range).
12.5.2 Signal Formation
The measuring principle is based on spectral interferometry. The signal from the object consists of many elementary waves emanating from different depths z. We neglect the dispersion in the object. The scattering amplitude of the elementary waves versus depth is aðzÞ. The object signal is superimposed with the plane reference wave aR. At the exit of the interferometer a spectrometer locally separates the different wavenumbers kð¼ 2=Þ. The interference signal IðkÞ is
342 |
|
|
Lindner et al. |
|
1 |
aðzÞ exp i2k½r þ nðzÞz& dz |
2 |
IðkÞ ¼ SðkÞ aRei2kr þ ð0 |
ð3Þ |
||
where |
|
|
|
2r |
= pathlength in the reference arm (Because we care only about path |
||
|
differences we define r ¼ 0, arbitrarily.) |
|
|
2ðr þ zÞ = pathlength in the object arm |
|
||
2z |
= pathlength in the object arm, measured from the reference plane |
||
|
(Table 1) |
|
|
z0 |
= offset distance between reference plane and object surface (Table 1) |
||
n= refractive index (n ¼ 1 for z < z0 and n 1:5 for longitudinal positions in the object z > z0)
aR |
= amplitude |
of the reference |
(For |
further investigations we set |
|||||
aðzÞ |
aR ¼ 1.) |
|
|
|
|
|
|
|
|
= backscattered amplitude of the object signal. With regard to the |
|||||||||
SðkÞ |
offset z0, aðzÞ is zero for z < z0 (see Fig. 4) |
|
|||||||
= spectral intensity distribution of the light source |
|
||||||||
With these assumptions the interference signal IðkÞ can be written as |
|
||||||||
|
1 |
aðzÞei2knzdz |
2 |
|
|
|
|
||
IðkÞ ¼ SðkÞ 1 þ ð0 |
|
|
|
|
ð4Þ |
||||
or |
|
|
|
|
|
|
|
|
|
IðkÞ ¼ SðkÞ 1 þ 2 |
|
1 |
|
|
1 |
1 |
|
|
|
ð0 |
aðzÞ cosð2knzÞ dz þ |
ð0 |
ð0 |
aðzÞaðz0Þe i2knðz z 0Þ dz dz 0 |
|||||
ð5Þ
It can be seen that IðkÞ is a sum of three terms. The first term is a constant offset. The second term encodes the depth information of the object; it is a sum of cosine functions, where the amplitude of each cosine is proportional to the scattering amplitude aðzÞ. The depth z of the scattering event is encoded in the frequency 2nz of the cosine function. This term describes the well-known Mu¨ller fringes in spectral interferometry [24]. It will be shown that aðzÞ can be acquired via a Fourier transformation of the interferogram [25]. The third(autocorrelation) term describes the mutual interference of all elementary waves.
We can get aðzÞ by Fourier transformation of IðkÞ under the assumption that aðzÞ is symmetrical with respect to z. Fortunately aðzÞ ¼ 0 for all z < z0. So we can replace aðzÞ by the symmetrical expansion a^ðzÞ ¼ aðzÞ þ að zÞ. After the Fourier transformation we have to restrict ourselves to z > z0, which gives us the depth information about the object
1 |
|
|
1 |
1 |
1 |
|
IðkÞ ¼ SðkÞ 1þð 1a^ðzÞ cosð2knzÞ dz þ |
|
ð 1 |
ð 1a^ðzÞa^ðz0Þe i2knðz z0Þ dz dz 0 |
ð6Þ |
||
4 |
||||||
þ1 |
1 |
|
þ1 |
|
|
|
IðkÞ ¼ SðkÞ 1 þ ð 1 a^ðzÞe i2knzdz þ |
|
ð 1 AC½a^ðzÞ&e i2knz dz |
ð7Þ |
|||
4 |
||||||
In this notation AC½a^ðzÞ& is the autocorrelation.
Spectral Radar |
|
|
|
|
|
|
|
|
|
|
|
|
343 |
|
IðkÞ ¼ SðkÞ 1 |
|
1 |
|
1 |
|
|
|
|
||||||
þ |
FOUz |
|
a^ðzÞ þ |
FOUz AC½a^ðzÞ& |
|
|
ð8Þ |
|||||||
|
|
|
|
|
||||||||||
2 |
|
8 |
|
|
||||||||||
Performing the inverse Fourier transformation we get |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
ð Þ |
|
FOU 1 IðkÞ |
¼ FOU 1 |
|
SðkÞ ½ ðzÞ& þ |
2 |
a^ðzÞ þ |
8 |
AC½a^ðzÞ& |
9 |
||||||
¼ A ðB þ C þ DÞ
where indicates convolution. From this result the symmetrized scattering amplitude a^ðzÞ and therefore aðzÞ can be deduced. In other words, we can see the strength of the scattering versus the depth. The main feature of all FDOCT sensors is that the total distribution of the scattering amplitude aðzÞ along one A-scan is measured at once. Light that is scattered back from each scatterer in the volume contributes to the interference signal during the whole measuring time.
Figure 3 demonstrates the signal evaluation for a mirror like object. Spectral radar measures the scattering amplitude in the Fourier domain [Eq. (5)], and the evaluation by Fourier transformation delivers the scattering amplitude in the spatial
Figure 3 Signal evaluation by spectral radar. (Left) Spectral radar measures the scattering amplitude (a,c for different distances z) in the Fourier domain. (Right) Evaluation by Fourier transformation delivers the scattering amplitude in the spatial domain (b,d). (a) Interference spectrum IðkÞ versus 1024 photodiodes of the detector array. Intensity is measured as number of electrons of each photodiode. (b) Fourier transformation gives the frequency of the cosine function, measured in Fourier periods FP (1 FP ¼ 7:2 m). (c,d) For a larger distance between the mirror and the reference plane, the cosine function displays a higher frequency.
344 Lindner et al.
domain [Eq. (9)]. The interference spectrum IðkÞ for a mirror at a distance z1 from the reference plane is shown in Fig. 3a. The signal contains the cosine function a cosð2kz1Þ, which is multiplied with SðkÞ. The Fourier transform gives the frequency of the cosine function (Fig. 3b). In the spatial domain the frequency is measured in Fourier periods ðFPÞ: 1 FP ¼ 7:2 m. For a larger distance z2 between the mirror and the reference plane the cosine function has a higher frequency (Figs 3c, 3d). The finite full width at half-maximum (FWHM) of the frequency peak of about 2 FP is caused by the convolution of the scattering amplitude with the correlogram of the light source [Eq. (9)].
Equation (9) contains the information about the scattering amplitude aðzÞ. However, besides the signal term C there are three further terms A, B, and D. These terms are well known from holography. As in holography we can get rid of them by an offset z0 of the reference plane with respect to the object surface.
In the first term, A B, we get the Fourier transform of the source spectrum (‘‘correlogram’’ ¼ A) located around z ¼ 0. To separate the correlogram A from the signal C we place the reference 200 m in front of the surface of the object. A distance of 200 m is sufficient for the separation, because the correlogram is much shorter (coherence length 35 m) if SðkÞ is a smooth function without ripples.
There is one more disturbing term: A D denotes the autocorrelation terms, which describe the mutual interference of all scattered elementary waves. In the strongly scattering skin the influence of D can be neglected because the autocorrelation term is much weaker than the signal term, which is weighted by the strong reference amplitude. Moreover, these terms are also located around z ¼ 0. Therefore, the center of the autocorrelation term is separated from the object signal aðzÞ even for a small offset z0 (Fig. 4). As mentioned earlier, the outer lobes of the AC terms are too weak to disturb us. If the object shows high backscattering from large depths, there is still the possibility to perform at each position a second measurement with no reference signal and subtract that signal from IðkÞ.
Finally, we have a convolution of the signal C with the correlogram of the light source. To achieve high resolution measurements the spectral characteristic of the light source has to be taken into account. Only if the light source has a broad and smooth spectrum without noise or ripple will the convolution peaks be sufficiently narrow. Otherwise we can overcome the influence of the convolution by dividing Eq.
(9) by the correlogram of the light source.
Figure 4 Sketch of the scattering amplitude aðzÞ and the autocorrelation terms. aðzÞ is zero up to the object surface that is located at z0. Due to the reference offset z0, the object reconstruction is separated from the AC terms, which are located around z ¼ 0.
Spectral Radar |
345 |
12.5.3 Measuring Range and Measuring Uncertainty
The measuring range z of spectral radar is basically limited by the resolution of the spectrometer [17]. A large difference between the object and reference optical paths will cause a high frequency in the spectrum. According to the sampling theorem the sample frequency of the photodiode array has to be twice as large as the highest occurring frequency in the spectrum. For z ¼ zMAX, the period of the cosine fringes is k ¼ =nzMAX. Therefore the spectrometer has to resolve at least k=2. With j kj ¼ 2=2, we get
z ¼ |
2 |
ð10Þ |
41n ! |
With the parameters in our setup the measuring range is z ¼ 2:4 mm (n ¼ 1:5). However, if we measure a mirror at variable distance z, we observe a decay of
the signal peak with increasing distance. This decay is caused by the finite width of the photodiodes (1=) of the spectrometer, although the reflected amplitude stays the same. To get the real scattering amplitude aðzÞ we must first compensate for this decay by normalization. Each measurement is normalized by dividing the measured signal by the measured decay.
The achievable spatial resolution depends on the coherence length of the light source and the scattering characteristics of the object. According to our experiments on the behavior of the spatial impulse response during signal propagation in volume scatterers with OCT methods [22], the minimum resolvable distance decreases with decreasing coherence length. However, it is finally limited by the scattering characteristics of the object. To determine the spatial resolution of spectral radar we measured a well-defined multilayer object with scattering coefficients adapted to human skin [26]. It turns out that the achievable spatial resolution within a range of 1000 m is about 10 m or less.
12.6EXPERIMENTAL RESULTS
12.6.1 Measurements on Biological Objects
Spectral radar was implemented as a fiber-optic system, and in vivo measurements were performed [4,18,19]. Figure 5 shows the head of the spectral radar including the collimation optics, the devices to adjust the reference plane and the depth of focus, and the lateral scanner.
Figure 6 displays two optograms comparing human skin at the forearm and the hand. We can distinguish the typical layer structure of human skin in both measurements. The scattering amplitudes along two A-scans (lines A and B) are shown at the bottom of the figure. Along each A-scan we see first a high peak obtained from direct reflection at the surface (position 100 m) and the strongly scattering stratum corneum (s.c.). The next layer of the skin is the boundary of the weakly scattering stratum germinativum (s.g.) followed by the strongly scattering stratum papillare (s.p.). In Fig. 6 the difference between different types of skin is evident. The structure in the forearm (e.g., thickness of the s.g.: 70–80 m) is thinner than the structure at the hand (thickness of the s.g.: 300–350 m).
346 |
Lindner et al. |
Figure 5 Spectral radar is used to check a patient’s arm.
The main aim of OCT is to investigate alterations of the skin. We investigated the influence of increased moisture content on the morphological structure of the skin by OCT. In the result the thickness of the epidermis with increased moisture content leads to an average extension of the epidermis of about 10% [4]. We can demonstrate the difference between healthy skin and a superficial spreading melanoma [4]. Further, we measured the morphological structure of a fingernail and that of the growing nail beneath the skin [19].
In Fig. 7, a malignant melanoma and Bowen’s disease, a preliminary stage of skin cancer, are compared. The melanoma causes an accumulation of melanin in the cells of the epidermis and therefore high backscattering. In addition, the layer structure of the healthy skin is destroyed. The epidermis is extended, and the backscattering caused by the melanin is very strong and homogeneous. The layer structure of the skin has vanished in the malignant melanoma as in the Bowen’s disease (Fig. 7). Differently from the melanoma, the scattering amplitude in the Bowen’s disease decays exponentially. This picture shows the possible first step of spectral radar to differentiate alterations in the skin. Many further measurements have to be done to confirm the result.
The optogram of Fig. 8 displays a tunnel in the epidermis built by a larva migrane transferred by a midge bite. The tunnel in the epidermis is characterized by an almost vanishing scattering amplitude due to tissue missing from the longitudinal position beneath 150 m. The tunnel was too deep for its end to be detected.
Our most recent measurements have been performed with a more powerful SLD (output power in the fiber P ¼ 10 mW; central wavelength ¼ 840 nm; FWHM ¼ 20 nm; coherence length lC ¼ 35 m). According to our considerations
Spectral Radar |
347 |
Figure 6 Optograms from two different parts of the human skin. Left: Skin at the forearm. Right: Skin at the hand.
about the dynamic range, the measuring range increases with increasing output power. In measurement with the former SLD (output power P ¼ 1:7 mW) we could achieve a depth of only about 600 m in the skin of the fingertip (Fig. 6). Now we are able to measure scattering amplitudes up to a depth of more than 1 mm with the new SLD. Figure 9 displays an optogram of human skin at the fingertip obtained with the new SLD. The epidermis (first two layers) has a thickness of about 250 m. The third layer with high signal due to backscattering and dark zones below is the stratum papillare. The surface of the skin follows the waves of papillary ridges in the depth. The lateral B-scan is performed over the range of 2 mm. The measuring time for the whole optogram is 21 s. (About half of the time is for scan and control of the spectrometer. The exposure time for one A-scan is 20 ms.)
The lateral boundary between normal skin and a malignant melanoma can be seen in Fig. 10. The melanoma causes an accumulation of melanin in the cells of the epidermis and therefore high backscattering (white areas at the bottom of the image). The epidermis is extended, and the typical layer structure is destroyed. At the top of the image you see healthy human skin with a three-layer structure. The result delivered in Fig. 10 corresponds well to the histology. The melanoma has an average
348 |
Lindner et al. |
Figure 7 Optograms of malignant melanoma (left) and Bowen’s disease (epithelioma) (right). Both skin alterations destroy the layer structure of the healthy skin. The epidermis is extended. The scattering amplitude in the malignant melanoma is strong and homogeneous. The scattering amplitude of the Bowen’s disease decays exponentially.
thickness of about 550 m. This depth was confirmed by further measurements. The measurement was performed with the 10 mW SLD.
12.6.2 Technical Applications
Although spectral radar was developed for the acquisition of 3-D data of biological tissue, it can be used for the examination of technological objects as well [18]. In Fig. 11 we show the structure of a multilayer printed circuit board. The board contains four layers of conducting lines, at depths of about 200, 400, 600, and 700 m. The first layer is on the surface. The last one is the copper baseplate. The conduction lines are also made of copper. Because the light cannot penetrate copper, no signal can be acquired from underneath the metal. Another application is the measurement of the thickness of paint. Three layers of paint on a car body can be detected and measured.
Spectral Radar |
349 |
Figure 8 Optogram of a tunnel in the epidermis built by a larva. Tunnel diameter 250 m:
Figure 9 Optogram of human skin at a fingertip. With the new SLD (output power P ¼ 100 mW) the A-scans shows the scattering amplitude in a depth of up to nearly 1 mm.
350 |
Lindner et al. |
Figure 10 Lateral boundary between normal skin and a malignant melanoma. The measurement was performed with the 10 mW SLD.
Figure 11 Structure of a multilayer printed circuit board. Fiber bundles and conduction lines in several layers are visible.
