Poly(vinyl alcohol) gels as photoacoustic breast phantoms revisited

A popular phantom in photoacoustic imaging is poly(vinyl alcohol) (PVA) hydrogel fabricated by freezing and thawing (F-T) aqueous solutions of PVA. The material possesses acoustic and optical properties similar to those of tissue. Earlier work characterized PVA gels in small test specimens where temperature distributions during F-T are relatively homogeneous. In this work, in breast-sized samples we observed substantial temperature differences between the shallow regions and the interior during the F-T procedure. We investigated whether spatial variations were also present in the acoustic and optical properties. The speed of sound, acoustic attenuation, and optical reduced scattering coefficients were measured on specimens sampled at various locations in a large phantom. In general, the properties matched values quoted for breast tissue. But while acoustic properties were relatively homogeneous, the reduced scattering was substantially different at the surface compared with the interior. We correlated these variations with gel microstructure inspected using scanning electron microscopy. Interestingly, the phantom's reduced scattering spatial distribution matches the optical properties of the standard two-layer breast model used in x ray dosimetry. We conclude that large PVA samples prepared using the standard recipe make excellent breast tissue phantoms.     

Mechanical property assessment of tissue-mimicking phantoms using remote palpation and optical read-out for amplitude of vibration and refractive index modulation

A coherent light beam is used to interrogate the focal region within a tissue-mimicking phantom insonified by an ultrasound transducer. The ultrasound-tagged photons exiting from the object carry with them information on local optical path length fluctuations caused by refractive index variations and medium vibration. Through estimation of the force distribution in the focal region of the ultrasound transducer, and solving the forward elastography problem for amplitude of vibration of tissue particles, we observe that the amplitude is directed along the axis of the transducer. It is shown that the focal region interrogated by photons launched along the transducer axis carries phase fluctuations owing to both refractive index variations and particle vibration, whereas the photons launched perpendicular to the transducer axis carry phase fluctuations arising mainly from the refractive index variations, with only smaller contribution from vibration of particles. Monte-Carlo simulations and experiments done on tissue-mimicking phantoms prove that as the storage modulus of the phantom is increased, the detected modulation depth in autocorrelation is reduced, significantly for axial photons and only marginally for the transverse-directed photons. It is observed that the depth of modulation is reduced to a significantly lower and constant value as the storage modulus of the medium is increased. This constant value is found to be the same for both axial and transverse optical interrogation. This proves that the residual modulation depth is owing to refractive index fluctuations alone, which can be subtracted from the overall measured modulation depth, paving the way for a possible quantitative reconstruction of storage modulus. Moreover, since the transverse-directed photons are not significantly affected by storage modulus variations, for a quantitatively accurate read-out of absorption coefficient variation, the interrogating light should be perpendicular to the focusing ultrasound transducer axis.     

Producing homogeneous cryogel phantoms for medical imaging: a finite-element approach

Tissue-mimicking phantoms with well-defined properties can help in identifying the potential weaknesses in medical imaging systems. Among the imaging systems, magnetic resonance elastography is a new noninvasive technique used to quantify the shear modulus of biological tissues, and therefore has shown promise in studying liver and brain pathologies. Polyvinyl alcohol (PVA) cryogel prepared by the freeze-thaw technique is a potential candidate for mimicking the mechanical properties of soft tissues and has been extensively used as a phantom material. However, large PVA cryogels suffer from variations in properties, partly due to the low thermal conductivity of PVA solution. The loss of homogeneity in cryogel phantoms is also attributed to inhomogeneous thawing rates during the freeze-thaw cycle. We have used a modified freeze-thaw process that imposes multiple isotherms so as to enhance the homogeneity of the produced cryogels. In addition, we have developed a finite-element modeling tool (a virtual controller) to optimize the temperature profile during the freeze-thaw cycle. Our experimental validations demonstrated the potential of the virtual controller in predicting the optimal temperature profile for the freeze-thaw process (phantom diameters: 60 and 100 mm). A robust simulation framework can fill the gap in the scientific literature with regard to phantom design for medical imaging and will help to reduce phantom development time and cost.     

