Glandular dose in breast computed tomography with synchrotron radiation

Virtual dosimetry study with three cone-beam breast computed tomography scanners using a fast GPU-based Monte Carlo code

Objective. To compare the dosimetric performance of three cone-beam breast computed tomography (BCT) scanners, using real-time Monte Carlo-based dose estimates obtained with the virtual clinical trials (VCT)-BREAST graphical processing unit (GPU)-accelerated platform dedicated to VCT in breast imaging. Approach. A GPU-based Monte Carlo (MC) code was developed for replicatingin silicothe geometric, x-ray spectra and detector setups adopted, respectively, in two research scanners and one commercial BCT scanner, adopting 80 kV, 60 kV and 49 kV tube voltage, respectively. Our cohort of virtual breasts included 16 anthropomorphic voxelized breast phantoms from a publicly available dataset. For each virtual patient, we simulated exams on the three scanners, up to a nominal simulated mean glandular dose of 5 mGy (primary photons launched, in the order of 1011-1012per scan). Simulated 3D dose maps (recorded for skin, adipose and glandular tissues) were compared for the same phantom, on the three scanners. MC simulations were implemented on a single NVIDIA GeForce RTX 3090 graphics card.Main results.Using the spread of the dose distribution as a figure of merit, we showed that, in the investigated phantoms, the glandular dose is more uniform within less dense breasts, and it is more uniformly distributed for scans at 80 kV and 60 kV, than at 49 kV. A realistic virtual study of each breast phantom was completed in about 3.0 h with less than 1% statistical uncertainty, with 109primary photons processed in 3.6 s computing time.Significance. We reported the first dosimetric study of the VCT-BREAST platform, a fast MC simulation tool for real-time virtual dosimetry and imaging trials in BCT, investigating the dose delivery performance of three clinical BCT scanners. This tool can be adopted to investigate also the effects on the 3D dose distribution produced by changes in the geometrical and spectrum characteristics of a cone-beam BCT scanner.

Multiscale X-ray phase-contrast tomography: From breast CT to micro-CT for virtual histology

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.

Optimization of a customized simultaneous algebraic reconstruction technique algorithm for phase-contrast breast computed tomography

Objective. To introduce the optimization of a customized GPU-based simultaneous algebraic reconstruction technique (cSART) in the field of phase-contrast breast computed tomography (bCT). The presented algorithm features a 3D bilateral regularization filter that can be tuned to yield optimal performance for clinical image visualization and tissues segmentation. Approach. Acquisitions of a dedicated test object and a breast specimen were performed at Elettra, the Italian synchrotron radiation (SR) facility (Trieste, Italy) using a large area CdTe single-photon counting detector. Tomographic images were obtained at 5 mGy of mean glandular dose, with a 32 keV monochromatic x-ray beam in the free-space propagation mode. Three independent algorithms parameters were optimized by using contrast-to-noise ratio (CNR), spatial resolution, and noise texture metrics. The results obtained with the cSART algorithm were compared with conventional SART and filtered back projection (FBP) reconstructions. Image segmentation was performed both with gray scale-based and supervised machine-learning approaches. Main results. Compared to conventional FBP reconstructions, results indicate that the proposed algorithm can yield images with a higher CNR (by 35% or more), retaining a high spatial resolution while preserving their textural properties. Alternatively, at the cost of an increased image ‘patchiness’, the cSART can be tuned to achieve a high-quality tissue segmentation, suggesting the possibility of performing an accurate glandularity estimation potentially of use in the realization of realistic 3D breast models starting from low radiation dose images. Significance. The study indicates that dedicated iterative reconstruction techniques could provide significant advantages in phase-contrast bCT imaging. The proposed algorithm offers great flexibility in terms of image reconstruction optimization, either toward diagnostic evaluation or image segmentation.

Motion artifacts assessment and correction using optical tracking in synchrotron radiation breast CT

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.

Normalized glandular dose coefficients for digital breast tomosynthesis systems with a homogeneous breast model

