• Nuclear Engineering and Technology
• /
• v.52 no.7
• /
• pp.1545-1556
• /
• 2020

Recently, the International Commission on Radiological Protection (ICRP) has developed the Mesh-type Reference Computational Phantoms (MRCPs) for adult male and female to overcome the limitations of the current Voxel-type Reference Computational Phantoms (VRCPs) described in ICRP Publication 110 due to the limited voxel resolutions and the nature of voxel geometry. In our previous study, the MRCPs were used to calculate the dose coefficients (DCs) for idealized external exposures of photons and electrons. The present study is an extension of the previous study to include three additional particles (i.e., neutrons, protons, and helium ions) into the DC library by conducting Monte Carlo radiation transport simulations with the Geant4 code. The calculated MRCP DCs were compared with the reference DCs of ICRP Publication 116 which are based on the VRCPs, to appreciate the impact of the new reference phantoms on the DC values. We found that the MRCP DCs of organ/tissue doses and effective doses were generally similar to the ICRP-116 DCs for neutrons, whereas there were significant DC differences up to several orders of magnitude for protons and helium ions due mainly to the improved representation of the detailed anatomical structures in the MRCPs over the VRCPs.

• Investigative Magnetic Resonance Imaging
• /
• v.23 no.1
• /
• pp.38-45
• /
• 2019

Purpose: To demonstrate the high-resolution numerical simulation of the respiration-induced dynamic $B_0$ shift in the head using generalized susceptibility voxel convolution (gSVC). Materials and Methods: Previous dynamic $B_0$ simulation research has been limited to low-resolution numerical models due to the large computational demands of conventional Fourier-based $B_0$ calculation methods. Here, we show that a recently-proposed gSVC method can simulate dynamic $B_0$ maps from a realistic breathing human body model with high spatiotemporal resolution in a time-efficient manner. For a human body model, we used the Extended Cardiac And Torso (XCAT) phantom originally developed for computed tomography. The spatial resolution (voxel size) was kept isotropic and varied from 1 to 10 mm. We calculated $B_0$ maps in the brain of the model at 10 equally spaced points in a respiration cycle and analyzed the spatial gradients of each of them. The results were compared with experimental measurements in the literature. Results: The simulation predicted a maximum temporal variation of the $B_0$ shift in the brain of about 7 Hz at 7T. The magnitudes of the respiration-induced $B_0$ gradient in the x (right/left), y (anterior/posterior), and z (head/feet) directions determined by volumetric linear fitting, were < 0.01 Hz/cm, 0.18 Hz/cm, and 0.26 Hz/cm, respectively. These compared favorably with previous reports. We found that simulation voxel sizes greater than 5 mm can produce unreliable results. Conclusion: We have presented an efficient simulation framework for respiration-induced $B_0$ variation in the head. The method can be used to predict $B_0$ shifts with high spatiotemporal resolution under different breathing conditions and aid in the design of dynamic $B_0$ compensation strategies.

• Journal of the Korean Society of Radiology
• /
• v.13 no.4
• /
• pp.565-572
• /
• 2019

In this study, shielding analysis of material and thickness of 3D printer filaments was performed for the manufacture of custom shielding by radiation workers during outdoor radiographic test. The shielding was attached to the ICRU Slab Phantom after selecting the voxel source $^{192}Ir$ and $^{75}Se$ through simulation using MCNPX, and the distance between the source and the slab Phantom was set at 100 cm. The 12 shielding materials were divided into 5 mm units up to 200 mm from the absence of shielding materials to evaluate the energy absorbed per unit mass of each shielding material. The results showed that the shielding effect was high in the order of ABS + Tungsten, ABS + Bismuth, PLA + Copper, PLA + Iron from all sources of radiographic test. However, compared to lead, the shielding effect was somewhat lower. Based on this study in the future, further study of the atomic number and the high density filament material is necessary.

• Journal of Biomedical Engineering Research
• /
• v.16 no.4
• /
• pp.395-402
• /
• 1995

