An automated technique for SPECT marker-based image registration in radionuclide therapy

SCreg: a registration-based platform to compare unicondylar knee arthroplasty SPECT/CT scans

Background

A combination of conventional computed tomography and single photon emitted computed tomography (SPECT/CT) provides simultaneous data on the intensity and location of osteoblastic activity. Currently, since SPECT/CT scans are not spatially aligned, scans following knee arthroplasty are compared by extracting average and maximal values of osteoblastic activity intensity from large subregions of the structure of interest, which leads to a loss of resolution, and hence, information. Therefore, this paper describes the SPECT/CT registration platform (SCreg) based on the principle of image registration to spatially align SPECT/CT scans following unicondylar knee arthroplasty (UKA) and allow full resolution intra-subject and inter-subject comparisons.

Methods

SPECT-CT scans of 20 patients were acquired before and 1 year after UKA. Firstly, scans were pre-processed to account for differences in voxel sizes and divided in volumes of interest. This was followed by optimization of registration parameters according to their volumetric agreement, and alignment using a combination of rigid, affine and non-rigid registration. Finally, radiotracer uptakes were normalized, and differences between pre-operative and post-operative activity were computed for each voxel. Wilcoxon signed rank sum test was performed to compare Dice similarity coefficients pre- and post-registration.

Results

Qualitative and quantitative validation of the platform assessing the correct alignment of SPECT/CT scans resulted in Dice similarity coefficient values over 80% and distances between predefined anatomical landmarks below the fixed threshold of (2;2;0) voxels. Locations of increased and decreased osteoblastic activity obtained during comparisons of osteoblastic activity before and after UKA were mainly consistent with literature.

Conclusions

Thus, a full resolution comparison performed on the platform could assist surgeons and engineers in optimizing surgical parameters in view of bone remodeling, thereby improving UKA survivorship.

Evaluation of sequential SPECT and CT for targeted radionuclide therapy dosimetry

PURPOSE

In targeted radionuclide therapy (TRT), a prior knowledge of the absorbed dose biodistribution is essential for pre-therapy treatment planning. Previously, we showed that non-rigid organ-by-organ registration in sequential quantitative SPECT images improved dose estimation. This study aims to investigate if sequential CT can further improve TRT dosimetric accuracy.

METHODS

We simulated SPECT/CT acquisitions at 1, 12, 24, 72 and 144 h In-111 Zevalin post-injection using an analytical MEGP projector, modeling attenuation, scatter and collimator-detector response. We later recruited a patient injected with 222 MBq In-111 DTPAOC imaged at 3 SPECT/CT sessions for clinical evaluations. Four registration schemes were evaluated: whole-body-based registration performed on sequential (1) SPECT (WB-SPECT) or (2) CT (WB-CT) images; organ-based registration applied on organs individually segmented from sequential (3) SPECT (O-SPECT) or (4) CT (O-CT) images. Voxel-by-voxel integration was performed followed by Y-90 voxel-S-kernel convolution. Organ-absorbed doses, iso-dose curves, dose-volume histograms (DVHs) were generated for targeted organs for analysis.

RESULTS

In simulation study, organ-absorbed dose errors were (- 8.66 ± 2.83)%, (- 2.51 ± 3.69)%, (- 9.23 ± 3.28)%, (- 7.17 ± 2.53)% for liver, (- 14.81 ± 4.91)%, (- 3.60 ± 4.37)%, (- 18.13 ± 4.44)%, (- 11.34 ± 4.22)% for spleen, for O-SPECT, O-CT, WB-SPECT and WB-CT registrations, respectively. For all organs, O-CT showed superior results. Results of iso-dose contour, DVHs were in accordance with the organ-absorbed doses. In clinical studies, the results were also consistent which showed O-CT method deviated the most from the result with no registration.

CONCLUSIONS

We conclude that if both sequential SPECT/CT scans are available, CT organ-based registration method can more effectively improve the 3D dose estimation. Sequential low-dose CT scans might be considered to be included in the standard TRT protocol.

MIRD pamphlet No. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy

In internal radionuclide therapy, a growing interest in voxel-level estimates of tissue-absorbed dose has been driven by the desire to report radiobiologic quantities that account for the biologic consequences of both spatial and temporal nonuniformities in these dose estimates. This report presents an overview of 3-dimensional SPECT methods and requirements for internal dosimetry at both regional and voxel levels. Combined SPECT/CT image-based methods are emphasized, because the CT-derived anatomic information allows one to address multiple technical factors that affect SPECT quantification while facilitating the patient-specific voxel-level dosimetry calculation itself. SPECT imaging and reconstruction techniques for quantification in radionuclide therapy are not necessarily the same as those designed to optimize diagnostic imaging quality. The current overview is intended as an introduction to an upcoming series of MIRD pamphlets with detailed radionuclide-specific recommendations intended to provide best-practice SPECT quantification-based guidance for radionuclide dosimetry.

