• Progress in Medical Physics
• /
• v.26 no.1
• /
• pp.28-35
• /
• 2015

In this study, we investigated the effects of an iterative reconstruction algorithm and an automatic exposure control (AEC) technique on image quality and radiation dose through phantom experiments with coronary computed tomography (CT) angiography protocols. We scanned the AAPM CT performance phantom using 320 multi-detector-row CT. At the tube voltages of 80, 100, and 120 kVp, the scanning was repeated with two settings of the AEC technique, i.e., with the target standard deviations (SD) values of 33 (the higher tube current) and 44 (the lower tube current). The scanned projection data were reconstructed also in two ways, with the filtered back projection (FBP) and with the iterative reconstruction technique (AIDR-3D). The image quality was evaluated quantitatively with the noise standard deviation, modulation transfer function, and the contrast to noise ratio (CNR). More specifically, we analyzed the influences of selection of a tube voltage and a reconstruction algorithm on tube current modulation and consequently on radiation dose. Reduction of image noise by the iterative reconstruction algorithm compared with the FBP was revealed eminently, especially with the lower tube current protocols, i.e., it was decreased by 46% and 38%, when the AEC was established with the lower dose (the target SD=44) and the higher dose (the target SD=33), respectively. As a side effect of iterative reconstruction, the spatial resolution was decreased by a degree that could not mar the remarkable gains in terms of noise reduction. Consequently, if coronary CT angiogprahy is scanned and reconstructed using both the automatic exposure control and iterative reconstruction techniques, it is anticipated that, in comparison with a conventional acquisition method, image noise can be reduced significantly with slight decrease in spatial resolution, implying clinical advantages of radiation dose reduction, still being faithful to the ALARA principle.

• Progress in Medical Physics
• /
• v.20 no.4
• /
• pp.298-307
• /
• 2009

By using the MR T2 map technique, this study intends, first, to measure the change of T2 values of cartilage between healthy people and patients with osteoarthritis and, second, to assess the form and the damage of cartilage in the knee-joint, through which this study would consider the utility of the T2 map technique. Thirty healthy people were selected based on their clinical history and current status and another thirty patients with osteoarthritis of the knee who were screened by simple X-ray from November 2007 to December 2008 were selected. Their T2 Spin Echo (SE hereafter) images for the cartilage of the knee joint were collected by using the T2 SE sequence, one of the multi-echo methods (TR: 1,000 ms; TE values: 6.5, 13, 19.5, 26, 32.5. 40, 45.5, 52). Based on these images, the changes in the signal intensity (SI hereafter) for each section of the cartilage of the knee joint were measured, which yielded average values of T2 through the Origin 7.0 Professional (Northampton, MA 01060 USA). With these T2s, the independent samples T-test was performed by SPSS Window version 12.0 to run the quantitative analysis and to test the statistical significance between the healthy group and the patient group. Closely looking at T2 values for each anterior and lateral articular cartilage of the sagittal plane and the coronal plane, in the sagittal plane, the average T2 of the femoral cartilage in the patient group with arthritis of the knee ($42.22{\pm}2.91$) was higher than the average T2 of the healthy group ($36.26{\pm}5.01$). Also, the average T2 of the tibial cartilage in the patient group ($43.83{\pm}1.43$) was higher than the average T2 in the healthy group ($36.45{\pm}3.15$). In the case of the coronal plane, the average T2 of the medial femoral cartilage in the patient group ($45.65{\pm}7.10$) was higher than the healthy group ($36.49{\pm}8.41$) and so did the average T2 of the anterior tibial cartilage (i.e., $44.46{\pm}3.44$ for the patient group vs. $37.61{\pm}1.97$ for the healthy group). As for the lateral femoral cartilage in the coronal plane, the patient group displayed the higher T2 ($43.41{\pm}4.99$) than the healthy group did ($37.64{\pm}4.02$) and this tendency was similar in the lateral tibial cartilage (i.e., $43.78{\pm}8.08$ for the patient group vs. $36.62{\pm}7.81$ for the healthy group). Along with the morphological MR imaging technique previously used, the T2 map technique seems to help patients with cartilage problems, in particular, those with the arthritis of the knee for early diagnosis by quantitatively analyzing the structural and functional changes of the cartilage.

• Progress in Medical Physics
• /
• v.21 no.1
• /
• pp.52-59
• /
• 2010

