Breast phantom with silicone implant for evaluation in conventional mammography

With the increased incidence of cancer and a similarly increased number of surgeries for insertion of silicone breast implants, it is necessary to assess the effect of such material within the breast tissue, particularly in mammography, because of the reduction in the power of breast cancer diagnosis. In this work, we introduce a breast phantom with silicone implants in order to evaluate the influence of the implant on the visibility of the main mammographic findings: fibers, microcalcifications and tumor masses. In this proposed phantom, the breast tissue was simulated using gel paraffin. In the optical density of phantom mammograms with implants, a reduction in breast tissue visibility was seen corresponding to 23% when compared to a phantom without silicone implants. This poor visibility was due to the X-ray beam scattering on silicone material; this effect produced a loss of visibility in the areas adjacent to the implant. It is expected that the proposed phantom model may be used as a device for the establishment of a technical standard for these types of procedures.

Three-dimensional, in-silico breast phantom for multimodality image simulations

We report on the construction of anatomically realistic three-dimensional in-silico breast phantoms with adjustable sizes, shapes and morphologic features. The concept of multiscale spatial resolution is implemented for generating breast tissue images from multiple modalities. Breast epidermal boundary and subcutaneous fat layer is generated by fitting an ellipsoid and 2nd degree polynomials to reconstructive surgical data and ultrasound imaging data. Intraglandular fat is simulated by randomly distributing and orienting adipose ellipsoids within a fibrous region immediately within the dermal layer. Cooper’s ligaments are simulated as fibrous ellipsoidal shells distributed within the subcutaneous fat layer. Individual ductal lobes are simulated following a random binary tree model which is generated based upon probabilistic branching conditions described by ramification matrices, as originally proposed by Bakic et al [3, 4]. The complete ductal structure of the breast is simulated from multiple lobes that extend from the base of the nipple and branch towards the chest wall. As lobe branching progresses, branches are reduced in height and radius and terminal branches are capped with spherical lobular clusters. Biophysical parameters are mapped onto the complete anatomical model and synthetic multimodal images (Mammography, Ultrasound, CT) are generated for phantoms of different adipose percentages (40%, 50%, 60%, and 70%) and are analytically compared with clinical examples. Results demonstrate that the in-silico breast phantom has applications in imaging performance evaluation and, specifically, great utility for solving image registration issues in multimodality imaging.

STIR, SPIR and SPAIR techniques in magnetic resonance of the breast: a comparative study

The amount of fat is a component that complicates the clinical evaluation and the differential diagnostic between benign and malign lesions in the breast MRI examinations. To overcome this problem, an effective erasing of the fat signal over the images acquisition process, is essentials. This study aims to compare three fat suppression techniques (STIR, SPIR, SPAIR) in the MR images of the breast and to evaluate the best image quality regarding its clinical usefulness. To mimic breast women, a breast phantom was constructed. First the exterior contour and, in second time, its content which was selected based on 7 samples with different components. Finally it was undergone to a MRI breast protocol with the three different fat saturation techniques. The examinations were performed on a 1.5 T MRI system (Philips®). A group of 5 experts evaluated 9 sequences, 3 of each with fat suppression techniques, in which the frequency offset and TI (Inversion Time) were the variables changed. This qualitative image analysis was performed according 4 parameters (saturation uniformity, saturation efficacy, detail of the anatomical structures and differentiation between the fibroglandular and adipose tissue), using a five-point Likert scale. The statistics analysis showed that anyone of the fat suppression techniques demonstrated significant differences compared to the others with (p > 0.05) and regarding each parameter independently. By Fleiss’ kappa coefficient there was a good agreement among observers P(e) = 0.68. When comparing STIR, SPIR and SPAIR techniques it was confirmed that all of them have advantages in the study of the breast MRI. For the studied parameters, the results through the Friedman Test showed that there are similar advantages applying anyone of these techniques.

Dielectric Studies of Corn Syrup for Applications in Microwave Breast Imaging

Permittivity and conductivity studies of corn syrup in various concentrations are performed using coaxial cavity perturbation technique over a frequency range of 250 MHz–3000 MHz. The results are utilized to estimate relaxation time and dipole moments of the samples. The stability of the material over the variations of time is studied. The measured specific absorption rate of the material complies with the microwave power absorption rate of biological tissues. This suggests the feasibility of using corn syrup as a suitable, cost effective coupling medium for microwave breast imaging. The material can also be used as an efficient breast phantom in microwave breast imaging studies.

Refraction corrected transmission ultrasound computed tomography for application in breast imaging

