A Novel Method to Model Image Creation Based on Mammographic Sensors Performance Parameters: A Theoretical Study.

BACKGROUND: Mammographic digital imaging is based on X-ray sensors with solid image quality characteristics. These primarily include (a) a response curve that yields high contrast and image latitude, (b) a frequency response given by the Modulation Transfer Function (MTF), which enables small detail imaging and (c) the Normalize Noise Power Spectrum (NNPS) that shows the extent of the noise effect on image clarity. METHODS: In this work, a methodological approach is introduced and described for creating digital phantom images based on the measured image quality properties of the sensor. For this purpose, a mathematical phantom, simulating breast tissue and lesions of blood, adipose, muscle, Ca and Ca(50%)-P(50%) was created by considering the corresponding X-ray attenuation coefficients. The simulated irradiation conditions of the phantom used four mammographic spectra assuming exponential attenuation. Published data regarding noise and blur of a commercial RadEye HR CMOS imaging sensor were used as input data for the resulting images. RESULTS: It was found that the Ca and Ca(50%)-P(50%) lesions were visible in all exposure conditions. In addition, the W/Rh spectrum at 28 kVp provided more detailed images than the corresponding Mo/Mo spectrum. CONCLUSIONS: The presented methodology can act complementarily to image quality measurements, leading to initial optimization of the X-ray exposure parameters per clinical condition.

Grain Size Distribution Analysis of Different Activator Doped Gd2O2S Powder Phosphors for Use in Medical Image Sensors.

The structural properties of phosphor materials, such as their grain size distribution (GSD), affect their overall optical emission performance. In the widely used gadolinium oxysulfide (Gd2O2S) host material, the type of activator is one significant parameter that also changes the GSD of the powder phosphor. For this reason, in this study, different phosphors samples of Gd2O2S:Tb, Gd2O2S:Eu, and Gd2O2S:Pr,Ce,F, were analyzed, their GSDs were experimentally determined using the scanning electron microscopy (SEM) technique, and thereafter, their optical emission profiles were investigated using the LIGHTAWE Monte Carlo simulation package. Two sets of GSDs were examined corresponding to approximately equal mean particle size, such as: (i) 1.232 µm, 1.769 µm and 1.784 µm, and (ii) 2.377 µm, 3.644 µm and 3.677 µm, for Tb, Eu and Pr,Ce,F, respectively. The results showed that light absorption was almost similar, for instance, 25.45% and 8.17% for both cases of Eu dopant utilizing a thin layer (100 µm), however, given a thicker layer (200 µm), the difference was more obvious, 22.82%. On the other hand, a high amount of light loss within the phosphor affects the laterally directed light quanta, which lead to sharper distributions and therefore to higher resolution properties of the samples.

A novel method for the optimization of positron emission tomography scanners imaging performance.

OBJECTIVE: The aim of the present study was to propose a comprehensive method for positron emission tomography (PET) scanners image quality assessment, by simulation of a thin layer chromatography (TLC) flood source with a previously validated Monte Carlo model. METHODS AND MATERIALS: We used the GATE Monte Carlo package (GEANT4 application for tomographic emission) and the reconstructed images were obtained using the software for tomographic image reconstruction (STIR), with cluster computing. The PET scanner used in this simulation study was the General Electric Discovery-ST (USA). The plane source that was used for the image quality assessment was a TLC plate, consisting of an aluminum (Al) foil, coated with a thin layer of silica and immersed in fluorodeoxyglucose (18F-FDG) bath solution (1 MBq). The influence of different scintillating crystals on PET scanner's image quality, in terms of the modulation transfer function (MTF), the normalized noise power spectrum (NNPS) and the detective quantum efficiency (DQE), were also investigated. Modulation transfer function was estimated from transverse slices of the plane source, whereas the NNPS from the corresponding coronal slices. Images were reconstructed by the commonly used 2D filtered back projection (FBP2D), the Kinahan and Rogers FPB3DRP and the maximum likelihood estimation (MLE)-OSMAPOSL algorithms. Images obtained using the OSMAPOSL algorithm were assessed by using 15 subsets and 3 iterations. RESULTS: The PET scanner configuration, equipped with LuAP crystals, exhibited the optimum MTF values in both 2D and 3D FBP image reconstruction, whereas the corresponding configuration with BGO crystals exhibited the optimum MTF values after the iterative algorithm. The scanner equipped with the BGO crystals was also found to exhibit overall the lowest noise levels and the highest DQE values after algorithms. These finding indicate that the GE Discovery ST PET scanner exhibits the optimum image quality parameters, in terms of MTF, NNPS and DQE, with BGO scintillating crystals. CONCLUSION: Our new method showed that the imaging performance of PET scanners can be fully characterized and further improved by investigation of the imaging chain components through Monte Carlo methods. To this aim, a TLC based plane source was used during the simulation, in order to assess the impact of the scintillating crystal material on PET image quality, with the application of a previously validated Monte Carlo model. The aforementioned plane source can be also useful for the further development of PET and SPET scanners through GATE simulations, for clinical applications.

