Temporal subtraction (TS) technique calculates a subtraction image between a pair of registered images acquired from the same patient at different times. Previous studies have shown that TS is effective for visualizing pathological changes over time; therefore, TS should be a useful tool for radiologists. However, artifacts caused by partial volume effects degrade the quality of thick-slice subtraction images, even with accurate image registration. Here, we propose a subtraction method for reducing artifacts in thick-slice images and discuss its implementation in high-speed processing. The proposed method is based on voxel matching, which reduces artifacts by considering gaps in discretized positions of two images in subtraction calculations. There are two different features between the proposed method and conventional voxel matching: (1) the size of a searching region to reduce artifacts is determined based on discretized position gaps between images and (2) the searching region is set on both images for symmetrical subtraction. The proposed method is implemented by adopting an accelerated subtraction calculation method that exploit the nature of liner interpolation for calculating the signal value at a point among discretized positions. We quantitatively evaluated the proposed method using synthetic data and qualitatively using clinical data interpreted by radiologists. The evaluation showed that the proposed method was superior to conventional methods. Moreover, the processing speed using the proposed method was almost unchanged from that of the conventional methods. The results indicate that the proposed method can improve the quality of subtraction images acquired from thick-slice images.
Keywords: Artifact reduction; Bone metastasis; CT; Temporal subtraction.