This work aims at calculating and releasing tabulated values of dose conversion coefficients, DgN_DBT, for mean glandular dose (MGD) estimates in digital breast tomosynthesis (DBT). The DgN_DBT coefficients are proposed as unique conversion coefficients for MGD estimates, in place of dose conversion coefficients in mammography (DgN_DM or c, g, s triad as proposed in worldwide quality assurance protocols) used together with the T correction factor. DgN_DBT is the MGD per unit incident air kerma measured at the breast surface for a 0° projection and the entire tube load used for the scan. The dataset of polyenergetic DgN_DBT coefficients was derived via a Monte Carlo software based on the Geant4 toolkit. Dose coefficients were calculated for a grid of values of breast characteristics (breast thickness in the range 20-90 mm and glandular fraction by mass of 1%, 25%, 50%, 75%, 100%) and the simulated geometries, scan protocols, irradiation geometries and typical spectral qualities replicated those of six commercial DBT systems (GE SenoClaire, Hologic Selenia Dimensions, GE Senographe Pristina, Fujifilm Amulet Innovality, Siemens Mammomat Inspiration and IMS Giotto Class). For given breast characteristics, target/filter combination, tube voltage and half value layer (HVL), two spectra with two HVL values have been simulated in order to permit MGD estimates from experimental HVL values via mathematical interpolation from tabulated values. The adopted breast model assumes homogenous composition of glandular and adipose tissues; it includes a 1.45 mm thick skin envelope in place of the 4-5 mm envelope commonly adopted in dosimetry protocols. The simulation code was validated versus AAPM Task group 195 Monte Carlo reference data sets (absolute differences not higher than 1.1%) and by comparison to relative dosimetry measurements with radiochromic film in a PMMA test object (differences within the maximum experimental uncertainty of 11%). The calculated coefficients show maximum relative deviations of -17.6% and +6.1% from those provided by the DBT dose coefficients adopted in the EUREF protocol and of 1.5%, on average, from data in the AAPM TG223 report. A spreadsheet is provided for interpolating the tabulated DgN_DBT coefficients for arbitrary values of HVL, compressed breast thickness and glandular fraction, in the corresponding investigated ranges, for each DBT unit modeled in this work.

Radiation dosimetry of a clinical prototype dedicated cone-beam breast CT system with offset detector

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 ( DgN CT ) 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 ( f g = 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.75 d ), aligned with the axis of rotation (AOR) located at 65 cm from the focal spot to determine the DgN CT . 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 DgN CT ranged 0.177 - 0.574  mGy/mGy and reduction in DgN CT 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 DgN CT ranged 0.221 - 0.581  mGy/mGy and reduction in DgN CT was observed for all breast diameters. The reduction in DgN CT 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 DgN CT with offset-detector geometries was not dependent on f g . Numerical fits of DgN CT d , l , f g were generated for each geometry.

CONCLUSION

The DgN CT and the numerical fit, D g N CT d , l , f g 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.

Experimental optimization of the energy for breast-CT with synchrotron radiation

Breast Computed Tomography (bCT) is a three-dimensional imaging technique that is raising interest among radiologists as a viable alternative to mammographic planar imaging. In X-rays imaging it would be desirable to maximize the capability of discriminating different tissues, described by the Contrast to Noise Ratio (CNR), while minimizing the dose (i.e. the radiological risk). Both dose and CNR are functions of the X-ray energy. This work aims at experimentally investigating the optimal energy that, at fixed dose, maximizes the CNR between glandular and adipose tissues. Acquisitions of both tissue-equivalent phantoms and actual breast specimens have been performed with the bCT system implemented within the Syrma-3D collaboration at the Syrmep beamline of the Elettra synchrotron (Trieste). The experimental data have been also compared with analytical simulations and the results are in agreement. The CNR is maximized at energies around 26–28 keV. These results are in line with the outcomes of a previously presented simulation study which determined an optimal energy of 28 keV for a large set of breast phantoms with different diameters and glandular fractions. Finally, a study on photon starvation has been carried out to investigate how far the dose can be reduced still having suitable images for diagnostics.

