Employing Multiple Low-Dose PET Images (at Different Dose Levels) as Prior Knowledge to Predict Standard-Dose PET Images.

The existing deep learning-based denoising methods predicting standard-dose PET images (S-PET) from the low-dose versions (L-PET) solely rely on a single-dose level of PET images as the input of deep learning network. In this work, we exploited the prior knowledge in the form of multiple low-dose levels of PET images to estimate the S-PET images. To this end, a high-resolution ResNet architecture was utilized to predict S-PET images from 6 to 4% L-PET images. For the 6% L-PET imaging, two models were developed; the first and second models were trained using a single input of 6% L-PET and three inputs of 6%, 4%, and 2% L-PET as input to predict S-PET images, respectively. Similarly, for 4% L-PET imaging, a model was trained using a single input of 4% low-dose data, and a three-channel model was developed getting 4%, 3%, and 2% L-PET images. The performance of the four models was evaluated using structural similarity index (SSI), peak signal-to-noise ratio (PSNR), and root mean square error (RMSE) within the entire head regions and malignant lesions. The 4% multi-input model led to improved SSI and PSNR and a significant decrease in RMSE by 22.22% and 25.42% within the entire head region and malignant lesions, respectively. Furthermore, the 4% multi-input network remarkably decreased the lesions' SUVmean bias and SUVmax bias by 64.58% and 37.12% comparing to single-input network. In addition, the 6% multi-input network decreased the RMSE within the entire head region, within the lesions, lesions' SUVmean bias, and SUVmax bias by 37.5%, 39.58%, 86.99%, and 45.60%, respectively. This study demonstrated the significant benefits of using prior knowledge in the form of multiple L-PET images to predict S-PET images.

Direct inference of Patlak parametric images in whole-body PET/CT imaging using convolutional neural networks.

PURPOSE: This study proposed and investigated the feasibility of estimating Patlak-derived influx rate constant (Ki) from standardized uptake value (SUV) and/or dynamic PET image series. METHODS: Whole-body 18F-FDG dynamic PET images of 19 subjects consisting of 13 frames or passes were employed for training a residual deep learning model with SUV and/or dynamic series as input and Ki-Patlak (slope) images as output. The training and evaluation were performed using a nine-fold cross-validation scheme. Owing to the availability of SUV images acquired 60 min post-injection (20 min total acquisition time), the data sets used for the training of the models were split into two groups: "With SUV" and "Without SUV." For "With SUV" group, the model was first trained using only SUV images and then the passes (starting from pass 13, the last pass, to pass 9) were added to the training of the model (one pass each time). For this group, 6 models were developed with input data consisting of SUV, SUV plus pass 13, SUV plus passes 13 and 12, SUV plus passes 13 to 11, SUV plus passes 13 to 10, and SUV plus passes 13 to 9. For the "Without SUV" group, the same trend was followed, but without using the SUV images (5 models were developed with input data of passes 13 to 9). For model performance evaluation, the mean absolute error (MAE), mean error (ME), mean relative absolute error (MRAE%), relative error (RE%), mean squared error (MSE), root mean squared error (RMSE), peak signal-to-noise ratio (PSNR), and structural similarity index (SSIM) were calculated between the predicted Ki-Patlak images by the two groups and the reference Ki-Patlak images generated through Patlak analysis using the whole acquired data sets. For specific evaluation of the method, regions of interest (ROIs) were drawn on representative organs, including the lung, liver, brain, and heart and around the identified malignant lesions. RESULTS: The MRAE%, RE%, PSNR, and SSIM indices across all patients were estimated as 7.45 ± 0.94%, 4.54 ± 2.93%, 46.89 ± 2.93, and 1.00 ± 6.7 × 10-7, respectively, for models predicted using SUV plus passes 13 to 9 as input. The predicted parameters using passes 13 to 11 as input exhibited almost similar results compared to the predicted models using SUV plus passes 13 to 9 as input. Yet, the bias was continuously reduced by adding passes until pass 11, after which the magnitude of error reduction was negligible. Hence, the predicted model with SUV plus passes 13 to 9 had the lowest quantification bias. Lesions invisible in one or both of SUV and Ki-Patlak images appeared similarly through visual inspection in the predicted images with tolerable bias. CONCLUSION: This study concluded the feasibility of direct deep learning-based approach to estimate Ki-Patlak parametric maps without requiring the input function and with a fewer number of passes. This would lead to shorter acquisition times for WB dynamic imaging with acceptable bias and comparable lesion detectability performance.

