Image analysis – Applications – Biomedical applications
Reexamination Certificate
2000-09-25
2001-12-18
Mancuso, Joseph (Department: 2621)
Image analysis
Applications
Biomedical applications
C378S004000
Reexamination Certificate
active
06332035
ABSTRACT:
FIELD OF THE INVENTION
This invention relates to imaging, and more particularly, to the high speed 3D Radon reprojection of tomographic images in three dimensions.
BACKGROUND OF THE INVENTION
Tomographic images are created from line or plane integral measurements of an unknown object at a variety of orientations. These integral measurements, which may represent measurements of density, reflectivity, etc., are then processed to yield an image that represents the unknown object. Projection data generated in this manner is collected into a sinogram, and the sinogram is processed and backprojected to create the image. Tomographic reconstruction is the technique underlying nearly all of the key diagnostic imaging modalities, including X-ray Computed Tomography (CT), Positron Emission Tomography (PET), Single Photon Emission Count Tomography (SPECT), certain acquisition methods for Magnetic Resonance Imaging (MRI), and newly emerging techniques such as electrical impedance tomography (EIT) and optical tomography.
The 3D Radon transform is the mathematical model that describes some of the true volumetric (i.e., 3D) imaging technologies. A single 3D Radon projection is a one dimensional function corresponding to a collection of integrals of the object function on planes at a given orientation and at different positions. The 3D Radon transform of the object is the collection of its 3D Radon projections at all orientations. In MRI, projection mode imaging collects samples of the 3D Radon transform. Projection mode imaging is considered a good candidate for imaging in the presence of patient motion. The 3D radon transform also underlies true 3D Synthetic Aperture Radar (SAR) imaging.
One of the most important applications of the 3D Radon transform (and of reprojection) is in the context of cone-beam tomography, which is an emerging imaging modality in medical imaging. Through specialized preprocessing methods, such as those described in U.S. Pat. Nos. 5,124,914, 5,257,183, 5,444,792 and 5,446,776, which are incorporated by reference in their entirety, cone-beam tomographic data can be mapped into samples of the 3D Radon transform. The problem of recovering the volume being imaged from these samples can then be solved in terms of 3D Radon transforms.
Regardless of how the data is collected, either by MRI, SAR, etc., the reconstruction problem for the 3D Radon transform often involves artifact correction. This artifact correction process involves the reprojection operation, which simulates the physical tomographic acquisition system. Furthermore, the complexity requirements are high for volumetric imaging, necessitating the use of fast algorithms. Indeed, the problem is so computationally demanding that true cone-beam tomography is not yet practical, although the hardware is becoming readily available.
Another important application of reprojection is in the realm of partial data problems. In many modalities, only partial data is available. Examples include 3D SAR and cone-beam tomography with any incomplete source orbit (such as single circle). The reconstruction problem in these cases requires the use of iterative reconstruction algorithms that require repeated reprojection operations, which also need significant acceleration to become practical. Reprojection can be used in conjunction with the hierarchical 3D Radon backprojection, as described in U.S. patent application Ser. No. 539,074, filed Mar. 30, 2000, in iterative reconstruction algorithms which require both reprojection and backprojection at each step.
Several methods of reprojection have been developed. Direct reprojection represents one method for reprojection in which the physical system is simulated using the defining equations for the 3D radon transform. This is a relatively slow method, which requires O (N
5
) computations for calculating O (N
2
) projections from an N×N×N sized volume.
The Separation method takes advantage of the fact that for certain distributions of the projection orientations, the 3D problem can be decoupled into a set of 2D problems. These 2D problems can then be solved using direct methods. In this way, a fast algorithm is constructed, that computes O (N
2
)projections from an N×N×N sized volume in O (N
4
) computations.
The Fourier Slice approach described in Swedish Patent Application No. 9700072-3, filed Jan. 14, 1997, works by the Fourier Slice Theorem, which states that the Fourier transform of the projections are “slices” of the Fourier transform of the object, and can thus be computed efficiently by interpolation from the 3D Fourier transform of the volume.
In the Linogram approach, a series of interpolations and resampling operations are used to transform the data into samples of the object's Fourier transform on a special grid. The object can then be recovered efficiently using Fast Fourier Transforms. This is also a fast algorithm which requires O (N
3
log
2
N) computations for computing O (N
2
) projections from an N
3
volume.
In all, there is a need for fast, efficient and accurate reprojection processes for 3D images.
SUMMARY OF THE INVENTION
In keeping with one aspect of this invention, a method for reprojecting a 3D image (3D array of numbers) into a 3D Radon sinogram includes the steps of dividing the three-dimensional image into sub-volumes as small as one volumetric picture element (voxel), and reprojecting the sub-volumes at a small number of orientations to form a plurality of sub-sinograms. These sub-sinograms are then successively aggregated and processed to form a full sinogram for the initial volume.
The method uses two processes to aggregate the sub-sinograms. In one process, aggregation is exact, and in the other process, aggregation is an approximation. The first process is accurate, but relatively slow, and the second process is faster, but less accurate. By performing some aggregations with the exact process and some aggregations with the approximate process, switching between the two processes in any of a number of suitable ways, an accurate result can be obtained quickly.
REFERENCES:
patent: 4042811 (1977-08-01), Brunnett et al.
patent: 4149247 (1979-04-01), Pavkovich et al.
patent: 4217641 (1980-08-01), Naparstek
patent: 4491932 (1985-01-01), Ruhman et al.
patent: 4616318 (1986-10-01), Crawford
patent: 4626991 (1986-12-01), Crawford et al.
patent: 4709333 (1987-11-01), Crawford
patent: 4714997 (1987-12-01), Crawford et al.
patent: 4718010 (1988-01-01), Fulii
patent: 4858128 (1989-08-01), Nowak
patent: 4930076 (1990-05-01), Meckley
patent: 4991093 (1991-02-01), Roberge et al.
patent: 5008822 (1991-04-01), Brunnett et al.
patent: 5136660 (1992-08-01), Flickner et al.
patent: 5224037 (1993-06-01), Jones et al.
patent: 5229934 (1993-07-01), Mattson et al.
patent: 5243664 (1993-09-01), Tuy
patent: 5253308 (1993-10-01), Johnson
patent: 5300782 (1994-04-01), Johnston et al.
patent: 5375156 (1994-12-01), Kuo-Petravic et al.
patent: 5396528 (1995-03-01), Hu et al.
patent: 5438602 (1995-08-01), Crawford et al.
patent: 5552605 (1996-09-01), Arata
patent: 5559335 (1996-09-01), Zeng et al.
patent: 5579358 (1996-11-01), Lin
patent: 5625190 (1997-04-01), Crandall
patent: 5654820 (1997-08-01), Lu et al.
patent: 5727041 (1998-03-01), Hsieh
patent: 5778038 (1998-07-01), Brandt et al.
patent: 5796803 (1998-08-01), Flohr et al.
patent: 5805098 (1998-09-01), McCorkle
patent: 5825031 (1998-10-01), Wong et al.
patent: 5848114 (1998-12-01), Kawai et al.
patent: 5862198 (1999-01-01), Samarasekera et al.
patent: 5901196 (1999-05-01), Sauer et al.
patent: 6028907 (2000-02-01), Adler et al.
patent: 6108007 (2000-08-01), Shochet
patent: 611181 (1994-08-01), None
G.L. Zeng; “A Rotating and Warping Projector/Backprojector for Fan-Beam and Cone-Beam Iterative Algorithm”,IEEE Transactions on Nuclear Science, vol. 41, No. 6, pp. 2807-2811; Dec. 1994.
Gary H. Glover et al.; “An Algorithm for the Reduction of Metal Clip Artifacts in CT Reconstructions”;Medical Physics, vol. 8, No. 6, pp. 799-807; Nov./Dec. 1981.
McCorkle et al.; “An Order N2log(N) Backprojector algorithm
Basu Samit
Bresler Yoram
Greer Burns & Crain Ltd.
Mancuso Joseph
Tabatabai A
The Board of Trustees of the University of Illinois
LandOfFree
Fast hierarchical reprojection algorithms for 3D radon... does not yet have a rating. At this time, there are no reviews or comments for this patent.
If you have personal experience with Fast hierarchical reprojection algorithms for 3D radon..., we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Fast hierarchical reprojection algorithms for 3D radon... will most certainly appreciate the feedback.
Profile ID: LFUS-PAI-O-2581086