X-ray characterisation of normal and neoplastic breast tissues

The correlation between mammographic densities and molecular pathology in breast cancer

This study aimed to analyze the correlation between mammographic density obtained by density analysis software (DAS)/radiologists visual (RV) classification with molecular subtype, and the expression levels of estrogen receptor (ER), progesterone receptor (PR), Ki67 antigen (Ki-67), p53 gene (p53), and human epidermal growth factor receptor-2 (HER2). A total of 688 breast cancer patients with digital mammography and complete molecular pathological results in Tianjin Medical University Cancer Institute and Hospital between February 2015 and February 2016 were collected. The DAS-density grade (DASD) and the radiologists visually classified density grade (RVD) were evaluated by 3 radiologists. The correlation between density grade and the expression levels of ER, PR, Ki-67, p53, HER2 and breast cancer molecular subtype (PMS) were analyzed. The agreement between DASD and RVD was explored. ER, PR and HER-2 positive rate were significantly different among patients with different RVD grades (P< 0.05). HER2 positive rates showed an increasing trend following RVD upgrading (P𝑡𝑟𝑒𝑛𝑑< 0.05). HER-2 positive rate in RVD D1 + D2 was 7.69%, which was higher than that in D3 + D4 (P< 0.05). The ER and Ki-67 expressions in patients were markedly different among DASD (P= 0.009 and 0.002) and RVD (P= 0.012 and 0.036) with different grades. The kappa value of each DASD to RVD was 0.31 (P< 0.01). The RVD 3 proportion was 14.58% (63/432) in HER2 Over-expressing subtype, which was apparently higher than RVD1 (2.43%, 1/41) (P< 0.05). Breast density may be partial correlated with molecular pathology in breast cancer.

Evaluation of human observer performance on lesion detectability in single-slice and multislice dedicated breast cone beam CT images with breast anatomical background

PURPOSE

We evaluate the lesion detectability using human and model observer studies in single-slice and multislice cone beam computed tomography (CBCT) images with a breast anatomical background. The purposes of this work are (a) to compare human observer detectability between single-slice and multislice images for different signal sizes and noise structures, (b) to investigate the effect of different multislice viewing modes (i.e., sequential and simultaneous) on the detectability by a human observer, and (c) to predict the detectability by a human observer in single-slice and multislice images using single-slice channelized Hotelling observer (ssCHO) and multislice CHO (msCHO), respectively.

METHODS

Breast anatomical background is modeled using a power law spectrum of mammograms and the lesion is modeled with a spherical signal. We conduct signal-known-exactly and background-known-statistically detection tasks on transverse and longitudinal images reconstructed using the Feldkamp-Davis-Kress algorithm with Hanning and Ram-Lak weighted ramp filters. The human observer study is conducted on three different viewing modes: single-slice, and sequential and simultaneous multislice. To predict the detectability by a human observer, we use ssCHO and msCHO with anthropomorphic channels (i.e., dense difference-of-Gaussian (D-DOG) and Gabor channels) and internal noise.

RESULTS

The detectability by a human observer increases for multislice images compared to single-slice images. For multislice images, the sequential viewing mode yields higher detectability than the simultaneous viewing mode. However, the relative rank of detectability by a human observer for different signal sizes, image planes, and reconstruction filters is not much different between the viewing modes. Detectability by CHO with internal noise shows good correlation with that of the human observer for all viewing modes.

CONCLUSIONS

Detectability by a human observer in CBCT images with breast anatomical background is affected by the image viewing mode, and the effect of the viewing mode depends on the signal size and noise structure. D-DOG and Gabor CHO with internal noise predict the detectability by a human observer well for both the single-slice and multislice image viewing modes.

Monochromatic mammography using scanning multilayer X-ray mirrors

A prototype system for breast imaging using monochromatic X-rays has been developed using a scanning multilayer X-ray mirror in combination with a conventional mammography tube and an imaging detector. The X-ray mirror produces a monochromatic fan beam tuned near 19 keV, with an energy bandpass of approximately 1.5 keV. Rotating the mirror about the tube's focal spot in synchronization with the X-ray generator and detector enables the acquisition of monochromatic X-ray images over large areas. The X-ray mirror also can be rotated completely out of the beam so that conventional polychromatic images can be acquired using a K-edge filter, facilitating direct comparison between the two modes of operation. The system was used to image synthetic, tissue-equivalent breast phantoms in order to experimentally quantify the improvements in image quality and dose that can be realized using monochromatic radiation. Nine custom phantoms spanning a range of thicknesses and glandular/adipose ratios, each containing both glandular- and calcification-equivalent features, were used to measure contrast and signal-difference-to-noise ratio (SDNR). Mean glandular dose (MGD) was computed from measured entrance exposure, and a figure-of-merit (FOM) was computed as FOM = SDNR²/MGD in each case. Monochromatic MGD ranges from 0.606 to 0.134 of polychromatic MGD for images having comparable glandular SDNR, depending on breast thickness and glandularity; relative monochromatic dose decreases with increasing glandularity for all thicknesses. Monochromatic FOM values are higher than the corresponding polychromatic FOM values in all but one case. Additionally, the monochromatic contrast for glandular features is higher than the polychromatic contrast in all but one case as well. These results represent important steps toward the realization of clinically practical monochromatic X-ray breast imaging systems having lower dose and better image quality, including those for digital mammography, digital breast tomosynthesis, contrast-enhanced spectral mammography and other modalities, for safer, more accurate breast cancer detection, diagnosis and staging.

Toward in vivo dosimetry in external beam radiotherapy using x-ray acoustic computed tomography: A soft-tissue phantom study validation

PURPOSE

To study, using phantoms made from biological tissues, the feasibility and practical challenges of monitoring the position of the radiation beam and the deposited dose by x-ray acoustic computed tomography (XACT) during external beam radiotherapy delivery.

MATERIAL AND METHODS

A prototype XACT system with a single immersion ultrasound transducer, which was positioned around the target sample driven by a motor-controlled rotation stage, was used to acquire the x-ray acoustic (XA) signals produced by a medical linear accelerator (Linac) to form an XACT image of the irradiated phantom. To investigate the feasibility of XACT in tracking the position of radiation dose, a large piece of veal liver with embedded fat tissue was imaged and beam misalignments were measured. Next, we explored the sensitivity of XACT in monitoring and quantifying the delivered dose, in which a block of porcine gel was embedded with equally spaced lard cylinders and imaged. The doses on the lard cylinders modulated by physical wedges were quantified from the XACT image and were verified by comparison to measurements from radiochromic films as the gold standard. Then, to simulate how XACT can perform in a targeted tissue in the human body, a porcine gel phantom with lard cylinders covered by different materials (bone, muscle, and air gap, respectively) was also imaged.

