Geometric models in dosimetry of thyroid remnant mass.

AIM: Absorbed dose to thyroid remnant tissue after 131I ablation becomes mass/size-dependent. This is a direct consequence of the small remnant size and radiation escape starts to be relevant. The self-absorbed fraction becomes mass/size-dependent. We have used Monte Carlo simulations to investigate the influence of the thyroid remnant shape upon the absorbed fraction calculation. METHODS: Thyroid residue was modeled using spherical, cylindrical and elliptical shapes. Uniform beta activity distribution and unit density medium (water) within a remnant was assumed. For each of the geometrical models beta self-absorbed fraction (varphi(beta)) was calculated using Monte Carlo codes, while the mean absorbed dose per unit cumulated activity (S(beta)) was calculated using MIRD formalism. RESULTS: For spherical objects varphi(mono) for mean beta energy (E = 0.182 MeV) of 131I is always greater than varphi(beta) calculated for the complete beta spectrum. For spheres having diameters 2-6 mm and assumption varphi(beta) = 1, S(beta) is over-estimated by 11-37%. For cylinder and prolate spheroid of the same length and thickness, S(beta) for cylinder is 30% smaller because of the greater mass. Similarly, elliptical cylinder and general ellipsoid of the same length and the same perpendicular dimensions (width and breadth), have similar varphi(beta), while S(beta) for elliptical cylinder is correspondingly smaller. CONCLUSION: For accurate dosimetry of thyroid remnants having masses <1 g and chordal lengths <1 cm it is necessary to calculate varphi(beta) for the full beta spectrum, or S(beta) will be overestimated. The shape of the remnant may also be important since elongated non-spherical objects may also have varphi(beta) < 1.

Using fast sequential asymmetric fanbeam transmission CT for attenuation correction of cardiac SPECT imaging.

UNLABELLED: The objective of this study was to determine the feasibility of using a fast (short-duration) transmission computed tomogram (TCT), acquired immediately before or after the emission CT, to correct for photon attenuation in cardiac SPECT. METHODS: The asymmetric fanbeam geometry with a 99mTc line source was used to acquire TCTs after conventional cardiac emission CT imaging on a triple-head SPECT system. The TCTs were reconstructed to generate patient-specific attenuation maps, which were used with an iterative maximum likelihood algorithm to reconstruct attenuation-corrected cardiac SPECT studies. The results of attenuation correction based on TCTs as short as 1 min were compared with long-duration transmission imaging for a phantom and several human studies. RESULTS: Attenuation correction based on asymmetric fanbeam TCT significantly improves the uniformity of images of a uniform tracer distribution in a cardiac-thorax phantom configured to simulate a large patient. By using a high-activity line source and a rapid camera rotation, a suitable attenuation map for this phantom can be obtained from a 4-min TCT. A similar result is obtained for patients with thorax widths of <40 cm. CONCLUSION: A sequential imaging protocol for acquiring a fast TCT can be used for attenuation correction of cardiac SPECT imaging. The sequential TCT can be acquired without significantly extending the duration of the imaging study. This method provides a way to perform attenuation correction on existing triple-head SPECT systems without extensively modifying the system.

Feasibility of dual radionuclide brain imaging with I-123 and Tc-99m.

A study was conducted to evaluate the feasibility of simultaneous dual radionuclide brain imaging with 123I and 99mTc using photopeak image subtraction techniques or offset photopeak image acquisition. The contribution of the photons from one radionuclide to a second radionuclide's photopeak energy window (crosstalk) was evaluated for SPECT and planar imaging of a brain phantom containing 123I and 99mTc for a range of activity levels and distribution properties approximating those in rCBF images of the adult human brain. Crosstalk was evaluated for 10% symmetrical energy windows centered on the 123I and 99mTc photopeaks and for 10% energy windows asymmetrically placed to the left and right of the center of the respective photopeaks. Major observations include: (1) in the centered photopeak windows, 99mTc crosstalk in the 123I window is 8.9% of the 99mTc seen in the 99mTc window and ranges from 37.5% to 75.0% of the 123I in the 123I window. 123I crosstalk is 37.8% of the 123I seen in the 123I window and ranges from 4.4% to 8.9% of the 99mTc seen in the 99mTc window; (2) the spatial distribution of a radionuclide's crosstalk photons differs from that observed in the radionuclide's photopeak window; (3) a 99mTc photopeak window offset to the left does not decrease 123I crosstalk, and the percentage of 99mTc scattered photons is significantly increased in the window.(ABSTRACT TRUNCATED AT 250 WORDS)

Asymmetric fan transmission CT on SPECT systems.

