WO2008017076A2 - Iterative methods for dose reduction and image enhancement in tomography - Google Patents
Iterative methods for dose reduction and image enhancement in tomography Download PDFInfo
- Publication number
- WO2008017076A2 WO2008017076A2 PCT/US2007/075220 US2007075220W WO2008017076A2 WO 2008017076 A2 WO2008017076 A2 WO 2008017076A2 US 2007075220 W US2007075220 W US 2007075220W WO 2008017076 A2 WO2008017076 A2 WO 2008017076A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- fourier
- recited
- space
- data
- projections
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 160
- 238000003325 tomography Methods 0.000 title abstract description 16
- 230000009467 reduction Effects 0.000 title abstract description 11
- 238000012986 modification Methods 0.000 claims abstract description 47
- 230000004048 modification Effects 0.000 claims abstract description 47
- 238000003384 imaging method Methods 0.000 claims abstract description 39
- 238000004422 calculation algorithm Methods 0.000 claims description 36
- 230000009466 transformation Effects 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000005855 radiation Effects 0.000 abstract description 18
- 238000007796 conventional method Methods 0.000 abstract description 15
- 238000013459 approach Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 20
- 238000002591 computed tomography Methods 0.000 description 19
- 230000000875 corresponding effect Effects 0.000 description 16
- 230000000153 supplemental effect Effects 0.000 description 13
- 230000005540 biological transmission Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 9
- 238000012804 iterative process Methods 0.000 description 8
- 230000015556 catabolic process Effects 0.000 description 7
- 238000006731 degradation reaction Methods 0.000 description 7
- 238000002603 single-photon emission computed tomography Methods 0.000 description 7
- 230000004907 flux Effects 0.000 description 6
- 239000002245 particle Substances 0.000 description 6
- 238000002600 positron emission tomography Methods 0.000 description 6
- 239000007787 solid Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 108060003552 hemocyanin Proteins 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 241000713772 Human immunodeficiency virus 1 Species 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000000844 transformation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004627 transmission electron microscopy Methods 0.000 description 2
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- 206010073306 Exposure to radiation Diseases 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000001493 electron microscopy Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 210000003463 organelle Anatomy 0.000 description 1
- 239000012857 radioactive material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009291 secondary effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Definitions
- This invention pertains generally to diagnostic and interventional imaging using beams and detectors, and more particularly to systems and methods for creating a three dimensional cross sectional image of an object by the reconstruction of its projections that have been iteratively refined through modification in object space and Fourier space.
- Tomography is an important diagnostic and interventional tool in medicine, biology, geophysics, materials science, archeology, astronomy and the other sciences. It is an imaging technique that produces a cross sectional 2D or 3D image of the internal structures of an object from its projections.
- the typical tomography apparatus contains a radiation source and a detector that is rotated around an axis extending perpendicularly from the plane of the examination table. Projections of the patient or object are conventionally taken at equal angle intervals such that the angle of the radiating source with respect to the isocenter of the scanner changes by a fixed amount from one projection to the next. Images have been produced from several different beam sources including x-rays, electrons, gamma rays, ions, neutrons, sound waves and others.
- a central problem in tomographic imaging is the radiation dose that is directed to the patient or biological specimen as the result of the imaging procedure.
- An additional problem in tomographic imaging is the degradation of resolution and other image quality parameters due to the problem of missing or incomplete sets of projection data. Missing projection data can arise from radiation dose restrictions or from practical mechanical limitations in the imaging procedure or imaging system.
- One example is the missing wedge problem that occurs in electron microscopy, i.e., specimens cannot be tilted beyond ⁇ 70° and the data in the remaining ⁇ 20° projections are missing.
- the filtered back projection scheme is computationally fast, but does not offer any solutions for the problem of excessive radiation dose exposure and the problem of image degradation due to missing projection data.
- FBP suffers from inaccuracies due to inherent interpolation problems that occur in the back projection process.
- images reconstructed with the FBP method often contain artifacts that degrade the resolution and overall image quality.
- other reconstruction algorithms exists that are not in general use because they are computationally expensive under practical imaging conditions and also suffer from the problem of interpolation, which reveals itself when the forward projection process is modeled into the system matrix.
- the present invention provides a system and method for use with any tomographic imaging system that reconstructs an object from its projections.
- These systems include variants of transmission CT (X-ray Computed
- CT/CAT Tomography
- X-ray MicroCT X-ray MicroCT
- Transmission Electron Microscopy variants of emission CT
- SPECT Single Photon Emission CT
- PET Positron Emission Tomography
- combined CTIPET combined CTIPET, CT/SPECT, SPECT/PET, CT/SPECT/PET systems.
- the present invention is an iterative reconstruction method capable of solving for a significant portion of the missing projection data.
- the invention presents a method to eliminate interpolations present in conventional tomography.
- the method has been experimentally shown to provide higher resolution and improved image quality parameters over existing approaches.
- a primary application of the method is dose reduction since the present invention can produce an image of a desired quality with a fewer number projections than conventional methods.
- embodiments of the invention make use of an iterative algorithm that solves for the missing projections and inaccessible data points during the reconstruction process. Once the experimental projection data from the imager is obtained, the projection data is reformulated so that it can be used in accordance of the Fourier Slice Theorem, and optionally, the oversampling method.
- the missing projection data may be filled in by assigning random or preconditioned values to the missing Fourier space data.
- the missing data is calculated and the image is reconstructed from the Fourier slices by using inverse and forward Fourier transforms to iterate back and forth between object space and Fourier space while modifying the data before each transform.
- the modifications in Fourier space may include making the resulting Fourier space data of the given iteration consistent with the transformed experimental projection data.
- the modifications in object space can consist of making the pixel (or voxel) values of the object data provided by the inverse Fourier space conform to a set of user defined values and/or conditions. The cycle of refinement will continue until a termination condition is satisfied and a final image is provided.
- Imaging of the raw projections begins by inputting the projection data. If necessary, the projection data is reformulated so that it can be utilized in accordance to the Fourier Slice Theorem. Optionally, the projections may be padded with a suitable number of zeros prior to placement on the Fourier space grid depending on whether oversampling is desired and the degree of oversampling. Padding with zeros will result in the object space data to be surrounded by mathematical zeros. The knowledge of these zeros can then be used to solve for the missing data in the iterative process.
- the projections are Fourier transformed and placed onto a grid in Fourier space.
- the grid and the Fourier transformation that are used are dependent on one another and may vary. Selection of the grid that the transformed projections are placed upon will influence the selection of a Fourier transform algorithm.
- the preferred computerized method for performing this step is to use a pseudopolar grid in Fourier space and then perform the Fourier transform (FT) using a Fractional Fast Fourier Transform (FrFFT) algorithm set to match the points of the pseudopolar grid at the given projection angle and slice location.
- FT Fourier transform
- FrFFT Fractional Fast Fourier Transform
- Other computerized implementations of the FT including the non-uniform Fast
- NUFFT NUFFT
- FFT Fast Fourier Transform
- the transformed projections are placed on the appropriate grid at the corresponding location in Fourier space. Placement of the projections on the Fourier space grid can also be accomplished by exact methods such as using the pseudopolar grid in conjunction with the FrFFT, exact or approximate mathematical transformations, or through the Gridding or other interpolation methods. [0020] Once placement on the Fourier space is complete, other optional supplemental transforms can be applied to correlate the slices. Such supplemental transform preferably includes a FT across all slices in the Fourier space.
- the locations and values of the Fourier transformed experimental data in Fourier space may also be stored. These values can be used to force the Fourier space data at these locations in future iterations to conform to the experimental data.
- the derived Fourier space data at these locations is replaced during each iteration by the transformed experimental values to ensure that the results of the reconstruction are consistent with the experimental data.
- the Fourier space data may be modified by assigning values to non- experimental locations and points beyond a resolution circle in Fourier space through preconditioning methods or by the random assignment of values.
- Preconditioning methods include assigning values to the missing Fourier space data based on the reconstruction of the projection data by some other method.
- Other modifications may include the enforcement of Fourier space conditions such as symmetry conditions on the phase or amplitude of the
- the modified Fourier space data is subject to an inverse Fourier
- Transform to provide an object space image. For example, this may be done slice by slice with a 2D algorithm or together with a 3D algorithm.
- the computerized implementation of this step will vary depending on the grid that is used.
- the Fourier space is on a pseudopolar grid while the object space is on a Cartesian grid
- the inverse Fourier Transform is an inverse Pseudopolar Fast Fourier Transform (IPPFFT) algorithm.
- IPPFFT inverse Pseudopolar Fast Fourier Transform
- Transform can be modified to make the object space pixel values conform to previously defined constraints.
- Constraints are restrictions on the object space image that can be used to aid in solving the reconstructed image.
- a positivity constraint could require the pixel values in object space to be positive.
- the pixel values that are negative are modified so that their modified values are closer to zero. This may be achieved by simply setting negative values to zero or pushing the values towards to zero.
- Other constraints that may be applied include a requirement that the object space image must be real. Complex values can be transformed to real values by taking the amplitude or real component of the given values. Furthermore, certain designated pixels can be driven towards values that are defined by the user.
- pixels in the object space can be modified based on general physical conditions which may restrict the pixel values to some maximum value. Pixels with values above the prescribed maximum values may be modified. In one embodiment, pixels with values that exceed a maximum are mathematically modified so that the modified value is closer to the maximum value than the unmodified value.
- the physical image region is preferably surrounded by a known number of pixels with values equal to zero.
- the pixels in this region referred to as the support, can be used as an additional constraint. Consequently, the object space can be further modified in the support region by pushing the non-zero pixel values in this region towards zero. This method aids the reconstruction by providing additional constraints by which the object space can be modified.
- the modified object space data is then subject to a Fourier transformation to provide Fourier space data. In one embodiment, a 2D or a 3D Fourier transformation is applied to the modified object space image.
- the Fourier space data is obtained by mathematically calculating the projections from the modified object space image and comparing the mathematically calculated projections with experimental projections.
- the mathematically calculated projections are also made compliant with the experimental projections. This may be achieved by replacing the points in the mathematically calculated projections at which corresponding experimental projections exist by the experimental projections.
- An appropriate 1 D Fourier Transform is applied to the calculated projections and the transformed projections are then placed onto the appropriate Fourier grid to arrive at the Fourier space data.
- the process of mathematically calculating the projections also includes a method by which the object space data is Fourier transformed using the 2D PPFFT algorithm on each slice.
- the values along the lines of the pseudopolar grid going through the origin of pseudopolar grid in Fourier space are extracted and then the mathematical projections are calculated using an inverse FrFFT algorithm in accordance to the Fourier Slice Theorem.
- the transformed Fourier space data can then be modified to complete the cycle.
- an additional supplemental transform may be applied to the Fourier data.
- an error function may be constructed to quantify the difference between the derived Fourier space data with the experimentally measured Fourier data at points where the measured data exist. If the error function reaches its best value as compared to the previous iteration, the Fourier space data corresponding for the iteration is stored.
- the Fourier space data corresponding to the best error is inverted to provide the best object space image according to the designated error function.
- the modification of the Fourier space data includes the replacement of the calculated Fourier slices at the locations of the experimentally measured slices with the stored experimentally measured Fourier slices at points inside the resolution circle. This method ensures that the reconstruction is consistent with the experimental data at each step of the iteration.
- the calculated slices are pushed towards the measured slices.
- the calculated slices are pushed towards the measured slices only if the calculated slices at a given location differ from the measured slices by a certain threshold percentage.
- the modification of the Fourier space can also include the application of other Fourier space conditions.
- the cycle can proceed until a termination condition is satisfied such as a predetermined number of iterations or a monitored error reaches a certain value.
- a final image can be produced from the results of the refined data.
- a method is provided that can be adapted to any imaging system that reconstructs an object from its projections and allows projection acquisition from parallel beam and non- parallel beam imagers.
- An aspect of the invention is to provide a method that allows a subset of tomographic slices to be reconstructed independently of the other slices and thus permitting parallel computing.
- Another aspect of the invention is to provide a method that produces images with improved image quality parameters than found with conventional methods.
- image quality parameters include resolution, contrast, signal to noise ratio and parameters specifying the degradation of an image by artifacts.
- Another aspect of the invention is to provide a method for reducing the dose that is imparted to the patient or the imaging object in tomographic imaging systems.
- An aspect of the invention is to provide a system and method that can provide detailed two and three dimensional images from a reduced number of projections as well as reduced radiation flux over existing methods to reduce the exposure of the subject to radiation.
- Another aspect of the invention is to provide a method for reducing the acquisition time of projection data in tomographic imaging systems. Fast acquisition times are particularly important in imaging an object that moves during the imaging procedures. Such objects may include the organs of a living organism that are in constant motion like the heart and lungs. Motion artifacts in tomographic imaging can be reduced with fast acquisition times.
- An aspect of the invention is to provide a method that allows a subset of tomographic slices to be reconstructed independently of the other slices and thus permitting parallel computing.
- An aspect of at least one embodiment of the invention includes the acquisition of a partial or full set of equally sloped (or pseudopolar) transmission or emission projections.
- An equally sloped or pseudopolar projection data set is defined as a transmission or emission projection data set that when converted from the source geometry to the corresponding parallel beam geometry (either exactly or by interpolation), the parallel projections at a given angle lie (within experimental error) along one of the lines of a chosen pseudopolar grid.
- a method for reconstructing an object from its projections comprising the steps of obtaining projection data from an imager; mapping Fourier transformed projection data in Fourier space; converting iteratively the projection data from
- a computer program product and system comprising a computer-readable storage medium having computer-readable program code portions stored therein and providing for associative relevancy knowledge profiling, the computer program product comprising a first program code configured for obtaining projection data; a second program code configured for mapping Fourier transformed projection data in Fourier space; a third program code configured for converting iteratively the projection data from Fourier space to object space; a fourth program code configured for modifying said projection data in Fourier space and in object space to produce progressively modified projection data; and a fifth program code configured for generating an image from the modified projection data.
- a system and method that permits efficient tomographic imaging by 1 ) converting projection data from source geometry to parallel beam geometry, if necessary. 2) The projections are padded with an appropriate number of zeros to match the oversampling of the pseudopolar grid. 3) The projections at each angle are mapped to their corresponding line on a pseudopolar grid in Fourier space using the Fractional Fast Fourier Transform in accordance to the Fourier Slice Theorem. 4) The values and locations corresponding to the transformed measured experimental data are stored. 5) Optionally, random or preconditioned values can be assigned to the pseudopolar grid points that are missing, unknown, inaccessible or have incomplete projection data.
- IPPFFT Inverse Pseudopolar Fast Fourier Transform
- the pseudopolar grid Since the pseudopolar grid is oversampled, the physical object should be surrounded by mathematical zeros in the support region. Consequently, this condition may be used as an additional constraint by pushing non-zero values in the support towards zero. Additionally, if information about pixels in any region is known and or if specific pixels are desired to be driven to desired values, these desired or known values can be used as constraints by pushing the original values to the desired or known values. 8) The forward Pseudopolar Fast Fourier Transform (PPFFT) is applied to the image to obtain the calculated frequency data. 9) An error function may be calculated from the difference of the calculated Fourier space data and the measured transformed experimental data at locations where experimental data exist.
- PPFFT forward Pseudopolar Fast Fourier Transform
- the calculated Fourier space data are replaced by the measured transformed experimental projections at locations where the measured data exist to make the reconstruction consistent with the experimental projections.
- the Fourier space data is further modified by any additional Fourier space constraints. 10) Steps 6 to 9 are repeated a termination condition is satisfied.
- step 9 of each iteration if the error function reaches its lowest (best) value as compared to the previous iteration, the Fourier space data corresponding for the iteration is stored.
- the Fourier space data corresponding to the best error is inverted by the IPPFFT, and cropped to exclude the support region, to yield the final reconstructed image.
- the termination condition includes a method by which the iterative algorithm is terminated if the error has not improved by a certain designated percentage in a certain number of iterations. Other termination conditions include satisfying a certain number of iterations or a certain designated error.
- FIG. 1 is a flow diagram of a method for image enhancement and dose reduction in tomography according to one embodiment of the present invention.
- FIG. 2 is a flow diagram of an alternative embodiment of a method for enhanced tomographic imaging and dose reduction according to the invention.
- FIG. 3 is a flow diagram illustrating a method of enhanced equally sloped tomographic imaging and dose reduction according to an alternative embodiment of the present invention.
- FIG. 4(a)-4(d) are cross correlated comparisons of 11 hemocyanin molecules iso-surface renderings imaged by conventional methods and with the methods according to the present invention.
- FIG. 5(A)-(D) are 3D image reconstructions of an HIV-1 virus like particle from cryo-EM data imaged by conventional methods and with the methods according to the present invention.
- FIGS. 1 through FIG. 5(D) the present invention is embodied in the compositions and methods generally shown in FIGS. 1 through FIG. 5(D). It will be appreciated that the apparatus may vary as to configuration and as to details of the parts, and that the methods may vary as to the specific steps and sequence, without departing from the basic concepts as disclosed herein.
- the present invention provides powerful new techniques for tomographic imaging that produces significantly higher spatial resolution, contrast, and signal to noise ratios than found with conventional methods in Computed Tomography (CT), for example. Problematic interpolations that are inherent in conventional methods of tomography are also eliminated by the invention. Furthermore, the present invention provides a method for significantly reducing the radiation dose that is required to be delivered to the patient or object. Additionally, the invention presents a method for performing high speed imaging by reducing the scan acquisition time. [0053] The method in general can be applied to any imaging system that reconstructs the cross section of an object from its projections.
- CT or CAT X-ray Computed Tomography
- PET Positron Emission Tomography
- SPECT Single Photon Emission Tomography
- Ultrasound Fluoroscopy
- FIG. 1 depicting one embodiment 10 of the invention.
- projection data from a parallel beam system is acquired at block 12.
- projection data can be acquired and then reformulated by a parallel-beam equivalent projection data set.
- Such reformulation at block 12 can be performed by mathematical transformations in an exact or approximate manner from the system geometry to an equivalent parallel beam geometry.
- the reformulation can be performed by existing 'rebinning' algorithms present in the art which specifically change the projection data from the given system to a parallel beam equivalent system in an exact or approximate method.
- Another related method includes an interpolation of the sinogram along the desired parallel projection angular location to arrive at a parallel beam projection.
- a further method for converting non-parallel cone or fan beam data to the parallel beam format comprises a) reconstructing the raw cone or fan beam projection data using standard algorithms for fan/cone beams known in the art, b) operation of a 2D Fourier transform (FT) across all slices, c) extracting the Fourier transformed projections at the desired location from the Fourier data with the use of the Fourier Slice Theorem to arrive at the transformed parallel beam projections, and, d) optionally, if an object space projection data set is desired, operating an inverse 1 D Fourier transform to arrive at the parallel beam projection data in object space.
- FT 2D Fourier transform
- a Fourier Transform of each projection is taken at block 14.
- the type of Fourier transform that is selected will depend on the type of grid in Fourier space that is desired.
- 'Fourier space' is defined as a space in which at least one Fourier transformation has been applied to a given data set.
- the Fourier space can be constructed from other transformations in addition to the Fourier transform. In some situations, the Fourier space used here is similar or equivalent to the phrase 'frequency domain' or 'k-space' in the art.
- the Fourier transform of the projections at a plurality of angles and slices at block 14 is a Fractional Fourier Transform set to match points and lines of a pseudopolar grid for each angle in accordance to the Fourier Slice Theorem.
- Other computerized implementations of the FT such as the non-uniform Fast Fourier Transform (NUFFT) and the conventional Fast Fourier Transform (FFT), can be utilized at this point in accordance with the Fourier Slice Theorem.
- NUFFT non-uniform Fast Fourier Transform
- FFT conventional Fast Fourier Transform
- the Fourier transforms of the projections are placed on a grid in Fourier space at block 16.
- a pseudopolar grid has points that sit at the intersections of linearly expanding concentric squares with angularly spaced rays.
- one embodiment of the present invention acquires the data in such a manner that the projections lay (when converted into parallel geometry) partially or fully along a selected number of equally sloped lines of a pseudopolar grid.
- the projections at each angle are then exactly mapped to their corresponding lines on the pseudopolar grid in the Fourier space, preferably using the Fractional Fast Fourier Transform (FrFFT).
- Fractional Fast Fourier Transform Fractional Fast Fourier Transform
- the placement of transformed projections onto a grid in Fourier space may also be performed with an exact or approximate mathematical transformation, or by interpolating the location and values of a transformed projection on the Fourier space grid.
- the placement of transformed projections onto a grid in Fourier space may also be performed by the
- the present invention is not restricted to a specific grid in Fourier space.
- Other grids that may be used for the Fourier space data include Cartesian type, polar type, 3D polar (i.e., spherical) type, 3D pseudopolar and others.
- the Fourier space data of the projections is then modified at block 18 of
- Modifications at block 18 may include filling in missing projection data by assigning values to non-experimental locations in Fourier space.
- the non- experimental locations can be visualized as regions in Fourier space that the transformed experimental projections do not cover when they are placed on the Fourier space grid.
- the non-experimental regions can be more generally defined as locations which the user or the algorithm designates as non-experimental.
- the locations and values of the Fourier transformed experimental data in Fourier space may also be stored. These stored values can be used to force the Fourier space data at these locations in future iterations to conform to the experimental data. For example, in one embodiment, during each iteration, the stored Fourier space experimental data are used to replace the derived
- the Fourier space data may also be modified by assigning values to non-experimental locations and points beyond a resolution circle in Fourier space through preconditioning methods or the random assignment of values.
- the resolution circle is a defined circular region in Fourier space which accordingly separates grid points into different regions as defined by being inside, on, or outside the resolution circle.
- Preconditioning methods include assigning values to the missing Fourier space data based on the reconstruction of the projection data by any method in the art.
- the Fourier data and grid may also be modified at this point by extending the grid and padding the Fourier space data with zeros. Other modifications may include the enforcement of Fourier space conditions such as symmetry conditions on the phase or amplitude of the Fourier data. If optional supplemental transforms were performed, they may be inverted at the final stage of this block.
- the modified Fourier space data begins an iterative process at block 20 that is able to solve for the missing/incomplete projection data and or experimentally inaccessible data points in Fourier space.
- These missing projection data typically arise in computed tomography due to the limitations in administered dose or system design or resources in practical applications. Shifting between Fourier space and object space while modifying the data in each space between shifts permits iterative refinements to the data until a final image is reconstructed. The resolution is improved as a result of removing the degradation of image quality that is present due to missing projection data.
- an inverse Fourier transform is performed upon the Fourier space data of block 20 to provide data in the real or object space at block 24.
- inverse Fourier transforms may be used at block 20 depending on the Fourier and object space grids that are used. For example, if the Fourier space is on a pseudopolar grid and the object space is on a Cartesian grid, the inverse pseudopolar fast Fourier transform (IPPFFT) may be used. Alternatively, if both grids are Cartesian, a conventional inverse FFT routine may be used. The inverse FT can be performed slice by slice or together at once or on grouped subsets of slices. The slice by slice method may be preferred as it easily allows for parallel computing by distributing subsets of slices across different computers.
- IPPFFT inverse pseudopolar fast Fourier transform
- termination of the process is evaluated to determine if a termination event has occurred.
- One termination event is the occurrence of a predetermined number of iterations.
- Another termination event can be the occurrence of a monitored error or image quality parameter function.
- Yet another termination event can be the occurrence of situation where the error or image quality parameter has not improved in a prescribed number of iterations.
- the final image may be produced by inverting the Fourier space data corresponding to the best error or image quality parameter, and optionally cropping the object space in accordance to the oversampling used.
- the final image does not necessarily have to be based upon the calculated data from the last iteration.
- the final image can be based on the data from previous iterations, or new data that is formed by combining data from previous iterations.
- the final image can be the average of the object space data from a number of different iterations.
- the termination event can be made to occur at any step in the iterative cycle.
- the final image can be based upon any step of the iterative process.
- the object space data at block 24 can be modified at block 28 to provide modified object space data at block 30.
- Object space data modifications in general make the object space pixel values conform to defined constraints and thereby provide an efficient method to iteratively solve the reconstruction problem in a manner consistent with the constraints.
- An object space pixel can be made to conform to an expected or desired value (i.e., a constraint) in any mathematical manner.
- a constraint i.e., a constraint
- One of the simplest and most common methods for achieving this is by simply replacing the pixel value by a value that is dictated by the constraint. For example, if a constraint dictates that a specific pixel should have the numerical value of zero, the value of the pixel may be changed to zero in the given iteration.
- Another method for making a specific pixel conform to a constraint is to modify it so that its value is not exactly the value of the constraint, but after the modification the value of the pixel is potentially closer to the value dictated by the constraint.
- This method is referred to as 'pushing' the value of a given pixel towards a different value.
- this method can be achieved in different ways.
- the equation for the modification can be any function of the current iteration pixel value and any number of the previous pixel values.
- the above modifications can be formulated in many different mathematical ways. However they all constitute the various implementations of the method presented in this invention.
- a given pixel value is only changed if it differs from the constraint by a certain percentage threshold.
- a positivity constraint is used to require the pixel values in object space to be positive.
- the pixel values that are negative are modified so that their modified values are closer to zero. This may be achieved by simply setting negative values to zero or pushing the values towards to zero.
- negative values are multiplied by a negative constant to make them a positive number.
- Other constraints can include conditions dictating that the object space image must be real. Complex values can be transformed to real values by taking the amplitude or real component of the given values.
- certain designated pixels can be driven towards defined desired values.
- pixels in the object space can be modified based on general physical conditions.
- a physical constraint may demand the pixel values to be below some maximum value. Pixels with values above the prescribed maximum values can be pushed towards the designated maximum.
- Other constraints may include conditions demanding that the object space data change continuously.
- Other constraints may include conditions on the mathematical derivative of the object space data.
- a given pixel may be modified so that the mathematical derivative along a certain direction is constrained to some numerical range.
- n is constant that can be taken to be equal the number of v's appearing in the abs parenthesis
- alpha and beta are arbitrary constant
- abs represents the absolute value function.
- Such a modification is useful for eliminating or de-amplifying variations that are not desired from one iteration to another.
- v 0.5
- the physical image region be surrounded by a known number of pixels (as determined by the degree of oversampling) with values equal to zeros.
- the pixels in this region referred to as the support, can be used as an additional constraint. Consequently, the object space can be further modified in the support region by pushing the nonzero pixel values in this region towards zero. This method aids the reconstruction by providing additional constraints by which the object space can be modified in that embodiment.
- a Fourier transform of the modified object space data at block 30 is applied at block 32 to provide Fourier space data at block 34.
- This transform converts the modified object space data into Fourier space data that will permit additional or different types of refinements to be made.
- the Fourier transform is a 2D Pseudopolar Fast Fourier Transform (PPFFT) applied to each slice.
- the Fourier transform is a 3D PPFFT algorithm applied to the totality or individual subsets of the object space data.
- the Fourier transform is a conventional 2D FFT applied to each slice.
- the Fourier transform is a 3D FFT algorithm applied to the totality or individual subsets of the object space data.
- projections from the new image at block 30 are calculated and made compliant with the observed experimental projections. Then at block 32, Fourier transforms of the projections are conducted and the projections are placed onto the appropriate Fourier grid in accordance to the Fourier Slice Theorem. This step can also include any additional desired transforms (such as a 1 D FT along the slice direction). [0077]
- the Fourier space data obtained at block 34 can be modified at block
- an additional supplemental transform may be applied to the Fourier data at this point if desired.
- an error function may be constructed to quantify the difference between the derived Fourier space data with the experimentally measured Fourier data at points at locations where the measured data exist.
- an error function can be constructed by calculating the normalized amplitude of the difference between the derived Fourier space data and the experimental Fourier data at the appropriate locations.
- the error function reaches its best value as compared to the previous iteration, the Fourier space data corresponding for the iteration with the best error is stored.
- the Fourier space data corresponding to the best error is inverted to provide the best object space image according to the designated error function.
- the modification of the Fourier space data includes the replacement of the calculated Fourier slices at the locations of the experimentally measured slices with the stored experimentally measured Fourier slices at points inside a resolution circle. This method ensures that the reconstruction is consistent with the experimental data at each step of the iteration.
- the calculated slices are pushed towards the measured slices. In another embodiment, the calculated slices are pushed towards the measured slices only if the calculated slices at a given location differ from the measured slices by a certain threshold percentage.
- the modification of the Fourier space includes the application of other Fourier space conditions. These constraints may include symmetry conditions on the phase and amplitude of the Fourier data. If the initial modification of Fourier space included an additional supplemental transformation, the supplemental transformation can be inverted at this point.
- the cycle begins again with the application of an inverse Fourier transform at block 22 and continues until a termination condition is satisfied at block 26 and a final image is created at block 38.
- a termination condition is satisfied at block 26 and a final image is created at block 38.
- the type of inverse and forward transforms and grids that are selected can change from cycle to cycle if desired.
- the types of modifications to the object and Fourier space data can be changed from cycle to cycle. This allows for the application of complimentary refinements to take place and missing or incomplete data points to be reliably filled.
- the raw projection data can be refined slice by slice independently. Independent slice computation permits parallel computing and greatly reduces the computing demand and the increases speed of image acquisition.
- FIG. 2 an alternative embodiment 40 of the general iterative reconstruction method is shown. Projections of a patient or object are obtained at block 42 and converted to parallel beam equivalents at block 44, if necessary. The projections are then padded with zeros at block 46 in this embodiment. The Fourier transform of the projections are taken in accordance of the Fourier Slice Theorem. As indicated previously, the computer implementation of the Fourier transform at this step may be varied in order to optimize or eliminate any interpolations in Fourier space.
- the Fourier grid is pseudopolar
- many different computer implementations of the Fourier transform such as the FFT or Fractional FFT (FrFFT) are possible.
- FFT Fractional FFT
- using the FrFFT potentially allows for the elimination of interpolation in Fourier space by varying the parameter of the
- FrFFT that controls the output spacing in such a manner that the sampling of the output matches the sampling of the pseudopolar grid in a given region.
- Transforms and grids can vary as described previously.
- additional transforms can be applied across the data in Fourier space.
- Such a transform may include applying a Fourier transform across different slices in Fourier space.
- the Fourier data can be transformed by extending the grid, and padding the Fourier data by zeros or other desired values.
- the transformed Fourier experimental values and locations are stored in block 54.
- a partial or full subset of the non-experimental region in Fourier space may be filled in by methods previously described.
- any desired Fourier space constraints i.e., conditions
- the resulting Fourier data from the above is symbolized by F',(k), where the subscript / denotes the iteration number. At this step the subscript may be taken to be 0 or 1 depending on the desired convention.
- the transformed and modified projections begin an iterative refinement cycle at block 60.
- the cycle shown in the embodiment of FIG. 2 can take two alternative paths (as shown by the dotted and solid arrows that extend out of block 70) or can produce a result influenced by a combination of both paths.
- Both cycles begin with the application of an inverse Fourier transform at block 62 of the Fourier data of block 60 to produce transformed data at block 64.
- transforms 62 are suitable depending on the grids used in Fourier space and object space as previously described.
- the transformed data at block 64 may be modified with a positivity constraint at block 66 if desired. If oversampling was used, the physical image region outside of the object should be surrounded by zeros in a region referred to as the support. Pixels in the support can then be pushed towards zero.
- any desired conditions i.e., constraints
- values may be limited to a certain range.
- the modified object space data at block 70 can go through to block 74 through different paths, as symbolized by the dotted and solid arrows.
- the algorithm only cycles through the steps depicted by the solid arrows.
- the algorithm only cycles through the path defined by the dashed lines in order to translate from block 70 to block 74.
- the algorithm chooses one or the other loop to translate from block 70 to 74.
- the algorithm produces two copies of the data in block 70, and then translates one copy through each path to block 74.
- a Fourier transform of the object space data is taken as described previously.
- the computerized algorithms for the Fourier transform can be used depending on the grids used in Fourier and object space.
- the data from block 70 is used to calculate mathematical projections at block 76.
- Various algorithms in the art that mathematically calculate projections for a given data set can be utilized.
- One method to calculate projections mathematically includes the use of the Fourier slice theorem. In this method, Fourier transforms of the object space data across all slices are taken. The transformed projections are extracted according to the Fourier Slice theorem and an inverse Fourier transform is then applied to extracted transformed projection data in accordance to the Fourier Slice
- Fourier space data at block 74 can come from one of the two cycles individually or can be a combination of the two cycles.
- Fourier space data of block 74 may be further modified and refined at block 84 by making the results compliant with stored transformed experimental values from block 54 as well as applying desired Fourier conditions at block 86 as previously described.
- the Fourier data can be transformed by extending the grid, and padding the Fourier data by zeros or other desired values.
- the modified Fourier data from block 86 may enter the cycle again at block 60.
- the cycle can be terminated at any stage of the iterative cycle to produce a final image. If the cycle is terminated at some point in Fourier space, the Fourier space data can be inverted to provide object space data as a basis for the final outputted image. Before final output, the data to be outputted can be processed in a number of ways. Such processing may include cropping the data, taking the real component of the data, taking the amplitude of data and others.
- FIG. 3 an alternative embodiment 90 of the invention is shown that is particularly preferred for use with a partial or full set of equally sloped transmission or emission parallel equivalent projections that lie (within experimental error) along one of the lines of a chosen pseudopolar grid.
- this embodiment uses a 1 D FrFFT transform in a manner consistent with the Fourier Slice Theorem and with the output sampling matching the grid sampling of the Fourier space grid at the given slice location and projection angle.
- a pseudopolar type grid in Fourier space is used that is preferably oversampled by a factor of two relative to the object space grid.
- the Fourier space grid specifically consists of a 2D pseudopolar grid at each slice location and uses a Cartesian grid in object space.
- Transformation from Fourier space to object space and back is accomplished with 2D Inverse and Forward Pseudopolar Fast Fourier Transform (IPPFFT) algorithm across all slices of the Fourier space.
- IPPFFT Inverse and Forward Pseudopolar Fast Fourier Transform
- no additional supplemental transformations are applied, and the reconstruction is performed slice by slice in a parallel computing environment.
- Projection data is obtained at block 92 of FIG. 3.
- the imaged subject, detectors, and/or the radiation source are rotated in manner that the majority of the resulting reformulated projection data, within experimental error, lies along isocenthc lines formed by grid points of a grid in Fourier space.
- An isocenthc line is a line intersecting a reasonable number of grid points as well as intersecting a central point of a slice or set of slices of a grid. Examples of such lines include the lines formed on 2D or 3D pseudopolar grids that intersect a central point in each slice of the grid or intersect the central point of a set of slices.
- Other grids include polar and spherical grids.
- One illustrative embodiment of the method includes rotating the sample in a transmission CT modality, such as the Transmission Electron Microscope, a SPECT or PET system or a medical x-ray CT/CAT scanner, in such a manner that resulting projections are parallel within experimental error to a partial set of lines defined by the Fourier space grid. If the preferred pseudopolar grid is used, the projections will reside along a partial set of pseudopolar lines.
- a transmission CT modality such as the Transmission Electron Microscope, a SPECT or PET system or a medical x-ray CT/CAT scanner
- the detectors as defined by the Gamma Camera or head(s) of the scanner can be rotated in such a manner that angles the head(s) make with the isocenter of the scanner correspond, within experimental error, to a partial set of angles made by the lines of pseudopolar grid that intersect the center of the grid.
- the detector elements can be arranged in such a way that a majority of the resulting reformulated projections are approximately parallel to a partial set of lines defined by the Fourier space grid. In the case of the pseudopolar grid, this can be achieved by arranging the detector elements along a square at each slice location as opposed to the conventional circular design.
- the projections that are acquired at block 92 and block 94 are optionally padded with zeros at block 96.
- a 1 D Fractional Fast Fourier Transform (FrFFT) is applied to the acquired projections at block 98.
- the parameter that controls the output spacing of the FrFFT should be set to match the spacing of the points of the lines of the Pseudopolar grid for the specific angle.
- the values of the transformed projections from block 98 are placed onto a pseudopolar grid at block 100.
- additional transforms such as a 1 D Fourier Transform along different slices may also be applied at block 102.
- the transformed Fourier experimental locations and values for each projection may be stored at 104 in the embodiment shown.
- values of non-experimental locations in Fourier space may be assigned and any Fourier space constrains may be applied at block 108 to produce the starting Fourier space data at block 110.
- a 2D Inverse Pseudopolar Fast Fourier Transform is performed on each two dimensional Fourier slice from the data obtained at block 110 to provide an object space image at block 114.
- Positivity constraints are applied to the support at block 116 and other constraints at block 118 to provide modified object space data at block 120.
- the modified object space data can be manipulated and can be translated to block 124 through alternative paths symbolized by the solid and dotted arrows or through a combination of paths as previously discussed in this document.
- a 2D forward Pseudopolar Fast Fourier Transform is applied to each slice of the image at block 122 in the embodiment shown in FIG. 3 to produce modified Fourier space data 124.
- projections can be calculated from the object space data from block 120 at block 126.
- the calculated projections are made to be compliant with experimental projections at block 128 as previously discussed.
- a 1 D FrFFT is then performed on the calculated projections using the appropriate output spacing parameter for a pseudopolar grid at block 130 and placed on to a pseudopolar grid at block 132.
- the modified Fourier space data at block 124 may be made compliant with transformed experimental values at block 134 and Fourier space conditions may be applied at block 136.
- the modified Fourier space data then enters the cycle again at block 110. The iterative cycle of transforms and modifications continues until a termination condition is met.
- a final image is then generated at block 138.
- this embodiment of the method resolves the interpolation problem typically found in conventional methods of tomography by matching the projections to isocentric lines of the pseudopolar grid and by matching the sampling of the FrFFT to that of the pseudopolar grid for each projection.
- the missing data can be solved and the process of solving for the missing data can be confirmed in computer simulations. Simulations generally show that after each iteration prior to the termination condition, there is a correlation between the reconstructed image with that of an ideal or perfect reconstruction that does not have a missing data increase.
- Embodiments of the present invention are described with reference to flowchart illustrations of methods and systems according to embodiments of the invention. These methods and systems can also be implemented as computer program products.
- each block or step of a flowchart, and combinations of blocks (and/or steps) in a flowchart can be implemented by various means, such as hardware, firmware, and/or software including one or more computer program instructions embodied in computer-readable program code logic.
- any such computer program instructions may be loaded onto a computer, including without limitation a general purpose computer or special purpose computer, or other programmable processing apparatus to produce a machine, such that the computer program instructions which execute on the computer or other programmable processing apparatus create means for implementing the functions specified in the block(s) of the flowchart(s).
- blocks of the flowcharts support combinations of means for performing the specified functions, combinations of steps for performing the specified functions, and computer program instructions, such as embodied in computer-readable program code logic means, for performing the specified functions.
- the computer program instructions may also be loaded onto a computer or other programmable processing apparatus to cause a series of operational steps to be performed on the computer or other programmable processing apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable processing apparatus provide steps for implementing the functions specified in the block(s) of the flowchart(s).
- the pseudopolar embodiment of the invention was applied to a 3D image experimental reconstruction of hemocyanin molecules.
- the hemocyanin molecule consists of a double-layered and hollow barrel complex about 30 nm in diameter and 35 nm in length.
- a data set was acquired by tilting the specimen from -69.4° to +69.4° with a total of 105 projections.
- Image reconstruction was carried out using the said method and the conventional method of FBP.
- the projection data from the imager was padded with zeros, Fourier transformed, placed onto a pseudopolar grid, and translated through the iterative process as described by the solid arrow loop of the diagram of
- FIG. 3 A positivity constraint was used in object space. In Fourier space, the data was updated and was made to conform to the experimental data. To perform a quantitative comparison, 11 individual hemocyanin molecules of the two reconstructions were chosen. The particles were compared with a model averaged from hundreds of hemocyanin molecules.
- FIG. 4(a) shows the cross-correlation values between the reconstructed particles and the model, which indicates that the said method improves the cross correlation values by about 40%.
- the iso-surface renderings shown in FIG. 4(b) through FIG. 4(d) further demonstrated that the said method is superior to the conventional reconstruction shown in FIG. 4(d).
- FIG. 5A through FIG. 5D image reconstruction of the HIV-1 like particle from cryo-EM data is shown.
- FIG. 6 (A) and FIG. 6(B) show sectioned images in the XY and YZ plane with 15-slice averaging and no denoising.
- the inset shows some of the gold fiducial markers reconstructed by conventional methods.
- the reconstruction comprised reformulating the projections and padding them with zeros. Thereafter a FrFFT was applied to the reformulated projections and the transformed projections were placed on a pseudopolar grid.
- the data was translated through the iterative process as depicted in the solid line loop of FIG. 3 and the reconstruction was terminated after 100 iterations and a final image was produced.
- the images from the application of the method in FIG. 5(A) and (B) show a better contrast and look sharper than the conventional images and potentially new detail can be seen in the reconstruction.
- the invention has been shown experimentally to yield higher resolution images than produced by conventional reconstruction methods such as
- FBP Filtered Back Projection
- one of the most important benefits of the invention is the reduction in exposure to radiation in tomographic imaging. It is possible to lower the radiation dose to the patient or object through either lowering the flux of radiation or by decreasing the number projections taken in transmission CT scans due to the increased accuracy of reconstructions and image quality parameters such as resolution, contrast, and signal to noise ratio provided by the present invention.
- the flux can be simply lowered by altering the current to the source (for example the mAs of an X-ray tube).
- a lower number of projections can be acquired in transmission CT thereby lowering the amount of radiation transmitted through the patient or object.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Mathematical Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Algebra (AREA)
- Biomedical Technology (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
Description
Claims
Priority Applications (9)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009523075A JP2009545409A (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image improvement in tomography |
CA002659467A CA2659467A1 (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image enhancement in tomography |
MX2009001136A MX2009001136A (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image enhancement in tomography. |
EP07840699A EP2049019A4 (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image enhancement in tomography |
BRPI0715216-7A BRPI0715216A2 (en) | 2006-08-03 | 2007-08-03 | Method for reconstructing an image representation of an object from its projections and for treating tomographic images and a computer program product |
AU2007281076A AU2007281076A1 (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image enhancement in tomography |
IL196694A IL196694A (en) | 2006-08-03 | 2009-01-25 | Iterative methods for dose reduction and image enhancement in tomography |
US12/363,079 US8270760B2 (en) | 2006-08-03 | 2009-01-30 | Iterative methods for dose reduction and image enhancement in tomography |
US12/978,093 US8442353B2 (en) | 2006-08-03 | 2010-12-23 | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US83555206P | 2006-08-03 | 2006-08-03 | |
US60/835,552 | 2006-08-03 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/363,079 Continuation US8270760B2 (en) | 2006-08-03 | 2009-01-30 | Iterative methods for dose reduction and image enhancement in tomography |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2008017076A2 true WO2008017076A2 (en) | 2008-02-07 |
WO2008017076A3 WO2008017076A3 (en) | 2008-07-03 |
Family
ID=38997903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2007/075220 WO2008017076A2 (en) | 2006-08-03 | 2007-08-03 | Iterative methods for dose reduction and image enhancement in tomography |
Country Status (11)
Country | Link |
---|---|
US (1) | US8270760B2 (en) |
EP (1) | EP2049019A4 (en) |
JP (1) | JP2009545409A (en) |
KR (1) | KR20090057972A (en) |
CN (1) | CN101495038A (en) |
AU (1) | AU2007281076A1 (en) |
BR (1) | BRPI0715216A2 (en) |
CA (1) | CA2659467A1 (en) |
IL (1) | IL196694A (en) |
MX (1) | MX2009001136A (en) |
WO (1) | WO2008017076A2 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009300393A (en) * | 2008-06-17 | 2009-12-24 | Shimadzu Corp | X-ray tomogram photographing apparatus |
WO2010011676A3 (en) * | 2008-07-23 | 2010-04-22 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
US8270760B2 (en) | 2006-08-03 | 2012-09-18 | The Regents Of The University Of California | Iterative methods for dose reduction and image enhancement in tomography |
CN102831629A (en) * | 2012-08-23 | 2012-12-19 | 汕头市东方超声科技有限公司 | Graphic processor based mammary gland CT (Computerized Tomography) image reconstruction method |
US8442353B2 (en) | 2006-08-03 | 2013-05-14 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
US8532350B2 (en) | 2008-01-30 | 2013-09-10 | The Regents Of The University Of California | Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints |
US8611626B2 (en) | 2009-04-30 | 2013-12-17 | The Regents Of The University Of California | System and methods for fast implementation of equally-sloped tomography |
EP2771771A4 (en) * | 2011-10-27 | 2015-06-17 | Flatfrog Lab Ab | Touch determination by tomographic reconstruction |
WO2017087203A1 (en) * | 2015-11-19 | 2017-05-26 | General Electric Company | Water equivalent diameter determination from scout images |
US20220254538A1 (en) * | 2012-10-30 | 2022-08-11 | California Institute Of Technology | Fourier ptychographic imaging systems, devices, and methods |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101660945B (en) * | 2008-08-27 | 2013-02-20 | 中国科学院光电技术研究所 | Rapid image reconstruction method |
US8983808B2 (en) * | 2008-10-15 | 2015-03-17 | The United States Of America, As Represented By The Secretary Of The Navy | Automated mesh generation and editing tools |
US8655924B2 (en) * | 2008-10-15 | 2014-02-18 | The United States Of America, As Represented By The Secretary Of The Navy | System and method for providing structured data to a structured or unstructured grid |
DE102009048302B4 (en) * | 2009-10-05 | 2011-07-07 | Siemens Aktiengesellschaft, 80333 | Correction of truncations in MR imaging |
FR2963145B1 (en) * | 2010-07-26 | 2013-06-21 | Gen Electric | MONITORING THE DOSE OF RADIATION ACCUMULATED BY A BODY |
US8379948B2 (en) | 2010-12-21 | 2013-02-19 | General Electric Company | Methods and systems for fast iterative reconstruction using separable system models |
US8781243B2 (en) * | 2011-01-07 | 2014-07-15 | Wisconsin Alumni Research Foundation | Method for constrained reconstruction of high signal-to-noise ratio images |
JP2013075035A (en) * | 2011-09-30 | 2013-04-25 | Canon Inc | Ophthalmic apparatus, ophthalmic image processing method, and recording medium |
JP5774447B2 (en) | 2011-10-28 | 2015-09-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT apparatus, dose calculation method and program |
US9208588B2 (en) | 2013-09-25 | 2015-12-08 | Wisconsin Alumni Research Foundation | Fast statistical imaging reconstruction via denoised ordered-subset statistically-penalized algebraic reconstruction technique |
US10278666B2 (en) | 2014-06-13 | 2019-05-07 | Palodex Group Oy | Systems and methods of automated dose control in x-ray imaging |
CN104240272B (en) * | 2014-07-16 | 2017-03-15 | 中国人民解放军信息工程大学 | A kind of puppet polar coordinate TV that is based on minimizes straight path CT image rebuilding methods |
KR101727836B1 (en) * | 2015-12-22 | 2017-04-17 | 한국과학기술원 | Method and apparatus for restoring image |
US10134155B2 (en) * | 2016-11-21 | 2018-11-20 | Elekta Limited | Systems and methods for real-time imaging |
US10332252B2 (en) * | 2016-12-29 | 2019-06-25 | General Electric Company | Slope constrained cubic interpolation |
US11381737B2 (en) * | 2018-03-28 | 2022-07-05 | Sony Corporation | Arithmetic device and arithmetic method |
US11195310B2 (en) * | 2018-08-06 | 2021-12-07 | General Electric Company | Iterative image reconstruction framework |
CN109949267B (en) * | 2019-01-21 | 2021-01-29 | 西安交通大学 | Rapid three-dimensional reconstruction method for multi-material fault data to entity |
Family Cites Families (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4506327A (en) * | 1981-11-23 | 1985-03-19 | General Electric Company | Limited-angle imaging using multiple energy scanning |
JPS59149555A (en) * | 1983-02-15 | 1984-08-27 | Toshiba Corp | Picture reprocessing device |
JPS59168839A (en) * | 1983-03-17 | 1984-09-22 | 横河メディカルシステム株式会社 | Calculator tomography apparatus |
US4616318A (en) * | 1983-06-07 | 1986-10-07 | Elscint, Inc. | System for reprojecting images using transform techniques |
JPS60111640A (en) * | 1983-11-21 | 1985-06-18 | 横河メディカルシステム株式会社 | Image re-constituting system |
US4752879A (en) * | 1985-01-30 | 1988-06-21 | Picker International, Inc. | Method and apparatus for medical imaging |
US4888693A (en) * | 1987-04-01 | 1989-12-19 | General Electric Company | Method to obtain object boundary information in limited-angle computerized tomography |
JPH01141646A (en) * | 1987-11-30 | 1989-06-02 | Yokogawa Medical Syst Ltd | Fourier transform type x-ray ct |
US5375156A (en) * | 1992-03-31 | 1994-12-20 | Siemens Medical Systems, Inc. | Method and apparatus for 3-D computer tomography |
US6304317B1 (en) * | 1993-07-15 | 2001-10-16 | Nikon Corporation | Projection apparatus and method |
JP3373720B2 (en) * | 1996-03-25 | 2003-02-04 | 株式会社日立メディコ | X-ray tomography equipment |
US5778038A (en) * | 1996-06-06 | 1998-07-07 | Yeda Research And Development Co., Ltd. | Computerized tomography scanner and method of performing computerized tomography |
US5937102A (en) * | 1996-10-09 | 1999-08-10 | California Institute Of Technology | Image reconstruction |
US5761267A (en) * | 1996-12-26 | 1998-06-02 | General Electric Company | Methods and apparatus for simplified filtering of scan data in an imaging system |
US6748098B1 (en) * | 1998-04-14 | 2004-06-08 | General Electric Company | Algebraic reconstruction of images from non-equidistant data |
JP2000090249A (en) * | 1998-09-14 | 2000-03-31 | Fujitsu Ltd | Image processor, image processing method and computer readable recording medium recorded with program for operating computer as the image processor |
US6744848B2 (en) * | 2000-02-11 | 2004-06-01 | Brandeis University | Method and system for low-dose three-dimensional imaging of a scene |
FR2810141B1 (en) * | 2000-06-07 | 2002-08-23 | Commissariat Energie Atomique | METHOD FOR ACCELERATED RECONSTRUCTION OF A THREE-DIMENSIONAL IMAGE |
US6366638B1 (en) * | 2001-03-16 | 2002-04-02 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for CT scout image processing |
US20030028386A1 (en) * | 2001-04-02 | 2003-02-06 | Zinser Richard L. | Compressed domain universal transcoder |
US7206459B2 (en) * | 2001-07-31 | 2007-04-17 | Ricoh Co., Ltd. | Enhancement of compressed images |
US6771733B2 (en) * | 2001-08-16 | 2004-08-03 | University Of Central Florida | Method of reconstructing images for spiral and non-spiral computer tomography |
US6873744B2 (en) * | 2002-04-17 | 2005-03-29 | Regents Of The University Of Minnesota | Image restoration from transformed component data |
CN100443917C (en) * | 2001-12-17 | 2008-12-17 | 皇家飞利浦电子股份有限公司 | Magnetic resonance method |
JP3950716B2 (en) * | 2002-03-06 | 2007-08-01 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT apparatus and control method thereof |
US7620440B2 (en) * | 2002-05-17 | 2009-11-17 | Case Western Reserve University | Direct temporal encoding of spatial information |
JP4467873B2 (en) * | 2002-09-24 | 2010-05-26 | 株式会社日立メディコ | Tomography equipment |
US6728331B1 (en) * | 2002-10-21 | 2004-04-27 | Ge Medical Systems Global Technology Company, Llc | Method and system for trauma application of CT imaging |
US20040215072A1 (en) * | 2003-01-24 | 2004-10-28 | Quing Zhu | Method of medical imaging using combined near infrared diffusive light and ultrasound |
JP4138558B2 (en) * | 2003-04-03 | 2008-08-27 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Image reconstruction device, image reconstruction method, and radiation tomography apparatus |
US7209535B2 (en) * | 2003-06-20 | 2007-04-24 | Wisconsin Alumni Research Foundation | Fourier space tomographic image reconstruction method |
US6862337B2 (en) * | 2003-06-25 | 2005-03-01 | General Electric Company | Linear track based digital tomosynthesis system and method |
US6990167B2 (en) * | 2003-08-29 | 2006-01-24 | Wisconsin Alumni Research Foundation | Image reconstruction method for divergent beam scanner |
US7315636B2 (en) * | 2003-09-18 | 2008-01-01 | Accuray, Inc. | Generation of reconstructed images |
JP2005218693A (en) * | 2004-02-06 | 2005-08-18 | Ge Medical Systems Global Technology Co Llc | Image reconstruction method and x-ray ct apparatus |
CN101237820B (en) * | 2005-04-25 | 2011-01-26 | 罗切斯特大学 | Method and apparatus of global de-noising for CT imaging |
US7702141B2 (en) * | 2005-06-29 | 2010-04-20 | General Electric Company | Method for quantifying an object in a larger structure using a reconstructed image |
EP1958161B1 (en) | 2005-08-31 | 2013-01-30 | GE Healthcare Limited | Method and system of multivariate analysis on normalized volume-wise data in the sinogram domain for improved quality in positron emission tomography studies |
JP4795830B2 (en) * | 2006-03-30 | 2011-10-19 | テルモ株式会社 | Diagnostic imaging apparatus and processing method thereof |
JP4760564B2 (en) * | 2006-06-20 | 2011-08-31 | 日本電気株式会社 | Pattern shape defect detection method and detection apparatus |
US7439739B2 (en) * | 2006-07-11 | 2008-10-21 | The Board Of Trustees Of The Leland Stanford Junior University | Anti-aliased magnetic resonance image reconstruction using partially parallel encoded data |
AU2007281076A1 (en) | 2006-08-03 | 2008-02-07 | The Regents Of The University Of California | Iterative methods for dose reduction and image enhancement in tomography |
US8442353B2 (en) * | 2006-08-03 | 2013-05-14 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
KR20100133950A (en) * | 2008-01-30 | 2010-12-22 | 더 리전츠 오브 더 유니버시티 오브 캘리포니아 | Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints |
EP2424438A4 (en) * | 2009-04-30 | 2013-04-03 | Univ California | System and methods for fast implementation of equally-sloped tomography |
-
2007
- 2007-08-03 AU AU2007281076A patent/AU2007281076A1/en not_active Abandoned
- 2007-08-03 BR BRPI0715216-7A patent/BRPI0715216A2/en not_active IP Right Cessation
- 2007-08-03 JP JP2009523075A patent/JP2009545409A/en active Pending
- 2007-08-03 CA CA002659467A patent/CA2659467A1/en not_active Abandoned
- 2007-08-03 MX MX2009001136A patent/MX2009001136A/en active IP Right Grant
- 2007-08-03 WO PCT/US2007/075220 patent/WO2008017076A2/en active Application Filing
- 2007-08-03 KR KR1020097002178A patent/KR20090057972A/en not_active Application Discontinuation
- 2007-08-03 CN CNA2007800286810A patent/CN101495038A/en active Pending
- 2007-08-03 EP EP07840699A patent/EP2049019A4/en not_active Withdrawn
-
2009
- 2009-01-25 IL IL196694A patent/IL196694A/en active IP Right Grant
- 2009-01-30 US US12/363,079 patent/US8270760B2/en not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
None |
See also references of EP2049019A4 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8442353B2 (en) | 2006-08-03 | 2013-05-14 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
US8270760B2 (en) | 2006-08-03 | 2012-09-18 | The Regents Of The University Of California | Iterative methods for dose reduction and image enhancement in tomography |
US8532350B2 (en) | 2008-01-30 | 2013-09-10 | The Regents Of The University Of California | Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints |
JP2009300393A (en) * | 2008-06-17 | 2009-12-24 | Shimadzu Corp | X-ray tomogram photographing apparatus |
WO2010011676A3 (en) * | 2008-07-23 | 2010-04-22 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
EP2317925A2 (en) * | 2008-07-23 | 2011-05-11 | The Regents of the University of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
EP2317925A4 (en) * | 2008-07-23 | 2012-10-24 | Univ California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
US8611626B2 (en) | 2009-04-30 | 2013-12-17 | The Regents Of The University Of California | System and methods for fast implementation of equally-sloped tomography |
EP2771771A4 (en) * | 2011-10-27 | 2015-06-17 | Flatfrog Lab Ab | Touch determination by tomographic reconstruction |
CN102831629A (en) * | 2012-08-23 | 2012-12-19 | 汕头市东方超声科技有限公司 | Graphic processor based mammary gland CT (Computerized Tomography) image reconstruction method |
US20220254538A1 (en) * | 2012-10-30 | 2022-08-11 | California Institute Of Technology | Fourier ptychographic imaging systems, devices, and methods |
WO2017087203A1 (en) * | 2015-11-19 | 2017-05-26 | General Electric Company | Water equivalent diameter determination from scout images |
US9895130B2 (en) | 2015-11-19 | 2018-02-20 | General Electric Company | Water equivalent diameter determination from scout images |
Also Published As
Publication number | Publication date |
---|---|
EP2049019A2 (en) | 2009-04-22 |
AU2007281076A2 (en) | 2009-03-26 |
IL196694A0 (en) | 2009-11-18 |
CA2659467A1 (en) | 2008-02-07 |
US8270760B2 (en) | 2012-09-18 |
JP2009545409A (en) | 2009-12-24 |
IL196694A (en) | 2013-04-30 |
EP2049019A4 (en) | 2011-11-30 |
US20090232377A1 (en) | 2009-09-17 |
BRPI0715216A2 (en) | 2013-06-18 |
AU2007281076A1 (en) | 2008-02-07 |
WO2008017076A3 (en) | 2008-07-03 |
MX2009001136A (en) | 2009-02-10 |
KR20090057972A (en) | 2009-06-08 |
CN101495038A (en) | 2009-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8270760B2 (en) | Iterative methods for dose reduction and image enhancement in tomography | |
US8611626B2 (en) | System and methods for fast implementation of equally-sloped tomography | |
US20110164799A1 (en) | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography | |
US8532350B2 (en) | Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints | |
US8805037B2 (en) | Method and system for reconstruction of tomographic images | |
US20060002509A1 (en) | Method and apparatus for direct reconstuction in tomosynthesis imaging | |
Badea et al. | 4D micro-CT for cardiac and perfusion applications with view under sampling | |
US9619902B2 (en) | Filter for tomographic reconstructions | |
JP2008006288A (en) | System and method for iterative image reconstruction | |
EP2317925A2 (en) | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography | |
NL2008891C2 (en) | Method and system for reconstruction of tomographic images. | |
US20120183194A1 (en) | Method for artifact reduction in cone-beam ct images | |
WO2012056361A1 (en) | Apparatus and method for hybrid reconstruction of an object from projection data. | |
US20140016847A1 (en) | Multi-phase computed tomography image reconstruction | |
WO2011011684A2 (en) | Selection of optimal views for computed tomography reconstruction | |
WO2016132880A1 (en) | Arithmetic device, x-ray ct device, and image reconstruction method | |
US12045916B2 (en) | Stochastic backprojection for 3D image reconstruction | |
US11270477B2 (en) | Systems and methods for tailored image texture in iterative image reconstruction | |
Konovalov et al. | Spatial resolution of few-view computed tomography using algebraic reconstruction techniques | |
Zhu | Improved compressed sensing algorithm for sparse-view CT | |
Fahimian et al. | Radiation dose reduction in medical CT through equally sloped tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200780028681.0 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 07840699 Country of ref document: EP Kind code of ref document: A2 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2007840699 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2659467 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2009/001136 Country of ref document: MX |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009523075 Country of ref document: JP Ref document number: 2007281076 Country of ref document: AU Ref document number: 557/CHENP/2009 Country of ref document: IN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 1020097002178 Country of ref document: KR |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2007281076 Country of ref document: AU Date of ref document: 20070803 Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: RU |
|
ENP | Entry into the national phase |
Ref document number: PI0715216 Country of ref document: BR Kind code of ref document: A2 Effective date: 20090130 |