Translator Disclaimer
1 September 2008 Determination of the optical properties of turbid media by measurements of the spatially resolved reflectance considering the point-spread function of the camera system
Author Affiliations +
Abstract
The radial dependence of the diffuse reflectance from a turbid medium that is due to a point source is basically influenced by the absorption and reduced scattering coefficients. A system consisting of a HeNe laser source and a CCD camera is described for making remote measurements of the spatially resolved diffuse reflectance. Liquid tissue phantoms were made of Intralipid and trypan blue to validate the experimental setup. We show that for the correct determination of the optical properties of the tissue phantoms, the point-spread function (PSF) of the camera system has to be considered. Convolving the PSF with a solution of the diffusion equation, the absorption and reduced scattering coefficients of the tissue phantoms could be determined with an average error of 8% in the absorption coefficient and 4% in the reduced scattering coefficient, whereas without convolution, the errors were considerably larger, especially for large optical parameters.

1.

Introduction

It is of the utmost importance to understand how light is distributed in biological tissues before applying optical diagnostic methods and laser treatments. Knowledge of the propagation of light may improve safety and can lead, in general, to a better understanding of medical treatments. Light transport in biological tissue, which is an optically turbid (i.e., scattering and absorbing) medium, can be described by the Boltzmann transport equation. This involves four optical parameters: the absorption coefficient μa , the scattering coefficient μs , the refraction index n , and the scattering phase function. In the diffusion approximation to the Boltzmann equation, the phase function is replaced by the mean cosine of the scattering angle g , which is combined with the scattering coefficient to give the reduced scattering coefficient μs=μs(1g) .

Several noninvasive methods exist to determine optical properties, for example, steady-state,1, 2, 3, 4 time-resolved,5 and frequency domain6, 7 reflectance methods. In spite of many progresses achieved both with theories and with experiments, the accurate determination of the optical properties is still a nontrivial task.

In this article, we present our approach to determine the optical properties of tissue phantoms using a CCD camera system. Such cameras have been widely used to simultaneously obtain an array of digitized radiometric values from the surface of an imaged turbid medium with high linearity and precision. One can easily interpret the image as the integral of the per-pixel distribution of radiance over an accepted solid angle, but as pointed out by Huang,8 Townshend,9 and Healey,10 as a consequence of many factors including the optics of the instrument, the detector, and electronics, a substantial portion of the signal of each pixel comes from surrounding areas. This means that the real radiance of a pixel should be regarded as a function of the signals of an area surrounding it, instead of just one pixel. Because of the impact of the point-spread function (PSF) on image acquisition, it is very important to characterize the PSF even for larger distances from the central point in order to obtain accurate intensity values. In this paper, we will introduce a method of measuring the PSF of a camera system with a laser source. The effect of the PSF is demonstrated by determining the optical properties of liquid tissue phantoms. It is shown that the PSF should be taken into account in any interpretation of data obtained with an imaging system.

2.

Materials and Methods

2.1.

Determination of the Point-Spread Function

Many techniques have been employed to account for the instrumental response of the imaging system11 and to obtain the PSF of a camera system.12, 13, 14 The most direct means of measuring the PSF is to take images of a point source. The CCD camera system used in this study was a TE/CCD 512 EFT (Princeton Instruments, Trenton, New Jersey) with a 16-bit analog-to-digital converter and an array of 512×512pixels corresponding to an area of the image of approximately 37×37mm2 . The collection optics consisted of a commercial camera zoom lens (Nikkor, 50-mm focal length, Nikon, Japan). The aperture of the camera was adjusted to a medium level. The camera chip was cooled by a Peltier cooling to 35°C to minimize dark noise. Additionally, we used water cooling for the camera.

Figure 1 shows the experimental setup. The beam of a HeNe laser emitting at 633nm (LHRR-0100, Laser 2000, Germany) was coupled into an optical single-mode fiber with a core diameter of 10μm . The fiber served as a point source so that the light of the laser appears in the center of a CCD pixel. The exposure time could be varied by means of a shutter. To minimize stray light, the experiment was performed in a dark room. Unwanted light strayed and reflected by parts of the experimental setup was carefully blocked with a black cloth. The end of the fiber was positioned in the image plane of the camera system. With this setup, the strongest signal generated by the laser on the CCD sensor occupied just one pixel. Before each measurement, a dark image was recorded and subtracted.

