Kristina Bliznakova

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.

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 51u2009keV) 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)u2009u2009=u2009u2009α/fu2009u2009 β . 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.

EUTEMPE-RX, an EC supported FP7 project for the training and education of medical physics experts in radiology

The core activity of the medical physics expert (MPE) is to ensure optimal use of ionising radiation in healthcare. It is essential that these healthcare professionals are trained to the highest level, defined as European Qualifications Framework for Lifelong Learning (EQF) level 8 by the European Commissions Radiation Protection Report 174 Guidelines on the MPE. The main objective of the EUTEMPE-RX project is to provide a model training scheme that allows the medical physicist in diagnostic and interventional radiology (D&IR) to reach this high level. A European network of partners was brought together in this FP7 EC project to ensure sufficient expertise in all aspects of the subject and to create a harmonised course programme. Targeted participants are medical physicists in D&IR in hospitals, engineers and scientists in medical device industries and officers working in regulatory authorities. Twelve course modules will be developed at EQF level 8, with radiation safety and diagnostic effectiveness being prevalent subjects. The modules will combine online with face-to-face teaching using a blended learning approach.

Breast tomosynthesis with monochromatic beams: a feasibility study using Monte Carlo simulations

The aim of this study is to investigate the impact on image quality of using monochromatic beams for lower dose breast tomosynthesis (BT). For this purpose, modeling and simulation of BT and mammography imaging processes have been performed using two x-ray beams: one at 28u2009kVp and a monochromatic one at 19u2009keV at different entrance surface air kerma ranging between 0.16 and 5.5u2009mGy. Two 4u2009cm thick computational breast models, in a compressed state, were used: one simple homogeneous and one heterogeneous based on CT breast images, with compositions of 50% glandular-50% adipose and 40% glandular-60% adipose tissues by weight, respectively. Modeled lesions, representing masses and calcifications, were inserted within these breast phantoms. X-ray transport in the breast models was simulated with previously developed and validated Monte Carlo application. Results showed that, for the same incident photon fluence, the use of the monochromatic beam in BT resulted in higher image quality compared to the one using polychromatic acquisition, especially in terms of contrast. For the homogenous phantom, the improvement ranged between 15% and 22% for calcifications and masses, respectively, while for the heterogeneous one this improvement was in the order of 33% for the masses and 17% for the calcifications. For different exposures, comparable image quality in terms of signal-difference-to-noise ratio and higher contrast for all features was obtained when using a monochromatic 19u2009keV beam at a lower mean glandular dose, compared to the polychromatic one. Monochromatic images also provide better detail and, in combination with BT, can lead to substantial improvement in visualization of features, and particularly better edge detection of low-contrast masses.

Power spectrum analysis of the x-ray scatter signal in mammography and breast tomosynthesis projections

PURPOSEnTo analyze the frequency domain characteristics of the signal in mammography images and breast tomosynthesis projections with patient tissue texture due to detected scattered x-rays.nnnMETHODSnAcquisitions of x-ray projection images of 19 different patient breasts were simulated using previously acquired volumetric patient images. Acquisition of these images was performed with a dedicated breast CT prototype system, and the images were classified into voxels representing skin, adipose, and glandular tissue with a previously validated automated algorithm. The classified three dimensional images then underwent simulated mechanical compression representing that which is performed during acquisition of mammography and breast tomosynthesis images. The acquisition of projection images of each patient breast was simulated using Monte Carlo methods with each simulation resulting in two images: one of the primary (non-scattered) signal and one of the scatter signal. To analyze the scatter signal for both mammography and breast tomosynthesis, two projections images of each patient breast were simulated, one with the x-ray source positioned at 0° (mammography and central tomosynthesis projection) and at 30° (wide tomosynthesis projection). The noise power spectra (NPS) for both the scatter signal alone and the total signal (primary + scatter) for all images were obtained and the combined results of all patients analyzed. The total NPS was fit to the expected power-law relationship NPS(f) = k/f β and the results were compared with those previously published on the power spectrum characteristics of mammographic texture. The scatter signal alone was analyzed qualitatively and a power-law fit was also performed.nnnRESULTSnThe mammography and tomosynthesis projections of three patient breasts were too small to analyze, so a total of 16 patient breasts were analyzed. The values of β for the total signal of the 0° projections agreed well with previously published results. As expected, the scatter power spectrum reflected a fast drop-off with increasing spatial frequency, with a reduction of four orders of magnitude by 0.1 lp/mm. The β values for the scatter signal were 6.14 and 6.39 for the 0° and 30° projections, respectively.nnnCONCLUSIONSnAlthough the low-frequency characteristics of scatter in mammography and breast tomosynthesis were known, a quantitative analysis of the frequency domain characteristics of this signal was needed in order to optimize previously proposed software-based x-ray scatter reduction algorithms for these imaging modalities.

