Organ and effective doses in pediatric patients undergoing helical multislice computed tomography examination

As multidetector computed tomography (CT) serves as an increasingly frequent diagnostic modality, radiation risks to patients became a greater concern, especially for children due to their inherently higher radiosensitivity to stochastic radiation damage. Current dose evaluation protocols include the computed tomography dose index (CTDI) or point detector measurements using anthropomorphic phantoms that do not sufficiently provide accurate information of the organ-averaged absorbed dose and corresponding effective dose to pediatric patients. In this study, organ and effective doses to pediatric patients under helical multislice computed tomography (MSCT) examinations were evaluated using an extensive series of anthropomorphic computational phantoms and Monte Carlo radiation transport simulations. Ten pediatric phantoms, five stylized (equation-based) ORNL phantoms (newborn, 1-year, 5-year, 10-year, and 15-year) and five tomographic (voxel-based) UF phantoms (9-month male, 4-year female, 8-year female, 11-year male, and 14-year male) were implemented into MCNPX for simulation, where a source subroutine was written to explicitly simulate the helical motion of the CT x-ray source and the fan beam angle and collimator width. Ionization chamber measurements were performed and used to normalize the Monte Carlo simulation results. On average, for the same tube current setting, a tube potential of 100 kVp resulted in effective doses that were 105% higher than seen at 80 kVp, and 210% higher at 120 kVp regardless of phantom type. Overall, the ORNL phantom series was shown to yield values of effective dose that were reasonably consistent with those of the gender-specific UF phantom series for CT examinations of the head, pelvis, and torso. However, the ORNL phantoms consistently overestimated values of the effective dose as seen in the UF phantom for MSCT scans of the chest, and underestimated values of the effective dose for abdominal CT scans. These discrepancies increased with increasing kVp. Finally, absorbed doses to select radiation sensitive organs such as the gonads, red bone marrow, colon, and thyroid were evaluated and compared between phantom types. Specific anatomical problems identified in the stylized phantoms included excessive pelvic shielding of the ovaries in the female phantoms, enhanced red bone marrow dose to the arms and rib cage for chest exams, an unrealistic and constant torso thickness resulting in excessive x-ray attenuation in the regions of the abdominal organs, and incorrect positioning of the thyroid within the stylized phantom neck resulting in insufficient shielding by clavicles and scapulae for lateral beam angles. To ensure more accurate estimates of organ absorbed dose in multislice CT, it is recommended that voxel-based phantoms, potentially tailored to individual body morphometry, be utilized in any future prospective epidemiological studies of medically exposed children.     

High temporal resolution for multislice helical computed tomography

Multislice helical computed tomography (CT) substantially reduces scanning time. However, the temporal resolution of individual images is still insufficient for imaging rapidly moving organs such as the heart and adjacent pulmonary vessels. It may, in some cases, be worse than with current single-slice helical CT. The purpose of this study is to describe a novel image reconstruction algorithm to improve temporal resolution in multislice helical CT, and to evaluate its performance against existing algorithms. The proposed image reconstruction algorithm uses helical interpolation followed by data weighting based on the acquisition time. The temporal resolution, the longitudinal (z-axis) spatial resolution, the image noise, and the in-plane image artifacts created by a moving phantom were compared with those from the basic multislice helical reconstruction (helical filter interpolation, HFI) algorithm and the basic single-slice helical reconstruction algorithm (180 degrees linear interpolation, 180LI) using computer simulations. Computer simulation results were verified with CT examinations of the heart and lung vasculature using a 0.5 second multislice scanner. The temporal resolution of HFI algorithm varies from 0.28 and 0.86 s, depending on helical pitch. The proposed method improves the resolution to a constant value of 0.29 s, independent of pitch, allowing moving objects to be imaged with reduced blurring or motion artifacts. The spatial (z) resolution was slightly worse than with the HFI algorithm; the image noise was worse than with the HFI algorithm but was comparable to axial (step-and-shoot) CT. The proposed method provided sharp images of the moving objects, portraying the anatomy accurately. The proposed algorithm for multislice helical CT allowed us to obtain CT images with high temporal resolution. It may improve the image quality of clinical cardiac, lung, and vascular CT imaging.     

