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.

Quantitative comparison between the commercial software STRATOS(®) by Philips and a homemade software for voxel-dosimetry in radiopeptide therapy.

BACKGROUND: Targeted radionuclide therapy is a rapidly growing modality. A few commercial treatment planning systems are entering the market. However, some in-house systems are currently developed for a more flexible and customized dosimetry calculation at voxel-level. For this purpose, we developed a novel software, VoxelMed, and performed a comparison with the software STRATOS. METHODS: The validation of both of them was undertaken using radioactive phantoms with different volume inserts. A cohort of 10 patients was also studied after a therapeutic administration of (177)Lu-labelled radiopeptides. The activity, number of disintegrations, absorbed dose and dose-volume histogram (DVH) were calculated for the phantoms and the kidneys in patients, which were the main critical organs at risk in this study. RESULTS: In phantoms the absorbed doses computed with VoxelMed and STRATOS agree within 5%. In patients at the voxel-level the absorbed dose to kidneys (VoxelMed: mean 0.66 Gy/GBq) showed a limited difference of 5%, but with a remarkable range (-40%, +60%) between the two software packages. Voxel-dosimetry allows to estimate the dose non-homogeneities in volumes, which may be evaluated through DVHs. CONCLUSION: This study demonstrates that a fully 3D voxel-dosimetry with multiple SPECT images is feasible by using home-made or commercial software package and absorbed dose results obtained are similar. The main difference between the studied tools was observed in the activity integration method (effective vs physical half-time to time activity curve tail). We believe that an effective half-time integration method produces a more accurate approximation of clinical uptake and resultant dosimetry.

A free database of radionuclide voxel S values for the dosimetry of nonuniform activity distributions.

The increasing availability of SPECT/CT devices with advanced technology offers the opportunity for the accurate assessment of the radiation dose to the biological target volume during radionuclide therapy. Voxel dosimetry can be performed employing direct Monte Carlo radiation transport simulations, based on both morphological and functional images of the patient. On the other hand, for voxel dosimetry calculations the voxel S value method can be considered an easier approach than patient-specific Monte Carlo simulations, ensuring a good dosimetric accuracy at least for anatomic regions which are characterized by uniform density tissue. However, this approach has been limited because of the lack of tabulated S values for different voxel dimensions and radionuclides. The aim of this work is to provide a free dataset of values which can be used for voxel dosimetry in targeted radionuclide studies. Seven different radionuclides (89Sr, 90Y, 131I, 153Sm, 177Lu, 186Re, 188Re), and 13 different voxel sizes (2.21, 2.33, 2.4, 3, 3.59, 3.9, 4, 4.42, 4.8, 5, 6, 6.8 and 9.28 mm) are considered. Voxel S values are calculated performing simulations of monochromatic photon and electron sources in two different homogeneous tissues (soft tissue and bone) with DOSXYZnrc code, and weighting the contributions on the basis of the radionuclide emission spectra. The outcomes are validated by comparison with Monte Carlo simulations obtained with other codes (PENELOPE and MCNP4c) performing direct simulation of the radionuclide emission spectra. The differences among the different Monte Carlo codes are of the order of a few per cent when considering the source voxel and the bremsstrahlung tail, whereas the highest differences are observed at a distance close to the maximum continuous slowing down approximation range of electrons. These discrepancies would negligibly affect dosimetric assessments. The dataset of voxel S values can be freely downloaded from the website www.medphys.it.

Physical and clinical comparison between a screen-film system and a dual-side reading mammography-dedicated computed radiography system.

BACKGROUND: Digital mammography systems, thanks to a physical performance better than conventional screen-film units, have the potential of reducing the dose to patients, without decreasing the diagnostic accuracy. PURPOSE: To achieve a physical and clinical comparison between two systems: a screen-film plate and a dual-side computed radiography system (CRM; FUJIFILM FCR 5000 MA). MATERIAL AND METHODS: A unique feature of the FCR 5000 MA system is that it has a clear support medium, allowing light emitted during the scanning process to be detected on the "back" of the storage phosphor plate, considerably improving the system's efficiency. The system's physical performance was tested by means of a quantitative analysis, with calculation of the modulation transfer function, detective quantum efficiency, and contrast-detail analysis; subsequently, the results were compared with those achieved using a screen-film system (SFM; Eastmann Kodak MinR-MinR 2000). A receiver operating characteristic (ROC) analysis was then performed on 120 paired clinical images obtained in a craniocaudal projection with the conventional SFM system under standard exposure conditions and also with the CRM system working with a dose reduced by 35% (average breast thickness: 4.3 cm; mean glandular dose: 1.45 mGy). CRM clinical images were interpreted both in hard copy and in soft copy. RESULTS: The ROC analysis revealed that the performances of the two systems (SFM and CRM with reduced dose) were similar (P>0.05): the diagnostic accuracy of the two systems, when valued in terms of the area underneath the ROC curve, was found to be 0.74 for the SFM, 0.78 for the CRM (hard copy), and 0.79 for the CRM (soft copy). CONCLUSION: The outcome obtained from our experiments shows that the use of the dual-side CRM system is a very good alternative to the screen-film system.

