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.

Monte Carlo and experimental evaluation of a Timepix4 compact gamma camera for coded aperture nuclear medicine imaging with depth resolution

PURPOSE

We designed a prototype compact gamma camera (MediPROBE4) for nuclear medicine tasks, including radio-guided surgery and sentinel lymph node imaging with a 99mTc radiotracer. We performed Monte Carlo (MC) simulations for image performance assessment, and first spectroscopic imaging tests with a 300 μm thick silicon detector.

METHODS

The hand-held camera (1 kg weight) is based on a Timepix4 readout circuit for photon-counting, energy-sensitive, hybrid pixel detectors (24.6 × 28.2 mm2 sensitive area, 55 μm pixel pitch), developed by the Medipix4 Collaboration. The camera design adopts a CdTe detector (1 or 2 mm thick) bump-bonded to a Timepix4 readout chip and a coded aperture collimator with 0.25 mm diameter round holes made of 3D printed 1-mm thick tungsten. Image reconstruction is performed via autocorrelation deconvolution.

RESULTS

Geant4 MC simulations showed that, for a 99mTc source in air, at 50 mm source-collimator distance, the estimated collimator sensitivity (4 × 10-4) is 292 times larger than that of a single hole in the mask; the system sensitivity is 0.22 cps/kBq (2 mm CdTe); the lateral spatial resolution is 1.7 mm FWHM. The estimated axial longitudinal resolution is 8.2 mm FWHM at 40 mm distance. First experimental tests with a 300 μm thick Silicon pixel detector bump-bonded to a Timepix4 chip and a high-resolution coded aperture collimator showed time-over-threshold and time-of-arrival capabilities with 241Am and 133Ba gamma-ray sources.

CONCLUSIONS

MC simulations and validation lab tests showed the expected performance of the MediPROBE4 compact gamma camera for gamma-ray 3D imaging.

Effect of the local energy distribution of x-ray beams generated through inverse Compton scattering in dual-energy imaging applications

X-ray sources based on the inverse Compton interaction between a laser and a relativistic electron beam are emerging as a promising compact alternative to synchrotron for the production of intense monochromatic and tunable radiation. The emission characteristics enable several innovative imaging techniques, including dual-energy K-edge subtraction (KES) imaging. The performance of these techniques is optimal in the case of perfectly monochromatic x-ray beams, and the implementation of KES was proven to be very effective with synchrotron radiation. Nonetheless, the features of inverse Compton scattering (ICS) sources make them good candidates for a more compact implementation of KES techniques. The energy and intensity distribution of the emitted radiation is related to the emission direction, which means different beam qualities in different spatial positions. In fact, as the polar angle increases, the average energy decreases, while the local energy bandwidth increases and the emission intensity decreases. The scope of this work is to describe the impact of the local energy distribution variations on KES imaging performance. By means of analytical simulations, the reconstructed signal, signal-to-noise ratio, and background contamination were evaluated as a function of the position of each detector pixel. The results show that KES imaging is possible with ICS x-ray beams, even if the image quality slightly degrades at the detector borders for a fixed collimation angle and, in general, as the beam divergence increases. Finally, an approach for the optimization of specific imaging tasks is proposed by considering the characteristics of a given source.

Patient-specific Monte Carlo-based organ dose estimates in spiral CT via optical 3D body scanning and adaptation of a voxelized phantom dataset: proof-of-principle