RESULTS

The reconstructed XACT images of the phantoms show congruence with the boundaries of the beam field and the interfaces between the different tissue materials. The beam displacement from the target was tracked properly by the reconstructed XACT images. An intensity difference as small as 2.9% in delivered dose region can be measured from XACT images P = 0.02. The intensities of XACT images were highly correlated to the film measurements with an R2 better than 0.986. The expected variances of dose delivered to different target regions as a result of the difference in attenuation were successfully captured by the XACT images.

CONCLUSIONS

This study validated the feasibility of XACT in accurately obtaining relative dose maps of tissue-mimicking phantoms. XACT offers a practical method for verifying the beam position against the target and quantifying the relative dose delivered to the target during external beam radiotherapy.

Mammographic breast density in recent and longer-standing ethiopian immigrants to israel

High breast density is associated with an increased risk of breast cancer development. Little is known concerning ethnic variations in breast density and its relevant contributing factors. We aimed to study breast density among Ethiopian immigrants to Israel in comparison with Israeli-born women and to determine any effect on breast density of the length of residency in the immigrant population. Mammographic breast density using the BI-RADS system was estimated and compared between 77 women of Ethiopian origin who live in Israel and 177 Israeli-born controls. Logistic regression analysis was performed to estimate the odds ratios (OR) for high density (BI-RADS score ≥ 3) vs low density (BI-RADS score < 3) cases, comparing the 2 origin groups. Ethiopian-born women had a crude OR of 0.15 (95% CI: 0.08-0.26) for high breast density compared with Israeli-born women. Adjustments for various cofounders did not affect the results. Time since immigration to Israel seemed to modify the relationship, with a stronger association for women who immigrated within 2 years prior to mammography (OR:0.07, 95% CI: 0.03-0.17) as opposed to women with a longer residency stay in Israel (OR:0.23, 95% CI:0.10-0.50). Adjustments of various confounders did not alter these findings. Breast density in Ethiopian immigrants to Israel is significantly lower than that of Israeli-born controls. Our study suggests a positive association between time since immigration and breast density. Future studies are required to define the possible effects of dietary change on mammographic density following immigration.

Human and model observer performance for lesion detection in breast cone beam CT images with the FDK reconstruction

We investigate the detectability of breast cone beam computed tomography images using human and model observers and the variations of exponent, β, of the inverse power-law spectrum for various reconstruction filters and interpolation methods in the Feldkamp-Davis-Kress (FDK) reconstruction. Using computer simulation, a breast volume with a 50% volume glandular fraction and a 2mm diameter lesion are generated and projection data are acquired. In the FDK reconstruction, projection data are apodized using one of three reconstruction filters; Hanning, Shepp-Logan, or Ram-Lak, and back-projection is performed with and without Fourier interpolation. We conduct signal-known-exactly and background-known-statistically detection tasks. Detectability is evaluated by human observers and their performance is compared with anthropomorphic model observers (a non-prewhitening observer with eye filter (NPWE) and a channelized Hotelling observer with either Gabor channels or dense difference-of-Gaussian channels). Our results show that the NPWE observer with a peak frequency of 7cyc/degree attains the best correlation with human observers for the various reconstruction filters and interpolation methods. We also discover that breast images with smaller β do not yield higher detectability in the presence of quantum noise.

The origins of breast cancer associated with mammographic density: a testable biological hypothesis

Background

Our purpose is to develop a testable biological hypothesis to explain the known increased risk of breast cancer associated with extensive percent mammographic density (PMD), and to reconcile the apparent paradox that although PMD decreases with increasing age, breast cancer incidence increases.

Methods

We used the Moolgavkar model of carcinogenesis as a framework to examine the known biological properties of the breast tissue components associated with PMD that includes epithelium and stroma, in relation to the development of breast cancer. In this model, normal epithelial cells undergo a mutation to become intermediate cells, which, after further mutation, become malignant cells. A clone of such cells grows to become a tumor. The model also incorporates changes with age in the number of susceptible epithelial cells associated with menarche, parity, and menopause. We used measurements of the radiological properties of breast tissue in 4454 healthy subjects aged from 15 to 80+ years to estimate cumulative exposure to PMD (CBD) in the population, and we examined the association of CBD with the age-incidence curve of breast cancer in the population.

Results

Extensive PMD is associated with a greater number of breast epithelial cells, lobules, and fibroblasts, and greater amounts of collagen and extracellular matrix. The known biological properties of these tissue components may, singly or in combination, promote the acquisition of mutations by breast epithelial cells specified by the Moolgavkar model, and the subsequent growth of a clone of malignant cells to form a tumor. We also show that estimated CBD in the population from ages 15 to 80+ years is closely associated with the age-incidence curve of breast cancer in the population.

Conclusions

These findings are consistent with the hypothesis that the biological properties of the breast tissue components associated with PMD increase the probability of the transition of normal epithelium to malignant cells, and that the accumulation of mutations with CBD may influence the age-incidence curve of breast cancer. This hypothesis gives rise to several testable predictions.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-0941-y) contains supplementary material, which is available to authorized users.

Characterization of CT Hounsfield Units for 3D acquisition trajectories on a dedicated breast CT system

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is important, especially when using multiple sources, novel detector and non-traditional trajectories. Dedicated clinical breast CT (BCT) systems therefore should be similarly evaluated. In this study, uniform cylindrical phantoms filled with various uniform density fluids were used to characterize differences in HU values between simple circular and complex 3D (saddle) orbits. Based on ACR recommendations, the HU accuracy, center-to-edge variability within a slice, and overall variability within the reconstructed volume were characterized for simple and complex acquisitions possible on a single versatile BCT system. Results illustrate the statistically significantly better performance of the saddle orbit, especially close to the chest and nipple regions of what would clinically be a pendant breast volume. The incomplete cone beam acquisition of a simple circular orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values throughout the reconstructed volumes. This study further establishes the advantages of using 3D acquisition trajectories for breast CT as well as other applications by demonstrating the robustness of HU values throughout large reconstructed volumes.