Radiochromic film dosimetry in synchrotron radiation breast computed tomography: a phantom study

This study relates to the INFN project SYRMA-3D for in vivo phase-contrast breast computed tomography using the SYRMEP synchrotron radiation beamline at the ELETTRA facility in Trieste, Italy. This peculiar imaging technique uses a novel dosimetric approach with respect to the standard clinical procedure. In this study, optimization of the acquisition procedure was evaluated in terms of dose delivered to the breast. An offline dose monitoring method was also investigated using radiochromic film dosimetry. Various irradiation geometries have been investigated for scanning the prone patient's pendant breast, simulated by a 14 cm-diameter polymethylmethacrylate cylindrical phantom containing pieces of calibrated radiochromic film type XR-QA2. Films were inserted mid-plane in the phantom, as well as wrapped around its external surface, and irradiated at 38 keV, with an air kerma value that would produce an estimated mean glandular dose of 5 mGy for a 14 cm-diameter 50% glandular breast. Axial scans were performed over a full rotation or over 180°. The results point out that a scheme adopting a stepped rotation irradiation represents the best geometry to optimize the dose distribution to the breast. The feasibility of using a piece of calibrated radiochromic film wrapped around a suitable holder around the breast to monitor the scan dose offline is demonstrated.

Flattening filter for Gaussian-shaped monochromatic X-ray beams: an application to breast computed tomography

The vertical intensity distribution of synchrotron-based X-ray beams usually has a Gaussian profile encompassing large intensity variations. For biomedical imaging applications this may entail sub-optimal dose distributions and large fluctuations in terms of image noise. Commonly, planar metallic filters coupled with absorbing slits systems are applied to adjust the delivered flux and to limit intensity variations, respectively. The latter results in a reduction of the effective beam size. A flattening filter that counterbalances the transverse inhomogeneity, while retaining a sufficient flux, has been developed in the context of a monochromatic phase-contrast breast computed tomography application, ongoing at the Elettra synchrotron facility. The implementation of the new filtration system results in homogeneous intensity (hence dose) distribution and signal-to-noise ratio across the imaged volume. Finally, and most importantly, it allows a wider portion of the beam to be used, directly translating into a major (∼40%) reduction of the overall scan time for samples requiring a field of view larger than the beam size (i.e. multiple translation steps).

Image quality comparison between a phase-contrast synchrotron radiation breast CT and a clinical breast CT: a phantom based study

In this study we compared the image quality of a synchrotron radiation (SR) breast computed tomography (BCT) system with a clinical BCT in terms of contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), noise power spectrum (NPS), spatial resolution and detail visibility. A breast phantom consisting of several slabs of breast-adipose equivalent material with different embedded targets (i.e., masses, fibers and calcifications) was used. Phantom images were acquired using a dedicated BCT system installed at the Radboud University Medical Center (Nijmegen, The Netherlands) and the SR BCT system at the SYRMEP beamline of Elettra SR facility (Trieste, Italy) based on a photon-counting detector. Images with the SR setup were acquired mimicking the clinical BCT conditions (i.e., energy of 30 keV and radiation dose of 6.5 mGy). Images were reconstructed with an isotropic cubic voxel of 273 µm for the clinical BCT, while for the SR setup two phase-retrieval (PhR) kernels (referred to as “smooth” and “sharp”) were alternatively applied to each projection before tomographic reconstruction, with voxel size of 57 × 57 × 50 µm³. The CNR for the clinical BCT system can be up to 2-times higher than SR system, while the SNR can be 3-times lower than SR system, when the smooth PhR is used. The peak frequency of the NPS for the SR BCT is 2 to 4-times higher (0.9 mm⁻¹ and 1.4 mm⁻¹ with smooth and sharp PhR, respectively) than the clinical BCT (0.4 mm⁻¹). The spatial resolution (MTF_10%) was estimated to be 1.3 lp/mm for the clinical BCT, and 5.0 lp/mm and 6.7 lp/mm for the SR BCT with the smooth and sharp PhR, respectively. The smallest fiber visible in the SR BCT has a diameter of 0.15 mm, while for the clinical BCT is 0.41 mm. Calcification clusters with diameter of 0.13 mm are visible in the SR BCT, while the smallest diameter for the clinical BCT is 0.29 mm. As expected, the image quality of the SR BCT outperforms the clinical BCT system, providing images with higher spatial resolution and SNR, and with finer granularity. Nevertheless, this study assesses the image quality gap quantitatively, giving indications on the benefits associated with SR BCT and providing a benchmarking basis for its clinical implementation. In addition, SR-based studies can provide a gold-standard in terms of achievable image quality, constituting an upper-limit to the potential clinical development of a given technique.