Fig. 1

Experimental setup for the determination of the PSF. The dashed line illustrates a beam scattered by the objective of the camera system, causing the PSF not to be a delta function.

054047_1_050805jbo1.jpg

To simplify the characterization of the PSF, we assumed that it is radially symmetric across the sensor. The exposure was divided into concentric rings, and pixel values in each ring were summed and divided by the number of pixels in the ring. Pixel values from different areas but with the same radial distance from the source were averaged. Radial averaging helps to increase the signal-to-noise ratio: As the signal decreases for larger distances to the illumination point, the number of pixels averaged increases because of the increasing radius of the ring. This leads to an improvement of the signal-to-noise ratio.

2.2.

Experimental Setup for Phantom Measurements

To experimentally examine the effect of the PSF on our CCD camera system and for the proper assessment of the experimental setup, measurements on tissue phantoms were performed. The phantom should mimic biological tissue and allow easy tuning of the optical properties. Several publications concerning the production of phantoms have been reviewed.15, 16 Studies of Lipovenös (Fresenius Kabi, Germany), an intravenous nutrient containing 20% soy oil, or similar fat solutions described measurements of its scattering properties and found that the absorbance can be neglected in many cases.17, 18 As an absorber, we used trypan blue (Biochrom, Germany), an organic dye, due to its negligible scattering. We determined the absorption properties of undiluted trypan blue by measuring the collimated transmission. Both substances were mixed with distilled water. Changing the concentrations of the added substances, the optical properties of the obtained solutions could be varied. We made up the tissue phantoms in a cylindrical black bin with a height of 18cm and a diameter of 10cm , so the phantoms could be regarded to be a semi-infinite half-space. The index of refraction of the phantom was assumed to be n=1.33 , because the phantom consists almost entirely of water. An analytic expression for the reflectance derived by means of the diffusion model was used for the determination of the optical properties.19

The experimental setup is provided in Fig. 2 . Again the same HeNe laser emitting at 633nm served as a light source. The laser beam was deflected onto the specimen by a mirror at an angle of incidence of 10degto15deg to avoid the detection of specularly reflected light. We confirmed by Monte Carlo simulations that this small angle had a negligible influence on the spatially resolved reflectance. The diffusely reflected light was detected by the same CCD camera mounted perpendicular to the specimen surface as described in Sec. 2.1. The distance between the front of the camera and the sample was fixed at 288mm . Neutral density filters were used to adjust the incident laser power to avoid saturation of the detector system. Again, we assumed that the reflectance is radially symmetric, and thus, we divided the exposure into concentric rings. Pixel values in each ring were summed and divided by the number of pixels in the ring to improve the signal-to-noise ratio.

Fig. 2

Experimental setup for phantom measurements.

054047_1_050805jbo2.jpg

When directing the beam into the middle of the objective of the camera and fixing the distance between the laser and the camera to 288mm , as mentioned in Sec. 2.1, we were able to determine the diameter of the laser beam. Referring to an intensity of (1e)37% of the maximum value, we received a beam diameter of 0.15mm .

2.3.

Data Analysis

Since the PSF is not just a delta function, effects of the PSF on the intensity values far away from the brightest pixel cannot be ignored without appreciable radiometric consequences. In a real optical system, images obtained are actually the convolution of the pixel intensity values with the normalized PSF of the system. We carried out a two-dimensional (2-D) convolution algorithm based on fast Fourier transform (FFT) and inverse FFT techniques. The Fourier transform of the convoluted raw image of the camera, I¯raw , equals the product of the Fourier transforms of the normalized PSF, I¯PSF , and the true image, I¯true ,

Eq. 1