Task-Based Modeling of a 5k Ultra-High-Resolution Medical Imaging System for Digital Breast Tomosynthesis

High-resolution, low-noise X-ray detectors based on CMOS active pixel sensor (APS) technology have demonstrated superior imaging performance for digital breast tomosynthesis (DBT). This paper presents a task-based model for a high-resolution medical imaging system to evaluate its ability to detect simulated microcalcifications and masses as lesions for breast cancer. A 3-D cascaded system analysis for a 50- [Formula: see text] pixel pitch CMOS APS X-ray detector was integrated with an object task function, a medical imaging display model, and the human eye contrast sensitivity function to calculate the detectability index and area under the ROC curve (AUC). It was demonstrated that the display pixel pitch and zoom factor should be optimized to improve the AUC for detecting small microcalcifications. In addition, detector electronic noise of smaller than 300 e^- and a high display maximum luminance (>1000 cd/cm ²) are desirable to distinguish microcalcifications of [Formula: see text] in size. For low contrast mass detection, a medical imaging display with a minimum of 12-bit gray levels is recommended to realize accurate luminance levels. A wide projection angle range of greater than ±30° in combination with the image gray level magnification could improve the mass detectability especially when the anatomical background noise is high. On the other hand, a narrower projection angle range below ±20° can improve the small, high contrast object detection. Due to the low mass contrast and luminance, the ambient luminance should be controlled below 5 cd/ [Formula: see text]. Task-based modeling provides important firsthand imaging performance of the high-resolution CMOS-based medical imaging system that is still at early stage development for DBT. The modeling results could guide the prototype design and clinical studies in the future.

Quantitative contrast-enhanced spectral mammography based on photon-counting detectors: A feasibility study

PURPOSE

To investigate the feasibility of accurate quantification of iodine mass thickness in contrast-enhanced spectral mammography.

MATERIALS AND METHODS

A computer simulation model was developed to evaluate the performance of a photon-counting spectral mammography system in the application of contrast-enhanced spectral mammography. A figure-of-merit (FOM), which was defined as the decomposed iodine signal-to-noise ratio (SNR) with respect to the square root of the mean glandular dose (MGD), was chosen to optimize the imaging parameters, in terms of beam energy, splitting energy, and prefiltrations for breasts of various thicknesses and densities. Experimental phantom studies were also performed using a beam energy of 40 kVp and a splitting energy of 34 keV with 3 mm Al prefiltration. A two-step calibration method was investigated to quantify the iodine mass thickness, and was validated using phantoms composed of a mixture of glandular and adipose materials, for various breast thicknesses and densities. Finally, the traditional dual-energy log-weighted subtraction method was also studied as a comparison. The measured iodine signal from both methods was compared to the known value to characterize the quantification accuracy and precision.

RESULTS

The optimal imaging parameters, which lead to the highest FOM, were found at a beam energy between 42 and 46 kVp with a splitting energy at 34 keV. The optimal tube voltage decreased as the breast thickness or the Al prefiltration increased. The proposed quantification method was able to measure iodine mass thickness on phantoms of various thicknesses and densities with high accuracy. The root-mean-square (RMS) error for cm-scale lesion phantoms was estimated to be 0.20 mg/cm² . The precision of the technique, characterized by the standard deviation of the measurements, was estimated to be 0.18 mg/cm² . The traditional weighted subtraction method also predicted a linear correlation between the measured signal and the known iodine mass thickness. However, the correlation slope and offset values were strongly dependent on the total breast thickness and density.

CONCLUSION

The results of this study suggest that iodine mass thickness for cm-scale lesions can be accurately quantified with contrast-enhanced spectral mammography. The quantitative information can potentially improve the differential power for malignancy.

Accuracy assessment and characterization of x-ray coded aperture coherent scatter spectral imaging for breast cancer classification

Although transmission-based x-ray imaging is the most commonly used imaging approach for breast cancer detection, it exhibits false negative rates higher than 15%. To improve cancer detection accuracy, x-ray coherent scatter computed tomography (CSCT) has been explored to potentially detect cancer with greater consistency. However, the 10-min scan duration of CSCT limits its possible clinical applications. The coded aperture coherent scatter spectral imaging (CACSSI) technique has been shown to reduce scan time through enabling single-angle imaging while providing high detection accuracy. Here, we use Monte Carlo simulations to test analytical optimization studies of the CACSSI technique, specifically for detecting cancer in ex vivo breast samples. An anthropomorphic breast tissue phantom was modeled, a CACSSI imaging system was virtually simulated to image the phantom, a diagnostic voxel classification algorithm was applied to all reconstructed voxels in the phantom, and receiver-operator characteristics analysis of the voxel classification was used to evaluate and characterize the imaging system for a range of parameters that have been optimized in a prior analytical study. The results indicate that CACSSI is able to identify the distribution of cancerous and healthy tissues (i.e., fibroglandular, adipose, or a mix of the two) in tissue samples with a cancerous voxel identification area-under-the-curve of 0.94 through a scan lasting less than 10 s per slice. These results show that coded aperture scatter imaging has the potential to provide scatter images that automatically differentiate cancerous and healthy tissue within ex vivo samples. Furthermore, the results indicate potential CACSSI imaging system configurations for implementation in subsequent imaging development studies.

Investigation of x-ray spectra for iodinated contrast-enhanced dedicated breast CT

Screening for breast cancer with mammography has been very successful, resulting in part to a reduction of breast cancer mortality by approximately 39% since 1990. However, mammography still has limitations in performance, especially for women with dense breast tissue. Iodinated contrast-enhanced, dedicated breast CT (BCT) has been proposed to improve lesion analysis and the accuracy of diagnostic workup for patients suspected of having breast cancer. A mathematical analysis to explore the use of various x-ray filters for iodinated contrast-enhanced BCT is presented. To assess task-based performance, the ideal linear observer signal-to-noise ratio (SNR) is used as a figure-of-merit under the assumptions of a linear, shift-invariant imaging system. To estimate signal and noise propagation through the BCT detector, a parallel-cascade model was used. The lesion model was embedded into a structured background and included a realistic level of iodine uptake. SNR was computed for 84,000 different exposure settings by varying the kV setting, x-ray filter materials and thickness, breast size, and composition and radiation dose. It is shown that some x-ray filter material/thickness combinations can provide up to 75% improvement in the linear ideal observer SNR over a conventionally used x-ray filter for BCT. This improvement in SNR can be traded off for substantial reductions in mean glandular dose.