Objective. We present a method for personalized organ dose estimates obtained before the CT exam, via 3D optical body scanning and Monte Carlo (MC) simulations.Approach. A voxelized phantom is derived by adapting a reference phantom to the body size and shape measured with a portable 3D optical scanner, which returns the 3D silhouette of the patient. This was used as an external rigid envelope for incorporating a tailored version of the internal body anatomy derived from a phantom dataset (National Cancer Institute, NIH, USA) matched for gender, age, weight, and height. The proof-of-principle was conducted on adult head phantoms. The Geant4 MC code provided estimates of the organ doses from 3D absorbed dose maps in voxelized body phantom.Main results. We applied this approach for head CT scanning using an anthropomorphic voxelized head phantom derived from 3D optical scans of mannequins. We compared the estimates of head organ doses with those provided by the NCICT3.0 software (NCI, NIH, USA). Head organ doses differed up to 38% using the proposed personalized estimate and MC code, with respect to corresponding estimates calculated for the standard (non-personalized) reference head phantom. Preliminary application of the MC code to chest CT scans is shown. Real-time pre-exam personalized CT dosimetry is envisaged with adoption of a GPU-based fast MC code.Significance. The developed procedure for personalized organ dose estimates before the CT exam, introduces a new approach for realistic description of size and shape of patients via voxelized phantoms specific for each patient.

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.

Fabrication of 3D printed patient-derived anthropomorphic breast phantoms for mammography and digital breast tomosynthesis: Imaging assessment with clinical X-ray spectra

PURPOSE

To design, fabricate and characterize 3D printed, anatomically realistic, compressed breast phantoms for digital mammography (DM) and digital breast tomosynthesis (DBT) x-ray imaging.

MATERIALS

We realized 3D printed phantoms simulating healthy breasts, via fused deposition modeling (FDM), with a layer resolution of 0.1 mm and 100% infill density, using a dual extruder printer. The digital models were derived from a public dataset of segmented clinical breast computed tomography scans. Three physical phantoms were printed in polyethylene terephthalate (PET), acrylonitrile-butadiene-styrene (ABS), or in polylactic-acid (PLA) materials, using ABS as a substitute for adipose tissue, and PLA or PET filaments for replicating glandular and skin tissues. 3D printed phantoms were imaged at three clinical centers with DM and DBT scanners, using typical spectra. Anatomical noise of the manufactured phantoms was evaluated via the estimates of the β parameter both in DM images and in images acquired via a clinical computed tomography (CT) scanner.

RESULTS

DM and DBT phantom images showed an inner texture qualitatively similar to the images of a clinical DM or DBT exam, suitably reproducing the glandular structure of their computational phantoms. β parameters evaluated in DM images of the manufactured phantoms ranged between 2.84 and 3.79; a lower β was calculated from the CT scan.

CONCLUSIONS

FDM 3D printed compressed breast phantoms have been fabricated using ABS, PLA and PET filaments. DM and DBT images with clinical x-ray spectra showed realistic textures. These phantoms appear promising for clinical applications in quality assurance, image quality and dosimetry assessments.

Comparisons of glandular breast dose between digital mammography, tomosynthesis and breast CT based on anthropomorphic patient-derived breast phantoms

PURPOSE

To evaluate the bias to the mean glandular dose (MGD) estimates introduced by the homogeneous breast models in digital breast tomosynthesis (DBT) and to have an insight into the glandular dose distributions in 2D (digital mammography, DM) and 3D (DBT and breast dedicated CT, BCT) x-ray breast imaging by employing breast models with realistic glandular tissue distribution and organ silhouette.

METHODS

A Monte Carlo software for DM, DBT and BCT simulations was adopted for the evaluation of glandular dose distribution in 60 computational anthropomorphic phantoms. These computational phantoms were derived from 3D breast images acquired via a clinical BCT scanner.

RESULTS

g·c·s·T conversion coefficients based on homogeneous breast model led to a MGD overestimate of 18% in DBT when compared to MGD estimated via anthropomorphic phantoms; this overestimate increased up to 21% for recently computed DgN_DBT conversion coefficients. The standard deviation of the glandular dose distribution in BCT resulted 60% lower than in DM and 55% lower than in DBT. The glandular dose peak - evaluated as the average value over the 5% of the gland receiving the highest dose - is 2.8 times the MGD in DM, this factor reducing to 2.6 and 1.6 in DBT and BCT, respectively.

CONCLUSIONS

