Frequency-domain reflectance for the determination of the scattering and absorption properties of tissue

Measurements of the phase and modulation of amplitudemodulated light diffusely reflected by turbid media can be used to deduce absorption and scattering coefficients.

Optical properties of normal and diseased human breast tissues in the visible and near infrared

The optical absorption and scattering coefficients have been determined for specimens of normal and diseased human breast tissues over the range of wavelengths from 500 to 1100 nm. Total attenuation coefficients were measured for thin slices of tissue cut on a microtome. The diffuse reflectance and transmittance were measured for 1.0 mm thick samples of these tissues, using standard integrating sphere techniques. Monte Carlo simulations were performed to derive the scattering and absorption coefficients, as well as the mean cosine of the scattering angle. The results indicate that scatter exceeds absorption by at least two orders of magnitude. Absorption is most significant at wavelengths below 600 nm. The scattering coefficients lie in the range 30-90 mm-1 at 500 nm, and fall smoothly with increasing wavelength to between 10 and 50 mm-1 at 1100 nm. The scattering coefficient for adipose tissue differs, in that it is invariant with wavelength over this spectral range. For all tissues examined, the scattered light is highly forward peaked, with the mean cosine of the scattering angle in the range 0.945-0.985. Systematic differences between the optical properties of some tissue types are demonstrated.

Experimental tests of the feasibility of singlet oxygen luminescence monitoring in vivo during photodynamic therapy

Singlet oxygen (1O2) is thought to be the cytotoxic agent in photodynamic therapy (PDT) with current photosensitizers. Direct monitoring of 1O2 concentration in vivo would be a valuable tool in studying biological response. Attempts were made to measure 1O2 IR luminescence during PDT of cell suspensions and two murine tumour models using the photosensitizers Photofrin II and aluminium chlorosulphonated phthalocyanine. Instrumentation was virtually identical to that devised by Parker in the one positive report of in vivo luminescence detection in the literature. Despite the fact that our treatments caused cell killing and tissue necrosis, we were unable to observe 1O2 emission under any conditions. We attribute this negative result to a reduction in 1O2 lifetime in the cellular environment. Quantitative calibration of our system allowed us to estimate that the singlet oxygen lifetime in tissue is less than 0.5 microsecond. Some technical improvements are suggested which would improve detector performance and perhaps make such measurements feasible.

In vivo tests of the concept of photodynamic threshold dose in normal rat liver photosensitized by aluminum chlorosulphonated phthalocyanine

In its simplest form, the photodynamic therapy (PDT) threshold dose model states that tissue necrosis due to PDT will occur if the number of photons absorbed by the photosensitizer per unit volume of tissue exceeds a critical value. This threshold is given by the product of photon fluence, photosensitizer concentration and specific absorption coefficient. To test the validity of this concept for PDT of normal rat liver sensitized with aluminum chlorosulphonated phthalocyanine (AISPC), all three of these parameters were varied by changing the injected AISPC dose, the wavelength of excitation and the irradiation geometry. The extent of necrosis caused by the treatment was consistent with the threshold model, except when the concentration of AISPC in the liver exceeded 20 micrograms g-1. For this animal model, we estimate the threshold to be (3.8 +/- 0.2) x 10(19) photons cm-3.

Similarity relations for the interaction parameters in radiation transport

It is possible to alter the interaction parameters characterizing neutral particle radiation transport without significantly altering the spatial distribution of the particle fluence rate. Practical mathematical relations specifying the constraints that such an alteration must satisfy are known as similarity relations. Similarity relations are derived in this work from integrated versions of the single energy neutral particle transport equation. The application of these relations in accelerating Monte Carlo deep penetration simulations is described and assessed. Computational reductions may exceed a factor of 10 in highly scattering media in which the scattering is highly forward peaked, such as applies to the propagation of red and near IR light through soft tissues.

Monte Carlo modeling of light propagation in highly scattering tissue--I: Model predictions and comparison with diffusion theory

Using optical interaction coefficients typical of mammalian soft tissues in the red and near infrared regions of the spectrum, calculations of fluence-depth distributions, effective penetration depths and diffuse reflectance from two models of radiative transfer, diffusion theory, and Monte Carlo simulation are compared for a semi-infinite medium. The predictions from diffusion theory are shown to be increasingly inaccurate as the albedo tends to zero and/or the average cosine of scatter tends to unity.