The intravoxel spin phases in magnetic resonance imaging (MRI) usually vary due to susceptibility differences of materials to be imaged. The phase variation in the voxel results in a reduction of the signal intensity. This signal intensity reduction is known as the susceptibility effect in MRI and has been studied extensively. In this paper, a new spectral decomposition technique Is proposed and the signal change due to the susceptibility effect can be analyzed. A pulse sequence for the spectral decomposition of the susceptibility was developed and applied to susceptibility imaging of venous blood possessing paramagnetic properties. The computer simulations and their corresponding experimental results obtained using both a phantom and human volunteers are reported. Key words : susceptibility effect in MRI : spectral decomposition of susceptibility effect.

• Nuclear Engineering and Technology
• /
• v.47 no.4
• /
• pp.472-478
• /
• 2015

The patient dose incurred from diagnostic procedures during advanced radiotherapy has become an important issue. Many researchers in medical physics are using computational simulations to calculate complex parameters in experiments. However, extended computation times make it difficult for personal computers to run the conventional Monte Carlo method to simulate radiological images with high-flux photons such as images produced by computed tomography (CT). To minimize the computation time without degrading imaging quality, we applied a deterministic adaptation to the Monte Carlo calculation and verified its effectiveness by simulating CT image reconstruction for an image evaluation phantom (Catphan; Phantom Laboratory, New York NY, USA) and a human-like voxel phantom (KTMAN-2) (Los Alamos National Laboratory, Los Alamos, NM, USA). For the deterministic adaptation, the relationship between iteration numbers and the simulations was estimated and the option to simulate scattered radiation was evaluated. The processing times of simulations using the adaptive method were at least 500 times faster than those using a conventional statistical process. In addition, compared with the conventional statistical method, the adaptive method provided images that were more similar to the experimental images, which proved that the adaptive method was highly effective for a simulation that requires a large number of iterations-assuming no radiation scattering in the vicinity of detectors minimized artifacts in the reconstructed image.

• Journal of Biomedical Engineering Research
• /
• v.41 no.1
• /
• pp.28-34
• /
• 2020

