Search

Dedicated breast CT is an imaging modality that provides true 3D imaging of the breast with many advantages over current conventional breast imaging modalities. The addition of intravascular contrast increases the sensitivity of breast CT substantially. As such, there are immediate potential applications in the clinical workflow. These include using breast CT to replace much of the traditional diagnostic workup when faced with indeterminate breast lesions. Contrast-enhanced breast CT may be appropriate as a supplemental screening tool for women at high risk of breast cancer, similar to breast MRI. In addition, emerging studies are demonstrating the utility of breast CT in neoadjuvant chemotherapy tumor response monitoring as well as planning for surgical treatment options. While short exam times and fully 3D imaging in a noncompressed position are advantages of this modality, limited coverage of chest wall/axilla due to prone positioning and use of ionizing radiation are drawbacks. To date, several studies have reported on the performance characteristics of this promising modality. [ABSTRACT FROM AUTHOR]

Background: X‐ray breast imaging modalities are commonly employed for breast cancer detection, from screening programs to diagnosis. Thus, dosimetry studies are important for quality control and risk estimation since ionizing radiation is used. Purpose: To perform multiscale dosimetry assessments for different breast imaging modalities and for a variety of breast sizes and compositions. The first part of our study is focused on macroscopic scales (down to millimeters). Methods: Nine anthropomorphic breast phantoms with a voxel resolution of 0.5 mm were computationally generated using the BreastPhantom software, representing three breast sizes with three distinct values of volume glandular fraction (VGF) for each size. Four breast imaging modalities were studied: digital mammography (DM), contrast‐enhanced digital mammography (CEDM), digital breast tomosynthesis (DBT) and dedicated breast computed tomography (BCT). Additionally, the impact of tissue elemental compositions from two databases were compared. Monte Carlo (MC) simulations were performed with the MC‐GPU code to obtain the 3D glandular dose distribution (GDD) for each case considered with the mean glandular dose (MGD) fixed at 4 mGy (to facilitate comparisons). Results: The GDD within the breast is more uniform for CEDM and BCT compared to DM and DBT. For large breasts and high VGF, the ratio between the minimum/maximum glandular dose to MGD is 0.12/4.02 for DM and 0.46/1.77 for BCT; the corresponding results for a small breast and low VGF are 0.35/1.98 (DM) and 0.63/1.42 (BCT). The elemental compositions of skin, adipose and glandular tissue have a considerable impact on the MGD, with variations up to 30% compared to the baseline. The inclusion of tissues other than glandular and adipose within the breast has a minor impact on MGD, with differences below 2%. Variations in the final compressed breast thickness alter the shape of the GDD, with a higher compression resulting in a more uniform GDD. Conclusions: For a constant MGD, the GDD varies with imaging modality and breast compression. Elemental tissue compositions are an important factor for obtaining MGD values, being a source of systematic uncertainties in MC simulations and, consequently, in breast dosimetry. [ABSTRACT FROM AUTHOR]

Background: Breast computed tomography (CT) is an emerging breast imaging modality, and ongoing developments aim to improve breast CT's ability to detect microcalcifications. To understand the effects of different parameters on microcalcification detectability, a virtual clinical trial study was conducted using hybrid images and convolutional neural network (CNN)‐based model observers. Mathematically generated microcalcifications were embedded into breast CT data sets acquired at our institution, and parameters related to calcification size, calcification contrast, cluster diameter, cluster density, and image display method (i.e., single slices, slice averaging, and maximum‐intensity projections) were evaluated for their influence on microcalcification detectability. Purpose: To investigate the individual effects and the interplay of parameters affecting microcalcification detectability in breast CT. Methods: Spherical microcalcifications of varying diameters (0.04, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40 mm) and native intensities were computer simulated to portray the partial volume effects of the imaging system. Calcifications were mathematically embedded into 109 patient breast CT volume data sets as individual calcifications or as clusters of calcifications. Six numbers of calcifications (1, 3, 5, 7, 10, 15) distributed within six cluster diameters (1, 3, 5, 6, 8, 10 mm) were simulated to study the effect of cluster density. To study the role of image display method, 2D regions of interest (ROIs) and 3D volumes of interest (VOIs) were generated using single slice extraction, slice averaging, and maximum‐intensity projection (MIP). 2D and 3D CNNs were trained on the ROIs and VOIs, and receiver operating characteristic (ROC) curve analysis was used to evaluate detection performance. The area under the ROC curve (AUC) was used as the primary performance metric. Results: Detection performance decreased with increasing section thickness, and peak detection performance occurred using the native section thickness (0.2 mm) and MIP display. The MIP display method, despite using a single slice, yielded comparable performance to the native section thickness, which employed 50 slices. Reduction in slices did not sacrifice detection accuracy and provided significant computational advantages over multi‐slice image volumes. Larger cluster diameters resulted in reduced overall detectability, while smaller cluster diameters led to increased detectability. Additionally, we observed that the presence of more calcifications within a cluster improved the overall detectability, while fewer calcifications decreased it. Conclusions: As breast CT is still a relatively new breast imaging modality, there is an ongoing need to identify optimal imaging protocols. This work demonstrated the utility of MIP presentation for displaying image volumes containing microcalcification clusters. It is likely that human observers may also benefit from viewing MIPs compared to individual slices. The results of this investigation begin to elucidate how model observers interact with microcalcification clusters in a 3D volume, and will be useful for future studies investigating a broader set of parameters related to breast CT. [ABSTRACT FROM AUTHOR]

Dedicated cone-beam breast computed tomography (CBBCT) is an emerging modality and provides fully three-dimensional (3D) images of the uncompressed breast at an isotropic voxel resolution. In an effort to translate this modality to breast cancer screening, advanced image reconstruction methods are being pursued. Since radiographic breast density is an established risk factor for breast cancer and CBBCT provides volumetric data, this study investigates the reproducibility of the volumetric glandular fraction (VGF), defined as the proportion of fibroglandular tissue volume relative to the total breast volume excluding the skin. Four image reconstruction methods were investigated: the analytical Feldkamp–Davis–Kress (FDK), a compressed sensing-based fast, regularized, iterative statistical technique (FRIST), a fully supervised deep learning approach using a multi-scale residual dense network (MS-RDN), and a self-supervised approach based on Noise-to-Noise (N2N) learning. Projection datasets from 106 women who participated in a prior clinical trial were reconstructed using each of these algorithms at a fixed isotropic voxel size of (0.273 mm3). Each reconstructed breast volume was segmented into skin, adipose, and fibroglandular tissues, and the VGF was computed. The VGF did not differ among the four reconstruction methods (p = 0.167), and none of the three advanced image reconstruction algorithms differed from the standard FDK reconstruction (p > 0.862). Advanced reconstruction algorithms developed for low-dose CBBCT reproduce the VGF to provide quantitative breast density, which can be used for risk estimation.

Dedicated cone-beam breast computed tomography (CBBCT) is an emerging modality and provides fully three-dimensional (3D) images of the uncompressed breast at an isotropic voxel resolution. In an effort to translate this modality to breast cancer screening, advanced image reconstruction methods are being pursued. Since radiographic breast density is an established risk factor for breast cancer and CBBCT provides volumetric data, this study investigates the reproducibility of the volumetric glandular fraction (VGF), defined as the proportion of fibroglandular tissue volume relative to the total breast volume excluding the skin. Four image reconstruction methods were investigated: the analytical Feldkamp–Davis–Kress (FDK), a compressed sensing-based fast, regularized, iterative statistical technique (FRIST), a fully supervised deep learning approach using a multi-scale residual dense network (MS-RDN), and a self-supervised approach based on Noise-to-Noise (N2N) learning. Projection datasets from 106 women who participated in a prior clinical trial were reconstructed using each of these algorithms at a fixed isotropic voxel size of (0.273 mm3). Each reconstructed breast volume was segmented into skin, adipose, and fibroglandular tissues, and the VGF was computed. The VGF did not differ among the four reconstruction methods (p = 0.167), and none of the three advanced image reconstruction algorithms differed from the standard FDK reconstruction (p > 0.862). Advanced reconstruction algorithms developed for low-dose CBBCT reproduce the VGF to provide quantitative breast density, which can be used for risk estimation. [ABSTRACT FROM AUTHOR]