Purpose: We present an iterative framework for CT reconstruction from transmission ultrasound data which accurately and efficiently models the strong refraction effects that occur in our target application: Imaging the female breast. Methods: Our refractive ray tracing framework has its foundation in the fast marching method (FNMM) and it allows an accurate as well as efficient modeling of curved rays. We also describe a novel regularization scheme that yields further significant reconstruction quality improvements. A final contribution is the development of a realistic anthropomorphic digital breast phantom based on the NIH Visible Female data set. Results: Our system is able to resolve very fine details even in the presence of significant noise, and it reconstructs both sound speed and attenuation data. Excellent correspondence with a traditional, but significantly more computationally expensive wave equation solver is achieved. Conclusions: Apart from the accurate modeling of curved rays, decisive factors have also been our regularization scheme and the high-quality interpolation filter we have used. An added benefit of our framework is that it accelerates well on GPUs where we have shown that clinical 3D reconstruction speeds on the order of minutes are possible.

Performance analysis of reconstruction algorithms for breast microwave imaging

The problem of localizing a scatterer, which represents a tumor, in a homogeneous circular domain, which represents a breast, is addressed. A breast imaging method based on microwaves is considered. The microwave imaging involves to several techniques for detecting, localizing and characterizing tumors in breast tissues. In all such methods an electromagnetic inverse scattering problem exists. For the scattering detection method, an algorithm based on a linear procedure solution, inspired by MUltiple SIgnal Classification algorithm (MUSIC) and Time Reversal method (TR), is implemented. The algorithm returns a reconstructed image of the investigation domain in which it is detected the scatterer position. This image is called pseudospectrum. A preliminary performance analysis of the algorithm vying the working frequency is performed: the resolution and the signal-to-noise ratio of the pseudospectra are improved if a multi-frequency approach is considered. The Geometrical Mean-MUSIC algorithm (GM- MUSIC) is proposed as multi-frequency method. The performance of the GMMUSIC is tested in different real life computer simulations. The performed analysis shows that the algorithm detects the scatterer until the electrical parameters of the breast are known. This is an evident limit, since, in a real life situation, the anatomy of the breast is unknown. An improvement in GM-MUSIC is proposed: the Eye-GMMUSIC algorithm. Eye-GMMUSIC algorithm needs no a priori information on the electrical parameters of the breast. It is an optimizing algorithm based on the pattern search algorithm: it searches the breast parameters which minimize the Signal-to-Clutter Mean Ratio (SCMR) in the signal. Finally, the GM-MUSIC and the Eye-GMMUSIC algorithms are tested on a microwave breast cancer detection system consisting of an dipole antenna, a Vector Network Analyzer and a novel breast phantom built at University of Bologna. The reconstruction of the experimental data confirm the GM-MUSIC ability to localize a scatterer in a homogeneous medium.

Tomographic 3D reconstruction through linear inversion for breast microwave imaging

In the present thesis we address the problem of detecting and localizing a small spherical target with characteristic electrical properties inside a volume of cylindrical shape, representing female breast, with MWI. One of the main works of this project is to properly extend the existing linear inversion algorithm from planar slice to volume reconstruction; results obtained, under the same conditions and experimental setup are reported for the two different approaches. Preliminar comparison and performance analysis of the reconstruction algorithms is performed via numerical simulations in a software-created environment: a single dipole antenna is used for illuminating the virtual breast phantom from different positions and, for each position, the corresponding scattered field value is registered. Collected data are then exploited in order to reconstruct the investigation domain, along with the scatterer position, in the form of image called pseudospectrum. During this process the tumor is modeled as a dielectric sphere of small radius and, for electromagnetic scattering purposes, it's treated as a point-like source. To improve the performance of reconstruction technique, we repeat the acquisition for a number of frequencies in a given range: the different pseudospectra, reconstructed from single frequency data, are incoherently combined with MUltiple SIgnal Classification (MUSIC) method which returns an overall enhanced image. We exploit multi-frequency approach to test the performance of 3D linear inversion reconstruction algorithm while varying the source position inside the phantom and the height of antenna plane. Analysis results and reconstructed images are then reported. Finally, we perform 3D reconstruction from experimental data gathered with the acquisition system in the microwave laboratory at DIFA, University of Bologna for a recently developed breast-phantom prototype; obtained pseudospectrum and performance analysis for the real model are reported.

Planification de traitements de curiethérapie du sein à l’aide de l’imagerie par résonance magnétique

