X-ray characterisation of normal and neoplastic breast tissues

X-ray source motion blur modeling and deblurring with generative diffusion for digital breast tomosynthesis

Objective. Digital breast tomosynthesis (DBT) has significantly improved the diagnosis of breast cancer due to its high sensitivity and specificity in detecting breast lesions compared to two-dimensional mammography. However, one of the primary challenges in DBT is the image blur resulting from x-ray source motion, particularly in DBT systems with a source in continuous-motion mode. This motion-induced blur can degrade the spatial resolution of DBT images, potentially affecting the visibility of subtle lesions such as microcalcifications.Approach. We addressed this issue by deriving an analytical in-plane source blur kernel for DBT images based on imaging geometry and proposing a post-processing image deblurring method with a generative diffusion model as an image prior.Main results. We showed that the source blur could be approximated by a shift-invariant kernel over the DBT slice at a given height above the detector, and we validated the accuracy of our blur kernel modeling through simulation. We also demonstrated the ability of the diffusion model to generate realistic DBT images. The proposed deblurring method successfully enhanced spatial resolution when applied to DBT images reconstructed with detector blur and correlated noise modeling.Significance. Our study demonstrated the advantages of modeling the imaging system components such as source motion blur for improving DBT image quality.

Knowledge of Non-Invasive Biomarkers of Breast Cancer, Risk Factors, and BSE Practices Among Nursing Undergraduates in Farasan Island, KSA

Background of the Study

Mammograms are sometimes met with issues of overdiagnosis and underdiagnosis; therefore, they are less reliable in identifying cancer in women with dense breasts. As a result, it is critical to be aware of other sensitive screening techniques for the early diagnosis of breast cancer.

Aim

The ultimate objective of this study was to assess the knowledge of nursing undergraduates regarding non-invasive biomarkers, such as volatile organic compounds in breath, nipple aspirate fluid, sweat, urine, and tears, for the early detection of breast cancer to help improve patient care, determine the risk factors, and encourage practice of breast self-examination.

Methods

Cross-sectional research was done in the Department of Nursing at Farasan campus using a self-structured questionnaire as the study tool. A total of 260 students willingly participated. The study tool had evaluation questions focused on the non-invasive biomarkers of breast cancer, risk factors, and breast self-examination practices to collect data. The data were subjected to descriptive and inferential statistics. The statistical significance was calculated at P < .05. Data analyses were done using Microsoft Excel (2013).

Results

A significant knowledge gap existed among the study participants about the non-invasive biomarkers of breast cancer. A lesser percentage of students (25%) stated that they do breast self-examination on a monthly basis. The most common reasons for not doing the breast self-examination were "not knowing how to do the breast self-examination" (77.3%), fear of a positive diagnosis (53.9%), thinking that they are not at risk as all were in their teens and hence not required (44.7%), and lack of time (48.7%). Age and frequency of breast self-examination were significantly associated (P < .05) as those few students (22.7%) who were doing breast self-examination practices every 2-4 months belonged to a higher study year. Furthermore, knowledge regarding incidence rates and health care expenditure by the government on breast cancer was also significantly low (P < .05).

Conclusions

Outcomes would help prioritize actions to help future nurses better understand breast cancer, allowing them to extend patient care in the best way possible.

Power-law spectrum-based objective function to train a generative adversarial network with transfer learning for the synthetic breast CT image

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.

Convolutional neural network-based model observer for signal known statistically task in breast tomosynthesis images

BACKGROUND

Since human observer studies are resource-intensive, mathematical model observers are frequently used to assess task-based image quality. The most common implementation of these model observers assume that the signal information is exactly known. However, these tasks cannot thoroughly represent situations where the signal information is not exactly known in terms of size and shape.

PURPOSE

Considering the limitations of the tasks for which signal information is exactly known, we proposed a convolutional neural network (CNN)-based model observer for signal known statistically (SKS) and background known statistically (BKS) detection tasks in breast tomosynthesis images.

METHODS

A wide parameter search was conducted from six different acquisition angles (i.e., 10°, 20°, 30°, 40°, 50°, and 60°) within the same dose level (i.e., 2.3 mGy) under two separate acquisition schemes: (1) constant total number of projections, and (2) constant angular separation between projections. Two different types of signals: spherical (i.e., SKE tasks) and spiculated (i.e., SKS tasks) were used. The detection performance of the CNN-based model observer was compared with that of the Hotelling observer (HO) instead of the IO. Pixel-wise gradient-weighted class activation mapping (pGrad-CAM) map was extracted from each reconstructed tomosynthesis image to provide an intuitive understanding of the trained CNN-based model observer.

RESULTS

The CNN-based model observer achieved a higher detection performance compared to that of the HO for all tasks. Moreover, the improvement in its detection performance was greater for SKS tasks compared to that for SKE tasks. These results demonstrated that the addition of nonlinearity improved the detection performance owing to the variation of the background and signal. Interestingly, the pGrad-CAM results effectively localized the class-specific discriminative region, further supporting the quantitative evaluation results of the CNN-based model observer. In addition, we verified that the CNN-based model observer required fewer images to achieve the detection performance of the HO.

CONCLUSIONS

In this work, we proposed a CNN-based model observer for SKS and BKS detection tasks in breast tomosynthesis images. Throughout the study, we demonstrated that the detection performance of the proposed CNN-based model observer was superior to that of the HO.

An in silico study on the detectability of field cancerization through parenchymal analysis of digital mammograms

BACKGROUND

Parenchymal analysis has shown promising performance for the assessment of breast cancer risk through the characterization of the texture features of mammography images. However, the working principles behind this practice are yet not well understood. Field cancerization is a phenomenon associated with genetic and epigenetic alterations in large volumes of cells, putting them on a path of malignancy before the appearance of recognizable cancer signs. Evidence suggests that it can induce changes in the biochemical and optical properties of the tissue.

PURPOSE

The aim of this work was to study whether the extended genetic mutations and epigenetic changes due to field cancerization, and the impact they have on the biochemistry of breast tissues are detectable in the radiological patterns of mammography images.

METHODS

An in silico experiment was designed, which implied the development of a field cancerization model to modify the optical tissue properties of a cohort of 60 voxelized virtual breast phantoms. Mammography images from these phantoms were generated and compared with images obtained from their non-modified counterparts, that is, without field cancerization. We extracted 33 texture features from the breast area to quantitatively assess the impact of the field cancerization model. We analyzed the similarity and statistical equivalence of texture features with and without field cancerization using the t-test, Wilcoxon sign rank test and Kolmogorov-Smirnov test, and performed a discrimination test using multinomial logistic regression analysis with lasso regularization.

RESULTS

With modifications of the optical tissue properties on 3.9% of the breast volume, some texture features started to fail to show equivalence (p < 0.05). At 7.9% volume modification, a high percent of texture features showed statistically significant differences (p < 0.05) and non-equivalence. At this level, multinomial logistic regression analysis of texture features showed a statistically significant performance in the discrimination of mammograms from breasts with and without field cancerization (AUC = 0.89, 95% CI: 0.75-1.00).

CONCLUSIONS

These results support the idea that field cancerization is a feasible underlying working principle behind the distinctive performance of parenchymal analysis in breast cancer risk assessment.

Technical note: Realization and uncertainty analysis for an adjustable 3D structured breast phantom in digital breast tomosynthesis

BACKGROUND

Projection imaging phantoms are often optimized for 2-dimensional image characteristics in homogeneous backgrounds. Therefore, evaluation of image quality in tomosynthesis (DBT) lacks accepted and established phantoms.