Estimation of effective doses to adult and pediatric patients from multislice computed tomography: A method based on energy imparted

The purpose of this study is to provide a method and required data for the estimation of effective dose (E) values to adult and pediatric patients from computed tomography (CT) scans of the head, chest abdomen, and pelvis, performed on multi-slice scanners. Mean section radiation dose (dm) to cylindrical water phantoms of varying radius normalized over CT dose index free-in-air (CTDIF) were calculated for the head and body scanning modes of a multislice scanner with use of Monte Carlo techniques. Patients were modeled as equivalent water phantoms and the energy imparted (epsilon) to simulated pediatric and adult patients was calculated on the basis of measured CTDI(F) values. Body region specific energy imparted to effective dose conversion coefficients (E/epsilon) for adult male and female patients were generated from previous data. Effective doses to patients aged newborn to adult were derived for all available helical and axial beam collimations, taking into account age specific patient mass and scanning length. Depending on high voltage, body region, and patient sex, E/epsilon values ranged from 0.008 mSv/mJ for head scans to 0.024 mSv/mJ for chest scans. When scanned with the same technique factors as the adults, pediatric patients absorb as little as 5% of the energy imparted to adults, but corresponding effective dose values are up to a factor of 1.6 higher. On average, pediatric patients absorb 44% less energy per examination but have a 24% higher effective dose, compared with adults. In clinical practice, effective dose values to pediatric patients are 2.5 to 10 times lower than in adults due to the adaptation of tube current. A method is provided for the calculation of effective dose to adult and pediatric patients on the basis of individual patient characteristics such as sex, mass, dimensions, and density of imaged anatomy, and the technical features of modern multislice scanners. It allows the optimum selection of scanning parameters regarding patient doses at CT.     

Organ and effective doses in infants undergoing upper gastrointestinal (UGI) fluoroscopic examination

To provide more detailed data on organ and effective doses in digital upper gastrointestinal (UGI) fluoroscopy studies of newborns and infants, the present study was conducted employing the time-sequence videotape-analysis technique used in a companion study of newborn and infant voiding cystourethrograms (VCUG). This technique was originally pioneered [O. H. Suleiman, J. Anderson, B. Jones, G. U. Rao, and M. Rosenstein, Radiology 178, 653-658 (1991)] for adult UGI examinations. Individual video frames were analyzed to include combinations of field size, field center, x-ray projection, image intensifier, and magnification mode. Additionally, the peak tube potential and the mA or mAs values for each segment/subsegment or digital photospot were recorded for both the fluoroscopic and radiographic modes of operation. The data from videotape analysis were then used in conjunction with a patient-scalable newborn tomographic computational phantom to report both organ and effective dose values via Monte Carlo radiation transport. The study includes dose estimates for five simulated UGI examinations representative of patients ranging from three to six months of age. Effective dose values for UGI examinations ranged from 1.17 to 6.47 mSv, with a mean of 3.14 mSv and a large standard deviation of 2.15 mSv. The colon, lungs, stomach, liver, and esophagus absorbed doses in sum were found to constitute between 63 and 75% of the effective dose in these UGI studies. Representing 23-30% of the effective dose, the lungs were found to be the most significant organ in the effective dose calculation. Approximately 80-95% of the effective dose is contributed by the dynamic fluoroscopy segments with larger percentages found in longer studies. The mean effective dose for newborn UGI examinations was not found to be statistically different from that seen in newborn VCUG examinations.     

Effective doses to patients undergoing thoracic computed tomography examinations