Differences among Monte Carlo codes in the calculations of voxel S values for radionuclide targeted therapy and analysis of their impact on absorbed dose evaluations.

Several updated Monte Carlo (MC) codes are available to perform calculations of voxel S values for radionuclide targeted therapy. The aim of this work is to analyze the differences in the calculations obtained by different MC codes and their impact on absorbed dose evaluations performed by voxel dosimetry. Voxel S values for monoenergetic sources (electrons and photons) and different radionuclides (90Y, 131I, and 188Re) were calculated. Simulations were performed in soft tissue. Three general-purpose MC codes were employed for simulating radiation transport: MCNP4C, EGSnrc, and GEANT4. The data published by the MIRD Committee in Pamphlet No. 17, obtained with the EGS4 MC code, were also included in the comparisons. The impact of the differences (in terms of voxel S values) among the MC codes was also studied by convolution calculations of the absorbed dose in a volume of interest. For uniform activity distribution of a given radionuclide, dose calculations were performed on spherical and elliptical volumes, varying the mass from 1 to 500 g. For simulations with monochromatic sources, differences for self-irradiation voxel S values were mostly confined within 10% for both photons and electrons, but with electron energy less than 500 keV, the voxel S values referred to the first neighbor voxels showed large differences (up to 130%, with respect to EGSnrc) among the updated MC codes. For radionuclide simulations, noticeable differences arose in voxel S values, especially in the bremsstrahlung tails, or when a high contribution from electrons with energy of less than 500 keV is involved. In particular, for 90Y the updated codes showed a remarkable divergence in the bremsstrahlung region (up to about 90% in terms of voxel S values) with respect to the EGS4 code. Further, variations were observed up to about 30%, for small source-target voxel distances, when low-energy electrons cover an important part of the emission spectrum of the radionuclide (in our case, for 131I). For 90Y and 188Re, the differences among the various codes have a negligible impact (within few percents) on convolution calculations of the absorbed dose; thus either one of the MC programs is suitable to produce voxel S values for radionuclide targeted therapy dosimetry. However, if a low-energy beta-emitting radionuclide is considered, these differences can affect also dose depositions at small source-target voxel distances, leading to more conspicuous variations (about 9% for 1311) when calculating the absorbed dose in the volume of interest.

Experimental validation of the filtering approach for dose monitoring in proton therapy at low energy.

The higher physical selectivity of proton therapy demands higher accuracy in monitoring of the delivered dose, especially when the target volume is located next to critical organs and a fractionated therapy is applied. A method to verify a treatment plan and to ensure the high quality of the hadrontherapy is to use Positron Emission Tomography (PET), which takes advantage of the nuclear reactions between protons and nuclei in the tissue during irradiation producing beta(+)-emitting isotopes. Unfortunately, the PET image is not directly proportional to the delivered radiation dose distribution; this is the reason why, at the present time, the verification of depth dose profiles with PET techniques is limited to a comparison between the measured activity and the one predicted for the planned treatment by a Monte Carlo model. In this paper we test the feasibility of a different scheme, which permits to reconstruct the expected PET signal from the planned radiation dose distribution along beam direction in a simpler and more direct way. The considered filter model, based on the description of the PET image as a convolution of the dose distribution with a filter function, has already demonstrated its potential applicability to beam energies above 70 MeV. Our experimental investigation provides support to the possibility of extending the same approach to the lower energy range ([40, 70] MeV), in the perspective of its clinical application in eye proton therapy.

An SVM classifier to separate false signals from microcalcifications in digital mammograms.

In this paper we investigate the feasibility of using an SVM (support vector machine) classifier in our automatic system for the detection of clustered microcalcifications in digital mammograms. SVM is a technique for pattern recognition which relies on the statistical learning theory. It minimizes a function of two terms: the number of misclassified vectors of the training set and a term regarding the generalization classifier capability. We compare the SVM classifier with an MLP (multi-layer perceptron) in the false-positive reduction phase of our detection scheme: a detected signal is considered either microcalcification or false signal, according to the value of a set of its features. The SVM classifier gets slightly better results than the MLP one (Az value of 0.963 against 0.958) in the presence of a high number of training data; the improvement becomes much more evident (Az value of 0.952 against 0.918) in training sets of reduced size. Finally, the setting of the SVM classifier is much easier than the MLP one.