Meei-Ling Jan

A three-dimensional registration method for automated fusion of micro PET-CT-SPECT whole-body images

Micro positron emission tomography (PET) and micro single-photon emission computed tomography (SPECT), used for imaging small animals, have become essential tools in developing new pharmaceuticals and can be used, among other things, to test new therapeutic approaches in animal models of human disease, as well as to image gene expression. These imaging techniques can be used noninvasively in both detection and quantification. However, functional images provide little information on the structure of tissues and organs, which makes the localization of lesions difficult. Image fusion techniques can be exploited to map the functional images to structural images, such as X-ray computed tomography (CT), to support target identification and to facilitate the interpretation of PET or SPECT studies. Furthermore, the mapping of two functional images of SPECT and PET on a structural CT image can be beneficial for those in vivo studies that require two biological processes to be monitored simultaneously. This paper proposes an automated method for registering PET, CT, and SPECT images for small animals. A calibration phantom and a holder were used to determine the relationship among three-dimensional fields of view of various modalities. The holder was arranged in fixed positions on the couches of the scanners, and the spatial transformation matrix between the modalities was held unchanged. As long as objects were scanned together with the holder, the predetermined matrix could register the acquired tomograms from different modalities, independently of the imaged objects. In this work, the PET scan was performed by Concordes microPET R4 scanner, and the SPECT and CT data were obtained using the Gamma Medicas X-SPECT/CT system. Fusion studies on phantoms and animals have been successfully performed using this method. For microPET-CT fusion, the maximum registration errors were 0.21 mm /spl plusmn/ 0.14 mm, 0.26 mm /spl plusmn/ 0.14 mm, and 0.45 mm /spl plusmn/ 0.34 mm in the X (right-left), Y (upper lower), and Z (rostral-caudal) directions, respectively; for the microPET-SPECT fusion, they were 0.24 mm /spl plusmn/ 0.14 mm, 0.28 mm /spl plusmn/ 0.15 mm, and 0.54 mm /spl plusmn/ 0.35 mm in the X, Y, and Z directions, respectively. The results indicate that this simple method can be used in routine fusion studies.

A novel image quality index using Moran I statistics

Measurement of image quality is very important for various applications such as image compression, restoration and enhancement. Conventional methods (e.g., mean squared error; MSE) use error summation to measure quality change pixel by pixel and do not correlate well with subjective quality measurement. This is due to the fact that human eyes extract structural information from the viewing field. In this study a new quality index using a Moran I statistics is proposed. The Moran statistic that measures the sharpness from a local area is a good index of quality as most image processing techniques alter the smoothness of the image. Preliminary results show that the new quality index outperforms the MSE significantly under various types of image distortions.

Quality Degradation in Lossy Wavelet Image Compression

The objective of this study was to develop a method for measuring quality degradation in lossy wavelet image compression. Quality degradation is due to denoising and edge blurring effects that cause smoothness in the compressed image. The peak Moran z histogram ratio between the reconstructed and original images is used as an index for degradation after image compression. The Moran test is applied to images randomly selected from each medical modality, computerized tomography, magnetic resonance imaging, and computed radiography and compressed using the wavelet compression at various levels. The relationship between the quality degradation and compression ratio for each image modality agrees with previous reports that showed a preference for mildly compressed images. Preliminary results show that the peak Moran z histogram ratio can be used to quantify the quality degradation in lossy image compression. The potential for this method is applications for determining the optimal compression ratio (the maximized compression without seriously degrading image quality) of an image for teleradiology.

Improvements on a patient-specific dose estimation system in nuclear medicine examination

The purpose of this paper is to develop a patient-specific dose estimation system in nuclear medicine examination. A dose deposition routine to store the deposited energy of the photons during their flights was embedded in the widely used SimSET Monte Carlo code and a user-friendly interface for reading PET and CT images was developed. Dose calculated on ORNL phantom was used to validate the accuracy of this system. The ratios of S value for (99m)Tc, (18)F and (131)I computed by this system to those obtained with OLINDA for various organs were ranged from 0.93 to 1.18, which were comparable to that obtained from MCNPX2.6 code (0.88-1.22). Our system developed provides opportunity for tumor dose estimation which cannot be known from the MIRD. The radiation dose can provide useful information in the amount of radioisotopes to be administered in radioimmunotherapy.