Mammographic Breast Density and Breast Cancer Risk: Implications of the Breast Density Legislation for Health Care Practitioners

Breast density has emerged as a critical phenotypic marker of increased breast cancer risk. The breast density legislation, passed in multiple states, requires patient notification of the implications of the breast density on breast cancer risk and screening. Supplemental screening may be suggested in the state regulation; however, there are limited data to guide conversations with patients. This article will review the current state of supplemental screening in women with dense breasts and discuss theories of the mechanism of action. Guidance is provided to assist in shared decision making and appropriate patient counseling.

Design and application of a structured phantom for detection performance comparison between breast tomosynthesis and digital mammography

This paper introduces and applies a structured phantom with inserted target objects for the comparison of detection performance of digital breast tomosynthesis (DBT) against 2D full field digital mammography (FFDM). The phantom consists of a 48 mm thick breast-shaped polymethyl methacrylate (PMMA) container filled with water and PMMA spheres of different diameters. Three-dimensionally (3D) printed spiculated masses (diameter range: 3.8-9.7 mm) and non-spiculated masses (1.6-6.2 mm) along with microcalcifications (90-250 µm) were inserted as targets. Reproducibility of the phantom application was studied on a single system using 30 acquisitions. Next, the phantom was evaluated on five different combined FFDM & DBT systems and target detection was compared for FFDM and DBT modes. Ten phantom images in both FFDM and DBT modes were acquired on these 5 systems using automatic exposure control. Five readers evaluated target detectability. Images were read with the four-alternative forced-choice (4-AFC) paradigm, with always one segment including a target and 3 normal background segments. The percentage of correct responses (PC) was assessed based on 10 trials of each reader for each object type, size and imaging modality. Additionally, detection threshold diameters at 62.5 PC were assessed via non-linear regression fitting of the psychometric curve. The reproducibility study showed no significant differences in PC values. Evaluation of target detection in FFDM showed that microcalcification detection thresholds ranged between 110 and 118 µm and were similar compared to the detection in DBT (range of 106-158 µm). In DBT, detection of both mass types increased significantly (p  =  0.0001 and p  =  0.0002 for non-spiculated and spiculated masses respectively) compared to FFDM, achieving almost 100% detection for all spiculated mass diameters. In conclusion, a structured phantom with inserted targets was able to show evidence for detectability differences between FFDM and DBT modes for five commercial systems. This phantom has potential for application in task-based assessment at acceptance and commissioning testing of DBT systems.

Characterization of X-ray scattering for various phantoms and clinical breast geometries using breast CT on a dedicated hybrid system

OBJECTIVE

The purpose of this study was to utilize a dedicated breast CT system using a 2D beam stop array to physically evaluate the scatter to primary ratios (SPRs) of different geometric phantoms and prospectively acquired clinical patient data.

METHODS

Including clinically unrealizable compositions of 100% glandular and 100% fat, projection images were acquired using three geometrically different phantoms filled with fluids simulating breast tissue. The beam stop array method was used for measuring scatter in projection space, and creating the scatter corrected primary images. 2D SPRs were calculated. Additionally, a new figure of merit, the 3D normalized scatter contribution (NSC) volumes were calculated.

RESULTS

The 2D SPR values (0.52-1.10) were primarily dependent on phantom geometry; a secondary dependence was due to their uniform density; 2D SPRs were low frequency and smoothly varying in the uniformly filled phantoms. SPRs of clinical patient data followed similar trends as phantoms, but with noticeable deviations and high frequency components due to the heterogeneous distribution of glandular tissue. The maximum measured patient 2D SPRs were all <0.6, even for the largest diameter breast. These results demonstrate modest scatter components with changing object geometries and densities; the 3D NSC volumes with higher frequency components help visualize scatter distribution throughout the reconstructed image volumes. Furthermore, the SPRs in the heterogeneous clinical breast cases were underestimated by the equivalent density, uniformly filled phantoms.

CONCLUSIONS

These results provide guidance on the use of uniformly distributed density and differently shaped phantoms when considering simulations. They also clearly demonstrate that results from patients can vary considerably from 2D SPRs of uniformly simulated phantoms.

Visual-search observers for assessing tomographic x-ray image quality

PURPOSE

Mathematical model observers commonly used for diagnostic image-quality assessments in x-ray imaging research are generally constrained to relatively simple detection tasks due to their need for statistical prior information. Visual-search (VS) model observers that employ morphological features in sequential search and analysis stages have less need for such information and fewer task constraints. The authors compared four VS observers against human observers and an existing scanning model observer in a pilot study that quantified how mass detection and localization in simulated digital breast tomosynthesis (DBT) can be affected by the number P of acquired projections.

METHODS

Digital breast phantoms with embedded spherical masses provided single-target cases for a localization receiver operating characteristic (LROC) study. DBT projection sets based on an acquisition arc of 60° were generated for values of P between 3 and 51. DBT volumes were reconstructed using filtered backprojection with a constant 3D Butterworth postfilter; extracted 2D slices were used as test images. Three imaging physicists participated as observers. A scanning channelized nonprewhitening (CNPW) observer had knowledge of the mean lesion-absent images. The VS observers computed an initial single-feature search statistic that identified candidate locations as local maxima of either a template matched-filter (MF) image or a gradient-template MF (GMF) image. Search inefficiencies that modified the statistic were also considered. Subsequent VS candidate analyses were carried out with (i) the CNPW statistical discriminant and (ii) the discriminant computed from GMF training images. These location-invariant discriminants did not utilize covariance information. All observers read 36 training images and 108 study images per P value. Performance was scored in terms of area under the LROC curve.

RESULTS

Average human-observer performance was stable for P between 7 and 35. In the absence of search inefficiencies, the VS models based on the GMF analysis provided the best correlation (Pearson ρ ≥ 0.62) with the human results. The CNPW-based VS observers deviated from the humans primarily at lower values of P. In this limited study, search inefficiencies allowed for good quantitative agreement with the humans for most of the VS observers.

CONCLUSIONS

The computationally efficient training requirements for the VS observer are suitable for high-resolution imaging, indicating that the observer framework has the potential to overcome important task limitations of current model observers for x-ray applications.

Dosimetry in x-ray-based breast imaging