To examine among patients with vertebral compression fracture the extent to which signal-to-noise ratio (SNR) and Apparent Diffusion Coefficient (ADC) values at the lumbar vertebral compression fracture site vary on diffusion-weighted MR images according to varying b values on the 1.5T MR device. Diffusion-weighted MR images of 30 patients with compression fracture due to chronic osteoporosis who underwent vertebral MRI from Jan. 2008 to Nov. 2009 were respectively obtained using a 1.5-T MR scanner with the b values increased from 400, 600, 800, 1,000 to $1,200\;s/mm^2$. For diffusion-weighted MR images with different b values, the signal-to-noise ratio (SNR) was assessed at three sites: the site of compression fracture of the lumbar vertebral body at L1 to L5, and both the upper and lower discs of the said fracture site, while for ADC map images with different b values, the SNR and ADC were respectively assessed at those three sites. As a quantitative analysis, diffusion-weighted MR images and ADC map images with b value of $400\;s/mm^2$ (the base b values) were respectively compared with the corresponding images with each different b value. As far as qualitative analysis is concerned, for both diffusion-weighted MR and ADC map images with b value of $400\;s/mm^2$, the extent to which signal intensity values obtained at the site of compression fracture of the lumbar vertebral body at L1 to L5 vary according to the increasing b values were examined. The quantitative analysis found that for both diffusion-weighted MR and ADC map images, as the b values increased, the SNR were relatively lowered at all the three sites, compared to the base b value. Also, it was found that as the b values increased, ADC valueswere relatively lowered at all the three sites on ADC map images. On the other hand, the qualitative analysis found that as the b values increased to more than $400\;s/mm^2$, the signal intensity gradually decreased at all the sites, while at the levels of more than $1,000\;s/mm^2$, severe image noises appeared at all of the three sites. In addition, higher signal intensity was found at the site of compression fracture of the lumbar vertebral body than at the discs. Findings showed that with the b value being increased, both the signal-to-noise ratio (SNR) and Apparent Diffusion Coefficient (ADC) values gradually decreased at all the sites of the lumbar vertebral compression fracture and both the upper and lower discs of the fracture site, suggesting that there is a possibility of a wider range of applications to assessment of various vertebral pathologies by utilizing multi b values in the diffusion-weighted MRI examination.

• Progress in Medical Physics
• /
• v.22 no.1
• /
• pp.42-51
• /
• 2011

Nuclear medicine images (SPECT, PET) were widely used tool for assessment of myocardial viability and perfusion. However it had difficult to define accurate myocardial infarct region. The purpose of this study was to investigate methodological approach for automatic measurement of rat myocardial infarct size using polar map with adaptive threshold. Rat myocardial infarction model was induced by ligation of the left circumflex artery. PET images were obtained after intravenous injection of 37 MBq $^{18}F$-FDG. After 60 min uptake, each animal was scanned for 20 min with ECG gating. PET data were reconstructed using ordered subset expectation maximization (OSEM) 2D. To automatically make the myocardial contour and generate polar map, we used QGS software (Cedars-Sinai Medical Center). The reference infarct size was defined by infarction area percentage of the total left myocardium using TTC staining. We used three threshold methods (predefined threshold, Otsu and Multi Gaussian mixture model; MGMM). Predefined threshold method was commonly used in other studies. We applied threshold value form 10% to 90% in step of 10%. Otsu algorithm calculated threshold with the maximum between class variance. MGMM method estimated the distribution of image intensity using multiple Gaussian mixture models (MGMM2, ${\cdots}$ MGMM5) and calculated adaptive threshold. The infarct size in polar map was calculated as the percentage of lower threshold area in polar map from the total polar map area. The measured infarct size using different threshold methods was evaluated by comparison with reference infarct size. The mean difference between with polar map defect size by predefined thresholds (20%, 30%, and 40%) and reference infarct size were $7.04{\pm}3.44%$, $3.87{\pm}2.09%$ and $2.15{\pm}2.07%$, respectively. Otsu verse reference infarct size was $3.56{\pm}4.16%$. MGMM methods verse reference infarct size was $2.29{\pm}1.94%$. The predefined threshold (30%) showed the smallest mean difference with reference infarct size. However, MGMM was more accurate than predefined threshold in under 10% reference infarct size case (MGMM: 0.006%, predefined threshold: 0.59%). In this study, we was to evaluate myocardial infarct size in polar map using multiple Gaussian mixture model. MGMM method was provide adaptive threshold in each subject and will be a useful for automatic measurement of infarct size.

• Progress in Medical Physics
• /
• v.27 no.1
• /
• pp.31-36
• /
• 2016