Dosimetry of gamma chamber blood irradiator using PAGAT gel dosimeter and Monte Carlo simulations.

Currently, the use of blood irradiation for inactivating pathogenic microbes in infected blood products and preventing graft-versus-host disease (GVHD) in immune suppressed patients is greater than ever before. In these systems, dose distribution and uniformity are two important concepts that should be checked. In this study, dosimetry of the gamma chamber blood irradiator model Gammacell 3000 Elan was performed by several dosimeter methods including thermoluminescence dosimeters (TLD), PAGAT gel dosimetry, and Monte Carlo simulations using MCNP4C code. The gel dosimeter was put inside a glass phantom and the TL dosimeters were placed on its surface, and the phantom was then irradiated for 5 min and 27 sec. The dose values at each point inside the vials were obtained from the magnetic resonance imaging of the phantom. For Monte Carlo simulations, all components of the irradiator were simulated and the dose values in a fine cubical lattice were calculated using tally F6. This study shows that PAGAT gel dosimetry results are in close agreement with the results of TL dosimetry, Monte Carlo simulations, and the results given by the vendor, and the percentage difference between the different methods is less than 4% at different points inside the phantom. According to the results obtained in this study, PAGAT gel dosimetry is a reliable method for dosimetry of the blood irradiator. The major advantage of this kind of dosimetry is that it is capable of 3D dose calculation.

Impact of rare earth element added filters on the X-ray beam spectra: a Monte Carlo approach.

The effectiveness of added filters including conventional and rare earth materials for dental radiography tasks was investigated using a simulation approach. Current study focuses on the combination of a range of various filters to investigate the reduction of radiation absorbed dose and improving the quality of a radiography image. To simulate the X-ray beam spectrum, a MCNP5 code was applied. Relative intensity, beam quality, and mean energy were investigated for a typical dental radiography machine. The impact of different rare-earth materials with different thicknesses and tube voltages on the X-ray spectrum was investigated. For Aluminum as a conventional filter, the modeled X-ray spectra and HVL values were in a good agreement with those reported by IPEM. The results showed that for a 70 kVp voltage, with an increase of the thickness and atomic number of a given added filters, an increase of HVL values were observed. However, with the increase of the attenuator thickness, X-ray beam intensity decreases. For mean energy, different results were observed. It was also found that rare earth made filters reduce high energy X-ray radiation due to k-edge absorption. This leads to an ideal beam for intra-oral radiography tasks. However, as a disadvantage of rare earth added filters, the reduction of the tube output levels should also be considered.

Design and implementation of continuous bed motion (CBM) in Xtrim preclinical PET scanner for whole-body Imaging: MC simulation and experimental measurements.

PURPOSE: Preclinical PET scanners often have limited axial field-of-view for whole-body (WB) scanning of the small-animal. Step-and-shoot(S&S) acquisition mode requires multiple bed positions (BPs) to cover the scan length. Alternatively, in Continuous Bed Motion(CBM) mode, data acquisition is performed while the bed is continuously moving. In this study, to reduce acquisition time and enhance image quality, the CBM acquisition protocol was optimized and implemented on the Xtrim-PET preclinical scanner for WB imaging. METHODS: The over-scan percentage(OS%) in CBM mode was optimized by Monte Carlo simulation. Bed movement speed was optimized considering ranges from 0.1 to 2.0 mm s-1, and absolute system sensitivities with the optimal OS% were calculated. The performance of the scanner in CBM mode was measured, and compared with S&S mode based on the NEMA-NU4 standard. RESULTS: The optimal trade-off between absolute sensitivity and uniformity of sensitivity profile was achieved at OS-50 %. In comparison to S&S mode with maximum ring differences (MRD) of 9 and 23, the calculated equivalent speeds in CBM(OS-50 %) mode were 0.3 and 0.14 mm s-1, respectively. In terms of data acquisition with equal sensitivity in both CBM(OS-50 %) and S&S(MRD-9) modes, the total scan time in CBM mode decreased by 25.9 %, 47.7 %, 54.7 %, and 58.2 % for scan lengths of 1 to 4 BPs, respectively. CONCLUSION: The CBM mode enhances WB PET scans for small-animals, offering rapid data acquisition, high system sensitivity, and uniform axial sensitivity, leading to improved image quality. Its efficiency and customizable scan length and bed speed make it a superior alternative.