The estimation of the mean glandular dose to the breast (MGD) for x-ray based imaging modalities forms an essential part of quality control and is needed for risk estimation and for system design and optimisation. This review considers the development of methods for estimating the MGD for mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (DBCT). Almost all of the methodology used employs Monte Carlo calculated conversion factors to relate the measurable quantity, generally the incident air kerma, to the MGD. After a review of the size and composition of the female breast, the various mathematical models used are discussed, with particular emphasis on models for mammography. These range from simple geometrical shapes, to the more recent complex models based on patient DBCT examinations. The possibility of patient-specific dose estimates is considered as well as special diagnostic views and the effect of breast implants. Calculations using the complex models show that the MGD for mammography is overestimated by about 30% when the simple models are used. The design and uses of breast-simulating test phantoms for measuring incident air kerma are outlined and comparisons made between patient and phantom-based dose estimates. The most widely used national and international dosimetry protocols for mammography are based on different simple geometrical models of the breast, and harmonisation of these protocols using more complex breast models is desirable.

Quantification of breast lesion compositions using low-dose spectral mammography: A feasibility study

PURPOSE

The positive predictive power for malignancy can potentially be improved, if the chemical compositions of suspicious breast lesions can be reliably measured in screening mammography. The purpose of this study is to investigate the feasibility of quantifying breast lesion composition, in terms of water and lipid contents, with spectral mammography.

METHODS

Phantom and tissue samples were imaged with a spectral mammography system based on silicon-strip photon-counting detectors. Dual-energy calibration was performed for material decomposition, using plastic water and adipose-equivalent phantoms as the basis materials. The step wedge calibration phantom consisted of 20 calibration configurations, which ranged from 2 to 8 cm in thickness and from 0% to 100% in plastic water density. A nonlinear rational fitting function was used in dual-energy calibration of the imaging system. Breast lesion phantoms, made from various combinations of plastic water and adipose-equivalent disks, were embedded in a breast mammography phantom with a heterogeneous background pattern. Lesion phantoms with water densities ranging from 0% to 100% were placed at different locations of the heterogeneous background phantom. The water density in the lesion phantoms was measured using dual-energy material decomposition. The thickness and density of the background phantom were varied to test the accuracy of the decomposition technique in different configurations. In addition, an in vitro study was also performed using mixtures of lean and fat bovine tissue of 25%, 50%, and 80% lean weight percentages as the background. Lesions were simulated by using breast lesion phantoms, as well as small bovine tissue samples, composed of carefully weighed lean and fat bovine tissues. The water densities in tissue samples were measured using spectral mammography and compared to measurement using chemical decomposition of the tissue.

RESULTS

The thickness of measured and known water contents was compared for various lesion configurations. There was a good linear correlation between the measured and the known values. The root-mean-square errors in water thickness measurements were 0.3 and 0.2 mm for the plastic phantom and bovine tissue backgrounds, respectively.

CONCLUSIONS

The results indicate that spectral mammography can be used to accurately characterize breast lesion composition in terms of their equivalent water and lipid contents.

Breast Cancer Risk and Mammographic Density Assessed with Semiautomated and Fully Automated Methods and BI-RADS

Purpose To compare three metrics of breast density on full-field digital mammographic (FFDM) images as predictors of future breast cancer risk. Materials and Methods This institutional review board-approved study included 125 women with invasive breast cancer and 274 age- and race-matched control subjects who underwent screening FFDM during 2004-2013 and provided informed consent. The percentage of density and dense area were assessed semiautomatically with software (Cumulus 4.0; University of Toronto, Toronto, Canada), and volumetric percentage of density and dense volume were assessed automatically with software (Volpara; Volpara Solutions, Wellington, New Zealand). Clinical Breast Imaging Reporting and Data System (BI-RADS) classifications of breast density were extracted from mammography reports. Odds ratios and 95% confidence intervals (CIs) were estimated by using conditional logistic regression stratified according to age and race and adjusted for body mass index, parity, and menopausal status, and the area under the receiver operating characteristic curve (AUC) was computed. Results The adjusted odds ratios and 95% CIs for each standard deviation increment of the percentage of density, dense area, volumetric percentage of density, and dense volume were 1.61 (95% CI: 1.19, 2.19), 1.49 (95% CI: 1.15, 1.92), 1.54 (95% CI: 1.12, 2.10), and 1.41 (95% CI: 1.11, 1.80), respectively. Odds ratios for women with extremely dense breasts compared with those with scattered areas of fibroglandular density were 2.06 (95% CI: 0.85, 4.97) and 2.05 (95% CI: 0.90, 4.64) for BI-RADS and Volpara density classifications, respectively. Clinical BI-RADS was more accurate (AUC, 0.68; 95% CI: 0.63, 0.74) than Volpara (AUC, 0.64; 95% CI: 0.58, 0.70) and continuous measures of percentage of density (AUC, 0.66; 95% CI: 0.60, 0.72), dense area (AUC, 0.66; 95% CI: 0.60, 0.72), volumetric percentage of density (AUC, 0.64; 95% CI: 0.58, 0.70), and density volume (AUC, 0.65; 95% CI: 0.59, 0.71), although the AUC differences were not statistically significant. Conclusion Mammographic density on FFDM images was positively associated with breast cancer risk by using the computer assisted methods and BI-RADS. BI-RADS classification was as accurate as computer-assisted methods for discrimination of patients from control subjects. ^© RSNA, 2016.

Effect of anatomical noise on the detectability of cone beam CT images with different slice direction, slice thickness, and volume glandular fraction

We investigate the effect of anatomical noise on the detectability of cone beam CT (CBCT) images with different slice directions, slice thicknesses, and volume glandular fractions (VGFs). Anatomical noise is generated using a power law spectrum of breast anatomy, and spherical objects with diameters from 1mm to 11mm are used as breast masses. CBCT projection images are simulated and reconstructed using the FDK algorithm. A channelized Hotelling observer (CHO) with Laguerre-Gauss (LG) channels is used to evaluate detectability for the signal-known-exactly (SKE) binary detection task. Detectability is calculated for various slice thicknesses in the transverse and longitudinal planes for 15%, 30% and 60% VGFs. The optimal slice thicknesses that maximize the detectability of the objects are determined. The results show that the β value increases as the slice thickness increases, but that thicker slices yield higher detectability in the transverse and longitudinal planes, except for the case of a 1mm diameter spherical object. It is also shown that the longitudinal plane with a 0.1mm slice thickness provides higher detectability than the transverse plane, despite its higher β value. With optimal slice thicknesses, the longitudinal plane exhibits better detectability for all VGFs and spherical objects.