The importance of the accuracy of image registration of SPECT images for 3D targeted radionuclide therapy dosimetry

In this paper, the importance of the accuracy of image registration of time-sequential SPECT images for 3D targeted radionuclide therapy dosimetry is studied. Image registration of a series of SPECT scans is required to allow the computation of the 3D absorbed dose distribution for both tumour sites and normal organs. Three simulated 4D datasets, based on patient therapy studies, were generated to allow the effect of mis-registration on the absorbed dose distribution to be investigated. The tumour sites studied range in size, shape and position, relative to the centre of the 3D SPECT scan. Randomly generated transformations along the x-, y- and z-axes and rotations around the z-axis were employed and the maximum and average absorbed dose distribution statistics, for the tumour sites present, were computed. It was shown that even small mis-registrations, translation of less than 9 mm and rotation of less than 5 degrees might cause differences in the absorbed dose statistics of up to 90%, especially when the size of the tumour is comparable to the induced mis-registration or when the tumour is situated close to the edge of the 3D dataset.

A generalized 4D image registration scheme for targeted radionuclide therapy dosimetry

An iterative, generalized four-dimensional (4D) method is presented in this paper that allows simultaneous registration of a series of single-photon emission computed tomography (SPECT) scans acquired in the course of a radionuclide therapy or pretherapy tracer study. The method combines temporal information with voxel-based similarity criteria to carry out simultaneous registration of the SPECT scans. A polynomial function was fitted to the maximum counts of each tumor site over the 4D study. Each tumor site was normalized to its maximum on the reference scan, and a template 4D dataset was generated, employing the polynomial fitting and the normalization map. Then, each 3D scan was registered to the corresponding simulated scan, using a 3D similarity criterion. The correlation coefficient (CC), the mutual information (MI), and the sum-of-absolute differences (SAD) similarity criteria were employed. Simulated data, based on a head-neck (131)I-MIBG study, were used to compare the proposed method for 4D registration with sequential 3D registration. Sequential 3D registration resulted in residual registration errors of 3.5 +/- 2.5, 3.2 +/- 2.0, and 7.0 +/- 3.5 mm for the CC, MI, and SAD criteria respectively, whereas the corresponding 4D method gave errors of 2.4 +/- 1.6, 1.9 +/- 1.1, and 5.3 +/- 2.9 mm for the CC, MI, and SAD criteria, respectively. The 4D method was applied to (186)Re HEDP SPECT patient studies and registration was verified by a dual-cursor display tool.

Quantitative analysis and effect of attenuation correction on lymph node staging of non-small cell lung cancer on SPECT and CT

OBJECTIVE

The purpose of our study was to assess quantitative indexes and the effect of attenuation correction on the evaluation of lymph node metastasis in the staging of non-small cell lung cancer (NSCLC) using fused thallium-201 SPECT/CT images.

MATERIALS AND METHODS

We evaluated 156 lymph nodes (66 metastatic, 90 nonmetastatic) from 29 patients with NSCLC. Using our combined SPECT/CT system, all patients underwent 201Tl SPECT and CT examinations immediately (early images) and 3 hr after (delayed images) the injection of 201Tl. SPECT images were reconstructed with and without attenuation correction. For the quantitative evaluation of lymph node metastasis, we calculated the early ratio, the delayed ratio, and the washout ratio for SPECT images and the short-axis diameter for CT images. Receiver operating characteristic (ROC) analysis was performed in each index for the differentiation between metastatic and nonmetastatic lymph nodes. Visual analysis was also performed by two experienced radiologists.

RESULTS

The area under the ROC curve (A(z)) showed that early ratio and delayed ratio were superior to short-axis diameter for the assessment of lymph node metastasis. In addition, early and delayed ratios on attenuation-corrected images were superior to those ratios on images without attenuation correction. However, the A(z) value for washout ratio was smaller than that for short-axis diameter. Early ratio on attenuation-corrected images was the most useful index (A(z) = 0.94). The sensitivity, specificity, and accuracy for early ratio on attenuation-corrected images were 78.8%, 94.4%, and 87.8% for the diagnosis of lymph node metastasis and 84.6%, 100%, and 93.1% for clinical staging (N0-N1 vs N2-N3), respectively. Fused images showed significantly higher diagnostic accuracy than CT images on visual analysis.