Monte Carlo modeling of light propagation in highly scattering tissues--II: Comparison with measurements in phantoms

Measurements of the fluence-depth distributions and of the diffuse reflectance of 633 nm light have been made in liquid media with optical properties similar to soft tissues. The results are compared with predictions of Monte Carlo computer calculations in order to test the adequacy of Monte Carlo modeling of light transport in tissue. Except at extremely high albedo, the experimental data and the Monte Carlo results agree well for the depth dependence of the fluence as a function of incident light beam diameter and optical absorption and scattering, and for the dependence of the diffuse reflectance on the albedo. The absolute experimental values for the fluence must be renormalized by a factor which varies with the albedo in order to match the model values, and the possible sources of this discrepancy are discussed.

Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties

When a picosecond light pulse is incident on biological tissue, the temporal characteristics of the light backscattered from, or transmitted through, the sample carry information about the optical absorption and scattering coefficients of the tissue. We develop a simple model, based on the diffusion approximation to radiative transfer theory, which yields analytic expressions for the pulse shape in terms of the interaction coefficients of a homogeneous slab. The model predictions are in good agreement with the results of preliminary in vivo experiments and Monte Carlo simulations.

Total attenuation coefficients and scattering phase functions of tissues and phantom materials at 633 nm

Measurements have been made of the total attenuation coefficient sigma t and the scattering phase function, S(theta), of 632.8 nm of light for a number of animal model tissues, blood, and inert scattering and absorbing media. Polystyrene microspheres of known size and refractive index, for which sigma t and S(theta) can be calculated by Mie theory, were used to test the experimental methods. The purpose of the study was to define typical ranges for these optical properties of tissues, as a contribution to the development of experimental and theoretical methods of light dosimetry in tissue, particularly related to photodynamic therapy of solid tumors. The results demonstrate that, for the representative tissues studied, the total attenuation coefficients are of the order of 10-100 mm-1, and that the scattering is highly forward peaked, with average cosine of scatter in the range 0.6-0.97.

Computer simulations of speckle in B-scan images

The granularity or speckle in medical ultrasound images tends to mask the presence of small lesions. As well, artifactual filling in of anechoic regions such as cysts, reduces the diagnostic potential of the images. These effects depend not only on the acoustic properties of the tissue but also are strongly influenced by the imaging system, especially the transducer geometry. To study the effect of the transducer on the final B-scan image, a computer model has been developed simulating the interaction of ultrasound with a simple scattering medium. This model, incorporating the position dependence of the point response of the transducer, is based on single scattering from a collection of points positioned randomly in a three-dimensional volume. Using this approach, B-scan images showing speckle have been generated for different transducer geometries. Inclusion of a 2.6 mm void mimicking a cyst within the three-dimensional scattering volume has allowed us to predict the cyst contrast in the image for the different transducer systems. Experimental B-scan images of a scattering phantom were obtained using the different pulse echo systems. Quantitative assessment using first and second order statistics of the images shows good agreement between experiment and theory.

The improvement and quantitative assessment of B-mode images produced by an annular array/cone hybrid

Hybrid ultrasound imaging systems, which combine spherical focusing on transmit with axicon focusing on receive, provide excellent resolution over a useful depth of field. This paper presents a new hybrid design with improved sensitivity, in which the axicon focusing is achieved by two conical mirrors and a PZT 5A disk out into 8 sectors. We have investigated two methods of processing the signals from the 8 sectors. In the first, phase insensitive sector addition (PISA), the B-scan is formed from the sum of the 8 demodulated signals. In the second, multiplicative processing (MP), the 8 rf waveforms are multiplied and the resultant is demodulated to form the image. Both techniques result in smoothed speckle but degraded lateral resolution. As well, MP decreases the off-axis sensitivity of the system and artifacts characteristic of axicon focusing. Quantitative assessment of the effects of PISA and MP was performed using a new approach called contrast-to-speckle ratio (CSR). The CSR data, which is a measure of the image contrast of cylindrical voids in a random scattering medium relative to contrast fluctuations due to speckle, shows the superiority of PISA and MP. This conclusion is supported by images of in vitro human breast tissue.

