This invention relates generally to processing and image reconstruction based on acquired raw image data. In particular, the present invention relates to increasing the processing performance with respect to the image reconstruction of diagnostic image data.
Raw image data from various diagnostic medical systems, such as Computed Tomography (CT) and Positron Emission Tomography (PET) systems, is acquired for diagnostic purposes. The CT and PET systems need to be able to support numerous scanning and reconstruction modes. The associated reconstruction algorithms are complex and computationally intensive. Users of the diagnostic medical systems desire an improvement in image quality, along with minimizing the time required to generate images based on raw image data and improving the reliability of the reconstruction process. By decreasing the amount of time needed to generate the desired images from raw image data, images can be evaluated sooner and patient through-put may be improved.
Previous diagnostic systems used different specialized processing units to accomplish specific tasks. That is, the reconstruction process was broken down according to the steps to be done. The processing units may operate in parallel or serially with respect to each other. In order to add processing capability, however, new processing units had to be added and the system design reconfigured and/or coordinated to integrate the new units and steps, both increasing the complexity and limiting the flexibility of the diagnostic system. Scalability and increased performance are thus difficult to achieve when adding additional processing units.
Thus, an apparatus and method are desired to reconstruct image data that addresses the problems noted above and others previously experienced.
In one embodiment, a method for generating images. Raw image data representative of an object of interest is acquired. The raw image data is decomposed into N subsets of raw image data. N is based on a number of available image generation processors. The N subsets of raw image data are processed to create processed image data. The image generation processors perform image processing on the image data in parallel with respect to each other.
In another embodiment, a method for increasing the performance of a system for processing raw image data. Raw image data is acquired which is representative of an object of interest. The raw image data is divided into substantially equal subsets of raw image data. At least one of the substantially equal subsets of raw image data is processed with an image generation processor.
In another embodiment, a scalable apparatus for processing raw image data. A data decomposer divides raw image data which is acquired by a data acquisition system. At least two image generation processors process the raw image data in parallel with respect to each other.
The x-ray source 12 and the two-dimensional detector array 20 are mounted on either side of a gantry 24 so as to rotate about an axis of rotation 26 generally positioned within the patient 18. The axis of rotation 26 forms the z-axis of a Cartesian coordinate system having its origin centered within the cone beam 14. The plane defined by the x and y axis of this coordinate system thus defines a plane of rotation, specifically the gantry plane 28 of the gantry 24.
The computer 38 is a general purpose minicomputer programmed to acquire and manipulate projection data. The computer 38 is connected to a data decomposer 46 and one or more image generation (IG) processors 40. The data decomposer 46 sends raw image data to the IG processors 40 which process raw image data as discussed below. In
The computer 38 receives commands and scanning parameters via operator console 42 which is generally a CRT display and keyboard that enables an operator to enter parameters for the CT scan and to display the reconstructed image. A mass storage device 44 provides a means for storing operating programs.
During data acquisition, the CT imaging system 10 functions as a conventional cone-beam system in gathering data. In the step-and-shoot acquisition method, the table 39 is held stationary as the x-ray emitter 12 and detector array 20 make a complete revolution around the gantry 24 about the axis of rotation 26. At each of a plurality of angular positions β, the attenuation data from all the detectors 22 in array 20 are stored in the mass memory 44. Upon completion of a full rotation, the computer 38 commands the table motor control 37 to advance the table 39 to another position along the z-axis 26 and another rotational scan of the patient 18 is performed. This process is repeated until the desired portion of the patient 18 has been fully scanned. Alternatively, the CT imaging system 10 may acquire data in the helical acquisition mode, wherein the table motor control 37 advances the table 39 as the x-ray emitter 12 and detector array 20 are rotated and scan data is being acquired.
The PET scanner system 100 includes a gantry 110 which supports a detector ring assembly 112 about a central bore which defines an imaging area 114. A patient table (not illustrated) is positioned in front of the gantry 110 and is aligned with the imaging area 114. A patient table controller (not shown) moves a table bed into the imaging area 114 in response to commands received from the operator work station 104 through a serial communications link 116.
A gantry controller 118 is mounted within the gantry 110 and is responsive to commands received from the operator work station 104 through the communication link 116 to operate the gantry 110. For example, the gantry 110 may perform a “coincidence timing calibration scan” to acquire corrective data, or an “emission scan” in which positron annihilation events are counted.
Returning to
Image reconstruction components 108 includes a sorter 134, a memory module 136, one or more IG processors 138, a data decomposer 162, an image CPU 140 and a backplane bus 142 which conforms to the VME standards and links all other processor components together. In
During a data acquisition, the coincidence counts are organized in memory 136 as a set of two-dimensional arrays, one for each axial image, and each having as one of its dimensions the projection angle θ and the other dimension the distance R. The θ by R array of detected events is called a histogram. Coincidence events occur at random and the sorter 134 quickly determines the θ and R values from the two crystal addresses in each coincidence data packet and increments the count of the corresponding sinogram array element. At the completion of an acquisition period, memory 136 stores the total number of annihilation events which occurred along each ray (R, θ) in the histogram.
Image CPU 140 controls the backplane bus 142 and links the image reconstruction components 108 to the communication link 116. The IG processors 138 also connect to the bus 142 and operate under the direction of the image CPU 140. The data decomposer 162 sends image data to the IG processors 138 which process the rash image data, or the histogram data, from the memory module 136 as discussed below. The resulting image array may be stored in memory module 136 and is output by the image CPU 140 to the operator work station 104.
The operator work station 104 includes a CPU 144, a CRT display 146 and a keyboard 148. The CPU 144 connects to the communications link 116 and scans the key board 148 for input information. Through the keyboard 148 and associated control panel switches, an operator can control calibration of the PET scanner system 100, its configuration, and the positioning of a patient table (not illustrated) during data acquisition.
Turning to
In step 202, the raw image data is sent to a data decomposer 250. The data decomposer 250 may comprise hardware and/or software. Therefore, the data decomposer 250 may be a separate component, or may be included within the computer 38 and mass storage 44 (
In step 204, the data decomposer 250 decomposes, or divides, the raw image data. The raw image data may be decomposed into three substantially equal portions, or subsets of image data. Continuing the example above wherein image data to form 900 data frames is acquired, each of three subsets of raw image data may comprise 300 data frames. It should be understood that the raw image data may not be decomposed based on data frames, and that the example is for illustrative purposes only. The first subset may comprise the first 300 data frames, the second subset may comprise the next 300 data frames, while the third subset may comprise the last 300 data frames. In another example, if the raw image data comprises data information for 1000 frames of data, the first and second subsets may comprise 333 data frames apiece, while the third subset comprises 334 data frames. The data decomposer 250 decomposes the raw image data into a number N of subsets based on the number of available IG processors 252-256. In other words, the data decomposer 250 determines how many IG processors 252-256 are available and automatically reconfigures itself to utilize the available IG processors 252-256.
In step 206, the data decomposer 250 sends the subsets of raw image data to the IG processors 252-256. For example, the first subset of raw image data may be sent to the IG processor 252, the second subset of raw image data may be sent to the IG processor 254, and the third subset of raw image data may be sent to the IG processor 256. The IG processors 252-256 each receive the respective subset of raw image data and begin the image generation process.
By way of example only, the data decomposer 250 may identify a predefined amount of raw image data for each IG processor 252-256. Therefore, once the image acquisition has been started and a portion of the raw image data has been acquired and received by the data decomposer 250, a predefined amount of raw image data is sent to each IG processor 252-256. The data decomposer 250 receives and holds raw image data as it is acquired. Then, when an IG processor 252-256 becomes available, the data decomposer 250 sends an amount of raw image data to the available IG processor 252-256. For example, the data decomposer 250 may send an amount of raw image data substantially equivalent to half of the total raw image data currently acquired and waiting to be processed. When the next IG processor 252-256 becomes available, the data decomposer 250 may again divide the raw image data currently acquired and waiting to be processed into two substantially equal subsets and send one subset to the next available IG processor 252-256 for processing. In this manner, the IG processors 252-256 are each engaged with raw image data to process, and the time any one IG processor 252-256 is idle is minimized.
Alternatively, the data decomposer 250 may decompose the raw image data into a number N of subsets where N is less than the number of total IG processors 252-256. For example, if raw image data from a previous scan remains to be processed, the data decomposer 250 may designate a portion of the total number of IG processors 252-256 to continue processing the previous scan, and divide the raw image data from the current scan between the remaining IG processors 252-256.
The IG processors 252-256 reconstruct the raw image data based on a reconstruction mode, which defines how the raw image data is processed. The reconstruction mode may be determined by the type of scan being performed, the anatomy being scanned, the desired image output, and the like. Many different reconstruction modes exist, and each mode comprises processing steps to be accomplished, optionally in a defined order, which may be different with respect to other reconstruction modes. Therefore, the steps 208-220 of
The IG processors 252-256 each process data independently with respect to each other. Thus, the IG processors 252-256 do not interact with each other, but rather operate in parallel. Also, as each of the IG processors 252-256 comprise the functional capability to process the raw image data independently, it is only necessary to have one IG processor 252-256. The IG processor 254 will be discussed below as an exemplary IG processor.
By way of example only, each of the IG processors 252-256 may comprise a “personal computer” or PC, having a motherboard, processors, and memory. Software processes are loaded in the memory and the processors process the raw image data according to the software processes. Therefore, each of the IG processors 252-256 may be a module, circuit board, or unit which is easily installed within the CT imaging system 10 and PET scanner system 100. As each of the IG processors 252-256 within an imaging system are substantially the same, more IG processors 252-256 may easily be added to the imaging system to increase the processing speed and capability.
In step 208, a preprocess or data correction process submodule 258 of the IG processor 254 processes the raw image data. For example, preprocessing may comprise normalizing the raw image data by applying corrections based on calibration data and offset data particular to the CT imaging system. Alternatively, the PET scanner system 100 may implement a data correction process in which calibration data or other knowledge of the PET scanner system 100 or patient 18 is used. For example, CT attenuation correction (CTAC) may be implemented for the PET scanner system 100, wherein CT images are used to obtain information about the patient's 18 anatomy and used to correct the image data acquired by the PET scanner system 100.
In step 210, a view-weight submodule 260 of the IG processor 254 processes the image data. An interpolation process is performed on views and rows to compensate for the scan acquisition mode (i.e. helical, cardiac for CT imaging system 10) wherein the desired weighting functions are based on the scan and reconstruction modes. The view-weight submodule 260 may generate view-weighting weights based on view-weighting parameters and applies the weights to the image data.
In step 212, the rebinning submodule 262 of the IG processor 254 performs rebinning of the image data. For example, the rebinning submodule 262 may perform an interpolation process on views and rows acquired by the CT imaging system 10 to transform image data from fan-beam data format to parallel-beam data format. Equal space rebinning may be applied by the rebinning submodule 262 to process either CT or PET image data, and Fourier rebinning may be implemented by the rebinning submodule 262 for PET image data to create a 2D data set from a 3D data set.
In step 214, a filter submodule 264 of the IG processor 254 filters the image data. The selected reconstruction filter may be based on the reconstruction mode or input by the user. The filter submodule 264 generates one or more Filters tables based on the Filter parameters and may access a preselected protocol or a reconstruction filter input by a user to generate the Filters tables. The filter submodule 264 may apply a mathematical filtration based on the Filters tables to the image data on a view basis in response to the selected reconstruction filter.
In step 216, a back projection or iterative image reconstruction submodule 266 of the IG processor 254 processes the image data. For example, in the CT imaging system 10, back projection of view data and summation into an image matrix may be performed. Alternatively, 2D or 3D iterative reconstruction methods for PET and CT may be implemented. In this case, one or more of the previous steps are optional, as filtering (step 214) is not required for iterative reconstruction techniques.
By way of example only, one or more additional submodules may be implemented if the desired reconstruction mode is the Iterative Bone Correction Algorithm. A reprojection submodule forms view data from image matrix data (step 216), then the iterative bone correction algorithm is applied.
In step 218, a post processing submodule 268 of the IG processor 254 performs additional processing on the image data, such as ring fix and image filters.
In step 220, the IG processor 254 sends the processed image data to the computer, such as computer 38 of
Therefore, it should be understood that each of the IG processors 252-256 comprise a copy of the entire image generation work flow which runs in parallel across the multiple IG processors 252-256. Performance is enhanced by reducing the image reconstruction times, and thus the overall time required for CT and PET studies. By adding additional IG processors 252-256, additional and/or more complex image reconstruction algorithms may be performed, improving image quality and/or quantity without negatively impacting the image reconstruction or exam times. Thus, the IG processors 252-256 provide improved performance, flexibility and reliability.
Alternatively, individual IG processors 252-256 may be directed, or configured via software, to support unique image generation process models, enabling simultaneous support for multiple image reconstruction modes. The reconfiguration of the IG processors 252-256 via software may be changed on the fly as the IG processors 252-256 are not dedicated to any particular image reconstruction mode.
For example, a first IG processor 252 may be directed to process raw image data from a first scan and a second IG processor 254 may be directed to process raw image data from a second scan wherein the IG processors 252 and 254 utilize different image generation process models. Alternatively, a first IG processor 252 may be directed to process raw image data using a first image reconstruction mode and a second IG processor 254 may be directed to process the same raw image data using a second image reconstruction mode. A third IG processor 256 may be directed to process raw image data in either of the first and second image reconstruction modes, operating in parallel with one of the first and second IG processors 254 and 256 on a subset of the same raw image data, directed to process raw image data according to a third image reconstruction mode, or to process raw image data from a different scan.
In addition, image reconstruction corrections specific to individual rows of data or processing multiple rows of data at one time can now be supported. Previously, processing speed may have been improved by processing a single row of data, which was then combined with other single rows into an image. By way of comparison, the IG processors 252-256 can process data from individual or multiple rows of data, providing more system flexibility and processing speed.
By utilizing the data decomposition model of
As discussed previously, the processing flow illustrated within the IG processors 252-256 in
Also, the image reconstruction processing flow may not utilize all of the submodules 258-268 described above. Some of the processing steps may be optional and thus any subset of the processing steps, in combination with additional steps if necessary, forms a valid image reconstruction process. Therefore, the image reconstruction processing flow may implement the submodules 258-268 in a different order, in combination with other submodules not specifically mentioned above, or also by excluding one or more of the submodules 258-268. For example, view-weight submodule 260 may be utilized only in specific applications for the PET image data. Volume CT may be accomplished by a first process flows of rebinning, filtering, and view-weighting; a second process flow of view-weighting, rebinning, and filtering; or a third process flow of filtering, rebinning, and view-weighting.
While the invention has been described in terms of various specific embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the claims.
The present application is a continuation application of U.S. patent application Ser. No. 10/999,347, filed Nov. 30, 2004, which is hereby incorporated by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
5970111 | Samarasekera et al. | Oct 1999 | A |
5987094 | Clarke et al. | Nov 1999 | A |
6044125 | Flohr et al. | Mar 2000 | A |
6081612 | Gutkowicz-Krusin et al. | Jun 2000 | A |
6169817 | Parker et al. | Jan 2001 | B1 |
7139416 | Vuylsteke | Nov 2006 | B2 |
20050207630 | Chan et al. | Sep 2005 | A1 |
Number | Date | Country | |
---|---|---|---|
20090169085 A1 | Jul 2009 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10999347 | Nov 2004 | US |
Child | 12367627 | US |