Medicine:Transconvolution

From HandWiki

The term transconvolution designates a numerical method used in medical imaging, in particular emission computed tomography. Transconvolution enables a subsequent manipulation of the Point spread function (PSF) in already recorded images.

Properties of an image such as the spatial resolution or the appearance of small objects are determined by the PSF of the imaging system used for image acquisition. Different imaging systems with different PSFs therefore provide slightly different images of one and the same object.

Starting from known PSFs of different tomographic systems, the transconvolution method allows an image recorded on a particular tomograph to be converted as if it had been acquired by another tomograph. The method can thus ensure the comparability of images that were originally recorded on different systems.

Definition

Given two different tomographs with different point spread functions [math]\displaystyle{ psf_1 }[/math] and [math]\displaystyle{ psf_2 }[/math] the imaging process can be defined in terms of convolution as

[math]\displaystyle{ obj * psf_1 = img_1 }[/math]
[math]\displaystyle{ obj * psf_2 = img_2 }[/math]

with "[math]\displaystyle{ * }[/math]" representing the convolution operator and [math]\displaystyle{ img_1 }[/math] and [math]\displaystyle{ img_2 }[/math] representing the two slightly different images of the same object [math]\displaystyle{ obj }[/math] as seen by the respective tomographs.

The two equations yield the relationship

[math]\displaystyle{ img_1 * psf_1^{-1} * psf_2 = img_2 }[/math]

with [math]\displaystyle{ psf_1^{-1} }[/math] representing the inverse function of the according point spread function [math]\displaystyle{ psf_1 }[/math].

The inverse point spreading function [math]\displaystyle{ psf_1^{-1} }[/math] diverges and can not be determined or handled numerically. But, within certain boundary conditions, the complete term [math]\displaystyle{ psf_1^{-1} * psf_2 }[/math] is approximately computable by numerical methods.

The transconvolution function [math]\displaystyle{ tf }[/math] is defined as

[math]\displaystyle{ tf = psf_1^{-1} * psf_2 }[/math]

which results in the formula

[math]\displaystyle{ img_1 * tf = img_2 }[/math]

With the PSFs of the respective tomographs known, it is thus possible to convert an image [math]\displaystyle{ img_1 }[/math] recorded by the first tomograph into an [math]\displaystyle{ img_2 }[/math] emulating an image as recorded by the second tomograph. Of course, the method is subject to certain limits, in particular the computed [math]\displaystyle{ img_2 }[/math] can not represent spatial frequencies not captured to at least some degree by [math]\displaystyle{ psf_1 }[/math]. Consequently, the spatial resolution of an image can not be increased arbitrarily.

Application in medical imaging

The second point spread function [math]\displaystyle{ psf_2 }[/math] does not have to represent a real tomograph, but can be purposely defined to represent a virtual tomograph with corresponding properties. Based on the definition of a standardized virtual tomograph and the determination of the imaging properties of different real tomographs, the transconvolution method allows a uniform and quantitatively comparable representation of the image data taken by the different tomographs or systems, as if all measurements were made consistently by the standardized virtual system. The method thus supports quantitative comparisons of images taken by different imaging systems and in particular by different clinical tomographs.

Another application of the transconvolution method in positron emission tomography allows to handle the diverse image blur caused by the according diverse positron range of the actual positron emitting radionuclide. In particular this allows to use differing radionuclides for calibration as opposed to the subsequent imaging.

Literature

  • T. Weitzel, F. Corminboeuf, B. Klaeser, T Krause: "Kreuzkalibrierung von Positronen-Emissions-Tomographen für multizentrische Studien: Festkörper-Phantom und Transconvolution." In: Bulletin der Schweizerische Gesellschaft für Strahlenbiologie und Medizinische Physik. Oct. 2010, S. 9ff. SGSMP-Bulletin 72
  • T. Weitzel, F. Corminboeuf, B. Klaeser, T. Krause, T Beyer: "Transconvolution and virtual PET: A new concept for quantification of PET in multi-center trials." In: Journal of Nuclear Medicine. 51, 2010, S. 115.
  • G. A. Prenosil, T. Weitzel, M. Hentschel, B. Klaeser, T. Krause: "Transconvolution and the virtual positron emission tomograph–a new method for cross calibration in quantitative PET/CT imaging." In: Medical physics. Vol 40, No. 6, June 2013, S. 062503, doi:10.1118/1.4805112, PMID 23718608.
  • G. A. Prenosil, M. Hentschel, M. Fürstner, T. Krause, T. Weitzel, B. Klaeser: "Technical Note: Transconvolution based equalization of positron energy effects for the use of 68Ge/68Ga phantoms in determining 18F PET recovery." In: Medical Physics., 44, 2017 S. 3761

External links