This article has Open Peer Review reports available.
Standardized volumetric 3D-analysis of SPECT/CT imaging in orthopaedics: overcoming the limitations of qualitative 2D analysis
© Hirschmann et al; licensee BioMed Central Ltd. 2012
Received: 11 February 2011
Accepted: 29 March 2012
Published: 29 March 2012
SPECT/CT combines high resolution anatomical 3D computerized tomography (CT) and single photon emission computerized tomography (SPECT) as functional imaging, which provides 3D information about biological processes into a single imaging modality. The clinical utility of SPECT/CT imaging has been recognized in a variety of medical fields and most recently in orthopaedics; however, clinical adoption has been limited due to shortcomings of analytical tools available. Specifically, SPECT analyses are mainly qualitative due to variation in overall metabolic uptake among patients. Furthermore, most analyses are done in 2D, although rich 3D data are available. Consequently, it is difficult to quantitatively compare the position, size, and intensity of SPECT uptake regions among patients, and therefore difficult to draw meaningful clinical conclusions.
We propose a method for normalizing orthopaedic SPECT/CT data that enables standardised 3D volumetric quantitative measurements and comparison among patients. Our method is based on 3D localisation using clinically relevant anatomical landmarks and frames of reference, along with intensity value normalisation using clinically relevant reference regions. Using the normalised data, we describe a thresholding technique to distinguish clinically relevant hot spots from background activity.
Using an exemplar comparison of two patients, we demonstrate how the normalised, 3D-rendered data can provide a richer source of clinical information and allow quantitative comparison of SPECT/CT measurements across patients. Specifically, we demonstrate how non-normalized SPECT/CT analysis can lead to different clinical conclusions than the normalized SPECT/CT analysis, and that normalized quantitative analysis can be a more accurate indicator of pathology.
Conventional orthopaedic frames of reference, 3D volumetric data analysis and thresholding are used to distinguish clinically relevant hot spots from background activity. Our goal is to facilitate a standardised approach to quantitative data collection and comparison of clinical studies using SPECT/CT, enabling more widespread clinical use of this powerful imaging tool.
SPECT/CT combines a high resolution anatomical 3D computerized tomography (CT) and a single photon emission computerized tomography (SPECT) as functional imaging, which provides 3D information about biological processes into a single imaging modality. Along with the fusion of these images the metabolic activity in a region of interest can be accurately allocated to specific anatomical areas [1–3]. Until recently, registration of the functional images (SPECT) with the structural images (CT) has been difficult due to SPECT's poor spatial resolution (3-10 mm) and its high variability in identifying anatomical landmarks [4, 5]. With the advent of integrated SPECT/CT technology, where the SPECT data is obtained in the same coordinate frame as the CT, image fusion has significantly improved .
Although SPECT/CT has proven useful as a diagnostic tool in a variety of medical fields and most recently in orthopaedics [1, 7–11], its use has a number of limitations. Firstly, SPECT tracer uptake values are only valid for the individual patient investigated and cannot easily be compared between patients. Therefore, it is difficult to draw clinically useful conclusions on intensity and size of SPECT regions. Several authors have introduced methods of quantification of SPECT data to minimize subjective tracer interpretation, however to date no quantitative analysis has been described in orthopaedic patients [12–16]. Secondly, although 3D-data are generally obtained, analysis is routinely performed only in 2D (ROI analysis) . Thirdly, accurate and reproducible localization of increased or decreased tracer uptake in 3D is difficult to achieve [10, 11]. Fourthly, to date clinical correlation of the biology, the biomechanics and the function of joint pathology has not been sufficiently studied in orthopaedics.
Enhanced tracer uptake regions
The pattern and intensity (increased, decreased) of tracer distribution in a SPECT scan gives information on bone and soft tissue pathologies. Radionuclide tracers are chosen that preferentially bind to a specific metabolic activity (e.g., technetium-labeled diphosphonate for osseous activity), and injected intravenously 3-4 hours before the SPECT and CT imaging. Regions of increased tracer uptake (which we term a "hotspot") may signal a metabolically active region and provide a diagnostic aid. Current analysis, however, is limited to this relative comparison (increased or decreased) as opposed to quantifying the intensity of the tracer uptake region in a clinically significant way.
A key hurdle is the significant inter- and intra-patient variation in the metabolic uptake of radionuclide tracers. Variation of more than 10 times in the average SPECT intensity is not uncommon between patients. Because of this variation, it is difficult to identify a clinically relevant SPECT intensity threshold to distinguish regions of interests from normal background variation.
The purpose of this paper is to propose a method to normalize the 1) location, 2) size, and 3) intensity of altered tracer uptake regions to allow quantitative comparison of orthopaedic SPECT/CT scans. Our method uses standardized orthopaedic frames of reference for orientation, and 3D volumetric data interpretation and thresholding to distinguish clinically relevant hot spots from background activity. Along with its introduction we strive to facilitate data collection and comparison of clinical studies using SPECT/CT. In this section we present the normalization method details. We then provide an exemplar patient comparison as a demonstration of the method's use in the Results.
The proposed normalization techniques and analysis can be carried out with any software that can reorient 3D data and quantify intensities of the SPECT voxels. We have implemented the software tools necessary for the exemplar patient comparison using a collection of open-source tools. The software was written in Python (v2.6, http://python.org/), using the Grassroots DICOM library (GDCM v2.0.12, http://gdcm.sourceforge.net/) to read the SPECT and CT DICOM files exported by the SPECT/CT machine, the Visualisation Toolkit (VTK v5.4.2, Kitware, http://vtk.org/) library for surface and volume rendering, and the NumPy (http://numpy.scipy.org/) Python library for data analysis. Note that a slight modification was required to the GDCM library to overcome an idiosyncrasy in the SPECT DICOM data for the specific SPECT/CT machine used (Symbia T16, Siemens).
Normalisation of position: orthopaedic reference frames
To reliably locate and describe the position of a region of altered metabolic activity it is essential that the reference system is well defined. Further, to overcome the relative position and orientation of the patient to the coordinate system of the scanner, the reference system should be based on anatomic landmarks that are identifiable in subsequent scans. The unique combination of anatomic and metabolic information available with SPECT/CT imaging allows the low resolution metabolic information from SPECT to be accurately localised using the high resolution anatomic information from CT.
We use well-recognised standard anatomic landmarks and corresponding frames of reference that relate to the biomechanics of the joint. The landmarks can be manually identified on a combination of CT slices and on the surface of a 3D bone reconstruction from the CT data. These landmarks define a coordinate system specific to the lower limb, independent of scanner-specific coordinates. We can then use this normalised coordinate system to localize the biomechanical and diagnostic SPECT regions of interest. Further, with this normalised coordinate system we can provide normalised views for the 3D bone reconstruction and 2D slices (i.e., "true" antero-posterior views and slices that will be the same perspective relative to the anatomy independent of the patient position and orientation in the scanner).
Normalisation of size and intensity: clinically neutral reference region
Here we propose a method for normalising the SPECT intensities between patients to allow for quantitative analysis, potentially increasing the value of SPECT imaging. Direct comparison of raw SPECT intensities is known to be problematic, due to significant variation among patients in overall uptake levels. Normalisation techniques that rely purely on the distribution of SPECT intensities, such as scaling by the maximum value, or normalising using the mean and standard deviation, are also problematic. These techniques can be biased by the fraction of tissue related voxels as compared to free air voxels (affected by patient size and field of view), the presence of an intense hotspot, or an overall inflammatory uptake response.
What we propose is to use the anatomic information provided by the CT component to establish, prior to analysis, a clinically relevant neutral region. For example, when analysing intensities relating to knee joint pain, we can define a volume in the middle of the femoral shaft (which is unrelated to the articular surface) as a reference region. The distribution of SPECT intensities from this region can be then be used to normalize other values. This approach avoids the biases associated with patient size, field of view, and normalising solely on the global distribution of SPECT intensities. It also takes advantage of the unique combination of anatomic and metabolic information offered by SPECT/CT imaging.
The optimal method to normalize hotspot intensities using reference region intensities, resulting in an analysis with the best correlation with clinical outcomes, remains an open research question and likely depends on the specific analysis conducted. Potential approaches include scaling intensities to the mean of the reference region, offsetting intensities with the mean of the reference region, using the standard deviation of the values in the reference region as a scaling parameter, or some combination of these techniques. Hotspots can then be compared in a statistical manner, using a t-test to compare the distribution of values in a hotspot with the values in the reference region. Intensity thresholds can also be defined in a straightforward way, such as more than n standard deviations away from the reference region mean.
Thresholding for volumetric 3D tracer uptake analysis and visualisation
Here we present an exemplar analysis of SPECT/CT scans of two patients presented to our hospital, demonstrating the establishment of orthopaedic frames of reference, the localisation of SPECT tracer uptake regions, the normalisation of the intensities (based on the intensities in a neutral region in the femoral shaft) and a quantitative comparison of the size and intensity of the SPECT voxels in patients with symptomatic patellofemoral joints.
Both patients presented to our hospital complaining about knee pain with the history and clinical examination suggestive of patellofemoral disease. We decided to demonstrate the proposed method on these two patients because most imaging modalities are limited when it comes to the patellofemoral joint. Both patients have been diagnosed with patellofemoral osteoarthritis, but show a different distribution of osteoarthritis. In one the cause is considerably a slight trochlear dysplasia and the other one is idiopathic.
For both patients, Tc-99m-HDP-SPECT/CT was performed using a hybrid system (Symbia T16, Siemens, Erlangen, Germany) that consists of a pair of low energy, high-resolution collimators and a dual-head gamma camera and an integrated 16x0.75-mm slice-thickness CT. A commercial 700 MBq Tc-99m HDP tracer (CIS Bio International Sur Yvette, France) was used. SPECT imaging was performed at 3-4 hours post injection, followed by a CT without moving the patient. In addition, 2D perfusion images were obtained in the initial 60 seconds post injection followed by 2D blood pool images 2-5 minutes later.
Non-normalized analysis and visualisation
Normalisation of position: orthopaedic reference frames
Normalisation of size and intensity: clinically neutral reference region
We have proposed a novel method of standardized quantitative volumetric 3D-analysis of SPECT/CT imaging to overcome the limitations of qualitative 2D analysis. A standardised approach to quantitative volumetric 3D data collection using SPECT/CT enabling analysis of the size, the intensity and the localisation of enhanced tracer uptake regions of SPECT/CT, that we put forward as "hotspot", was introduced. The major contributions of our method are the following:
1. Standardised 3D-CT reconstructions using relevant surgical anatomical landmarks are used to define orthopaedic frames of reference for biomechanical orientation . Only after orientation of the data in relation to the mechanical axis of the lower limb can reliable and meaningful information be drawn for the diagnosis and treatment of patients with musculoskeletal disorders [11, 19–21]. To our knowledge the proposed method offers for the first time the chance to relate the information of joint mechanics (CT) with clinical hotspots as biologic/metabolic data (SPECT).
Current imaging and image processing techniques are not sufficiently advanced to investigate certain pathologies and compartments of joints. The patellofemoral disease is one such condition, which is difficult to differentiate from the other two compartments of the knee. In addition, it is a typical example for the symbiosis of joint mechanics and biology. This is also relevant for patients following joint replacement surgery [11, 18, 20, 21]. This 3D localisation method enables the surgeon to more clearly identify the clinically relevant hotspots and aid in understanding of the origin of the patient's symptoms.
2. It is widely known that individual SPECT hotspots vary greatly among patients. The normalizing techniques described here allow comparison of quantitative measurements amongst different patients. It may contribute to the understanding of pathologies and clinical value in clinical trials using SPECT/CT [11, 20, 22, 23]. Using our proposed method of image analysis and processing we hope to distinguish clinically relevant hotspots from normal background activity. Normalization approaches are biased by the patient size, the field of view and metabolic activity in areas outside the region of clinical interest . The establishment of clinically neutral reference regions prior to analysis addresses these problems. The limitations of the conventional approach to SPECT data normalisation, which relies on the distribution of SPECT intensities, scaling by the maximum values or normalisation using the mean or standard deviation are also avoided. However, the optimal method to normalize hotspot intensities using reference region intensities is an open research question. It may depend on different methods of data analysis [11, 20, 22, 23]. The optimal approach to scaling of intensities has to be investigated in clinical studies. We suggest studies looking into the normalisation of hotspots and offsetting intensities with the mean of the reference region. Using the standard deviation of the values in the reference region as a scaling parameter or some combination of these techniques might also be helpful to improve diagnostics further . Hotspots can then be compared in a statistical manner to compare the distribution of values within a hotspot with that of the reference region. Intensity thresholds can also be defined as a factor of the standard deviation from the reference region mean.
This proposed method of quantitative 3D analysis is novel and contrasts previous clinical studies using SPECT or SPECT/CT in orthopaedics, where the SPECT intensity and size of hotspots was graded using a Likert scale (0-10) or descriptively analysed [11, 20–23]. To date, most authors consider the area of maximum tracer activity as the area of interest analysing the peak values only [11, 21]. However, this method neglects the lower intensity SPECT values and may lead to a decrease in sensitivity and specificity in specific clinical conditions. It may also be case that the pattern of uptake provides the surgeon with more clinically relevant information leading to a more accurate diagnosis.
We propose a standardized method of analysing the location, size and the intensity of SPECT/CT tracer uptake regions ("hotspots"). Conventional orthopaedic frames of reference, 3D volumetric data analysis and thresholding are used distinguish clinically relevant hot spots from background activity. We demonstrated how the normalised, 3D-rendered data may provide a richer source of clinical information and allow quantitative comparison of SPECT/CT measurements across patients. This article gives a detailed description on the proposed method, but its clinical utility will be investigated in further clinical studies.
We would like to thank the Deutsche Arthrose-Hilfe e.V., Germany and the University of Basel for supporting our research activities with a project grant.
- Bybel B, Brunken RC, DiFilippo FP, Neumann DR, Wu G, Cerqueira MD: SPECT/CT imaging: clinical utility of an emerging technology. Radiographics. 2008, 28: 1097-1113.View ArticlePubMedGoogle Scholar
- Patton JA, Turkington TG: SPECT/CT physical principles and attenuation correction. J Nucl Med Technol. 2008, 36: 1-10.View ArticlePubMedGoogle Scholar
- Schillaci O: Hybrid SPECT/CT: a new era for SPECT imaging?. Eur J Nucl Med Mol Imaging. 2005, 32: 521-524.View ArticlePubMedGoogle Scholar
- O'Connor MK, Kemp BJ: Single-photon emission computed tomography/computed tomography: basic instrumentation and innovations. Semin Nucl Med. 2006, 36: 258-266.View ArticlePubMedGoogle Scholar
- Bernstein MP, Caldwell CB, Antonyshyn OM, Cooper PW, Ehrlich LE: Spatial and temporal registration of CT and SPECT images: development and validation of a technique for in vivo three-dimensional semiquantitative analysis of bone. J Nucl Med. 2000, 41: 1075-1081.PubMedGoogle Scholar
- Bocher M, Balan A, Krausz Y, Shrem Y, Lonn A, Wilk M, Chisin R: Gamma camera-mounted anatomical X-ray tomography: technology, system characteristics and first images. Eur J Nucl Med. 2000, 27: 619-627.View ArticlePubMedGoogle Scholar
- Agriantonis DJ, Hall L, Wilson MA: Utility of SPECT/CT as an adjunct to planar whole body I-131 imaging: liver metastasis from papillary thyroid cancer. Clin Nucl Med. 2009, 34: 247-248.View ArticlePubMedGoogle Scholar
- Akincioglu C, Belhocine T, Gambhir S, Romsa J, Chouraiki G, Kribs S, Elliot JA, Vezina W, Urbain JL: Complementary roles of low-dose SPECT-CT and high-resolution volume CT for detection of coronary artery disease. Clin Nucl Med. 2008, 33: 285-287.View ArticlePubMedGoogle Scholar
- Even-Sapir E, Keidar Z, Bar-Shalom R: Hybrid Imaging (SPECT/CT and PET/CT)-improving the diagnostic accuracy of functional/metabolic and anatomic imaging. Semin Nucl Med. 2009, 39: 264-275.View ArticlePubMedGoogle Scholar
- Hirschmann MT, Iranpour F, Davda K, Rasch H, Hugli R, Friederich NF: Combined single-photon emission computerized tomography and conventional computerized tomography (SPECT/CT): clinical value for the knee surgeons?. Knee Surg Sports Traumatol Arthrosc. 2009, 18: 341-345.View ArticlePubMedGoogle Scholar
- Hirschmann MT, Iranpour F, Konala P, Kerner A, Rasch H, Cobb JP, Friederich NF: A novel standardized algorithm for evaluating patients with painful total knee arthroplasty using combined single photon emission tomography and conventional computerized tomography. Knee Surg Sports Traumatol Arthrosc. 2010, 18: 939-44.View ArticlePubMedGoogle Scholar
- Blake GM, Frost ML, Fogelman I: Quantitative radionuclide studies of bone. J Nucl Med. 2009, 50: 1747-1750.View ArticlePubMedGoogle Scholar
- Blake GM, Moore AE, Fogelman I: Quantitative studies of bone using (99m)Tc-methylene diphosphonate skeletal plasma clearance. Semin Nucl Med. 2009, 39: 369-379.View ArticlePubMedGoogle Scholar
- Filippi L, Manni C, Pierantozzi M, Brusa L, Danieli R, Stanzione P, Schillaci O: 123I-FP-CIT semi-quantitative SPECT detects preclinical bilateral dopaminergic deficit in early Parkinson's disease with unilateral symptoms. Nucl Med Commun. 2005, 26: 421-426.View ArticlePubMedGoogle Scholar
- Sacchetti GM, Ghisellini F, Brambilla M, De Consoli A, Fornara P, Rizzo E, Rudoni M, Inglese E, Cannas M: Quantitative scintigraphic evaluation of total knee arthroplasties: a feasibility study. Clin Orthop Relat Res. 1996, 325: 181-189.View ArticlePubMedGoogle Scholar
- Shiraishi S, Tomiguchi S, Utsunomiya D, Kawanaka K, Awai K, Morishita S, Okuda T, Yokotsuka K, Yamashita Y: Quantitative analysis and effect of attenuation correction on lymph node staging of non-small cell lung cancer on SPECT and CT. AJR Am J Roentgenol. 2006, 186: 1450-1457.View ArticlePubMedGoogle Scholar
- Pennock GR, Clark KJ: An anatomy-based coordinate system for the description of the kinematic displacements in the human knee. J Biomech. 1990, 23: 1209-1218.View ArticlePubMedGoogle Scholar
- Hirschmann MT, Konala P, Amsler F, Iranpour F, Friederich NF, Cobb JP: Comparison of inter- and intra-observer reliability of component position measurements in total knee arthroplasty on conventional radiographs, transverse 2D-CT slices and 3D-CT. J Bone Joint Surg Br. 2010, 93: 629-633.Google Scholar
- Henckel J, Richards R, Lozhkin K, Harris S, Baena FM, Barrett AR, Cobb JP: Very low-dose computed tomography for planning and outcome measurement in knee replacement. The imperial knee protocol. J Bone Joint Surg Br. 2006, 88: 1513-1518.View ArticlePubMedGoogle Scholar
- Hirschmann MT, Davda K, Iranpour F, Rasch H, Friederich NF: Combined single photon emission computerised tomography and conventional computerised tomography (SPECT/CT) in patellofemoral disorders: a clinical review. Int Orthop. 2011, 35: 675-680.View ArticlePubMedGoogle Scholar
- Hirschmann MT, Konala P, Iranpour F, Kerner A, Rasch H, Friederich NF: Clinical value of SPECT-CT for evaluation of patients with painful knees after total knee arthroplasty - a new dimension of diagnostics. BMC Musculoskelet Disord. 2011, 12: 36-View ArticlePubMedPubMed CentralGoogle Scholar
- Hirschmann MT, Adler T, Rasch H, Hugli RW, Friederich NF, Arnold MP: Painful knee joint after ACL reconstruction using biodegradable interference screws- SPECT/CT a valuable diagnostic tool? A case report. Sports Med Arthrosc Rehabil Ther Technol. 2: 24-Google Scholar
- Konala P, Iranpour F, Kerner A, Rasch H, Friederich NF, Hirschmann MT: Clinical benefit of SPECT/CT for follow-up of surgical treatment of osteochondritis dissecans. Ann Nucl Med. 24: 621-624.Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/12/5/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.