Perturbation of TG-43 parameters of the brachytherapy sources under insufficient scattering materials.

In the recommendations of Task Group #43 from American Association of Physicists in Medicine (AAPM TG43), methods of brachytherapy source dosimetry are recommended, under full scattering conditions. However, in actual brachytherapy procedures, sources may not be surrounded by full scattering tissue in all directions. Clinical examples include high-dose-rate (HDR) brachytherapy of the breast or low-dose-rate (LDR) brachytherapy of ocular melanoma using eye plaque treatment with 125I and 103Pd. In this work, the impact of the missing tissue on the TG-43-recommended dosimetric parameters of different brachytherapy sources was investigated. The impact of missing tissue on the TG-43-recommended dosimetric parameters of 137Cs, 192Ir, and 103Pd brachytherapy sources was investigated using the MCNP5 Monte Carlo code. These evaluations were performed by placing the sources at different locations inside a 30 × 30 × 30 cm3 cubical water phantom and comparing the results with the values of the source located at the center of the phantom, which is in a full scattering condition. The differences between the thickness of the overlying tissues for different source positions and the thickness of the overlying tissue in full scattering condition is referred to as missing tissue. The results of these investigations indicate that values of the radial dose function and 2D anisotropy function vary as a function of the thickness of missing tissue, only in the direction of the missing tissue. These changes for radial dose function were up to 5%, 11%, and 8% for 137Cs, 192Ir, and 103Pd, respectively. No significant changes are observed for the values of the dose rate constants. In this project, we have demonstrated that the TG-43 dosimetric parameters may only change in the directions of the missing tissue. These results are more practical than the published data by different investigators in which a symmetric effect of the missing tissue on the dosimetric parameters of brachytherapy source are being considered, regardless of the implant geometry in real clinical cases.

Monte Carlo Dosimetric Study of Percutaneous Vertebroplasty and Brachytherapy for the Treatment of Spinal Metastases.

Background: Percutaneous vertebroplasty employs bone cement for injecting into the fractured vertebral body (VB) caused by spinal metastases. Radioactive bone cement and also brachytherapy seeds have been utilized to suppress the tumor growth in the VB. Objective: This study aims to investigate the dose distributions of low-energy brachytherapy seeds, and to compare them to those of radioactive bone cement, by Monte Carlo simulation. Material and Methods: In this simulation study, nine CT scan images were imported in Geant4. For the simulation of brachytherapy, I-125, Cs-131, or Pd-103 seeds were positioned in the VB, and for the simulation of vertebroplasty, the VB was filled by a radioactive cement loaded by P-32, Ho-166, Y-90, or Sm-153 radioisotopes. The dose-volume histograms of the VB, and the spinal cord (SC) were obtained after segmentation, considering that the reference dose is the minimum dose covered 95% of the VB. Results: The SC sparing was improved by using beta-emitting cement because of their steep gradient dose distribution. I-125 seeds and Y-90 radioisotope showed better VB coverage for brachytherapy and vertebroplasty techniques, respectively. Pd-103 seeds and P-32 radioisotope showed better SC sparing for brachytherapy and vertebroplasty, respectively. The minimum mean doses that covered 100% of the VB were 62.0%, 56.5%, and 45.0% for I-125, Cs-131, and Pd-103 seeds, and 28.3%, 28.6%, 32.9%, and 17.7%, for P-32, Ho-166, Y-90, and Sm-153 sources, respectively. Conclusion: I-125 and Cs-131 seeds may be useful for large tumors filling the entire VB, and also for the extended tumors invading multiple vertebrae. Beta-emitting bone cement is recommended for tumors located near the SC.

QCard-NM: Developing a semiautomatic segmentation method for quantitative analysis of the right ventricle in non-gated myocardial perfusion SPECT imaging.

