Z. Kolitsi

A three-dimensional breast software phantom for mammography simulation

This paper presents a methodology for three-dimensional (3D) computer modelling of the breast, using a combination of 3D geometrical primitives and voxel matrices that can be further subjected to simulated x-ray imaging, to produce synthetic mammograms. The breast phantom is a composite model of the breast and includes the breast surface, the duct system and terminal ductal lobular units. Coopers ligaments, the pectoral muscle, the 3D mammographic background and breast abnormalities. A second analytical x-ray matter interaction modelling module is used to generate synthetic images from monoenergetic fan beams. Mammographic images of various synthesized breast models differing in size, shape and composition were produced. A preliminary qualitative assessment performed by three radiologists and a quantitative evaluation study using fractal and grey-level histogram analysis were conducted. A comparative study of extracted features with published data has also been performed. The evaluation results indicated good correlation of characteristics between synthetic and actual radiographs. Applications foreseen are not only in the area of breast imaging experimentation but also in education and training.

A multiple projection method for digital tomosynthesis

A new method of optimized efficiency for the retrospective reconstruction of tomograms is presented. The method has been developed for use with isocentric fluoroscopic units and is capable of performing digital tomosynthesis of anatomical planes of user selected orientation and distance from the isocenter. Optimization of efficiency has been achieved by segmenting the reconstruction process into discrete transformations that are specific to groups of pixels, rather than performing pixel by pixel operations. These involve a number of projections of the acquired image matrices as well as parallel translations and summing. Application of this method has resulted in a significant reduction of computing time. The proposed algorithm has been experimentally tested on a radiotherapy simulator unit with the use of a phantom and the obtained results are reported and discussed.

Dual-energy mammography: simulation studies

This paper presents a mammography simulator and demonstrates its applicability in feasibility studies in dual-energy (DE) subtraction mammography. This mammography simulator is an evolution of a previously presented x-ray imaging simulation system, which has been extended with new functionalities that are specific for DE simulations. The new features include incident exposure and dose calculations, the implementation of a DE subtraction algorithm as well as amendments to the detector and source modelling. The system was then verified by simulating experiments and comparing their results against published data. The simulator was used to carry out a feasibility study of the applicability of DE techniques in mammography, and more precisely to examine whether this modality could result in better visualization and detection of microcalcifications. Investigations were carried out using a 3D breast software phantom of average thickness, monoenergetic and polyenergetic beam spectra and various detector configurations. Dual-shot techniques were simulated. Results showed the advantage of using monoenergetic in comparison with polyenergetic beams. Optimization studies with monochromatic sources were carried out to obtain the optimal low and high incident energies, based on the assessment of the figure of merit of the simulated microcalcifications in the subtracted images. The results of the simulation study with the optimal energies demonstrated that the use of the DE technique can improve visualization and increase detectability, allowing identification of microcalcifications of sizes as small as 200 microm. The quantitative results are also verified by means of a visual inspection of the synthetic images.

An integrated research tool for X-ray imaging simulation

This paper presents a software simulation package of the entire X-ray projection radiography process including beam generation, absorber structure and composition, irradiation set up, radiation transport through the absorbing medium, image formation and dose calculation. Phantoms are created as composite objects from geometrical or voxelized primitives and can be subjected to simulated irradiation process. The acquired projection images represent the two-dimensional spatial distribution of the energy absorbed in the detector and are formed at any geometry, taking into account energy spectrum, beam geometry and detector response. This software tool is the evolution of a previously presented system, with new functionalities, user interface and an expanded range of applications. This has been achieved mainly by the use of combinatorial geometry for phantom design and the implementation of a Monte Carlo code for the simulation of the radiation interaction at the absorber and the detector.

Image quality in extended arc filtered digital tomosynthesis.

PURPOSE To study image quality in filtered digital tomosynthesis (FDTS) tomograms as a function of their reconstruction arc, using isocentrically acquired, fluoroscopic projection data. MATERIAL AND METHODS Both digital tomosynthesis (DTS) and cone beam CT (CBCT) reconstruction algorithms are based on backprojection and use cone beam projection data as input. Under limited angle conditions, CBCT is reduced to FDTS, where only a subset of projection data are used for reconstruction. The effect of the reconstruction arc on the spatial resolution, slice thickness, contrast sensitivity, shape distortion and artifacts, was also experimentally studied. The investigation was performed using both simulated and actual fluoroscopic images. RESULTS AND CONCLUSION Image quality in terms of spatial resolution, slice thickness, shape distortion and artifacts, improved with increasing reconstruction arc and was optimized at 180 degrees, while contrast continued to improve as the arc was increased to 360 degrees. However, DTS was determined to be the technique of choice when reconstruction arcs of less than 40 degrees were used. Consequently, FDTS may be successfully implemented in applications involving extended arc reconstructions, in the range between 40 degrees delimiting the DTS domain and 360 degrees corresponding to CBCT.

A wavelet-based method for removal of out-of-plane structures in digital tomosynthesis

Reconstructed images in digital tomosynthesis (DTS) are affected by artifacts due to blur from planes other than the fulcrum plane. A wavelet-based method has been developed for the discrimination and subsequent removal of unrelated structures from the reconstructed plane. The approach exploits both the specific pattern of noise in DTS and the spatial locality of the wavelet transformation. The technique was implemented on a DTS clinical protoype system. Experimental evaluation on angiographic types of images demonstrated excellent noise differentiation and elimination. The method is therefore particularly useful for certain medical imaging applications such as vascular DTS imaging.