The purpose of this study was to investigate how x-ray technique factors and effective doses vary with patient size in chest CT examinations. Technique factors (kVp, mAs, section thickness, and number of sections) were recorded for 44 patients who underwent a routine chest CT examination. Patient weights were recorded together with dimensions and mean Hounsfield unit values obtained from representative axial CT images. The total mass of directly irradiated patient was modeled as a cylinder of water to permit the computation of the mean patient dose and total energy imparted for each chest CT examination. Computed values of energy imparted during the chest CT examination were converted into effective doses taking into account the patient weight. Patient weights ranged from 4.5 to 127 kg, and half the patients in this study were children under 18 years of age. All scans were performed at 120 kVp with a 1 s scan time. The selected tube current showed no correlation with patient weight (r2=0.06), indicating that chest CT examination protocols do not take into account for the size of the patient. Energy imparted increased with increasing patient weight, with values of energy imparted for 10 and 70 kg patients being 85 and 310 mJ, respectively. The effective dose showed an inverse correlation with increasing patient weight, however, with values of effective dose for 10 and 70 kg patients being 9.6 and 5.4 mSv, respectively. Current CT technique factors (kVp/mAs) used to perform chest CT examinations result in relatively high patient doses, which could be reduced by adjusting technique factors based on patient size.     

The effect of z overscanning on patient effective dose from multidetector helical computed tomography examinations

z overscanning in multidetector (MD) helical CT scanning is prerequisite for the interpolation of acquired data required during image reconstruction and refers to the exposure of tissues beyond the boundaries of the volume to be imaged. The aim of the present study was to evaluate the effect of z overscanning on the patient effective dose from helical MD CT examinations. The Monte Carlo N-particle radiation transport code was employed in the current study to simulate CT exposure. The validity of the Monte Carlo simulation was verified by (a) a comparison of calculated and measured standard computed tomography dose index (CTDI) dosimetric data, and (b) a comparison of calculated and measured dose profiles along the z axis. CTDI was measured using a pencil ionization chamber and head and body CT phantoms. Dose profiles along the z axis were obtained using thermoluminescence dosimeters. A commercially available mathematical anthropomorphic phantom was used for the estimation of effective doses from four standard CT examinations, i.e., head and neck, chest, abdomen and pelvis, and trunk studies. Data for both axial and helical modes of operation were obtained. In the helical mode, z overscanning was taken into account. The calculated effective dose from a CT exposure was normalized to CTDI(free in air). The percentage differences in the normalized effective dose between contiguous axial and helical scans with pitch = 1, may reach 13.1%, 35.8%, 29.0%, and 21.5%, for head and neck, chest, abdomen and pelvis, and trunk studies, respectively. Given that the same kilovoltage and tube load per rotation were used in both axial and helical scans, the above differences may be attributed to z overscanning. For helical scans with pitch = 1, broader beam collimation is associated with increased z overscanning and consequently higher normalized effective dose value, when other scanning parameters are held constant. For a given beam collimation, the selection of a higher value of reconstructed image slice width increases the normalized effective dose. In conclusion, z overscanning may significantly affect the patient effective dose from CT examinations performed on MD CT scanners. Therefore, an estimation of the patient effective dose from MD helical CT examinations should always take into consideration the effect of z overscanning.     

General surface reconstruction for cone-beam multislice spiral computed tomography

A new family of cone-beam reconstruction algorithm, the General Surface Reconstruction (GSR), is proposed and formulated in this paper for multislice spiral computed tomography (CT) reconstructions. It provides a general framework to allow the reconstruction of planar or nonplanar surfaces on a set of rebinned short-scan parallel beam projection data. An iterative surface formation method is proposed as an example to show the possibility to form nonplanar reconstruction surfaces to minimize the adverse effect between the collected cone-beam projection data and the reconstruction surfaces. The improvement in accuracy of the nonplanar surfaces over planar surfaces in the two-dimensional approximate cone-beam reconstructions is mathematically proved and demonstrated using numerical simulations. The proposed GSR algorithm is evaluated by the computer simulation of cone-beam spiral scanning geometry and various mathematical phantoms. The results demonstrate that the GSR algorithm generates much better image quality compared to conventional multislice reconstruction algorithms. For a table speed up to 100 mm per rotation, GSR demonstrates good image quality for both the low-contrast ball phantom and thorax phantom. All other performance parameters are comparable to the single-slice 180 degrees LI (linear interpolation) algorithm, which is considered the "gold standard." GSR also achieves high computing efficiency and good temporal resolution, making it an attractive alternative for the reconstruction of next generation multislice spiral CT data.     