Ce mémoire présente l’étude de la faisabilité de la planification de traitements pour la curiethérapie interstitielle du sein en utilisant l’imagerie par résonance magnétique (IRM) seule au lieu de l’imagerie par tomodensitométrie (CT). L'imagerie CT étant la référence, la mesure des différences observables sur la distribution de doses provenant des deux types d’imagerie a été effectuée. Des fantômes de seins ont été fabriqués et utilisés, ainsi que l’imagerie de patients. La taille des fantômes en fonction du positionnement dans l’appareil d’IRM et la longueur reconstruite des cathéters ont été analysées. Les différences dans les distributions de doses de fantômes et de patients ont été calculées en s’assurant que la reconstruction des cathéters provenant des images CT et IRM est la seule variable. La différence dans les critères de doses à respecter est plus grande lorsque la taille du fantôme et/ou un déplacement latéral dans l’IRM sont plus grands. La longueur reconstruite des cathéters est comparable entre les deux techniques d’imagerie. Pour le petit fantôme des différences <2% ont été observées pour tous les critères de dose. Pour le grand fantôme et pour les patients, une valeur maximale de 5% est observée pour les critères sur la cible, mais peut atteindre 19% pour le critère Externe V150%/V100% pour le grand fantôme et 33% pour les patients. Par contre, le seuil clinique de ce critére est toujours respecté. Ceci nous indique que pour la plupart des patients, la zone à traiter serait bien couverte en utilisant les images IRM uniquement pour planifier.

A peripheral dose calculation model for estimating second cancer risk after breast radiotherapy