Scatter correction for 3D PET using beam stoppers combined with dual-energy window acquisition: a feasibility study.

Fully three-dimensional (3D) positron emission tomography (PET) can achieve high sensitivity of coincidence events, but the absence of inter-slice septa inevitably leads to increased scattered events. The scattered events can represent as much as 50% of the total detected events. In this research, we proposed a scatter correction method for 3D PET based on beam stoppers and dual-energy window acquisition. The beam stoppers were placed surrounding the object to attenuate primary beams. The scatter fractions were directly estimated at those blocked lines of response and then the entire scatter fraction distribution was recovered using the dual-energy window ratio as reference. The performance was evaluated by using Monte Carlo simulations of various digital phantoms. For the Utah phantom study, the proposed method accurately estimated the scatter fraction distribution, and improved image contrast and quantification based on four different quality indices as performance measures. For the non-homogeneous Zubal phantom, the simulated results also demonstrated that the proposed method achieved a better restoration of image contrast than the dual-energy window method. We conclude that the proposed scatter correction method could effectively suppress various kinds of scattered events, including multiple scatter and scatter from outside the field of view.

Determination of beam intensity in a single step for IMRT inverse planning.

In intensity modulated radiotherapy (IMRT), targets are treated by multiple beams at different orientations each with spatially-modulated beam intensities. This approach spreads the normal tissue dose to a greater volume and produces a higher dose conformation to the target. In general, inverse planning is used for IMRT treatment planning. The inverse planning requires iterative calculation of dose distribution in order to optimize the intensity profile for each beam and is very computation intensive. In this paper, we propose a single-step method utilizing a figure of merit (FoM) to estimate the beam intensities for IMRT treatment planning. The FoM of a ray is defined as the ratio between the delivered tumour dose and normal tissue dose and is a good index for the dose efficacy of the ray. To maximize the beam utility, it is natural to irradiate the tumour with intensity of each ray proportional to the value of the FoM. The nonuniform beam intensity profiles are then fixed and the weights of the beam are determined iteratively in order to yield a uniform tumour dose. In this study, beams are employed at equispaced angles around the patient. Each beam with its field size that just covers the tumour is divided into a fixed number of beamlets. The FoM is calculated for each beamlet and this value is assigned to be the beam intensity. Various weighting factors are incorporated in the FoM computation to accommodate different clinical considerations. Two clinical datasets are used to test the feasibility of the algorithm. The resultant dose-volume histograms of this method are presented and compared to that of conformal therapy. Preliminary results indicate that this method reduces the critical organ doses at a small expense of uniformity in tumour dose distribution. This method estimates the beam intensity in one single step and the computation time is extremely fast and can be finished in less than one minute using a regular PC.

Noninvasive measurement of radiopharmaceutical time–activity data using external thermoluminescent dosimeters (TLDs)

In this study, we present a new method for estimating the time-activity data using serial timely measurements of thermoluminescent dosimeters (TLDs). The approach is based on the combination of the measurement of surface dose using TLD and Monte Carlo (MC) simulation to estimate the radiopharmaceutical time-activity data. It involves four steps: (1) identify the source organs and outline their contours in computed tomography images; (2) compute the S values on the body surface for each source organ using a MC code; (3) obtain a serial measurement of the dose with numerous TLDs placed on the body surface; (4) solve the dose-activity equation to generate organ cumulative activity for each period of measurement. The activity of each organ at the time of measurement is simply the cumulative activity divided by the timespan between measurements. The usefulness of this method was studied using a MC simulation based on an Oak Ridge National Laboratory mathematical phantom with 18F-FDG filled in six source organs. Numerous TLDs were placed on different locations of the surface and were repeatedly read and replaced. The time-activity curves (TACs) of all organs were successfully reconstructed. Experiments on a physical phantom were also performed. Preliminary results indicate that it is an effective, robust, and simple method for assessing the TAC. The proposed method holds great potential for a range of applications in areas such as targeted radionuclide therapy, pharmaceutical research, and patient-specific dose estimation.

Recovering the triple coincidence of non-pure positron emitters in preclinical PET*