BACKGROUND: Recent studies have shown that the right ventricular (RV) quantitative analysis in myocardial perfusion imaging (MPI) SPECT can be beneficial in the diagnosis of many cardiopulmonary diseases. This study proposes a new algorithm for right ventricular 3D segmentation and quantification. METHODS: The proposed Quantitative Cardiac analysis in Nuclear Medicine imaging (QCard-NM) algorithm provides RV myocardial surface estimation and creates myocardial contour using an iterative 3D model fitting method. The founded contour is then used for quantitative RV analysis. The proposed method was assessed using various patient datasets and digital phantoms. First, the physician's manually drawn contours were compared to the QCard-NM RV segmentation using the Dice similarity coefficient (DSC). Second, using repeated MPI scans, the QCard-NM's repeatability was evaluated and compared with the QPS (quantitative perfusion SPECT) algorithm. Third, the bias of the calculated RV cavity volume was analyzed using 31 digital phantoms using the QCard-NM and QPS algorithms. Fourth, the ability of QCard-NM analysis to diagnose coronary artery diseases was assessed in 60 patients referred for both MPI and coronary angiography. RESULTS: The average DSC value was 0.83 in the first dataset. In the second dataset, the coefficient of repeatability of the calculated RV volume between two repeated scans was 13.57 and 43.41 ml for the QCard-NM and QPS, respectively. In the phantom study, the mean absolute percentage errors for the calculated cavity volume were 22.6% and 42.2% for the QCard-NM and QPS, respectively. RV quantitative analysis using QCard-NM in detecting patients with severe left coronary system stenosis and/or three-vessel diseases achieved a fair performance with the area under the ROC curve of 0.77. CONCLUSION: A novel model-based iterative method for RV segmentation task in non-gated MPI SPECT is proposed. The precision, accuracy, and consistency of the proposed method are demonstrated by various validation techniques. We believe this preliminary study could lead to developing a framework for improving the diagnosis of cardiopulmonary diseases using RV quantitative analysis in MPI SPECT.

Characteristics of miniature electronic brachytherapy x-ray sources based on TG-43U1 formalism using Monte Carlo simulation techniques.

PURPOSE: The goal of this study is to determine a method for Monte Carlo (MC) characterization of the miniature electronic brachytherapy x-ray sources (MEBXS) and to set dosimetric parameters according to TG-43U1 formalism. TG-43U1 parameters were used to get optimal designs of MEBXS. Parameters that affect the dose distribution such as anode shapes, target thickness, target angles, and electron beam source characteristics were evaluated. Optimized MEBXS designs were obtained and used to determine radial dose functions and 2D anisotropy functions in the electron energy range of 25-80 keV. METHODS: Tungsten anode material was considered in two different geometries, hemispherical and conical-hemisphere. These configurations were analyzed by the 4C MC code with several different optimization techniques. The first optimization compared target thickness layers versus electron energy. These optimized thicknesses were compared with published results by Ihsan et al. [Nucl. Instrum. Methods Phys. Res. B 264, 371-377 (2007)]. The second optimization evaluated electron source characteristics by changing the cathode shapes and electron energies. Electron sources studied included; (1) point sources, (2) uniform cylinders, and (3) nonuniform cylindrical shell geometries. The third optimization was used to assess the apex angle of the conical-hemisphere target. The goal of these optimizations was to produce 2D-dose anisotropy functions closer to unity. An overall optimized MEBXS was developed from this analysis. The results obtained from this model were compared to known characteristics of HDR (125)I, LDR (103)Pd, and Xoft Axxent™ electronic brachytherapy source (XAEBS) [Med. Phys. 33, 4020-4032 (2006)]. RESULTS: The optimized anode thicknesses as a function of electron energy is fitted by the linear equation Y (µm) = 0.0459X (keV)-0.7342. The optimized electron source geometry is obtained for a disk-shaped parallel beam (uniform cylinder) with 0.9 mm radius. The TG-43 distribution is less sensitive to the shape of the conical-hemisphere anode than the hemispherical anode. However, the optimized apex angle of conical-hemisphere anode was determined to be 60°. For the hemispherical targets, calculated radial dose function values at a distance of 5 cm were 0.137, 0.191, 0.247, and 0.331 for 40, 50, 60, and 80 keV electrons, respectively. These values for the conical-hemisphere targets are 0.165, 0.239, 0.305, and 0.412, respectively. Calculated 2D anisotropy functions values for the hemispherical target shape were F(1 cm, 0°) = 1.438 and F(1 cm, 0°) = 1.465 for 30 and 80 keV electrons, respectively. The corresponding values for conical-hemisphere targets are 1.091 and 1.241, respectively. CONCLUSIONS: A method for the characterizations of MEBXS using TG-43U1 dosimetric data using the MC MCNP4C has been presented. The effects of target geometry, thicknesses, and electron source geometry have been investigated. The final choices of MEBXS design are conical-hemisphere target shapes having an apex angle of 60°. Tungsten material having an optimized thickness versus electron energy and a 0.9 mm radius of uniform cylinder as a cathode produces optimal electron source characteristics.