A robust method, based on a novel source, for performance and diagnostic capabilities assessment of the positron emission tomography system.

The aim of our work was to provide a robust method for evaluating imaging performance of positron emission tomography (PET) systems and particularly to estimate the modulation transfer function (MTF) using the line spread function (LSF) method. A novel plane source was prepared using thin layer chromatography (TLC) of a fluorine-18-fluorodeoxyglucose ((18)F-FDG) solution. The source was placed within a phantom, and imaged using the whole body (WB) two dimensional (2D) and three dimensional (3D) standard imaging protocols in a GE Discovery ST hybrid PET/CT scanner. Modulation transfer function was evaluated by determining the LSF, for various reconstruction methods and filters. The proposed MTF measurement method was validated against the conventional method, based on point spread function (PSF). Higher MTF values were obtained with 3D scanning protocol and 3D iterative reconstruction algorithm. All MTF obtained using 3D reconstruction algorithms showed better preservation of higher frequencies than the 2D algorithms. They also exhibited better contrast and resolution. MTF derived from LSF were more precise compared with those obtained from PSF since their reproducibility was better in all cases, providing a mean standard deviation of 0.0043, in contrary to the PSF method which gave 0.0405. In conclusion, the proposed method is novel and easy to implement for characterization of the signal transfer properties and image quality of PET/computed tomography (CT) systems. It provides an easy way to evaluate the frequency response of each kernel available. The proposed method requires cheap and easily accessible materials, available to the medical physicist in the nuclear medicine department. Furthermore, it is robust to aliasing and since this method is based on the LSF, is more resilient to noise due to greater data averaging than conventional PSF-integration techniques.

Towards the optimization of nuclear medicine procedures for better spatial resolution, sensitivity, scan image quality and quantitation measurements by using a new Monte Carlo model featuring PET imaging.

The aim of this study was to develop a Monte-Carlo model that can be used for the optimization of positron emission tomography (PET) procedures and image quality metrics. This model was developed using the Monte Carlo package of Geant4 application for tomographic emission (GATE) and the software for tomographic image reconstruction (STIR) with cluster computing to obtain reconstructed images. The PET scanner used in this study was the General Electric Discovery-ST (US). The GATE model was validated by comparing results obtained in accordance with the National Electrical Manufacturers Association NEMA-NU-2-2001 protocol [Mawlawi et al (2004) and Bettinardi et al (2004)]. All images were reconstructed with the commonly used 2D filtered back projection and the 3D reprojection algorithms. We found that the simulated spatial resolution in terms of full width at half maximum (FWHM) agreed within less than 3.29% in 2D and less than 2.51% in 3D with published data of others, respectively. The 2D values for the sensitivity, scatter fraction and count-rate were found to agree within less than 0.46%, 4.59% and 7.86%, respectively with these published values. Accordingly, our study showed that the corresponding 3D values were found to agree to less than 1.62%, 2.85% and 9.13%, respectively with Mawlawi et al (2004) published values. Sensitivity, which was also estimated without the presence of attenuation material by simulating an ideal source, showed differences between the extrapolated and the ideal source values (with and without attenuation) ranging in 2D from 0.04% to 0.82% (radial location R=0cm) and 0.52% to 0.67% in 3D mode (radial locations R=10cm). The simulated noise equivalent count rate was found to be 94.31kcps in 2D and 66.9kcps in 3D at 70 and 15kBq/mL respectively, compared to 94.08kcps in 2D and 70.88kcps in 3D at 54.6kBq/mL and 14kBq/mL respectively, from the published by others values. The simulated image quality was found in excellent agreement with these published values. In conclusion, our study showed that our Monte Carlo model can be used to assess, optimize, simplify and reduce the simulation time for the quality control procedure of PET scanners. By using this model, sensitivity can be obtained in a more simplified procedure. Reconstructed images by STIR can be also used to obtain radiopharmaceutical distribution of images and direct dose maps, quite useful to nuclear medicine practitioners.

