N. Pallikarakis

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.

Influence of the bicarbonate pool and on the occurrence of 13CO2 in exhaled air

SummaryIn 13C02 breath tests, based on 13C:12C ratio measurements, the appearance of 13C in exhaled C02 was monitored after the administration of a 13C-la-belled compound. Independently of the substrate used, the existence of a bicarbonate pool into which the C02 produced enters before being exhaled, imposes a delay on the appearance of changes in the 13C:12C ratio. To estimate the nature and magnitude of this delay, we applied a two-compartment model to describe the kinetics of the body bicarbonate pool and we evaluated the 13C:12C ratio of C02 entering that pool from the measured 13C:12C ratio in the exhaled C02 after an oral intake of “naturally labelled” 13C-glucose. Our results demonstrated that discrepancies between total and exogenous glucose oxidation in relation to the peak occurrence time, as well as the absolute quantities, could be adequately explained by the interference of the bicarbonate stores.

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.

Evaluation of an improved algorithm for producing realistic 3D breast software phantoms: Application for mammography

PURPOSE This work presents an improved algorithm for the generation of 3D breast software phantoms and its evaluation for mammography. METHODS The improved methodology has evolved from a previously presented 3D noncompressed breast modeling method used for the creation of breast models of different size, shape, and composition. The breast phantom is composed of breast surface, duct system and terminal ductal lobular units, Coopers ligaments, lymphatic and blood vessel systems, pectoral muscle, skin, 3D mammographic background texture, and breast abnormalities. The key improvement is the development of a new algorithm for 3D mammographic texture generation. Simulated images of the enhanced 3D breast model without lesions were produced by simulating mammographic image acquisition and were evaluated subjectively and quantitatively. For evaluation purposes, a database with regions of interest taken from simulated and real mammograms was created. Four experienced radiologists participated in a visual subjective evaluation trial, as they judged the quality of the simulated mammograms, using the new algorithm compared to mammograms, obtained with the old modeling approach. In addition, extensive quantitative evaluation included power spectral analysis and calculation of fractal dimension, skewness, and kurtosis of simulated and real mammograms from the database. RESULTS The results from the subjective evaluation strongly suggest that the new methodology for mammographic breast texture creates improved breast models compared to the old approach. Calculated parameters on simulated images such as beta exponent deducted from the power law spectral analysis and fractal dimension are similar to those calculated on real mammograms. The results for the kurtosis and skewness are also in good coincidence with those calculated from clinical images. Comparison with similar calculations published in the literature showed good agreement in the majority of cases. CONCLUSIONS The improved methodology generated breast models with increased realism compared to the older model as shown in evaluations of simulated images by experienced radiologists. It is anticipated that the realism will be further improved using an advanced image simulator so that simulated images may be used in feasibility studies in mammography.

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 software data generator for radiographic imaging investigations

A software data generation tool, intended to be used in radiographic applications, has been developed. The application integrates a phantom design module and an imaging simulator. Phantoms can be described either as a set of geometrical objects or voxels, or as contours drawn on multiple tomographic slices. Radiographic projections of the phantoms are formed on the basis of a simulated irradiation process with selectable imaging parameters. A comparison between actual projection images from a physical and a simulated phantom shows good correspondence. The application was used for digital tomosynthesis (DTS) investigations and has proven to be a useful tool in the study of tomographic noise. Further development is expected to expand the use of the application to more areas of radiological imaging research.

A novel simulation algorithm for soft tissue compression.

This paper presents a novel general approach to simulation of soft tissue compression. A theoretical framework of the compression algorithm has been developed and implemented, based on the concept of a simple spring. The volume subjected to compression is divided into a number of “model elements”, each one consisting of 27 nodes. The nodes are connected with springs. The mechanical properties of the tissues are assumed to be linear and isotropic. The compressed volume remains constant due to the introduced concept of spring variable equilibrium lengths. Initial settings for compression simulation are introduced in order that the algorithm converges faster. The developed compression algorithm was used to model breast compression during a standard mammography examination. Specifically, the method was applied to a high-resolution three-dimensional software breast phantom, composed to have a medium glandularity and calcification abnormalities. The compression was set to 50%. Results showed that the abnormalities maintain their shape and dimensions during the compression, while the surrounding breast tissues undergo considerable deformation and displacement. A “decompression” algorithm was also applied to test the reversibility of the model.

Glucose Utilization During Exercise in Normal and Diabetic Subjects: The Role of Insulin

Due to selective isotopic effects occurring during photosynthesis, certain natural sugars are enriched in 13C Using such “naturally labeled 13C-glucose,” we studied glucose oxidation during exercise in seven normal volunteers and in six insulin-dependent diabetics after an overnight fast. In the diabetics, blood glucose was monitored the night before the test and adjusted to about 100 mg/dl by intravenous insulin infusion. The insulin infusion was withheld 15 min before exercise in four diabetics and maintained at 0.9 U/h for 2 h; then it was maintained at 0.6 U/h for 2 h in five diabetics. Three patients underwent both tests. All subjects exercised on a treadmill for 4 h at about 45% of their max. After 15 min adaptation, all received 100 g 13C-labeled glucose orally. Total glucose oxidation was derived from non-protein RQ and exogenous glucose oxidation evaluated as previously described. The diabetics had no residual B-cell function as indicated by negligible plasma C-peptide values and a lack of Cpeptide response to the oral glucose challenge. Total glucose oxidation averaged 230 ± 14 g/4 h in the normal subjects. It was similar (238 ± 19 g/4 h) in the diabetics receiving an intravenous insulin infusion, but decreased to 176 ± 14 g/4 h when no insulin was infused. Exogenous glucose oxidation was 92 ± 3 g/4 h and 84 ± 8 g/4 h (not statistically different) in the controls and in the insulin-infused diabetics, respectively. It was 43 ± 11 g/4 h in the diabetics exercising without being infused with insulin. We conclude that (1) in well-insulinized diabetic patients, prolonged muscular exercise can be performed under metabolic conditions which are basically similar to those of normal subjects; (2) during prolonged exercise, well-insulinized diabetic patients are able to oxidize up to 85–90% of a 100-g exogenous glucose load given orally and oral glucose can thus be ingested during prolonged exercise in well-controlled juvenile insulintreated diabetics; (3) even in the absence of insulin administration during exercise, juvenile diabetics who start exercising when blood glucose is near normal are able to perform a 4-h exercise at 45–50% of their max. Under these conditions, however, they are unable to utilize more than 40–45% of a 100-g glucose load given orally. They rely more upon lipid stores than the normal subjects or the well-insulinized diabetics.