PURPOSE

We describe a 3D breast phantom with a structured, variable background. The phantom is an adaptable and advanced version of the L1 phantom by Cockmartin et al. Phantom design and its use for quality assurance measurements for DBT devices are described. Four phantoms were compared to assess the objectivity.

METHODS

The container size was increased to a diameter of 24 cm and a total height of 53.5 mm. Spiculated masses were replaced by five additional non-spiculated masses for higher granularity in threshold diameter resolution. These patterns are adjustable to the imaging device. The masses were printed in one session with a base layer using two-component 3D printing. New materials compared to the L1 phantom improved the attenuation difference between the lesion models and the background. Four phantoms were built and intra-human observer, inter-human observer and inter-phantom variations were determined. The latter assess the reproducibility of the phantom production. Coefficients of variance (V) were calculated for all three variations.

RESULTS

The difference of the attenuation coefficients between the lesion models and the background was 0.20 cm-1 (with W/Al at 32 kV, equivalent to 19-20 keV effective energy) compared to 0.21 cm-1 for 50/50 glandular/adipose breast tissue and cancerous lesions. PMMA equivalent thickness of the phantom was 47.0 mm for the Siemens Mammomat Revelation. For the masses, theV i n t r a $V_{intra}$ for the intra-observer variation was 0.248, the averaged inter-observer variation,V ¯ i n t e r $\overline{V}_{inter}$ was 0.383.V p h a n t o m $V_{phantom}$ for phantom variance was 0.321. For the micro-calcifications,V i n t r a $V_{intra}$ was 0.0429,V ¯ i n t e r = $\overline{V}_{inter}=$ 0.0731 andV p h a n t o m = $V_{phantom}=$ 0.0759.

CONCLUSIONS

Position, orientation and shape of the masses are reproducible and attenuation differences appropriate. The phantom presented proved to be a candidate test object for quality control.

Breast Density Evaluation According to BI-RADS 5th Edition on Digital Breast Tomosynthesis: AI Automated Assessment Versus Human Visual Assessment

Background: The assessment of breast density is one of the main goals of radiologists because the masking effect of dense fibroglandular tissue may affect the mammographic identification of lesions. The BI-RADS 5th Edition has revised the mammographic breast density categories, focusing on a qualitative evaluation rather than a quantitative one. Our purpose is to compare the concordance of the automatic classification of breast density with the visual assessment according to the latest available classification. Methods: A sample of 1075 digital breast tomosynthesis images from women aged between 40 and 86 years (58 ± 7.1) was retrospectively analyzed by three independent readers according to the BI-RADS 5th Edition. Automated breast density assessment was performed on digital breast tomosynthesis images with the Quantra software version 2.2.3. Interobserver agreement was assessed with kappa statistics. The distributions of breast density categories were compared and correlated with age. Results: The agreement on breast density categories was substantial to almost perfect between radiologists (κ = 0.63–0.83), moderate to substantial between radiologists and the Quantra software (κ = 0.44–0.78), and the consensus of radiologists and the Quantra software (κ = 0.60–0.77). Comparing the assessment for dense and non-dense breasts, the agreement was almost perfect in the screening age range without a statistically significant difference between concordant and discordant cases when compared by age. Conclusions: The categorization proposed by the Quantra software has shown a good agreement with the radiological evaluations, even though it did not completely reflect the visual assessment. Thus, clinical decisions regarding supplemental screening should be based on the radiologist’s perceived masking effect rather than the data produced exclusively by the Quantra software.

Beyond Breast Density: Risk Measures for Breast Cancer in Multiple Imaging Modalities

Breast density is an independent risk factor for breast cancer. In digital mammography and digital breast tomosynthesis, breast density is assessed visually using the four-category scale developed by the American College of Radiology Breast Imaging Reporting and Data System (5th edition as of November 2022). Epidemiologically based risk models, such as the Tyrer-Cuzick model (version 8), demonstrate superior modeling performance when mammographic density is incorporated. Beyond just density, a separate mammographic measure of breast cancer risk is parenchymal textural complexity. With advancements in radiomics and deep learning, mammographic textural patterns can be assessed quantitatively and incorporated into risk models. Other supplemental screening modalities, such as breast US and MRI, offer independent risk measures complementary to those derived from mammography. Breast US allows the two components of fibroglandular tissue (stromal and glandular) to be visualized separately in a manner that is not possible with mammography. A higher glandular component at screening breast US is associated with higher risk. With MRI, a higher background parenchymal enhancement of the fibroglandular tissue has also emerged as an imaging marker for risk assessment. Imaging markers observed at mammography, US, and MRI are powerful tools in refining breast cancer risk prediction, beyond mammographic density alone.

A voxel-by-voxel method for mixing two filaments during a 3D printing process for soft-tissue replication in an anthropomorphic breast phantom

Objective. In this study, a novel voxel-by-voxel mixing method is presented, according to which two filaments of different material are combined during the three dimensional (3D) printing process.Approach. In our approach, two types of filaments were used for the replication of soft-tissues, a polylactic acid (PLA) filament and a polypropylene (PP) filament. A custom-made software was used, while a series of breast patient CT scan images were directly associated to the 3D printing process. Each phantom´s layer was printed twice, once with the PLA filament and a second time with the PP filament. For each material, the filament extrusion rate was controlled voxel-by-voxel and was based on the Hounsfield units (HU) of the imported CT images. The phantom was scanned at clinical CT, breast tomosynthesis and micro CT facilities, as the major processing was performed on data from the CT. A side by side comparison between patient´s and phantom´s CT slices by means of profile and histogram comparison was accomplished. Further, in case of profile comparison, the Pearson´s coefficients were calculated.Main results. The visual assessment of the distribution of the glandular tissue in the CT slices of the printed breast anatomy showed high degree of radiological similarity to the corresponding patient´s glandular distribution. The profile plots´ comparison showed that the HU of the replicated and original patient soft tissues match adequately. In overall, the Pearson´s coefficients were above 0.91, suggesting a close match of the CT images of the phantom with those of the patient. The overall HU were close in terms of HU ranges. The HU mean, median and standard deviation of the original and the phantom CT slices were -149, -167, ±65 and -121, -130, ±91, respectively.Significance. The results suggest that the proposed methodology is appropriate for manufacturing of anthropomorphic soft tissue phantoms for x-ray imaging and dosimetry purposes, since it may offer an accurate replication of these tissues.

Performance evaluation of digital breast tomosynthesis systems: physical methods and experimental data

Digital breast tomosynthesis (DBT) has become a well-established breast imaging technique, whose performance has been investigated in many clinical studies, including a number of prospective clinical trials. Results from these studies generally point to non-inferiority in terms of microcalcification detection and superior mass-lesion detection for DBT imaging compared to digital mammography (DM). This modality has become an essential tool in the clinic for assessment and ad-hoc screening but is not yet implemented in most breast screening programmes at a state or national level. While evidence on the clinical utility of DBT has been accumulating, there has also been progress in the development of methods for technical performance assessment and quality control of these imaging systems. DBT is a relatively complicated ‘pseudo-3D’ modality whose technical assessment poses a number of difficulties. This paper reviews methods for the technical performance assessment of DBT devices, starting at the component level in part one and leading up to discussion of system evaluation with physical test objects in part two. We provide some historical and basic theoretical perspective, often starting from methods developed for DM imaging. Data from a multi-vendor comparison are also included, acquired under the medical physics quality control protocol developed by EUREF and currently being consolidated by a European Federation of Organisations for Medical Physics working group. These data and associated methods can serve as a reference for the development of reference data and provide some context for clinical studies.