Objectives: Breast density is considered an independent risk factor for the development of breast cancer. This study aimed to quantitatively assess the percent breast density (PBD) and the mammary glands volume (MGV) according to the patient's age and breast quadrant. We propose a regression model to estimate PBD and MGV as a function of the patient's age. Methods: The breast composition in 1027 spiral breast CT (BCT) datasets without soft tissue masses, calcifications, or implants from 517 women (57 ± 8 years) were segmented. The breast tissue volume (BTV), MGV, and PBD of the breasts were measured in the entire breast and each of the four quadrants. The three breast composition features were analyzed in the seven age groups, from 40 to 74 years in 5-year intervals. A logarithmic model was fitted to the BTV, and a multiplicative inverse model to the MGV and PBD as a function of age was established using a least-squares method. Results: The BTV increased from 545 ± 345 to 676 ± 412 cm3, and the MGV and PBD decreased from 111 ± 164 to 57 ± 43 cm3 and from 21 ± 21 to 11 ± 9%, respectively, from the youngest to the oldest group (p < 0.05). The average PBD over all ages were 14 ± 13%. The regression models could predict the BTV, MGV, and PBD based on the patient's age with residual standard errors of 386 cm3, 67 cm3, and 13%, respectively. The reduction in MGV and PBD in each quadrant followed the ones in the entire breast. Conclusions: The PBD and MGV computed from BCT examinations provide important information for breast cancer risk assessment in women. The study quantified the breast mammary gland reduction and density decrease over the entire breast. It established mathematical models to estimate the breast composition features—BTV, MGV, and PBD, as a function of the patient's age. [ABSTRACT FROM AUTHOR]

Background: Radiopaque breast markers cause artifacts in dedicated spiral breast-computed tomography (SBCT). This study investigates the extent of artifacts in different marker types and the feasibility of reducing artifacts through a metal artifact reduction (MAR) algorithm. Methods: The pilot study included 18 women who underwent contrast-enhanced SBCT. In total, 20 markers of 4 different types were analyzed for artifacts. The extent of artifacts with and without MAR was measured via the consensus of two readers. Image noise was quantitatively evaluated, and the effect of MAR on the detectability of breast lesions was evaluated on a 3-point Likert scale. Results: Breast markers caused significant artifacts that impaired image quality and the detectability of lesions. MAR decreased artifact size in all analyzed cases, even in cases with multiple markers in a single slice. The median length of in-plain artifacts significantly decreased from 31 mm (range 11–51 mm) in uncorrected to 2 mm (range 1–5 mm) in corrected images (p ≤ 0.05). Artifact size was dependent on marker size. Image noise in slices affected by artifacts was significantly lower in corrected (13.6 ± 2.2 HU) than in uncorrected images (19.2 ± 6.8 HU, p ≤ 0.05). MAR improved the detectability of lesions affected by artifacts in 5 out of 11 cases. Conclusion: MAR is feasible in SBCT and improves the image quality and detectability of lesions. [ABSTRACT FROM AUTHOR]

PurposeTo develop and evaluate the diagnostic performance of an algorithm for multi-marker radiomic-based classification of breast masses in dedicated breast computed tomography (bCT) images.MethodsOver 1000 radiomic descriptors aimed at quantifying mass and border heterogeneity, morphology, and margin sharpness were developed and implemented. These included well-established texture and shape feature descriptors, which were supplemented with additional approaches for contour irregularity quantification, spicule and lobe detection, characterization of degree of infiltration, and differences in peritumoral compartments. All descriptors were extracted from a training set of 202 bCT masses (133 benign and 69 malignant), and their individual diagnostic performance was investigated in terms of area under the receiver operating characteristics (ROC) curve (AUC) of single-feature-based linear discriminant analysis (LDA) classifiers. Subsequently, the most relevant descriptors were selected through a multiple-step feature selection process (including stability analysis, statistical significance, evaluation of feature interaction, and dimensionality reduction), and used to develop a final LDA radiomic model for classification of benign and malignant masses, which was then tested on an independent test set of 82 cases (45 benign and 37 malignant).ResultsThe majority of the individual radiomic descriptors showed, on the training set, an AUC value deriving from a linear decision boundary higher than 0.65, with the lower limit of the associated 95% confidence interval (C.I.) not overlapping with random chance (AUC = 0.5). The final LDA radiomic model resulted in a test set AUC of 0.90 (95% C.I. 0.80-0.96).ConclusionsThe proposed multi-marker radiomic approach achieved high diagnostic accuracy in bCT mass classification, using a radiomic signature based on different feature types. While future studies with larger datasets are needed to further validate these results, quantitative radiomics applied to bCT shows potential to improve the breast cancer diagnosis pipeline.

Background: In breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior studies provided evidence for a new image acquisition design, dubbed Narrow Beam Breast CT (NB‐bCT), in preventing scatter acquisition. Purpose: Here, we report the design, implementation, and initial characterization of the first NB‐bCT prototype. Methods: The imaging system's apparatus is composed of two primary assemblies: a dynamic Fluence Modulator (collimator) and a photon‐counting line detector. The design of the assemblies enables them to operate in lockstep during image acquisition, converting sourced x‐rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. The assemblies are each comprised of a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. To study the impact of fluence modulation on the signal received in the detector, three physical breast phantoms, representative of small, average, and large size breasts, were developed and imaged, and acquired projections analyzed. The scatter acquisition in each projection as a function of breast phantom size was investigated. The imaging system's spatial resolution at the center and periphery of the field of view was measured. Results: Minimal acquisition of scattered rays occurs during image acquisition with NB‐bCT; results in minimal scatter to primary ratios in small, average, and large breast phantoms imaged were 0.05, 0.07, and 0.9, respectively. System spatial resolution of 5.2 lp/mm at 10% max MTF and 2.9 lp/mm at 50% max MTF at the center of the field of view was achieved, with minimal loss with the shift toward the corner (5.0 lp/mm at 10% max MTF and 2.5 lp/mm at 50% max MTF). Conclusion: The disclosed development, implementation, and characterization of a physical NB‐bCT prototype system demonstrates a new method of CT‐based image acquisition that yields high spatial resolution while minimizing scatter‐components in acquired projections. This methodology holds promise for high‐resolution CT‐imaging applications in which reduction of scatter contamination is desirable. [ABSTRACT FROM AUTHOR]

PurposeClinical use of dedicated breast computed tomography (bCT) requires relatively short scan times necessitating systems with high frame rates. This in turn impacts the x-ray tube operating range. We characterize the effects of tube voltage, beam filtration, dose, and object size on contrast and noise properties related to soft tissue and iodine contrast agents as a way to optimize imaging protocols for soft tissue and iodine contrast at high frame rates.MethodsThis study design uses the signal-difference-to-noise ratio (SDNR), noise-equivalent quanta (NEQ), and detectability (d´) as measures of imaging performance for a prototype breast CT scanner that utilizes a pulsed x-ray tube (with a 4 ms pulse width) at 43.5 fps acquisition rate. We assess a range of kV, filtration, breast phantom size, and mean glandular dose (MGD). Performance measures are estimated from images of adipose-equivalent breast phantoms machined to have a representative size and shape of small, medium, and large breasts. Water (glandular tissue equivalent) and iodine contrast (5 mg/ml) were used to fill two cylindrical wells in the phantoms.ResultsAir kerma levels required for obtaining an MGD of 6 mGy ranged from 7.1 to 9.1 mGy and are reported across all kV, filtration, and breast phantom sizes. However, at 50 kV, the thick filters (0.3 mm of Cu or Gd) exceeded the maximum available mA of the x-ray generator, and hence, these conditions were excluded from subsequent analysis. There was a strong positive association between measurements of SDNR and d' (R2  > 0.97) within the range of parameters investigated in this work. A significant decrease in soft tissue SDNR was observed for increasing phantom size and increasing kV with a maximum SDNR at 50 kV with 0.2 mm Cu or 0.2 mm Gd filtration. For iodine contrast SDNR, a significant decrease was observed with increasing phantom size, but a decrease in SDNR for increasing kV was only observed for 70 kV (50 and 60 kV were not significantly different). Thicker Gd filtration (0.3 mm Gd) resulted in a significant increase in iodine SDNR and decrease in soft tissue SDNR but requires significantly more tube current to deliver the same MGD.ConclusionsThe choice of 60 kV with 0.2 mm Gd filtration provides a good trade-off for maximizing both soft tissue and iodine contrast. This scanning technique takes advantage of the ~50 keV Gd k-edge to produce contrast and can be achieved within operating range of the x-ray generator used in this work. Imaging at 60 kV allows for a greater range in dose delivered to the large breast sizes when uniform image quality is desired across all breast sizes. While imaging performance metrics (i.e., detectability index and SDNR) were shown to be strongly correlated, the methodologies presented in this work for the estimation of NEQ (and subsequently d') provides a meaningful description of the spatial resolution and noise characteristics of this prototype bCT system across a range of beam quality, dose, and object sizes.

