Development of stopping rule methods for the MLEM and OSEM algorithms used in PET image reconstruction

• Κύριος συγγραφέας: Γαϊτάνης, Αναστάσιος
• Άλλοι συγγραφείς: Παναγιωτάκης, Γεώργιος
• Μορφή: Thesis
• Γλώσσα: English
• Έκδοση: 2011
• Θέματα: Positron emission tomography (PET); Medical image reconstruction; Monte Carlo methods; Iterative reconstruction algorithms; OSEM; MLEM; Stopping rules; Updating coefficients; Τομογραφία εκπομπής ποζιτρονίου; Ανακατασκευή ιατρικής εικόνας; Μέθοδοι Monte Carlo; Επαναληπτικοί αλγόριθμοι; Κριτήρια παύσης; Συντελεστές αναβάθμισης; 616.075 75

The aim of this Thesis is the development of stopping rule methods for the MLEM and OSEM algorithms used in image reconstruction positron emission tomography (PET). The development of the stopping rules is based on the study of the properties of both algorithms. Analyzing their mathematical expressions, it can be observed that the pixel updating coefficients (PUC) play a key role in the upgrading process of the reconstructed image from iteration k to k+1. For the analysis of the properties of the PUC, a PET scanner geometry was simulated using Monte Carlo methods. For image reconstruction using iterative techniques, the calculation of the transition matrix is essential. And it fully depends on the geometrical characteristics of the PET scanner. The MLEM and OSEM algorithms were used to reconstruct the projection data. In order to compare the reconstructed and true images, two figures of merit (FOM) were used; a) the Normalized Root Mean Square Deviation (NRMSD) and b) the chi-square χ2. The behaviour of the PUC C values for a zero and non-zero pixel in the phantom image was analyzed and it has been found different behavior for zero and non-zero pixels. Based on this assumption, the vector of all C values was analyzed for all non-zero pixels of the reconstructed image and it was found that the histograms of the values of the PUC have two components: one component around C(i)=1.0 and a tail component, for values C(i)<1.0. In this way, a vector variable has been defined, where I is the total number of pixels in the image and k is the iteration number. is the minimum value of the vector of the pixel updating coefficients among the non-zero pixels of the reconstructed image at iteration k. Further work was performed to find out the dependence of Cmin on the image characteristics, image topology and activity level. The analysis shows that the parameterization of Cmin is reliable and allows the establishment of a robust stopping rule for the MLEM algorithm. Furthermore, following a different approach, a new stopping rule using the log-likelihood properties of the MLEM algorithm has been developed. The two rules were evaluated using the independent Digimouse phantom. The study revealed that both stopping rules produce reconstructed images with similar properties. The same study was performed for the OSEM algorithm and a stopping rule for the OSEM algorithm dedicated to each number of subset was developed.