I¯raw(x,y)=I¯true(x,y)I¯PSF(x,y).
In practice, because we did not have an exact formula for the PSF but only a chart, the PSF and the intensity values of the camera image were all expressed as M×M matrices. We expanded the intensity values calculated according to Sec. 2.1 to two dimensions instead of using the measured intensity values, because this approach clearly improves the signal-to-noise ratio. The origin of the radial distance r was set at the (0, 0) element of the M×M matrix.

The absorption and the reduced scattering coefficients were obtained by (1) fitting the experimental data to an analytical solution of the diffusion equation for a homogenous infinite turbid medium according to Kienle and Patterson,19 and (2) fitting the 2-D convolution of the preceding solution of the diffusion equation and of the PSF to the experimental data. We used a convolution algorithm because a deconvolution process tends to increase the noise level in the image and so enhances the contrast of both the signal and the noise. The extrapolated boundary condition was regarded. The fitting range contained only distances larger than 1μs because the diffusion equation is not valid for distances shorter than one transport mean free path.19 It included three orders of magnitude from this starting distance. For large μs , the fitting range included only distances larger than 0.8mm because of the finite beam diameter used in the experiments. If the optical parameters of the phantoms were roughly known, the starting parameters were chosen to vary at about 30 to 40% from the expected values. If the optical parameters were completely unknown, we started with μa=0.01mm1 and μs=1.0mm1 . As the obtained optical parameters depend on the fitting range and the fitting range—as described earlier—depends on the starting value of μs , the fitting algorithm sometimes had to be accomplished twice or triply. This was especially the case if the starting parameter of μs differed strongly from the real value.

The best fit was reached with a Levenberg-Marquardt algorithm by varying μa and μs , a multiplicative constant due to the unknown conversion of pixel counts to reflectance values, and an additive constant because of the drifting of the camera background. In the fitting procedure, we assumed a Poisson-like distribution of noise, so the weight of every value was calculated by its inverse intensity value measured by the camera. The optical properties could be found most accurately in this way.

3.

Results and Discussion

3.1.

Point-Spread Function of the CCD Camera System

The first measurements were performed without using a neutral density filter. The exposure time was set to obtain maximum counts in the central bright spot without saturating the pixel. We used t=0.01s in this case. The bright spot had digital counts on the order of 3104 . Pixels immediately surrounding the bright spot consisted of intensities of 103 digital counts or lower. The challenge was that outside of the few bright pixels, the signal level was low, which led to a low signal-to-noise ratio. It was, therefore, necessary to raise digital counts in that area. Thus, another measurement had to be accomplished with a longer exposure time (t=1s) , which resulted in saturation of the central pixels. This led to an improvement of the signal-to-noise ratio for the pixels far away from the point source but also increased the intensity signal of the whole curve. For further interpretation, we took only the part of the curve with intensity values below saturation into account.

Because of the uncertain effective exposure time, i.e., the unknown opening time of the shutter, it is inevitable to determine the exact factor between the opening times of both measurements. This can be achieved by comparing two nonsaturated measurements with different shutter opening times. As a result, we received the PSF with a high signal-to-noise ratio even for large distances to the illumination point.

We determined the PSF in this way with the point source illuminating several different positions of the CCD chip. Small variations between the single PSFs were found only for large radial distances, i.e., distances between the image of the point source and the other pixels larger than 20mm , resulting from the drifting of the camera’s dark noise. Despite the fact that we took an associating dark image immediately before the real measurement, it is not possible to determine the low intensity values exactly. Figure 3 shows the average of five different PSFs, which have been determined for five different positions of the CCD chip. As can be seen, relatively large error bars appear for distances far away from the illumination point, whereas for small distances, the errors are smaller.

Fig. 3

Mean PSF measured for five different positions on the CCD chip. Error bars show standard deviation. Errors bars are not centered because of logarithmic scaling.

054047_1_050805jbo3.jpg