Fourier-based approach to interpolation in single-slice helical computed tomography

It has recently been shown that longitudinal aliasing can be a significant and detrimental presence in reconstructed single-slice helical computed tomography (CT) volumes. This aliasing arises because the directly measured data in helical CT are generally undersampled by a factor of at least 2 in the longitudinal direction and because the exploitation of the redundancy of fanbeam data acquired over 360 degrees to generate additional longitudinal samples does not automatically eliminate the aliasing. In this paper we demonstrate that for pitches near 1 or lower, the redundant fanbeam data, when used properly, can provide sufficient information to satisfy a generalized sampling theorem and thus to eliminate aliasing. We develop and evaluate a Fourier-based algorithm, called 180FT, that accomplishes this. As background we present a second Fourier-based approach, called 360FT, that makes use only of the directly measured data. Both Fourier-based approaches exploit the fast Fourier transform and the Fourier shift theorem to generate from the helical projection data a set of fanbeam sinograms corresponding to equispaced transverse slices. Slice-by-slice reconstruction is then performed by use of two-dimensional fanbeam algorithms. The proposed approaches are compared to their counterparts based on the use of linear interpolation-the 360LI and 180LI approaches. The aliasing suppression property of the 180FT approach is a clear advantage of the approach and represents a step toward the desirable goal of achieving uniform longitudinal resolution properties in reconstructed helical CT volumes.     

Evaluation of coccygeal bone variability, intercoccygeal and lumbo-sacral angles in asymptomatic patients in multislice computed tomography

The coccyx is a highly variable structure in the human caudal spine. Previous studies have revealed a significant correlation between coccyx shape and the pain syndrome coccygodynia. The aim of this study was to carry out a complex morphological evaluation of the coccyx in a group of asymptomatic patients of different sex and age examined by multislice computed tomography (MSCT) of the pelvis for different clinical reasons. MSCT pelvis examinations from various nontraumatic clinical conditions from consecutive adult patients (250 males and 250 females of comparable age, mean 54.9 ± 14.8 years) were used. Based on middle sagittal plane reconstructions: coccyx configuration (types I–IV according to Postacchini and Massobrio classification, each successive type characterized by a more pronounced anterior position of coccyx), number of segments, length and angles (intercoccygeal and lumbo-sacral) were measured. The results obtained were analyzed statistically. The following types of coccyx were observed in the study group: type I in 16.2 %, type II 40.0 %, type III 32.4 %, and type IV 11.4 % cases. In most cases (50.8 %), three segments were noted. Lumbo-sacral angle varied from 15.6° to 66.4° (average 41.6° ± 7.7°), and intercoccygeal angle from 0° to 107° (average 51° ± 23.3°). A significant negative correlation between age and number of segments as well as age and intercoccygeal angle was observed. In males, the coccyx was significantly longer, while in females the intercoccygeal angle was significantly wider. Type I was significantly more frequent in males, while type IV was found more often in females. The results obtained differ from other results in the literature. Our research could be useful to determine population standards, and help (together with clinical history) future studies of associations between idiopathic coccygodynia and coccyx morphology.     

Survey of computed tomography scanners in Taiwan: dose descriptors, dose guidance levels, and effective doses