Poly(vinyl alcohol) gels for use as tissue phantoms in photoacoustic mammography

Materials for solid photoacoustic breast phantoms, based on poly(vinyl alcohol) hydrogels, are presented. Phantoms intended for use in photoacoustics must possess both optical and acoustic properties of tissue. To realize the optical properties of tissue, one approach was to optimize the number of freezing and thawing cycles of aqueous poly(vinyl alcohol) solutions, a procedure which increases the turbidity of the gel while rigidifying it. The second approach concentrated on forming a clear matrix of the rigid poly(vinyl alcohol) gel without any scattering, so that appropriate amounts of optical scatterers could be added at the time of formation, to tune the optical properties as per requirement. The relevant optical and acoustic properties of such samples were measured to be close to the average properties of human breast tissue. Tumour simulating gel samples of suitable absorption coefficient were created by adding appropriate quantities of dye at the time of formation; the samples were then cut into spheres. A breast phantom embedded with such 'tumours' was developed for studying the applicability of photoacoustics in mammography.     

Optical scattering, absorption, and polarization of healthy and neovascularized human retinal tissues

The optical scattering, absorption, and polarization properties of human retinal tissues are investigated for a number of laser wavelengths in the visible range. The indices of refraction of these tissues are determined by applying Brewster's law. The inverse adding doubling method based on the diffusion approximation and radiative transport theory is applied to the measured values of total diffuse transmission, total diffuse reflection, and index of refraction to determine the optical absorption, scattering, and scattering anisotropy coefficients of the intact retinal tissues from healthy and diseased (neovascularized) human eyes. The polarization studies show that the retinal tissues possess significant intrinsic polarization characteristics, that are more pronounced in diseased tissues than in healthy tissues.     

Measurement of visco-elastic properties of breast-tissue mimicking materials using diffusing wave spectroscopy

Diffusing wave spectroscopy (DWS), without the use of tracer particles, has been used to study the internal dynamics of polyvinyl alcohol (PVA) phantoms, which mimic the properties of normal and malignant breast tissues. From the measured intensity autocorrelations, the mean square displacement (MSD) of phantom meshing is estimated, leading to the storage and loss moduli of the medium covering frequencies up to 10 KHz. These are verified with independent measurements from a dynamic mechanical analyzer (DMA) at low frequencies. We thus prove the usefulness of DWS to extract visco-elastic properties of the phantom and its possible application in detecting malignancy in soft tissues.     

Model-based elastography: a survey of approaches to the inverse elasticity problem

Elastography is emerging as an imaging modality that can distinguish normal versus diseased tissues via their biomechanical properties. This paper reviews current approaches to elastography in three areas--quasi-static, harmonic and transient--and describes inversion schemes for each elastographic imaging approach. Approaches include first-order approximation methods; direct and iterative inversion schemes for linear elastic; isotropic materials and advanced reconstruction methods for recovering parameters that characterize complex mechanical behavior. The paper's objective is to document efforts to develop elastography within the framework of solving an inverse problem, so that elastography may provide reliable estimates of shear modulus and other mechanical parameters. We discuss issues that must be addressed if model-based elastography is to become the prevailing approach to quasi-static, harmonic and transient elastography: (1) developing practical techniques to transform the ill-posed problem with a well-posed one; (2) devising better forward models to capture the complex mechanical behavior of soft tissues and (3) developing better test procedures to evaluate the performance of modulus elastograms.     

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.     

Light scattering properties of natural and artificially demineralized dental enamel at 1310 nm