Optimization of the energy for Breast monochromatic absorption X-ray Computed Tomography

The limits of mammography have led to an increasing interest on possible alternatives such as the breast Computed Tomography (bCT). The common goal of all X-ray imaging techniques is to achieve the optimal contrast resolution, measured through the Contrast to Noise Ratio (CNR), while minimizing the radiological risks, quantified by the dose. Both dose and CNR depend on the energy and the intensity of the X-rays employed for the specific imaging technique. Some attempts to determine an optimal energy for bCT have suggested the range 22 keV-34 keV, some others instead suggested the range 50 keV-60 keV depending on the parameters considered in the study. Recent experimental works, based on the use of monochromatic radiation and breast specimens, show that energies around 32 keV give better image quality respect to setups based on higher energies. In this paper we report a systematic study aiming at defining the range of energies that maximizes the CNR at fixed dose in bCT. The study evaluates several compositions and diameters of the breast and includes various reconstruction algorithms as well as different dose levels. The results show that a good compromise between CNR and dose is obtained using energies around 28 keV.

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.

Advancements towards the implementation of clinical phase-contrast breast computed tomography at Elettra

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.

Geant4 physics list comparison for the simulation of phase-contrast mammography (XPulse project)

PURPOSE

Breast cancer is the most frequent cancer in women. Early and accurate detection of the disease is a major factor in patient survival. To this end, phase-contrast imaging has gained significant interest in recent years. The aim of this work was to validate the physics models of a Geant4 mammography imaging simulation (in the context of the XPulse project) by comparing to EGSnrc results.

METHODS

We used three Geant4 electromagnetic physics lists of the version 10.4 of the toolkit: Standard, Livermore and Penelope. We calculated energy distributions in homogeneous and inhomogeneous phantoms and breast doses in DICOM images. The simulations used photon beams of energies 20-100 keV. The Geant4 calculations were compared with EGSnrc/DOSXYZnrc simulations.

RESULTS

We found a very good agreement between the Standard Electromagnetic option 4 and Livermore Physics Lists (within 1% for all beam energies). Larger differences were found between Standard Electromagnetic option 4 and Penelope Physics Lists (about 4%). The agreement of longitudinal energy distributions between Geant4 Standard Electromagnetic option 4 and EGSnrc was good in water and light biological materials, but important discrepancies were found in heavy elements. We confirmed with both codes that dose to the breast is minimal at beam energy around 60 keV.

CONCLUSIONS

Overall, we found good agreement between the option 4 of the Standard Electromagnetic physics list and Livermore physics lists of Geant4, as well as EGSnrc for materials relevant to mammography screening. Further investigations are needed for the case of heavier materials.

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.

Large-area single-photon-counting CdTe detector for synchrotron radiation computed tomography: a dedicated pre-processing procedure