Breast dosimetry in alternative X-ray-based imaging modalities used in current clinical practices

In X-ray breast imaging, Digital Mammography (DM) and Digital Breast Tomosynthesis (DBT), are the standard and largely used techniques, both for diagnostic and screening purposes. Other techniques, such as dedicated Breast Computed Tomography (BCT) and Contrast Enhanced Mammography (CEM) have been developed as an alternative or a complementary technique to the established ones. The performance of these imaging techniques is being continuously assessed to improve the image quality and to reduce the radiation dose. These imaging modalities are predominantly used in the diagnostic setting to resolve incomplete or indeterminate findings detected with conventional screening examinations and could potentially be used either as an adjunct or as a primary screening tool in select populations, such as for women with dense breasts. The aim of this review is to describe the radiation dosimetry for these imaging techniques, and to compare the mean glandular dose with standard breast imaging modalities, such as DM and DBT.

Mammographic breast density and the risk of breast cancer: A systematic review and meta-analysis

OBJECTIVES

Mammographic density is a well-defined risk factor for breast cancer and having extremely dense breast tissue is associated with a one-to six-fold increased risk of breast cancer. However, it is questioned whether this increased risk estimate is applicable to current breast density classification methods. Therefore, the aim of this study was to further investigate and clarify the association between mammographic density and breast cancer risk based on current literature.

METHODS

Medline, Embase and Web of Science were systematically searched for articles published since 2013, that used BI-RADS lexicon 5th edition and incorporated data on digital mammography. Crude and maximally confounder-adjusted data were pooled in odds ratios (ORs) using random-effects models. Heterogeneity regarding breast cancer risks were investigated using I² statistic, stratified and sensitivity analyses.

RESULTS

Nine observational studies were included. Having extremely dense breast tissue (BI-RADS density D) resulted in a 2.11-fold (95% CI 1.84-2.42) increased breast cancer risk compared to having scattered dense breast tissue (BI-RADS density B). Sensitivity analysis showed that when only using data that had adjusted for age and BMI, the breast cancer risk was 1.83-fold (95% CI 1.52-2.21) increased. Both results were statistically significant and homogenous.

CONCLUSIONS

Mammographic breast density BI-RADS D is associated with an approximately two-fold increased risk of breast cancer compared to having BI-RADS density B in general population women. This is a novel and lower risk estimate compared to previously reported and might be explained due to the use of digital mammography and BI-RADS lexicon 5th edition.

Attenuation coefficient in the energy range 14–36 keV of 3D printing materials for physical breast phantoms

Objective. To measure the monoenergetic x-ray linear attenuation coefficient, μ, of fused deposition modeling (FDM) colored 3D printing materials (ABS, PLAwhite, PLAorange, PET and NYLON), used as adipose, glandular or skin tissue substitutes for manufacturing physical breast phantoms. Approach. Attenuation data (at 14, 18, 20, 24, 28, 30 and 36 keV) were acquired at Elettra synchrotron radiation facility, with step-wedge objects, using the Lambert–Beer law and a CCD imaging detector. Test objects were 3D printed using the Ultimaker 3 FDM printer. PMMA, Nylon-6 and high-density polyethylene step objects were also investigated for the validation of the proposed methodology. Printing uniformity was assessed via monoenergetic and polyenergetic imaging (32 kV, W/Rh). Main results. Maximum absolute deviation of μ for PMMA, Nylon-6 and HD-PE was 5.0%, with reference to literature data. For ABS and NYLON, μ differed by less than 6.1% and 7.1% from that of adipose tissue, respectively; for PET and PLAorange the difference was less than 11.3% and 6.3% from glandular tissue, respectively. PLAorange is a good substitute of skin (differences from −9.4% to +1.2%). Hence, ABS and NYLON filaments are suitable adipose tissue substitutes, while PET and PLAorange mimick the glandular tissue. PLAwhite could be printed at less than 100% infill density for matching the attenuation of glandular tissue, using the measured density calibration curve. The printing mesh was observed for sample thicknesses less than 60 mm, imaged in the direction normal to the printing layers. Printing dimensional repeatability and reproducibility was less 1%. Significance. For the first time an experimental determination was provided of the linear attenuation coefficient of common 3D printing filament materials with estimates of μ at all energies in the range 14–36 keV, for their use in mammography, breast tomosynthesis and breast computed tomography investigations.

The association between breast density and breast cancer pathological response to neoadjuvant chemotherapy

Purpose

Mammographic Density (MD) refers to the amount of fibroglandular breast tissue present in the breast and is an established risk factor for developing breast cancer. The ability to evaluate treatment response dynamically renders neoadjuvant chemotherapy (NACT) the preferred treatment option in many clinical scenarios. Previous studies have suggested that MD can predict patients likely to achieve a pathological complete response (pCR) to NACT. We aimed to determine whether there is a causal relationship between BI-RADS breast composition categories for breast density at diagnosis and the pCR rate and residual cancer burden score (RCB) by performing a retrospective review on consecutive breast cancer patients who received NACT in a tertiary referral centre from 2015 to 2021.

Methods

The Mann–Whitney U Test was used to test for differences between two independent groups (i.e. those who achieved pCR and those who did not). A binary logistic regression model was used to estimate odds ratios (OR) and corresponding 95% confidence intervals (CI) for an association between the independent variables of molecular subtype, MD, histological grade and FNA positivity and the dependant variable of pCR. Statistical analysis was conducted with SPSS (IBM SPSS for Mac, Version 26.0; IBM Corp).

Results

292 patients were included in the current study. There were 124, 155 and 13 patients in the BI-RADS MD category b, c and d, respectively. There were no patients in the BI-RADS MD category a. The patients with less dense breast composition (MD category b) were significantly older than patients with denser breast composition (MD category c, d) (p = 0.001) and patients who had a denser breast composition (MD category d) were more likely to have ER+ tumours. There was no significant difference in PgR status, HER2 status, pathological complete response (pCR), FNA positivity, or RCB class dependent upon the three MD categories. A binary logistic regression revealed that patients with HER2-enriched breast cancer and triple-negative breast cancer are more likely to achieve pCR with an OR of 3.630 (95% CI 1.360–9.691, p = 0.010) and 2.445 (95% CI 1.131–5.288, p = 0.023), respectively.

Conclusion

Whilst dense MD was associated with ER positivity and these women were less likely to achieve a pCR, MD did not appear to independently predict pCR post-NACT.

Review of Thermal and Physiological Properties of Human Breast Tissue

Electromagnetic thermal therapies for cancer treatment, such as microwave hyperthermia, aim to heat up a targeted tumour site to temperatures within 40 and 44 °C. Computational simulations used to investigate such heating systems employ the Pennes’ bioheat equation to model the heat exchange within the tissue, which accounts for several tissue properties: density, specific heat capacity, thermal conductivity, metabolic heat generation rate, and blood perfusion rate. We present a review of these thermal and physiological properties relevant for hyperthermia treatments of breast including fibroglandular breast, fatty breast, and breast tumours. The data included in this review were obtained from both experimental measurement studies and estimated properties of human breast tissues. The latter were used in computational studies of breast thermal treatments. The measurement methods, where available, are discussed together with the estimations and approximations considered for values where measurements were unavailable. The review concludes that measurement data for the thermal and physiological properties of breast and tumour tissue are limited. Fibroglandular and fatty breast tissue properties are often approximated from those of generic muscle or fat tissue. Tumour tissue properties are mostly obtained from approximating equations or assumed to be the same as those of glandular tissue. We also present a set of reliable data, which can be used for more accurate modelling and simulation studies to better treat breast cancer using thermal therapies.