A new approach to overcome the inconsistency between SPECT and the anatomical map in maximum A-posterior expectation-maximization reconstruction algorithm.

Noise reduction while preserving spatial resolution is one of the most important challenges in the reconstructing of emission tomography images. One of the resolving methods is the Bowsher maximum a-posteriori expectation-maximization reconstruction (MAPEM) algorithm. This method considers a binary selection of the neighbors of each voxel based on the prior anatomical values to use in the regularization function. This method is particularly susceptible to imposing the wrong data into the reconstructed image due to the spatial or functional inconsistencies between the anatomical image and the actual activity distribution. Because of the poor spatial resolution of single-photon emission tomography (SPECT) images and the different nature of emission and anatomical imaging, there is not enough certainty of inconsistency with anatomical images. Therefore, we proposed a new weighted Bowsher method that can overcome this weakness while the image quality indexes, especially the spatial resolution, are almost preserved. In the proposed method, each of the neighbors of a specific voxel takes a constant weight depending on the order of its value and independent of its intensity quantity. The proposed method was evaluated using some different physical phantoms and a patient scan. The results show that the proposed method has superiority in the presence of inconsistency; moreover, the proposed method gives nearly similar results to the regular Bowsher MAPEM in case of consistency. In conclusion, we show that using a suitable constant weighting factor in Bowsher MAPEM, one can operatively reduce the image noise while preserving the image quality parameters where the emission tomography images are either consistent or inconsistent with the prior anatomical map.

Impact of the vaginal applicator and dummy pellets on the dosimetry parameters of Cs-137 brachytherapy source.

In this study, dose rate distribution around a spherical 137Cs pellet source, from a low-dose-rate (LDR) Selectron remote afterloading system used in gynecological brachytherapy, has been determined using experimental and Monte Carlo simulation techniques. Monte Carlo simulations were performed using MCNP4C code, for a single pellet source in water medium and Plexiglas, and measurements were performed in Plexiglas phantom material using LiF TLD chips. Absolute dose rate distribution and the dosimetric parameters, such as dose rate constant, radial dose functions, and anisotropy functions, were obtained for a single pellet source. In order to investigate the effect of the applicator and surrounding pellets on dosimetric parameters of the source, the simulations were repeated for six different arrangements with a single active source and five non-active pellets inside central metallic tubing of a vaginal cylindrical applicator. In commercial treatment planning systems (TPS), the attenuation effects of the applicator and inactive spacers on total dose are neglected. The results indicate that this effect could lead to overestimation of the calculated F(r,θ), by up to 7% along the longitudinal axis of the applicator, especially beyond the applicator tip. According to the results obtained in this study, in a real situation in treatment of patients using cylindrical vaginal applicator and using several active pellets, there will be a large discrepancy between the result of superposition and Monte Carlo simulations.

The effects of a transverse magnetic field on the dose enhancement of nanoparticles in a proton beam: a Monte Carlo simulation.

High-Z nanoparticles (NP) as radio-sensitization agents provide the feasibility of dose localization within the tumor in radiotherapy. Dose enhancement of NPs in the presence of a magnetic field (MF) could be challenged when magnetic resonance imaging (MRI) systems are used as an image-guided system. The MF can influence dose enhancement of NPs at their interfaces and surrounding medium and affect their dose deposition behavior. In the TOPAS Monte Carlo code, gold nanoparticle (GNP) and superparamagnetic iron oxide nanoparticle (SPION) were irradiated using 70 and 150 MeV proton beams, in presence of transverse MF strengths with 0, 1, 3, and 7 T. The changes in the liberated secondary electrons from NPs and their dose enhancement ratio (DER), magnetic dose enhancement ratio (MDER), and angular dose distribution in 10 nm shell thicknesses up to 500 nanometers from their centers were measured. The central plane of NPs was considered as a scorer. Its thickness was 2 nm and divided into 6-degree sectors with 10 nm radial length. The dose deposition in this voxelated scorer was calculated. The values of the deposited doses around NPs decrease rapidly while the DERs resulted from the secondary electrons are increased. MDERs are changed within [Formula: see text] and [Formula: see text] for 20 and 50 nm radius NPs, respectively. The variation in the angular dose distribution around a singular NP was not considerable when different MF strengths were applied. The dose values in the voxelated central plane show very similar results for the same NPs types in the different MF strengths. The typically used MF in the MRI systems would not considerably affect the energy deposition behavior of the secondary electrons produced in the interaction of proton beam with NPs, at least in the near vicinity of NPs. The DERs of NPs in a water medium resulted from emerged secondary electrons, experience a low degree of perturbation in the presence of an MF. The results of this study show that the NPs as dose enhancement agents can also be used in an MF without pronounced modification in their efficacy.