CONCLUSION

Quantitative assessment using fused SPECT/CT images is useful for the diagnosis of lymph node metastasis in patients with NSCLC.

Evaluation of a Rigid Registration Method of Lung Perfusion SPECT and Thoracic CT

OBJECTIVE

The objective of our study was to evaluate a rigid registration method in lung perfusion SPECT using thoracic CT as a standard.

MATERIALS AND METHODS

The reproducibility of markers selection and the robustness of the method were assessed on a torso phantom. The accuracy of registration regarding the number and location of markers and the breathing state during CT was evaluated on eight patients using 10 external markers placed around the thorax before SPECT and CT acquisitions. The accuracy of registration was assessed using the mean errors (ME) between 10 markers after registration.

RESULTS

Registration using external markers on a phantom was accurate (ME, < 3 mm) when rotation was less than 40 degrees (p = 0.02). The accuracy of registration improved markedly from four to six markers for phantom (5.5-3.6 mm) and patients (11.2-9.5 mm) and then remained constant up to 10 markers. The ME was less when using markers that well encompassed the thorax for phantom and patients (p = 0.02 and p = 0.05, respectively). The use of four anatomic markers was not accurate (ME, 20 mm).

CONCLUSION

The registration method is reproducible and accurate, and six external markers were required to get an ME of less than 10 mm in patients.

Prediction of radiation-induced normal tissue complications in radiotherapy using functional image data

The aim of this study has been to explicitly include the functional heterogeneity of an organ as a factor that contributes to the probability of complication of normal tissues following radiotherapy. Situations for which the inclusion of this information can be advantageous to the design of treatment plans are then investigated. A Java program has been implemented for this purpose. This makes use of a voxelated model of a patient, which is based on registered anatomical and functional data in order to enable functional voxel weighting. Using this model, the functional dose-volume histogram (fDVH) and the functional normal tissue complication probability (fNTCP) are then introduced as extensions to the conventional dose-volume histogram (DVH) and normal tissue complication probability (NTCP). In the presence of functional heterogeneity, these tools are physically more meaningful for plan evaluation than the traditional indices, as they incorporate additional information and are anticipated to show a better correlation with outcome. New parameters m(f), n(f) and TD(50f) are required to replace the m, n and TD(50) parameters. A range of plausible values was investigated, awaiting fitting of these new parameters to patient outcomes where functional data have been measured. As an example, the model is applied to two lung datasets utilizing accurately registered computed tomography (CT) and single photon emission computed tomography (SPECT) perfusion scans. Assuming a linear perfusion-function relationship, the biological index mean perfusion weighted lung dose (MPWLD) has been extracted from integration over outlined regions of interest. In agreement with the MPWLD ranking, the fNTCP predictions reveal that incorporation of functional imaging in radiotherapy treatment planning is most beneficial for organs with a large volume effect and large focal areas of dysfunction. There is, however, no additional advantage in cases presenting with homogeneous function. Although presented for lung radiotherapy, this model is general. It can also be applied to positron emission tomography (PET)-CT or functional magnetic resonance imaging (fMRI)-CT registered data and extended to the functional description of tumour control probability.

A novel four-dimensional image registration method for radionuclide therapy dosimetry

A novel method for registering sequential SPECT scans (4DRRT) is described, whereby all sequential scans acquired in the course of a therapy or a pre-therapy tracer study may be registered in one pass. The method assumes that a monoexponential decay function can be fitted to the series of sequential SPECT scans. Multiple volumes, presenting with different decay rates, are fitted with different mono-exponential functions. The MSSE (mean sum of squared errors in the least-squares fit algorithm), over the volume used for registration, is the cost function minimized at registration. Simulated data were used to assess the effect of thresholding, smoothing, noise and the multi-exponential nature of the four-dimensional (4D) SPECT studies on the performance of 4DRRT, resulting in three-dimensional (3D) residual registration errors <3.5 mm. The 4DRRT method was then compared to the following 3D registration methods: the correlation coefficient, the sum of absolute differences, the variance of image ratios and the mutual information. The comparisons, using both simulated and clinical data, were based on the standard deviation of the effective decay time distribution, generated from the registered 4D dataset, and showed that image registration using 4DRRT is simpler and more robust compared to the 3D techniques, especially when multiple tumour sites with different decay rates are present.

