Signal attenuation and localization in optical coherence tomography studied by Monte Carlo simulation.

A Monte Carlo model has been developed for optical coherence tomography (OCT). A geometrical optics implementation of the OCT probe with low-coherence interferometric detection was combined with three-dimensional stochastic Monte Carlo modelling of photon propagation in the homogeneous sample medium. Optical properties of the sample were selected to simulate intralipid and blood, representing moderately (g = 0.7) and highly (g = 0.99) anisotropic scattering respectively. For shallow optical depths in simulated intralipid (<3 scattering mean free path (mfp) units), the number of detected backscattered photons followed the extinction-single-backscatter model, and OCT was found to detect only minimally scattered photons. Within this depth range the backscatter positions of detected photons corresponded well with the nominal focus position of the probe. For propagation to deeper positions in intralipid, localization of backscattering was quickly lost due to detection of stray photons, and the number of detected photons remained constant with increasing depth in the non-absorbing medium. For strongly forward-directed scattering in simulated blood, the number of detected photons approached the extinction-single-backscatter model only for very shallow depths (<2 mfp units). However, backscattering positions for detected photons correlated well with the nominal focus position of the probe even for optical depths greater than 40 mfp units.

Accuracy and noise in optical Doppler tomography studied by Monte Carlo simulation.

A Monte Carlo model has been developed for optical Doppler tomography (ODT) within the framework of a model for optical coherence tomography (OCT). A phantom situation represented by blood flowing in a horizontal 100 microm diameter vessel placed at 250 microm axial depth in 2% intralipid solution was implemented for the Monte Carlo simulation, and a similar configuration used for experimental ODT measurements in the laboratory. Simulated depth profiles through the centre of the vessel of average Doppler frequency demonstrated an accuracy of 3-4% deviation in frequency values and position localization of flow borders, compared with true values. Stochastic Doppler frequency noise was experimentally observed as a shadowing in regions underneath the vessel and also seen in simulated Doppler frequency depth profiles. By Monte Carlo simulation, this Doppler noise was shown to represent a nearly constant level over an investigated 100 microm interval of depth underneath the vessel. The noise level was essentially independent of the numerical aperture of the detector and angle between the flow velocity and the direction of observation, as long as this angle was larger than 60 degrees. Since this angle determines the magnitude of the Doppler frequency for backscattering from the flow region, this means that the signal-to-noise ratio between Doppler signal from the flow region to Doppler noise from regions underneath the flow is improved by decreasing the angle between the flow direction and direction of observation. Doppler noise values from Monte Carlo simulations were compared with values from statistical analysis.

Wavelengths for port wine stain laser treatment: influence of vessel radius and skin anatomy.

Recent Monte Carlo computations in realistic port wine stain (PWS) models containing numerous uniformly distributed vessels suggest equal depth of vascular injury at wavelengths of 577 and 585 nm. This finding contradicts clinical experience and previous theory. From a skin model containing normal and PWS vessels in separate dermal layers, we estimate analytically the average volumetric heat production in the deepest targeted PWS vessel. The fluence rate distribution is approximated by Beer's law, which depends upon the tissue's effective attenuation coefficient, and includes a homogeneous fractional volumetric blood concentration corrected for finite-size blood vessels. The model predicts 585-587 nm wavelengths are optimal in adult PWSs containing at least one layer of small-radius blood vessels. In superficial PWSs, typically in young children with small-radius vessels, 577-580 nm wavelengths are optimal. Wavelength-independent results similar to those from Monte Carlo models are valid in single-layered PWSs of large-radius vessels. In conclusion, the volumetric heat production in the deepest targeted PWS blood vessel can be maximized on an individual patient basis. However, absorption of 585-587 nm wavelengths is sufficiently high in superficial lesions, so we hypothesize that these wavelengths may be considered adequate for the treatment of any PWS.

Doppler grid surface scanning applications for pulmonary subsurface parenchymal perfusion assessment.

Subsurface perfusion to lung parenchyma underlying the pleura is difficult to assess in live ventilated animals. The purpose of this study was to assess applicability of a newly developed laser Doppler grid scanning imaging technology that measures perfusion of pleural subsurface lung regions in intact normal and abnormal animal lungs. Eighty-six Doppler grid perfusion measurements were performed in five New Zealand White Rabbits (3-5 kg); four with unilateral bullous lung disease, one normal control. Left upper lobe lung surface was exposed to 10 1-sec spot Nd:YAG exposures (70 W/cm2). One week following laser exposure, all rabbits underwent sequential bilateral open thoracotomy. Unaffected left lower lobes in these animals and all four lobes of a previously untreated rabbit were used as controls. Pleural subsurface perfusion measurements were recorded over a contiguous 900-pixel square surface grid using quantitative noncontact laser Doppler imaging during open thoracotomy procedures. Scans were obtained in a normal volume ventilation mode, at 30 cm of inspiratory hold airway pressure, and postinflation. A perfusion-pressure response curve was obtained in normal lung at 10-, 20-, and 30-cm static airway pressure. Post mortem measurements were used as 0 flow controls. Normal lung tissue was found to have relatively high pleural subsurface perfusion (1362 +/- 328 corrected units on a scale of 0-4095). Areas of atelectasis had decreased perfusion (659 +/- 512 U., 48.4 +/- 12.5% compared to normal lung, p < 0.02), but returned to normal levels after inflation of the lung (1253 +/- 363 U., p = 0.21 compared to normal). Pleural subsurface perfusion decreased uniformly and progressively as lung inflation pressure increased (p < 0.0001). Perfusion increased immediately to supranormal values following release of high inspiratory inflation pressure holds (1603 +/- 626 U., 117 +/- 18% compared to normal lung, p = 0.03). Bullae had markedly decreased perfusion (541 +/- 68 U.) that was not further reduced by increased inflation pressures. Noncontact laser Doppler grid perfusion imaging appears to provide a new tool for measuring pleural subsurface perfusion over a large area of lung surface in clinical experimental settings. Results are rapid, reproducible, and consistent. Sampling errors inherent in current point sampling Doppler flow techniques are reduced by the multiple contiguous measurements. We have used this technique to demonstrate inspiratory pressure-related reduction in pleural subsurface perfusion in normal lung, reversible decreased perfusion in atelectatic regions, and reduced perfusion in bullous and laser-treated lung regions.