Large-area CdTe single-photon-counting detectors are becoming more and more attractive in view of low-dose imaging applications due to their high efficiency, low intrinsic noise and absence of a scintillating screen which affects spatial resolution. At present, however, since the dimensions of a single sensor are small (typically a few cm²), multi-module architectures are needed to obtain a large field of view. This requires coping with inter-module gaps and with close-to-edge pixels, which generally show a non-optimal behavior. Moreover, high-Z detectors often show gain variations in time due to charge trapping: this effect is detrimental especially in computed tomography (CT) applications where a single tomographic image requires hundreds of projections continuously acquired in several seconds. This work has been carried out at the SYRMEP beamline of the Elettra synchrotron radiation facility (Trieste, Italy), in the framework of the SYRMA-3D project, which aims to perform the world's first breast-CT clinical study with synchrotron radiation. An ad hoc data pre-processing procedure has been developed for the PIXIRAD-8 CdTe single-photon-counting detector, comprising an array of eight 30.7 mm × 24.8 mm modules tiling a 246 mm × 25 mm sensitive area, which covers the full synchrotron radiation beam. The procedure consists of five building blocks, namely dynamic flat-fielding, gap seaming, dynamic ring removal, projection despeckling and around-gap equalization. Each block is discussed and compared, when existing, with conventional approaches. The effectiveness of the pre-processing is demonstrated for phase-contrast CT images of a human breast specimen. The dynamic nature of the proposed procedure, which provides corrections dependent upon the projection index, allows the effective removal of time-dependent artifacts, preserving the main image features including phase effects.

Synchrotron radiation external beam rotational radiotherapy of breast cancer: proof of principle