Stereotactic radiosurgery planning with ictal SPECT images

This paper is motivated by a clinical requirement to utilise ictal SPECT images for target localisation in stereotactic radiosurgery treatment planning using the xknife system which only supports CT and MRI images. To achieve this, the SPECT images were converted from raw (pixel data only) format into a part 10 compliant DICOM CT fileset. The minimum requirements for the recasting of a raw format image as DICOM CT or MRI data set are described in detail. The method can be applied to the importation of raw format images into any radiotherapy treatment planning system that supports CT or MRI import. It is demonstrated that the combination of the low spatial resolution SPECT images, depicting functional information, with high spatial resolution MRI images, which show the structural information, is suitable for stereotactic radiosurgery treatment planning.

Software for image registration: algorithms, accuracy, efficacy

Image registration is finding increased clinical use both in aiding diagnosis and guiding therapy. There are numerous algorithms for registration, which all involve maximizing a measure of similarity between a transformed floating image and a fixed reference image. The choice of the similarity measure depends, to some extent, on the application. Methods based on the use of the joint intensity histogram have become popular because of their flexibility and robustness. A distinction is made between rigid-body and non-rigid transformations. The latter are needed for inter-subject registration or intra-subject registration in cases where the region of the body of interest is not considered rigid. Non-rigid transformation is normally achieved using a global model of the deformation but can also be defined by a set of locally rigid transformations, each constrained to a small block in the image. There is scope for further research on the incorporation of appropriate constraints, especially for the application of non-rigid transformations to nuclear medicine studies. Most of the initial practical concerns regarding image registration have been overcome and there is increasing availability of commercial software. There are several approaches to the validation of registration software, with validation of non-rigid algorithms being particularly difficult. Studies have demonstrated the accuracy on the order of half a pixel for both intra- and inter-modality registration (typically 2 to 3 mm). Although hardware-based registration has now become possible by using dual-modality instruments, software-based registration will continue to play an important role in nuclear medicine.

RMDP: a dedicated package for 131I SPECT quantification, registration and patient-specific dosimetry

The limitations of traditional targeted radionuclide therapy (TRT) dosimetry can be overcome by using voxel-based techniques. All dosimetry techniques are reliant on a sequence of quantitative emission and transmission data. The use of (131)I, for example, with NaI or mIBG, presents additional quantification challenges beyond those encountered in low-energy NM diagnostic imaging, including dead-time correction and additional photon scatter and penetration in the camera head. The Royal Marsden Dosimetry Package (RMDP) offers a complete package for the accurate processing and analysis of raw emission and transmission patient data. Quantitative SPECT reconstruction is possible using either FBP or OS-EM algorithms. Manual, marker- or voxel-based registration can be used to register images from different modalities and the sequence of SPECT studies required for 3-D dosimetry calculations. The 3-D patient-specific dosimetry routines, using either a beta-kernel or voxel S-factor, are included. Phase-fitting each voxel's activity series enables more robust maps to be generated in the presence of imaging noise, such as is encountered during late, low-count scans or when there is significant redistribution within the VOI between scans. Error analysis can be applied to each generated dose-map. Patients receiving (131)I-mIBG, (131)I-NaI, and (186)Re-HEDP therapies have been analyzed using RMDP. A Monte-Carlo package, developed specifically to address the problems of (131)I quantification by including full photon interactions in a hexagonal-hole collimator and the gamma camera crystal, has been included in the dosimetry package. It is hoped that the addition of this code will lead to improved (131)I image quantification and will contribute towards more accurate 3-D dosimetry.

Automated CT marker segmentation for image registration in radionuclide therapy

In this paper a novel, automated CT marker segmentation technique for image registration is described. The technique, which is based on analysing each CT slice contour individually, treats the cross sections of the external markers as protrusions of the slice contour. Knowledge-based criteria, using the shape and dimensions of the markers, are defined to enable marker identification and segmentation. Following segmentation, the three-dimensional (3D) markers' centroids are localized using an intensity-weighted algorithm. Finally, image registration is performed using a least-squares fit algorithm. The technique was applied to both simulated and patient studies. The patients were undergoing 131I-mIBG radionuclide therapy with each study comprising several 99mTc single photon emission computed tomography (SPECT) scans and one CT marker scan. The mean residual 3D registration errors (+/- 1 SD) computed for the simulated and patient studies were 1.8 +/- 0.3 mm and 4.3 +/- 0.5 mm respectively.