Image Reconstruction Methods Based on Block Circulant System Matrices

摘要

An iterative image reconstruction method used with an imaging system that generates projection data, the method comprises: collecting the projection data; choosing a polar or cylindrical image definition comprising a polar or cylindrical grid representation and a number of basis functions positioned according to the polar or cylindrical grid so that the number of basis functions at different radius positions of the polar or cylindrical image grid is a factor of a number of in-plane symmetries between lines of response along which the projection data are measured by the imaging system; obtaining a system probability matrix that relates each of the projection data to each basis function of the polar or cylindrical image definition; restructuring the system probability matrix into a block circulant matrix and converting the system probability matrix in the Fourier domain; storing the projection data into a measurement data vector; providing an initial polar or cylindrical image estimate; for each iteration; recalculating the polar or cylindrical image estimate according to an iterative solver based on forward and back projection operations with the system probability matrix in the Fourier domain; and converting the polar or cylindrical image estimate into a Cartesian image representation to thereby obtain a reconstructed image.