WO2002067223A2 - System and method for fast parallel cone-beam reconstruction using one or more microprocessors - Google Patents
System and method for fast parallel cone-beam reconstruction using one or more microprocessors Download PDFInfo
- Publication number
- WO2002067223A2 WO2002067223A2 PCT/US2002/004183 US0204183W WO02067223A2 WO 2002067223 A2 WO2002067223 A2 WO 2002067223A2 US 0204183 W US0204183 W US 0204183W WO 02067223 A2 WO02067223 A2 WO 02067223A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- processing unit
- data
- computer
- cone
- projection
- Prior art date
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/428—Real-time
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Definitions
- the present invention is directed to a system and method for cone-beam reconstruction in medical imaging or the like and more particularly to such a system and method implemented on one or more microprocessors.
- the present invention is also useful for nondestructive testing, single photon emission tomography and CT-based explosive detection, micro CT or micro cone beam volume CT, etc.
- Cone-beam reconstruction has attracted much attention in the medical imaging , community. Examples of cone-beam reconstruction are found in the commonly assigned
- CT computed tomography
- IR iterative reconstruction
- the filtered backprojection is more often discussed because it is accurate and amenable to fast implementation.
- the filtered backprojection can be implemented as an exact reconstruction method or as an approximate reconstruction method, both based on the Radon transform and/or the Fourier transform.
- the cone beam reconstruction process is time-consuming and needs a lot of computing operation. Currently, the cone beam reconstruction process is prohibitively long for clinical and other practical applications.
- GFLOPS gigaflops
- the "link” method has been extended to 3D cone-beam FBP; after rebinning the projection data in each row, the same method as in 2D can be applied to rebinning data, and data processing time can be brought down to O(N 3 logN) complexity for cone beam reconstruction.
- Another fast algorithm has been presented, using Fast Hierarchical Backprojection (FHBP) algorithms for 2D FBP, which address some of the shortcomings of existing fast algorithms.
- FHBP algorithms are based on a hierarchical decomposition of the Radon transform and need O(N logN) computing complexity for reconstruction.
- a customized backprojection hardware engine having parallelism and pipelining of various kinds can push the execution speed to the very limit.
- the hardware can be an FPGA based module or an ASIC module, a customized mask-programmable gate array, a cell-based IC and field programmable logic device or an add-in board with high speed RISC or DSP processors. Those boards usually use high-speed multi-port buffer memory or a DMA controller to increase data exchanging speed between boards. Some techniques, like vector computing and pre-interpolating projection data, are used with the customized engine to decrease reconstruction operation. Most of the customized hardware is built for 2D FBP reconstruction applications. No reconstruction engine-based a single or multiple microprocessors that is specially designed for fast cone beam reconstruction is commercially available.
- a multi-processor computer or a multi-computer system can be used to accelerate the cone beam reconstruction algorithm.
- Many large-scale parallel computers have tightly coupled processors interconnected by high-speed data paths.
- the multi-processor computer can be a shared memory computer or a distributed memory computer.
- Much work has been done on the large-scale and extremely expensive parallel computer.
- Most of that work uses an algorithm based on the 3D Radon transform.
- the Feldkamp algorithm and two iterative algorithms, 3D ART and SIRT have been implemented on large-scale computers such as Cray-3D, Paragon and SP1.
- the local data partition is used for the Feldkamp algorithm and the SIRT algorithm
- the global data partition is used for the ART algorithm.
- the implementation is voxel driven.
- processors The communication speed between processors is important to the reconstruction time, and the Feldkamp implementation can gain best performance in Multiple Instruction Multiple Data (MIMD) computers.
- MIMD Multiple Instruction Multiple Data
- Parallel 2D FBP has been implemented on Intel Paragon and CM5 computers. Using customized accelerating hardware or a large-scale parallel computer is not a cost-effective fast reconstruction solution, and it is not convenient to modify or add a new algorithm for research work.
- processors can gain data and operation parallelism with some micro-architecture techniques.
- Instruction-level Parallelism is a family of processor and compiler design techniques that speed up execution by causing individual machine operations to execute in parallel. Modern processors can divide instruction executing into several stages; some techniques such as pipeline and branch prediction permit the execution of multiple instructions simultaneously.
- SBVID single instruction multiple data
- Such processors include Intel's IA-32 architecture with MMXTM and SSE/SSE2, Motorola's PowerPCTM with AltNecTM and AMD Athlon with 3DnowTM.
- SBVID single instruction multiple data
- Such processors include Intel's IA-32 architecture with MMXTM and SSE/SSE2, Motorola's PowerPCTM with AltNecTM and AMD Athlon with 3DnowTM.
- the present invention is directed to a practical implementation for high-speed CBR on a commercially available PC based on hybrid computing (HC).
- Feldkamp CBR is implemented with multi-level acceleration, performing HC utilizing single instruction multiple data (SIMD) and making execution units (EU) in the processor work effectively.
- SIMD single instruction multiple data
- EU execution units
- the multi-thread and fiber support in the operating system can be exploited, which automatically enable the reconstruction parallelism in a multi-processor environment and also make data I/O to the hard disk more effective.
- Memory and cache access are optimized by proper data partitioning.
- the present invention can decrease filtering time by more than 75% for 288 projections each having 512 data points and can save more than 60% of the reconstruction time for 512 cube, while maintaining good precision with less than 0.08% average error.
- the resulting system is cost-effective and highspeed.
- An effective reconstruction engine can be built with a commercially available Symmetric Multi-processor (SMP) computer, which is easy and inexpensive to upgrade along with newer PC processors and memory with higher access speed.
- SMP Symmetric Multi-processor
- the Feldkamp algorithm cone beam reconstruction can achieve high speed with good precision.
- the test environment is an Intel Pentium III 500 Mhz with 640 MB 100 Mhz memory.
- the result shows that the reconstruction for a 512 3 cube with 288 projections can be finished in less than 20 minutes and maintains good precision, while the old implementation required more than 100 minutes.
- Several simulated phantoms have been used to test the precision of the HC FACBR. Comparing the reconstructed image with a simulated phantom image and images reconstructed by the traditional method shows less than a 0.04% average error compared to traditional method images and good precision to computer-simulated phantoms.
- a linear attenuation coefficient distribution of a three-dimensional object can be reconstructed quickly and accurately.
- a higher speed SSE-2 enabled Pentium IV and a 2- or 4-processor PC are expected to permit 512 3 cube FACBR in a few minutes in the future.
- FACBR is implemented with multilevel acceleration and hybrid computing utilizing the SIMD and ILP technology.
- the memory and cache access are optimized by proper data partition.
- the present invention is cost-effective and high-speed.
- a market available SMP computer provides an effective reconstruction engine which is easy and inexpensive to be upgraded along with newer PC processors. By contrast, custom built hardware is expensive and very difficult to upgrade.
- a high-speed implementation will be disclosed for FACBR on a PC. Techniques for hybrid execution (HE) and hybrid data (HD) will also be disclosed.
- a high speed Feldkamp implementation can be implemented on a general purpose PC with a high performance to price ratio.
- the HD and HE can also be applied to implementation on other hardware platforms to improve the FACBR perfo ⁇ nance.
- higher clock frequency processors and an inexpensive market available SMP PC it is possible to gain good performance as done by expensive, inconvenient customized hardware.
- As a commercial market available PC is used to achieve high performance, it is convenient to design new algorithms and a new system for cone beam reconstruction, and it is useful to integrate an image grab system and 3D rendering system, in a single system which is easy to configure and upgrade.
- the present invention implements parallel processing on a single microprocessor or multiple processors.
- the use of hybrid computing both fixed and floating point calculation) accelerates the cone-beam reconstruction without reducing the accuracy of the reconstruction and without increasing image noise. Those characteristics are particularly important for the reconstruction of soft tissue, e.g., cancer detection.
- Fig. 1 shows a cone-beam coordinate system used in reconstruction in the preferred embodiment
- Fig. 2 shows the architecture of an Intel 86 processor
- Fig. 3 shows an UML diagram of a known FACBR implementation
- Fig. 4 shows a UML diagram of hybrid execution according to the preferred embodiment
- Fig. 5 shows a data partition scheme used in the preferred embodiment
- Fig. 4 shows a UML diagram of hybrid execution according to the preferred embodiment
- Fig. 5 shows a data partition scheme used in the preferred embodiment
- Fig. 5 shows a data partition scheme used in the preferred embodiment
- FIG. 6 shows a block diagram of a system on which the preferred embodiment of the present invention can be implemented.
- the O-XYZ is the world coordinate system.
- the X-Y-Z axis gives the physical coordinates for the reconstructed voxels.
- the Z-axis is the rotation axis.
- the t-s axis is the rotated gantry X-Y coordinate system. The s-axis always passes through the x-ray source and is perpendicular to the detector plane.
- the Feldkamp algorithm falls into the class of filtered backprojection algorithms.
- the implementation of the Feldkamp algorithm contains following steps: a) Apply weight and ramp filter to the projections data; this is done by applying a weight to
- (t,s) is the coordinate in gantry system, which is the rotation transform of
- the reconstruction volume N x ⁇ N y ⁇ N z voxels in the x, y, and z directions.
- the Intel processor 200 has multiple execution units (EU's) to do integer and float operations simultaneously, hi an Intel Pentium III, there are two integer unit (ALUO 202 and ALU1 204), one float unit (FPU) 206, one MMX unit 208 to process 8-bit and 16-bit integers in parallel, and one Streaming SIMD Executing unit (SSE) to process four 32-bit single-precision float point data in parallel 210. Also present are an address generation unit 212, a memory load 214, a store address calculation unit 216, a memory store 218 and a reservation station 220.
- the SIMD instruction in the SSE enables four float integer operations at one instruction.
- the Pentium III processor has five pipelines 222 to exploit the parallelism in instruction execution. In each
- the processor core may dispatch zero or one ⁇ op on a port to any of the five
- Level one (LI) cache is the on-chip cache subsystem and consists of two 16-Kbyte four- ay set associative caches with a cache line length of 32bytes for instruction and data.
- the data cache has eight banks interleaved on four-byte boundaries.
- Level two (L2) cache is off-chip but in the same processor package. It usually has a size from 128Kbytes to 1Mbyte. L2 usually has a latency from 4 to 10 cycles for data access. When the processor needs to fetch instructions or data, LI is much faster than L2, and L2 is faster than access to main memory.
- the Instruction- Level Parallelism has two basic kinds: Individual instructions are overlapped (executed at the same time) in the processor, a given instruction is decomposed into sub operations and the sub operations are overlapped. As described in Feldkamp algorithm, a set of M projections is used, each projection having a size NxN pixels, to reconstruct anN cube. Each projection requires N loop calculations to do backprojection. projections require M*N loop calculations. Usually, M should be on the same level as Nto get a better result.
- the total actual reconstruction time can be written as:
- the reconstruction is finished, and the data are saved or rendered in a 3D display.
- the filtering time is about 1/15 to 1/30 of the backprojection time or less.
- Equation (2) shows that (u,t,s) depends only on (x,y), so that the projection map needs only O(N 3 ) computation time; thus both k and t unit . can be decreased.
- b) Use some a priori knowledge to generate some boundary as a sphere or cylinder; the computation can be skipped for some voxels which are outside the boundary and unable to be reconstructed, thereby providing a smaller k. If the reconstructed voxels are visualized as a cube with N length, then the full number of voxels is on the order of N 3 ,
- multi-processor computer works by multithread implementation and carefully allocates the tasks among the processors. Operating systems capable of controlling a multiprocessor computer in such a manner are known in the art, as noted above. For a single processor, the context switching will sacrifice the CPU time and so may actually decrease the performance, so it is contemplated that the multi-thread method will be used only when SMP is available.
- HD is used to decrease and make ALU units work in parallel.
- the SSE unit is independent from the FPU unit and the MMX unit, the SSE unit can work with the ALU unit, the MMX unit and even the FPU unit simultaneously, thus allowing a hybrid execute mode for either PF data or HD.
- the map data and some intermediate results can be processed by the ALU in fixed point data format, and the reconstruction data and finally output results can be processed in floating point format. That hybrid data format for different data and stages can improve the EU's efficiency.
- the best method is to use the MMX unit to adjust the data address and map data, and to use the SSE to do the backprojection calculation.
- the MMX can process data address and map data for two or more points, while the ALU can deal with only one point.
- the HE method for PF can be shown as a UML activity diagram in Figure 4. Since the MMX unit in a Pentium III processor can only process 8-bit and 16-bit integer multiplication, it is not so effective to do HE for HD data as to do HE for PF data. However, with new processor techniques such as SSE2 in the Pentium IN processor, the hybrid execution with HD will gain more improvements on speed.
- the second parallelism consideration is the data partition schema.
- the reconstruction data are partitioned into different sub-units.
- a data partition scheme is shown in Figure 5.
- Data are stored in memory as a one-dimensional array, in which the index of each data point increases for z, then for x and last for y.
- Data are processed in z-lines because the same projection data u value can be used for one z-line
- the projection data used to do backprojection for voxels in one z-line are actually in two adjacent w-lines, since four adjacent points (two in each of the adjacent w-lines) are used to interpolate the data for one voxel in the z-line.
- the reconstruction accuracy of the implementation will be determined using computer-simulated phantom.
- the reconstruction error noise level and uniformity of the reconstructed images are quantified using both pure float point implementation and HD computing implementation, and the reconstruction results from the two implementations are compared with both simulated phantom and experimental phantom data.
- the speedup of MC implementation is evaluated compared to normal pure float-point computing reconstruction.
- Experimental phantom data are also used to evaluate the effectiveness of the implementation in the real world.
- the Shepp Logan phantom is used as a general precision error compare reference.
- the cylinder phantom is used to compare the precision error at different z positions. Normally, the Feldkamp Algorithm has the best result at center slice, and the precision error increases for the slices at two ends.
- the cylinder phantom is used to check whether the HD and PF precision error varies with z-distance to center. Table 2
- the total FACBR time contains the filtering time and backprojection time; the filtering time takes only a small part in the total time.
- the backprojection process is the most time-consuming part of FACBR.
- the acceleration of five implementations has been tested; the results are shown in Table 4 below.
- the tests used 288 512 2 projections to reconstruct 512 3 data. All the reconstructions are run with a cylinder boundary.
- the program runs in Windows NT 4.0 and takes 95% to 98% of the processor time.
- the first column is the traditional PF method with boundary
- the second is PF with data partition
- the third one is HD method
- the forth is HD data with HE
- the fifth one is PF with SSE acceleration
- the last one is PF with HE.
- HD provides a 3 to 3.5 speed-up over the traditional implementation
- HE-HD provides a 4 to 5 speed-up, which is ahnost same as PF with SSE.
- This result declares two points: first, the HE-HD does not involve the SSE unit, so that cheaper processor like the Celeron can be used to get almost the same performance with SSE-PF; second, a higher speed-up can be obtained by using a functional unit which works with fixed-point data in the same way in which the SSE works with floating-point data; such functionality already appears in the Pentium IN processor.
- the HE-PF is the most efficient method in a Pentium III processor.
- Table 5 below shows the effective t umt of different slices for the HE-HD and HE-PF methods. Since the program runs in a multi-processor operating system, the processor time resource varies over time, so that the effective t U nit also varies over time. Basically, t un i t becomes stable as the slice number increases. When the slice number is less than 4 or the data are not 16-bytes aligned, the processor is unable to use SSE, and then the t unauti t is greater than when SSE is available. Therefore, the time for a single slice can be greater than for other slices. Table 5
- the error ratio will be calculated for each pixel. Then, the average error ratio will be calculated for the whole comparing Region of Interest (ROI).
- ROI Region of Interest
- HE-PF Since HE-PF works with floating-point data, it will not introduce an extra precision error compared to tradition PF method. The greatest concern is whether HD computing will bring more precision error or not. If the relative precision error between a PF reconstructed image and a phantom image is Epp, the relative error between a HD reconstructed image and a phantom image is E HD , and the relative error between a HD reconstructed image and a PF reconstructed image is E HP , the ratio of the hybrid computing error to the whole precision error is defined as:
- the precision errpr has been determined for a HD reconstructed image relative to a simulated phantom image and a PF reconstructed image.
- the precision error between the HD image and the PF image is less 0.03%; for the cylinder phantom, the E HP is less than 0.02%.
- the HD image keeps a good precision compared to the PF image.
- the E HP contributes less than 5% of the total error percentage to the simulated phantom image. This means that the algorithm introduces more than 95% of the total error.
- the HD images have enough precision and are comparable to PF images.
- FIG. 6 An apparatus on which the invention can be implemented is shown in Fig. 6, which is reproduced from Fig. 9 of the above-referenced U.S. Patent No. 5,999,587.
- a standard CT a 3-D reconstruction is obtained by stacking a series of slices.
- a volume CT a direct reconstruction of an object can be obtained.
- FIG. 6 it is shown how the cone-beam tomography system 600 of the present invention can be used to obtain a direct 3- D reconstruction of an object.
- the cone beam volume CT scanning apparatus 600 is illustrated in a simplified block diagram form. The invention may preferably be employed in conjunction with such a cone beam volume CT scanning apparatus to generate a 3-D reconstruction matrix of the object.
- a cone beam volume CT scanning apparatus examines a body P using a cone shaped radiation beam 604 which traverses a set of paths across the body.
- an x- ray source 610 and a 2-D detector 611 are mounted on a gantry frame 602 that rotates around the body P being examined.
- the operating voltage for the x-ray source is obtained from a conventional high- voltage generator 608 in such a manner that the x-ray source 610 produces the desired cone-shaped beam of radiation when the high- voltage is applied to it.
- the high- voltage generator 608 is energized by means of a power source 618, through a switch 616.
- a contrast solution injector 640 can be used as needed.
- a first motor 612 is also powered by the power source 618- such that it drives the gantry frame 602 in its orbit about the body, for example, in a clockwise direction as shown by the arrows adjacent to the frame.
- the power source 618 is turned on by means of switch 620 or other conventional control devices, in order to initiate a measurement sequence.
- a speed control circuit 614 is used to control the speed of rotation of the gantry frame 602 and to provide an output control signal which indicates when the speed of the motor 712 is at the desired level for taking measurements.
- the output from the rotational control 614 may also be utilized to operate the switch 616 such that the high- voltage generator 608 is only turned on when the gantry frame 602 is driven at the desired speed for making measurements.
- a tilt control 615 is
- the gantry frame tilt motor 613 means of the gantry frame tilt motor 613. That tilting allows the acquisition of arc projection data on the perpendicular arc. Such geometry results in a complete set of data for an object with a 25-40 cm diameter corresponding to a 37-60 cm field at the detectors 611 with a magnification of 1.5.
- the tilting of the gantry 602 is generally available in a standard CT gantry, to acquire arc projections, the minimal modification of a standard CT gantry has to be made such that the tilting of the gantry, the x-ray exposure timing and the projection acquisition are synchronized by the system control computer 624 as shown in FIG. 6.
- the circle-plus-arc geometry can be implemented in one of the following two ways. In the first
- the gantry 602 is tilted to a small angle ( ⁇ 15° to ⁇ 30°)
- the gantry 602 will be tilted toward -15° to -30°. Another half set of arc projections will be acquired only when the x-ray tube 610 and the 2-D detector 611 are at
- the second alternative method is to mechanically modify a standard CT gantry such that two short arc orbits are added to the gantry, and the x-ray tube 610 and the 2-D detector 611 can be moved on the arc to acquire the arc projections and on the circle to acquire the circle projections.
- One arc constitutes the orbit of the x-ray tube 610 and the other arc is the
- the two arc orbits are mounted 180° apart from each other.
- x-ray tube 610 and the 2-D detector 611 are synchronously moved on the arc orbits to acquire arc projections. Then, the x-ray tube 610 and the 2-D detector 611 are rotated on the gantry to acquire circle projections.
- a 2-D detector 611 Mounted on the gantry frame 602 opposite the x-ray source 610 is a 2-D detector 611 which has a dynamic range equal to or greater than 1000:1 and an image lag of less than 10%, for example a selenium thin film transistor (STFT) array or a silicon STFT array, in order to provide 2-D projections that correspond to an x-ray attenuation signal pattern.
- STFT selenium thin film transistor
- the x- ray source 610 and the 2-D detector 611 are mounted on the gantry frame 602 in such a manner that they both move synchronously.
- the cone-shaped beam of radiation 604 generated by the x-ray source 610 is projected through the body or object under test.
- the 2-D detector cone measures the radiation transmitted along the set of beam paths across the cone.
- a continuous series of two-dimensional detectors can be fixedly mounted proximate to the gantry frame 602 and the x-ray source 610 is mounted to the gantry frame such that, upon rotation of the gantry frame, the cone-shaped radiation beam 604 is projected through the body P under test and sequentially received by each of the series of detectors.
- a 2-D projection acquisition control and A/D conversion unit 626 under control of the scanning pulses sequentially obtained from the system control computer 624, which includes the clock 622, receives a sequence of outputs corresponding to different lines of the 2-D detector 611.
- Each line of the 2-D detector consists of many detection cells (at least • >100).
- the output of each detector cell represents a line integral of attenuation values measurable along one of the respective beam paths.
- the cone-shaped beam 604 subtends a cone angle sufficient to include the entire region of interest of the body.
- the analog-to-digital conversion unit 626 serves to digitize the projection signals and to save them in the 3-D image reconstruction array processor 628 and storage device 630.
- the method employed by the 3-D image reconstruction array processor 628 is the invented algorithm and method described in this application.
- the 3-D image reconstruction array processor 628 serves to transform the digitized projection signals into x-ray attenuation data vectors.
- the x-ray attenuation data matrix corresponds to x-ray attenuation at spaced grid locations within the body trunk being examined. Each data element of the matrix represents an x-ray attenuation value and the location of the element corresponds to a respective 3-D grid location within the body.
- a display processor 632 obtains the data stored as 3-D x-ray attenuation signal patterns in the memory storage 630, processes that data as previously described, and then the desired 3-D images are displayed on a 3-D display device 634.
- the 3-D image reconstruction array processor 632 may, for example, be a computer as described above with one or more Intel or Intel- compatible 86-class microprocessors. However, any processor or processors capable of the same or substantially the same parallel operation can be used.
- the present invention specific to x86 processors; instead, the invention can be used with any processor capable of implementing the algorithms described above and has particular utility with any processor that has a floating-point unit that can process more than one single- precision 32-bit datum within one instruction set and a fixed-point unit that can process more than one 16- or 32-bit data within one instruction set. Therefore, the present invention should be construed as limited only by the appended claims.
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2002251922A AU2002251922B2 (en) | 2001-02-16 | 2002-02-13 | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
CA002438387A CA2438387A1 (en) | 2001-02-16 | 2002-02-13 | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
EP02720960.0A EP1366469B1 (en) | 2001-02-16 | 2002-02-13 | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/784,331 US6477221B1 (en) | 2001-02-16 | 2001-02-16 | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
US09/784,331 | 2001-02-16 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2002067223A2 true WO2002067223A2 (en) | 2002-08-29 |
WO2002067223A3 WO2002067223A3 (en) | 2002-10-24 |
Family
ID=25132101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2002/004183 WO2002067223A2 (en) | 2001-02-16 | 2002-02-13 | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
Country Status (6)
Country | Link |
---|---|
US (1) | US6477221B1 (en) |
EP (1) | EP1366469B1 (en) |
CN (1) | CN1284122C (en) |
AU (1) | AU2002251922B2 (en) |
CA (1) | CA2438387A1 (en) |
WO (1) | WO2002067223A2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1579372A2 (en) * | 2002-12-02 | 2005-09-28 | Visiongate, Inc. | Method and apparatus for three-dimensional imaging in the fourier domain |
CN102877828A (en) * | 2012-09-09 | 2013-01-16 | 山西山地物探技术有限公司 | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land |
Families Citing this family (50)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6987831B2 (en) | 1999-11-18 | 2006-01-17 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US20020169680A1 (en) * | 2001-05-10 | 2002-11-14 | International Business Machines Corporation | Method and apparatus for building commercial distributed computing networks via computer cost subsidization |
US6771733B2 (en) * | 2001-08-16 | 2004-08-03 | University Of Central Florida | Method of reconstructing images for spiral and non-spiral computer tomography |
US6638226B2 (en) | 2001-09-28 | 2003-10-28 | Teratech Corporation | Ultrasound imaging system |
JP3870105B2 (en) * | 2002-02-22 | 2007-01-17 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Back projection method and X-ray CT apparatus |
AU2003301341A1 (en) * | 2002-10-15 | 2004-05-04 | Digitome Corporation | Ray tracing kernel |
US6904117B2 (en) * | 2002-10-30 | 2005-06-07 | Toshiba Corporation | Tilted gantry helical cone-beam Feldkamp reconstruction for multislice CT |
DE10304662A1 (en) * | 2003-02-05 | 2004-08-19 | Siemens Ag | Method for generating images in computer tomography using a 3D image reconstruction method |
EP1627361A1 (en) * | 2003-05-27 | 2006-02-22 | Clean Earth Technologies, LLC | Method for fast image reconstruction with compact radiation source and detector arrangement using computerized tomography |
US7134036B1 (en) * | 2003-12-12 | 2006-11-07 | Sun Microsystems, Inc. | Processor core clock generation circuits |
US7362843B2 (en) * | 2004-09-23 | 2008-04-22 | General Electric Company | System and method for reconstruction of cone beam tomographic projections with missing data |
WO2006116316A2 (en) | 2005-04-22 | 2006-11-02 | University Of Chicago | Open source trajectory method and apparatus for interior imaging |
CA2608119A1 (en) | 2005-05-11 | 2006-11-16 | Optosecurity Inc. | Method and system for screening luggage items, cargo containers or persons |
US7991242B2 (en) | 2005-05-11 | 2011-08-02 | Optosecurity Inc. | Apparatus, method and system for screening receptacles and persons, having image distortion correction functionality |
US7492858B2 (en) | 2005-05-20 | 2009-02-17 | Varian Medical Systems, Inc. | System and method for imaging and treatment of tumorous tissue in breasts using computed tomography and radiotherapy |
US7646842B2 (en) * | 2005-09-23 | 2010-01-12 | General Electric Company | Methods and apparatus for reconstructing thick image slices |
US20070132754A1 (en) * | 2005-12-12 | 2007-06-14 | Intel Corporation | Method and apparatus for binary image classification and segmentation |
CN101297325B (en) * | 2005-12-29 | 2013-04-24 | 英特尔公司 | Method and device for radio tracking |
WO2007098284A2 (en) * | 2006-02-27 | 2007-08-30 | University Of Rochester | Method and apparatus for cone beam ct dynamic imaging |
AU2007221086A1 (en) * | 2006-02-27 | 2007-09-07 | University Of Rochester | Phase contrast cone-beam CT imaging |
WO2007126932A1 (en) * | 2006-03-28 | 2007-11-08 | Xoran Technologies, Inc. | Ct scanner with automatic determination of volume of interest |
WO2007124338A1 (en) * | 2006-04-19 | 2007-11-01 | Xoran Technologies, Inc. | Ct scanner with untracked markers |
DE102007020879A1 (en) | 2006-05-10 | 2009-04-02 | Gachon University Of Medicine & Science Industry-Academic Cooperation Foundation | Three dimensional image e.g. positron emission tomography image, reconstructing method, involves back-projecting sinogram data with projection angles for producing image data using symmetry properties, and projecting image data |
US7899232B2 (en) | 2006-05-11 | 2011-03-01 | Optosecurity Inc. | Method and apparatus for providing threat image projection (TIP) in a luggage screening system, and luggage screening system implementing same |
CN100386779C (en) * | 2006-06-02 | 2008-05-07 | 清华大学 | Orthographic projection and back projection methods for body to be measured based on general image display card |
US8494210B2 (en) | 2007-03-30 | 2013-07-23 | Optosecurity Inc. | User interface for use in security screening providing image enhancement capabilities and apparatus for implementing same |
DE102006036327A1 (en) * | 2006-08-03 | 2008-02-14 | Siemens Ag | Method for providing 3D image data and system for taking X-ray images |
US8217937B2 (en) * | 2007-03-28 | 2012-07-10 | The Aerospace Corporation | Isosurfacial three-dimensional imaging system and method |
GB2465726A (en) | 2007-08-23 | 2010-06-02 | Fischer Medical Technologies Inc | Improved computed tomography breast imaging and biopsy system |
US8023767B1 (en) | 2008-03-10 | 2011-09-20 | University Of Rochester | Method and apparatus for 3D metal and high-density artifact correction for cone-beam and fan-beam CT imaging |
US7940891B2 (en) | 2008-10-22 | 2011-05-10 | Varian Medical Systems, Inc. | Methods and systems for treating breast cancer using external beam radiation |
US8285971B2 (en) * | 2008-12-16 | 2012-10-09 | International Business Machines Corporation | Block driven computation with an address generation accelerator |
US8327345B2 (en) * | 2008-12-16 | 2012-12-04 | International Business Machines Corporation | Computation table for block computation |
US8458439B2 (en) * | 2008-12-16 | 2013-06-04 | International Business Machines Corporation | Block driven computation using a caching policy specified in an operand data structure |
US8407680B2 (en) * | 2008-12-16 | 2013-03-26 | International Business Machines Corporation | Operand data structure for block computation |
US8281106B2 (en) * | 2008-12-16 | 2012-10-02 | International Business Machines Corporation | Specifying an addressing relationship in an operand data structure |
WO2010093357A1 (en) * | 2009-02-11 | 2010-08-19 | Tomotherapy Incorporated | Target pedestal assembly and method of preserving the target |
US7949095B2 (en) | 2009-03-02 | 2011-05-24 | University Of Rochester | Methods and apparatus for differential phase-contrast fan beam CT, cone-beam CT and hybrid cone-beam CT |
CN101664583B (en) * | 2009-09-09 | 2012-05-09 | 深圳市海博科技有限公司 | Dosage calculation optimization method and system based on CUDA |
CA2849398C (en) | 2011-09-07 | 2020-12-29 | Rapiscan Systems, Inc. | X-ray inspection system that integrates manifest data with imaging/detection processing |
CN104540451B (en) | 2012-03-05 | 2019-03-08 | 罗切斯特大学 | Method and apparatus for differential phase contrast Cone-Beam CT and mixing Cone-Beam CT |
US9364191B2 (en) | 2013-02-11 | 2016-06-14 | University Of Rochester | Method and apparatus of spectral differential phase-contrast cone-beam CT and hybrid cone-beam CT |
CN105027227B (en) | 2013-02-26 | 2017-09-08 | 安科锐公司 | Electromagnetically actuated multi-diaphragm collimator |
KR102075191B1 (en) * | 2013-06-05 | 2020-02-07 | 에베 그룹 에. 탈너 게엠베하 | Measuring device and method for ascertaining a pressure map |
CN105806858B (en) * | 2014-12-31 | 2019-05-17 | 北京固鸿科技有限公司 | CT detection method and CT equipment |
WO2017146930A1 (en) | 2016-02-22 | 2017-08-31 | Rapiscan Systems, Inc. | Systems and methods for detecting threats and contraband in cargo |
US10482632B2 (en) * | 2017-04-28 | 2019-11-19 | Uih America, Inc. | System and method for image reconstruction |
CN109157215B (en) * | 2018-08-29 | 2021-09-28 | 中国医学科学院生物医学工程研究所 | Magnetic induction magnetoacoustic conductivity image reconstruction method based on system matrix |
US20230175987A1 (en) * | 2020-11-18 | 2023-06-08 | Jed Co., Ltd | X-ray inspection apparatus |
CN116543088B (en) * | 2023-07-07 | 2023-09-19 | 有方(合肥)医疗科技有限公司 | CBCT image reconstruction method and device |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5777620A (en) * | 1992-10-02 | 1998-07-07 | Canon Kabushiki Kaisha | 3D graphics system grouping surface primitives with and without specularity |
US6049343A (en) * | 1997-01-20 | 2000-04-11 | Hitachi, Ltd. | Graphics processing unit and graphics processing system |
US6078638A (en) * | 1998-09-30 | 2000-06-20 | Siemens Corporate Research, Inc. | Pixel grouping for filtering cone beam detector data during 3D image reconstruction |
EP1061476A2 (en) * | 1999-06-18 | 2000-12-20 | Marconi Medical Systems, Inc. | Volumetric image reconstruction |
EP1077430A2 (en) * | 1999-08-16 | 2001-02-21 | Analogic Corporation | Apparatus and method for reconstruction of volumetric images in a computed tomography system |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5253171A (en) * | 1990-09-21 | 1993-10-12 | General Electric Company | Parallel processing method and apparatus based on the algebra reconstruction technique for reconstructing a three-dimensional computerized tomography (CT) image from cone beam projection data |
US5257183A (en) | 1990-12-21 | 1993-10-26 | General Electric Company | Method and apparatus for converting cone beam X-ray projection data to planar integral and reconstructing a three-dimensional computerized tomography (CT) image of an object |
US5170439A (en) | 1991-06-11 | 1992-12-08 | Picker International, Inc. | Cone beam reconstruction using combined circle and line orbits |
US5365560A (en) | 1991-07-29 | 1994-11-15 | General Electric Company | Method and apparatus for acquiring a uniform distribution of radon data sufficiently dense to constitute a complete set for exact image reconstruction of an object irradiated by a cone beam source |
US5333164A (en) * | 1991-12-11 | 1994-07-26 | General Electric Company | Method and apparatus for acquiring and processing only a necessary volume of radon data consistent with the overall shape of the object for efficient three dimensional image reconstruction |
US5390226A (en) | 1992-07-02 | 1995-02-14 | General Electric Company | Method and apparatus for pre-processing cone beam projection data for exact three dimensional computer tomographic image reconstruction of a portion of an object |
US5517602A (en) | 1992-12-03 | 1996-05-14 | Hewlett-Packard Company | Method and apparatus for generating a topologically consistent visual representation of a three dimensional surface |
US5278884A (en) | 1992-12-18 | 1994-01-11 | General Electric Company | Complete 3D CT data acquisition using practical scanning paths on the surface of a sphere |
US5461650A (en) | 1993-10-18 | 1995-10-24 | General Electric Company | Method and system for pre-processing cone beam data for reconstructing free of interpolation-induced artifacts a three dimensional computerized tomography image |
US5400255A (en) | 1994-02-14 | 1995-03-21 | General Electric Company | Reconstruction of images from cone beam data |
US6002738A (en) * | 1995-07-07 | 1999-12-14 | Silicon Graphics, Inc. | System and method of performing tomographic reconstruction and volume rendering using texture mapping |
US5671265A (en) | 1995-07-14 | 1997-09-23 | Siemens Corporate Research, Inc. | Evidential reconstruction of vessel trees from X-ray angiograms with a dynamic contrast bolus |
JPH09149902A (en) | 1995-12-01 | 1997-06-10 | Hitachi Medical Corp | Tomography and tomograph |
-
2001
- 2001-02-16 US US09/784,331 patent/US6477221B1/en not_active Expired - Lifetime
-
2002
- 2002-02-13 AU AU2002251922A patent/AU2002251922B2/en not_active Ceased
- 2002-02-13 EP EP02720960.0A patent/EP1366469B1/en not_active Expired - Lifetime
- 2002-02-13 CA CA002438387A patent/CA2438387A1/en not_active Abandoned
- 2002-02-13 WO PCT/US2002/004183 patent/WO2002067223A2/en not_active Application Discontinuation
- 2002-02-13 CN CNB028050894A patent/CN1284122C/en not_active Expired - Lifetime
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5777620A (en) * | 1992-10-02 | 1998-07-07 | Canon Kabushiki Kaisha | 3D graphics system grouping surface primitives with and without specularity |
US6049343A (en) * | 1997-01-20 | 2000-04-11 | Hitachi, Ltd. | Graphics processing unit and graphics processing system |
US6078638A (en) * | 1998-09-30 | 2000-06-20 | Siemens Corporate Research, Inc. | Pixel grouping for filtering cone beam detector data during 3D image reconstruction |
EP1061476A2 (en) * | 1999-06-18 | 2000-12-20 | Marconi Medical Systems, Inc. | Volumetric image reconstruction |
EP1077430A2 (en) * | 1999-08-16 | 2001-02-21 | Analogic Corporation | Apparatus and method for reconstruction of volumetric images in a computed tomography system |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1579372A2 (en) * | 2002-12-02 | 2005-09-28 | Visiongate, Inc. | Method and apparatus for three-dimensional imaging in the fourier domain |
EP1579372A4 (en) * | 2002-12-02 | 2006-09-20 | Visiongate Inc | Method and apparatus for three-dimensional imaging in the fourier domain |
CN102877828A (en) * | 2012-09-09 | 2013-01-16 | 山西山地物探技术有限公司 | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land |
Also Published As
Publication number | Publication date |
---|---|
CA2438387A1 (en) | 2002-08-29 |
US6477221B1 (en) | 2002-11-05 |
AU2002251922B2 (en) | 2008-01-03 |
EP1366469A2 (en) | 2003-12-03 |
EP1366469B1 (en) | 2017-10-04 |
WO2002067223A3 (en) | 2002-10-24 |
US20020154727A1 (en) | 2002-10-24 |
CN1491404A (en) | 2004-04-21 |
CN1284122C (en) | 2006-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6477221B1 (en) | System and method for fast parallel cone-beam reconstruction using one or more microprocessors | |
AU2002251922A1 (en) | System and method for fast parallel cone-beam reconstruction using one or more microprocessors | |
Myagotin et al. | Efficient volume reconstruction for parallel-beam computed laminography by filtered backprojection on multi-core clusters | |
Scherl et al. | Evaluation of state-of-the-art hardware architectures for fast cone-beam CT reconstruction | |
US20080095300A1 (en) | System and method for iterative reconstruction using parallel processing | |
US5901196A (en) | Reduction of hitlist size in spiral cone beam CT by use of local radon origins | |
US7209535B2 (en) | Fourier space tomographic image reconstruction method | |
WO2005071601A1 (en) | Improved methods and apparatus for back-projection and forward-projection | |
US20080085040A1 (en) | System and method for iterative reconstruction using mask images | |
Keck et al. | GPU-accelerated SART reconstruction using the CUDA programming environment | |
Chen et al. | A hybrid architecture for compressive sensing 3-D CT reconstruction | |
Scherl et al. | Implementation of the FDK algorithm for cone-beam CT on the cell broadband engine architecture | |
Lu et al. | Cache-aware GPU optimization for out-of-core cone beam CT reconstruction of high-resolution volumes | |
Kachelrieß et al. | Hyperfast perspective cone--beam backprojection | |
Li | Design of an FPGA-based computing platform for real-time three-dimensional medical imaging | |
Goddard et al. | High-speed cone-beam reconstruction: An embedded systems approach | |
Zeng et al. | A fast CT reconstruction scheme for a general multi-core PC | |
Yu et al. | High-speed cone-beam reconstruction on PC | |
US20060039525A1 (en) | Method and apparatus for exact cone beam computed tomography | |
Yang et al. | Parallel implementation of Katsevich's FBP algorithm | |
Zheng et al. | A distributed multi-GPU system for high speed electron microscopic tomographic reconstruction | |
Okitsu et al. | Accelerating cone beam reconstruction using the CUDA-enabled GPU | |
Xu et al. | Mapping iterative medical imaging algorithm on cell accelerator | |
Kingswood et al. | Image reconstruction using the transputer | |
Shih et al. | Fast algorithm for X-ray cone-beam microtomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 02805089.4 Country of ref document: CN |
|
AK | Designated states |
Kind code of ref document: A2 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG UZ VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A2 Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
AK | Designated states |
Kind code of ref document: A3 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG UZ VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A3 Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2002251922 Country of ref document: AU |
|
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2438387 Country of ref document: CA |
|
REEP | Request for entry into the european phase |
Ref document number: 2002720960 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2002720960 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 2002720960 Country of ref document: EP |
|
REG | Reference to national code |
Ref country code: DE Ref legal event code: 8642 |
|
NENP | Non-entry into the national phase |
Ref country code: JP |
|
WWW | Wipo information: withdrawn in national office |
Ref document number: JP |
|
ENP | Entry into the national phase |
Ref document number: 2002251922 Country of ref document: AU Date of ref document: 20020213 Kind code of ref document: B |