Visibility of microcalcification in cone beam breast CT: effects of X-ray tube voltage and radiation dose

Mammography is the only technique currently used for detecting microcalcification (MC) clusters, an early indicator of breast cancer. However, mammographic images superimpose a three-dimensional compressed breast image onto two-dimensional projection views, resulting in overlapped anatomical breast structures that may obscure the detection and visualization of MCs. One possible solution to this problem is the use of cone beam computed tomography (CBCT) with a flat-panel (FP) digital detector. Although feasibility studies of CBCT techniques for breast imaging have yielded promising results, they have not shown how radiation dose and x-ray tube voltage affect the accuracy with which MCs are detected by CBCT experimentally. We therefore conducted a phantom study using a FP-based CBCT system with various mean glandular doses and kVp values. An experimental CBCT scanner was constructed with a data acquisition rate of 7.5 frames/s. 10.5 and 14.5 cm diameter breast phantoms made of gelatin were used to simulate uncompressed breasts consisting of 100% glandular tissue. Eight different MC sizes of calcium carbonate grains, ranging from 180-200 microm to 355-425 microm, were used to simulate MCs. MCs of the same size were arranged to form a 5 x 5 MC cluster and embedded in the breast phantoms. These MC clusters were positioned at 2.8 cm away from the center of the breast phantoms. The phantoms were imaged at 60, 80, and 100 kVp. With a single scan (360 degrees), 300 projection images were acquired with 0.5 x, 1x, and 2x mean glandular dose limit for 10.5 cm phantom and with 1x, 2x, and 4x for 14.5 cm phantom. A Feldkamp algorithm with a pure ramp filter was used for image reconstruction. The normalized noise level was calculated for each x-ray tube voltage and dose level. The image quality of the CBCT images was evaluated by counting the number of visible MCs for each MC cluster for various conditions. The average percentage of the visible MCs was computed and plotted as a function of the MGD, the kVp, and the average MC size. The results showed that the MC visibility increased with the MGD significantly but decreased with the breast size. The results also showed that the x-ray tube voltage affects the detection of MCs under different circumstances. With a 50% threshold, the minimum detectable MC sizes for the 10.5 cm phantom were 348(+/-2), 288(+/-7), 257(+/-2) microm at 3, 6, and 12 mGy, respectively. Those for the 14.5 cm phantom were 355 (+/-1), 307 (+/-7), 275 (+/-5) microm at 6, 12, and 24 mGy, respectively. With a 75% threshold, the minimum detectable MC sizes for the 10.5 cm phantom were 367 (+/-1), 316 (+/-7), 265 (+/-3) microm at 3, 6, and 12 mGy, respectively. Those for the 14.5 cm phantom were 377 (+/-3), 334 (+/-5), 300 (+/-2) microm at 6, 12, and 24 mGy, respectively.

Measurement of breast-tissue x-ray attenuation by spectral mammography: solid lesions

Knowledge of x-ray attenuation is essential for developing and evaluating x-ray imaging technologies. For instance, techniques to distinguish between cysts and solid tumours at mammography screening would be highly desirable to reduce recalls, but the development requires knowledge of the x-ray attenuation for cysts and tumours. We have previously measured the attenuation of cyst fluid using photon-counting spectral mammography. Data on x-ray attenuation for solid breast lesions are available in the literature, but cover a relatively wide range, likely caused by natural spread between samples, random measurement errors, and different experimental conditions. In this study, we have adapted a previously developed spectral method to measure the linear attenuation of solid breast lesions. A total of 56 malignant and 5 benign lesions were included in the study. The samples were placed in a holder that allowed for thickness measurement. Spectral (energy-resolved) images of the samples were acquired and the image signal was mapped to equivalent thicknesses of two known reference materials, which can be used to derive the x-ray attenuation as a function of energy. The spread in equivalent material thicknesses was relatively large between samples, which is likely to be caused mainly by natural variation and only to a minor extent by random measurement errors and sample inhomogeneity. No significant difference in attenuation was found between benign and malignant solid lesions. The separation between cyst-fluid and tumour attenuation was, however, significant, which suggests it may be possible to distinguish cystic from solid breast lesions, and the results lay the groundwork for a clinical trial. In addition, the study adds a relatively large sample set to the published data and may contribute to a reduction in the overall uncertainty in the literature.

BI-RADS Density Classification From Areometric and Volumetric Automatic Breast Density Measurements

RATIONALE AND OBJECTIVES

The aim of our study was to classify breast density using areometric and volumetric automatic measurements to best match Breast Imaging-Reporting and Data System (BI-RADS) density scores, and determine which technique best agrees with BI-RADS. Second, this study aimed to provide a set of threshold values for areometric and volumetric density to estimate BI-RADS categories.

MATERIALS AND METHODS

We randomly selected 537 full-field digital mammography examinations from a population-based screening program. Five radiologists assessed breast density using BI-RADS with all views available. A commercial program calculated areometric and volumetric breast density automatically. We compared automatically calculated density to all BI-RADS density thresholds using area under the receiver operating characteristic curve, and used Youden's index to estimate thresholds in automatic densities, with matching sensitivity and specificity. The 95% confidence intervals were estimated by bootstrapping.

RESULTS

Areometric density correlated well with volumetric density (r(2) = 0.76, excluding outliers, n = 2). For the BI-RADS threshold between II and III, areometric and volumetric assessment showed about equal area under the curve (0.94 vs. 0.93). For the threshold between I and II, areometric assessment was better than volumetric assessment (0.91 vs. 0.86). For the threshold between III and IV, volumetric assessment was better than areometric assessment (0.97 vs. 0.92).

CONCLUSIONS

Volumetric assessment is equal to or better than areometric assessment for the most clinically relevant thresholds (ie, between scattered fibroglandular and heterogeneously dense, and between heterogeneously dense and extremely dense breasts). Thresholds found in this study can be applied in daily practice to automatic measurements of density to estimate BI-RADS classification.

Model predictions for the wide-angle x-ray scatter signals of healthy and malignant breast duct biopsies