To assess the density values of breast lesions and breast tissue using non-contrast spiral breast CT (nc-SBCT) imaging. In this prospective study women undergoing nc-SBCT between April-October 2023 for any purpose were included in case of: histologically proven malignant lesion (ML); fibroadenoma (FA) with histologic confirmation or stability > 24 months (retrospectively); cysts with ultrasound correlation; and women with extremely dense breast (EDB) and no sonographic findings. Three regions of interest were placed on each lesion and 3 different area of EDB. The evaluation was performed by two readers (R1 and R2). Kruskal-Wallis test, intraclass correlation (ICC) and ROC analysis were used. 40 women with 12 ML, 10 FA, 15 cysts and 9 with EDB were included. Median density values and interquartile ranges for R1 and R2 were: 60.2 (53.3–67.3) and 62.5 (55.67–76.3) HU for ML; 46.3 (41.9–59.5) and 44.5 (40.5–59.8) HU for FA; 35.3 (24.3–46.0) and 39.7 (26.7–52.0) HU for cysts; and 28.7 (24.2–33.0) and 33.3 (31.7–36.8) HU for EDB. For both readers, densities were significantly different for ML versus EDB (p < 0.001) and cysts (p < 0.001) and for FA versus EDB (p=/<0.003). The AUC was 0.925 (95 %CI 0.858–0.993) for R1 and 0.942 (0.884–1.00) for R2 when comparing ML versus others and 0.792 (0.596–0.987) and 0.833 (0.659–1) when comparing ML versus FA. The ICC showed an almost perfect inter-reader (0.978) and intra-reader agreement (>0.879 for both readers). In nc-SBCT malignant lesions have higher density values compared to normal tissue and measurements of density values are reproducible between different readers. [ABSTRACT FROM AUTHOR]

Objectives: Breast density is considered an independent risk factor for the development of breast cancer. This study aimed to quantitatively assess the percent breast density (PBD) and the mammary glands volume (MGV) according to the patient’s age and breast quadrant. We propose a regression model to estimate PBD and MGV as a function of the patient’s age. Methods: The breast composition in 1027 spiral breast CT (BCT) datasets without soft tissue masses, calcifications, or implants from 517 women (57 ± 8 years) were segmented. The breast tissue volume (BTV), MGV, and PBD of the breasts were measured in the entire breast and each of the four quadrants. The three breast composition features were analyzed in the seven age groups, from 40 to 74 years in 5-year intervals. A logarithmic model was fitted to the BTV, and a multiplicative inverse model to the MGV and PBD as a function of age was established using a least-squares method. Results: The BTV increased from 545 ± 345 to 676 ± 412 cm3, and the MGV and PBD decreased from 111 ± 164 to 57 ± 43 cm3 and from 21 ± 21 to 11 ± 9%, respectively, from the youngest to the oldest group (p < 0.05). The average PBD over all ages were 14 ± 13%. The regression models could predict the BTV, MGV, and PBD based on the patient’s age with residual standard errors of 386 cm3, 67 cm3, and 13%, respectively. The reduction in MGV and PBD in each quadrant followed the ones in the entire breast. Conclusions: The PBD and MGV computed from BCT examinations provide important information for breast cancer risk assessment in women. The study quantified the breast mammary gland reduction and density decrease over the entire breast. It established mathematical models to estimate the breast composition features—BTV, MGV, and PBD, as a function of the patient’s age.

Background: Radiopaque breast markers cause artifacts in dedicated spiral breast-computed tomography (SBCT). This study investigates the extent of artifacts in different marker types and the feasibility of reducing artifacts through a metal artifact reduction (MAR) algorithm. Methods: The pilot study included 18 women who underwent contrast-enhanced SBCT. In total, 20 markers of 4 different types were analyzed for artifacts. The extent of artifacts with and without MAR was measured via the consensus of two readers. Image noise was quantitatively evaluated, and the effect of MAR on the detectability of breast lesions was evaluated on a 3-point Likert scale. Results: Breast markers caused significant artifacts that impaired image quality and the detectability of lesions. MAR decreased artifact size in all analyzed cases, even in cases with multiple markers in a single slice. The median length of in-plain artifacts significantly decreased from 31 mm (range 11–51 mm) in uncorrected to 2 mm (range 1–5 mm) in corrected images (p ≤ 0.05). Artifact size was dependent on marker size. Image noise in slices affected by artifacts was significantly lower in corrected (13.6 ± 2.2 HU) than in uncorrected images (19.2 ± 6.8 HU, p ≤ 0.05). MAR improved the detectability of lesions affected by artifacts in 5 out of 11 cases. Conclusion: MAR is feasible in SBCT and improves the image quality and detectability of lesions.