Three-dimensional localisation based on projectional and tomographic image correlation: an application for digital tomosynthesis

Accurate three-dimensional tumor localisation in Radiotherapy, is critical to the treatment outcome, particularly when high dose gradients are present. A number of techniques have been proposed for the localisation of anatomical structures or markers. The present study proposes an approach to a concurrent maximisation of localisation accuracy and efficiency by correlation of tomographic and projectional images. The method introduces an element of direct verification and interactive optimisation of the process. Tomographic images are used for the identification of a point of interest. Its position is computed within the treatment co-ordinate system and verification of this position is achieved by obtaining the beams eye view of the identified point on two projection radiographs. The key element of the approach is that all images used should be part of one single image data set. The implementation of this localisation method, as part of the functionality of a Digital Tomosynthesis prototype, has provided an integrated facility for localisation, of optimised accuracy and precision, while easy and efficient to use. The considerations are general and apply in principle to any imaging system that can augment tomographic images with projections.

Experimental validation of a radiographic simulation code using breast phantom for X-ray imaging

Computer models and simulations of X-ray imaging systems are becoming a very precious tool during the development and evaluation of new X-ray imaging techniques. To provide, however, a faithful simulation of a system, all components must be accurately modelled and tested, followed by verification through experimental measurements. This paper presents a validation study of the XRayImagingSimulator, an in-house developed X-ray imaging simulator, which is extensively used as a basic tool in carrying out complex breast imaging simulations. The approach followed compares results obtained via an experimental setup for breast phantom (CIRS 011A) imaging, using synchrotron radiation (SYRMEP beamline at ELETTRA), with those from its simulated setup under the same conditions. The study demonstrated a very good agreement between experimental and simulated images compared both in terms of subjective and objective criteria. The combination of the XRayImagingSimulator with our BreastSimulator provides a powerful tool for in silico testing of new X-ray breast imaging approaches.

Evaluation of digital breast tomosynthesis reconstruction algorithms using synchrotron radiation in standard geometry

PURPOSE In this article, the image quality of reconstructed volumes by four algorithms for digital tomosynthesis, applied in the case of breast, is investigated using synchrotron radiation. METHODS An angular data set of 21 images of a complex phantom with heterogeneous tissue-mimicking background was obtained using the SYRMEP beamline at ELETTRA Synchrotron Light Laboratory, Trieste, Italy. The irradiated part was reconstructed using the multiple projection algorithm (MPA) and the filtered backprojection with ramp followed by hamming windows (FBR-RH) and filtered backprojection with ramp (FBP-R). Additionally, an algorithm for reducing the noise in reconstructed planes based on noise mask subtraction from the planes of the originally reconstructed volume using MPA (MPA-NM) has been further developed. The reconstruction techniques were evaluated in terms of calculations and comparison of the contrast-to-noise ratio (CNR) and artifact spread function. RESULTS It was found that the MPA-NM resulted in higher CNR, comparable with the CNR of FBP-RH for high contrast details. Low contrast objects are well visualized and characterized by high CNR using the simple MPA and the MPA-NM. In addition, the image quality of the reconstructed features in terms of CNR and visual appearance as a function of the initial number of projection images and the reconstruction arc was carried out. Slices reconstructed with more input projection images result in less reconstruction artifacts and higher detail CNR, while those reconstructed from projection images acquired in reduced angular range causes pronounced streak artifacts. CONCLUSIONS Of the reconstruction algorithms implemented, the MPA-NM and MPA are a good choice for detecting low contrast objects, while the FBP-RH, FBP-R, and MPA-NM provide high CNR and well outlined edges in case of microcalcifications.

Estimation of the heart respiratory motion with applications for cone beam computed tomography imaging: a simulation study

Computed tomography (CT) reconstruction methods assume imaging of static objects; object movement during projection data acquisition causes tomogram artifacts. The continuously moving heart, therefore, represents a complicated imaging case. The associated problems due to the heart beating can be overcome either by using very short projection acquisition times, during which the heart may be considered static, or by ECG-gated acquisition. In the latter case, however, the acquisition of a large number of projections may not be completed in a single breath hold, thus heart displacement occurs as an additional problem. This problem has been addressed by applying heart motion models in various respiratory motion compensation algorithms. Our paper focuses on cone beam computed tomography (CBCT), performed in conjunction with isocentric, fluoroscopic equipment, and continuous ECG and respiratory monitoring. Such equipment is used primarily for in-theater three-dimensional (3-D) imaging and benefits particularly from the recent developments in flat panel detector technologies. The objectives of this paper are: (i) to develop a model for the motion of the heart due to respiration during the respiratory cycle; (ii) to apply this model to the tomographic reconstruction algorithm, in order to account for heart movement due to respiration in the reconstruction; and (iii) to initially evaluate this method by means of simulation studies. Based on simulation studies, we were able to demonstrate that heart displacement due to respiration can be estimated from the same projection data, required for a CBCT reconstruction. Our paper includes semiautomatic segmentation of the heart on the X-ray projections and reconstruction of a convex 3-D-heart object that performs the same motion as the heart during respiration, and use of this information into the CBCT reconstruction algorithm. The results reveal significant image quality improvements in cardiac image reconstruction.