AbstractBreast cancer is one of the most common cancers affecting one in eight women during their lives. Survival rates have increased steadily thanks to early diagnosis with mammography screening and more efficient treatment strategies. Post-operative radiation therapy is a standard of care in the management of breast cancer and has been shown to reduce efficiently both local recurrence rate and breast cancer mortality. Radiation therapy is however associated with some late effects for long-term survivors. Radiation-induced secondary cancer is a relatively rare but severe late effect of radiation therapy. Currently, radiotherapy plans are essentially optimized to maximize tumor control and minimize late deterministic effects (tissue reactions) that are mainly associated with high doses (» 1 Gy). With improved cure rates and new radiation therapy technologies, it is also important to evaluate and minimize secondary cancer risks for different treatment techniques. This is a particularly challenging task due to the large uncertainties in the dose-response relationship.In contrast with late deterministic effects, secondary cancers may be associated with much lower doses and therefore out-of-field doses (also called peripheral doses) that are typically inferior to 1 Gy need to be determined accurately. Out-of-field doses result from patient scatter and head scatter from the treatment unit. These doses are particularly challenging to compute and we characterized it by Monte Carlo (MC) calculation. A detailed MC model of the Siemens Primus linear accelerator has been thoroughly validated with measurements. We investigated the accuracy of such a model for retrospective dosimetry in epidemiological studies on secondary cancers. Considering that patients in such large studies could be treated on a variety of machines, we assessed the uncertainty in reconstructed peripheral dose due to the variability of peripheral dose among various linac geometries. For large open fields (> 10x10 cm2), the uncertainty would be less than 50%, but for small fields and wedged fields the uncertainty in reconstructed dose could rise up to a factor of 10. It was concluded that such a model could be used for conventional treatments using large open fields only.The MC model of the Siemens Primus linac was then used to compare out-of-field doses for different treatment techniques in a female whole-body CT-based phantom. Current techniques such as conformai wedged-based radiotherapy and hybrid IMRT were investigated and compared to older two-dimensional radiotherapy techniques. MC doses were also compared to those of a commercial Treatment Planning System (TPS). While the TPS is routinely used to determine the dose to the contralateral breast and the ipsilateral lung which are mostly out of the treatment fields, we have shown that these doses may be highly inaccurate depending on the treatment technique investigated. MC shows that hybrid IMRT is dosimetrically similar to three-dimensional wedge-based radiotherapy within the field, but offers substantially reduced doses to out-of-field healthy organs.Finally, many different approaches to risk estimations extracted from the literature were applied to the calculated MC dose distribution. Absolute risks varied substantially as did the ratio of risk between two treatment techniques, reflecting the large uncertainties involved with current risk models. Despite all these uncertainties, the hybrid IMRT investigated resulted in systematically lower cancer risks than any of the other treatment techniques. More epidemiological studies with accurate dosimetry are required in the future to construct robust risk models. In the meantime, any treatment strategy that reduces out-of-field doses to healthy organs should be investigated. Electron radiotherapy might offer interesting possibilities with this regard.RésuméLe cancer du sein affecte une femme sur huit au cours de sa vie. Grâce au dépistage précoce et à des thérapies de plus en plus efficaces, le taux de guérison a augmenté au cours du temps. La radiothérapie postopératoire joue un rôle important dans le traitement du cancer du sein en réduisant le taux de récidive et la mortalité. Malheureusement, la radiothérapie peut aussi induire des toxicités tardives chez les patients guéris. En particulier, les cancers secondaires radio-induits sont une complication rare mais sévère de la radiothérapie. En routine clinique, les plans de radiothérapie sont essentiellement optimisées pour un contrôle local le plus élevé possible tout en minimisant les réactions tissulaires tardives qui sont essentiellement associées avec des hautes doses (» 1 Gy). Toutefois, avec l'introduction de différentes nouvelles techniques et avec l'augmentation des taux de survie, il devient impératif d'évaluer et de minimiser les risques de cancer secondaire pour différentes techniques de traitement. Une telle évaluation du risque est une tâche ardue étant donné les nombreuses incertitudes liées à la relation dose-risque.Contrairement aux effets tissulaires, les cancers secondaires peuvent aussi être induits par des basses doses dans des organes qui se trouvent hors des champs d'irradiation. Ces organes reçoivent des doses périphériques typiquement inférieures à 1 Gy qui résultent du diffusé du patient et du diffusé de l'accélérateur. Ces doses sont difficiles à calculer précisément, mais les algorithmes Monte Carlo (MC) permettent de les estimer avec une bonne précision. Un modèle MC détaillé de l'accélérateur Primus de Siemens a été élaboré et validé avec des mesures. La précision de ce modèle a également été déterminée pour la reconstruction de dose en épidémiologie. Si on considère que les patients inclus dans de larges cohortes sont traités sur une variété de machines, l'incertitude dans la reconstruction de dose périphérique a été étudiée en fonction de la variabilité de la dose périphérique pour différents types d'accélérateurs. Pour de grands champs (> 10x10 cm ), l'incertitude est inférieure à 50%, mais pour de petits champs et des champs filtrés, l'incertitude de la dose peut monter jusqu'à un facteur 10. En conclusion, un tel modèle ne peut être utilisé que pour les traitements conventionnels utilisant des grands champs.Le modèle MC de l'accélérateur Primus a été utilisé ensuite pour déterminer la dose périphérique pour différentes techniques dans un fantôme corps entier basé sur des coupes CT d'une patiente. Les techniques actuelles utilisant des champs filtrés ou encore l'IMRT hybride ont été étudiées et comparées par rapport aux techniques plus anciennes. Les doses calculées par MC ont été comparées à celles obtenues d'un logiciel de planification commercial (TPS). Alors que le TPS est utilisé en routine pour déterminer la dose au sein contralatéral et au poumon ipsilatéral qui sont principalement hors des faisceaux, nous avons montré que ces doses peuvent être plus ou moins précises selon la technTque étudiée. Les calculs MC montrent que la technique IMRT est dosimétriquement équivalente à celle basée sur des champs filtrés à l'intérieur des champs de traitement, mais offre une réduction importante de la dose aux organes périphériques.Finalement différents modèles de risque ont été étudiés sur la base des distributions de dose calculées par MC. Les risques absolus et le rapport des risques entre deux techniques de traitement varient grandement, ce qui reflète les grandes incertitudes liées aux différents modèles de risque. Malgré ces incertitudes, on a pu montrer que la technique IMRT offrait une réduction du risque systématique par rapport aux autres techniques. En attendant des données épidémiologiques supplémentaires sur la relation dose-risque, toute technique offrant une réduction des doses périphériques aux organes sains mérite d'être étudiée. La radiothérapie avec des électrons offre à ce titre des possibilités intéressantes.

A nonlinear elasticity phantom containing spherical inclusions

The strain image contrast of some in vivo breast lesions changes with increasing applied load. This change is attributed to differences in the nonlinear elastic properties of the constituent tissues suggesting some potential to help classify breast diseases by their nonlinear elastic properties. A phantom with inclusions and long-term stability is desired to serve as a test bed for nonlinear elasticity imaging method development, testing, etc. This study reports a phantom designed to investigate nonlinear elastic properties with ultrasound elastographic techniques. The phantom contains four spherical inclusions and was manufactured from a mixture of gelatin, agar and oil. The phantom background and each of the inclusions have distinct Young's modulus and nonlinear mechanical behavior. This phantom was subjected to large deformations (up to 20%) while scanning with ultrasound, and changes in strain image contrast and contrast-to-noise ratio between inclusion and background, as a function of applied deformation, were investigated. The changes in contrast over a large deformation range predicted by the finite element analysis (FEA) were consistent with those experimentally observed. Therefore, the paper reports a procedure for making phantoms with predictable nonlinear behavior, based on independent measurements of the constituent materials, and shows that the resulting strain images (e. g., strain contrast) agree with that predicted with nonlinear FEA.

