The invention concerns the reproducible and accurate measurement of radio-tracer distribution in a body from the results of functional medical scanning procedures such as Positron Emission Tomography (PET) or Single Photon Emission Computed Tomography (SPECT).
Reconstruction is a step in a functional medical scanning procedure where an image is generated from the data acquired during the scan, the image depicting the localisation and concentration of radioisotope within the plane of the organ that was scanned.
For example, during a PET procedure, a subject is injected with a radioisotope such as [18F]-labeled 2-deoxyglucose (FDG) and, in time, the distribution of tracer throughout the body is measured by detecting the radiation produced when positrons emitted by the 18F radiolabel collide with electrons in neighboring atoms. This energy (in the form of gamma rays) is emitted as a pair of photons at approximately 180 degrees to each other and detection of such a pair allows a coincidence line to be deduced, along which the source of the radiation is assumed to lie.
Reconstruction algorithms use the information derived from a large number of such measurements at a large number of angular and linear positions to generate the image depicting the tracer concentration.
The reconstructed image provides a convenient tool for assessing the results of the scan and drawing inferences concerning metabolic reactions of the subject.
General purpose reconstruction techniques introduce algorithm specific artefacts. For example, Filtered Back Projection (FBP) (Kak and Slaney 1988) blurs the data uniformly, introducing strong correlations between neighbouring voxels.
Maximum A Posteriori (MAP) (Green 1990) and Maximum Likelihood Expectation Maximum (ML-EM) (Shepp and Vardi 1982) (Lange and Carson 1984) approaches both introduce non-uniform blurring and correlations that are spatially variant (Barrett, Wilson et al. 1994) and, in the case of dynamic data, temporally variant.
This causes sub-optimality in terms of task-specific operations performed by the user, including regional quantification, lesion detection, computer aided analysis and automated comparison against normal databases.
Thus far, work on addressing these problems has exploited the quantitative nature of the reconstructed PET/SPECT data, as illustrated generally with phantom experiments. Specific methods have been developed for the case of MAP reconstruction algorithms for achieving approximations of either uniform variance or uniform bias (Stayman 2002). In addition, effort has been directed towards the problem of how to control the smoothing by the use of informative priors onto the cost function (Comtat, Kinahan et al. 2002).
For ML-EM reconstruction, effort has been focussed on establishing the variance/bias at different points of convergence and finding optimal post-reconstruction smoothing parameters (Nuyts 2001). For FBP, effort has been directed towards identifying non-uniform blurring kernels to try and preserve edges.
All of these approaches involve developing new reconstruction algorithms, each of which introduces its own inaccuracies and variability, thus compounding the problems associated with inter-study comparisons.
Following reconstruction, post-analysis tasks are performed: these include lesion detection, region of interest (ROI) analysis and comparison with databases.
The invention reduces the problems associated with general purpose reconstruction by providing an approach that is optimized for the task under consideration which reduces inter-study variability.
A number of techniques are known that are applicable to reconstruction of specific regions of interest, see for example D. C. Schottlander, T. Kadir, J. M. Declerck, and M. Brady, “Unbiased quantification of tomographic data by projecting continuous regions of interest,” presented at the IEEE Nuclear Science Symposium and Medical Imaging Conference, San Diego, 2006; R. H. Huesman, “A new fast algorithm for the evaluation of regions of interest and statistical uncertainty in computed tomography,” Physics in Medicine and Biology, vol. 29, no. 5, pp. 543-52, 1984, UK Journal-Paper English 0031-9155. 92, 98 and R. E. Carson, “A maximum likelihood method for region-of-interest evaluation in emission tomography,” Journal of Computer Assisted Tomography, vol. 10, pp. 654-63, 1986.
According to the invention, a method of processing data acquired during the performance of a functional medical scan comprises the steps set out in claim 1 attached hereto.
The invention will now be described with reference to:
Referring to
In step 2, regions of interest in the reconstructed image are identified. This could be done using any of a number of methods including computer aided detection; registration with an atlas; visual inspection and qualities such as contrast or brightness.
In step 3, task specific reconstruction of the regions of interest identified at step 2 is performed. Possible approaches here include parametric map reconstruction using a model based approach such as compartment modelling; reconstruction to maximize lesion detectability; reconstruction of average image intensity values and associated confidence for each ROI and reconstruction of ROI homogeneity index.
The results of the task specific reconstruction and/or confidence values, can be displayed or presented as a report or graphic overlay.
Quantititative cardiac PET imaging requires accurate data measurements from the myocardium as a precursor to performing weighted model based fitting to a physiological model of myocardial blood flow. One way of achieving this is to reconstruct the data using an iterative algorithm, executed to, or close to, convergence. This implies running the algorithm for many hundreds of iterations, which is prohibitively expensive from a computational point of view. Task specific reconstruction provides a means of achieving this by using the following steps:
Number | Date | Country | Kind |
---|---|---|---|
0710428.4 | Jun 2007 | GB | national |