Effect of optical blurring of X-ray source on breast tomosynthesis image quality: Modulation transfer function, anatomical noise power spectrum, and signal detectability perspectives

We investigated the effect of the optical blurring of X-ray source on digital breast tomosynthesis (DBT) image quality using well-designed DBT simulator and table-top experimental systems. To measure the in-plane modulation transfer function (MTF), we used simulated sphere phantom and Teflon sphere phantom and generated their projection data using two acquisition modes (i.e., step-and-shoot mode and continuous mode). After reconstruction, we measured the in-plane MTF using reconstructed sphere phantom images. In addition, we measured the anatomical noise power spectrum (aNPS) and signal detectability. We constructed simulated breast phantoms with a 50% volume glandular fraction (VGF) of breast anatomy using the power law spectrum and inserted spherical objects with 1 mm, 2 mm, and 5 mm diameters as breast masses. Projection data were acquired using two acquisition modes, and in-plane breast images were reconstructed using the Feldkamp-Davis-Kress (FDK) algorithm. For the experimental study, we used BR3D breast phantom with 50% VGF and obtained projection data using a table-top experimental system. To compare the detection performance of the two acquisition modes, we calculated the signal detectability using the channelized Hotelling observer (CHO) with Laguerre-Gauss (LG) channels. Our results show that spatial resolution of in-plane image in continuous mode was degraded due to the optical blurring of X-ray source. This blurring effect was reflected in aNPS, resulting in large β values. From a signal detectability perspective, the signal detectability in step-and-shoot mode is higher than that in continuous mode for small spherical signals but not large spherical signals. Although the step-and-shoot mode has disadvantage in terms of scan time compared to the continuous mode, scanning in step-and-shoot mode is better for detecting small signals, indicating that there is a tradeoff between scan time and image quality.

Two-phase learning-based 3D deblurring method for digital breast tomosynthesis images

In digital breast tomosynthesis (DBT) systems, projection data are acquired from a limited number of angles. Consequently, the reconstructed images contain severe blurring artifacts that might heavily degrade the DBT image quality and cause difficulties in detecting lesions. In this study, we propose a two-phase learning approach for artifact compensation in a coarse-to-fine manner to mitigate blurring artifacts effectively along all viewing directions of the DBT image volume (i.e., along the axial, coronal, and sagittal planes) to improve the detection performance of lesions. The proposed method employs a convolutional neural network model comprising two submodels/phases, with Phase 1 performing three-dimensional (3D) deblurring and Phase 2 performing additional 2D deblurring. To investigate the effects of loss functions on the proposed model’s deblurring performance, we evaluated several loss functions, such as the pixel-based loss function, adversarial-based loss function, and perception-based loss function. Compared with the DBT image, the mean squared error of the image and the root mean squared errors of the gradient of the image decreased by 82.8% and 44.9%, respectively, and the contrast-to-noise ratio increased by 183.4% in the in-focus plane. We verified that the proposed method sequentially restored the missing frequency components as the DBT images were processed through the Phase 1 and Phase 2 steps. These results indicate that the proposed method performs effective 3D deblurring, significantly reducing the blurring artifacts in the in-focus plane and other planes of the DBT image, thus improving the detection performance of lesions.

Biological Mechanisms and Therapeutic Opportunities in Mammographic Density and Breast Cancer Risk

Mammographic density is an important risk factor for breast cancer; women with extremely dense breasts have a four to six fold increased risk of breast cancer compared to women with mostly fatty breasts, when matched with age and body mass index. High mammographic density is characterised by high proportions of stroma, containing fibroblasts, collagen and immune cells that suggest a pro-tumour inflammatory microenvironment. However, the biological mechanisms that drive increased mammographic density and the associated increased risk of breast cancer are not yet understood. Inflammatory factors such as monocyte chemotactic protein 1, peroxidase enzymes, transforming growth factor beta, and tumour necrosis factor alpha have been implicated in breast development as well as breast cancer risk, and also influence functions of stromal fibroblasts. Here, the current knowledge and understanding of the underlying biological mechanisms that lead to high mammographic density and the associated increased risk of breast cancer are reviewed, with particular consideration to potential immune factors that may contribute to this process.

DIR-DBTnet: Deep iterative reconstruction network for three-dimensional digital breast tomosynthesis imaging

PURPOSE

The goal of this study is to develop a three-dimensional (3D) iterative reconstruction framework based on the deep learning (DL) technique to improve the digital breast tomosynthesis (DBT) imaging performance.

METHODS

In this work, the DIR-DBTnet is developed for DBT image reconstruction by mapping the conventional iterative reconstruction (IR) algorithm to the deep neural network. By design, the DIR-DBTnet learns and optimizes the regularizer and the iteration parameters automatically during the network training with a large amount of simulated DBT data. Numerical, experimental, and clinical data are used to evaluate its performance. Quantitative metrics such as the artifact spread function (ASF), breast density, and the signal difference to noise ratio (SDNR) are measured to assess the image quality.

RESULTS

Results show that the proposed DIR-DBTnet is able to reduce the in-plane shadow artifacts and the out-of-plane signal leaking artifacts compared to the filtered backprojection (FBP) and the total variation (TV)-based IR methods. Quantitatively, the full width half maximum (FWHM) of the measured ASF from the clinical data is 27.1% and 23.0% smaller than those obtained with the FBP and TV methods, while the SDNR is increased by 194.5% and 21.8%, respectively. In addition, the breast density obtained from the DIR-DBTnet network is more accurate and consistent with the ground truth.

CONCLUSIONS

In conclusion, a deep iterative reconstruction network, DIR-DBTnet, has been proposed for 3D DBT image reconstruction. Both qualitative and quantitative analyses of the numerical, experimental, and clinical results demonstrate that the DIR-DBTnet has superior DBT imaging performance than the conventional algorithms.

In Silico Phase-Contrast X-Ray Imaging of Anthropomorphic Voxel-Based Phantoms

Propagation-based phase-contrast X-ray imaging is an emerging technique that can improve dose efficiency in clinical imaging. In silico tools are key to understanding the fundamental imaging mechanisms and develop new applications. Here, due to the coherent nature of the phase-contrast effects, tools based on wave propagation (WP) are preferred over Monte Carlo (MC) based methods. WP simulations require very high wave-front sampling which typically limits simulations to small idealized objects. Virtual anthropomorphic voxel-based phantoms are typically provided with a resolution lower than imposed sampling requirements and, thus, cannot be directly translated for use in WP simulations. In the present paper we propose a general strategy to enable the use of these phantoms for WP simulations. The strategy is based on upsampling in the 3D domain followed by projection resulting in high-resolution maps of the projected thickness for each phantom material. These maps can then be efficiently used for simulations of Fresnel diffraction to generate in silico phase-contrast X-ray images. We demonstrate the strategy on an anthropomorphic breast phantom to simulate propagation-based phase-contrast mammography using a laboratory micro-focus X-ray source.

Comprehensive characterization of ExacTrac stereoscopic image guidance system using Monte Carlo and Spektr simulations