Assessment of spring water quality and associated health risks in a high-level natural radiation area, North Iran.

In this study, the spring water quality of a high natural background radiation region in North Iran was evaluated by measuring hydrochemical characteristics and concentration of potentially toxic elements (PTEs) and activity concentration of radioactive elements. The carcinogenicity potential from exposure to PTEs and radioactive elements was also investigated using probabilistic approach. The hydrochemical properties of water samples revealed that there were two different water types in the study area: (1) non-thermal Ca-HCO3 type and (2) thermal Na-Cl type. The concentrations of Al, Cd, Cr, Cu, Pb, Mo, Ni, Sb, Zn, and 40K were within the recommended water quality standard set by the World Health Organization (WHO). Elevated concentrations of As, Mn, Hg, 226Ra, and 232Th occur in Na-Cl water type while elevated concentrations of Fe occur in Ca-HCO3 water type. Also, health complications of dermal contact (via balneology or bathing) are within the safe limits. The major concern regarding the contaminated springs is the possibility of soil and groundwater contamination through uncontrolled runoff and spa effluents. Preventing the spread of toxic constituents in the study area via high-risk spring water requires periodic monitoring, and applying control measures where necessary.

Superior X-ray Radiation Shielding Effectiveness of Biocompatible Polyaniline Reinforced with Hybrid Graphene Oxide-Iron Tungsten Nitride Flakes.

X-ray radiation is a harmful carcinogenic electromagnetic source that can adversely affect the health of living species and deteriorate the DNA of cells, thus it's vital to protect vulnerable sources from them. To address this flaw, the conductive polymeric structure of polyaniline (PANi) was reinforced with diverse filler loadings (i.e., 25 wt % and 50 wt %) of hybrid graphene oxide-iron tungsten nitride (ITN) flakes toward attenuation of X-ray beams and inhabitation of microorganisms' growth. Primary characterizations confirmed the successful decoration of graphene oxide (GO) with interconnected and highly dense structure of iron tungsten nitride with a density of about 24.21 g.cm⁻3 and reinforcement of PANi with GO-ITN. Additionally, the outcome of evaluations showed the superior performance of developed shields, where a shield with 1.2 mm thickness containing 50 wt % GO-ITN showed 131.73 % increase in the electrical conductivity (compared with neat PANi) along with 78.07%, 57.12%, and 44.99% decrease in the amplitude of the total irradiated X-ray waves at 30, 40, and 60 kVp tube voltages, respectively, compared with control X-ray dosage. More importantly, the developed shields not only showed non-toxic nature and improved the viability of cells, but also completely removed the selected microorganisms at a concentration of 1000 µg.mL-1.

Quantification of contrast agent materials using a new image- domain multi material decomposition algorithm based on dual energy CT.

OBJECTIVE: Dual-Energy CT (DECT) is an imaging modality in which the objects are scanned by two different energy spectra. Using these two measurements, two type of materials can be separated and density image pairs can be generated as well. Decomposing more than two materials is necessary in both clinical and industrial CT applications. METHODS: In our MMD, barycentric coordinates were chosen using an innovative local clustering method. Local clustering increases precision in the barycentric coordinates assignment by decreasing search domain. Therefore the algorithm can be run in parallel. For optimizing coordinates selection, a fast bi-directional Hausdorff distance measurement is used. To deal with the significant obstacle of noise, we used Doubly Local Wiener Filter Directional Window (DLWFDW) algorithm. RESULTS: Briefly, the proposed algorithm separates blood and fat ROIs with errors of less than 2 and 9 % respectively on the clinical images. Also, the ability to decompose different materials with different concentrations is evaluated employing the phantom data. The highest accuracy obtained in separating different materials with different concentrations was 93 % (for calcium plaque) and 97.1 % (for iodine contrast agent) respectively. The obtained results discussed in detail in the following results section. CONCLUSION: In this study, we propose a new material decomposition algorithm. It improves the MMD work flow by employing tools which are easy to implement. Furthermore, in this study, an effort has been made to turn the MMD algorithm into a semi-automatic algorithm by employing clustering concept in material coordinate's assignment. The performance of the proposed method is comparable to existing methods from qualitative and quantitative aspects. ADVANCES IN KNOWLEDGE: All decomposition methods have their own specific problems. Image- domain decomposition also has barriers and problems, including the need for a predetermined table for the separation of different materials with specified coordinates. In the present study, it attempts to solve this problem by using clustering methods and relying on the intervals between different materials in the attenuation domain.