The IAEA and the ICRP recommended dose guidance levels for the most frequent computed tomography (CT) examinations to promote strategies for the optimization of radiation dose to CT patients. A national survey, including on-site measurements and questionnaires, was conducted in Taiwan in order to establish dose guidance levels and evaluate effective doses for CT. The beam quality and output and the phantom doses were measured for nine representative CT scanners. Questionnaire forms were completed by respondents from facilities of 146 CT scanners out of 285 total scanners. Information on patient, procedure, scanner, and technique for the head and body examinations was provided. The weighted computed tomography dose index (CTDI(w)), the dose length product (DLP), organ doses and effective dose were calculated using measured data, questionnaire information and Monte Carlo simulation results. A cost-effective analysis was applied to derive the dose guidance levels on CTDI(w) and DLP for several CT examinations. The mean effective dose +/- standard deviation distributes from 1.6 +/- 0.9 mSv for the routine head examination to 13 +/- 11 mSv for the examination of liver, spleen, and pancreas. The surveyed results and the dose guidance levels were provided to the national authorities to develop quality control standards and protocols for CT examinations.     

In vivo dosimetry for estimation of effective doses in multislice CT coronary angiography

In vivo dosimetry represents a technique that has been widely employed to evaluate the dose to the patient mainly in radiotherapy. Considering the increment in dose to the population due to new high-dose multislice CT examinations, such as coronary angiography, it is becoming important to more accurately know the dose to the patient. The desire to know patient dose extends even to radiological examinations. Thermoluminescent dosimeters are considered the gold standard for in vivo dosimetry, but their use is time consuming. A rapid, less labor-intensive method has been developed to perform in vivo dosimetry using radiochromic film positioned next to the patient's skin. Multislice CT scanners allow the estimation of the effective dose to the patient from the dose length product (DLP) parameter, the value of which is displayed on the acquisition console, simply multiplying the DLP by published conversion factors. The method represents only an approximation based on standard size circular phantoms and neglects the actual size of the patient. More accurate evaluations can be carried out using software-based Monte Carlo simulations. However, these methods do not consider possible dose reduction techniques, such as automatic tube-current modulation. For 22 patients effective doses measured by in vivo dosimetry and calculated by software were compared. The technique of using in vivo dosimetry measured with radiochromic film appears a promising procedure for improving the assessment of the effective dose to the patient.     

Assessment of coronary artery bypass graft patency by multislice computed tomography

The recently developed multislice computed tomography (MSCT) is capable of rapid imaging of cardiac structures, including coronary arteries, during a single breath-hold. We evaluated coronary artery bypass graft (CABG) patency by comparing MSCT results to those of contrast angiography. MSCT and contrast angiography were performed in 39 patients (10 women, 29 men and mean age 60.0 +/- 7.8 years) with a total of 115 bypass grafts including 36 left internal mammary arteries, 4 right internal mammary arteries, 19 radial arteries, 2 gastroepiploic arteries and 54 vein grafts. Patients were investigated for an average of 14 +/- 27 months (range 1 - 108 months) after CABG surgery. Contrast angiography showed a patency rate of 87.0% (100/115). Ninety-nine of these 100 patent grafts by contrast angiography and 14 of the remaining 15 occluded grafts were correctly classified by MSCT (93.3% sensitivity and 99.0% specificity for bypass graft occlusion). The positive and negative predictive values for bypass graft occlusion were 93.3% and 99%, respectively, with an overall diagnostic accuracy of 98.3% (97.2% for left internal mammary artery, 100% for radial artery, 98.1% for vein graft and 100% for other grafts). In conclusion, MSCT is a useful and accurate diagnostic tool for the evaluation of bypass graft patency.     

The cause of the artifact in 4-slice helical computed tomography