The principle of rotational summation of the absorbed dose for breast cancer treatment with orthovoltage X-ray beams was proposed by J. Boone in 2012. Here, use of X-ray synchrotron radiation for image guided external beam rotational radiotherapy treatment of breast cancer is proposed. Tumor irradiation occurs with the patient in the prone position hosted on a rotating bed, with her breast hanging from a hole in the bed, which rotates around a vertical axis passing through the tumor site. Horizontal collimation of the X-ray beam provides for whole breast or partial breast irradiation, while vertical translation of the bed and successive rotations allow for irradiation of the full tumor volume, with dose rates which permit also hypofractionated treatments. In this work, which follows a previous preliminary report, results are shown of a full series of measurements on polyethylene and acrylic cylindrical phantoms carried out at the Australian Synchrotron, confirmed by Geant4 Monte Carlo simulations, intended to demonstrate the proof of principle of the technique. Dose measurements were carried out with calibrated ion chambers, radiochromic films and thermoluminescence dosimeters. The photon energy investigated was 60 keV. Image guidance may occur with the transmitted beam for contrast-enhanced breast computed tomography. For a horizontal beam collimation of 1.5 cm and rotation around the central axis of a 14 cm-diameter polyethylene phantom, a periphery-to-center dose ratio of 14% was measured. The simulations showed that under the same conditions the dose ratio decreases with increasing photon energy down to 10% at 175 keV. These values are comparable with those achievable with conventional megavoltage radiotherapy of breast cancer with a medical linear accelerator. Dose painting was demonstrated with two off-center `cancer foci' with 1.3 Gy and 0.6 Gy target doses. The use of a radiosensitizing agent for dose enhancement is foreseen.

A proposal for a quality control protocol in breast CT with synchrotron radiation

Background

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.

Materials and methods

A set of projections is acquired and used to perform a CT reconstruction of two selected portions of the phantom. Such portions contain a uniform layer of water and a set of radiographic inserts, respectively. Together, they allow to perform several image quality measurements, namely CT number linearity, reconstruction accuracy, uniformity, noise, and low contrast resolution. All measurements are repeated at different beam energies in the range of interest, and at two different dose values.

Results

Measurements show a good linearity in the soft tissue range, paired to a high accuracy of the CT number reconstruction. Uniformity and noise measurements show that reconstruction inhomogeneities are bound to a few percent of the average pixel values. However, low contrast detectability is limited to the higher portion of the explored energy range.

Conclusions

The results of the measurements are satisfactory in terms of their quality, feasibility and reproducibility. With minimal modifications, the phantom is promising to allow a set of image quality measurements to be used in the upcoming clinical trial.

Evaluation of the BreastSimulator software platform for breast tomography

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.

A Monte Carlo model for mean glandular dose evaluation in spot compression mammography

PURPOSE

To characterize the dependence of normalized glandular dose (DgN) on various breast model and image acquisition parameters during spot compression mammography and other partial breast irradiation conditions, and evaluate alternative previously proposed dose-related metrics for this breast imaging modality.

METHODS

Using Monte Carlo simulations with both simple homogeneous breast models and patient-specific breasts, three different dose-related metrics for spot compression mammography were compared: the standard DgN, the normalized glandular dose to only the directly irradiated portion of the breast (DgNv), and the DgN obtained by the product of the DgN for full field irradiation and the ratio of the mid-height area of the irradiated breast to the entire breast area (DgNM ). How these metrics vary with field-of-view size, spot area thickness, x-ray energy, spot area and position, breast shape and size, and system geometry was characterized for the simple breast model and a comparison of the simple model results to those with patient-specific breasts was also performed.

RESULTS

The DgN in spot compression mammography can vary considerably with breast area. However, the difference in breast thickness between the spot compressed area and the uncompressed area does not introduce a variation in DgN. As long as the spot compressed area is completely within the breast area and only the compressed breast portion is directly irradiated, its position and size does not introduce a variation in DgN for the homogeneous breast model. As expected, DgN is lower than DgNv for all partial breast irradiation areas, especially when considering spot compression areas within the clinically used range. DgNM underestimates DgN by 6.7% for a W/Rh spectrum at 28 kVp and for a 9 × 9 cm2 compression paddle.

CONCLUSION

As part of the development of a new breast dosimetry model, a task undertaken by the American Association of Physicists in Medicine and the European Federation of Organizations of Medical Physics, these results provide insight on how DgN and two alternative dose metrics behave with various image acquisition and model parameters.

Monte Carlo simulation of radiation dose distribution in X-ray imaging at Shanghai Synchrotron Radiation Facility

The aim of this study is to investigate the dose distribution in the process of X-ray imaging at Shanghai Synchrotron Radiation Facility. We used the Monte Carlo software EGSnrc based on the statistics of the beamline BL13W to simulate the actual irradiation environment and calculate the dose distribution. Three different sizes of sources and two different tomographic models were created to predict the dose distribution at different regions of interest (ROI) with different fields of view (FOV). The simulation results showed that the dose in ROI is higher. The dose accumulated in the ROI with smaller FOV was more concentrated and relatively smaller compared with the larger FOV. This calculated dose distribution provides guidance for in-vivo animal experiments and other applications such as microbeam therapy.

Towards breast cancer rotational radiotherapy with synchrotron radiation

PURPOSE

We performed the first investigations, via measurements and Monte Carlo simulations on phantoms, of the feasibility of a new technique for synchrotron radiation rotational radiotherapy for breast cancer (SR³T).

METHODS

A Monte Carlo (MC) code based on Geant4 toolkit was developed in order to simulate the irradiation with the SR³T technique and to evaluate the skin sparing effect in terms of centre-to-periphery dose ratio at different energies in the range 60-175keV. Preliminary measurements were performed at the Australian Synchrotron facility. Radial dose profiles in a 14-cm diameter polyethylene phantom were measured with a 100-mm pencil ionization chamber for different beam sizes and compared with the results of MC simulations. Finally, the dose painting feasibility was demonstrated with measurements with EBT3 radiochromic films in a phantom and collimating the SR beam at 1.5cm in the horizontal direction.

RESULTS

MC simulations showed that the SR³T technique assures a tumour-to-skin absorbed dose ratio from about 7:1 (at 60keV photon energy) to about 10:1 (at 175keV), sufficient for skin sparing during radiotherapy. The comparison between the results of MC simulations and measurements showed an agreement within 5%. Two off-centre foci were irradiated shifting the rotation centre in the horizontal direction.

CONCLUSIONS

The SR³T technique permits to obtain different dose distributions in the target with multiple rotations and can be guided via synchrotron radiation breast computed tomography imaging, in propagation based phase-contrast conditions. Use of contrast agents like iodinated solutions or gold nanoparticles for dose enhancement (DE-SR³T) is foreseen and will be investigated in future work.

Towards breast tomography with synchrotron radiation at Elettra: first images

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 μm)(3) CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR.