|
1.IntroductionThere are several optical noninvasive methods currently used in different fields of knowledge to study object displacements and deformations. For instance, automotive and aerodynamics industries, as well as medical and biological fields use, on a daily basis, optical noninvasive methods for remote inspection, design processes, measurement of mechanical parameters (such as strain/stress), and eye correction, among many other applications. Particularly, there are two such methods that have been applied extensively and have proved their effectiveness as problem solvers: electronic speckle pattern interferometry (ESPI) and digital holography (DH). For instance, the latter has been used to measure deformations on biological tissues from a pig and on a human hand,1 or for deformation analysis and shape measurement.2 A variation of the former based on the speckle pattern contrast has been used to acquire optical tomography images for the detection of inhomogeneities.3 Recently, DH and other techniques have been applied to measure depth position of inhomogeneities,4 and the detection and analysis of defects.5, 6 In what follows, the research done on 3-D pulsed digital holography (3-D PDH) for the detection and depth quantification of inhomogeneities within semisolid objects, also known as phantoms, is presented. The object is opaque and the inhomogeneity is placed well inside the material so that it is not visible and does not cause a dent on the object surface. Previous reported research work7 used an out-of-plane sensitive PDH system to detect surface microdisplacements on phantoms. The optical setup for this out-of-plane sensitivity measures the three orthogonal components of surface displacement, with the largest component being out of the plane of the object surface, e.g., the component of a flat surface. The other two components, , , are there, however small, and cannot be separated from the former, i.e., the three displacement components are mixed. The inhomogeneities consisted in human breast tumors. Sinusoidal sound waves were used to scan frequencies from to to find the gel resonant modes response. At this frequency range, two resonant modes were found, the first at and the second at . The latter frequency was chosen for the experiments, since it produced the largest surface displacement, with an input sound power of approximately , equivalent to a pressure of . The digital hologram taken for the inhomogeneity is compared with the digital hologram taken from a phantom that does not contain an inhomogeneity. The result shows as a set of nonconcentric rings on the phantom surface. However, this is an out-of-plane sensitive system that does not render precise and reliable information regarding the actual inhomogeneity depth position, thus a 3-D PDH system was devised for its quantification. 3-D PDH employs a setup where the phantom is illuminated from three different directions.8 This is necessary to have from each illumination direction an equation containing the mixed surface displacement components , , and , thus having three unknowns with the corresponding three equations needed to solve the problem. The usual pulsed digital hologram subtraction is performed separately: a digital hologram pair from each illumination direction. Each pair is formed with the holograms taken before (first laser pulse), e.g., at the peak of the vibration cycle, and after surface displacement (second pulse), after the first, giving as a result a phase map directly related with the direction and magnitude of the displacement. The second pulse is fired as needed according to the mechanical wave propagating through the gel and that reaches its surface. Moving on with the process, each phase map is unwrapped to get the surface displacement data, vector in 3-D, via the following equation, where is the unitary sensitivity vector,with and being unitary vectors along the direction of object illumination and observation, respectively. Consider an object point located at the origin of an orthogonal coordinate system that serves to measure in the experimental set up the location coordinates for the charge-coupled device (CCD) camera and the three sources of illumination (see Fig. 1 ). Each phase map for the corresponding object illumination direction is represented by , in such a way that they are combined through Eq. 1, leaving the only unknown, which is the object vector displacement , in all three directions , , and . From these data it is now possible to individually obtain displacement information along each individual axis and thus measure the depth of the inhomogeneity within the phantom.The characteristic unwrapped phase map for a phantom with and without an inhomogeneity is such that the concentric rings obtained for the latter are destroyed for the former, and this may be immediately noticed by simple observation of the phase map. The quantitative evaluation of such unwrapped phase maps through Eq. 1 gives the actual depth position of the inhomogeneity, a unique feature of 3-D PDH. The observed surface displacements are all in the region of micrometers. 2.Experimental MethodThree hologram pairs were obtained sequentially for five locations of the inhomogeneity: 10, 14, 16, 18 and below the gel surface and to the outer edge of the inhomogeneity. The inhomogeneities used were tumor tissues of approximately in diameter. The optical arrangement is shown in Fig. 2 , where a high resolution CCD is used, with at 12 bits. A beamsplitter divides the pulsed laser beam, at , 15-ns pulse width, and of coherence length, in two: a reference and object beam. By using a rotating mirror (RM), the latter is directed sequentially, i.e., for each hologram pair, to mirrors , , and , to illuminate the object from three different positions. The reference beam is conveyed to the CCD sensor via a single-mode optical fiber. Care was taken to match the optical path lengths by pairs to be within the laser coherence length. This arrangement is now sensitive for each object illumination direction to all , , and axes. An electronic circuit was designed to synchronize and control the sinusoidal sound waves, the firing of the laser pulse, and the CCD camera trigger in such a way that the first laser pulse was fired at the peak of the sound vibration cycle, i.e., the camera shutter opens a few microseconds before this, and closes a few microseconds after the laser pulse was fired, acquiring the first digital hologram. A second laser pulse comes after the first one, so the second digital hologram is obtained. As is commonly done in digital holography, these two holograms are subtracted after the usual Fourier algorithm routine is employed, obtaining a phase map for an illumination direction. This procedure is repeated for the other two illumination directions. Thus, three phase maps are obtained containing the direction and magnitude of the surface displacement in 3-D. The phantom is a semisphere with an 8.4 cm diameter and 4 cm height, subject to sinusoidal sound waves at by means of an off-the-shelf speaker. The phantom center is located at the origin of the Cartesian coordinate system. The axis data are used to evaluate the inhomogeneity depth position. To validate and make sure that the results were repeatable, several holograms were taken for each of the phantoms at the positions given before. Since with 3-D PDH it is possible to separate individually displacement components along each axis, it is then possible to calculate normal and tangential vectors to the surface. With this in mind, it is now possible to probe that particular line that passes through the largest surface displacement. This procedure is repeated with all the available data previously taken for each object and its depth position, calibrating the 3-D PDH system in this manner. 3.ResultsSeveral experiments were performed to verify the repeatability of the technique. Figure 3 shows three unwrapped phase maps for each illumination direction, chosen at random from many taken. It is seen that the phase pattern changes according to the illumination direction, and is clearly nonsymmetric, taken for a tumor inside the phantom. These data are used in Eq. 1 to evaluate the object surface displacement along the , , and axes. Figure 4a shows the result in 3-D, for all three displacement components, for a 1 cm-diam tumor that gives a total phantom surface displacement of . Figure 4b corresponds to a phantom without any inhomogeneity and any sound present: the very small surface displacement was due to mechanical disturbances on the optics table, a feature that shows how the surface deforms under noncontrolled arbitrary conditions. The images in this figure may be compared, and by simple observation, or by looking at the color scale at its bottom, it would be possible to decide which sample has an inhomogeneity in it and which one does not. Data, converted to a pixel size, along the axis were used to find the depth location of the inhomogeneity. Figure 5 shows images chosen at random for an inhomogeneity placed at the same depth. Several lines were scanned around the area where the inhomogeneity is located on the image, and the one corresponding to the largest surface displacement is chosen, shown in white. Vertical scans were also performed to corroborate this line data and position. The procedure was repeated for all acquired image data and for an inhomogeneity 1.2 cm diameter located below the surface. It is pointed out here that to calibrate the system, all of the previous was done for all cited locations. Figure 6 shows the depth location, along the axis, of the inhomogeneity with respect to the surface. The maximum surface displacement for this tumor size and position is about . Blue, pink, and brown lines correspond to Figs. 5a, 5b, 5c, respectively. The actual value to determine the depth position is averaged, taking into consideration the system 0.297-mm/pixel resolution for 25 data points/line and a phantom container diameter of . The same procedure was repeated for all images at the locations given. Figure 7 shows the displacement results for tumor tissue with a 1-cm diameter placed at about below the phantom surface. The maximum surface displacement, averaged, is about . Figure 8 shows the error found between the optical measurements and the minimum square fit to the data, which shows a straight line that serves to predict the inhomogeneity depth location. The horizontal axis shows the real, physically measured depth of the inhomogeneity immersed in the gel, while the measurements on the vertical axis are displacement averages obtained with the optical noninvasive method dealt with here. The line scans indicate that the displacement measurements performed decrease as the depth of the inhomogeneity increases. The optical measured values are in very good agreement with the actual depth position, within an error of about 0.03 %. The plot serves as a calibration chart, and thus can be used to predict from the optical data the depth of an inhomogeneity immersed in the phantom. 4.Concluding RemarksWithin the scope of this research, it is shown that 3-D pulsed digital holography is an optical noninvasive method capable of measuring with great accuracy the depth position of human tumor tissue, breast tumors, immersed in gel phantoms. Separation of the three displacement components, , , and , from unwrapped phase maps quantifies their depth. The technique was tested for different inhomogeneity known positions up to depths around below the gel surface. The results show that the optical data can be related to the depth location via a straight line fitted through the experimental results. Future research will be oriented toward the application of the method in live animals to later transfer the gained knowledge to live human tissue. The pulsed laser used delivers pulses of width at of energy for the experiments conducted here. Since the illuminating object beam is expanded, the live tissue will be exposed to an average power of , which is well within the limits of human tissue exposure to this type of laser. Furthermore, the sound power that may be used in live tissues will be around that used here, which is also well below the audible decibel limits. It is envisaged that problems will arise in the placement of the human breast in the optical setup, and the finding of the optimum resonant mode frequency directly related to the type of tissue properties, like its density and fat. AcknowledgmentsThe authors would like to acknowledge partial financial support from Consejo Nacional de Ciencia y Tecnología, grants 42971 and 48177. ReferencesS. Schedin,
G. Pedrini, and
H. J. Tiziani,
“Pulsed digital holography for deformation measurements on biological tissue,”
Appl. Opt., 39
(16), 2853
–2857
(2000). 0003-6935 Google Scholar
G. Pedrini,
H. J. Tiziani, and
M. E. Gusev,
“Pulsed digital holographic interferometry with 694- and wavelengths,”
Appl. Opt., 39
(2), 246
–249
(2000). 0003-6935 Google Scholar
J. Li,
G. Ku, and
L. V. Wang,
“Ultrasound-modulated optical tomography of biological tissue by use of contrast of laser speckles,”
Appl. Opt., 41
(28), 6030
–6035
(2002). 0003-6935 Google Scholar
B. C. Forget,
F. Ramaz,
M. Atlan,
J. Selb, and
A. C. Boccara,
“High-contrast fast Fourier transform acousto-optical tomography of phantom tissues with a frequency-chirp modulation of the ultrasound,”
Appl. Opt., 42
(7), 1379
–1383
(2003). 0003-6935 Google Scholar
S. Schedin,
G. Pedrini,
H. J. Tiziani,
A. K. Aggarwal, and
M. E. Gusev,
“Highly sensitive pulsed digital holography for built-in defect analysis with a laser excitation,”
Appl. Opt., 40
(1), 100
–103
(2001). 0003-6935 Google Scholar
M. Bashkansky,
M. D. Duncan,
M. Kahn, D. Lewis III, J. Reintjes,
“Subsurface defect detection in ceramics by high-speed high-resolution optical coherent tomography,”
Opt. Lett., 22
(1), 61
–63
(1997). 0146-9592 Google Scholar
M. del Socorro Hernández-Montes,
C. Pérez-López, and
F. Mendoza Santoyo,
“Detection of biological tissue in gels using pulsed digital holography,”
Opt. Express, 12 853
–858
(2004). https://doi.org/10.1364/OPEX.12.000853 1094-4087 Google Scholar
G. Pedrini,
F. Mendoza Santoyo,
S. Schedin,
P. Fröning, and
H. J. Tiziani,
“Whole 3D-digital holographic measurements of vibrating objects,”
Proc. SPIE, 3823 53
–63
(1999). 0277-786X Google Scholar
|