A fundamental understanding of how near-IR light propagates through sound and carious dental hard tissues is essential for the development of clinically useful optical diagnostic systems, since image contrast is based on changes in the optical properties of these tissues on demineralization. During the caries (decay) process, micropores are formed in the lesion due to partial dissolution of the individual mineral crystals. Such small pores behave as scattering centers, strongly scattering visible and near-IR light. The optical properties of enamel can be quantitatively described by the absorption and scattering coefficients, and the scattering phase function. Our aim is to measure the optical scattering behavior of natural and artificial enamel caries. Near-IR attenuation measurements and angular-resolved goniometer measurements coupled with Monte Carlo simulations are used to determine changes in the scattering coefficient and the scattering anisotropy on demineralization at 1310 nm. An ultra-high resolution digital microradiography system is used to quantify the lesion severity by measurement of the relative mineral loss for comparison with optical scattering measurements. The scattering coefficient increases exponentially with increasing mineral loss. Natural and artificial demineralization increases the scattering coefficient more than two orders of magnitude at 1310 nm, and the scattering is highly forward directed.     

Acoustical properties of selected tissue phantom materials for ultrasound imaging

This note summarizes the characterization of the acoustic properties of four materials intended for the development of tissue, and especially breast tissue, phantoms for the use in photoacoustic and ultrasound imaging. The materials are agar, silicone, polyvinyl alcohol gel (PVA) and polyacrylamide gel (PAA). The acoustical properties, i.e., the speed of sound, impedance and acoustic attenuation, are determined by transmission measurements of sound waves at room temperature under controlled conditions. Although the materials are tested for application such as photoacoustic phantoms, we focus here on the acoustic properties, while the optical properties will be discussed elsewhere. To obtain the acoustic attenuation in a frequency range from 4 MHz to 14 MHz, two ultrasound sources of 5 MHz and 10 MHz core frequencies are used. For preparation, each sample is cast into blocks of three different thicknesses. Agar, PVA and PAA show similar acoustic properties as water. Within silicone polymer, a significantly lower speed of sound and higher acoustical attenuation than in water and human tissue were found. All materials can be cast into arbitrary shapes and are suitable for tissue-mimicking phantoms. Due to its lower speed of sound, silicone is generally less suitable than the other presented materials.     

Refractive index of carcinogen-induced rat mammary tumours

Near-infrared optical techniques for clinical breast cancer screening in humans are rapidly advancing. Based on the computational inversion of the photon diffusion process through the breast, these techniques rely on optical tissue models for accurate image reconstruction. Recent interest has surfaced regarding the effect of refractive index variations on these reconstructions. Although many data exist regarding the scattering and absorption properties of normal and diseased tissue, no measurements of refractive index appear in the literature. In this paper, we present near-infrared refractive index data acquired from N-methyl-N-nitrosourea-induced rat mammary tumours, which are similar in pathology and disease progression to human ductal carcinoma. Eight animals, including one control, were employed in this study, yielding data from 32 tumours as well as adjacent adipose and connective tissues.     

Measurement of optical scattering properties with low-coherence enhanced backscattering spectroscopy

Low-coherence enhanced backscattering (LEBS) is a depth selective technique that allows noninvasive characterization of turbid media such as biological tissue. LEBS provides a spectral measurement of the tissue reflectance distribution as a function of distance between incident and reflected ray pairs through the use of partial spatial coherence broadband illumination. We present LEBS as a new depth-selective technique to measure optical properties of tissue in situ. Because LEBS enables measurements of reflectance due to initial scattering events, LEBS is sensitive to the shape of the phase function in addition to the reduced scattering coefficient (μ(s) (*)). We introduce a simulation of LEBS that implements a two parameter phase function based on the Whittle-Mate?rn refractive index correlation function model. We show that the LEBS enhancement factor (E) primarily depends on μ(s) (*), the normalized spectral dependence of E (S(n)) depends on one of the two parameters of the phase function that also defines the functional type of the refractive index correlation function (m), and the LEBS peak width depends on both the anisotropy factor (g) and m. Three inverse models for calculating these optical properties are described and the calculations are validated with an experimental measurement from a tissue phantom.     

Elastic moduli of normal and pathological human breast tissues: an inversion-technique-based investigation of 169 samples