Conventional conversion coefficients for MGD estimates based on homogeneous breast models overestimate MGD by 18%, when compared to MGD estimated via anthropomorphic phantoms. The ratio between the peak glandular dose and the MGD is 2.8 in DM. This ratio is 8% and 75% higher than in DBT and BCT, respectively.

Breast radiotherapy with kilovoltage photons and gold nanoparticles as radiosensitizer: An in vitro study

Purpose

We investigated the dose enhancement and internalization of gold nanoparticles (AuNPs) used as a radiosensitizer agent for rotational radiotherapy of breast cancer using a kilovoltage (kV) X‐ray beam.

Methods

Human breast cancer cells MDA‐MB‐231 were incubated with or without 100 μg/mL (4.87 nM) or 200 μg/mL (9.74 nM) 15 nm AuNPs and irradiated with 100 kV, 190 kV, or 6 MV X‐rays. To assess the toxicity of the AuNPs, we performed a Sulforhodamine B assay. Using atomic absorption spectroscopy, scanning electron microscopy, transmission electron microscopy, and time‐lapse optical microscopy (rate of 2 frames per minute), we carried out a quantitative assessment of the amount of gold internalized by MDA‐MB‐231 cells and a characterization of the static and dynamical aspects of this internalization process.

Results

No effect of AuNPs alone was shown on cell viability. Time‐lapse optical microscopy showed for the first time AuNPs cellular uptake and the dynamics of AuNPs internalization. Electron microscopy demonstrated AuNPs localization in endosomal vesicles, preferentially in the perinuclear region. After irradiation at doses up to 2 Gy, cell survival fraction curves showed increased mortality with AuNPs, with respect to irradiation without AuNPs. The highest effect of radioenhancement by AuNPs (at 9.74 nM AuNPs concentration) was observed at 190 kV showing a dose enhancement factor of 1.33 ± 0.06 (1.34 ± 0.02 at 100 kV), while at 6 MV it was 1.14 ± 0.06.

Conclusions

The observed radio‐sensitization effect is promising for future radio‐enhanced kV radiotherapy of breast cancer and quantitatively in the order of previous observations for 15 nm AuNPs. These results of a significant dose enhancement were obtained at 15 nm AuNPs concentration as low as several nanomolar units, at dose levels typical of a single dose fraction in a radiotherapy session. Dynamical behavior of the 3D spatial distribution of 15 nm AuNPs outside the nucleus of single breast cancer cell was observed, with possible implications for future models of AuNPs sensitization.

Radiomics software for breast imaging optimization and simulation studies

BACKGROUND AND OBJECTIVE

The development, control and optimisation of new x-ray breast imaging modalities could benefit from a quantitative assessment of the resulting image textures. The aim of this work was to develop a software tool for routine radiomics applications in breast imaging, which will also be available upon request.

METHODS

The tool (developed in MATLAB) allows image reading, selection of Regions of Interest (ROI), analysis and comparison. Requirements towards the tool also included convenient handling of common medical and simulated images, building and providing a library of commonly applied algorithms and a friendly graphical user interface. Initial set of features and analyses have been selected after a literature search. Being open, the tool can be extended, if necessary.

RESULTS

The tool allows semi-automatic extracting of ROIs, calculating and processing a total of 23 different metrics or features in 2D images and/or in 3D image volumes. Computations of the features were verified against computations with other software packages performed with test images. Two case studies illustrate the applicability of the tool - (i) features on a series of 2D 'left' and 'right' CC mammograms acquired on a Siemens Inspiration system were computed and compared, and (ii) evaluation of the suitability of newly proposed and developed breast phantoms for x-ray-based imaging based on reference values from clinical mammography images. Obtained results could steer the further development of the physical breast phantoms.

CONCLUSIONS

A new image analysis toolbox was realized and can now be used in a multitude of radiomics applications, on both clinical and test images.

Dataset of patient-derived digital breast phantoms for in silico studies in breast computed tomography, digital breast tomosynthesis, and digital mammography

PURPOSE