For proper attenuation correction of SPECT images, a set of 3D attenuation maps specific to the imaging slices is needed. Among the many different approaches for deriving the attenuation maps, fan beam transmission CT (FBTCT), performed on the same SPECT system as emission imaging, has many promising and clinically practical features. The major problem of FBTCT is that the current SPECT systems do not have a large enough field of view (FOV) to cover the typical cross-sectional size of patients. To address this problem, we have developed a novel asymmetric fan (AsF) sampling scheme to extend the FOV to practical sizes for clinical TCT imaging on existing SPECT systems. This AsF scheme samples only half of the intended FOV in each projection; the other half would be sampled in an opposing projection after detector rotation. We have implemented the AsF sampling on a three-head SPECT system through a specially designed source-collimator assembly. We have modified the conventional convolution backprojection algorithm to facilitate simple and fast image reconstruction. The feasibility of the approach is confirmed by the quality of the derived TCT images of various phantoms and human subjects. The AsF sampling scheme could also have applications in other general transmission CT systems.

Triple-head gamma camera PET: system overview and performance characteristics.

Positron emission tomography (PET) is currently performed using either a dedicated PET scanner or scintillation gamma camera equipped with electronic circuitry for coincidence detection of 511 keV annihilation quanta (gamma camera PET system). Although the resolution limits of these two instruments are comparable, the sensitivity and count rate performance of the gamma camera PET system are several times lower than that of the PET scanner. Most gamma camera PET systems are manufactured as dual-detector systems capable of performing dual-head coincidence imaging. One possible step towards the improvement of the sensitivity of the gamma camera PET system is to add another detector head. This work investigates the characteristics of one such triple-head gamma camera PET system capable of performing triple-head coincidence imaging. The following performance characteristics of the system were assessed: spatial resolution, sensitivity, count rate performance. The spatial resolution, expressed as the full width at half-maximum (FWHM), at 1 cm radius is 5.9 mm; at 10 cm radius, the transverse radial resolution is 5.3 mm, whilst the transverse tangential and axial resolutions are 8.9 mm and 13.3 mm, respectively. The sensitivity for a standard cylindrical phantom is 255 counts.s(-1).MBq*(-1)), using a 30% width photopeak energy window. An increase of 35% in the PET sensitivity is achievable by opening an additional 30% width energy window in the Compton region. The count rate in coincidence mode, at the upper limit of the systems optimal performance, is 45 kc.s(-1) (kc=kilocounts) using the photopeak energy window only, and increases to 60 kc.s(-1) using the photopeak + Compton windows. Sensitivity results are compared with published data for a similar dual-head detector system.

Determination of small objects from pinhole scintigrams.

This study assessed the possibility of measuring the linear dimensions of small structures using pinhole scintigraphy. A number of glass objects were made with a spherical, cylindrical or conical shape. Their maximum dimensions (diameters and heights) were 3.5-22.5 mm. These glass objects were filled with 131I, placed inside a plastic neck phantom and imaged using a gamma camera equipped with a pinhole collimator. The source-to-collimator distance was varied from 2 to 12 cm. An algorithm for image segmentation (threshold selection) was used to divide the image into object and background. On the segmented image, the number of non-zero pixels in the direction of the principal axes was multiplied by the appropriate calibration factor to obtain the linear dimensions of the object. Spatial resolution of the pinhole collimator, expressed as the full-width at half-maximum (FWHM), varied from 8 to 10 mm for the range of source-to-collimator distances examined. We found that, for dimensions up to 1.5 x FWHM, finite spatial resolution affects the accuracy of measurement. Non-linear correlation between true and calculated dimensions was used to take the latter into account. Our results are now being used to improve quantitation of remnant thyroid tissue masses for the calculation of radioiodine ablation doses.

3-D image analysis of intra-cerebral brain hemorrhage from digitized CT films.

A new 3-D technique for the segmentation and quantification of human spontaneous intra-cerebral brain hemorrhage (ICH) is presented in this paper. The algorithm for ICH primary region segmentation uses the spatially weighted K-means histogram-based clustering algorithm. The ICH edema region segmentation algorithm employs an iterative morphological processing of the ICH brain data. A volume rendering technique is used for the effective 3-D visualization of ICH segmented regions. A computer program is developed for use in the human spontaneous ICH study involving a large number of patients. Experimental measurements and visualization results are presented which were computed on real ICH patient brain data.

Web-based virtual endoscopy.

A web-based application for computed tomography (CT) image analysis and visualization of the bronchial airways is presented in the paper. The technique used in this application provides a noninvasive way to examine the interior of the bronchial tubes. Many different properties of the tubes such as presence of the foreign objects or structure abnormalities of the tubes can be detected this way. This method can also be used for other similar structures, such as blood vessels. The application uses the CT images obtained by spiral computed tomography as input. Amplitude thresholding is used for segmentation of the airways. The medial axis of the airways has been used as the path for animation. The procedure produces a fly-through animation in the MPEG format, which is obtained by volume rendering of the original CT images along the extracted path. The user can choose the appropriate path using the Virtual Reality Modeling Language (VRML) model, which is sufficient for visualization of basic structures. The algorithms and the software have been developed on the SUN Ultra 1 workstation.