A phantom based evaluation of the clinical imaging performance of electronic portal imaging devices.

Purpose: In this study an evaluation of the imaging performance of an electronic portal imaging device (EPID) is presented. The evaluation performed employing the QC-3V image quality phantom. Methods: An EPID system of a 6 MV LINAC, was used to obtain images of a QC-3V EPID phantom. The X-ray source to phantom distance was 100 cm and the field size was 15x15 cm2. The irradiation conditions comprised Dose Rates (DR) of 200, 400 and 600 for a 2 MU-100 MU range. The Contrast Transfer Function (CTF), the Noise Power Spectrum (NPS), the Normalized Noise Power Spectrum (NNPS) and the Contrast-to-Noise Ratio (CNR) were studied. In addition, an alternative factor showing a frequency related output signal-to-noise ratio (SNR), the Signal-to-Noise-Frequency Response (SNFR), has been introduced. SNFR is a comprehensive quality index, easily determined in clinical environment. Results: The CTF curves were found comparable to each other. The lowest values were measured at 2 MU and 200 MU/min. Concerning the NPS and NNPS graphs it was found that the values decrease up to approximately 0.3 lp/mm and demonstrate a white noise shape afterwards. SNFR values were found reducing with spatial frequency. Highest CNR were found between the region 7 and 11 of the phantom. Conclusions: The influence of MU and DR on EPID performance were investigated. Image quality was assessed using the QC-3V phantom. The presented results can lead to image quality amelioration and act supportively to current image quality control routine protocols.

MTF of columnar phosphors with a homogenous part: an analytical approach.

A method for the theoretical estimation of the MTF of columnar phosphors with a homogeneous part at the end used in X-ray imaging has been developed. This method considers the light transport inside the scintillator through an analytical modelling, the optical photon beams distribution on the scintillator-optical sensor interface, and uses the definition of the PSF and a Gauss fitted LSF to estimate the MTF of an indirect detector. This method was applied to a columnar CsI:Tl scintillator and validated against experimental results found in literature, and a good agreement was observed. It was found that, by increasing the pixel size of the optical detector and the thickness of the scintillator, the MTF decreased as expected. This method may be used in evaluating the performance of the columnar phosphors used in medical imaging, given their physical and geometrical characteristics.Graphical abstract (a) Side view of a part of the scintillator where five crystal columns with homogeneous ends attached to an optical sensor is shown. (b) Propagation of two random optical photon beams emitted from point K with different angles of emission is shown. All the symbols are explained analytically in the text.

Monte Carlo study on the imaging performance of powder Lu2SiO5:Ce phosphor screens under x-ray excitation: comparison with Gd2O2S:Tb screens.

Lu2SiO5: Ce (LSO) scintillator is a relatively new luminescent material which has been successfully applied in positron emission tomography systems. Since it has been recently commercially available in powder form, it could be of value to investigate its performance for use in x-ray projection imaging as both physical and scintillating properties indicate a promising material for such applications. In the present study, a custom and validated Monte Carlo simulation code was used in order to examine the performance of LSO, under diagnostic radiology (mammography and general radiography) conditions. The Monte Carlo code was based on a model using Mie scattering theory for the description of light attenuation. Imaging characteristics, related to image brightness, spatial resolution and noise of LSO screens were predicted using only physical parameters of the phosphor. The overall performance of LSO powder phosphor screens was investigated in terms of the: (i) quantum detection efficiency (ii) emitted K-characteristic radiation (iii) luminescence efficiency (iv) modulation transfer function (v) Swank factor and (vi) zero-frequency detective quantum efficiency [DQE(0)]. Results were compared to the traditional rare-earth Gd2O2S:Tb (GOS) phosphor material. The relative luminescence efficiency of LSO phosphor was found inferior to that of GOS. This is due to the lower intrinsic conversion efficiency of LSO (0.08 instead of 0.15 of GOS) and the relatively high light extinction coefficient mext of this phosphor (0.239 mircom(-1) instead of 0.218 /microm(-1) for GOS). However, the property of increased light extinction combined with the rather sharp angular distribution of scattered light photons (anisotropy factor g=0.624 for LSO instead of 0.494 for GOS) reduce lateral light spreading and improve spatial resolution. In addition, LSO screens were found to exhibit better x-ray absorption as well as higher signal to noise transfer properties in the energy range from 18 keV up to 50.2 keV (e.g. DQE(0)=0.62 at 18 keV and for 34 mg/cm2, instead of 0.58 for GOS). The results indicate that certain optical properties of LSO (optical extinction coefficient, scattering anisotropy factor) combined with the relatively high x-ray coefficients, make this material a promising phosphor which, under appropriate conditions, could be considered for use in x-ray projection imaging detectors.