The purpose of this work is to develop accurate computational methods to comprehensively characterize and model the clinical ExacTrac imaging system, which is used as an image guidance system for stereotactic treatment applications. The Spektr toolkit was utilized to simulate the spectral and imaging characterization of the system. Since Spektr only simulates the primary beam (ignoring scatter), a full model of ExacTrac was also developed in Monte Carlo (MC) to characterize the imaging system. To ensure proper performance of both simulation models, Spektr and MC data were compared to the measured spectral and half value layers (HVLs) values. To validate the simulation results, x-ray spectra of the ExacTrac system were measured for various tube potentials using a CdTe spectrometer with multiple added narrow collimators. The raw spectra were calibrated using a ⁵⁷Co source and corrected for the escape peaks and detector efficiency. HVLs in mm of Al for various energies were measured using a calibrated RaySafe detector. Spektr and MC HVLs were calculated and compared to the measured values. The patient surface dose was calculated for different clinical imaging protocols from the measured air kerma and HVL values following the TG-61 methodology. The x-ray focal spot was measured by slanted edge technique using gafchromic films. ExacTrac imaging system beam profiles were simulated for various energies by MC simulation and the results were benchmarked by experimentally acquired beam profiles using gafchromic films. The effect of 6D IGRT treatment couch on beam hardening, dynamic range of the flat panel detector and scatter effect were determined using both Spektr simulation and experimental measurements. The measured and simulated spectra (of both MC and Spektr) for various kVps were compared and agreed within acceptable error. As another validation, the measured HVLs agreed with the Spektr and MC simulated HVLs on average within 1.0% for all kVps. The maximum and minimum patient surface doses were found to be 1.06 mGy for shoulder (high) and 0.051 mGy for cranial (low) imaging protocols, respectively. The MC simulated beam profiles were well matched with experimental results and replicated the penumbral slopes, the heel effect, and out-of-field regions. Dynamic range of detector (in terms of air kerma at detector surface) was found to be in the range of [6.1 × 10^-6, 5.3 × 10^-3] mGy. Accurate MC and Spektr models of the ExacTrac image guidance system were successfully developed and benchmarked via experimental validation. While patient surface dose for available imaging protocols were reported in this study, the established MC model may be used to obtain 3D imaging dose distribution for real patient geometries.

Microwave Breast Imaging Using Rotational Bistatic Impulse Radar for the Detection of Breast Cancer: Protocol for a Prospective Diagnostic Study

Background

Mammography is the standard examination for breast cancer screening; however, it is associated with pain and exposure to ionizing radiation. Microwave breast imaging is a less invasive method for breast cancer surveillance. A bistatic impulse radar–based breast cancer detector has recently been developed.

Objective

This study aims to present a protocol for evaluating the diagnostic accuracy of the novel microwave breast imaging device.

Methods

This is a prospective diagnostic study. A total of 120 participants were recruited before treatment administration and divided into 2 cohorts: 100 patients diagnosed with breast cancer and 20 participants with benign breast tumors. The detector will be directly placed on each breast, while the participant is in supine position, without a coupling medium. Confocal images will be created based on the analyzed data, and the presence of breast tumors will be assessed. The primary endpoint will be the diagnostic accuracy, sensitivity, and specificity of the detector for breast cancer and benign tumors. The secondary endpoint will be the safety and detectability of each molecular subtype of breast cancer. For an exploratory endpoint, the influence of breast density and tumor size on tumor detection will be investigated.

Results

Recruitment began in November 2018 and was completed by March 2020. We anticipate the preliminary results to be available by summer 2021.

Conclusions

This study will provide insights on the diagnostic accuracy of microwave breast imaging using a rotational bistatic impulse radar. The collected data will improve the diagnostic algorithm of microwave imaging and lead to enhanced device performance.

Trial Registration

Japan Registry of Clinical Trials jRCTs062180005; https://jrct.niph.go.jp/en-latest-detail/jRCTs062180005

International Registered Report Identifier (IRRID)

DERR1-10.2196/17524

On the correlation between second order texture features and human observer detection performance in digital images

Image texture, the relative spatial arrangement of intensity values in an image, encodes valuable information about the scene. As it stands, much of this potential information remains untapped. Understanding how to decipher textural details would afford another method of extracting knowledge of the physical world from images. In this work, we attempt to bridge the gap in research between quantitative texture analysis and the visual perception of textures. The impact of changes in image texture on human observer's ability to perform signal detection and localization tasks in complex digital images is not understood. We examine this critical question by studying task-based human observer performance in detecting and localizing signals in tomographic breast images. We have also investigated how these changes impact the formation of second-order image texture. We used digital breast tomosynthesis (DBT) an FDA approved tomographic X-ray breast imaging method as the modality of choice to show our preliminary results. Our human observer studies involve localization ROC (LROC) studies for low contrast mass detection in DBT. Simulated images are used as they offer the benefit of known ground truth. Our results prove that changes in system geometry or processing leads to changes in image texture magnitudes. We show that the variations in several well-known texture features estimated in digital images correlate with human observer detection-localization performance for signals embedded in them. This insight can allow efficient and practical techniques to identify the best imaging system design and algorithms or filtering tools by examining the changes in these texture features. This concept linking texture feature estimates and task based image quality assessment can be extended to several other imaging modalities and applications as well. It can also offer feedback in system and algorithm designs with a goal to improve perceptual benefits. Broader impact can be in wide array of areas including imaging system design, image processing, data science, machine learning, computer vision, perceptual and vision science. Our results also point to the caution that must be exercised in using these texture features as image-based radiomic features or as predictive markers for risk assessment as they are sensitive to system or image processing changes.

Age-Specific Breast Density Changes in Taiwanese Women: A Cross-Sectional Study

Breast density is a risk factor for breast cancer. This study explored distribution of mammographic density quantitatively and qualitatively in a wide age range of Taiwanese women. Subjects with negative and benign mammographic findings were included. According to the Breast Imaging Reporting and Data System, the proportion of extremely dense breasts declined from 58.0% in women < 30 years to 1.9% in women > 74 years. More than 80% of mammograms in women < 55 years old were classified as extremely or heterogeneously dense, while the proportion of dense breasts was still high in women aged 60–64 years (59.3%). The absolute dense area of the breast declined from 35.8% in women < 30 years to 18.5% in women > 74 years. The correlation between breast density and age was significant, with and without controlling for the effect of body composition (p < 0.001), implying that the relationship between breast density and age was not wholly related to body composition. In conclusion, the higher breast density in Taiwanese women aged 60–64 years was comparable to that of Western women aged 40–44 years in the literature. This suggests that breast cancer screening using mammography may be more challenging for Asian women than for Western women of the same age.

Effects of kV, filtration, dose, and object size on soft tissue and iodine contrast in dedicated breast CT

PURPOSE

Clinical 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.

METHODS

This 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.

RESULTS

Air 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' (R²  > 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.

CONCLUSIONS

The 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.

Multimodal Breast Phantoms for Microwave, Ultrasound, Mammography, Magnetic Resonance and Computed Tomography Imaging

