|
1.Introduction1.1.Photoacoustic Tomography and Photoacoustic Tomography PhantomsBreast cancer is the second leading cause of cancer-related deaths in American women.1 X-ray mammography, the main screening method, requires harmful radiation and has low sensitivity for cancer detection in dense breasts.2 Ultrasound can improve breast cancer detection in dense breasts3 and also distinguish benign from cancerous masses, but results in unnecessary biopsies due to low specificity.4 A possible remedy for this problem, the emerging imaging modality photoacoustic tomography (PAT), holds strong potential for mammography applications.5 PAT combines optical tissue excitation with acoustic detection to enable imaging of light-absorbing chromophores (melanin, hemoglobin, water, and lipids). Short-duration light pulses are delivered to tissue, where light absorption results in rapid thermoelastic expansion and subsequent emission of ultrasonic pressure waves. These waves can be detected at the tissue surface using acoustic transducers, such as clinical ultrasound probes.6 Beamforming algorithms are then used to backproject detected acoustic waves to their sources (light-absorbing tissue regions). Because sound is attenuated much more weakly than light in tissue, PAT can produce images of optical absorption in deep tissue (up to several centimeters).7 The primary endogenous source of PAT image contrast is hemoglobin found in blood, and several groups have proposed using multiwavelength PAT imaging to generate maps of oxygen saturation.8,9 These images could be used to identify the abnormal vascular structure seen in cancers, resulting in higher specificity for cancer diagnosis than ultrasound alone. While there are currently no clinical PAT systems approved or cleared for use in the United States, several commercial systems, e.g., Seno Medical (San Antonio, Texas),10 iThera Acuity (München, Germany)11, and Canon (Tokyo, Japan),12 and research systems, e.g., University of Twente (Enschede, The Netherlands)13 and Washington University (St. Louis, Missouri),14 have been applied in patients. However, due to the nascent status of PAT, there are no well-established benchtop methods for predicting in vivo PAT system performance. The development of such approaches is essential for improving the understanding of complex tissue interactions and their impact on real-world clinical imaging system performance. Consensus standards and best-practice documents for mature medical imaging modalities (e.g., mammography, magnetic resonance imaging, and ultrasound) identify tissue-simulating phantoms as essential tools for preclinical performance testing.15,16 Standardized phantoms often implement an idealized geometry to enable uniform evaluation of the effect of one type of parameter on PAT image quality while keeping other parameters the same. A wide variety of tissue-mimicking materials (TMMs) have been studied for use as PAT phantoms, including traditional gelatin hydrogel,17,18 poly(vinyl alcohol) cryogel,19,20 and poly(vinyl chloride) plastisol (PVCP).21–24 Most PAT phantoms presented in the literature use a homogeneous background medium to simulate healthy tissue. However, breast tissue is a highly heterogeneous mixture of fat and fibroglandular tissue. The fibroglandular tissue is surrounded by, admixed with, and covered by a layer of fatty tissue.25 These two types of tissue have been reported to have substantially different optical and acoustic properties, and spatial variations of these properties in the breast may cause errors in image reconstruction or changes in image quality. In ultrasound imaging, the heterogeneity in acoustic properties between the fat and fibroglandular tissue results in significant degradation of ultrasound image quality and spatial resolution. Using acoustic detection, PAT also needs to evaluate this effect of acoustic heterogeneity on PAT image quality. Homogeneous phantoms with uniform optical and acoustic properties may be suitable for basic performance testing, but it is likely that heterogeneous phantoms are needed to simulate complex tissue interactions that affect in vivo imaging performance. To realize such phantoms, the spatial distribution of TMMs and their optical and acoustic properties need to mimic those reported for complex breast tissue. 1.2.Optical and Acoustic Properties of Breast TissueDefining appropriate representative optical and acoustic properties of breast tissue is a challenging problem. The main concern is that both optical and acoustic properties vary significantly due to high variations in intrinsic breast tissue composition. Breast density, or the relative content of fatty tissue and fibroglandular tissue, is a major driver of this variation and may depend on many factors, including body mass index, age, and menopausal status.26 In the field of mammography, breast density is used as a coarse estimate of the amount of fibroglandular tissue present in a breast. The Breast Imaging and Reporting Data System (BI-RADS)27 is a well-accepted standardized format for reporting mammography cases. BI-RADS includes a categorization structure for mammographic density with four types: a (almost entirely fat), b (scattered fibroglandular densities), c (heterogeneously dense), and d (extremely dense). The wide variations in subject breast density have thus led to a broad variation in breast property measurements reported in the literature.28–32 The most important optical properties governing PAT system performance are the optical absorption and scattering coefficients, which have been extensively studied in breasts for the past three decades. There is a significant volume of literature available on the optical absorption26,28,33–35 and scattering properties26,34 of in vivo breast tissue measured using diffuse reflectance spectroscopy and diffuse optical tomography. One important limitation is that these measurements are averaged over large tissue volumes, and it is difficult to isolate the optical properties of pure fatty tissue from pure fibroglandular tissue. Therefore, the majority of optical properties reported are not well categorized according to BI-RADS breast tissue types. Some reported values for pure fatty or fibroglandular breast tissue are available from ex vivo measurements, but these may not be completely representative of in vivo tissue due to effects such as dehydration or exsanguination. Several studies have reported variations of in vivo breast optical properties with breast density (relative fat/fibroglandular tissue composition). However, these data show high variability due to other factors such as age and menopausal state, which may contribute to changes in breast density and relative water and lipid content.26,28 Optical property ranges of BI-RADS breast tissue types at 750 nm are listed in Table 1. The values reported using the combination of the middle two types are labeled as b + c. Optical properties for additional wavelengths will be presented in Sec. 3.1. Table 1Reported optical and acoustic properties as functions of BI-RADS breast tissue type.
The critical acoustic properties affecting PAT system performance are speed of sound, frequency-dependent acoustic attenuation, and backscattering coefficient. Speed of sound describes the propagation velocity of acoustic pressure waves in the medium and is a key input parameter to image reconstruction algorithms, such as a delay-and-sum beamformer, where the received ultrasound signals are delayed to account for time-of-flight differences across the acoustic detector surface and summed at each reconstruction point. Both in vitro measurements31,39 on excised breast tissues maintained at 37°C and in vivo measurements32,36,40,41 reported sound speeds of less than in fat and sound speeds more than in fibroglandular tissue over the frequency range of 2 to 13 MHz. Speed of sound values of BI-RADS breast tissue types at 8 MHz are listed in Table 1 because this frequency is close to the central frequency, 7.5 MHz, of our imaging probe and there is one reference at this frequency. Note: the speeds of sound reported36 at 8 MHz in Table 1 are lower compared to the corresponding values at the other frequencies in the figures in Sec. 3.2. However, the trend in speed of sound versus breast tissue types is similar. The spatial variation in breast speed of sound due to the heterogeneous distribution of fatty and fibroglandular tissue can degrade the spatial resolution and accuracy of spatial registration of ultrasound and photoacoustic images when a single uniform speed of sound is assumed in the image reconstruction.42,43 Breast acoustic attenuation has been extensively characterized by several groups, particularly as a function of breast composition. Ex vivo31,44–46 and in vivo measurements38 have demonstrated lower acoustic attenuation in fatty tissue and higher attenuation in fibroglandular tissue over the frequency range of 1 to 10 MHz. Attenuation of BI-RADS a and b + c breast tissue types at 8 MHz is shown in Table 1. As with optical properties, both the speed of sound and acoustic attenuation of the whole breast also vary with menopausal status and age.30 Because many PAT systems presented in the literature are bimodal,47–49 allowing for overlay of coregistered PAT and B-mode ultrasound images, complex PAT tissue phantoms should be suitable for imaging in both domains. Thus, it is important to account for the acoustic backscattering coefficient, which describes the relative amount of energy reflected back toward an acoustic emitter. This mechanism represents the primary source of contrast in B-mode ultrasound images. It is well known that breast fatty and fibroglandular tissues have substantially different backscattering properties, as fibroglandular tissues produce high intensity in B-mode ultrasound images while fatty tissue is relatively hypoechoic. Ex vivo and in vivo studies have demonstrated that breast fibroglandular tissue may possess an order of magnitude higher backscattering coefficient compared to fatty tissue.37,38,44 Backscattering of BI-RADS a and d breast tissue types at 8 MHz is shown in Table 1. 1.3.Study ObjectivesStated succinctly, the overall objective of this research was to advance clinical translation of mammographic PAT through the development of a realistic tissue-simulating phantom and its implementation to elucidate practical issues impacting PAT image quality in the breast. The specific goals of this work included: (1) developing custom TMMs incorporating realistic optical and acoustic properties of breast fat and mixed fatty-fibroglandular tissue, (2) fabricating two-layer PVCP phantoms composed of these TMMs arranged in a tissue-relevant layered geometry with an undulating boundary, and (3) using this innovative breast phantom to assess the effect of tissue properties and reconstruction algorithm parameters on PAT image quality. By achieving these goals, we can take a major step toward filling critical gaps in the literature and facilitate standardization of this innovative technology. 2.Methods2.1.Phantom FabricationWe have previously developed custom PVCP formulations for use in making PAT phantoms with tunable, biologically relevant optical and acoustic properties.24 PVCP is a suspension of PVC resin in liquid plasticizers, which fuses into a gel after heating to high temperatures (). In this study, we utilize the same base materials, including a dispersion-grade PVC resin (Geon 121A, Mexichem Specialty Resins, Inc., Avon Lake, Ohio), benzyl butyl phthalate (BBP) (Sigma-Aldrich, St Louis, Missouri), and di(2-ethylhexyl) adipate (DEHA) (Sigma-Aldrich). For the first aim of the study, we were able to develop two TMM formulations simulating BI-RADS a and c type breast tissues, respectively. The acoustic properties of the two TMMs were kept within desired acoustic property ranges over frequencies from 4 to 9 MHz of the corresponding tissue types, while their optical properties were kept similarly within desired optical properties of BI-RADS b + c type breast tissue over a wavelength range of 650 to 850 nm. These two types of TMM were labeled as a-fat and c-fibrogland.+fat, respectively (Table 2). The speed of sound and acoustic attenuation were tuned by modifying the volume ratio between the two plasticizers (BBP and DEHA) and the mass/mass PVC concentration. Acoustic backscattering was adjusted by adding soda lime glass microspheres (Spheriglass A, Potter Industries LLC, Malvern, Pennsylvania). Microspheres with diameters of 38 to and 63 to were used in the formulations for a-fat and c-fibrogland.+fat tissue, respectively. Optical scattering was induced by adding anatase titanium dioxide (Sigma-Aldrich). Once the optical and acoustic properties fell within the desired range of the tissue type, the tuning process was stopped. There are two limitations to using the combination of BBP, DEHA, and PVC. One is that the current TMM is not able to achieve a speed of sound higher than . However, the two TMM formulations utilized in this study are suitable because their acoustic properties are biologically relevant and present sufficient differences in acoustic properties to affect image quality. The other is that optical absorption of the current TMM cannot be further lowered to match that of fatty tissue. Since we previously showed this TMM to be optically tunable, we could have imparted distinct sets of optical properties to each layer, mimicking different breast tissue composition types. However, to specifically investigate the effect of acoustical heterogeneity on photoacoustic image quality, we tuned the two formulas to have similar, but still biologically relevant, optical properties. Table 2TMM formulations for stimulating breast fat (a-fat) and fatty-fibroglandular (c-fibrogland.+fat) tissue.
PVCP phantoms were fabricated as described previously.24 Briefly, a large stock PVCP solution was prepared by mixing a binary mixture of plasticizers with 1% v/v calcium-zinc heat stabilizer (M-F Manufacturing Co., Fort Worth, Texas), adding PVC resin, then stirring for 30 min followed by degassing for 60 min. A 75-mL volume of PVCP solution was mixed with , sonicated at 40°C for 40 min, and then reintroduced to the original PVCP stock solution and stirred for 5 min. Glass microspheres were added to a 100-mL batch of PVCP solution and stirred for 5 min. This batch volume was then poured into a 200-mL round-bottom flask immersed in a magnetically stirred oil bath maintained at 190°C using a thermocouple. The flask, which contained a stir bar, was evacuated and stirred at for 15 min. The flask was then removed from the oil bath and stirred magnetically at a gradually reducing speed until the flask temperature was . After reaching this temperature, PVCP was poured into aluminum molds. This whole cooling protocol was used to increase PVCP viscosity at time of pouring, which would minimize settling of the glass microspheres. Samples for phantom material characterization were first produced, including 5-mm-thick, 38-mm-diameter disks for speed of sound, optical property measurements (Sec. 2.2), and acoustic attenuation (Sec. 2.3). A 20-mm-thick, 38-mm-diameter disk was produced for backscatter coefficient measurements (Sec. 2.3). 2.2.Optical CharacterizationAs described in Ref. 24, PVCP disks were placed between 1-mm-thick glass slides (refractive index = 1.51), and diffuse transmittance and reflectance measurements were made over 600 to 1000 nm using an integrating sphere spectrophotometer (Lambda 1050, PerkinElmer, Waltham, Massachusetts). NIST-traceable Spectralon standards were used to normalize measurements. Optical absorption coefficients () and reduced scattering coefficients () were calculated using the inverse adding-doubling method,50 which requires a priori knowledge of the anisotropy factor and refractive index of the sample. Using the Lorentz–Lorenz mixture rule for refractive index,51 a-fat and c-fibrogland.+fat TMMs are expected to have a refractive index of 1.485 and 1.528, respectively. The scattering anisotropy factor was assumed to equal 0.7 based on the Mie scattering theory of (using open-source MATLAB® software52) for both TMMs. 2.3.Acoustic CharacterizationThe sound speed of PVCP-based TMM was measured on 5-mm-thick, 38-mm-diameter disks1,24,53 based on through-transmission53 using a pair of identical broadband transducers (V320, Panametrics, Waltham, Massachusetts) with focal length of 3.81 cm, diameter of 1.27 cm, and center frequency of 7.5 MHz. The transducers faced each other in a water bath, with one transmitting, the other receiving, and the PVCP disk under test positioned at their focus. They were connected to a pulser/receiver (Model 5800PR, Panametrics), and the received ultrasound signals were digitized (8 bit, 50 MHz) using an oscilloscope (9310C, Teledyne LeCroy, Chestnut Ridge, New York). The speed of sound in a sample disk, , was calculated as53 where is the speed of sound in water, is the sample thickness, and is the time delay between pulses traveling through a water path with the sample and a water-only path. Four spatial locations in the sample were measured. For each location, 60 acquisitions were averaged before estimating time delay and sound speed. Mean and standard deviation were calculated for the four estimates.Achieved concurrently with the sound speed measurement, the frequency-dependent acoustic attenuation coefficient was derived as1 where is the acoustic power spectrum measured through a water-only path and is the power spectrum measured through the sample, both calculated over 4 to 9 MHz. It was then fitted to the power-law relationship , where and are the fitting parameters.Backscatter coefficients were measured using the same setup, except in reflection mode involving only one transducer.54–56 The central plane of the 20-mm-thick TMM disk under test was aligned with the focal plane of the transducer. Pulse-echo signals were recorded at 30 spatial positions across the face of the disk by laterally translating the TMM disk at 0.5-mm steps. At each position, 20 acquisitions were averaged. The signals were time gated to isolate echoes of 1.5 mm in depth centered at the focal plane, and the mean, , of the corresponding power spectra was calculated. A reference spectrum, , was obtained using the same approach but replacing the TMM disk with a planar reflector made of low-density polyethylene at the transducer focus. This type of material was chosen to avoid signal saturation during the measurements. The normalized power spectrum, , was calculated as56 where compensates for attenuation effects and is calculated as54,55 where is the frequency-dependent attenuation coefficient in for the sample over the gated volume, the gated length is 1.5 mm, and is the distance from the TMM surface close to the transducer to the center of the gated signal. The backscatter coefficients were then calculated after properly compensating for the transducer geometry from the normalized power spectrum as56 where is the aperture area of the transducer with radius , is the pressure reflection coefficient of the planar reflector (in this case, ), is the transducer focal length, is the pressure focusing gain of the transducer, is the wave number, and and are the zeroth- and first-order Bessel functions.56 The coefficient was then fitted to the power-law relationship , where and are the fitting parameters. Mean and standard deviation were calculated for the set of 30 fitted backscatter spectra for each sample.2.4.Homogeneous and Heterogeneous PhantomsA phantom (, 180.6 mL) with homogeneous background (no undulating boundary) was constructed by pouring PVCP (c-fibrogland.+fat formulation) into a mold containing six 0.5-mm-diameter steel wires (A-M Systems, Sequim, Washington) spaced 5 mm apart vertically and horizontally, as shown in the left side of Fig. 1(a). Metallic wires were chosen because they are visible in both ultrasound and PAT images and easily aligned due to their high tensile strength. The diameter of the wire is similar to the size of a brachytherapy seed (metallic implant used for lesion identification and treatment) and in the diameter range of a blood vessel.57,58 Since the phantom was larger than our 100-mL batch volume, the phantom was produced by sequential layer pouring. Prior to pouring the second layer, a heat gun was applied to the surface of the first layer to produce a thin layer of remelted PVCP. As the second layer was poured, this remelting prevented formation of air gaps and poor PVCP bonding, which would cause photoacoustic artifacts at the boundary. Both layers used the c-fibrogland.+fat formulation to mimic dense breast. A heterogeneous phantom with the same layout was created using the same mold, but pouring the a-fat PVCP formulation for the first layer and the c-fibrogland.+fat formulation for the second layer. This phantom also had an undulating boundary between its two layers to mimic the boundary between superficial fatty tissue and deeper fibroglandular tissue in breast. An in vivo ultrasound B-mode image is shown in Fig. 1(b) as an illustration of this tissue structure. To create this feature, an aluminum block with machined concave surfaces [Fig. 1(d)] was inserted into the mold [Fig. 1(c)]. PVCP solution for the first (fatty) layer was poured onto the plate, and after cooling, the plate was removed, resulting in convexities at the surface of the layer onto which the deep (fibroglandular) layer of TMM was poured. The boundary layer consisted of four segments whose radii of curvature were designed to reflect curvatures found by analysis of the fat-glandular tissue boundary in the ultrasound B-mode images of 10 patients acquired for a different research study (an IRB-approved ultrasound elasticity study conducted at George Washington University). All the images included in this analysis present a clear boundary that can be approximately fitted with an arc. Based on the fitting results, we then modeled the boundary between superficial and deep breast layers as a cascade of convexities covering a reasonable range of radii of curvature. In Fig. 1(a), the boundary convexities (from left to right) have radii of curvature of 30, 24, 16, and 16 mm and boundary angles of 58 deg, 66 deg, 81 deg, and 41 deg. Imaging was performed using a previously developed custom PAT system composed of a cart-based tunable near-infrared pulsed laser (Phocus Mobile, Opotek, Inc., Carlsbad, California) and a research-grade ultrasound system (Vantage 128, Verasonics, Inc., Kirkland, Washington). Phantoms were imaged at 750 nm and a radiant exposure of . Acoustic sensing was performed using a 128-channel ultrasound linear array transducer with a 7.5-MHz center frequency, 7.0-MHz bandwidth, and 38.1-mm length (L11-4v, Verasonics). A significant clutter artifact encountered using this approach was caused by photoacoustic generation at the transducer face resulting from high fluence at the phantom surface and optical absorption of transducer surface material. To reduce this effect, the transducer surface was covered with aluminum foil and coupled to the transducer surface with a thin layer of acoustic coupling gel. While the foil can still cause clutter artifacts, the overall image quality is significantly improved. The transducer was then brought into contact with the phantom top surface using a thin layer of water as the couplant. One-way delay-sum beamforming was applied to the received radio-frequency ultrasound signals to create photoacoustic images with dynamic receiving focusing. Data were apodized using Hanning windows, and the number was set at 0.9. The reconstruction speed of sound was varied between 1437 and in the heterogeneous phantom and between 1440 and in the homogeneous phantom. The spatial size of the metal wires with the same diameter at different depths across the phantom was estimated and compared in the reconstructed PAT imaging domain. The axial and lateral spatial sizes of these imaged targets were defined as the full-width half-maximum (FWHM). All PAT images generated in this study were normalized using the same maximum intensity value. B-mode ultrasound images were also acquired using the same PAT system, applying two-way delay-sum beamforming to received data from plane wave transmission. Gray levels for all ultrasound images were normalized to a single maximum intensity value. Target signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were also computed using uncompressed image intensity data to study the effects of the choice of reconstruction speed of sound on these two imaging quality parameters. To calculate these parameters, a rectangular region of interest (ROI) was first centered at each target. A masked subset of this ROI using an FWHM () was selected as a target ROI. The background ROI was the same size and immediately to the right of the target ROI. Target SNR values were calculated as the mean of the target ROI divided by the standard deviation of the background ROI. Target CNR was calculated as the difference between the mean of the target ROI and background ROI over the standard deviation of the background ROI. 2.5.Undulation Effects in Heterogeneous PhantomsTo evaluate the effects of the boundary undulation between superficial (fatty) and deeper (fibroglandular) layers on PAT image quality, a second complex phantom was constructed using the same process as the heterogeneous phantom described in Sec. 2.4. Instead of wire inclusions at different depths, a horizontal array of six polytetrafluoroethylene (PTFE) tubes (Component Supply Co., Fort Meade, Florida) was positioned at a depth of 20 mm, as shown in the right side of Fig. 1. These tubes have an inner diameter of 0.559 mm and a wall thickness of 0.150 mm and were filled with an India ink solution. We used a 1% aqueous solution of India ink (Speedball, Statesville, North Carolina) with an optical absorption of at 750 nm. This value was chosen to simulate venous blood with a hemoglobin concentration of and oxygen saturation of 70%.59 Similarly, as with the target size measurements in Sec. 2.4, the spatial extent of these photoacoustic signals is also estimated using their FWHM in the corresponding direction. There are signals close to the top and bottom of the tubes. Those close to the top wall were used for spatial size and intensity estimation. Axial size was calculated by doubling the half width at half maximum of the rising part of the bell-shaped signal closest to the surface of the phantom since the falling portion of the bell-shaped signal was affected by the signal close to the bottom wall of the tubing. 3.Results and Discussion3.1.Optical PropertiesOptical properties from the literature have been classified according to four BI-RADS categories. Optical properties of b and d breast tissue types are included to delineate the property range of the other two breast tissue types because the change of optical properties follows the change of breast density. The same rule is also applied to acoustic properties reported from the literature. Sandell and Zhu35 summarized the optical properties of fatty breast tissue measured with four in vivo experimental methods. The resulting 95% confidence intervals are shown in Fig. 2. Jacques34 focuses on the optical properties of mixed breast tissue, which can be represented by the fitted curves. The mean and standard deviation of these curves are also shown in Fig. 2 as shaded areas with the mean at the center. The range of values in Jacques’s review paper overlapped with the lower values in Sandell and Zhu’s review paper. The large variation might be due to the variation in a range of biological factors, such as breast density across subjects, and the measurement variation among experimental methods. Optical properties for four tissue types measured in vivo on human subjects are reported by Taroni et al.26 The values shown in Fig. 2 are calculated from the supplementary material of Taroni et al.’s60 paper provided online. The standard deviations were shown only for tissue type c for clarity. The optical properties of our two phantom formulations were designed to be similar to each other, within the optical property range of breast tissue, and also capable of elucidating the effect of heterogeneous acoustic properties on image quality. As shown in Fig. 2(a), the optical absorptions of two formulas are similar to each other except for the wavelength range 900 to 940 nm. These values are slightly higher than the results reported by Taroni et al. over 670- to 920-nm wavelengths, but still in the range of values reported by Sandell and Zhu. The reduced scattering coefficients of both phantom materials are similar and consistent with the values of BI-RADS c type breast tissue reported by Taroni et al. over the wavelength range of 600 to 1000 nm. 3.2.Acoustic PropertiesSpeed-of-sound values reported in the literature31,32,36,39–41 are categorized into four BI-RADS tissue types (a, b, c, and d) according to the description of tissue types in the literature (fat, fatty-fibroglandular, less fatty-fibroglandular, and fibroglandular tissue in order), as shown in Fig. 3. These values are reported as either or as a range of values. There are large variations in speed of sound for each tissue type. The general trend of increasing speed of sound is reported with increasing amounts of fibroglandular tissue. Our formulations are designed to achieve a speed of sound of for a-fat and for c-fibrogland.+fat. These values fall in the sound speed range of tissue types of a and c and are consistent with the references. Acoustic attenuation values reported in the literature31,38,44,45 are shown in Fig. 4. Shaded areas represent either the range of values or standard deviation around the plotted mean values. The trend of increasing sound speed with increasing glandular tissue content is also observed in acoustic attenuation values, with lower values in fat and higher ones in fatty-fibroglandular tissue. The acoustic attenuation of the phantom formulations is consistent with those in the references and is representative of tissue types a and c. Acoustic backscatter coefficients reported in the literature37,38,44 are shown in Fig. 5 for tissue types a and d. Tissue type d has significantly higher backscatter coefficients than type a. The a-fat formulation gives a backscatter coefficient consistent with that of type a breast tissue reported. The c-fibrogland.+fat formulation gives backscatter coefficients slightly lower than the reported values for d, which is consistent with the dependence of acoustic properties on the proportion of fibroglandular tissue present. 3.3.Homogeneous and Heterogeneous PhantomsThe TMM formulations we have developed enable the construction of homogenous and heterogeneous phantoms. Figures 6(a) and 6(b) are the reconstructed ultrasound and PAT images of a heterogeneous phantom with six wires using the speed of sound of the top layer, . Figures 6(c) and 6(d) are the reconstructed ultrasound and PAT images of a homogeneous phantom with six wires using a speed of sound of . The heterogeneous phantom can be used to evaluate the image quality of the dual-modality photoacoustic/ultrasound system. The ultrasound images shown in Fig. 6(a) provide valuable anatomical information mimicking the clinical imaging environment. This is due to a realistic difference in backscatter coefficient for the two layers. While we simulated one specific breast heterogeneity case with a specific set of layer properties, this TMM and phantom method could be readily tuned/adapted to simulate different tissue composition ratios and geometries. 3.4.Effects of Reconstruction Speed of Sound on Photoacoustic Tomography Image QualityThe heterogeneous phantom described above provides a realistic imaging environment to evaluate the effects of boundary undulation and reconstruction speed of sound on PAT target size. Our conventional delay-and-sum reconstruction algorithm assumed a homogeneous speed of sound, and the mismatch between this assumed speed of sound value and the tissue’s true sound speed will lead to overestimation in PAT target size due to incorrect signal alignment. Using a homogeneous phantom, Fig. 7(a) demonstrates that axial FWHM of targets with the same diameter at different depths is not sensitive to the reconstruction speed of sound. The axial FWHM of the first three targets is smaller than the wire diameter; because of the high optical absorption coefficient of the metal target, target light penetration is shallow and thus photoacoustic signals are generated only near the target surface.57 Figure 7(b) shows that lateral FWHM is sensitive to reconstruction sound speed. This is because incorrect speed of sound leads to residual phase modulation across presummation signals even after beamforming delays are applied, and such a phase term increases lateral FWHM. Reconstruction sound speeds higher than the true values result in greater degradation than do sound speeds lower than the true value. For example, underestimating sound speed by (, green curve) results in approximately one-third of the lateral resolution degradation compared with overestimating sound speed by (, gold curve) for depths between 10 and 25 mm. At all speed-of-sound values, both axial and lateral FWHMs increase as the depth increases. This is mainly due to the TMM’s nonlinear acoustic attenuation frequency spectrum,61,62 which causes a decrease in the spectral bandwidth as shown in Fig. 7(c) and a reduction of central frequency as shown in Fig. 7(d). The heterogeneous phantom contains a top layer with a sound speed of and a bottom layer with a sound speed of . The choice of reconstruction speed of sound affects the estimation of spatial size of the imaging targets in the two layers differently. As in the homogeneous phantom, axial resolution is similar across all reconstruction speeds in the heterogeneous phantom [Fig. 8(a)]. Lateral resolution is more severely degraded along the depth for the first three targets (located in the top layer) than for the fifth and sixth targets (in the bottom layer) [Fig. 8(b)]. Two factors contribute to this phenomenon. First, applying overestimated reconstruction sound speed in the top layer can cause more degradation in lateral FWHM calculation of imaging targets than applying underestimated reconstruction sound speed in the bottom layer. Second, the underestimated sound path length in the deeper layer can be partially compensated by the overestimated path length in the superficial layer when the reconstruction speed falls between the true speeds of the two layers. The lateral FWHM of the fourth target located in the deeper layer shows deviation from an expected monotonic trend and is typically larger than the deeper fifth and sixth targets; this may be due to undulating boundary effects, which will be discussed in Sec. 3.5. Target SNR values are shown in Fig. 9(a) for a homogeneous phantom over a reconstruction speed of sound from 1450 to . As expected, target SNR decreases with depth. The highest SNR is achieved using the true speed of sound at . For heterogeneous phantom images reconstructed over a reconstruction speed of sound from 1437 to , shown in Fig. 9(b), target SNR decreases with the increase of the depth, except for the first target, which may have been due to strong ringing effects for the first target. Setting the reconstruction speed of sound close to that of the top layer provides higher target SNR, except for the first target. Target CNRs have similar plots and are not shown here. As an example, target CNRs at in the homogeneous phantom are provided here. The values are 46.0, 45.3, 43.5, 36.8, 28.8 in dB and decrease as the depth of the targets increases. Mismatch between reconstruction sound speed and true speed of sound also causes incorrect estimation of target depth. Figure 10(a) shows that reconstruction sound speeds higher than the true sound speed cause overestimation of target depth as expected (because reconstruction depth is the product of reconstruction sound speed and time delay), whereas values lower than the true sound speed cause underestimation of target depth in a homogenous phantom. The error in reconstruction depth is directly proportional to the depth. In this analysis, we considered target depths reconstructed using the true speed-of-sound value as ground truth, and the depth of the first target was considered as a reference point. Similar depth overestimation can also be seen for the first three targets in the top layer of the heterogeneous phantom, as shown in Fig. 10(b). Underestimated depths for the second target may be due to slight misalignment during the construction process. The estimation error might be partially offset by the accumulative error from the top layer. The estimation errors for the deepest targets (fourth and fifth) are due to a combination of three factors: reconstruction speed mismatch, compensation from the overestimated path length in the superficial layer, and undulating boundary effects. 3.5.Boundary Undulation Effects on Photoacoustic Tomography Image QualityFigure 11 shows ultrasound (a) and PAT (b) images of a heterogeneous phantom containing six PTFE tubes at the same depth as the fourth wire in the above heterogeneous phantom. Assuming that a two-layer phantom with a flat boundary will result in equal size and intensity of each tube, a phantom containing an undulating boundary can illustrate the changes of the dimensions and intensity of the imaged targets as a function of lateral position relative to the undulation pattern. Figure 12(a) indicates that, as expected, axial FWHM values are not significantly different from each other for all reconstruction sound speed values. Figure 12(b) shows that lateral FWHM values of the targets at the same depth depend on the reconstruction speed of sound. When reconstruction speed of sound is between 1437 and , the lateral FWHM values of the middle four targets are similar across the transducer surface. When the reconstruction speed of sound is between 1460 and , the lateral FWHM values are related to the locations of the targets relative to the boundary undulation pattern. The first and sixth targets from left to right have the worst lateral FWHM. This is mainly because (1) reconstruction is performed with dynamic subaperture, (2) only half of the normal subaperture can be used for the targets close to the edge, and (3) lateral FWHM near the focal plane is inversely proportional to subaperture width, which is 23 mm for a focal ratio ( number) of 0.9 at the imaging depth of about 21 mm. Across the different speeds of sound, the second and fifth targets show similar dependences of lateral FWHM values on reconstruction sound speed. This might be due to the similar positions of the two targets relative to the boundary undulations. The fourth target FWHM is more sensitive to the reconstruction sound speed. This might be due to this target being directly under one of the undulation valleys, which has a focusing effect and caused the worst signal alignment in image reconstruction due to the refraction of sound across this valley. According to Snell’s law, when the incident angle changes from 0 deg to 60 deg in the deeper layer with a speed of sound of , the corresponding refraction angle in the top layer with a speed of sound of changes from 0 deg to 56 deg. In other words, a beam can be bent inward by 0 deg to 4 deg for these incident angles. Such bending was ignored in reconstruction of the PAT image, which was based on delay-and-sum and assumed constant speed of sound and wave propagation along a straight line, causing inaccurate delay compensation before beam summation across the probe surface and a subsequent increase of the lateral FWHM of a target. To further illustrate the effects of sound speed mismatch, we simulated a region (, ) in Fig. 11(a) centered around the fourth target and calculated the beamformer delay error profiles across 64 ultrasound probe elements for this target. For reconstruction sound speed between 1437 to , the low delay error shown in Fig. 12(c) explains the small lateral FWHMs shown in Fig. 12(b), while the largest delay error from using correctly predicts the largest lateral FWHM of this target. In addition to sound speed variation, the delay error also depends on geometry. Figure 12(d) shows beamformer delay error profiles for the same simulated region with the convex boundary replaced by a planar boundary positioned at 12 mm below the probe surface. The profiles are different from those in Fig. 12(c) as well as those (not shown) generated from the same convex boundary at the same location as the planar boundary. In summary, the speed difference across a boundary, the shape of the boundary, the location of the boundary relative to the ultrasound probe, and the light absorber all play a role in signal misalignments in image reconstruction. This finding also suggests that the undulation also contributes to the large variation in lateral FWHM estimation for the fourth wire at the same depth in the first heterogeneous phantom, as discussed in Sec. 3.4. Figure 12(e) shows that the peak intensity of the six target signals at the same depth is also sensitive to both the choice of reconstruction speed of sound and acoustic attenuation distribution across the undulating boundary. For the former, for example, the peak intensities of the fourth target decreased with the increase of reconstruction speed of sound. This trend is the opposite of that of the lateral FWHM values versus reconstruction speed of sound. This opposite relationship may signify that the decrease in peak intensity is due to averaging effects in PAT imaging resulting from the corresponding increased target’s lateral size due to mismatched reconstruction sound speed. The second dependence can be observed crossing the middle four targets at one fixed reconstruction sound speed, such as . Even though the lateral and axial FWHMs of the fourth and fifth targets were similar to each other at this speed, the photoacoustic signal intensity of the fifth target reduced to about 83% of that of the fourth target. This decrease may be explained by the relatively larger portion of highly attenuating fibroglandular TMM in the signal propagation path for the fifth target. For example, averaging 1-mm extra fibroglandular TMM thickness can cause about a 12% intensity decrease using the acoustic attenuation difference between two PVCP formulations at the central frequency of 7.5 MHz. The actual variation of PAT intensity due to acoustic attenuation can be more complex due to the acoustic heterogeneity caused by the undulating boundary. This complexity might also explain the different PAT intensities between the second and fifth targets even though they have similar geometries. As expected, the target SNR shown in Fig. 12(f) has a similar trend as the peak intensity for the middle four targets. 4.ConclusionIn this study, we have developed two customizable PVCP formulations simulating breast fat and mixed fatty-fibroglandular tissue. Using these formulations, we have designed and constructed two-layer heterogeneous phantoms with similar tissue-relevant optical properties and proper acoustic heterogeneity between two layers that provide an imaging environment that is similar to the structure of the human breast for evaluating photoacoustic imaging system performance. With these phantoms, we have studied reconstruction errors using an assumed homogeneous sound speed, as well as the effects of tissue boundary undulations on PAT target size and PAT signal intensity. Our findings indicate that phantoms incorporating realistic properties and geometric features can provide useful tools for generating quantitative insights into the relationship between device parameters, tissue properties, and image quality. DisclosuresThe mention of commercial products, their sources, or their use in connection with material reported herein is not to be construed as either an actual or implied endorsement of such products by the Department of Health and Human Services. No conflicts of interest, financial or otherwise, are declared by the authors. AcknowledgmentsThe authors gratefully acknowledge funding support from the FDA Medical Countermeasures Initiative, the FDA Office of Women’s Health, and the FDA Critical Path Initiative, as well as the ORISE fellowship program through Oak Ridge Associated Universities. The authors also wish to thank Rachel Brem of George Washington University for patient data acquisition, Dr. S. Kaisar Alam for helpful discussions, and Randolph L. Bidinger of FDA for fabricating phantom molds. ReferencesK. A. Wear,
“Cancellous bone analysis with modified least squares Prony’s method and chirp filter: phantom experiments and simulation,”
J. Acoust. Soc. Am., 128
(4), 2191
–2203
(2010). http://dx.doi.org/10.1121/1.3478779 JASMAN 0001-4966 Google Scholar
D. Euhus, P. A. Di Carlo and N. F. Khouri,
“Breast cancer screening,”
Surg. Clin., 95
(5), 991
–1011
(2015). http://dx.doi.org/10.1016/j.suc.2015.05.008 Google Scholar
C. M. Rumack, S. R. Wilson and J. W. Charboneau, Diagnostic Ultrasound, 541
–564 Mosby, Philadelphia, Pennsylvania
(2001). Google Scholar
H. Zhao et al.,
“Limitations of mammography in the diagnosis of breast diseases compared with ultrasonography: a single-center retrospective analysis of 274 cases,”
Eur. J. Med. Res., 20
(1), 49
(2015). http://dx.doi.org/10.1186/s40001-015-0140-6 EJMRDJ 0379-8399 Google Scholar
M. Heijblom et al.,
“Photoacoustic image patterns of breast carcinoma and comparisons with magnetic resonance imaging and vascular stained histopathology,”
Sci. Rep., 5
(1), 11778
(2015). http://dx.doi.org/10.1038/srep11778 Google Scholar
M. Pramanik et al.,
“Design and evaluation of a novel breast cancer detection system combining both thermoacoustic (TA) and photoacoustic (PA) tomography,”
Med. Phys., 35
(6), 2218
–2223
(2008). http://dx.doi.org/10.1118/1.2911157 MPHYA6 0094-2405 Google Scholar
M. H. Xu and L. H. V. Wang,
“Photoacoustic imaging in biomedicine,”
Rev. Sci. Instrum., 77
(4), 041101
(2006). http://dx.doi.org/10.1063/1.2195024 RSINAK 0034-6748 Google Scholar
X. Q. Li et al.,
“High resolution functional photoacoustic tomography of breast cancer,”
Med. Phys., 42
(9), 5321
–5328
(2015). http://dx.doi.org/10.1118/1.4928598 MPHYA6 0094-2405 Google Scholar
K. E. Wilson et al.,
“Multiparametric spectroscopic photoacoustic imaging of breast cancer development in a transgenic mouse model,”
Theranostics, 4
(11), 1062
–1071
(2014). http://dx.doi.org/10.7150/thno.9922 Google Scholar
J. Zalev et al.,
“Opto-acoustic image fusion technology for diagnostic breast imaging in a feasibility study,”
Proc. SPIE, 9419 941909
(2015). http://dx.doi.org/10.1117/12.2081679 PSISDG 0277-786X Google Scholar
A. Taruttis et al.,
“Optoacoustic imaging of human vasculature: feasibility by using a handheld probe,”
Radiology, 281
(1), 256
–263
(2016). http://dx.doi.org/10.1148/radiol.2016152160 RADLAX 0033-8419 Google Scholar
E. Fakhrejahani et al.,
“Clinical report on the first prototype of a photoacoustic tomography system with dual illumination for breast cancer imaging,”
PLoS One, 10
(10), e0139113
(2015). http://dx.doi.org/10.1371/journal.pone.0139113 POLNCL 1932-6203 Google Scholar
A. Garcia-Uribe et al.,
“Dual-modality photoacoustic and ultrasound imaging system for noninvasive sentinel lymph node detection in patients with breast cancer,”
Sci. Rep., 5 15748
(2015). http://dx.doi.org/10.1038/srep15748 SRCEC3 2045-2322 Google Scholar
M. Heijblom, W. Steenbergen and S. Manohar,
“Clinical photoacoustic breast imaging: the Twente experience,”
IEEE Pulse, 6
(3), 42
–46
(2015). http://dx.doi.org/10.1109/MPUL.2015.2409102 Google Scholar
T. J. Pfefer and A. Agrawal,
“A review of consensus test methods for established medical imaging modalities and their implications for optical coherence tomography,”
Proc. SPIE, 8215 82150D
(2012). http://dx.doi.org/10.1117/12.912371 PSISDG 0277-786X Google Scholar
J. P. B. O’Connor et al.,
“Imaging biomarker roadmap for cancer studies,”
Nat. Rev. Clin. Oncol., 14
(3), 169
–186
(2016). http://dx.doi.org/10.1038/nrclinonc.2016.162 Google Scholar
J. R. Cook, R. R. Bouchard and S. Y. Emelianov,
“Tissue-mimicking phantoms for photoacoustic and ultrasonic imaging,”
Biomed. Opt. Express, 2
(11), 3193
–3206
(2011). http://dx.doi.org/10.1364/BOE.2.003193 BOEICL 2156-7085 Google Scholar
P. Ephrat et al.,
“Localization of spherical lesions in tumor-mimicking phantoms by 3D sparse array photoacoustic imaging,”
Med. Phys., 37
(4), 1619
–1628
(2010). http://dx.doi.org/10.1118/1.3352785 MPHYA6 0094-2405 Google Scholar
S. Manohar et al.,
“Photoacoustic mammography laboratory prototype: imaging of breast tissue phantoms,”
J. Biomed. Opt., 9
(6), 1172
–1181
(2004). http://dx.doi.org/10.1117/1.1803548 JBOPFO 1083-3668 Google Scholar
A. Kharine et al.,
“Poly(vinyl alcohol) gels for use as tissue phantoms in photoacoustic mammography,”
Phys. Med. Biol., 48
(3), 357
–370
(2003). http://dx.doi.org/10.1088/0031-9155/48/3/306 PHMBA7 0031-9155 Google Scholar
M. Fonseca et al.,
“Characterisation of a PVCP based tissue-mimicking phantom for quantitative photoacoustic imaging,”
Proc. SPIE, 9539 953911
(2015). http://dx.doi.org/10.1117/12.2183684 PSISDG 0277-786X Google Scholar
S. E. Bohndiek et al.,
“Development and application of stable phantoms for the evaluation of photoacoustic imaging instruments,”
PLoS One, 8
(9), e75533
(2013). http://dx.doi.org/10.1371/journal.pone.0075533 POLNCL 1932-6203 Google Scholar
G. M. Spirou et al.,
“Optical and acoustic properties at 1064 nm of polyvinyl chloride-plastisol for use as a tissue phantom in biomedical optoacoustics,”
Phys. Med. Biol., 50
(14), N141
–N153
(2005). http://dx.doi.org/10.1088/0031-9155/50/14/N01 PHMBA7 0031-9155 Google Scholar
W. C. Vogt et al.,
“Biologically relevant photoacoustic imaging phantoms with tunable optical and acoustic properties,”
J. Biomed. Opt., 21
(10), 101405
(2016). http://dx.doi.org/10.1117/1.JBO.21.10.101405 JBOPFO 1083-3668 Google Scholar
K. L. Bontrager and J. P. Lampignano, Textbook of Radiographic Positioning and Related Anatomy, 8th ed.Mosby, St. Louis, Missouri
(2013). Google Scholar
P. Taroni et al.,
“Breast tissue composition and its dependence on demographic risk factors for breast cancer: non-invasive assessment by time domain diffuse optical spectroscopy,”
PLoS One, 10
(6), e0128941
(2015). http://dx.doi.org/10.1371/journal.pone.0128941 POLNCL 1932-6203 Google Scholar
E. A. Sickles et al.,
“ACR BI-RADS® mammography,”
ACR BI-RADS® Atlas, Breast Imaging Reporting and Data System, American College of Radiology, Reston, Virginia
(2013). http://www.acr.org/quality-safety/resources/birads/mammography Google Scholar
N. Shah et al.,
“Spatial variations in optical and physiological properties of healthy breast tissue,”
J. Biomed. Opt., 9
(3), 534
–540
(2004). http://dx.doi.org/10.1117/1.1695560 JBOPFO 1083-3668 Google Scholar
L. Spinelli et al.,
“Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography,”
J. Biomed. Opt., 9
(6), 1137
–1142
(2004). http://dx.doi.org/10.1117/1.1803546 JBOPFO 1083-3668 Google Scholar
I. Katz-Hanani et al.,
“Age-related ultrasonic properties of breast tissue in vivo,”
Ultrasound Med. Biol., 40
(9), 2265
–2271
(2014). http://dx.doi.org/10.1016/j.ultrasmedbio.2014.03.034 USMBA3 0301-5629 Google Scholar
J. C. Bamber,
“Ultrasonic propagation properties of the breast,”
Ultrasonic Examination of the Breast, John Wiley & Sons Ltd., Chichester, United Kingdom
(1983). Google Scholar
G. Kossoff, E. K. Fry and J. Jellins,
“Average velocity of ultrasound in the human female breast,”
J. Acoust. Soc. Am., 53
(6), 1730
–1736
(1973). http://dx.doi.org/10.1121/1.1913526 Google Scholar
P. Taroni et al.,
“Absorption of collagen: effects on the estimate of breast composition and related diagnostic implications,”
J. Biomed. Opt., 12
(1), 014021
(2007). http://dx.doi.org/10.1117/1.2699170 JBOPFO 1083-3668 Google Scholar
S. L. Jacques,
“Optical properties of biological tissues: a review,”
Phys. Med. Biol., 58
(11), R37
–R61
(2013). http://dx.doi.org/10.1088/0031-9155/58/11/R37 PHMBA7 0031-9155 Google Scholar
J. L. Sandell and T. C. Zhu,
“A review of in-vivo optical properties of human tissues and its impact on PDT,”
J. Biophotonics, 4
(11–12), 773
–787
(2011). http://dx.doi.org/10.1002/jbio.v4.11/12 Google Scholar
J. F. Greenleaf and R. C. Bahn,
“Clinical imaging with transmissive ultrasonic computerized tomography,”
IEEE Trans. Biomed. Eng., BME-28
(2), 177
–185
(1981). http://dx.doi.org/10.1109/TBME.1981.324789 IEBEAX 0018-9294 Google Scholar
M. E. Anderson, M. S. Soo and G. E. Trahey,
“In vivo breast tissue backscatter measurements with 7.5-and 10-MHz transducers,”
Ultrasound Med. Biol., 27
(1), 75
–81
(2001). http://dx.doi.org/10.1016/S0301-5629(00)00310-0 USMBA3 0301-5629 Google Scholar
H. G. Nasief et al.,
“Acoustic properties of breast fat,”
J. Ultrasound Med., 34
(11), 2007
–2016
(2015). http://dx.doi.org/10.7863/ultra.14.07039 JUMEDA 0278-4297 Google Scholar
F. S. Foster, M. Strban and G. Austin,
“The ultrasound macroscope: initial studies of breast tissue,”
Ultrason. Imaging, 6
(3), 243
–261
(1984). http://dx.doi.org/10.1177/016173468400600301 ULIMD4 0161-7346 Google Scholar
P. L. Carson et al.,
“Breast imaging in coronal planes with simultaneous pulse echo and transmission ultrasound,”
Science, 214
(4525), 1141
–1143
(1981). http://dx.doi.org/10.1126/science.7302585 SCIEAS 0036-8075 Google Scholar
G. H. Glover,
“Characterization of in vivo breast tissue by ultrasonic time-of-flight computed tomography,”
Ultrasonic Tissue Characterization II, 221
–225 Washington, D.C.(1979). Google Scholar
J. Jose et al.,
“Speed-of-sound compensated photoacoustic tomography for accurate imaging,”
Med. Phys., 39
(12), 7262
–7271
(2012). http://dx.doi.org/10.1118/1.4764911 MPHYA6 0094-2405 Google Scholar
M. Ye et al.,
“Adaptive speed of sound correction with photoacoustic tomography for imaging quality optimization,”
Chin. Opt. Lett., 14
(8), 081701
(2016). http://dx.doi.org/10.3788/COL CJOEE3 1671-7694 Google Scholar
F. T. D’Astous and F. S. Foster,
“Frequency dependence of ultrasound attenuation and backscatter in breast tissue,”
Ultrasound Med. Biol., 12
(10), 795
–808
(1986). http://dx.doi.org/10.1016/0301-5629(86)90077-3 USMBA3 0301-5629 Google Scholar
F. S. Foster and J. W. Hunt,
“Transmission of ultrasound beams through human tissue-focusing and attenuation studies,”
Ultrasound Med. Biol., 5
(3), 257
–268
(1979). http://dx.doi.org/10.1016/0301-5629(79)90017-6 Google Scholar
P. D. Edmonds et al.,
“Ultrasound tissue characterization of breast biopsy specimens,”
Ultrason. Imaging, 13
(2), 162
–185
(1991). http://dx.doi.org/10.1177/016173469101300204 ULIMD4 0161-7346 Google Scholar
T. N. Erpelding et al.,
“Sentinel lymph nodes in the rat: noninvasive photoacoustic and US imaging with a clinical US system,”
Radiology, 256
(1), 102
–110
(2010). http://dx.doi.org/10.1148/radiol.10091772 RADLAX 0033-8419 Google Scholar
H. S. Salehi et al.,
“Coregistered photoacoustic and ultrasound imaging and classification of ovarian cancer: ex vivo and in vivo studies,”
J. Biomed. Opt., 21
(4), 046006
(2016). http://dx.doi.org/10.1117/1.JBO.21.4.046006 JBOPFO 1083-3668 Google Scholar
Y. T. Wang et al.,
“Toward in vivo biopsy of melanoma based on photoacoustic and ultrasound dual imaging with an integrated detector,”
Biomed. Opt. Express, 7
(2), 279
–286
(2016). http://dx.doi.org/10.1364/BOE.7.000279 BOEICL 2156-7085 Google Scholar
S. A. Prahl, M. J. C. van Gemert and A. J. Welch,
“Determining the optical properties of turbid media by using the adding-doubling method,”
Appl. Opt., 32
(4), 559
–568
(1993). http://dx.doi.org/10.1364/AO.32.000559 APOPAI 0003-6935 Google Scholar
W. Heller,
“Remarks on refractive index mixture rules,”
J. Phys. Chem., 69
(4), 1123
–1129
(1965). http://dx.doi.org/10.1021/j100888a006 JPCHAX 0022-3654 Google Scholar
S. A. Prahl,
“Mie scattering,”
(2016) http://omlc.org/software/mie/ September ). 2016). Google Scholar
K. A. Wear,
“Measurements of phase velocity and group velocity in human calcaneus,”
Ultrasound Med. Biol., 26
(4), 641
–646
(2000). http://dx.doi.org/10.1016/S0301-5629(99)00172-6 USMBA3 0301-5629 Google Scholar
M. O’Donnell and J. G. Miller,
“Quantitative broadband ultrasonic backscatter: an approach to nondestructive evaluation in acoustically inhomogeneous materials,”
J. Appl. Phys., 52
(2), 1056
–1065
(1981). http://dx.doi.org/10.1063/1.328803 JAPIAU 0021-8979 Google Scholar
M. L. Oelze and W. D. O’Brien,
“Frequency-dependent attenuation-compensation functions for ultrasonic signals backscattered from random media,”
J. Acoust. Soc. Am., 111
(5 Pt. 1), 2308
–2319
(2002). http://dx.doi.org/10.1121/1.1452743 JASMAN 0001-4966 Google Scholar
R. J. Lavarello, G. Ghoshal and M. L. Oelze,
“On the estimation of backscatter coefficients using single-element focused transducers,”
J. Acoust. Soc. Am., 129
(5), 2903
–2911
(2011). http://dx.doi.org/10.1121/1.3557036 JASMAN 0001-4966 Google Scholar
T. Mitcham et al.,
“Modulation of photoacoustic signal generation from metallic surfaces,”
J. Biomed. Opt., 18
(5), 056008
(2013). http://dx.doi.org/10.1117/1.JBO.18.5.056008 JBOPFO 1083-3668 Google Scholar
M. A. L. Bell et al.,
“In vivo visualization of prostate brachytherapy seeds with photoacoustic imaging,”
J. Biomed. Opt., 19
(12), 126011
(2014). http://dx.doi.org/10.1117/1.JBO.19.12.126011 JBOPFO 1083-3668 Google Scholar
S. A. Prahl,
“Tabulated molar extinction coefficient for hemoglobin in water,”
(2016) http://omlc.org/spectra/hemoglobin/summary.html September ). 2016). Google Scholar
P. Taroni et al.,
“Breast tissue composition and its dependence on demographic risk factors for breast cancer: non-invasive assessment by time domain diffuse optical spectroscopy,”
PLoS One, 10
(6), e0128941
(2015). http://dx.doi.org/10.1371/journal.pone.0128941 POLNCL 1932-6203 Google Scholar
P. A. Narayana and J. Ophir,
“A closed form method for the measurement of attenuation in nonlinearly dispersive media,”
Ultrason. Imaging, 5
(1), 17
–21
(1983). http://dx.doi.org/10.1177/016173468300500102 ULIMD4 0161-7346 Google Scholar
J. Ophir and P. Jaeger,
“Spectral shifts of ultrasonic propagation through media with non-linear dispersive attenuation,”
Ultrason. Imaging, 4
(3), 282
–289
(1982). http://dx.doi.org/10.1177/016173468200400304 ULIMD4 0161-7346 Google Scholar
BiographyCongxian Jia received her BS and MS degrees in mechanical engineering from Beijing University, China, in 1999 and 2002, respectively. She also received her MS degree in biomechanics of aerospace and mechanical engineering from Boston University in 2004 and her PhD in biomedical engineering from the University of Michigan in 2010. Currently, she works as a research fellow at the FDA, and her research interests include ultrasound elasticity imaging and photoacoustic imaging. William C. Vogt received his BS degree in mechanical engineering from the University of Massachusetts Amherst in 2009 and his PhD in biomedical engineering from Virginia Polytechnic Institute and State University in 2013. Since 2013, he has been conducting photoacoustic imaging research as a research fellow at the FDA in the Office of Science and Engineering Laboratories. His research interests include photoacoustic imaging, phantom-based image quality testing, and biophotonic medical device characterization and evaluation. Keith A. Wear received his BA degree in applied physics from the University of California, San Diego, and his MS and PhD degrees in applied physics from Stanford University. He is the FDA Acoustics Laboratory leader. He is an associate editor of Journal of the Acoustical Society of America, Ultrasonic Imaging, and IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control. He is a fellow of the Acoustical Society of America, American Institute for Medical and Biological Engineering, and American Institute of Ultrasound in Medicine. T. Joshua Pfefer received his BS degree in mechanical engineering from Northwestern University, his MS degree in mechanical engineering, and his PhD in biomedical engineering from the University of Texas at Austin and was a research fellow at the Wellman Laboratories of Photomedicine. In 2000, he joined the FDA, where he is currently the leader of the Optical Diagnostic Devices Laboratory. His group’s research focuses on safety and effectiveness in emerging clinical biophotonic spectroscopy and imaging technologies. Brian S. Garra was trained at the University of Washington and the University of Utah and currently practices radiology at the Washington DC Veterans Affairs Medical Center. He also does research and medical device evaluation at the FDA as an associate director-clinical in the Division of Imaging, Diagnostics, and Software Reliability. Currently, he is working on elastography, phantoms for photoacoustic system/elastographic system evaluation, the quantitative imaging biomarker alliance, and ultrasound outreach in Peru. |