Many computer aided diagnosis (CADx) tools for breast cancer begin by fully or semiautomatically segmenting a given breast lesion and then classifying the lesion's likelihood of malignancy using quantitative features extracted from the image. It is often assumed that better segmentation will result in better classification. However, this has not been thoroughly evaluated. The purpose of this study is to evaluate the relationship between computer segmentation performance and computer classification performance. We used 85 breast lesions (32 benign, 56 malignant) from breast computed tomography (CT) cases of 82 women. We prepared one smooth and one sharp iterative image reconstructions (IIR) and a clinical reconstruction for each of the 82 breast CT scans. For each reconstruction, we created 15 segmentation outcomes by applying 15 different segmentation algorithms. Specifically, we simulated 15 segmentation algorithms by changing parameters in a single segmentation algorithm. We then created 15 classification outcomes by conducting quantitative image feature analysis on the segmented image results. Using a 10-fold cross-validation, we evaluated the relationship between segmentation and classification performances. We found a low positive correlation between segmentation and classification performances for the smooth IIR (median Pearson's rho = 0.18), while a moderate positive correlation (median Pearson's rho = 0.4-0.43) was found between the two performances for the sharp IIR and clinical reconstruction. However, we found large variations in both segmentation and classification performances for the sharp IIR and clinical reconstruction. There were cases where segmentation algorithms resulted in similar segmentation performances, but the corresponding classification performances were different. These results indicate that an improvement in segmentation performance does not guarantee an improvement in the corresponding classification performance. Computer segmentation is an indirect variable affecting the computer classification. As better segmentation does not guarantee better classification, we should report both segmentation and classification performances when comparing segmentation algorithms.

We proposed the neutrosophic approach for segmenting breast lesions in breast computed tomography (bCT) images. The neutrosophic set considers the nature and properties of neutrality (or indeterminacy). We considered the image noise as an indeterminate component while treating the breast lesion and other breast areas as true and false components. We iteratively smoothed and contrast-enhanced the image to reduce the noise level of the true set. We then applied one existing algorithm for bCT images, the RGI segmentation, on the resulting noise-reduced image to segment the breast lesions. We compared the segmentation performance of the proposed method (named as NS-RGI) to that of the regular RGI segmentation. We used 122 breast lesions (44 benign and 78 malignant) of 111 noncontrast enhanced bCT cases. We measured the segmentation performances of the NS-RGI and the RGI using the Dice coefficient. The average Dice values of the NS-RGI and RGI were 0.82 and 0.80, respectively, and their difference was statistically significant ([Formula: see text]). We conducted a subsequent feature analysis on the resulting segmentations. The classifier performance for the NS-RGI ([Formula: see text]) improved over that of the RGI ([Formula: see text], [Formula: see text]).

Purpose: Breast cancer is the most common malignancy in women. Unfortunately, current breast imaging techniques all suffer from certain limitations: they are either not fully three dimensional, have an insufficient resolution or low soft‐tissue contrast. Grating interferometry breast computed tomography (GI‐BCT) is a promising X‐ray phase contrast modality that could overcome these limitations by offering high soft‐tissue contrast and excellent three‐dimensional resolution. To enable the transition of this technology to clinical practice, dedicated data‐processing algorithms must be developed in order to effectively retrieve the signals of interest from the measured raw data. Methods: This article proposes a novel denoising algorithm that can cope with the high‐noise amplitudes and heteroscedasticity which arise in GI‐BCT when operated in a low‐dose regime to effectively regularize the ill‐conditioned GI‐BCT inverse problem. We present a data‐driven algorithm called INSIDEnet, which combines different ideas such as multiscale image processing, transform‐domain filtering, transform learning, and explicit orthogonality to build an Interpretable NonexpanSIve Data‐Efficient network (INSIDEnet). Results: We apply the method to simulated breast phantom datasets and to real data acquired on a GI‐BCT prototype and show that the proposed algorithm outperforms traditional state‐of‐the‐art filters and is competitive with deep neural networks. The strong inductive bias given by the proposed model's architecture allows to reliably train the algorithm with very limited data while providing high model interpretability, thus offering a great advantage over classical convolutional neural networks (CNNs). Conclusions: The proposed INSIDEnet is highly data‐efficient, interpretable, and outperforms state‐of‐the‐art CNNs when trained on very limited training data. We expect the proposed method to become an important tool as part of a dedicated plug‐and‐play GI‐BCT reconstruction framework, needed to translate this promising technology to the clinics. [ABSTRACT FROM AUTHOR]

The aim of this work was the evaluation of the software BreastSimulator, a breast x-ray imaging simulation software, as a tool for the creation of 3D uncompressed breast digital models and for the simulation and the optimization of computed tomography (CT) scanners dedicated to the breast. Eight 3D digital breast phantoms were created with glandular fractions in the range 10%-35%. The models are characterised by different sizes and modelled realistic anatomical features. X-ray CT projections were simulated for a dedicated cone-beam CT scanner and reconstructed with the FDK algorithm. X-ray projection images were simulated for 5 mono-energetic (27, 32, 35, 43 and 51 keV) and 3 poly-energetic x-ray spectra typically employed in current CT scanners dedicated to the breast (49, 60, or 80 kVp). Clinical CT images acquired from two different clinical breast CT scanners were used for comparison purposes. The quantitative evaluation included calculation of the power-law exponent, β, from simulated and real breast tomograms, based on the power spectrum fitted with a function of the spatial frequency, f, of the form S(f)  =  α/f   β . The breast models were validated by comparison against clinical breast CT and published data. We found that the calculated β coefficients were close to that of clinical CT data from a dedicated breast CT scanner and reported data in the literature. In evaluating the software package BreastSimulator to generate breast models suitable for use with breast CT imaging, we found that the breast phantoms produced with the software tool can reproduce the anatomical structure of real breasts, as evaluated by calculating the β exponent from the power spectral analysis of simulated images. As such, this research tool might contribute considerably to the further development, testing and optimisation of breast CT imaging techniques.

Conventional mammographic dosimetry has been developed over the past 40 years. Prior to the availability of high-resolution three-dimensional breast images, certain assumptions about breast anatomy were required. These assumptions were based on the information evident on two-dimensional mammograms; they included assumptions of thick skin, a uniform mixture of glandular and adipose tissue, and a median breast density of 50%. Recently, the availability of high-resolution breast CT studies has provided more accurate data about breast anatomy, and this, in turn, has provided the opportunity to update mammographic dosimetry. Based on hundreds of data sets on breast CT volume, a number of studies were performed and reported which have shed light on the basic breast anatomy specific to dosimetry in mammography. It was shown that the average skin thickness of the breast was approximately 1.5 mm, instead of the 4 or 5 mm in the past. In another study, 3-D breast CT data sets were used for validation of the 2-D algorithm developed at the University of Toronto, leading to data suggesting that the overall average breast density is of the order of 16-20%, rather than the previously assumed 50%. Both of these assumptions led to normalized glandular dose (DgN) coefficients which are higher than those of the past. However, a comprehensive study on hundreds of breast CT data sets confirmed the findings of other investigators that there is a more centralized average location of glandular tissue within the breast. Combined with Monte Carlo studies for dosimetry, when accurate models of the distribution of glandular tissue were used, a 30% reduction in the radiation dose (as determined by the DgN coefficient) was found as an average across typical molybdenum and tungsten spectra used clinically. The 30% average reduction was found even when the thinner skin and the lower average breast density were considered. The article reviews three specific anatomic observations made possible based on high-resolution breast CT data by several different research groups. It is noted that, periodically, previous assumptions pertaining to dosimetry can be updated when new information becomes available, so that more accurate dosimetry is achieved. Dogmatic practices typically change slowly, but it is hoped that the medical physics community will continue to evaluate changes in DgN coefficients such that they become more accurate.

PurposeThe purpose of this study is to determine the optimal representative reconstruction and quantitative image feature set for a computer-aided diagnosis (CADx) scheme for dedicated breast computer tomography (bCT).MethodWe used 93 bCT scans that contain 102 breast lesions (62 malignant, 40 benign). Using an iterative image reconstruction (IIR) algorithm, we created 37 reconstructions with different image appearances for each case. In addition, we added a clinical reconstruction for comparison purposes. We used image sharpness, determined by the gradient of gray value in a parenchymal portion of the reconstructed breast, as a surrogate measure of the image qualities/appearances for the 38 reconstructions. After segmentation of the breast lesion, we extracted 23 quantitative image features. Using leave-one-out-cross-validation (LOOCV), we conducted the feature selection, classifier training, and testing. For this study, we used the linear discriminant analysis classifier. Then, we selected the representative reconstruction and feature set for the classifier with the best diagnostic performance among all reconstructions and feature sets. Then, we conducted an observer study with six radiologists using a subset of breast lesions (N = 50). Using 1000 bootstrap samples, we compared the diagnostic performance of the trained classifier to those of the radiologists.ResultThe diagnostic performance of the trained classifier increased as the image sharpness of a given reconstruction increased. Among combinations of reconstructions and quantitative image feature sets, we selected one of the sharp reconstructions and three quantitative image feature sets with the first three highest diagnostic performances under LOOCV as the representative reconstruction and feature set for the classifier. The classifier on the representative reconstruction and feature set achieved better diagnostic performance with an area under the ROC curve (AUC) of 0.94 (95% CI = [0.81, 0.98]) than those of the radiologists, where their maximum AUC was 0.78 (95% CI = [0.63, 0.90]). Moreover, the partial AUC, at 90% sensitivity or higher, of the classifier (pAUC = 0.085 with 95% CI = [0.063, 0.094]) was statistically better (P-value

PurposeCompare conspicuity of ductal carcinoma in-situ (DCIS) to benign calcifications on unenhanced (bCT), contrast-enhanced dedicated breast CT (CEbCT) and mammography (DM).Methods and materialsThe institutional review board approved this HIPAA-compliant study. 42 women with Breast Imaging Reporting and Data System 4 or 5 category micro-calcifications had breast CT before biopsy. Three subjects with invasive disease at surgery were excluded. Two breast radiologists independently compared lesion conspicuity scores (CS) for CEbCT, to bCT and DM. Enhancement was measured in Hounsfield units (HU). Mean CS ± standard deviations are shown. Receiver operating characteristic analysis (ROC) measured radiologists' discrimination performance by comparing CS to enhancement alone. Statistical measurements were made using ANOVA F-test, Wilcoxon rank-sum test and robust linear regression analyses.Results39 lesions (17 DCIS, 22 benign) were analyzed. DCIS (8.5 ± 0.9, n=17) was more conspicuous than benign micro-calcifications (3.6 ± 2.9, n=22; p

PurposeCurrent dosimetry methods in mammography assume that the breast is comprised of a homogeneous mixture of glandular and adipose tissues. Three-dimensional (3D) dedicated breast CT (bCT) data sets were used previously to assess the complex anatomical structure within the breast, characterizing the statistical distribution of glandular tissue in the breast. The purpose of this work was to investigate the effect of bCT-derived heterogeneous glandular distributions on dosimetry in mammography.MethodsbCT-derived breast diameters, volumes, and 3D fibroglandular distributions were used to design realistic compressed breast models comprised of heterogeneous distributions of glandular tissue. The bCT-derived glandular distributions were fit to biGaussian functions and used as probability density maps to assign the density distributions within compressed breast models. The MCNPX 2.6.0 Monte Carlo code was used to estimate monoenergetic normalized mean glandular dose "DgN(E)" values in mammography geometry. The DgN(E) values were then weighted by typical mammography x-ray spectra to determine polyenergetic DgN (pDgN) coefficients for heterogeneous (pDgNhetero) and homogeneous (pDgNhomo) cases. The dependence of estimated pDgN values on phantom size, volumetric glandular fraction (VGF), x-ray technique factors, and location of the heterogeneous glandular distributions was investigated.ResultsThe pDgNhetero coefficients were on average 35.3% (SD, 4.1) and 24.2% (SD, 3.0) lower than the pDgNhomo coefficients for the Mo-Mo and W-Rh x-ray spectra, respectively, across all phantom sizes and VGFs when the glandular distributions were centered within the breast phantom in the coronal plane. At constant breast size, increasing VGF from 7.3% to 19.1% lead to a reduction in pDgNhetero relative to pDgNhomo of 23.6%-27.4% for a W-Rh spectrum. Displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width in the superior and inferior directions, resulted in a 37.3% and a -26.6% change in the pDgNhetero coefficient, respectively, relative to the centered distribution for the Mo-Mo spectrum. Lateral displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width, resulted in a 1.5% change in the pDgNhetero coefficient relative to the centered distribution for the W-Rh spectrum.ConclusionsIntroducing bCT-derived heterogeneous glandular distributions into mammography phantom design resulted in decreased glandular dose relative to the widely used homogeneous assumption. A homogeneous distribution overestimates the amount of glandular tissue near the entrant surface of the breast, where dose deposition is exponentially higher. While these findings are based on clinically measured distributions of glandular tissue using a large cohort of women, future work is required to improve the classification of glandular distributions based on breast size and overall glandular fraction.

Breast CT provides image volumes with isotropic resolution in high contrast, enabling detection of small calcification (down to a few hundred microns in size) and subtle density differences. Since breast is sensitive to x-ray radiation, dose reduction of breast CT is an important topic, and for this purpose, few-view scanning is a main approach. In this article, we propose a Deep Efficient End-to-end Reconstruction (DEER) network for few-view breast CT image reconstruction. The major merits of our network include high dose efficiency, excellent image quality, and low model complexity. By the design, the proposed network can learn the reconstruction process with as few as (9(N) parameters, where N is the side length of an image to be reconstructed, which represents orders of magnitude improvements relative to the state-of-the-art deeplearning-based reconstruction methods that map raw data to tomographic images directly. Also, validated on a cone-beam breast CT dataset prepared by Koning Corporation on a commercial scanner, our method demonstrates a competitive performance over the state-of-the-art reconstruction networks in terms of image quality. The source code of this paper is available at: https://github.com/HuidongXie/DEER.

Purpose: We aim to determine the combination of X-ray spectrum and detector scintillator thickness that maximizes the detectability of microcalcification clusters in dedicated cone-beam breast CT., Approach: A cascaded linear system analysis was implemented in the spatial frequency domain and was used to determine the detectability index using numerical observers for the imaging task of detecting a microcalcification cluster with 0.17 mm diameter calcium carbonate spheres. The analysis considered a thallium-doped cesium iodide scintillator coupled to a complementary metal-oxide semiconductor detector and an analytical filtered-back-projection reconstruction algorithm. Independent system parameters considered were the scintillator thickness, applied X-ray tube voltage, and X-ray beam filtration. The combination of these parameters that maximized the detectability index was considered optimal., Results: Prewhitening, nonprewhitening, and nonprewhitening with eye filter numerical observers indicate that the combination of 0.525 to 0.6 mm thick scintillator, 70 kV, and 0.25 to 0.4 mm added copper filtration maximized the detectability index at a mean glandular dose (MGD) of 4.5 mGy., Conclusion: Using parallel cascade systems' analysis, the combination of parameters that could maximize the detection of microcalcifications was identified. The analysis indicates that a harder beam than that used in current practice may be beneficial for the task of detecting microcalcifications at an MGD suitable for breast cancer screening., (© 2024 Society of Photo-Optical Instrumentation Engineers (SPIE).)

Purpose: The SYRMA‐3D collaboration is setting up a breast computed tomography (bCT) clinical program at the Elettra synchrotron radiation facility in Trieste, Italy. Unlike the few dedicated scanners available at hospitals, synchrotron radiation bCT requires the patient's rotation, which in turn implies a long scan duration (from tens of seconds to few minutes). At the same time, it allows the achievement of high spatial resolution. These features make synchrotron radiation bCT prone to motion artifacts. This article aims at assessing and compensating for motion artifacts through an optical tracking approach. Methods: In this study, patients' movements due to breathing have been first assessed on seven volunteers and then simulated during the CT scans of a breast phantom and a surgical specimen, by adding a periodic oscillatory motion (constant speed, 1 mm amplitude, 12 cycles/minute). CT scans were carried out at 28 keV with a mean glandular dose of 5 mGy. Motion artifacts were evaluated and a correction algorithm based on the optical tracking of fiducial marks was introduced. A quantitative analysis based on the structural similarity (SSIM) index and the normalized mean square error (nMSE) was performed on the reconstructed CT images. Results: CT images reconstructed through the optical tracking procedure were found to be as good as the motionless reference image. Moreover, the analysis of SSIM and nMSE demonstrated that an uncorrected motion of the order of the system's point spread function (around 0.1 mm in the present case) can be tolerated. Conclusions: Results suggest that a motion correction procedure based on an optical tracking system would be beneficial in synchrotron radiation bCT. [ABSTRACT FROM AUTHOR]

PURPOSE:The purpose of this study is to measure the effectiveness of local curvature measures as novel image features for classifying breast tumors. METHODS:A total of 119 breast lesions from 104 noncontrast dedicated breast computed tomography images of women were used in this study. Volumetric segmentation was done using a seed-based segmentation algorithm and then a triangulated surface was extracted from the resulting segmentation. Total, mean, and Gaussian curvatures were then computed. Normalized curvatures were used as classification features. In addition, traditional image features were also extracted and a forward feature selection scheme was used to select the optimal feature set. Logistic regression was used as a classifier and leave-one-out cross-validation was utilized to evaluate the classification performances of the features. The area under the receiver operating characteristic curve (AUC, area under curve) was used as a figure of merit. RESULTS:Among curvature measures, the normalized total curvature (CT) showed the best classification performance (AUC of 0.74), while the others showed no classification power individually. Five traditional image features (two shape, two margin, and one texture descriptors) were selected via the feature selection scheme and its resulting classifier achieved an AUC of 0.83. Among those five features, the radial gradient index (RGI), which is a margin descriptor, showed the best classification performance (AUC of 0.73). A classifier combining RGI and CT yielded an AUC of 0.81, which showed similar performance (i.e., no statistically significant difference) to the classifier with the above five traditional image features. Additional comparisons in AUC values between classifiers using different combinations of traditional image features and CT were conducted. The results showed that CT was able to replace the other four image features for the classification task. CONCLUSIONS:The normalized curvature measure contains useful information in classifying breast tumors. Using this, one can reduce the number of features in a classifier, which may result in more robust classifiers for different datasets.

Purpose: A clinical‐prototype, dedicated, cone‐beam breast computed tomography (CBBCT) system with offset detector is undergoing clinical evaluation at our institution. This study is to estimate the normalized glandular dose coefficients (DgNCT) that provide air kerma‐to‐mean glandular dose conversion factors using Monte Carlo simulations. Materials and methods: The clinical prototype CBBCT system uses 49 kV x‐ray spectrum with 1.39 mm 1st half‐value layer thickness. Monte Carlo simulations (GATE, version 8) were performed with semi‐ellipsoidal, homogeneous breasts of various fibroglandular weight fractions (fg=0.01,0.15,0.5,1), chest wall diameters (d=8,10,14,18,20 cm), and chest wall to nipple length (l=0.75d), aligned with the axis of rotation (AOR) located at 65 cm from the focal spot to determine the DgNCT. Three geometries were considered – 40×30‐cm detector with no offset that served as reference and corresponds to a clinical CBBCT system, 30×30‐cm detector with 5 cm offset, and a 30×30‐cm detector with 10 cm offset. Results: For 5 cm lateral offset, the DgNCT ranged 0.177‐0.574 mGy/mGy and reduction in DgNCT with respect to reference geometry was observed only for 18 cm (6.4%±0.23%) and 20 cm (9.6%±0.22%) diameter breasts. For the 10 cm lateral offset, the DgNCT ranged 0.221‐0.581 mGy/mGy and reduction in DgNCT was observed for all breast diameters. The reduction in DgNCT was 1.4%±0.48%, 7.1%±0.13%, 17.5%±0.19%, 25.1%±0.15%, and 27.7%±0.08% for 8, 10, 14, 18, and 20 cm diameter breasts, respectively. For a given breast diameter, the reduction in DgNCT with offset‐detector geometries was not dependent on fg. Numerical fits of DgNCTd,l,fg were generated for each geometry. Conclusion: The DgNCT and the numerical fit, DgNCTd,l,fg would be of benefit for current CBBCT systems using the reference geometry and for future generations using offset‐detector geometry. There exists a potential for radiation dose reduction with offset‐detector geometry, provided the same technique factors as the reference geometry are used, and the image quality is clinically acceptable. [ABSTRACT FROM AUTHOR]

The SYRMA-3D collaboration is setting up the first clinical trial of phase-contrast breast CT with synchrotron radiation at the Elettra synchrotron facility in Trieste, Italy. In this communication, a quality control protocol for breast CT is proposed, and a first test of image quality measurements is performed by means of a custom-made radiographic phantom.

Purpose: Contrast‐enhanced imaging of the breast is frequently used in breast MRI and has recently become more common in mammography. The purpose of this study was to make single‐scan contrast‐enhanced imaging feasible for photon‐counting breast CT (pcBCT) and to assess the spectral performance of a pcBCT scanner by evaluating iodine maps and virtual non‐contrast (VNC) images. Methods: We optimized the settings of a pcBCT to maximize the signal‐to‐noise ratio between iodinated contrast agent and breast tissue. Therefore, an electronic energy threshold dividing the x‐ray spectrum used into two energy bins was swept from 23.17 keV to 50.65 keV. Validation measurements were performed by placing syringes with contrast agent (2.5 mg/ml to 40 mg/ml) in phantoms with 7.5 cm and 12 cm in diameter. Images were acquired at different tube currents and reconstructed with 300 μm isotropic voxel size. Iodine maps and VNC images were generated using image‐based material decomposition. Iodine concentrations and CT values were measured for each syringe and compared to the known concentrations and reference CT values. Results: Maximal signal‐to‐noise ratios were found at a threshold position of 32.59 keV. Accurate iodine quantification (average root mean square error of 0.56 mg/ml) was possible down to a concentration of 2.5 mg/ml for all tube currents investigated. The enhancement has been sufficiently removed in the VNC images, so they can be interpreted as unenhanced CT images. Only minor changes of CT values compared to a conventional CT scan were observed. Noise was increased by the decomposition by a factor of 2.62 and 4.87 (7.5 cm and 12 cm phantoms) but did not compromise the accuracy of the iodine quantification. Conclusions: Accurate iodine quantification and generation of VNC images can be achieved using contrast‐enhanced pcBCT from a single CT scan in the absence of temporal or spatial misalignment. Using iodine maps and VNC images, pcBCT has the potential to reduce dose, shorten examination and reading time, and to increase cancer detection rates. [ABSTRACT FROM AUTHOR]

The goal of this work is to determine the optimal projection to detect breast masses on Chest CT. Institutional Review Board (HIPPA compliant) approval was obtained with a waiver of consent. 10 image pairs of Chest CT images containing breast masses were selected for review by 10 chest radiologists: the pairs consisted of axial, sagittal, coronal and axial MIP images (MIP images) with each projection compared to a MIP and with one another. For each pair, the image where the mass was most conspicuous was recorded. MIPs were preferred to any cross sectional projection 82% of the time; sagittal (63%) or coronal (63%) images were preferred to the axial projection. When sagittal and coronal images were compared there was no preference. MIP images should be obtained and reviewed for breast pathology; sagittal or coronal projections may provide additional information. • Denser breast mass is more conspicuous on the MIP images. • Radiologists preferred the MIP images for breast mass detection. • Sagittal and coronal images were found to be equally preferable to axial. [ABSTRACT FROM AUTHOR]

BACKGROUND: High-resolution, low-noise detectors with minimal dead-space at chest-wall could improve posterior coverage and microcalcification visibility in the dedicated cone-beam breast CT (CBBCT). However, the smaller field-of-view necessitates laterally-shifted detector geometry to enable optimizing the air-gap for x-ray scatter rejection. OBJECTIVE: To evaluate laterally-shifted detector geometry for CBBCT with clinical projection datasets that provide for anatomical structures and lesions. METHODS: CBBCT projection datasets (n = 17 breasts) acquired with a 40×30 cm detector (1024×768-pixels, 0.388-mm pixels) were truncated along the fan-angle to emulate 20.3×30 cm, 22.2×30 cm and 24.1×30 cm detector formats and correspond to 20, 120, 220 pixels overlap in conjugate views, respectively. Feldkamp-Davis-Kress (FDK) algorithm with 3 different weighting schemes were used for reconstruction. Visual analysis for artifacts and quantitative analysis of root-mean-squared-error (RMSE), absolute difference between truncated and 40×30 cm reconstructions (Diff), and its power spectrum (PSDiff) were performed. RESULTS: Artifacts were observed for 20.3×30 cm, but not for other formats. The 24.1×30 cm provided the best quantitative results with RMSE and Diff (both in units of μ, cm-1) of 4.39×10-3±1.98×10-3 and 4.95×10-4±1.34×10-4, respectively. The PSDiff (>0.3 cycles/mm) was in the order of 10-14μ2mm3 and was spatial-frequency independent. CONCLUSIONS: Laterally-shifted detector CBBCT with at least 220 pixels overlap in conjugate views (24.1×30 cm detector format) provides quantitatively accurate and artifact-free image reconstruction. [ABSTRACT FROM AUTHOR]

This work demonstrates the potential for using a deformable mapping method to register lesions between dedicated breast computed tomography (bCT) and both automated breast ultrasound (ABUS) and digital breast tomosynthesis (DBT) images (craniocaudal [CC] and mediolateral oblique [MLO] views). Two multi-modality breast phantoms with external fiducial markers attached were imaged by the three modalities. The DBT MLO view was excluded for the second phantom. The automated deformable mapping algorithm uses biomechanical modeling to determine corresponding lesions based on distances between their centers of mass (dCOM) in the deformed bCT model and the reference model (DBT or ABUS). For bCT to ABUS, the mean dCOM was 5.2 ± 2.6 mm. For bCT to DBT (CC), the mean dCOM was 5.1 ± 2.4 mm. For bCT to DBT (MLO), the mean dCOM was 4.7 ± 2.5 mm. This application could help improve a radiologist's efficiency and accuracy in breast lesion characterization, using multiple imaging modalities. [ABSTRACT FROM AUTHOR]

Dedicated cone-beam breast computed tomography (CBBCT) is an emerging modality and provides fully three-dimensional (3D) images of the uncompressed breast at an isotropic voxel resolution. In an effort to translate this modality to breast cancer screening, advanced image reconstruction methods are being pursued. Since radiographic breast density is an established risk factor for breast cancer and CBBCT provides volumetric data, this study investigates the reproducibility of the volumetric glandular fraction (VGF), defined as the proportion of fibroglandular tissue volume relative to the total breast volume excluding the skin. Four image reconstruction methods were investigated: the analytical Feldkamp-Davis-Kress (FDK), a compressed sensing-based fast, regularized, iterative statistical technique (FRIST), a fully supervised deep learning approach using a multi-scale residual dense network (MS-RDN), and a self-supervised approach based on Noise-to-Noise (N2N) learning. Projection datasets from 106 women who participated in a prior clinical trial were reconstructed using each of these algorithms at a fixed isotropic voxel size of (0.273 mm 3 ). Each reconstructed breast volume was segmented into skin, adipose, and fibroglandular tissues, and the VGF was computed. The VGF did not differ among the four reconstruction methods ( p = 0.167), and none of the three advanced image reconstruction algorithms differed from the standard FDK reconstruction ( p > 0.862). Advanced reconstruction algorithms developed for low-dose CBBCT reproduce the VGF to provide quantitative breast density, which can be used for risk estimation.

Breast computed tomography (BCT) is an emerging application of X‐ray tomography in radiological practice. A few clinical prototypes are under evaluation in hospitals and new systems are under development aiming at improving spatial and contrast resolution and reducing delivered dose. At the same time, synchrotron‐radiation phase‐contrast mammography has been demonstrated to offer substantial advantages when compared with conventional mammography. At Elettra, the Italian synchrotron radiation facility, a clinical program of phase‐contrast BCT based on the free‐space propagation approach is under development. In this paper, full‐volume breast samples imaged with a beam energy of 32 keV delivering a mean glandular dose of 5 mGy are presented. The whole acquisition setup mimics a clinical study in order to evaluate its feasibility in terms of acquisition time and image quality. Acquisitions are performed using a high‐resolution CdTe photon‐counting detector and the projection data are processed via a phase‐retrieval algorithm. Tomographic reconstructions are compared with conventional mammographic images acquired prior to surgery and with histologic examinations. Results indicate that BCT with monochromatic beam and free‐space propagation phase‐contrast imaging provide relevant three‐dimensional insights of breast morphology at clinically acceptable doses and scan times. [ABSTRACT FROM AUTHOR]

Purpose: The purpose of this work was to generate uncompressed heterogeneous breast phantom models using size‐dependent fibroglandular distributions derived from a large cohort of breast CT (bCT) datasets, and to compare differences in normalized glandular dose coefficients for bCT "DgNCT" when the realistic heterogeneous model is considered relative to the simple, homogeneous model used in the past. Methods: A cohort of 274 segmented bCT datasets were used to quantify the fibroglandular tissue distribution within the breast parenchyma. Each dataset was interpolated to an isotropic voxel size of 0.25 mm and the breast center‐of‐mass was aligned for all coronal slices. Each aligned dataset was converted into two binarized volumes representing voxels containing only glandular tissue "G(x,y,z)" and voxels containing glandular or adipose tissue "AG(x,y,z)". The datasets were classified by volume in accordance with previously reported size‐dependent, breast‐shaped phantoms. Within the five groups — each containing on average 55 datasets, all G(x,y,z) and AG(x,y,z) volumes were summed separately representing the cumulative distribution of glandular tissue or breast parenchyma (glandular + adipose), respectively. G(x,y,z) was divided by AG(x,y,z) on a voxel‐by‐voxel basis resulting in a glandular fraction "GF(x,y,z)" distribution for each phantom size. The GF(x,y,z) distributions were used to construct heterogeneous mathematical phantoms for the small, median, and large breast sizes with a 1.5 mm skin thickness — based on previously reported measurements from bCT, and a 5 mm skin thickness for comparison with outdated assumptions of skin thickness. A subset of 15 bCT datasets from the cohort (five for each breast size) were used to construct voxelized patient models for validation of the heterogeneous phantom models. Monte Carlo techniques were used to estimate monoenergetic DgN(E)CT values for photon energies from 9 to 70 keV (in 1 keV intervals) using the mathematical phantoms composed of either heterogeneous or homogeneous breast parenchyma. Polyenergetic (pDgNCT) coefficients were determined by weighting the DgN(E)CT values by x‐ray spectra tuned to the beam characteristic of breast CT. Dose coefficients were compared between the two breast compositions for each volume class, breast density, and skin thickness. Results: For photon energies ≲45 keV, the homogeneous model overestimates DgN(E) values relative to the realistic heterogeneous model sorted into five volume classes. The 5 mm skin thickness underestimates DgN(E) values relative to the realistic 1.5 mm thickness for lower energies and the differences diminish up to 70 keV. Averaged across all phantom sizes the homogeneous model overestimates pDgNCT by 5.7% and 23.3% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively. The heterogeneous model was also found to be in agreement with the voxelized bCT patient models with pDgNCT differences less than 2.3% and 5.2% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively, across all phantom sizes. Conclusion: Anatomically accurate heterogeneous phantom models were developed using bCT image‐derived fibroglandular tissue distributions. These new models improve the accuracy of bCT dosimetry, and in conjunction with previous models for mammography, may help in providing a more universally accepted breast dosimetry model. [ABSTRACT FROM AUTHOR]

Highlights • Modeling a breast compressed in MLO view instead of in CC view does not produce any significant differences in MGD estimates. • Skin model presents a noticeable influence on T-factor in DBT. • Correction factor for taking into account the ion-chamber model in MGD estimates in BCT was calculated. Abstract Purpose To provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT). Methods The MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements. Results In mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT. Conclusions A series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models. [ABSTRACT FROM AUTHOR]

• Phase-contrast CT addresses poor soft-tissue contrast in X-ray breast imaging. • Interactive 3D breast virtual visualizations available at different resolutions. • Spatial resolution bridging from radiological level to virtual histology level. • 2D and 3D images may be combinedly used in a multidisciplinary clinical context. • Discernibility of breast structures and microcalcifications is enhanced. Phase-contrast imaging techniques address the issue of poor soft-tissue contrast encountered in traditional X-ray imaging. This can be accomplished with the propagation-based phase-contrast technique by employing a coherent photon beam, which is available at synchrotron facilities, as well as long sample-to-detector distances. This study demonstrates the optimization of propagation-based phase-contrast computed tomography (CT) techniques for multiscale X-ray imaging of the breast at the Elettra synchrotron facility (Trieste, Italy). Two whole breast mastectomy samples were acquired with propagation-based breast-CT using a monochromatic synchrotron beam at a pixel size of 60 µm. Paraffin-embedded blocks sampled from the same tissues were scanned with propagation-based micro-CT imaging using a polychromatic synchrotron beam at a pixel size of 4 µm. Images of both methodologies and of the same sample were spatially registered. The resulting images showed the transition from whole breast imaging with propagation-based breast-CT methodology to virtual histology with propagation-based micro-CT imaging of the same sample. Additionally, conventional histological images were matched to virtual histology images. Phase-contrast images offer a high resolution with low noise, which allows for a highly precise match between virtual and conventional histology. Furthermore, those techniques allow a clear discernment of breast structures, lesions, and microcalcifications, being a promising clinically-compatible tool for breast imaging in a multiscale approach, to either assist in the detection of cancer in full volume breast samples or to complement structure identification in paraffin-embedded breast tissue samples. [ABSTRACT FROM AUTHOR]

Objective. This paper proposes a new objective function to improve the quality of synthesized breast CT images generated by the GAN and compares the GAN performances on transfer learning datasets from different image domains. Approach. The proposed objective function, named beta loss function, is based on the fact that x-ray-based breast images follow the power-law spectrum. Accordingly, the exponent of the power-law spectrum (beta value) for breast CT images is approximately two. The beta loss function is defined in terms of L1 distance between the beta value of synthetic images and validation samples. To compare the GAN performances for transfer learning datasets from different image domains, ImageNet and anatomical noise images are used in the transfer learning dataset. We employ styleGAN2 as the backbone network and add the proposed beta loss function. The patient-derived breast CT dataset is used as the training and validation dataset; 7355 and 212 images are used for network training and validation, respectively. We use the beta value evaluation and Fréchet inception distance (FID) score for quantitative evaluation. Main results. For qualitative assessment, we attempt to replicate the images from the validation dataset using the trained GAN. Our results show that the proposed beta loss function achieves a more similar beta value to real images and a lower FID score. Moreover, we observe that the GAN pretrained with anatomical noise images achieves better equality than ImageNet for beta value evaluation and FID score. Finally, the beta loss function with anatomical noise as the transfer learning dataset achieves the lowest FID score. Significance. Overall, the GAN using the proposed beta loss function with anatomical noise images as the transfer learning dataset provides the lowest FID score among all tested cases. Hence, this work has implications for developing GAN-based breast image synthesis methods for medical imaging applications., (© 2023 Institute of Physics and Engineering in Medicine.)

Purpose: The purpose of this study was to develop a realistic patient‐based 4D digital breast phantom including time‐varying contrast enhancement for simulation of dedicated breast CT perfusion imaging. Methods: A 3D static phantom is first created by segmenting a breast CT image from a healthy patient into skin, fibroglandular tissue, adipose tissue, and vasculature. For the creation of abnormal cases, a breast lesion model was developed and can be added to the phantom. After defining the necessary perfusion parameters for each tissue (e.g., arterial input function for vasculature, blood volume and blood flow for the other normal tissues) based on contrast‐enhanced dynamic breast MRI data, the corresponding time‐enhancement curves are computed for each voxel in the phantom, according to tissue type. These curves are calculated by convolution between the arterial input function and a shifted exponential function. This exponential depends on the perfusion parameters associated with each tissue voxel, and, to incorporate normal biological variability, a uniform random distribution is used to vary the perfusion parameters on a voxel‐basis. Finally, a 4D array is produced by sampling the continuous time‐enhancement curves at the desired sampling rate. Beside modeling different enhancement dynamics according to the given input perfusion parameters, the phantom also includes the possibility to realistically simulate different spatial enhancement patterns for the breast parenchyma, taking into account the arterial sources supplying the breast. Finally, different patterns of contrast medium uptake can also be simulated for the tumor models (homogeneous and rim enhancement). Results: As an example, a typical 4D phantom has dimensions of 426 × 421 × 260 × 559 (x, y, z, t), with a voxel size of 273 μm and a sampling time of 1 s. The characteristics of the tumor model can be modified at will to evaluate perfusion in different types of breast lesions. Results show the expected enhancement of tissues, consistent with the given input parameters. Moreover, the tumor models evaluated in this work show different enhancement dynamics according to the tumor type (defined by different input perfusion parameters), and also present a higher enhancement compared to the other healthy tissues, as expected. Conclusions: The proposed digital phantom can model the breast tissue perfusion during 4D breast CT image acquisition, displaying the different enhancement dynamics that could be found in a real patient breast. This phantom can be used during the development of dynamic contrast‐enhanced dedicated breast CT imaging, for optimization of image acquisition, image reconstruction, and image analysis. This modality could provide functional information of the breast, resulting in detection, diagnosis, and treatment improvements of breast cancer with breast CT. [ABSTRACT FROM AUTHOR]

Purpose: With the advent of three‐dimensional (3D) breast imaging modalities such as digital breast tomosynthesis (DBT) and dedicated breast CT (bCT), research into new anthropomorphic breast phantoms has accelerated. These breast phantoms are important for the optimization of new breast imaging systems, assessing new regulatory submissions to prove safety and effectiveness, and for developing new approaches to acceptance and constancy testing of 3D breast imaging systems. This paper provides a review of current research investigating both digital and physical breast phantom development for use in x‐ray based imaging. Methods: Two approaches for designing anthropomorphic, digital breast phantoms are discussed, procedural model‐based phantom generation, where breast features are expressed using mathematical models, and patient‐based generation, where breast structures from tissue specimens or patient‐based breast MR or CT volumes are segmented. Following this discussion, a review of physical anthropomorphic phantoms is given, with emphasis on the advantages and disadvantages present with each approach. Conclusions: This paper provides a summary of the state‐of‐the‐art in anthropomorphic breast phantom development for x‐ray breast imaging. The primary advantage of model‐based digital phantoms is that an unlimited number of phantoms with varying size, shape, and density can be generated. Current research on model‐based breast phantoms is producing more and more realistic breast models; however, they probably are not yet able to pass the so‐called "fool the radiologist" visualization test. Empirical patient‐based breast phantoms are typically based on clinical breast CT data and look more realistic. However, clinical breast CT images have limited spatial resolution and thus do not always portray the finer details in the breast. A number of innovative solutions have been proposed for fabricating physical anthropomorphic breast phantoms based on digital phantom models; however, a number of challenges remain, including realistic modeling of x‐ray attenuation properties and accurately representing high‐frequency structures within breast. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Purpose In cone-beam computed tomography dedicated to the breast (BCT), the mean glandular dose (MGD) is the dose metric of reference, evaluated from the measured air kerma by means of normalized glandular dose coefficients (DgN CT ). This work aimed at computing, for a simple breast model, a set of DgN CT values for monoenergetic and polyenergetic X-ray beams, and at validating the results vs. those for patient specific digital phantoms from BCT scans. Methods We developed a Monte Carlo code for calculation of monoenergetic DgN CT coefficients (energy range 4.25–82.25 keV). The pendant breast was modelled as a cylinder of a homogeneous mixture of adipose and glandular tissue with glandular fractions by mass of 0.1%, 14.3%, 25%, 50% or 100%, enveloped by a 1.45 mm-thick skin layer. The breast diameter ranged between 8 cm and 18 cm. Then, polyenergetic DgN CT coefficients were analytically derived for 49-kVp W-anode spectra (half value layer 1.25–1.50 mm Al), as in a commercial BCT scanner. We compared the homogeneous models to 20 digital phantoms produced from classified 3D breast images. Results Polyenergetic DgN CT resulted 13% lower than most recent published data. The comparison vs. patient specific breast phantoms showed that the homogeneous cylindrical model leads to a DgN CT percentage difference between −15% and +27%, with an average overestimation of 8%. Conclusions A dataset of monoenergetic and polyenergetic DgN CT coefficients for BCT was provided. Patient specific breast models showed a different volume distribution of glandular dose and determined a DgN CT 8% lower, on average, than homogeneous breast model. [ABSTRACT FROM AUTHOR]

Rationale and Objectives: Dedicated breast computed tomography (DBCT) is an emerging and promising modality for breast lesions. The objective of this study was to evaluate the potential use of applying the BI-RADS Mammography Atlas 5th Edition for reporting and assessing breast lesions on DBCT. Currently, no atlas exists for DBCT.Materials and Methods: Four radiologists trained in breast imaging were recruited in this institutional review board-approved, Health Insurance Portability and Accountability Act-compliant study. The enrolled radiologists, who were blinded to mammographic and histopathologic findings, individually reviewed 30 randomized DBCT cases that contained marked lesions. Thirty-four lesions were included in this study: 24 (70.6%) masses, 7 (20.6%) calcifications, and 3 (8.8%) architectural distortions. Eight (23.5%) lesions were malignant and 26 (76.5%) were benign. The reader was asked to specify according to the BI-RADS Mammography Atlas for each marked DBCT lesion: primary findings, features, breast density, and final assessment. We calculated readers' diagnostic performances for differentiating between benign and malignant lesions and interobserver variability for reporting and assessing lesions using a generalized estimating equation and the Fleiss kappa (κ) statistic.Results: The estimated overall sensitivity of the readers was 0.969, and the specificity was 0.529. There were no significant differences in the sensitivity and the specificity between lesion types. For reporting the presence of a primary finding, the overall substantial agreement (κ = 0.70) was seen. In assigning the breast density and the final assessment, the overall agreement was moderate (κ = 0.53) and fair (κ = 0.30).Conclusion: The use of the BI-RADS Mammography Atlas 5th Edition for DBCT showed high performance and good agreement among readers. [ABSTRACT FROM AUTHOR]