The causes of the image artifacts in a 4-slice helical computed tomography have been discussed as follows: (1) changeover in pairs of data used in z interpolation, (2) sampling interval in z, and (3) the cone angle. This study analyzes the first two causes of the artifact and describes how the current algorithm [K. Taguchi and H. Aradate, Radiology 205P, 390 (1997); 205P, 618 (1997); Med. Phys. 25, 550-561 (1998); H. Hu, ibid. 26, 5-18 (1999); S. Schaller et al., IEEE Trans. Med. Imaging 19, 822-834 (2000); K. Taguchi, Ph.D. thesis, University of Tsukuba, 2002] solves the problem. An interpolated sinogram for a slice at the edge of a ball phantom shows discontinuity caused by the changeover. If we extend the streak artifact in the reconstructed image, it crosses the focus orbit at the corresponding projection angle. Applying z filtering can reduce such causes by its feathering effect and mixing data obtained by different cone angles; the best results are provided when z filtering is applied to densely sampled helical data.     

Usefulness of helical computed tomography in diagnosing pulmonary vein stenosis in infants

We investigated the usefulness of helical computed tomography(CT)in the morphological diagnosis of pulmonary vein stenosis, particularly that in infants and small children. In total, 20 helical CT examinations were performed in 10 post-operative cases of Total Anomalous Pulmonary Venous Drainage(TAPVD), 3 cases of single right ventricle, and 1 case of single left ventricle. In all cases, distinct morphological imaging was possible. Pulmonary vein stenosis could be categorized into three types: (1)stenosis from the anastomosis of the common pulmonary vein (CPV)-the left atrium (LA) to the peripheral pulmonary vein; (2) stenosis only at the anastomosis of CPV-LA; and (3) stenosis due to compression by nearby organs. Coronal views by multiplanar reconstruction (MPR) provided morphological information along the up-down direction of the body axis. Morphological diagnosis of pulmonary vein stenosis is important in deciding prognosis and therapeutic regimens, and helical CT was considered useful for such diagnosis in our 14 young patients.     

PI-line-based image reconstruction in helical cone-beam computed tomography with a variable pitch

Current applications of helical cone-beam computed tomography (CT) involve primarily a constant pitch where the translating speed of the table and the rotation speed of the source-detector remain constant. However, situations do exist where it may be more desirable to use a helical scan with a variable translating speed of the table, leading a variable pitch. One of such applications could arise in helical cone-beam CT fluoroscopy for the determination of vascular structures through real-time imaging of contrast bolus arrival. Most of the existing reconstruction algorithms have been developed only for helical cone-beam CT with constant pitch, including the backprojection-filtration (BPF) and filtered-backprojection (FBP) algorithms that we proposed previously. It is possible to generalize some of these algorithms to reconstruct images exactly for helical cone-beam CT with a variable pitch. In this work, we generalize our BPF and FBP algorithms to reconstruct images directly from data acquired in helical cone-beam CT with a variable pitch. We have also performed a preliminary numerical study to demonstrate and verify the generalization of the two algorithms. The results of the study confirm that our generalized BPF and FBP algorithms can yield exact reconstruction in helical cone-beam CT with a variable pitch. It should be pointed out that our generalized BPF algorithm is the only algorithm that is capable of reconstructing exactly region-of-interest image from data containing transverse truncations.     

Three-dimensional imaging and osteometry of adult human skulls using helical computed tomography

To make a digital image database of human craniology, we optimized the three-dimensional (3-D) images of 29 dried human skull specimens by helical computed tomography (CT). For the verification of the quantitative exactitude of these image data, we manually measured nine items of direct distances between standard anthropologic points on each skull and the corresponding distances projected on the CT monitor by specifying the respective points. The results obtained by the two methods of manual and CT measurements were compared and statistically analyzed. The CT measurements were so exact that the lower limit of correlation coefficients (95% of the confidence interval) between the two results was more than 0.8 in six items; i.e., maximal cranial length and breadth, minimal frontal breadth, bizygomatic breadth, distance between ectomolares and nasion-basion length. In contrast, the CT results were less well correlated with the manual measurements of three items; i.e., distance between bilateral mastoidales, total facial height, and nasal breadth. We concluded that the qualitative representation of 3-D CT images was adequate, although some quantitative data may be incorrect. The inaccuracy is suspected to be due to the difficulty in specifying the standard points on the CT images, and due to the differences in measurement procedures between the direct and projected distances.