The aim of this work was to develop multimodal anthropomorphic breast phantoms suitable for evaluating the imaging performance of a recently-introduced Microwave Imaging (MWI) technique in comparison to the established diagnostic imaging modalities of Magnetic Resonance Imaging (MRI), Ultrasound (US), mammography and Computed Tomography (CT). MWI is an emerging technique with significant potential to supplement established imaging techniques to improve diagnostic confidence for breast cancer detection. To date, numerical simulations have been used to assess the different MWI scanning and image reconstruction algorithms in current use, while only a few clinical trials have been conducted. To bridge the gap between the numerical simulation environment and a more realistic diagnostic scenario, anthropomorphic phantoms which mimic breast tissues in terms of their heterogeneity, anatomy, morphology, and mechanical and dielectric characteristics, may be used. Key in this regard is achieving realism in the imaging appearance of the different healthy and pathologic tissue types for each of the modalities, taking into consideration the differing imaging and contrast mechanisms for each modality. Suitable phantoms can thus be used by radiologists to correlate image findings between the emerging MWI technique and the more familiar images generated by the conventional modalities. Two phantoms were developed in this study, representing difficult-to-image and easy-to-image patients: the former contained a complex boundary between the mammary fat and fibroglandular tissues, extracted from real patient MRI datasets, while the latter contained a simpler and less morphologically accurate interface. Both phantoms were otherwise identical, with tissue-mimicking materials (TMMs) developed to mimic skin, subcutaneous fat, fibroglandular tissue, tumor and pectoral muscle. The phantoms’ construction used non-toxic materials, and they were inexpensive and relatively easy to manufacture. Both phantoms were scanned using conventional modalities (MRI, US, mammography and CT) and a recently introduced MWI radar detection procedure called in-coherent Multiple Signal Classification (I-MUSIC). Clinically realistic artifact-free images of the anthropomorphic breast phantoms were obtained using the conventional imaging techniques as well as the emerging technique of MWI.

Human observer performance on in-plane digital breast tomosynthesis images: Effects of reconstruction filters and data acquisition angles on signal detection

For digital breast tomosynthesis (DBT) systems, we investigate the effects of the reconstruction filters for different data acquisition angles on signal detection. We simulated a breast phantom with a 30% volume glandular fraction (VGF) of breast anatomy using the power law spectrum and modeled the breast mass as a spherical object with a 1 mm diameter. Projection data were acquired using two different data acquisition angles and numbers of projection view pairs, and in-plane breast images were reconstructed using the Feldkamp-Davis-Kress (FDK) algorithm with three different reconstruction filter schemes. To measure the ability to detect a signal, we conducted the human observer study with a binary detection task and compared the signal detectability of human to that of channelized Hotelling observer (CHO) with Laguerre-Gauss (LG) channels and dense difference-of-Gaussian (D-DOG) channels. We also measured the contrast-to-noise ratio (CNR), signal power spectrum (SPS), and β values of the anatomical noise power spectrum (NPS) to show the association between human observer performance and these traditional metrics. Our results show that using a slice thickness (ST) filter degraded the signal detection performance of human observers at the same data acquisition angle. This could be predicted by D-DOG CHO with internal noise, but the correlation between the traditional metrics and signal detectability was not observed in this work.

Three-layer heterogeneous mammographic phantoms for Monte Carlo simulation of normalized glandular dose coefficients in mammography

Normalized glandular dose (DgN) coefficients obtained using homogeneous breast phantoms are commonly used in breast dosimetry for mammography. However, glandular tissue is heterogeneously distributed in the breast. This study aimed to construct three-layer heterogeneous mammographic phantoms (THEPs) to examine the effect of glandular distribution on DgN coefficient. Each layer of THEPs was set to 25%, 50%, or 75% glandular fraction to emulate heterogeneous glandular distribution. Monte Carlo simulation was performed to attain mean glandular dose (MGD) and air kerma at 22-36 kVp and W/Al, W/Rh, and W/Ag target-filter combinations. The heterogeneous DgN coefficient was calculated as functions of the mean glandular fraction (MGF), breast thickness, tube voltage, and half-value layer. At 50% MGF, the heterogeneous DgN coefficients for W/Al, W/Rh, and W/Ag differed by 40.3%, 36.7%, and 31.2%. At 9-cm breast thickness, the DgN values of superior and inferior glandular distributions were 25.4% higher and 29.2% lower than those of uniform distribution. The proposed THEPs can be integrated with conventional breast dosimetry to consider the heterogeneous glandular distribution in clinical practice.

A fluence modulation and scatter shielding apparatus for dedicated breast CT: Theory of operation

PURPOSE

To introduce an auxiliary apparatus of fluence modulation and scatter shielding for dedicated breast computed tomography (bCT) and the corresponding patient-specific method of image acquisition.

METHODS

The apparatus is composed of two assemblies, referred to herein as the "Dynamic Fluence Gate" (FG) and "Scatter Shield" (SS), that work in synchrony to form a narrow beam sweeping the entire fan angle coverage of the imaging system during a projection. The apparatus, as a whole, is referred to as FG-SS. FG and SS are pre-patient and post-patient units, respectively. Each is composed of a rotating drum, on top of which are installed two sheets of high x-ray attenuating material, a sensory system, and the constituent robotics. The sheets of each unit are positioned such that an opening - a window Fluence Modulation and Scatter Shielding is formed between them. The rotations of the drums and positioning of the sheets are synchronized and adjusted such that a line of sight is created between the source, FG window, the breast, and the SS window. With line of sight achieved, the narrow beam transitions from the source to the detector. The fluence of the narrow beam during a projection depends on the size, shape, and positioning of the breast. The FG-SS method of imaging is discussed mathematically and demonstrated using computer simulations. A series of Monte Carlo simulations are conducted to evaluate the performance of the system as relates to its impact on the imager's dynamic range, dose distribution to the breast, noise inhomogeneity in reconstructed images, and scatter buildup in projections within small, medium, and large breasts composed of homogeneous medium breast tissue.

RESULTS

Implementation of FG-SS results in near scatter-free projections, reduction in both dose and the required dynamic range of the imager, and equalization of the quantum noise distribution in the reconstructed image. Using the disclosed design, the dynamic range was reduced by factors ranging from 1.6 to 5.5 across the range of breast sizes studied. A reduction in the acquisition of the scattered rays, varying between the factors of 6.1 (in the small breast) and 9.8 (in that large breast) was achieved and consequently, shading artifacts were suppressed. Noise in the CT image was equalized by reducing the overall spatial variation from 29% to <5% in small breast and from 45% to 14% in the large breast. An overall reduction in deposited dose to the breast was achieved - between 26% and 39% depending on the breast size.

CONCLUSIONS

Utilization of the FG-SS apparatus and technique was demonstrated via simulations to result in: significant reductions in dose to the breast, reductions in scatter uptake in projections, reduced required dynamic range of the imager, and homogenizing of quantum noise throughout the reconstructed image.

Field cancerization in the understanding of parenchymal analysis of mammograms for breast cancer risk assessment

In recent years, mammographic image analysis has shown great potential for breast cancer risk assessment. The aim of risk assessment is to predict how likely a woman is to develop breast cancer in the future. Several studies suggest that computerized parenchymal analysis of mammograms can be utilized as an independent imaging biomarker of breast cancer. Parenchymal analysis consists of the quantitative assessment of visual texture patterns in mammograms to infer the level of risk. In spite of substantial evidence of the association between parenchymal patterns and breast cancer risk, its biological foundations remain poorly understood. In this work, we draw a hypothesis that links the field cancerization (FC) with breast cancer risk assessment based on the parenchymal analysis. In the literature, the FC is interpreted as a biochemical anomaly amplification in otherwise healthy cells due to the effect of pre-cancerous transformed cells in surrounding regions. Our hypothesis is that these biochemical anomaly amplifications change the cellular micro-environment which, in turn, alter tissue responses to X-ray radiation. As a result, it is reasonable to think that these changes influence the interaction of X-rays with parenchymal - the functional - breast tissue thus enabling cancer prediction by analyzing X-ray images of the breast. We believe that our hypothesis provides an actionable explanation as to how computerized parenchymal analysis of apparently normal mammograms can be successfully utilized for the stratification of breast cancer risk.

Cascaded systems analysis of anatomic noise in digital mammography and dual-energy digital mammography