Therefore, careful characterization of the camera’s PSF is essential because the accuracy of the PSF will affect the accuracy of the received optical coefficients. Each optical system will have a characteristic PSF that must be measured precisely. The PSF tends to spread energy everywhere and lowers the contrast. Images and boundaries of bright objects will be seriously affected by the PSF. This can make per-pixel interpretation erroneous because pixels that map places on the phantom surface where no light is remitted can be increased by a bias. The larger the extension of the source is, the worse is the bias. One way to remove the artificial radiance variation is to deconvolute the image with the system’s PSF and return the intensity values back to the true values.

3.2.

Optical Parameters of Phantom Measurements

As described earlier, we assumed that in liquid tissue phantoms, Lipovenös is the only scatterer and trypan blue is the only absorber; however, for very small concentrations of trypan blue, the absorption of water, which is μa0.0003mm1 at 633nm , has to be regarded. We investigated phantoms with optical properties between μa0.0030.1mm1 and μs×=0.552.0mm1 , which are typical for biological tissues. Table 1 shows the results for the fit of the experimental data with the solution of the diffusion equation for an infinite homogenous turbid medium both when using and without using a convolution algorithm.

Table 1

Expected and calculated optical parameters of tissue phantoms without and with a convolution algorithm. The absorption coefficient of water, which is ≈0.0003mm−1 at 633nm , has been considered

Expected ParametersFit without ConvolutionFit with Convolution
μa(mm1) μs(mm1) μa(mm1) μs(mm1) μa(mm1) μs(mm1)
0.003350.550.003290.560.003240.56
0.006060.550.005410.570.005540.58
0.01080.550.01130.530.01050.55
0.03120.550.03180.550.03240.56
0.05720.550.05980.550.05870.57
0.09780.550.09150.600.09280.59
0.003101.100.002661.110.002531.06
0.005871.100.005931.080.005701.09
0.009861.100.009471.110.009431.10
0.02991.100.02711.140.02861.10
0.05981.100.05271.140.05811.11
0.09911.100.08261.190.09561.14
0.003142.000.002562.290.002941.97
0.005742.000.004352.380.006171.98
0.01062.000.007952.390.01131.99
0.03162.000.02062.450.02992.01
0.05922.000.04092.510.06042.02
0.10032.000.06292.610.09912.02

In general, we found that the reduced scattering coefficient can be determined more accurately than the absorption coefficient for the range stated earlier. For both small absorption and reduced scattering coefficients, we found good agreement between the expected values and the optical parameters received with both fitting algorithms (rms error in μa6% , rms error in μs ).

On the other hand, as can be seen in Table 1, for large optical parameters, the fit without the convolution algorithm leads to systematically decreased absorption coefficients and systematically increased reduced scattering coefficients. For μa=0.1003mm1 and μs=2.00mm1 , the errors were 38% in μa and 31% in μs . Figure 4 shows a comparison of the measured reflectance and the intensities for both fitting algorithms. As can be seen, especially the part for distances larger than 8mm was not fitted well without using the convolution algorithm. The additive constant was overestimated in this case, which led to a constant intensity for distances larger than 11mm . So the bending of the curve cannot be described accurately in that way. However, as the measured intensity decreases slightly for larger distances, the systematically falsified optical parameters are explainable.

Fig. 4

Measured reflectance and fit to measured data with and without using a convolution algorithm for μa0.10mm1 , μs2.0mm1 . The fitting range starts at 0.8mm and includes three orders of magnitude.

054047_1_050805jbo4.jpg

Only when considering the specific influence of the objective by convoluting the data with the PSF could a significant improvement be achieved, and systematic discrepancies could be eliminated. The measured data could be fitted quite exactly (see Fig. 4) The convoluted curve now agrees with the experimental data both for small and for larger distances to the central beam. The slight decrease of the intensity for distances larger than 10mm is recorded quite well.

Again four fitting parameters had to be used. Despite the consideration of light that is scattered by the objective of the camera with the help of the convolution algorithm, we found that the use of an additive constant is essential for correct results of the optical parameters, although the additive constant takes over a fairly small value in this case. The dark noise of the CCD camera is not stable temporally despite cooling down the camera to 35°C . Although before all measurements, a camera background was determined and these background data were subtracted from the immediately following measurements, the drift of the dark noise led to nonsatisfying results. Thus, the additive fitting parameter accounts for the instabilities of the dark noise. In conclusion, the fitting algorithm with both a variable additive and a variable multiplicative constant allows accurate determination of the optical properties for tissue phantoms in the range of typical biological values.