The paper discusses radiation dose of dual energy CT on which copper modulation layer, is mounted in order to improve diagnostic performance of the dual energy CT. The radiation dose is estimated using MCNPX and its results are compared with that of the conventional dual energy CT system. CT X-ray spectra of 80 and 120 kVp, which are usually used for thorax, abdominal, head, and neck CT scans, were generated by the SPEC78 code and were used for the source specification 'SDEF' card for MCNPX dose modeling. The copper modulation layer was located 20 cm away from a source covering half of the X-ray window. The radiation dose was measured as changing its thickness from 0.5 to 2.0 mm at intervals of 0.5 mm. Since the MCNPX tally provides only normalized values to a single particle, the dose conversion coefficients of F6 tally for the modulation layer-based dual energy CBCT should be calculated for matching the modeling results into the actual dose. The dose conversion coefficient is $7.2*10^4cGy/output$ that is obtained from dose calibration curve between F6 tally and experimental results in which GAFCHORMIC EBT3 films were exposed by an already known source. Consequently, the dose of the modulation layer-based dual energy cone beam CT is 33~40% less than that of the single energy CT system. On the basis of the results, it is considered that scattered dose produced by the copper modulation layer is very small. It shows that the modulation layer-based dual energy CBCT system can effectively reduce radiation dose, which is the major disadvantage of established dual energy CT.

• Progress in Medical Physics
• /
• v.18 no.3
• /
• pp.179-185
• /
• 2007

The purpose of this study is to investigate the spectra of a magnetic resonance spectroscopy (MRS) in accordance with the variance of TE and the volumes of metabolites in a localized voxel for the quality assurance using a designed single voxel spectroscopy QA phantom. Because a cone-shade phantom is designed as the volume of metabolite in a localized voxel is changeable, we try to analyze the peaks of each metabolite (NAA, Cr, Cho, Lac, etc.) in accordance with metabolite volume in a localized voxel as well as echo time (TE). All data were obtained using a 3T MRI/MRS machine and analyzed using $jMRUI^{(R)}$. The results of this study show that TE is in inverse proportion to the noise of MRS and the longer TE and the less metabolite volume in the localized voxel, the peak intensities of each metabolite decrease. In case of the lactate, its peak was observed on the all TE only if the greatest metabolite is included in the localized voxel. Then, the intensity of a metabolite is more sensitive to the metabolite volume in the localized voxel than the TE. These obtained in vitro MRS data is provide the guideline that is important for in vivo metabolite quantification. But, in the edge of cone-shape vial air bubbles were observed and spectrum could not obtained. Therefore our cone-shape MRS phantom needs to be modified in order to solve these problems.

• Progress in Medical Physics
• /
• v.20 no.4
• /
• pp.269-276
• /
• 2009

Radiation treatment techniques using photon beam such as three-dimensional conformal radiation therapy (3D-CRT) as well as intensity modulated radiotherapy treatment (IMRT) demand accurate dose calculation in order to increase target coverage and spare healthy tissue. Both jaw collimator and multi-leaf collimators (MLCs) for photon beams have been used to achieve such goals. In the Pinnacle3 treatment planning system (TPS), which we are using in our clinics, a set of model parameters like jaw collimator transmission factor (JTF) and MLC transmission factor (MLCTF) are determined from the measured data because it is using a model-based photon dose algorithm. However, model parameters obtained by this auto-modeling process can be different from those by direct measurement, which can have a dosimetric effect on the dose distribution. In this paper we estimated JTF and MLCTF obtained by the auto-modeling process in the Pinnacle3 TPS. At first, we obtained JTF and MLCTF by direct measurement, which were the ratio of the output at the reference depth under the closed jaw collimator (MLCs for MLCTF) to that at the same depth with the field size $10{\times}10\;cm^2$ in the water phantom. And then JTF and MLCTF were also obtained by auto-modeling process. And we evaluated the dose difference through phantom and patient study in the 3D-CRT plan. For direct measurement, JTF was 0.001966 for 6 MV and 0.002971 for 10 MV, and MLCTF was 0.01657 for 6 MV and 0.01925 for 10 MV. On the other hand, for auto-modeling process, JTF was 0.001983 for 6 MV and 0.010431 for 10 MV, and MLCTF was 0.00188 for 6 MV and 0.00453 for 10 MV. JTF and MLCTF by direct measurement were very different from those by auto-modeling process and even more reasonable considering each beam quality of 6 MV and 10 MV. These different parameters affect the dose in the low-dose region. Since the wrong estimation of JTF and MLCTF can lead some dosimetric error, comparison of direct measurement and auto-modeling of JTF and MLCTF would be helpful during the beam commissioning.

