Estimating Center of Rotation of Single-Photon Emission Computerized Tomography Projection Images Using MATLAB for Standardization and Calibration

Abstract

Objective: The objective of the study was to develop a Personal Computer (PC) based tool to estimate the center of rotation (COR) offsets from COR projection datasets using methods mentioned in IAEA-TECDOC-602. Materials and Methods: Twenty-four COR studies were acquired on Discovery NM 630 Dual head gamma camera fitted with parallel hole collimator, and COR offsets were estimated with the software available at the terminal for processing a COR study. These COR projection images were exported in DICOM. A MATLAB script (software program) was written to estimate COR offset using Method A (using opposite pair of projections) and Method B (using curve fitting method) as mentioned in IAEA-TECDOC-602. Our program read the COR study (in DICOM) and estimated COR offsets using Method A and Method B. The accuracy of the program was verified using simulated projection dataset of a point source object acquired at 6° interval in the range of 0°–360° angle. Bland Altman plot was used for analyzing the agreement between the COR offsets estimated using (1) Method A and Method B mentioned in IAEA-TECDOC-602, and (2) Our program and vendor program available at Discovery NM 630 acquisition terminal. Results: On simulated data, offset from center of gravity (COG) in X direction (COGX) and COG in Y direction (COGY) estimated using Method A was constant (same) at each pair of angles while using Method B, it was found to be in the range (−2 × 10 −10 , 1 × 10 −10) which is negligible. Most of the differences (23 out of 24) between the result of Method A and Method B, and between the results of our program and vendor program was found to be within 95% confidence interval (mean ± 1.96 standard deviation). Conclusions: Our PC-based tool to estimate COR offsets from COR projection datasets using methods mentioned in IAEA-TECDOC-602 was found to be accurate and provides results in agreement with vendor's program. It can be used as an independent tool to estimate COR offset for standardization and calibration purposes.