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. Cooper's 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.

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.

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.

The quality of the reconstructed images in Digital Tomosynthesis is often limited by the presence of artifacts due to blur from planes other than the fulcrum plane. A technique has been developed for the separation and subsequent removal of unrelated structures from the reconstructed plane. The method involves the reconstruction of the blur originating in user-selected "noisy" planes as it appears on the plane of interest. This is achieved by projecting the reconstructed images of the selected plane on the image formation plane for all viewing angles, and subsequently, synthesizing its blurred image on the plane of interest. There are no restrictions as to the identity of the planes to be removed. Reproduction of the noise is performed using the tomosynthesis algorithm itself, thus the technique can be modified to suit any reconstruction algorithm. The proposed technique was implemented on the Multiple Projection Algorithm and was experimentally evaluated using a radiotherapy simulator unit.