Tumors in the contralateral breast following radiation therapy for breast cancer

The proportional distribution of independent malignant tumors in the contralateral breast following treatment for breast cancer was investigated to assess the influence of scattered radiation as a cause of these tumors. In a population of 172 patients the proportion of contralateral tumors in each quadrant and the center (the nipple-areolar complex) was compared with the expected, or natural, distribution found in the general population, in the absence of radiation. The observed/expected ratio for contralateral tumors was 1.43 for the upper-inner quadrant; 0.97, lower-inner quadrant; 1.51, center; 0.76, upper-outer quadrant; and 0.64, lower-outer quadrant. In each quadrant, except the lower-inner, the observed/expected ratio differed from 1.00 with statistical significance at the 5% level (one-tail). The same analysis, stratified by age and menopausal status, showed a similar shift of tumors, with more than expected in the inner quadrants and center and less than expected in the outer quadrants, although the results did not show statistical significance at the 5% level for all strata. For each patient the mean absorbed radiation dose for each quadrant and center of the breast was estimated, based on measurements in a tissue-equivalent phantom. Among patients the doses ranged from 0.5 to 8 Gy; within individuals, doses to the inner quadrants typically were a factor of three times higher than doses to the outer quadrants. The results suggest that radiation may be a risk factor for contralateral breast tumors and warrants further investigation. ^

FDTD study of an UWB radar technique for breast tumor detection

A finite difference time domain (FDTD) method is applied to investigate capabilities of an ultra-wide band (UWB) radar system to detect a breast tumor. The first part of the investigations concerns FDTD simulations of a phantom formed by a plastic container with liquid and a small reflecting target. The second part focuses on a three-dimensional numerical breast model with a small tumor. FDTD simulations are carried out assuming a planar incident wave. Various time snap shots of the electromagnetic field are recorded to learn about the physical phenomenon of reflection and scattering in different layers of the phantom.

Clinical translation of a hand-held optical imager for breast cancer diagnostics: In vitro and in vivo tomography studies

Optical imaging is an emerging technology towards non-invasive breast cancer diagnostics. In recent years, portable and patient comfortable hand-held optical imagers are developed towards two-dimensional (2D) tumor detections. However, these imagers are not capable of three-dimensional (3D) tomography because they cannot register the positional information of the hand-held probe onto the imaged tissue. A hand-held optical imager has been developed in our Optical Imaging Laboratory with 3D tomography capabilities, as demonstrated from tissue phantom studies. The overall goal of my dissertation is towards the translation of our imager to the clinical setting for 3D tomographic imaging in human breast tissues. A systematic experimental approach was designed and executed as follows: (i) fast 2D imaging, (ii) coregistered imaging, and (iii) 3D tomographic imaging studies. (i) Fast 2D imaging was initially demonstrated in tissue phantoms (1% Liposyn solution) and in vitro (minced chicken breast and 1% Liposyn). A 0.45 cm3 fluorescent target at 1:0 contrast ratio was detectable up to 2.5 cm deep. Fast 2D imaging experiments performed in vivo with healthy female subjects also detected a 0.45 cm3 fluorescent target superficially placed ∼2.5 cm under the breast tissue. (ii) Coregistered imaging was automated and validated in phantoms with ∼0.19 cm error in the probe’s positional information. Coregistration also improved the target depth detection to 3.5 cm, from multi-location imaging approach. Coregistered imaging was further validated in-vivo , although the error in probe’s positional information increased to ∼0.9 cm (subject to soft tissue deformation and movement). (iii) Three-dimensional tomography studies were successfully demonstrated in vitro using 0.45 cm3 fluorescence targets. The feasibility of 3D tomography was demonstrated for the first time in breast tissues using the hand-held optical imager, wherein a 0.45 cm3 fluorescent target (superficially placed) was recovered along with artifacts. Diffuse optical imaging studies were performed in two breast cancer patients with invasive ductal carcinoma. The images showed greater absorption at the tumor cites (as observed from x-ray mammography, ultrasound, and/or MRI). In summary, my dissertation demonstrated the potential of a hand-held optical imager towards 2D breast tumor detection and 3D breast tomography, holding a promise for extensive clinical translational efforts.