In-line phase-contrast breast tomosynthesis: a phantom feasibility study at a synchrotron radiation facility.

The major objective is to adopt, apply and test developed in-house algorithms for volumetric breast reconstructions from projection images, obtained in in-line phase-contrast mode. Four angular sets, each consisting of 17 projection images obtained from four physical phantoms, were acquired at beamline ID17, European Synchroton Radiation Facility, Grenoble, France. The tomosynthesis arc wasu2009u2009±32°. The physical phantoms differed in complexity of texture and introduced features of interest. Three of the used phantoms were in-house developed, and made of epoxy resin, polymethyl-methacrylate and paraffin wax, while the fourth phantom was the CIRS BR3D. The projection images had a pixel size of 47 µmu2009u2009×u2009u200947 µm. Tomosynthesis images were reconstructed with standard shift-and-add (SAA) and filtered backprojection (FBP) algorithms. It was found that the edge enhancement observed in planar x-ray images is preserved in tomosynthesis images from both phantoms with homogeneous and highly heterogeneous backgrounds. In case of BR3D, it was found that features not visible in the planar case were well outlined in the tomosynthesis slices. In addition, the edge enhancement index calculated for features of interest was found to be much higher in tomosynthesis images reconstructed with FBP than in planar images and tomosynthesis images reconstructed with SAA. The comparison between images reconstructed by the two reconstruction algorithms shows an advantage for the FBP method in terms of better edge enhancement. Phase-contrast breast tomosynthesis realized in in-line mode benefits the detection of suspicious areas in mammography images by adding the edge enhancement effect to the reconstructed slices.

Evaluation of the effect of silicone breast inserts on X-ray mammography and breast tomosynthesis images: A Monte Carlo simulation study.

PURPOSEnBreast augmentation is one of the most popular cosmetic surgeries worldwide. The aim of this study is to investigate the effect of breast implant insertion on the detectability and visibility of lesions on mammography and breast tomosynthesis (BT) images.nnnMATERIALS AND METHODSnThree software phantoms, composed of a homogeneous background with embedded silicone gel structures, and two types of breast abnormalities, microcalcifications (µCs) and masses, were generated. Two X-ray breast imaging modalities were simulated: mammography and BT with six incident monochromatic X-ray beams with energies in the interval between 20 and 30u2009keV. Projection images were generated using an in-house developed Monte Carlo simulator. The detectability of mammographic findings adjacent to the implant material and the influence of the incident beam energy and implant thickness on the feature detection were studied.nnnRESULTSnIt was found that implants thicker than 26u2009mm for the case of mammography and 14u2009mm for the case of BT obscured the visibility of underlying structures. Although BT demonstrated a lack of contrast, this modality was able to visualize µCs under considerable depths of implant. Increasing the incident beam energy led to better visualization of small µCs, while in the case of breast masses, their detectability was limited.nnnCONCLUSIONSnSilicone gel implants introduce a limitation in the image quality of mammograms resulting in low detectability of features. In addition, silicone gel implants obscure partially or totally parts of the image, depending on the size and the thickness of the implant as well the energy of the X-rays used.

Evaluation of a breast software model for 2D and 3D X-ray imaging studies of the breast