Non-pure positron emitters, with their long half-lives, allow for the tracing of slow biochemical processes which cannot be adequately examined by the commonly used short-lived positron emitters. Most of these isotopes emit high-energy cascade gamma rays in addition to positron decay that can be detected and create a triple coincidence with annihilation photons. Triple coincidence is discarded in most scanners, however, the majority of the triple coincidence contains true photon pairs that can be recovered. In this study, we propose a strategy for recovering triple coincidence events to raise the sensitivity of PET imaging for non-pure positron emitters. To identify the true line of response (LOR) from a triple coincidence, a framework utilizing geometrical, energy and temporal information is proposed. The geometrical criterion is based on the assumption that the LOR with the largest radial offset among the three sub pairs of triple coincidences is least likely to be a true LOR. Then, a confidence time window is used to test the valid LOR among those within triple coincidence. Finally, a likelihood ratio discriminant rule based on the energy probability density distribution of cascade and annihilation gammas is established to identify the true LOR. An Inveon preclinical PET scanner was modeled with GATE (GEANT4 application for tomographic emission) Monte Carlo software. We evaluated the performance of the proposed method in terms of identification fraction, noise equivalent count rates (NECR), and image quality on various phantoms. With the inclusion of triple coincidence events using the proposed method, the NECR was found to increase from 11% to 26% and 19% to 29% for I-124 and Br-76, respectively, when 7.4-185 MBq of activity was used. Compared to the reconstructed images using double coincidence, this technique increased the SNR by 5.1-7.3% for I-124 and 9.3-10.3% for Br-76 within the activity range of 9.25-74 MBq, without compromising the spatial resolution or contrast. We conclude that the proposed method can improve the counting statistics of PET imaging for non-pure positron emitters and is ready to be implemented on current PET systems.

Scatter fraction performance tests for positron imaging system with dual plane geometry

NEMA NU 2 and IEC 61675–1 standards are widely used for characterizing the performance of PET scanners. These documents specify procedures for acquiring and analyzing data by using standard phantoms and radio-sources. However, some dedicated PET systems, such as breast and prostate PET systems, have non-ring geometric design may not meet these standards. In this work, we proposed certain modifications to the NEMA NU 2 for analyzing the scatter fraction (SF) of a dual-plane positron imaging system reasonably. A projection-map with geometric compensation was performed instead of sinogram to get analyzed profile. Three methods to determine analyzed window (4 × FWHM, 40 mm width, and 14 mm width) were used to estimate SF. The effects of diameter and location of line source were also considered. Three diameters (0.6, 2, 3.2 mm) and nine locations of line source were set for GATE simulation. The equivalent SF of the dual-plane system was calculated as the area-weighted average of nine-position SF values. The results reveal that the presented method here can estimate SF value better, compared with the methods suggested by the standards.

Application of the Intraoperative Dual Photon Emission Computed Tomography System in Sentinel Lymph Node Detection: A Simulation Study

The sentinel lymph node (SLN) hypothesis is applied as part of the standard procedure for identifying early-stage breast cancer. Thus, an imaging system that can locate SLNs in operating rooms is desired. Many 2-D probe imaging systems and a freehand single-photon emission-computed tomography (fhSPECT) system have been proposed. However, 2-D probe imaging systems are affected by shine-through and shadowing effects. Here, we propose an alternative to 3-D imaging systems, i.e., a dual-photon emission computed tomography (DuPECT) system, which integrates both preoperative and intraoperative information to locate SLNs using cascade isotopes such as Se-75. The system consists of a LaBr3-based probe and planar head, a collimation system, and a coincidence circuit. For each disintegration, the slat and parallel-hole collimator define a plane and a line, respectively, which represent the possible flight paths of each photon. SLNs can be located using the line-plane intersection. Here, the performance is evaluated using Monte Carlo software developed in our laboratory, integrated with SimSET and GATE software. A measurement study indicates that the randoms rate increases with increased initial activities, while the scatter rate is lower than 1.2 count/s for various activities. In a simulated imaging study, four injection sites and two LNs placed at various depths are minimally distinguishable. However, the LNs are clearly identifiable in the absence of injection sites. Our results indicate that the proposed three-dimensional imaging system has the potential to identify injection sites and various SLNs. However, difficulties with low sensitivity for LN detection, especially in the presence of activity from injection sites, and the choice of appropriate radioisotope must be overcome for its clinical usage.