• Progress in Medical Physics
• /
• v.19 no.4
• /
• pp.209-218
• /
• 2008

According to improved radiation therapy technology such as IMRT and proton therapy, the accuracy of patient alignment system is more emphasized and IGRT is dominated research field in radiation oncology. We proposed to study the feasibility of cone-beam CT system using simple x-ray imaging systems for image guided proton therapy at National Cancer Center. 180 projection views ($2,304{\times}3,200$, 14 bit with 127 ${\mu}m$ pixel pitch) for the geometrical calibration phantom and humanoid phantoms (skull, abdomen) were acquired with $2^{\circ}$ step angle using x-ray imaging system of proton therapy gantry room ($360^{\circ}$ for 1 rotation). The geometrical calibration was performed for misalignments between the x-ray source and the flat-panel detector, such as distances and slanted angle using available algorithm. With the geometrically calibrated projection view, Feldkamp cone-beam algorithm using Ram-Lak filter was implemented for CBCT reconstruction images for skull and abdomen phantom. The distance from x-ray source to the gantry isocenter, the distance from the flat panel to the isocenter were calculated as 1,517.5 mm, 591.12 mm and the rotated angle of flat panel detector around x-ray beam axis was considered as $0.25^{\circ}$. It was observed that the blurring artifacts, originated from the rotation of the detector, in the reconstructed toomographs were significantly reduced after the geometrical calibration. The demonstrated CBCT images for the skull and abdomen phantoms are very promising. We performed the geometrical calibration of the large gantry rotation system with simple x-ray imaging devices for CBCT reconstruction. The CBCT system for proton therapy will be used as a main patient alignment system for image guided proton therapy.

• Progress in Medical Physics
• /
• v.20 no.3
• /
• pp.132-138
• /
• 2009

Studies on target motion in 4-dimensional radiotherapy are being world-widely conducted to enhance treatment record and protection of normal organs. Prediction of tumor motion might be very useful and/or essential for especially free-breathing system during radiation delivery such as respiratory gating system and tumor tracking system. Neural network is powerful to express a time series with nonlinearity because its prediction algorithm is not governed by statistic formula but finds a rule of data expression. This study intended to assess applicability of neural network method to predict tumor motion in 4-dimensional radiotherapy. Scaled Conjugate Gradient algorithm was employed as a learning algorithm. Considering reparation data for 10 patients, prediction by the neural network algorithms was compared with the measurement by the real-time position management (RPM) system. The results showed that the neural network algorithm has the excellent accuracy of maximum absolute error smaller than 3 mm, except for the cases in which the maximum amplitude of respiration is over the range of respiration used in the learning process of neural network. It indicates the insufficient learning of the neural network for extrapolation. The problem could be solved by acquiring a full range of respiration before learning procedure. Further works are programmed to verify a feasibility of practical application for 4-dimensional treatment system, including prediction performance according to various system latency and irregular patterns of respiration.

• Progress in Medical Physics
• /
• v.9 no.4
• /
• pp.227-245
• /
• 1998

A most appropriate model of 3-D conformal radiotherapy has been induced by clinical evaluation and animal study, and therapeutic gains were evaluated by numerical equation of tumor control probability(TCP) and normal tissue complication probability (NTCP). The radiation dose to the tumor and the adjacent normal organs was accurately evaluated and compared using the dose volume histogram(DVH). The TCP and NTCP was derived from the distribution of given dosage and irradiated volume, and these numbers were used as the biological index for the assessment of the treatment effects. Ten patients with liver disease have been evaluated and 3 dogs were sacrificed for this study. Based on the 3-D images of the tumor and adjacent organs, the optimum radiation dose and the projection direction which could maximize the radiation effect while minimizing the effects to the adjacent organs could be decided. 3). The most effective collimation for the normal adjacent organs was made through the beams eye view with the use of multileaf collimator. When the dose was increased from 50Gy to 70Gy, the TCP for the conventional 2-port radiation and the 5-port multidimensional therapy was 0.982 and 0.995 respectively, while the NTCP was 0.725 and 0.142 respectively, suggesting that the 3-D conformal radiotherapy might be the appropriate therapy to apply sufficient radiation dose to the tumor while minimizing the damages to the normal areas of the liver. Positive correlation was observed between the NTCP and the actual complication of the normal liver in the animal study. The present study suggest that the use of 3-D conformal radiotherapy and the application of the mathematical models of TCP and NTCP may provide the improvements in the treatment of hepatoma with enhanced results.