Breast imaging with a conical transducer/annular array hybrid scanner

A hybrid conical scanner consisting of a 37.5 degrees, 4 MHz conical receiver and a five element 5 cm dia. 4 MHz annular array transmitter has been developed for breast imaging. The ultrasound fields and imaging characteristics of three interesting geometries are compared. These are: (1) the annular array used in conventional pulse-echo mode, (2) the conical transducer used in combination with the full annular array and (3) the conical transducer used in combination with the center element of the annular array (simple conical scanner). Both conical systems provide extremely narrow (0.3 mm) beamwidths at -6dB but the cone/annular array combination has significantly improved sidelobe characteristics compared with the simple conical scanner. The combination of high resolution and reasonable sidelobe levels for the cone/annular array hybrid results in images with very fine speckle patterns and good definition of small structures. Breast images from young normal volunteers show that this approach has considerable potential for clinical application.

Improved method for the design of tissue compensators

A program, WBEAM, is described which calculates dose distributions in planes perpendicular to the beam axis, taking into account both field shape and patient contour. WBEAM can be used to design compensators which when placed in the beam will produce uniform dose distributions in the plane of calculation. The program was tested in three different situations: a 30 x 30 cm field with a flat contour, a "mantle" field with a flat contour, a "mantle" field with a human-like contour. The program performs as designed: the dose distributions are accurate, and the compensators flatten the radiation beams to the specified limits.

Characteristics of an 18 MV photon beam from a Therac 20 Medical Linear Accelerator

The 18 MV photon beam characteristics of a Therac 20 Medical Linear Accelerator manufactured by Atomic Energy of Canada Ltd, are presented. Tissue phantom ratios (TRP's) and percent depth dose data are given; for a 10 x 10 cm field, the percent depth dose at a depth of 10 cm is 78.5 (SSD 100 cm). The relative dose factors (RDF'S) are given and are analyzed to elucidate the relative contributions from phantom scatter, collimator scatter, and backscatter from the top of the collimators into the monitor chambers. The effect of field size and depth on the penumbra is described. Crossplots of the beam at a depth of 5 cm indicate that the flattening filter could be improved; there are hot spots of 108% near the corners of 40 x 40 fields.

Spontaneous human lymphocyte-mediated cytotoxicity against tumor target cells. IX. The quantitation of natural killer cell activity

On analysis of in vitro assays of human natural killer (NK) cell function the inadequacy of commonly used methods of expressing lytic activity was apparent. A comparison was made of the data obtained using modifications of two equations-the simple exponential fit and the von Krogh equations. Both of these equations were found to satisfy the following essential criteria for use in these assays. First, the majority of the results obtained in the chromium-release assay could be used in data reduction; second, the resultant "dose-response" curve was reduced to linearity; and third, a single numerical expression was obtained which was directly proportional to the cytotoxic activity. Of the two methods the more conventional exponential fit was found to be the simpler to use. The closeness of fit of the experimentally derived data to the ideal curves did not support the possibility that normal lymphocyte preparations contain suppressor cells capable of inhibiting NK activity. Data have also been presented showing that NK-sensitive targets could be categorized with respect to their susceptibility by comparing the slopes of the target cell survival curves obtained using the exponential fit equation. These observations are relevant to the accurate assessment of NK activity in patient populations and to the determination of the effects of disease and its treatment on this activity.

The conical scanner: a two transducer ultrasound scatter imaging technique

A new two transducer ultrasound scatter imaging system has been developed. In operation, a large conical transducer generates an ultrasound beam that converges to a sharp line focus in the tissue being imaged. The second transducer is aimed along the axis of the cone and detects scattered ultrasound as a function of time, converting this information into a high resolution image. Excellent focusing properties are maintained over a large depth of field with no increase in the complexity of standard B-scanner electronics. As well, no increase in time is required to generate the image.

The performance of a 45° conical scanner is examined theoretically and experimentally. The FWHM of the lateral point response of the scanner remains 0.25 mm over a 5 cm depth of field. Images of phantoms and in vitro breast tissue demonstrate a great improvement in image quality compared with conventional pulse-echo imaging.