An analytical approach to the light transport in columnar phosphors. Detector Optical Gain, angular distribution and the CsI:Tl paradigm.

PURPOSE: An analytical model has been developed for the light propagation in columnar phosphors, based on the optical photon propagation physical and geometrical principles. METHODS: This model accounts for the multiple reflections on the sides of the crystal column, as well as for the infinite forward and backward reflections of the propagated optical photon beams created in the crystal bulk. Additionally it considers the lateral propagated optical photon beams after multiple refractions from the neighbor columns and the optical photon attenuation inside the scintillator. The model was used to predict the Detector Optical Gain (DOG), and the angular distribution, of the columnar CsI:Tl scintillators, used in medical imaging. RESULTS: The model was validated against CsI:Tl optical photon transmission published results and good agreement was observed. It was, also, found that the DOG is affected by the length of the columns, as well as the incident X-ray energy spectrum. The results of the angular distribution are in accordance with the theory that the longer crystal columns have more directional light distribution. CONCLUSIONS: The results of DOG are in accordance with the use of short crystal columns for lower energies (mammography) and the use of long crystal columns for higher energies (general radiology). Angular distribution was found more directive for long crystal columns.

Detective quantum efficiency (DQE) in PET scanners: A simulation study.

The aim of the present study is to introduce the detective quantum efficiency (DQE) for the image quality assessment of positron emission tomography (PET) scanners. For this purpose, a thin layer chromatography (TLC) plane source was simulated using a previously validated, scanner and source geometry, Monte Carlo (MC) model. The model was developed with the Geant4 application for tomographic emission (GATE) MC package and reconstructed images obtained with the software for tomographic image reconstruction (STIR), with cluster computing. The GE Discovery ST PET scanner was simulated by using a previously validated code. A plane source consisting of a TLC plate, was simulated by a layer of silica gel on aluminum (Al) foil substrate, immersed in 18F-FDG bath solution (1MBq). Image quality was assessed in terms of the modulation transfer function (MTF) and the normalized noise power spectrum (NNPS) in order to obtain the detective quantum efficiency (DQE). MTF curves were estimated from transverse reconstructed images of the plane source, whereas the NNPS data were estimated from the corresponding coronal images. Images were reconstructed by the maximum likelihood estimation ordered subsets maximum a posteriori one step late (MLE)-OS-MAP-OSL algorithm, by using various subsets 1-21) and iterations 1-20). MTF values were found to increase up to the 12th iteration whereas remain almost constant thereafter. However, the range of the increase in the MTF is limited as the number of subsets increases. The noise levels were found to increase with the corresponding increase of both the number of iterations and subsets. The maximum NNPS value (0.517mm2) was observed for the 420 MLEM-equivalent iterations reconstructed image at 0cycles/mm. Finally DQE values were found to increase for spatial frequencies up to 0.038cycles/mm and to decrease thereafter with the corresponding increase in both number of iterations and subsets. The maximum DQE value (0.48 at 0.038cycles/mm) was obtained for the 8 MLEM-equivalent iterations image. The simulated PET evaluation method based on the TLC plane source can be useful in the quality control and in the further development of PET and SPECT scanners though GATE simulations.

Modeling granular phosphor screens by Monte Carlo methods.