Understanding and quantifying the mechanical properties of breast tissues has been a subject of interest for the past two decades. This has been motivated in part by interest in modelling soft tissue response for surgery planning and virtual-reality-based surgical training. Interpreting elastography images for diagnostic purposes also requires a sound understanding of normal and pathological tissue mechanical properties. Reliable data on tissue elastic properties are very limited and those which are available tend to be inconsistent, in part as a result of measurement methodology. We have developed specialized techniques to measure tissue elasticity of breast normal tissues and tumour specimens and applied them to 169 fresh ex vivo breast tissue samples including fat and fibroglandular tissue as well as a range of benign and malignant breast tumour types. Results show that, under small deformation conditions, the elastic modulus of normal breast fat and fibroglandular tissues are similar while fibroadenomas were approximately twice the stiffness. Fibrocystic disease and malignant tumours exhibited a 3-6-fold increased stiffness with high-grade invasive ductal carcinoma exhibiting up to a 13-fold increase in stiffness compared to fibrogalndular tissue. A statistical analysis showed that differences between the elastic modulus of the majority of those tissues were statistically significant. Implications for the specificity advantages of elastography are reviewed.     

Application of ultrasound-tagged photons for measurement of amplitude of vibration of tissue caused by ultrasound: theory, simulation, and experiments

We investigate the modulation of an optical field caused by its interaction with an ultrasound beam in a tissue mimicking phantom. This modulation appears as a modulation in the intensity autocorrelation, which is measured by a photon counting correlator. The factors contributing to the modulation are: 1. amplitude of vibration of the particles of the tissue, 2. refractive index modulation, and 3. absorption coefficient in the region of the tissue intercepted by the ultrasound beam and light. We show in this work that a significant part of the contribution to this modulation comes from displacement of the tissue particles, which in turn is governed by the elastic properties of the tissue. We establish, both through simulations and experiments using an optical elastography phantom, the effects of the elasticity and absorption coefficient variations on the modulation of intensity autocorrelation. In the case where there is no absorption coefficient variation, we suggest that the depth of modulation can be calibrated to measure the displacement of tissue particles that, in turn, can be used to measure the tissue elasticity.     

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.     

High-resolution tensor MR elastography for breast tumour detection

MR elastography is a novel imaging technique for the visualization of elastic properties of tissue. It is expected that this method will have diagnostic value for the clarification of suspicious breast lesions. Low-frequency mechanical waves are coupled into the tissue and visualized via an MR sequence which is phase-locked to the mechanical excitation. Commonly, elasticity is assumed to be isotropic and reconstruction is performed in only two dimensions. The technique is extended to three dimensions such that the entire symmetric elasticity tensor is assessed. This is achieved by measuring different phases of the mechanical wave during one oscillatory cycle. Thereby it is possible to provide information about the anisotropy of the elasticity tensor. Finite-element simulations as well as phantom experiments are performed to demonstrate the feasibility of the method. Initial clinical results of a breast carcinoma are presented. The analysis of the eigenvalues of the elasticity tensor support the hypothesis that breast carcinoma might exhibit an anisotropic elasticity distribution. The surrounding benign tissue appears isotropic. Thereby new and additional diagnostic information is provided which might help in distinguishing between benign and malignant breast diseases.     

Reconstruction of optical properties of phantom and breast lesion in vivo from paraxial scanning data

We report on the reconstruction of absorption and reduced scattering coefficients of breast tissue in vivo of a patient with mastopathic disease. Distributions of times of flight of photons through the compressed breast were recorded by paraxial scanning. From data measured at four different source-detector offsets optical properties were reconstructed within the linear Rytov approximation by a fast inverse Fourier space method. Low-pass filtering in Fourier space was employed to remove excessive noise from high spatial frequency components and to reduce the computational efforts by a factor of 3, typically. The mammograms displaying reconstructed absorption and reduced scattering coefficients were compared with projection mammograms either obtained by time-window analysis of experimental data or based on average absorption and reduced scattering coefficients which were derived from measured temporal point spread functions within a simple homogeneous model. All inhomogeneities which were visible in the projection mammograms and which could be associated with specific breast tissue compartments could be correlated with inhomogeneities in the reconstructed absorption coefficient. In particular, the mastopathic disease was detected in the reconstructed absorption mammogram. In order to assess reliability of optical properties reconstructed from data obtained by paraxial scanning, corresponding phantom experiments and reconstructions of phantom optical properties were carried out. Because of the limited angular range sampled by the in vivo and phantom measurements, considerable blurring of the absorption coefficient occurs along the compression direction, compromising longitudinal resolution.     

