An embodiment relates in general to imaging by computed tomography, and more specifically to acquisition of projection data for generating an image without motion artifacts for a moving imaging subject.
Computed Tomography (CT) images may be created by combining X-ray images taken while rotating a gantry holding an x-ray source and detector around an imaging subject. For some CT imaging systems, the best CT images may be produced when all parts of the imaging subject are stationary with respect to the imaging system. Motions by the imaging subject may cause degradation of a CT image formed from the individual X-ray images. CT image degradations from subject motion may be referred to as motion artifacts. Examples of motion artifacts from subject motion include loss of image contrast, loss of image resolution, blurring of some or all of an image, and other undesirable effects.
CT imaging has many important medical applications but has proven challenging when imaging parts of a body that undergo movement. A motion artifact in a CT image may obscure a clinically important aspect of the image or may be mistaken for an indication of the condition of the imaging subject. Diagnostic accuracy may suffer as a result of the presence of motion artifacts in CT images. For medical imaging, examples of image subject motions which may cause motion artifacts include motions of heart muscles, coronary vessels, the lungs and other parts of the chest cavity moving in response to breathing, opening and closing of valves and orifices in the circulatory and digestive systems, peristaltic waves, and other voluntary and involuntary motions of the imaging subject. For industrial imaging, examples of image subject motions which may cause motion artifacts include rotating and reciprocating parts of pumps, compressors, generators, motors, engines, linear and rotary actuators, solenoids, switches, relays, and other devices having moving parts.
Techniques and instruments have been developed to improve CT imaging of subjects having moving parts. For example, cardiac CT imaging may use an electrocardiogram (ECG or EKG) signal to determine the time interval during a heartbeat cycle when view data for a CT image will be acquired. The EKG signal may be used to predict the portion of a heartbeat when heart motions are at or near a minimum, as may occur for example during diastole for low heart rates and systole for high heart rates, where a heart rate is an expression of the number of complete heart cycles per unit time. Techniques have been developed to acquire view data for forming a CT image from a single heart cycle and from multiple heart cycles. For example, view data may be collected throughout the duration of a complete heart cycle and then CT images may be made retrospectively at selected times during the heart cycle to find the time relative to the start of a cycle, also referred to as a phase of a heart cycle, at which the best image for diagnostic purposes can be made.
Several problems are associated with using an EKG signal to control acquisition and selection of CT images. Placement of sensors and connecting wires to detect EKG signals that can be compared to diagnostic references requires skilled, experience staff who may not be available at the time a CT image is to be made. Electrical noise picked up by the EKG sensors or wires may interfere with interpretation of EKG signals. The imaging subject's heart rate may be so high that a sufficiently stable period for CT imaging cannot be located, possibly requiring the administration of drugs to reduce heart rate. Using an EKG to control CT imaging may therefore interfere with workflow in a CT imaging facility or increase the cost of CT imaging. The person responsible for selecting and interpreting CT images must know which EKG sensor, and which portion of an EKG signal from a selected sensor, may produce a CT image of the selected phase that has the fewest number of motion artifacts, and must be able to distinguish between a motion artifact and an actual medical condition of the imaging subject.
Improper visualization from motion artifacts may lead to errors in identification of pathology and possibly to misdiagnosis. For example, a physician may need to closely and carefully examine many CT images at different phases of a heart cycle to find an image that is sufficiently free of motion artifacts for accurate diagnosis. The skill of the operator of the CT instrument may affect the quality of CT images and the usefulness of the images for diagnostic purposes. The difficulty is increased when the clinically relevant portions of an imaging subject are small in size, for example when looking for the presence of plaque or blockage in the circulatory system. The time and experience needed required of clinical staff to create and sort through CT images has impeded the use of CT imaging and increased the cost of CT imaging.
Algorithms have been developed to perform pre-processing of view data in the projection or image domain to find a phase of least motion for capturing images of the heart and vessels. However, these algorithms are unable to guarantee that the resulting CT images will be free from motion artifacts, for reasons including the variability of heart rate in some patients, premature heart beats, and errors in feature detection algorithms caused by noise and low image contrast.
Post-acquisition image processing methods have been introduced to reduce motion artifacts from images of coronary vessels. These methods may depend on accurate identification of a coronary vessel tree from image data to detect and compensate for motions of the vessels, but may fail to eliminate the effects of large imaging subject motions and high noise levels from CT images.
Previously known methods have been unable to produce high-resolution CT images of rapidly moving imaging subjects without using complex gating signals such as an EKG. Selection of the wrong EKG lead, or selecting the wrong part of an EKG signal on a correctly chosen lead, may result in a CT image with motion artifacts that obscure or misrepresent diagnostically important aspects of the imaging subject. Producing high-contrast, high-resolution CT images that are free from motion artifacts, without depending on the skill and experience of an imaging system operator to interpret the physiological and imaging implications of a gating signal, are problems that have remained unsolved by previously known methods.
An example of a method embodiment includes determining a time duration tR of a data record to reconstruct exactly one motion-corrected computed tomography (CT) image of an imaging subject having a moving part; rotating a gantry holding a detector and a source through a full revolution; and while rotating the gantry, capturing a sequence of views of the imaging subject for the time duration tR. The example of a method embodiment further includes saving the sequence of views in the data record; forming a motion-corrected image of the imaging subject by partial angle reconstruction of the sequence of views from the data record; and displaying the motion-corrected CT image.
An embodiment may optionally include acquiring view data asynchronously to a heart cycle. Alternatively, an embodiment may optionally include acquiring view data synchronously with a pulse signal representative of a heart rate.
Accompanying drawings show one or more embodiments; however, the accompanying drawings should not be taken to limit the invention to only the embodiments shown. Various aspects and advantages will become apparent upon review of the following detailed description and upon reference to the drawings in which:
An exemplary embodiment improves acquisition of view data for computed tomography (CT) images by collecting a data record, and in some embodiments exactly one data record, having sufficient information to produce a motion-corrected image. A data record in accord with an embodiment may be collected asynchronously with respect to a cyclical motion of an imaging subject. A data record collected asynchronously with respect to the motions of the imaging subject may be used to generate motion-corrected CT images without using or referencing any signal representing motions of the imaging subject, such as an EKG signal or a timing pulse measured from the imaging subject. However, when a signal having components which are synchronous with motions of the imaging subject is available, embodiments may use the synchronization signal to reduce an amount of data presented to an image reconstruction system for formation of a motion-corrected CT image, and may optionally use the synchronization signal to target image acquisition at a selected time during a cyclical motion of the imaging subject.
Embodiments are advantageous for collecting projection data, also referred to as view data, to be combined into a data record submitted to an image reconstruction system for forming motion-compensated CT images. In contrast to previously known methods for CT imaging of moving imaging subjects, embodiments eliminate the screening of multiple CT images or the associated data records, either by a person or by a screening algorithm, to identify an image with motion artifacts that are too small to interfere with diagnostic accuracy.
The time duration of the one data record depends on whether the acquisition is performed synchronously or asynchronously relative to cyclical motions of the imaging subject, as will be discussed in more detail below. Synchronous acquisition refers to data records having a start time deliberately aligned in time with a specific event in a cyclical motion of the imaging subject, for example a selected heart phase or the time of a synchronization pulse. Asynchronous acquisition refers to a data record that is not deliberately time-aligned with an event in the cyclical motion of the imaging subject, and in which events occurring during the cyclical motion may appear at any time during the acquisition of the data record.
For discussion purposes herein, a cyclical motion is a motion which repeats at approximately regular intervals. A period of the cyclical motion is a time duration from an event in one cycle to the corresponding event in the next cycle. For the embodiments disclosed herein, it is not necessary for the cyclical motion to repeat with perfect uniformity. An image free from motion artifacts may still be reconstructed even when one period of cyclical motion does not exactly match another period of cyclical motion. This contrasts with some previously known methods which may collect data for reconstruction of one image by combining image data from more than one period of cyclical motion of the imaging subject, where any variation in subject motion from one period to another may lead to image degradation. Embodiments preferably collect a minimum amount of data needed to reconstruct exactly one motion-corrected CT image and preferably avoid collecting more than the minimum amount of data needed for image reconstruction. A motion-corrected CT image in accord with an embodiment is a CT image formed by partial angle reconstruction in which any motion artifacts that may be present are too small to affect the accuracy of a diagnosis of the imaging subject based on the reconstructed CT image.
CT images in accord with an embodiment are preferably formed by partial angle reconstruction. Partial angle reconstruction is capable of forming high-resolution CT images in which the subject motion is averaged over a much shorter time period than other methods of CT imaging. For example, at a heart rate of 75 beats per minute, the duration of one heart cycle is 800 milliseconds (ms). For a gantry rotation rate of 240 revolutions per minute, partial angle reconstruction is capable of producing a CT image with subject motion averaged over an interval of about 175 microseconds, a duration over which heart motion is insignificant for most diagnostic imaging purposes, except at the most rapid portion of a heart cycle represented on an EKG by an interval including the peak R in a QRS complex. Some embodiments operate to provide sufficient data for reconstruction of a motion-corrected image even when the QRS complex occurs during collection of view data. Other embodiments collect view data away from the relatively short portion of a heart cycle associated with maximum heart motion. In contrast to partial angle reconstruction techniques used by an embodiment, other previously known CT imaging methods may average subject motion over an interval of about 60 ms, long enough for many of the heart motions that occur during a heart cycle to cause motion artifacts that are large enough to interfere with diagnostic interpretation of CT images.
To perform partial angle reconstruction, a CT imaging system acquires a set of views corresponding to the fan angle of an X-ray beam and another set of views with the gantry rotated 180° from the first set of views. The two sets of views represent conjugate views of the imaging subject. The conjugate views are compared to one another to determine which parts of the imaging subject may have moved from one view to another, and a motion correction estimate is formed based on the analysis of the sequence of views. The motion correction estimate is applied to the views to compensate for subject motion, and the CT image is reconstructed from motion-compensated data to form a motion-corrected image. At present, partial angle reconstruction is capable of compensating for heart motions everywhere in a heart cycle but very close to the peak R of a QRS complex. Heart motions during other parts of a heart cycle are sufficiently long compared to the temporal resolution of the partial angle reconstruction technique that motions away from the peak of the QRS complex for a normal sinus rhythm do not cause motion artifacts that could affect diagnostic accuracy. Motions of other organs, for example respiratory motions, are substantially slower than the fastest heart motions than can be imaged and may therefore be imaged accurately without motion artifacts by partial angle reconstruction.
The image reconstruction system 140 may activate output from the source 102 through a source control and interface 134. Projection data 138 from the detector 104 may be received by the CPU 142 from a detector control and interface 136. In some embodiments, a pulse signal 150 synchronized with a cyclical motion of the moving object in the imaging subject may be detected by an optional pulse sensor 148 and communicated to the CPU 142. Examples of a pulse sensor 148 for cardiac imaging include, but are not limited to, a pulse Doppler probe, an acoustic sensor, a pressure sensor, and a pulse oximetry probe, any of which may output a pulse signal useful for synchronous acquisition of view data.
To produce a conjugate image pair for partial angle reconstruction, the gantry may be swept through a gantry rotation angle θG determined by equation (1):
θG=180°+θF+θ (1)
where θF is the fan angle of the beam, θ is an optional incremental angle used by some embodiments, and 0°<θG<360°. The time required for the gantry to traverse θG at a specified gantry rotation rate RG determines the minimum time duration tN 154 of a data record 164 having sufficient data to form a motion-corrected CT image by partial angle reconstruction.
As noted earlier, imaging by an embodiment preferably avoids the period of maximum heart motion near the peak R 160 of the QRS complex 162 in the EKG signal example 172. When a data record 164 is acquired asynchronously, at some times during the cycle the data record may by chance be collected concurrently with the maximum heart motion and at other times may by chance be collected when the heart motion is within the motion correction capabilities of partial angle reconstruction. For the examples of data records 164 in
t
RA=2×tN+tM (2)
where tM 194 is the duration of a time interval during which subject motion may not be well-corrected by partial angle reconstruction. An embodiment performs asynchronous acquisition by rotating the gantry at a specified rate of rotation RG and determining the time duration tN for the gantry to rotate through the gantry rotation angle θG at the specified rate of rotation RG.
Acquiring view data for the time duration tRA may be performed without concurrent reference to an EKG signal or any other timing signal related to the cyclical motion of the imaging subject, for example a signal representative of a heartbeat. For cardiac imaging, the interval tM 194 may correspond to the duration of the QRS complex 162, tQRS 166. A data record 170 of duration tRA 168 will always contain projection data 138 in a continuous block with time duration tN 154 to guarantee that the image reconstruction system can form a motion-corrected CT image, no matter when the interval of maximum subject motion tM 194 occurs relative to the acquisition of projection data. Only the one data record 170 of duration tRA need be processed by the image reconstruction system to form a motion-corrected image 158, without screening or sorting of multiple CT images or data records.
Heart motions sufficiently large to cause motion artifacts in previously known CT imaging systems may occur at phases of the heart cycle other than the QRS complex. For example, heart motion may be sufficient to cause undesirable motion artifacts in a CT image formed by previously known CT systems from projection data 164C that includes data from the interval tY 193 in
Some embodiments may optionally coordinate acquisition of view data with a synchronization signal to reduce an amount of view data acquired for CT image reconstruction. In the example of
When the detected time tsync 198 of the synchronization pulse 176 on the synchronization signal 180 occurs too close to the QRS complex 162, or when it is desirable to target a specific heart phase for imaging, an embodiment may perform synchronous data acquisition in accord with the example of
t
start
=|t
C
−t
D|+(tP−(tN/2) (3)
and where |tC−tD| represents the absolute value of the difference of tC and tD. Collection of view data continues for a time duration of tN 154 and the resulting data record 164 is submitted to the image reconstruction system for formation of one motion-corrected image centered at the selected heart phase 184.
Although embodiments do not require an EKG signal to form motion-corrected CT images, if an EKG signal is available it may be used as the synchronization signal 180. For example, the R peak of a QRS complex may be used to image at a specified target phase, according the operations performed with respect to
Some of the operations described herein may be performed in a different order than implied by the sequence of blocks in
Unless expressly stated otherwise herein, ordinary terms have their corresponding ordinary meanings within the respective contexts of their presentations, and ordinary terms of art have their corresponding regular meanings.
This application claims priority to U.S. Provisional Patent Application No. 62/167,126 filed May 27, 2015, titled “System and Method for Motion Free Computed Tomography Images”, incorporated herein by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
62167126 | May 2015 | US |