The intrinsic phosphor properties are of significant importance for the performance of phosphor screens used in medical imaging systems. In previous analytical-theoretical and Monte Carlo studies on granular phosphor materials, values of optical properties, and light interaction cross sections were found by fitting to experimental data. These values were then employed for the assessment of phosphor screen imaging performance. However, it was found that, depending on the experimental technique and fitting methodology, the optical parameters of a specific phosphor material varied within a wide range of values, i.e., variations of light scattering with respect to light absorption coefficients were often observed for the same phosphor material. In this study, x-ray and light transport within granular phosphor materials was studied by developing a computational model using Monte Carlo methods. The model was based on the intrinsic physical characteristics of the phosphor. Input values required to feed the model can be easily obtained from tabulated data. The complex refractive index was introduced and microscopic probabilities for light interactions were produced, using Mie scattering theory. Model validation was carried out by comparing model results on x-ray and light parameters (x-ray absorption, statistical fluctuations in the x-ray to light conversion process, number of emitted light photons, output light spatial distribution) with previous published experimental data on Gd2O2S: Tb phosphor material (Kodak Min-R screen). Results showed the dependence of the modulation transfer function (MTF) on phosphor grain size and material packing density. It was predicted that granular Gd2O2S: Tb screens of high packing density and small grain size may exhibit considerably better resolution and light emission properties than the conventional Gd2O2S: Tb screens, under similar conditions (x-ray incident energy, screen thickness).

Measuring scatter radiation in diagnostic X rays for radiation protection purposes.

During the last decades, radiation protection and dosimetry in medical X-ray imaging practice has been extensively studied. The purpose of this study was to measure secondary radiation in a conventional radiographic room, in terms of ambient dose rate equivalent H*(10) and its dependence on the radiographic exposure parameters such as X-ray tube voltage, tube current and distance. With some exceptions, the results indicated that the scattered radiation was uniform in the space around the water cylindrical phantom. The results also showed that the tube voltage and filtration affect the dose rate due to the scatter radiation. Finally, the scattered X-ray energy distribution was experimentally calculated.

Evolution of microbial populations during traditional Feta cheese manufacture and ripening.

In three different dairies (A, B and C) located in Peloponess region (Southern Greece), traditional Feta cheese trials took place February to March using mixtures of sheep's and goat's milk. Only small variations in the evolution of microbial groups were observed during the whole ripening period. The main groups, such as thermophilic cocci, mesophilic lactococci, thermophilic lactobacilli, nonstarter lactic acid bacteria (NSLAB), presumptive Leuconostoc, enterococci and micrococci, reached their highest levels during the first 16 days, and then declined approximately 1-2 log units until the end of ripening. The remaining groups investigated, comprising yeasts, coliforms and Escherichia coli, were highest at day 4. The yeasts remained constant, while coliforms and E. coli decreased sharply and were not detectable after 120 days of ripening. A number of 146 isolates (dairy A) taken from all stages of the manufacturing and ripening process were purified and studied. Lactobacillus plantarum (58/146) and isolates of related species Lactobacillus pentosus and Lactobacillus paraplantarum (16/146) were the most common microorganisms found during cheese ripening. Streptococcus thermophilus (23/146) and Lactobacillus delbrueckii subsp. bulgaricus (20/146) were detected in high levels up to 20 days, and then gradually reduced. Enterococcus faecium (29/146) was found in all manufacturing and ripening stages.

A new PET resolution measurement method through Monte-Carlo simulations.

OBJECTIVE: The aim of this study was to propose a novel method for image quality assessment in PET scanners through estimation of the modulation transfer function (MTF) of a plane source. The simulation was implemented using the previously validated Monte-Carlo model. A comparison of the proposed method with the more traditional technique, based on a line source, was also performed. MATERIALS AND METHODS: The Geant4 application for tomographic emission (GATE) Monte-Carlo package was used for model development, and reconstructed images were obtained using software for tomographic image reconstruction (STIR) with cluster computing. A novel plane source consisting of a radioactive ((18)F-fluorodeoxyglucose) thin-layer chromatography plate was simulated (total source activity: 44.4 MBq) to assess image quality through the MTF. All images were reconstructed with the three-dimensional filtered back projection (FBP3DRP) and ordered-subsets expectation maximization (OSEM) reprojection algorithms. RESULTS: The MTFs obtained using ordered-subsets expectation maximization show, in all cases, that higher frequencies are preserved compared with those obtained using the FBP3DRP. In addition, the plane source method is less prone to noise than the conventional line source method (SD=0.0031 and 0.0203, respectively). CONCLUSION: The thin-layer chromatography-based plane source presented requires materials commonly found in a clinical environment and could be used to assess image quality in nuclear medicine departments and to further develop PET and single-photon emission computed tomography scanners through Monte-Carlo simulations.