It should be added that if a camera system with a stable background exists, it would be possible to accomplish the fitting algorithm with only three fitting parameters, neglecting the variable additive constant. If so, the convolution of the measured intensity values with the exact PSF would be even more important because the fitting algorithm is less flexible with only three fitting parameters.

With our approach, we found good agreement for the whole range of optical parameters with rms errors in μa of 8% and rms errors in μs of 4% (see Table 1). Relatively large errors appearing, e.g., for μa=0.0978mm1 , μs=0.55mm1 are probably due to diffusion theory, which requires μaμs . This assumption is not fulfilled here.

We also accomplished determination of the optical properties of tissue phantoms with slightly varying intensities of the PSF for larger distances to the illumination point but within the standard deviation (see Fig. 3). The calculated optical parameters differ only marginally with an rms error in μa of 11% and an rms error in μs of 5%.

To check the consistency of the method, used we accomplished a blind study with 12 tissue phantoms and their optical parameters extending the whole range of Table 1. When running the fitting algorithm a few times (this had to be done because of the completely unknown reduced scattering coefficient, the invalidity of the diffusion equation for distances smaller than 1μs , and the dependence of the beginning of the fitting range on the value chosen for μs at the beginning; see Sec. 2.3), the results were similar with errors in μa of 10% and errors in μs of 5%.

4.

Conclusion

In conclusion, we have successfully determined the PSF of a CCD camera system and showed that the optical parameters of liquid tissue phantoms can be determined accurately in the whole parameter range only with the help of the PSF. A fit using a convolution algorithm technique based on fast Fourier transform/inverse fast Fourier transform was described. For the fitting algorithm, four fitting parameters had to be used. We have seen that the characterization of the PSF, even for large distances to the illumination point, is an indispensable part of measurement utilizing a camera system for large μa and μs values. So the whole system provides a basis for checking and determining optical properties of unknown scattering media.

References

1. 

T. J. Farrell, M. S. Patterson, and B. Wilson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys., 19 879 –888 (1992). https://doi.org/10.1118/1.596777 0094-2405 Google Scholar

2. 

A. Kienle, L. Lilge, M. S. Patterson, R. Hibst, R. Steiner, and B. C. Wilson, “Spatially resolved absolute diffuse reflectance measurements for noninvasive determination of optical scattering and absorption coefficients of biological tissue,” Appl. Opt., 35 2304 –2314 (1996). 0003-6935 Google Scholar

3. 

M. G. Nichols, E. L. Hull, and T. H. Foster, “Design and testing of a white-light, steady-state diffuse reflectance spectrometer for determination of optical properties of highly scattering systems,” Appl. Opt., 36 93 –104 (1997). 0003-6935 Google Scholar

4. 

N. Dögnitz and G. Wagnierés, “Determination of tissue optical properties by steady-state spatial frequency-domain reflectometry,” Lasers Med. Sci., 13 55 –65 (1998). 0268-8921 Google Scholar

5. 

M. S. Patterson, B. Chance, and B. C. Wilson, “Time-resolved reflectance and transmittance for the noninvasive measurement of tissue optical properties,” Appl. Opt., 28 2331 –2336 (1989). 0003-6935 Google Scholar

6. 

M. S. Patterson, J. D. Moulton, B. C. Wilson, K. W. Berndt, and J. R. Lackowicz, “Frequency-domain reflectance for determination of the scattering and absorption properties of tissue,” Appl. Opt., 30 4474 –4476 (1991). 0003-6935 Google Scholar

7. 

S. R. Arridge, M. Cope, and D. T. Delpy, “The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis,” Phys. Med. Biol., 37 1531 –1560 (1992). https://doi.org/10.1088/0031-9155/37/7/005 0031-9155 Google Scholar