INTRODUCTIONnIn X-ray imaging, test objects reproducing breast anatomy characteristics are realized to optimize issues such as image processing or reconstruction, lesion detection performance, image quality and radiation induced detriment. Recently, a physical phantom with a structured background has been introduced for both 2D mammography and breast tomosynthesis. A software version of this phantom and a few related versions are now available and a comparison between these 3D software phantoms and the physical phantom will be presented.nnnMETHODSnThe software breast phantom simulates a semi-cylindrical container filled with spherical beads of different diameters. Four computational breast phantoms were generated with a dedicated software application and for two of these, physical phantoms are also available and they are used for the side by side comparison. Planar projections in mammography and tomosynthesis were simulated under identical incident air kerma conditions. Tomosynthesis slices were reconstructed with an in-house developed reconstruction software. In addition to a visual comparison, parameters like fractal dimension, power law exponent β and second order statistics (skewness, kurtosis) of planar projections and tomosynthesis reconstructed images were compared.nnnRESULTSnVisually, an excellent agreement between simulated and real planar and tomosynthesis images is observed. The comparison shows also an overall very good agreement between parameters evaluated from simulated and experimental images.nnnCONCLUSIONnThe computational breast phantoms showed a close match with their physical versions. The detailed mathematical analysis of the images confirms the agreement between real and simulated 2D mammography and tomosynthesis images. The software phantom is ready for optimization purpose and extrapolation of the phantom to other breast imaging techniques.

Computer Aided Preoperative Evaluation of the Residual Liver Volume Using Computed Tomography Images

Major hepatectomy causes a risk of postoperative liver dysfunction, failure, and infections like surgical site infection. Preoperative assessment of the liver volume and function of the remnant liver is a mandatory prerequisite before performing such surgery. The aim of this work is to develop and test a software application for evaluation of the residual function of the liver prior to the intervention of the surgeons. For this purpose, a technique for evaluation of liver volume from computed tomography (CT) images has been developed. Furthermore, the methodology algorithms were implemented and incorporated within a software tool with three basic functionalities: volume determination based on segmentation of liver from CT images, virtual tumour resection and estimation of the residual liver function and 3D visualisation. Forty-one sets of abdominal CT images consisting of different number of tomographic slice images were used to test and evaluate the proposed approach. Volumes that were obtained after manual tracing by two surgeon experts showed a relative difference of 3.5xa0%. The suggested methodology was encapsulated within an application with user-friendly interface that allows surgeons interactively to perform virtual tumour resection, to evaluate the relative residual liver and render the final result. Thereby, it is a tool in the surgeons’ hands that significantly facilitates their duties, saves time, and allows them to objectively evaluate the situation and take the right decisions. At the same time, the tool appears to be appropriate educational instrument for virtual training of young surgeon specialists.

A software platform for phase contrast x-ray breast imaging research

PURPOSEnTo present and validate a computer-based simulation platform dedicated for phase contrast x-ray breast imaging research.nnnMETHODSnThe software platform, developed at the Technical University of Varna on the basis of a previously validated x-ray imaging software simulator, comprises modules for object creation and for x-ray image formation. These modules were updated to take into account the refractive index for phase contrast imaging as well as implementation of the Fresnel-Kirchhoff diffraction theory of the propagating x-ray waves. Projection images are generated in an in-line acquisition geometry. To test and validate the platform, several phantoms differing in their complexity were constructed and imaged at 25 keV and 60 keV at the beamline ID17 of the European Synchrotron Radiation Facility. The software platform was used to design computational phantoms that mimic those used in the experimental study and to generate x-ray images in absorption and phase contrast modes.nnnRESULTSnThe visual and quantitative results of the validation process showed an overall good correlation between simulated and experimental images and show the potential of this platform for research in phase contrast x-ray imaging of the breast. The application of the platform is demonstrated in a feasibility study for phase contrast images of complex inhomogeneous and anthropomorphic breast phantoms, compared to x-ray images generated in absorption mode.nnnCONCLUSIONSnThe improved visibility of mammographic structures suggests further investigation and optimisation of phase contrast x-ray breast imaging, especially when abnormalities are present. The software platform can be exploited also for educational purposes.