Figure of image quality and information capacity in digital mammography.

OBJECTIVES: In this work, a simple technique to assess the image quality characteristics of the postprocessed image is developed and an easy to use figure of image quality (FIQ) is introduced. This FIQ characterizes images in terms of resolution and noise. In addition information capacity, defined within the context of Shannon's information theory, was used as an overall image quality index. MATERIALS AND METHODS: A digital mammographic image was postprocessed with three digital filters. Resolution and noise were calculated via the Modulation Transfer Function (MTF), the coefficient of variation, and the figure of image quality. In addition, frequency dependent parameters such as the noise power spectrum (NPS) and noise equivalent quanta (NEQ) were estimated and used to assess information capacity. RESULTS: FIQs for the "raw image" data and the image processed with the "sharpen edges" filter were found 907.3 and 1906.1, correspondingly. The information capacity values were 60.86 × 10(3) and 78.96 × 10(3) bits/mm(2). CONCLUSION: It was found that, after the application of the postprocessing techniques (even commercial nondedicated software) on the raw digital mammograms, MTF, NPS, and NEQ are improved for medium to high spatial frequencies leading to resolving smaller structures in the final image.

Effect of high-pressure treatment at various temperatures on indigenous proteolytic enzymes and whey protein denaturation in bovine milk.

The objective of the present study was to determine the effect of high pressure (HP) processing (200, 450 and 650 MPa) at various temperatures (20, 40 and 55 degrees C) on the total plasmin plus plasminogen-derived activity (PL), plasminogen activator(s) (PA) and cathepsin D activities and on denaturation of major whey proteins in bovine milk. Data indicated that transfer of both PL and PA from the casein micelles to milk serum occurred at all pressures utilized at room temperature (20 degrees C). In addition to the transfer of PL and PA from micelles, there were reductions in activities of PL (16-18%) and PA (38-62%) for the pressures 450 and 650 MPa, at room temperature. There were synergistic negative effects between pressure and temperature on residual PL activity at 450 and 650 MPa and on residual PA activity only at 450 MPa. Cathepsin D activity in the acid whey from HP-treated milk was in general baroresistant at room temperature. The residual activity of cathepsin D decreased significantly at 650 MPa and 40 degrees C and at the pressures 450 and 650 MPa at 55 degrees C. Synergistic negative effects on the amount of native beta-lactoglobulin were observed at 450 and 650 MPa and on the amount of native alpha-lactalbumin at 650 MPa. There were significant correlations between enzymatic activities (PL, PA and cathepsin D) and the residual native beta-lactoglobulin and alpha-lactalbumin in bovine milk. In conclusion, HP significantly affected the activity of indigenous proteolytic enzymes and whey protein denaturation in bovine milk. Reduction in activity of indigenous enzymes (PL, PA and cathepsin D) and transfer of PL and PA from the casein to milk serum induced by HP is expected to have a profound effect on cheese yield, proteolysis during cheese ripening and quality of UHT milk during storage.

Purification and characterization of chymosin and pepsin from kid.

The objective of this work was to study the characteristics of the gastric aspartic proteinases chymosin and pepsin which are constituents of the kid rennet. The two enzymes were extracted from abomasal tissue of one kid from a local indigenous breed, separated from each other by DEAE-cellulose chromatography and then were purified by gel filtration and anion-exchange chromatography. The molecular weights of the purified kid chymosin and pepsin as determined by gel filtration were 36 kDa and 40 kDa respectively. The isoelectric point of kid chymosin was as multiple forms of 3-6 zones at pH 4.6-5.1, while that of kid pepsin was at pH < or =3.0. Kid pepsin contained 0.37 molecules phosphorous per molecule and was totally inhibited by 5 muM pepstatin A, being more sensitive than kid chymosin. Both enzymes were almost equally as proteolytic as calf chymosin on total casein at pH 5.6. Kid pepsin activity was more pH and temperature dependent than kid chymosin activity. In comparison with the calf chymosin temperature sensitivity, the order of increased sensitivity was: calf chymosin