In x-ray based imaging of the breast, contrast between fibroglandular (Fg) tissue and adipose (Ad) tissue is a source of anatomic noise. The goal of this work was to validate by simulation and experiment a mathematical framework for modelling the Fg component of anatomic noise in digital mammograpy (DM) and dual-energy (DE) DM. Our mathematical framework unifies and generalizes existing approaches. We compared mathematical predictions directly with empirical measurements of the anatomic noise power spectrum of the CIRS BR3D structured breast phantom using two clinical mammography systems and four beam qualities. Our simulation and experimental results showed agreement with mathematical predictions. As a demonstration of utility, we used our mathematical framework in a theoretical spectral optimization of DM for the task of detecting breast masses. Our theoretical optimization showed that the optimal tube voltage for DM may be higher than that based on predictions that do not account for anatomic noise, in agreement with recent theoretical findings. Additionally, our theoretical optimization predicts that filtering tungsten-anode x-ray spectra with rhodium has little influence on lesion detectability, in contrast with previous findings. The mathematical methods validated in this work can be incorporated easily into cascaded systems analysis of breast imaging systems and will be useful when optimizating novel techniques for x-ray-based imaging of the breast.

Mammographic density: intersection of advocacy, science, and clinical practice

Purpose

Here we aim to review the association between mammographic density, collagen structure and breast cancer risk.

Findings

While mammographic density is a strong predictor of breast cancer risk in populations, studies by Boyd show that mammographic density does not predict breast cancer risk in individuals. Mammographic density is affected by age, parity, menopausal status, race/ethnicity, and body mass index (BMI).New studies normalize mammographic density to BMI may provide a more accurate way to compare mammographic density in women of diverse race and ethnicity. Preclinical and tissue-based studies have investigated the role collagen composition and structure in predicting breast cancer risk. There is emerging evidence that collagen structure may activate signaling pathways associated with aggressive breast cancer biology.

Summary

Measurement of film mammographic density does not adequately capture the complex signaling that occurs in women with at-risk collagen. New ways to measure at-risk collagen potentially can provide a more accurate view of risk.

Quantitative characterization of breast tissues with dedicated CT imaging

A quantitative characterization of the soft tissues composing the human breast is achieved by means of a monochromatic CT phase-contrast imaging system, through accurate measurements of their attenuation coefficients within the energy range of interest for breast CT clinical examinations. Quantitative measurements of linear attenuation coefficients are performed on tomographic reconstructions of surgical samples, using monochromatic x-ray beams from a synchrotron source and a free space propagation setup. An online calibration is performed on the obtained reconstructions, in order to reassess the validity of the standard calibration procedure of the CT scanner. Three types of healthy tissues (adipose, glandular, and skin) and malignant tumors, when present, are considered from each sample. The measured attenuation coefficients are in very good agreement with the outcomes of similar studies available in the literature, although they span an energy range that was mostly neglected in the previous studies. No globally significant differences are observed between healthy and malignant dense tissues, although the number of considered samples does not appear sufficient to address the issue of a quantitative differentiation of tumors. The study assesses the viability of the proposed methodology for the measurement of linear attenuation coefficients, and provides a denser sampling of attenuation data in the energy range useful to breast CT.

Development of breast lesions models database

PURPOSE

We present the development and the current state of the MaXIMA Breast Lesions Models Database, which is intended to provide researchers with both segmented and mathematical computer-based breast lesion models with realistic shape.

METHODS

The database contains various 3D images of breast lesions of irregular shapes, collected from routine patient examinations or dedicated scientific experiments. It also contains images of simulated tumour models. In order to extract the 3D shapes of the breast cancers from patient images, an in-house segmentation algorithm was developed for the analysis of 50 tomosynthesis sets from patients diagnosed with malignant and benign lesions. In addition, computed tomography (CT) scans of three breast mastectomy cases were added, as well as five whole-body CT scans. The segmentation algorithm includes a series of image processing operations and region-growing techniques with minimal interaction from the user, with the purpose of finding and segmenting the areas of the lesion. Mathematically modelled computational breast lesions, also stored in the database, are based on the 3D random walk approach.

RESULTS

The MaXIMA Imaging Database currently contains 50 breast cancer models obtained by segmentation of 3D patient breast tomosynthesis images; 8 models obtained by segmentation of whole body and breast cadavers CT images; and 80 models based on a mathematical algorithm. Each record in the database is supported with relevant information. Two applications of the database are highlighted: inserting the lesions into computationally generated breast phantoms and an approach in generating mammography images with variously shaped breast lesion models from the database for evaluation purposes. Both cases demonstrate the implementation of multiple scenarios and of an unlimited number of cases, which can be used for further software modelling and investigation of breast imaging techniques. The created database interface is web-based, user friendly and is intended to be made freely accessible through internet after the completion of the MaXIMA project.

CONCLUSIONS

The developed database will serve as an imaging data source for researchers, working on breast diagnostic imaging and on improving early breast cancer detection techniques, using existing or newly developed imaging modalities.

Multi-Energy Computed Tomography Breast Imaging with Monte Carlo Simulations: Contrast-to-Noise-Based Image Weighting

Context:

Photon-counting detectors and breast computed tomography imaging have been an active area of research. With these detectors, photons are assigned an equal weight and weighting schemes can be enabled. More weight can be assigned to lower energies, resulting in an increase in the contrast-to-noise ratio (CNR).

Aims:

The aim of this study is to develop and evaluate an energy weighting imaging technique to improve the CNR of simulated breast phantoms and to improve tumour detection.

Materials and Methods:

Breast phantoms consisting of adipose, glandular, malignant tissues and iodine contrast were constructed with BreastSimulator software. The phantoms were used in egs_cbct simulations for energies ranging between 20 and 65 keV from which multiple images were reconstructed. A new CNR-based image weighting method was proposed based on the CNR values obtained from the images. This method improves on previous methods and can be applied to complicated phantoms since no structural information is needed.

Results:

An increase in the CNR can be seen for lower energies. A sharp increase in the CNR is seen just above the K-edge for the phantoms with the iodine contrast. The CNR-based image weighting leads to a 68.47% (1.68-fold) increase in the CNR for the malignant tissue without iodine. For the malignant tissue with iodine contrast, the increase in the CNR was 96.14% (1.96-fold).

Conclusions:

The new proposed CNR-based image weighting scheme is easy to implement and can be used for complicated phantoms with varying structures. A large increase in the CNR is seen with or without the use of iodine contrast.

Three-dimensionally-printed anthropomorphic physical phantom for mammography and digital breast tomosynthesis with custom materials, lesions, and uniform quality control region

Anthropomorphic breast phantoms mimic patient anatomy in order to evaluate clinical mammography and digital breast tomosynthesis system performance. Our goal is to create a modular phantom with an anthropomorphic region to allow for improved lesion and calcification detection as well as a uniform region to evaluate standard quality control (QC) metrics. Previous versions of this phantom used commercial photopolymer inkjet three-dimensional printers to recreate breast anatomy using four surfaces that were fabricated with commercial materials spanning only a limited breast density range of 36% to 64%. We use modified printers to create voxelized, dithered breast phantoms with continuous gradations between glandular and adipose tissues. Moreover, the new phantom replicates the low-end density (representing adipose tissue) using third party material, Jf Flexible, and increases the high-end density to the density of glandular tissue and beyond by either doping Jf Flexible with salts and nanoparticles or using a new commercial resin, VeroPureWhite. An insert design is utilized to add masses, calcifications, and iodinated objects into the phantom for increased utility. The uniform chest wall region provides a space for traditional QC objects such as line pair patterns for measuring resolution and scale bars for measuring printer linearity. Incorporating these distinct design modules enables us to create an improved, complete breast phantom to better evaluate clinical mammography systems for lesion and calcification detection and standard QC performance evaluation.