3-D image analysis of abdominal aortic aneurysm.

In this paper we propose a technique for 3-D segmentation of abdominal aortic aneurysm (AAA) from computed tomography (CT) angiography images. Output data form the proposed method can be used for measurement of aortic shape and dimensions. Knowledge of aortic shape and size is very important for selection of appropriate stent graft device for treatment of AAA. The technique is based on a 3-D deformable model and utilizes the level-set algorithm for implementation of the method. The method performs 3-D segmentation of CT images and extracts a 3-D model of aortic wall. Once the 3-D model of aortic wall is available it is easy to perform all required measurements for appropriate stent graft selection. The method proposed in this paper uses the level-set algorithm instead of the classical active contour algorithm developed by Kass et al. The main advantage of the level set algorithm is that it enables easy segmentation surpassing most of the drawbacks of the classical approach. In the level-set approach for shape modeling, a 3-D surface is represented by a real 3-D function or equivalent 4-D surface. The 4-D surface is then evolved through an iterative process of solving the differential equation of surface motion. Surface motion is defined by velocity at each point. The velocity is a sum of constant velocity and curvature-dependent velocity. The stopping criterion is calculated based on image gradient. The algorithm has been implemented in MATLAB and C languages. Experiments have been performed using real patient CT angiography images and have shown good results.

Real-time 3D position reconstruction of guidewire for monoplane X-ray.

We present a novel real-time method for the 3D reconstruction of the guidewire using a monoplane X-ray. The method consists of two steps: (1) the backprojection step to reconstruct a 3D surface that contains the guidewire and (2) the optimization step to select a curve on the surface that is the best match under the pre-specified constraints. The proposed method utilizes a priori knowledge in the form of a volume that indicates positions of the blood vessels and thus restricts the reconstruction. The reconstruction precision is limited by the local thickness of the vessels. The method is quantitatively evaluated on five phantom datasets and qualitatively on two patient datasets. For the phantom datasets the average reconstruction error is resolution limited to 1-2 voxels and is biased in the depth direction. The worst-case reconstruction error for any point, including the guidewire tip, is not larger than the local vessel thickness. A visual inspection of results for the patient datasets shows the guidewire is always placed in the proper vessel and is aligned with the 2D image, which is sufficient for the guidewire navigation. The developed implementation achieves the processing speed of 12 fps using Core™i7 CPU 920 at 2.67 GHz.

Comparison of transaxial resolution in 180 degrees and 360 degrees SPECT with a rotating scintillation camera.

Using circular 180 degrees and 360 degrees SPECT acquisition modes the transaxial resolution of line sources in air and water were measured at different positions in the field of view. With the 180 degrees acquisition mode, all line sources in air located off the axis of camera rotation (AoR) showed an oval distortion. This distortion was systematically related to the starting point of the rotating detector. On axis line sources in air were undistorted, regardless of the 180 degrees acquisition starting angle. The 360 degrees acquisition images of the line sources in air showed a similar effect but in a very mild form. In water, transaxial reconstructions of line sources (off axis) showed an enhancement of the oval distortion for both the 180 degrees and 360 degrees acquisitions. Computer simulations of the line source measurements were performed and correlated well with the experimental data. The line source results are explainable by the inherent depth dependent response of the scintillation camera. In clinical SPECT studies, distortions of this nature will be most appreciable with 180 degrees imaging of small organs that are located off the AoR.

Short-axis circumferential profiles of the heart in healthy subjects: comparison of T1-201 SPECT and two-dimensional echocardiography.

To correlate the findings at thallium-201 single photon emission computed tomography (SPECT) with left ventricular anatomy, circumferential wall thickness was measured in 12 healthy subjects with two-dimensional echocardiography. At the midpapillary level, eight anatomic structures were identified with echocardiography, as were four maxima (papillary muscles and junctions of the left ventricular myocardium with the right ventricle, both anteriorly and posteriorly) and the intervening minima (including the septum). Tl-201 SPECT was performed in the same subjects. The normalized circumferential count profiles of the short-axis sections that included the papillary muscles showed the same basic pattern as that on echocardiograms. In most cases, the posterior papillary muscle and the posterior junction with the right ventricle were not distinguishable from each other, which produced three circumferential profile maxima instead of four. When values from all subjects were averaged, left ventricular anatomy was less evident in the circumferential profile. Left ventricular anatomic structures were reflected to a similar degree with either 180 degrees or 360 degrees data sampling.