To present a dataset of computational digital breast phantoms derived from high-resolution three-dimensional (3D) clinical breast images for the use in virtual clinical trials in two-dimensional (2D) and 3D x-ray breast imaging.

ACQUISITION AND VALIDATION METHODS

Uncompressed computational breast phantoms for investigations in dedicated breast CT (BCT) were derived from 150 clinical 3D breast images acquired via a BCT scanner at UC Davis (California, USA). Each image voxel was classified in one out of the four main materials presented in the field of view: fibroglandular tissue, adipose tissue, skin tissue, and air. For the image classification, a semi-automatic software was developed. The semi-automatic classification was compared via manual glandular classification performed by two researchers. A total of 60 compressed computational phantoms for virtual clinical trials in digital mammography (DM) and digital breast tomosynthesis (DBT) were obtained from the corresponding uncompressed phantoms via a software algorithm simulating the compression and the elastic deformation of the breast, using the tissue's elastic coefficient. This process was evaluated in terms of glandular fraction modification introduced by the compression procedure. The generated cohort of 150 uncompressed computational breast phantoms presented a mean value of the glandular fraction by mass of 12.3%; the average diameter of the breast evaluated at the center of mass was 105 mm. Despite the slight differences between the two manual segmentations, the resulting glandular tissue segmentation did not consistently differ from that obtained via the semi-automatic classification. The difference between the glandular fraction by mass before and after the compression was 2.1% on average. The 60 compressed phantoms presented an average glandular fraction by mass of 12.1% and an average compressed thickness of 61 mm.

DATA FORMAT AND ACCESS

The generated digital breast phantoms are stored in DICOM files. Image voxels can present one out of four values representing the different classified materials: 0 for the air, 1 for the adipose tissue, 2 for the glandular tissue, and 3 for the skin tissue. The generated computational phantoms datasets were stored in the Zenodo public repository for research purposes (http://doi.org/10.5281/zenodo.4529852, http://doi.org/10.5281/zenodo.4515360).

POTENTIAL APPLICATIONS

The dataset developed within the INFN AGATA project will be used for developing a platform for virtual clinical trials in x-ray breast imaging and dosimetry. In addition, they will represent a valid support for introducing new breast models for dose estimates in 2D and 3D x-ray breast imaging and as models for manufacturing anthropomorphic physical phantoms.

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.

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.

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.

Dose intercomparison at Italian hadrontherapy centers

PURPOSE

To perform the first dosimetric intercomparison for proton beams in Italy using ionization chambers, according to the IAEA TRS-398 code of practice.

METHODS

Measurement sites included: National Center for Oncological Hadron Therapy (CNAO, Pavia), Center for Proton Therapy (CTP, Trento) and Center for Hadron Therapy and for advanced Nuclear Applications (CATANA, Catania). For comparison we also included a 6 MV photon beam produced at Istituti Clinici Scientifici Maugeri (ICSM, Pavia). For proton beams, both single pseudo-monoenergetic layers (in order to obtain a planned dose of 2 Gy at the reference depth of 2 cm in a water phantom) and Spread-out Bragg peaks (SOBP) have been delivered. Measurements were performed with a PTW Farmer 30010-1 and a PTW Advanced Markus type 34,045 ionization chamber.

RESULTS

Data obtained at CATANA, CNAO and CPT in terms of absorbed dose to water depth show good consistency within the experimental uncertainties, with a weighted mean of 1.99 ± 0.01 Gy and a standard error of 0.003 Gy, with reference to a nominal dose of 2 Gy as designed by the treatment planning system.

CONCLUSIONS

The results showed a standard deviation of less than 1% for single layer and SOBP beams, for all chambers and a percent deviation less than 1.5% for single layer measurements. The weighted means of the absorbed doses for clinical proton beams (118.19 MeV and 173.61 MeV) are consistent within less than 1%. These results agree within the 1.5% difference considered acceptable for national dose intercomparison.

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.

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.

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.

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.