Preliminary investigation of mammographic density among women in Riyadh: association with breast cancer risk factors and implications for screening practices

PURPOSE

Improved mammographic sensitivity is associated with reduced mammographic density. This study aims to: provide a preliminary report on mammographic density among women in Riyadh; identify risk factors associated with mammographic density; and consider the potential implications for screening practices.

METHODS

Based on a cross-sectional design, we examined a total of 792 women using an area-based mammographic density method (LIBRA). Spearman's correlation, Mann-Whitney U, Kruskal-Wallis and binary logistic regression were used for analyses.

RESULTS

The study population had a mean age of 49.6 years and a high proportion of participants were overweight or obese (90%). A large number of women had low mammographic density, with a mean dense breast area of 19.1 cm² and percent density of 10.3 cm². Slightly more than half of the variations in the dense breast area and percent density models were explained by BMI. In the adjusted analyses, BMI, menopausal status, age at menarche and number of children remained statistically significant predictors.

CONCLUSION

Given the high proportion of women with low mammographic density, our findings suggest that women living in Riyadh may not require additional imaging strategies beyond mammography to detect breast cancers. The high proportion of obese women reported here and across Saudi Arabia suggests that mammographic density is less likely to have an adverse impact on mammographic sensitivity. Thus and to improve clinical outcomes among Saudi women, annual mammography and commencing screening at a younger age are suggested. Additional studies are required to shed further light on our findings.

Non-stationary model of oblique x-ray incidence in amorphous selenium detectors: I. Point spread function

PURPOSE

In previous work, a theoretical model of the point spread function (PSF) for oblique x-ray incidence in amorphous selenium (a-Se) detectors was proposed. The purpose of this paper is to develop a complementary model that includes two additional features. First, the incidence angle and the directionality of ray incidence are calculated at each position, assuming a divergent x-ray beam geometry. This approach allows the non-stationarity of the PSF to be modeled. Second, this paper develops a framework that is applicable to a digital system, unlike previous work which did not model the presence of a thin-film transistor (TFT) array.

METHODS

At each point on the detector, the incidence angle and the ray incidence direction are determined using ray tracing. Based on these calculations, an existing model for the PSF of the x-ray converter (Med Phys. 1995;22:365-374) is generalized to a non-stationary model. The PSF is convolved with the product of two rectangle functions, which model the sampling of the TFT array. The rectangle functions match the detector element (del) size in two dimensions.

RESULTS

It is shown that the PSF can be calculated in closed form. This solution is used to simulate a digital mammography (DM) system at two x-ray energies (20 and 40 keV). Based on the divergence of the x-ray beam, the direction of ray incidence varies with position. Along this direction, the PSF is broader than the reference rect function matching the del size. The broadening is more pronounced with increasing obliquity. At high energy, the PSF deviates more strongly from the reference rect function, indicating that there is more blurring. In addition, the PSF is calculated along the polar angle perpendicular to the ray incidence direction. For this polar angle, the shape of the PSF is dependent upon whether the ray incidence direction is parallel with the sides of the detector. If the ray incidence direction is parallel with either dimension, the PSF is a perfect rectangle function, matching the del size. However, if the ray incidence direction is at an oblique angle relative to the sides of the detector, the PSF is not rectangular. These results illustrate the non-stationarity of the PSF.

CONCLUSIONS

This paper demonstrates that an existing model of the PSF of a-Se detectors can be generalized to include the effects of non-stationarity and digitization. The PSF is determined in closed form. This solution offers the advantage of shorter computation time relative to approaches that use numerical methods. This model is a tool for simulating a-Se detectors in future work, such as in virtual clinical trials with computational phantoms.

Monochromatic breast computed tomography with synchrotron radiation: phase-contrast and phase-retrieved image comparison and full-volume reconstruction

A program devoted to performing the first in vivo synchrotron radiation (SR) breast computed tomography (BCT) is ongoing at the Elettra facility. Using the high spatial coherence of SR, phase-contrast (PhC) imaging techniques can be used. The latest high-resolution BCT acquisitions of breast specimens, obtained with the propagation-based PhC approach, are herein presented as part of the SYRMA-3D collaboration effort toward the clinical exam. Images are acquired with a 60 - μ m   pixel dead-time-free single-photon-counting CdTe detector. The samples are imaged at 32 and 38 keV in a continuous rotating mode, delivering 5 to 20 mGy of mean glandular dose. Contrast-to-noise ratio (CNR) and spatial resolution performances are evaluated for both PhC and phase-retrieved images, showing that by applying the phase-retrieval algorithm a five-time CNR increase can be obtained with a minor loss in spatial resolution across soft tissue interfaces. It is shown that, despite having a poorer CNR, PhC images can provide a sharper visualization of microcalcifications, thus being complementary to phase-retrieved images. Furthermore, the first full-volume scan of a mastectomy sample ( 9 × 9 × 3    cm 3 ) is reported. This investigation into surgical specimens indicates that SR BCT in terms of CNR, spatial resolution, scan duration, and scan volume is feasible.

Measurement of breast-tissue x-ray attenuation by spectral imaging: fresh and fixed normal and malignant tissue

Knowledge of x-ray attenuation is essential for developing and evaluating x-ray imaging technologies. In mammography, measurement of breast density, dose estimation, and differentiation between cysts and solid tumours are example applications requiring accurate data on tissue attenuation. Published attenuation data are, however, sparse and cover a relatively wide range. To supplement available data we have previously measured the attenuation of cyst fluid and solid lesions using photon-counting spectral mammography. The present study aims to measure the attenuation of normal adipose and glandular tissue, and to measure the effect of formalin fixation, a major uncertainty in published data. A total of 27 tumour specimens, seven fibro-glandular tissue specimens, and 15 adipose tissue specimens were included. Spectral (energy-resolved) images of the samples were acquired and the image signal was mapped to equivalent thicknesses of two known reference materials, from which x-ray attenuation as a function of energy can be derived. The spread in attenuation between samples was relatively large, partly because of natural variation. The variation of malignant and glandular tissue was similar, whereas that of adipose tissue was lower. Formalin fixation slightly altered the attenuation of malignant and glandular tissue, whereas the attenuation of adipose tissue was not significantly affected. The difference in attenuation between fresh tumour tissue and cyst fluid was smaller than has previously been measured for fixed tissue, but the difference was still significant and discrimination of these two tissue types is still possible. The difference between glandular and malignant tissue was close-to significant; it is reasonable to expect a significant difference with a larger set of samples. We believe that our studies have contributed to lower the overall uncertainty of breast tissue attenuation in the literature due to the relatively large sample sets, the novel measurement method, and by clarifying the difference between fresh and fixed tissue.

Optimizing illumination for full field imaging at high brilliance hard X-ray synchrotron sources

A new technique is presented to overcome beam size limitation in full field imaging at high brilliance synchrotron sources using specially designed refractive X-ray optics. These optics defocus the incoming beam in vertical direction and reshape the intensity distribution from a Gaussian to a more desirable top-hat-shaped profile at the same time. With these optics X-ray full-field imaging of extended objects becomes possible without having to stack several scans or applying a cone beam geometry in order to image the entire specimen. For in situ experiments in general and for diffraction limited sources in particular this gain in field of view and the optimization of the intensity distribution is going to be very beneficial.