Wide-angle x-ray scatter (WAXS) could potentially be used to diagnose ductal carcinoma in situ (DCIS) in breast biopsies. The regions of interest were assumed to consist of fibroglandular tissue and epithelial cells and the model assumed that biopsies with DCIS would have a higher concentration of the latter. The scattered number of photons from a 2-mm diameter column of tissue was simulated using a 110-kV beam and selectively added in terms of momentum transfer. For a 1-min exposure, specificities and sensitivities of unity were obtained for biopsies 2- to 20-mm thick. The impact of sample and tumor cell layer thicknesses was studied. For example, a biopsy erroneously estimated to be 8 mm would be correctly diagnosed if its actual thickness was between 7.3 and 8.7 mm. An 8-mm thick malignant biopsy can be correctly diagnosed provided the malignant cell layer thickness is [Formula: see text]. WAXS methods could become a diagnostic tool for DCIS within breast biopsies.

A computational model to generate simulated three-dimensional breast masses

PURPOSE

To develop algorithms for creating realistic three-dimensional (3D) simulated breast masses and embedding them within actual clinical mammograms. The proposed techniques yield high-resolution simulated breast masses having randomized shapes, with user-defined mass type, size, location, and shape characteristics.

METHODS

The authors describe a method of producing 3D digital simulations of breast masses and a technique for embedding these simulated masses within actual digitized mammograms. Simulated 3D breast masses were generated by using a modified stochastic Gaussian random sphere model to generate a central tumor mass, and an iterative fractal branching algorithm to add complex spicule structures. The simulated masses were embedded within actual digitized mammograms. The authors evaluated the realism of the resulting hybrid phantoms by generating corresponding left- and right-breast image pairs, consisting of one breast image containing a real mass, and the opposite breast image of the same patient containing a similar simulated mass. The authors then used computer-aided diagnosis (CAD) methods and expert radiologist readers to determine whether significant differences can be observed between the real and hybrid images.

RESULTS

The authors found no statistically significant difference between the CAD features obtained from the real and simulated images of masses with either spiculated or nonspiculated margins. Likewise, the authors found that expert human readers performed very poorly in discriminating their hybrid images from real mammograms.

CONCLUSIONS

The authors' proposed method permits the realistic simulation of 3D breast masses having user-defined characteristics, enabling the creation of a large set of hybrid breast images containing a well-characterized mass, embedded within real breast background. The computational nature of the model makes it suitable for detectability studies, evaluation of computer aided diagnosis algorithms, and teaching purposes.

Contrast-Enhanced X-ray Detection of Microcalcifications in Radiographically Dense Mammary Tissue Using Targeted Gold Nanoparticles

Breast density reduces the accuracy of mammography, motivating methods to improve sensitivity and specificity for detecting abnormalities within dense breast tissue, but preclinical animal models are lacking. Therefore, the objectives of this study were to investigate a murine model of radiographically dense mammary tissue and contrast-enhanced X-ray detection of microcalcifications in dense mammary tissue by targeted delivery of bisphosphonate-functionalized gold nanoparticles (BP-Au NPs). Mammary glands (MGs) in the mouse mammary tumor virus - polyomavirus middle T antigen (MMTV-PyMT or PyMT) model exhibited greater radiographic density with age and compared with strain- and age-matched wild-type (WT) controls at 6-10 weeks of age. The greater radiographic density of MGs in PyMT mice obscured radiographic detection of microcalcifications that were otherwise detectable in MGs of WT mice. However, BP-Au NPs provided enhanced contrast for the detection of microcalcifications in both radiographically dense (PyMT) and WT mammary tissues as measured by computed tomography after intramammary delivery. BP-Au NPs targeted microcalcifications to enhance X-ray contrast with surrounding mammary tissue, which resulted in improved sensitivity and specificity for detection in radiographically dense mammary tissues.

Characterization of structural-prior guided optical tomography using realistic breast models derived from dual-energy x-ray mammography