Fast Padé transform for increasing the signal to noise ratio of spectra provided by STEAM pulse sequence.

BACKGROUND/OBJECTIVE: There are two routine pulse-sequences for single voxel spectroscopy (SVS), point resolved spectroscopy (PRESS) and stimulated echo acquisition mode (STEAM). Although STEAM has several advantages in comparison to PRESS, signal/noise ratio (SNR) superiority of PRESS makes it the first choice for SVS. Application of fast Padé transform (FPT) instead of Fast Furrier transform (FFT) might increase the SNR of the signal produced by STEAM pulse-sequence and therefore allows the benefits of its advantages. We aimed to evaluate and compare the noise root mean square (RMS) and SNR provided by STEAM pulse-sequence using both FPT and FFT. MATERIALS AND METHOD: A gelatin-based phantom was constructed in a 19-cm acrylic cylinder. The phantom had two normal/tumoral parts. The SVS was performed using a 3T MRI scanner. STEAM pulse-sequence were used with the following parameters: TR = 2000 ms, TM = 10 ms, and three TEs of 20, 135 and 270 ms with two data-points of 1024 and 512 and voxel-size of 1 cm3. The raw data were extracted and processed using both FFT and FPT estimators to produce the spectrum. The noise RMS and SNR of Cho and Cr metabolites were assessed. RESULTS: According to the results, noise RMS of spectra provided by FPT were decreased between 3619.01-14252.94% in comparison to FFT (p< 0.00001). The SNR of Cr1 and Cho peaks of the spectra provided by FPT were increased more than 96.80 and 97.18, respectively (0.00006

A method for cranial target delineation in radiotherapy treatment planning aided by single-voxel magnetic resonance spectroscopy: evaluation using a custom-designed gel-based phantom and simulations.

OBJECTIVE: Magnetic resonance spectroscopy (MRS) has been useful in radiotherapy treatment planning (RTP) especially in tumor delineation. Routinely, 2D/3D MRSI data are used for this application. However, not all centers have access to 2D/3D MRSI. The objective of this study was to introduce a method of using single-voxel spectroscopy (SVS) data in target delineation and assess its reliability. METHODS: A gel-based phantom containing Creatine (Cr), N-acetyl-l-aspartic-acid (NAA), and Choline (Cho) was designed and built. The metabolite ratios simulate the normal and tumoral part of the brain. The jMRUI software (v. 6.0) was used to simulate a 1.5 T GE MRI scanner. The metabolite spectra provided by different time of echos (TE)s of the Point-RESolved Spectroscopy pulse-sequence (PRESS), different data-points, and post-processings were quantized by jMRUI. PseudoMRSI maps of Cho/Cr, NAA/Cr, and Cho + Cr/NAA were created. A conformity index (CI) was used to determine which metabolite-ratio isolines are more appropriate for tumor delineation. RESULTS: The simulation accuracy was verified. There were no differences > 4% between the measured and simulated spectra in peak regions. The pseudoMRSI map of Cho + Cr/NAA smoothly followed the complicated geometry of the tumor inside the gel-based phantom. The results showed that the single-voxel spectra produced by the PRESS pulse sequence with the TE of 144 ms, 512 data-points, and minimum post-processings of water suppression, eddy current correction, and baseline correction can be used for target delineation. CONCLUSION: This study suggests that SVS data can be used to aid target delineation by using a mathematical approach. This can enable a wider use of MR-derived information in radiotherapy. ADVANCES IN KNOWLEDGE: To the best of our knowledge, until now, 2D or 3D MRSI data provided from 3T MRI scanners have been used for MRS-based radiotherapy treatment planning. However, there are a lot of centers that are equipped to 1.5 T MRI scanners and some of them just equipped to SVS. This study introduces a mathematical approach to help these centers to take the benefits of MRS-based treatment planning.

Effects of the attenuation correction and reconstruction method parameters on conventional cardiac dynamic SPECT.

Nuclear cardiology has not witnessed development of new tracers or hardware for many years. Hence there is a need for the development of improvised techniques. Dynamic cardiac single photon emission computed tomography (SPECT) is one such technique that has a potential to overcome the limitations of conventional myocardial SPECT including the absolute quantification of blood flow. The main goal of this study is to evaluate the effect of attenuation correction (AC) on estimation of the washout parameters extracted from dynamic SPECT using a conventional protocol. The effect of the postprocessing on quantitative evaluation of dynamic SPECT is also assessed.A physical phantom was employed to physically simulate the dynamic behavior of a heart in the thorax. Using a dual detector SPECT system, 180° tomographic data in every 90 seconds were acquired. The SPECT data were reconstructed using ordered subset expectation maximization (OSEM) method while different iterations and a Butterworth filter with different cut-off frequencies were applied. Estimated washout parameter of the time activity curves (TACs) was compared with applying AC or without it.Results show that AC can improve the bias of computed washout parameter in normal regions (average bias reduction in normal ROI: 7%). Moreover, the postreconstruction filtering and reducing the number of iterations in reconstructing phase can reduce the variance of the computed washout values in normal regions (from 3.99% for cut-off frequency 0.5 cycle/cm and 32 times update in OSEM to 2.05% for cut-off frequency 0.35 cycle/cm and 16 times update in OSEM). They also reduce the actual size of the defect region (13% reduction in defect extent for above change in reconstruction parameters).According to the results, the AC and postprocessing filtration can directly affect the standard deviation of washout value acquired by cardiac dynamic SPECT. These parameters also showed a direct effect on the defect extent in final results. The study showed that the AC may partly improve the bias of calculated normal washout value. The effect of attenuation correction on the defective washout value could not be answered comprehensively in this paper.

The effect of tandem-ovoid titanium applicator on points A, B, bladder, and rectum doses in gynecological brachytherapy using 192Ir.

PURPOSE: The dosimetry procedure by simple superposition accounts only for the self-shielding of the source and does not take into account the attenuation of photons by the applicators. The purpose of this investigation is an estimation of the effects of the tandem and ovoid applicator on dose distribution inside the phantom by MCNP5 Monte Carlo simulations. MATERIAL AND METHODS: In this study, the superposition method is used for obtaining the dose distribution in the phantom without using the applicator for a typical gynecological brachytherapy (superposition-1). Then, the sources are simulated inside the tandem and ovoid applicator to identify the effect of applicator attenuation (superposition-2), and the dose at points A, B, bladder, and rectum were compared with the results of superposition. The exact dwell positions, times of the source, and positions of the dosimetry points were determined in images of a patient and treatment data of an adult woman patient from a cancer center. The MCNP5 Monte Carlo (MC) code was used for simulation of the phantoms, applicators, and the sources. RESULTS: The results of this study showed no significant differences between the results of superposition method and the MC simulations for different dosimetry points. The difference in all important dosimetry points was found to be less than 5%. CONCLUSIONS: According to the results, applicator attenuation has no significant effect on the calculated points dose, the superposition method, adding the dose of each source obtained by the MC simulation, can estimate the dose to points A, B, bladder, and rectum with good accuracy.

Developing an Optimum Protocol for Thermoluminescence Dosimetry with GR-200 Chips using Taguchi Method.

Thermoluminescence dosimetry (TLD) is a powerful technique with wide applications in personal, environmental and clinical dosimetry. The optimum annealing, storage and reading protocols are very effective in accuracy of TLD response. The purpose of this study is to obtain an optimum protocol for GR-200; LiF: Mg, Cu, P, by optimizing the effective parameters, to increase the reliability of the TLD response using Taguchi method. Taguchi method has been used in this study for optimization of annealing, storage and reading protocols of the TLDs. A number of 108 GR-200 chips were divided into 27 groups, each containing four chips. The TLDs were exposed to three different doses, and stored, annealed and read out by different procedures as suggested by Taguchi Method. By comparing the signal-to-noise ratios the optimum dosimetry procedure was obtained. According to the results, the optimum values for annealing temperature (°C), Annealing Time (s), Annealing to Exposure time (d), Exposure to Readout time (d), Pre-heat Temperature (°C), Pre-heat Time (s), Heating Rate (°C/s), Maximum Temperature of Readout (°C), readout time (s) and Storage Temperature (°C) are 240, 90, 1, 2, 50, 0, 15, 240, 13 and -20, respectively. Using the optimum protocol, an efficient glow curve with low residual signals can be achieved. Using optimum protocol obtained by Taguchi method, the dosimetry can be effectively performed with great accuracy.