In this paper, we introduce motion artifact reduction algorithm for interleaved MRI using an advanced 3D approximation algorithm. The motion artifact framework of this paper is data corrected by post-processing with a new 3-D approximation algorithm which uses data structure for each voxel. In this study, we simulate and evaluate our algorithm using Shepp-Logan phantom and T1-MRI template for both scattered dataset and uniform dataset. We generated motion artifact using random generated motion parameters for the interleaved MRI. In simulation, we use image coregistration by SPM12 (https://www.fil.ion.ucl.ac.uk/spm/) to estimate the motion parameters. The motion artifact correction is done with using full dataset with estimated motion parameters, as well as use only one half of the full data which is the case when the half volume is corrupted by severe movement. We evaluate using numerical metrics and visualize error images.

• Progress in Medical Physics
• /
• v.14 no.3
• /
• pp.203-210
• /
• 2003

PET (positron emission tomography) permits the investigation of physiological and biochemical processes in vivo. The accuracy of quantifying PET data is affected by its finite spatial resolution, which causes partial volume effects. In this study, we developed a method for partial volume correction using Hoffman phantom PET and MR data, and applied various FWHM (full width at half maximum) levels. We also applied this method to PET images of normal controls and tested for the possibility of clinical application. $^{18}$ F-PET Hoffman phantom images were co-registered to MR slices. The gray matter and white matter regions were then segmented into binary images. Each binary image was convolved by 4, 8, 12, 16 mm FWHM levels. These convolved images of gray and white matter were merged corresponding to the same level of FWHM. The original PET images were then divided by the convolved binary images voxel-by-voxel. These corrected PET images were multiplied by binary images. The corrected PET images were evaluated by analyzing regions of interests, which were drawn on the gray and white matter regions of the original MR image slices. We calculated the ratio of white to gray matter. We also applied this method to the PET images of normal controls. On analyzing the corrected PET images of Hoffman phantom, the ratios of the corrected images increased more than that of the uncorrected images. With the normal controls, the ratio of the corrected images increased more than that of the uncorrected images. The ratio increase of the corrected PET images was lower than that of the corrected phantom PET images. In conclusion, the method developed for partial volume correction in PET data may be clinically applied, although further study may be required for optimal correction.

• Nuclear Engineering and Technology
• /
• v.49 no.7
• /
• pp.1495-1504
• /
• 2017

In a previous study, a set of polygon-mesh (PM)-based skin models including a $50-{\mu}m-thick$ radiosensitive target layer were constructed and used to calculate skin dose coefficients (DCs) for idealized external beams of electrons. The results showed that the calculated skin DCs were significantly different from the International Commission on Radiological Protection (ICRP) Publication 116 skin DCs calculated using voxel-type ICRP reference phantoms that do not include the thin target layer. The difference was as large as 7,700 times for electron energies less than 1 MeV, which raises a significant issue that should be addressed subsequently. In the present study, therefore, as an extension of the initial, previous study, skin DCs for three other particles (photons, protons, and helium ions) were calculated by using the PM-based skin models and the calculated values were compared with the ICRP-116 skin DCs. The analysis of our results showed that for the photon exposures, the calculated values were generally in good agreement with the ICRP-116 values. For the charged particles, by contrast, there was a significant difference between the PM-model-calculated skin DCs and the ICRP-116 values. Specifically, the ICRP-116 skin DCs were smaller than those calculated by the PM models-which is to say that they were under-estimated-by up to ~16 times for both protons and helium ions. These differences in skin dose also significantly affected the calculation of the effective dose (E) values, which is reasonable, considering that the skin dose is the major factor determining effective dose calculation for charged particles. The results of the current study generally show that the ICRP-116 DCs for skin dose and effective dose are not reliable for charged particles.

• Journal of Biomedical Engineering Research
• /
• v.24 no.2
• /
• pp.77-82
• /
• 2003

Virtual colonoscopy is a non-invasive computerized procedure to detect polyps by examining the colon from a CT data set. To fly through the inside of colons. the extraction of a suitable flight-path is necessary to Provide the viewpoint and view direction of a virtual camera. However. manual path extraction by Picking Points is a very time-consuming and difficult task due 1,c, the long and complex shape of colon. Also, existing automatic methods are computationally complex. and tend to generate an improper and/or discontinuous path for complicated regions. In this paper, we propose a fast flight-path generation algorithm using the distance and order maps. The order map Provides all Possible directions of a path. The distance map assigns the Euclidean distance value from each inside voxel to the nearest background voxel. By jointly using these two maps. we can obtain a proper centerline regardless of thickness and curvature of an object. Also, we Propose a simple smoothing technique that guarantees not to collide with the surface of an object. The phantom and real colon data are used for experiments. Experimental results show that for a set of human colon data, the proposed algorithm can provide a smoothened and connected flight-path within a minute on an 800MHz PC. And it is proved that the obtained flight-Path provides successive volume-rendered images satisfactory for virtual navigation.

• Progress in Medical Physics
• /
• v.25 no.3
• /
• pp.139-142
• /
• 2014

The purpose of this study was to develop a system of clinical application of reconstructed dose that includes dose reconstruction, reconstructed dose registration between fractions of treatment, and dose-volume-histogram generation and to demonstrate the system on a deformable prostate phantom. To achieve this purpose, a deformable prostate phantom was embedded into a 20 cm-deep and 40 cm-wide water phantom. The phantom was CT scanned and the anatomical models of prostate, seminal vesicles, and rectum were contoured. A coplanar 4-field intensity modulated radiation therapy (IMRT) plan was used for this study. Organ deformation was simulated by inserting a "transrectal" balloon containing 20 ml of water. A new CT scan was obtained and the deformed structures were contoured. Dose responses in phantoms and electronic portal imaging device (EPID) were calculated by using the XVMC Monte Carlo code. The IMRT plan was delivered to the two phantoms and integrated EPID images were respectively acquired. Dose reconstruction was performed on these images using the calculated responses. The deformed phantom was registered to the original phantom using an in-house developed software based on the Demons algorithm. The transfer matrix for each voxel was obtained and used to correlate the two sets of the reconstructed dose to generate a cumulative reconstructed dose on the original phantom. Forwardly calculated planning dose in the original phantom was compared to the cumulative reconstructed dose from EPID in the original phantom. The prescribed 200 cGy isodose lines showed little difference with respect to the "prostate" and "seminal vesicles", but appreciable difference (3%) was observed at the dose level greater than 210 cGy. In the rectum, the reconstructed dose showed lower volume coverage by a few percent than the plan dose in the dose range of 150 to 200 cGy. Through this study, the system of clinical application of reconstructed dose was successfully developed and demonstrated. The organ deformation simulated in this study resulted in small but observable dose changes in the target and critical structure.