Multi-spectral near-infrared diffuse optical tomography (DOT) is capable of providing functional tissue assessment that can complement structural mammographic images for more comprehensive breast cancer diagnosis. To take full advantage of the readily available sub-millimeter resolution structural information in a multi-modal imaging setting, an efficient x-ray/optical joint image reconstruction model has been proposed previously to utilize anatomical information from a mammogram as a structural prior. In this work, we develop a complex digital breast phantom (available at http://openjd.sf.net/digibreast) based on direct measurements of fibroglandular tissue volume fractions using dual-energy mammographic imaging of a human breast. We also extend our prior-guided reconstruction algorithm to facilitate the recovery of breast tumors, and perform a series of simulation-based studies to systematically evaluate the impact of lesion sizes and contrasts, tissue background, mesh resolution, inaccurate priors, and regularization parameters, on the recovery of breast tumors using multi-modal DOT/x-ray measurements. Our studies reveal that the optical property estimation error can be reduced by half by utilizing structural priors; the minimum detectable tumor size can also be reduced by half when prior knowledge regarding the tumor location is provided. Moreover, our algorithm is shown to be robust to false priors on tumor location.

Task-based optimization of dedicated breast CT via Hotelling observer metrics

PURPOSE

The purpose of this work is to develop and demonstrate a set of practical metrics for CT systems optimization. These metrics, based on the Hotelling observer (HO) figure of merit, are task-based. The authors therefore take the specific example of optimizing a dedicated breast CT system, including the reconstruction algorithm, for two relevant tasks, signal detection and Rayleigh discrimination.

METHODS

A dedicated breast CT system is simulated using specifications in the literature from an existing prototype. The authors optimize configuration and image reconstruction algorithm parameters for two tasks: the detection of simulated microcalcifications and the discrimination of two adjacent, high-contrast signals, known as the Rayleigh discrimination task. The effects on task performance of breast diameter, signal location, image grid size, projection view number, and reconstruction filter were all investigated. Two HO metrics were evaluated: the percentage of correct decisions in a two-alternative forced choice experiment (equivalent to area under the ROC curve or AUC), and the HO efficiency, defined as the squared ratio of HO signal-to-noise ratio (SNR) in the reconstructed image to HO SNR in the projection data.

RESULTS

The ease and efficiency of the HO metric computation allows a rapid high-resolution survey of many system parameters. Optimization of a range of system parameters using the HO results in images that subjectively appear optimal for the tasks investigated. Further, the results of assessment through the HO reproduce closely many existing results in the literature regarding the impact of parameter selection on image quality.

CONCLUSIONS

This study demonstrates the utility of a task-based approach to system design, evaluation, and optimization. The methodology presented is equally applicable to determining the impact of a wide range of factors, including patient parameters, system and acquisition design, and the reconstruction algorithm. The results demonstrate the versatility of the proposed HO formalism by not only generating a set of parameters that are optimal for a given task but also by qualitatively reproducing many existing results from the breast CT literature. Meanwhile, the implementation of the proposed methodology is straightforward and entirely simulation-based. This is an attractive feature for many system optimization problems, where the goal is to analyze the individual system components such as the image reconstruction algorithm. Final assessment of the system as a whole should be based also on real data studies.

Mammographic breast density: comparison of methods for quantitative evaluation

PURPOSE

To evaluate the results from two software tools for measurement of mammographic breast density and compare them with observer-based scores in a large cohort of women.

MATERIALS AND METHODS

Following written informed consent, a data set of 36 281 mammograms from 8867 women were collected from six United Kingdom centers in an ethically approved trial. Breast density was assessed by one of 26 readers on a visual analog scale and with two automated density tools. Mean differences were calculated as the mean of all the individual percentage differences between each measurement for each case (woman). Agreement in total breast volume, fibroglandular volume, and percentage density was assessed with the Bland-Altman method. Association with observer's scores was calculated by using the Pearson correlation coefficient (r).

RESULTS

Correlation between the Quantra and Volpara outputs for total breast volume was r = 0.97 (P < .001), with a mean difference of 43.5 cm(3) for all cases representing 5.0% of the mean total breast volume. Correlation of the two measures was lower for fibroglandular volume (r = 0.86, P < .001). The mean difference was 30.3 cm(3) for all cases representing 21.2% of the mean fibroglandular tissue volume result. Quantra gave the larger value and the difference tended to increase with volume. For the two measures of percentage volume density, the mean difference was 1.61 percentage points (r = 0.78, P < .001). Comparison of observer's scores with the area-based density given by Quantra yielded a low correlation (r = 0.55, P < .001). Correlations of observer's scores with the volumetric density results gave r values of 0.60 (P < .001) and 0.63 (P < .001) for Quantra and Volpara, respectively.

CONCLUSION

Automated techniques for measuring breast density show good correlation, but these are poorly correlated with observer's scores. However automated techniques do give different results that should be considered when informing patient personalized imaging. (©) RSNA, 2015 Clinical trial registration no. ISRCTN 73467396.

Evaluating the sensitivity of the optimization of acquisition geometry to the choice of reconstruction algorithm in digital breast tomosynthesis through a simulation study

Due to the limited number of views and limited angular span in digital breast tomosynthesis (DBT), the acquisition geometry design is an important factor that affects the image quality. Therefore, intensive studies have been conducted regarding the optimization of the acquisition geometry. However, different reconstruction algorithms were used in most of the reported studies. Because each type of reconstruction algorithm can provide images with its own image resolution, noise properties and artifact appearance, it is unclear whether the optimal geometries concluded for the DBT system in one study can be generalized to the DBT systems with a reconstruction algorithm different to the one applied in that study. Hence, we investigated the effect of the reconstruction algorithm on the optimization of acquisition geometry parameters through carefully designed simulation studies. Our results show that using various reconstruction algorithms, including the filtered back-projection, the simultaneous algebraic reconstruction technique, the maximum-likelihood method and the total-variation regularized least-square method, gave similar performance trends for the acquisition parameters for detecting lesions. The consistency of system ranking indicates that the choice of the reconstruction algorithm may not be critical for DBT system geometry optimization.

Region of interest based Hotelling observer for computed tomography with comparison to alternative methods

We compare several approaches to estimation of Hotelling observer (HO) performance in x-ray computed tomography (CT). We consider the case where the signal of interest is small so that the reconstructed image can be restricted to a small region of interest (ROI) surrounding the signal. This reduces the dimensionality of the image covariance matrix so that direct computation of HO metrics within the ROI is feasible. We propose that this approach is directly applicable to systems optimization in CT; however, many alternative approaches exist, which make computation of HO performance tractable through a range of approximations, assumptions, or estimation strategies. Here, we compare several of these methods, including the use of Laguerre-Gauss channels, discrete Fourier domain computation of the HO (which assumes noise stationarity), and two approaches to HO estimation through samples of noisy images. Since our method computes HO performance exactly within an ROI, this allows us to investigate the validity of the assumptions inherent in various common approaches to HO estimation, such as the stationarity assumption in the case of the discrete Fourier transform domain method.

Three-dimensional x-ray fluorescence mapping of a gold nanoparticle-loaded phantom

PURPOSE

X-ray fluorescence (XRF) is a promising technique with sufficient specificity and sensitivity for identifying and quantifying features in small samples containing high atomic number (Z) materials such as iodine, gadolinium, and gold. In this study, the feasibility of applying XRF to early breast cancer diagnosis and treatment is studied using a novel approach for three-dimensional (3D) x-ray fluorescence mapping (XFM) of gold nanoparticle (GNP)-loaded objects in a physical phantom at the technical level.

METHODS

All the theoretical analysis and experiments are conducted under the condition of using x-ray pencil beam and a compactly integrated x-ray spectrometer. The penetrability of the fluorescence x-rays from GNPs is first investigated by adopting a combination of BR12 with 70 mm/50 mm in thickness on the excitation/emission path to mimic the possible position of tumor goldin vivo. Then, a physical phantom made of BR12 is designed to translate in 3D space with three precise linear stages and subsequently the step by step XFM scanning is performed. The experimental technique named as background subtraction is applied to isolate the gold fluorescence from each spectrum obtained by the spectrometer. Afterwards, the attenuations of both the incident primary x-ray beam with energies beyond the gold K-edge energy (80.725 keV) and the isolated gold Kα fluorescence x-rays (65.99 -69.80 keV) acquired after background subtraction are well calibrated, and finally the unattenuated Kα fluorescence counts are used to realize mapping reconstruction and to describe the linear relationship between gold fluorescence counts and corresponding concentration of gold solutions.

RESULTS

The penetration results show that the goldKα fluorescence x-rays have sufficient penetrability for this phantom study, and the reconstructed mapping results indicate that both the spatial distribution and relative concentration of GNPs within the designed BR12 phantom can be well identified and quantified.

CONCLUSIONS

Although the XFM method in this investigation is still studied at the technical level and is not yet practical for routinein vivo mapping tasks with GNPs, the current penetrability measurements and phantom study strongly suggest the feasibility to establish and develop a 3D XFM system.