8. 

C. Huang, J. R. G. Townshend, S. N. V. Kalluri, and R. S. de Fries, “Impact of sensor’s point spread function on land cover characterization: assessment and deconvolution,” Remote Sens. Environ., 80 203 –212 (2002). 0034-4257 Google Scholar

9. 

J. R. G. Townshend, C. Huang, S. N. Kalluri, R. S. de Fries, S. Liang, and K. Yang, “Beware of per-pixel characterization of land cover,” Int. J. Remote Sens., 21 839 –843 (2000). 0143-1161 Google Scholar

10. 

G. E. Healey and R. Kondepudy, “Radiometric CCD camera calibration and noise estimation,” IEEE Trans. Pattern Anal. Mach. Intell., 16 267 –276 (1994). https://doi.org/10.1109/34.276126 0162-8828 Google Scholar

11. 

T. H. Pham, C. Eker, A. Durkin, B. J. Tromberg, and S. Andersson-Engels, “Quantifying the optical properties and chromophore concentrations of turbid media by chemometric analysis of hyperspectral diffuse reflectance data collected using a Fourier interferometric imaging system,” Appl. Spectrosc., 55 1035 –1045 (2001). https://doi.org/10.1366/0003702011952938 0003-7028 Google Scholar

12. 

A. L. Weickenmeier, W. Nüchter, and J. Mayer, “Quantitative characterization of point-spread function and detection quantum efficiency for a YAG scintillator slow scan CCD camera,” Optik (Jena), 99 147 –154 (1995). 0030-4026 Google Scholar

13. 

J. Lehr, J. B. Sibarita, and J. M. Chassery, “Image restoration in x-ray microscopy: PSF determination and biological applications,” IEEE Trans. Image Process., 7 258 –263 (1998). https://doi.org/10.1109/83.661006 1057-7149 Google Scholar

14. 

J. Zandhuis, D. Pycock, S. Quigley, and P. Webb, “Sub-pixel nonparametric PSF estimation for image enhancement,” IEE Proc. Vision Image Signal Process., 14 285 –292 (1997). 1350-245X Google Scholar

15. 

R. Cubeddu, A. Pifferi, P. Taroni, A. Torricelli, and G. Valentini, “A solid tissue phantom for photon migration studies,” Phys. Med. Biol., 42 1971 –1979 (1997). https://doi.org/10.1088/0031-9155/42/10/011 0031-9155 Google Scholar

16. 

B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt., 11 041102 (2006). https://doi.org/10.1117/1.2335429 1083-3668 Google Scholar

17. 

S. T. Flock, S. L. Jacques, B. C. Wilson, W. M. Star, and M. J. C. van Gemert, “Optical properties of Intralipid: a phantom medium for light propagation studies,” Lasers Surg. Med., 12 510 –519 (1992). https://doi.org/10.1002/lsm.1900120510 0196-8092 Google Scholar

18. 

H. J. van Staveren, C. J. M. Moes, J. van Marle, S. A. Prahl, and M. J. C. van Gemert, “Light scattering in Intralipid-10% in the wavelength range of 4001100nm,” Appl. Opt., 30 4507 –4514 (1991). 0003-6935 Google Scholar

19. 

A. Kienle and M. S. Patterson, “Improved solutions of the steady-state and the time-resolved diffusion equations for the reflectance from a semi-infinite turbid medium,” J. Opt. Soc. Am. A, 14 246 –254 (1996). https://doi.org/10.1364/JOSAA.14.000246 0740-3232 Google Scholar
©(2008) Society of Photo-Optical Instrumentation Engineers (SPIE)
Marco Pilz, Sibylle Honold, and Alwin Kienle "Determination of the optical properties of turbid media by measurements of the spatially resolved reflectance considering the point-spread function of the camera system," Journal of Biomedical Optics 13(5), 054047 (1 September 2008). https://doi.org/10.1117/1.2983675
Published: 1 September 2008
JOURNAL ARTICLE
6 PAGES


SHARE
Advertisement
Advertisement
Back to Top