Optical and acoustic properties at 1064 nm of polyvinyl chloride-plastisol for use as a tissue phantom in biomedical optoacoustics

A novel optoacoustic phantom made of polyvinyl chloride-plastisol (PVCP) for optoacoustic studies is described. The optical and acoustic properties of PVCP were measured. Titanium dioxide (TiO2) powder and black plastic colour (BPC) were used to introduce scattering and absorption, respectively, in the phantoms. The optical absorption coefficient (mua) at 1064 nm was determined using an optoacoustic method, while diffuse reflectance measurements were used to obtain the optical reduced scattering coefficient (mu's). These optical properties were calculated to be mua = (12.818 +/- 0.001)ABPC cm(-1) and mu's = (2.6 +/- 0.2)S(TiO2) + (1.4 +/- 0.1) cm(-1), where ABPC is the BPC per cent volume concentration, and S(TiO2) is the TiO2 volume concentration (mg mL(-1)). The speed of sound in PVCP was measured to be (1.40 +/- 0.02) x 10(3) m s(-1) using the pulse echo transmit receive method, with an acoustic attenuation of (0.56 +/- 1.01) f(1.51+/-0.06)MHz (dB cm(-1)) in the frequency range of 0.61-1.25 MHz, and a density, calculated by measuring the displacement of water, of 1.00 +/- 0.04 g cm(-3). The speed of sound and density of PVCP are similar to tissue, and together with the user-adjustable optical properties, make this material well suited for developing tissue-equivalent phantoms for biomedical optoacoustics.     

Magnetic resonance acoustic radiation force imaging

Acoustic radiation force impulse imaging is an elastography method developed for ultrasound imaging that maps displacements produced by focused ultrasound pulses systematically applied to different locations. The resulting images are "stiffness weighted" and yield information about local mechanical tissue properties. Here, the feasibility of magnetic resonance acoustic radiation force imaging (MR-ARFI) was tested. Quasistatic MR elastography was used to measure focal displacements using a one-dimensional MRI pulse sequence. A 1.63 or 1.5 MHz transducer supplied ultrasound pulses which were triggered by the magnetic resonance imaging hardware to occur before a displacement-encoding gradient. Displacements in and around the focus were mapped in a tissue-mimicking phantom and in an ex vivo bovine kidney. They were readily observed and increased linearly with acoustic power in the phantom (R2=0.99). At higher acoustic power levels, the displacement substantially increased and was associated with irreversible changes in the phantom. At these levels, transverse displacement components could also be detected. Displacements in the kidney were also observed and increased after thermal ablation. While the measurements need validation, the authors have demonstrated the feasibility of detecting small displacements induced by low-power ultrasound pulses using an efficient magnetic resonance imaging pulse sequence that is compatible with tracking of a dynamically steered ultrasound focal spot, and that the displacement increases with acoustic power. MR-ARFI has potential for elastography or to guide ultrasound therapies that use low-power pulsed ultrasound exposures, such as drug delivery.     

Development of tissue-simulating optical phantoms: poly-N-isopropylacrylamide solution entrapped inside a hydrogel

The average turbid optical properties of the N-isopropylacrylamide (NIPA) polymer solution entrapped inside a polyacrylamide hydrogel (called an NIPA/PAAM gel system) were studied using a multiwavelength oblique-incidence reflectometer. The turbidity of such a system can be drastically changed by simply switching the temperature from below the low critical solution temperature of the NIPA, around 33 degrees C, to above. The absorption coefficient and the reduced scattering coefficient were obtained as a function of wavelength for samples with selected NIPA and blue dextran concentrations. It is found that the scattering of the optical phantom comes from the NIPA polymer chains and the absorption from the blue dextran. The turbid optical properties of an NIPA/PAAM gel system can be tuned to simulate biological tissues at a specific wavelength by varying compositions of NIPA and blue dextran and further modified by controlling the temperature.