1. Field of the Invention
The present invention relates generally to the field of computed tomography (CT). In particular, the present invention provides a method of reconstructing a three-dimensional image of an object's volume of interest from a set of conical-beam, intensity-modulated projections of this object. Also provided are image processing devices, examination apparatus, as well as a computer-readable medium and a program element that is adapted and configured to perform aspects of the methods disclosed herein.
2. Description of Related Art
Compared to traditional slice-at-a-time tomographic instruments, the cone-beam (CB) CT offers faster scans, higher patient throughput, significant reduction in X-ray dosage, and isotropic resolution. It has a great potential to be applied to a wide range of medical and industrial applications.
Cone-beam computerized tomography (CBCT) with flat panel detectors (FPD) has become a prevalent three-dimensional (3-D) imaging system for clinical application. Its short time acquisition and reconstruction with mobile C-arm system allow the integration of 3-D imaging into the interventional procedures such as endovascular techniques. In comparison with image intensifiers, the flat panel detectors provide no geometrical distortion, a higher dynamic range, and homogeneous signal distribution1. The Feldkamp, Davis, Kress (FDK) algorithm2 and its modifications in the filtered-backprojection (FBP) family have been implemented for reconstructing approximate 3-D images from circular CB projection data. Correction algorithms such as scatter correction, beam hardening correction, truncation correction, and ring correction have been developed for removing artifacts of reconstructed image prior to and after FDK reconstruction algorithm3.
Because an X-ray tube generates X-rays in all directions, beam restrictors such as collimators or aperture systems are set up in the X-ray tube for absorbing unnecessary X-ray beams outside field of view (FOV). Satisfactory CBCT systems employing collimators to reduce radiation exposure, however, have not been implemented. As CBCT scans become more common, patients are exposed to greater amounts of X-ray radiation. This exposure carries an increased risk for developing cancer, especially for children and for individuals who require multiple scans. In most cases, the region of interest for diagnosis is smaller than the window of exposure to radiation, meaning that healthy surrounding tissue receives a full dose of radiation despite not being used in the diagnosis of the patient's condition. Current CT scanning technology, unlike certain other scanning technologies, does not allow for shaping the window of exposure because of the complexity of interpreting less than the full radiation dosage across the imaging detectors.
The present invention overcomes these and many other limitations inherent in the prior art by providing, inter alia, intensity-modulated conical-beam CT imaging methods that significantly lower the overall radiation exposure during a computed tomography (CT) scan. The disclosed methods also provide a mechanism for interpreting a collimated, intensity-modulated beam, thus enabling a medical practitioner to narrow the window of radiation exposure to encompass only the particular region of interest to be imaged.
Unlike conventional CT instrumentation that uses constant beam intensity for every projection angle to produce the final image, the present invention provides an aperture, such as a collimator, to modulate the intensity of the X-ray beam, thereby reducing the amount of radiation impinging on the object of interest. In the case of medical CT imaging, the methods result is significantly lower background radiation, and a reduced level of radiation exposure to surrounding tissues, while still maintaining a sufficient intensity for imaging the target tissue. CT instrumentation employing aperture-controlled CB energy provides overall reduction in radiation exposure, which is expected to lessen the risk of adverse effects, including carcinogenesis, in the patient undergoing CBCT imaging and examination.
In illustrative embodiments, the method is accomplished using one or more aperture-controlling and/or beam intensity-restricting or limiting devices, such as, for example, a plurality of collimators, preferably placed between the radiation source and the object, to modulate the intensity of the energy beam. Modulation of the beam intensity using the intensity modulating device is accomplished such that one or more regions of interest (ROI) (alternatively referred to as a “volume of interest” [VOI]) for which analysis and imaging is desired, receive a dose of radiation sufficient to effect the desired imaging characteristics, while the surrounding tissue outside the ROI/VOI receives much lower (i.e., attenuated) ancillary radiation. By narrowing the target “window” to encompass substantially only the ROI/VOI, the present invention significantly improves the quality of diagnostic imaging, while also reducing the risk(s) that are inherent to a patient undergoing imaging using modalities (such as CT) that employ ionizing radiation.
The purposes for applying collimation to CBCT are to reduce exposing areas of patients that need not be imaged and to moderate the scatter effect, which degrades image quality4. Collimators can be designed as variable diaphragms composed of movable pieces of metal. They include two movable pieces of metal in the longitudinal direction (perpendicular to the plane of X-ray source trajectory) and two in the transverse direction (tangential to the circle of X-ray source trajectory). The relation between scatter effect and collimation on different organs of a phantom has been previously demonstrated5. The experimental results showed that scatter-to-primary ratio decreases when FOV in the longitudinal direction is reduced. Moore et al. calculated the contrast-to-noise ratio (CNR) of a targeted area using aluminum-made filters and lead-made collimators with varying thickness and with different fan angles of FOV in the longitudinal direction6. The present invention demonstrates significant improvements in CNR can be achieved when using collimated beam intensity-modulated imaging versus traditional non-collimated CT imaging methodologies and instrumentation.
In an overall and general sense, CT systems capable of providing data for reconstructing a three-dimensional (3-D) object using one or more rotations with an intensity-modulatable beam of radiation, typically include a source of X-rays distributed along a line parallel to an axis of rotation, an aperture control positioned perpendicular to the axis of rotation and near the x-ray sources so to limit the X-rays illuminating the object to contain only X-rays that travel substantially along lines perpendicular to the axis of rotation, and a detector for detecting and measuring transmitted X-rays emitted by the source. In illustrative embodiments, the aperture for controlling the size of the CB may include one or more collimators to restrict the X-rays primarily to impinge primarily upon the ROI/VOI in the object being imaged. Exemplary collimating devices may include one or more arrays or pluralities of collimators, positioned to modulate the intensity of the radiation source. In certain embodiments, the collimators may include a plurality of flat collimator plates, blades, lamellas, gates, vanes, irises, pinholes, or windows, or any combination thereof, including without limitation, those composed of a piece of dense plastic, composite, alloy, or metal, including, without limitation, aluminum, lead, or tungsten, that effectively block at least a first portion of the beam from reaching the object.
In a first embodiment, the invention provides a method for operating computed tomographic imaging using a radiation source and a plurality of detectors to generate an image of an object. In an overall and general sense, this method involves one or more steps including a) defining desired image characteristics; b) performing calculations to determine the modulation intensity to be applied to a radiation source by at least a first collimator to generate the desired image characteristics; and c) modulating the radiation source using the at least a first aperture system or collimator device to generate a desired pattern of fluence between the beam source and the object to be imaged. Exemplary desired image characteristics include, without limitation, desired levels of contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), or a combination thereof. The method may involve providing at least one desired image quality in at least one defined region of interest, or alternatively, providing at least a desired plurality or distribution of images or image qualities. In certain applications, the step of performing the calculations will include solving an inverse problem using one or more iterative solutions. Such methods may also further optionally include a step of defining at least a first region of interest from a library of population modules, or from at least one previously acquired image of the object.
In the practice of these methods upon animal (and preferably, human) subjects, the total radiation dose to the organism will be lower (and in some cases substantially lower) than that typically required performing a similar method in the absence of intensity modulation, and/or without using at least one modulatable aperture or collimator device to modulate the intensity and fluence of the radiation beam emanating from the radiation source.
In certain embodiments, the method will include providing at least a first temporal modulation of the radiation source, at least a first spatial modulation of the radiation source, or a combination of both spatial and temporal modulation of the radiation source.
In exemplary embodiments, the collimator will include a plurality of individual elements adapted to absorb radiation, and in some instances, this plurality of individual elements is composed of lead or other such like element, compound, or composition.
In another embodiment, the invention provides an imaging system adapted and configured to perform one or more of the methods described herein. In the practice of the invention, the imaging system may include at least a first collimator that is composed of two or more individual elements, each of which is substantially impervious to the radiation, and being substantially movable between an open position and a closed position, whereby open positions of the individual elements within the collimator define an aperture that permits passage of the beam from the radiation source to impinge upon the object being imaged. Such imaging systems include, without limitation, CT systems, and preferably, CBCT systems, including, without limitation, those CBCT systems in which a computer forms at least a first image using at least a first mathematical algorithm adapted and configured for obtaining and/or analyzing an image of an object of interest.
In certain embodiments, the imaging system will further optionally include one or more post-processing modules that may be useful in the enhancement of at least a first three-dimensional model of at least a first region of interest from on, or substantially within one or more portions of the object being imaged by the system. Preferably, the plurality of X-ray projections will originate from a circular, or even helical, trajectory of one or more X-ray radiation sources around, or in substantial proximity to, the object of interest being imaged.
In another embodiment, the invention provides a record carrier on which a computer program for the generation of a three-dimensional model of at least a first region of interest of an object from a plurality of collimated X-ray projections is stored. The record carrier preferably contains a computer program that is adapted to execute at least one, and preferably, substantially all of the iterative steps of one or more of the methods disclosed herein.
The invention also provides a computer-readable medium, onto which a computer program for reconstructing one or more three-dimensional image(s) of an object's ROI/VOI from a set of collimated CB X-ray projections of the object is stored. Preferably, the computer program, when being executed by a processor, is adapted to execute at least one, and substantially all, of the steps of one or more of the image acquisition and/or image analysis methods disclosed herein.
The invention, in another embodiment, also provides a program element for reconstructing a three-dimensional image of an object's ROI or VOI from at least a first set of CB X-ray projections of the object, which, when being executed by a processor, is adapted to execute at least one step, and substantially all, of the steps of one or more methods as disclosed herein.
In another embodiment, the invention provides an examination apparatus that includes (a) an X-ray device for the generation of X-ray projections of an object's region or volume of interest from at least two different directions, wherein the projections are obtained from at least two different samplings of a collimated beam of X-rays generated from the device; and (b) at least a first imaging acquisition or analysis device as disclosed herein. In one embodiment, the examination is characterized as a CBCT system, and is preferably adapted and/or configured to provide at least a first intensity-modulated (i.e., aperture-limited or collimated) beam of X-rays for imaging at least a first region of interest of an object that is undergoing examination in the apparatus by one or more technicians, operators, medical practitioners, or such like.
In exemplary embodiments, the examination apparatus as disclosed herein may be adapted or configured as a baggage inspection apparatus, a medical diagnostic apparatus, a material testing apparatus, or a materials science analytic apparatus, or such like.
The invention also further herein provides methods for generating three-dimensional images of an object using a plurality of CB projections that are passed through the object and attenuated thereby. In an overall and general sense, such methods include at least the steps of a) positioning a source at a position on a predetermined scan path; b) passing a projection of CB X-ray radiation comprising a plurality of projection rays from the source through an object, wherein the CB projection is attenuated by at least one or more collimators adapted and configured to attenuate the projection; c) detecting the radiation intensity of the CB projection passing through the object of interest onto an area of at least a first detector; d) obtaining a two-dimensional attenuation image of the CB projection from the detected radiation intensity; e) generating an intermediate, locally reconstructed, three-dimensional image with constant values assigned along each projection ray; f) at least once, repeating steps (d)-(e); and g) summing the plurality of intermediate, locally reconstructed, three-dimensional images obtained for the plurality of CB projections to obtain an ultimate, reconstructed, three-dimensional image of the object.
The invention also further provides a CBCT apparatus that includes: (a) at least a first radiation source, (b) at least a first radiation detector, (c) at least a first support for an object of interest, or alternatively, at least a first region of interest from such an object, to be illuminated by a collimated beam of radiation projected from the radiation source, and (d) at least one computer-readable storage medium for storing computer-executable software to generate a reconstruction of the CB radiation attenuation in the illuminated object. In such applications, the software will preferably include: (a) computer-executable code for obtaining and/or generating at least a first intermediate, locally reconstructed, three-dimensional image; (b) computer-executable code for summing a plurality of intermediate, locally reconstructed, three-dimensional images obtained for the CB projections to obtain an ultimate, reconstructed, three-dimensional image of the object; and (c) computer-executable code for displacing the source and detector relative to the support in a predetermined scan path for radiation transmitted from the source, through an object positioned by the support, and to the detector.
The invention also provides methods and systems for forming an image of an object that generally involves) exposing the object to a cone beam of radiation rendered spatially coherent by its passage through at least a first collimator; (b) projecting the spatially coherent collimated conical beam of radiation onto the object and collecting the radiation which has passed through the object in at least a first detector to produce detected data; (c) deriving, from the detected data, at least a first image; and (d) repeating step (c) at least once to form a plurality of images of at least a first region of interest within the object. Exemplary systems include, without limitation, CBCT system adapted and configured for medical diagnostic imaging and/or therapy.
According to another exemplary embodiment of the present invention, a computer-readable medium may be provided, in which a computer program for reconstructing an image from a set of projections of an object of interest with an examination apparatus is stored which, when being executed by a processor, is adapted to carry out the above-mentioned method steps.
The present invention also relates to a program element of reconstructing an image from a projection data set of an object of interest, which, when being executed by a processor, is adapted to carry out the above-mentioned method steps. The program element may preferably be loaded into the working memory of a data processor. The data processor may thus be equipped to carry out exemplary embodiments of the methods of the present invention. The computer program may be written in any suitable programming language, and may be stored on a computer-readable medium, such as a computer hard drive, a CD-ROM, DVD-ROM, or the like. In addition, the computer program may be available from a network, from which it may be downloaded into image processing devices or processors, or any suitable computer(s).
For promoting an understanding of the principles of the invention, reference will now be made to the embodiments, or examples, illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended. Any alterations and further modifications in the described embodiments, and any further applications of the principles of the invention as described herein are contemplated as would normally occur to one of ordinary skill in the art to which the invention relates.
The following drawings form part of the present specification and are included to demonstrate certain aspects of the present invention. The invention may be better understood by reference to the following description taken in conjunction with the accompanying drawings, in which like reference numerals identify like elements, and in which:
Illustrative embodiments of the invention are described below. In the interest of clarity, not all features of an actual implementation are described in this specification. It will of course be appreciated that in the development of any such actual embodiment, numerous implementation-specific decisions must be made to achieve the developers' specific goals, such as compliance with system-related and business-related constraints, which will vary from one implementation to another. Moreover, it will be appreciated that such a development effort might be complex and time-consuming, but would be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure.
In CBCT, a set of X-ray cone-beam projections of an object's ROI is acquired while an X-ray source moves along some source trajectory around the object. Using a CBCT scanner, a 3-D image of the object's ROI can be reconstructed from this set of cone-beam projections. Conventional CBCT scanners are equipped with a point-like X-ray source and a large-area X-ray detector. Typically, the detector is flat and subdivided into a two-dimensional (2-D) array of small detector elements. Such a flat detector naturally defines a 2-D plane in 3-D space, i.e., the detector plane, which is rigidly attached to the detector and moves as the detector moves. The cone-beam is formed by those X-rays that emanate from the source and intercept the detector. A 3-D image is reconstructed from the set of cone-beam projections by means of an image-processing device, which is typically a computer that executes a computer program that implements a reconstruction algorithm. The reconstructed image provides a sampled estimate of the 3-D map of the X-ray attenuation coefficient within the ROI of the imaged object.
CBCT scanners may be realized in various ways, which are known to the person of ordinary skill in the medical imaging arts. U.S. Pat. Nos. 6,582,120 and 5,124,914 (each of which is specifically incorporated herein in its entirety by express reference thereto) illustrate examples of conventional arm-based and gantry-based CBCT apparatus. C-arm based scanners are normally equipped with a flat detector, while gantry-based scanners have either flat or non-flat detectors. Flat detectors have natural detector planes associated with them that are usually subdivided into a 2-D array of small detector elements. In such a case, each CB projection of the acquired set of CB projections will be sampled on an equidistant Cartesian grid in the detector plane. In non-flat detector models, a virtual detector plane can be associated with the detector in which each CB projection can be sampled in the virtual detector plane. In these detectors, the grid of sampling points is be planar, but not necessarily equidistant Cartesian.
CBCT scanners include one or more image processing devices that reconstruct the acquired images by executing one or more particular reconstruction algorithms, and typically include one or more viewing consoles for displaying reconstructed images, and is configured and adapted to allow an operator, such as a CT technician or radiologist, to operate the device, as well as obtain and/or record one or more images obtained by the scanner.
A number of CT systems, including CBCT, and methods for collimating, processing and forming an image using such systems are known in the art, and exemplified, without limitation, in U.S. Pat. Nos. 7,424,090, 7,418,082, 7,412,029, 7,397,887, 7,396,162, 7,386,088, 7,379,529, 7,375,338, 7,372,937, 7,339,174, 7,317,819, 7,310,411, 7,308,072, 7,305,063, 7,292,717, 7,292,674, 7,289,599, 7,249,886, 7,245,697, 7,233,644, 7,209,535, 7,203,272, 7,187,750, 7,176,466, 7,151,817, 7,147,372, 7,145,981, 7,142,633, 7,116,749, 7,113,570, 7,103,137, 7,101,078, 7,099,428, 7,085,345, 7,076,023, 7,072,444, 7,072,436, 7,058,159, 7,050,534, 7,046,757, 7,042,975, 7,027,558, 7,023,956, 7,020,232, 7,015,477, 7,015,460, 6,865,246, 6,842,502, and 6,735,277; as well as U.S. Patent Appl. Publ. Nos. 20080253516, 20080226018, 20080212860, 20080212735, 20080205727, 20080205587, 20080205585, 20080192886, 20080187195, 20080165923, 20080165918, 20080144904, 20080137802, 20080118032, 20080116386, 20080112531, 20080104530, 20080095305, 20080095301, 20080089468, 20080061241, 20080031507, 20070269001, 20070172026, 20070140410, 20070104310, 20070003010, 20060269049, 20060262898, 20090292898, 20060067464, 20060002506, 20030235265, 20030202637, and 20030043957, the contents of each of which is specifically incorporated herein in its entirety by express reference thereto.
Collimators in CBCT systems are designed to acquire customized field of view (FOV) for a smaller object and to reduce radiation dose and scatter effect. In most applications, only a region of interest (ROI) within a large body needs to be imaged and acquired for reducing radiation exposure to patients. However, the reconstructed ROI image from incomplete projections may contain severe artifacts.
In the example that follows, an imaging method is described, which customizes a small ROI within an object by applying adjustable collimators and a simple and practical reconstruction by truncation correction in order to reduce a large amount of unnecessary radiation dose and to obtain tolerable image quality of the reconstructed image compared with that without ROI imaging.
In the studies described herein, ROI imaging was performed by setting transversely and longitudinally movable lead collimators between X-ray source and the object. The ROI is constrained to a small area within an object. Entrance dose has been utilized for estimating the radiation dose absorbed in the body and conduct experiments for measuring radiation dose with ROI imaging. Two reconstruction techniques for possibility of overcoming truncation problems were investigated: one was the truncation correction technique with the Feldkamp, Davis, Kress (FDK) algorithm and the other was π-line reconstruction algorithm.
The truncation correction technique was applied with FDK to reconstruct the image of cylindrical polystyrene with water as background. The imaged regions are small circles with diameters of 70- and 94-mm. From the radiation dose measurement, it was found that at least 60% and 70% of radiation dose was reduced in each of ROI imaging. Furthermore, the image quality was still acceptable with little variation of image gray values; however, it was not possible using the analytical reconstruction algorithms employed to fully reconstruct a ROI that was totally encompassed by the object.
A truncation correction technique was applied to ROI imaging with ROI within the object, and a surprising discovery was made: at least half of radiation dose could be reduced without sacrificing significant image quality. These results demonstrated that the use of collimators to reduce radiation dose during CBCT imaging represents a fundamental change in the way medical CT imaging may be performed in the future.
The detection, acquisition, analysis, and display of CT images collected using the disclosed methods and apparatus may conveniently be implemented using a conventional general purpose computer or micro-processor programmed according to the teachings of the present invention, as will be apparent to those skilled in the computer arts. Appropriate software may readily be prepared by programmers of ordinary skill based on the teachings of the present disclosure, as will be apparent to those skilled in the software arts. Exemplary computers may include a motherboard that contains a CPU, memory (e.g., DRAM, ROM, EPROM, EEPROM, SRAM, SDRAM, Flash-RAM, or the like), and other optional special purpose logic devices (e.g., ASICS, and the like) or configurable logic devices (e.g., GAL, reprogrammable FPGA, and the like). The computer also includes plural input devices, including, without limitation, a keyboard and a mouse, one or more display/graphics card or adaptors for controlling the output to one or more monitor or display devices. Additionally, the computer may include removable media devices, such as, without limitation, a compact disc, magnetic tape, removable magneto-optical media, and the like); and one or more fixed, or high-density media drives including, without limitation, a “hard disk,” connected using an appropriate device bus (e.g., a SCSI bus, an Enhanced-IDE bus, an Ultra-DMA bus, and the like). The computer may also include a compact disc reader, a compact disc reader/writer unit, which may be connected to the same device bus or to another device bus.
Examples of computer readable media associated with the present invention include, without limitation, compact discs, hard disks, magnetic tapes, magneto-optical disks, programmed read-only memory (e.g., without limitation, EPROM, EEPROM, Flash-EPROM, and the like), DRAM, SRAM, SDRAM, Rambus-DRAM, etc. and one or more combinations thereof. Stored on any one or on a combination of these computer readable media, the present invention includes software for controlling both the hardware of the computer and for enabling the computer to interact with a human user. Such software may include, without limitation, one or more device drivers, one or more operating systems and one or more user applications, including development tools and the like. Computer program products of the present invention include any computer readable medium which stores computer program instructions (e.g., computer code devices) which when executed by a computer causes the computer to perform the method of the present invention. The acquired data may be digitized, if not already in digital form, and alternatively, the source of image data being obtained and processed may be a memory storing data produced by an image acquisition device, and the memory may be local or remote, in which case a data communication network may be used to access the image data for processing according to the present invention.
The following example is included to demonstrate illustrative embodiments of the invention. It should be appreciated by those of skill in the art that the techniques disclosed in the example that follows, represent techniques discovered by the inventors to function well in the practice of the invention, and thus can be considered to constitute preferred modes for its practice. However, those of ordinary skill in the art should, in light of the present disclosure, appreciate that many changes can be made in the specific embodiments which are disclosed and still obtain a like or similar result without departing from the spirit and scope of the invention.
When objects or bodies extend outside scan FOV, there is abrupt discontinuity at the edges of detector, and projection data is truncated since missing data outside the detector is theoretically assumed to be zero during CBCT reconstruction. In practice, nearly all of the reconstruction algorithms in CBCT scanners are based on the FDK algorithm. Performing one-dimensional (1-D) filtering on truncated projection data leads to bright band artifacts extending inside scan FOV and incorrect reconstruction around the edges of scan FOV7.
Wf=p, (1)
where f represents the image vector, p the projection vector, and W the weighting matrix. Equation (1) is expressed as a linear equation from line integral (Radon transform). The image vector f can be obtained by solving the inverse problem in least square sense. As described previously, biased but informative solutions can be iteratively computed10. However, it lacks computation speed, and due to severe truncation, little information in p is available to increase accuracy of the reconstructed image, f.
A number of truncation correction techniques commonly implemented on mildly truncated projections as the case in
When changes of 3-D images are confined to a small region during interventions, ROI imaging is applied during the second scan in order to save radiation dose. Projections are generally truncated at all views as in
In this example, ROI imaging technique uses movable pieces of lead as collimators. The collimators make the scan FOV reduced such that the imaged region is as small as a circle having a diameter of approximately 70-mm, and all projections were under severe truncation at all views as shown in
ROI Imaging with Collimators
ROI Imaging Truncation Correction
First, the truncation correction technique is as designed as follows: letting Pk(u,v) be the k-th projection image with u-v coordinate shown in
where ur and ul are determined by the right and left edges of the truncated projection PFOV(u,v). Assuming the lengths for extrapolations at the both sides are Nr and Nl with 0<Nl,Nr<ur−ul, based on the symmetric mirroring extrapolation technique11, the extrapolated projection (
By assuming that the image variation in ROI is small compared with the magnitudes PFOV(ur,v) and PFOV(ul,v), the assumption, in general, is held for small ROI imaging of a body or head. Therefore, the extrapolated projection Pex(u,v) is positive under this assumption:
P
FOV(ur,v)>|PFOV(ur,v)−PFOV(−u+2ur,v)|
P
FOV(ul,v)>|PFOV(ul,v)−PFOV(−u+2ul,v)|
To obtain smooth transition of the extrapolated data to zero, Pex(u,v) is weighted by Equation (4):
The weighted extrapolation is illustrated in
π-Line Reconstruction Algorithm
The π-line reconstruction algorithm is described below:
Assuming that the π-line segment
r=r(θA)+zruz+λuπ (6)
where λε[0,1], and uz=(0,0,1). The reconstruction has two steps: first, back-project the projections into the π-line segments, and secondly, perform Hilbert transformation along each of π-line segments. The 3-D volume is union of π-line segments and, for simplicity, we describe the practical reconstruction of 3-D images only on a π-line segment based on theoretical results previously described21.
Considering the π-line segment,
where d(r,θ) was defined as the distance between r and r(θ). Pθ(ur,vr) is the line-integral which starts from the source at the rotational angle θ, passes through the point r, and ends at a point (ur,vr) on the detector.
Finally, the reconstructed image along the π-line segment
where c is the midpoint of
Calculation of Radiation Dose Reduction
An “entrance dose” of X-rays is the dose absorbed at surface of the skin where X-ray beams enter. An X-ray beam enters the body from the direction of the X-ray tube. A small share of the beam exits from the body on the opposite side, where it exposes the detector. The share of the beam that never exits from the body is absorbed as extra energy by the body's internal organs and bones. Between the collimators and the patient, imagine that there is a virtual plane parallel to the detector and close to the top of the patient illustrated in
Letting X-ray intensity be a unit of energy per unit area, and assuming that the energy which is obtained by integrating the X-ray beam intensity with a small sphere surrounding the X-ray source is Es in
where 4πd2 is the surface area of a sphere with radius d. Assuming that the coordinate of the point Q is (x,y), then the intensity IQ at the point Q is
From Equation (9) and Equation (10), IQ becomes
The obliquity effect occurs since surface of the detector is not perpendicular to the direction of X-ray beams propagation. In
IQ′=IQ cos θ (12)
where cos
The combination of inverse square law and obliquity effect is multiplicative. From Equation (11) and Equation (12), the overall intensity at the point Q is given by:
Letting the area projected on the virtual plane within FOV be AFOV shown in
Letting AFOV,m be the maximal area on the virtual plane when full-field imaging is applied, and denoting the entrance dose as DE,m, entrance dose was defined relative to that in the full field imaging in Equation (15) as:
i.e., D% is the entrance dose relative to the maximal entrance dose DE,m.
Measurement of Radiation Dose
Dose measurement in the study was made using a cylindrical CT dose-index (CTDI) phantom with a diameter of 160 mm. This is the standard Perspex® CTDI head phantom with a central hole and four peripheral holes 10-mm below the surface.
CTDI100=0.876×2×R (16)
where 2 is a modified factor for the specific CT chamber in the measurement, and 0.876 is a factor converting Roentgen units into rads in air. The weighted dose CTDIw was also calculated from Equation (17) as:
CTDIw=⅓CTDI100,c+⅔CTDI100,p, (17)
where CTDI100,c and CTDI100,p are measured doses at the center and at the peripheral, respectively.
The CBCT imaging of the described studies were based on the C-arm flat panel (FP) CT system (Siemens Medical Solutions, Forchheim, Germany), and the parameters of CBCT system are given in Table 1. The collimators were made of lead with a thickness of 5 mm, which was enough to attenuate the intensity of X-ray beams close to zero. Two collimator pieces are movable in the transverse direction, while the two remaining collimator pieces are movable in the longitudinal direction as illustrated in
FDK & π-Line Reconstruction Under Severe Truncation Correction
Using a head phantom with tissue- and bone-like materials as an object, two studies were conducted: the first was “without collimation” imaging, and the other was “with collimation” imaging with Ct=130 mm and Cl=130 mm. The reconstructions were the FDK algorithm and π-line BPF algorithm, each of which is applied in both studies. The projection images were stored as integers within the range [32,4096], and were scaled to [0,1].
ROI Reconstruction
Cylindrical polystyrene (about −30 HU) was used with a diameter of 30 mm and height of 25 mm as the object in the ROI and water as the background. The cylindrical polystyrene was stacked above a cylindrical plastic water phantom (about 30 HU) in a tank of water. These materials and background provide constant attenuation coefficient that clearly demonstrated the variation of ROI images in the “with collimation” imaging. The projection of the cylindrical polystyrene is approximately in the isocenter of the detector at all views.
Three studies were performed: one “without collimation” imaging, another “with collimation” imaging with Ct=150 mm and Cl=83 mm and the third “with collimation” imaging where Ct=110 mm and Cl=77 mm. The later study “with collimation” imaging was conducted with the smallest FOV, which enabled the projection of the cylindrical polystyrene to exactly fit within FOV at all views. The CBCT 3-D reconstruction used the FDK algorithm with ring correction, scatter correction, and overexposure correction in the CBCT system. In the “with collimation” imaging, the truncation correction described above was applied on projection images before FDK reconstruction. The extrapolation lengths Nr and Nl were assumed to be one-half Ct.
Simulation for Entrance Dose
In these simulations, the calculation of entrance dose was based on parameters of the C-arm flat panel (FP) CT system shown in Table 1. The estimated source-to-virtual-plane distance was 650 mm, and the estimated AFOV,m 35213.75 mm2 (162.5×216.7 mm2), which was AFOV when the “without collimation” imaging is applied.
Radiation Dose Measurement with Collimation
The radiation doses were measured with varying open areas of transverse and longitudinal collimators. The size Ct×Cl of full field (“without collimation”) was 385×295 mm2, and its measured doses, tube voltages, and tube current at the center and at the peripheral were 2.35 Roentgen, 79 kV, and 251 mA; and 2.03 Roentgen, 79 kV, and 251 mA, respectively. Two studies were conducted: one only varied the length Ct with the fixed Cl=295 mm, while the other varied both length Ct and width Cl. The doses, tube voltages, and tube currents at the center and peripheral were measured in these two studies as shown in Table 2 and Table 3.
It is worth noticing that the doses were not lower at the smallest open areas (Cl=25 mm) of collimators in Table 2 and Table 3. The tube voltages and currents for Cl=25 mm in Table 2 and Table 3 were above 90 kV and below 220 mA compared with the others at about 80 kV and 250 mA. This can be explained by the automatic exposure control in the CBCT system. The tube voltage of X-ray beams would be automatically increased for fewer photons detected on the detector. Due to the exposure control mechanism, it was not possible to use the same tube voltage and current for all the measurements. It was also not possible to scale the doses to those with fixed tube voltage and current due to lack of relationship of dose and tube parameters under collimated imaging. For the other measurements, the variation of tube voltages and currents was within 5% relative to 80 kV and 250 mA, and the measured doses were still reliable in this study.
The dose for full-field (i.e., “without collimation”) was denoted as 100%.
The main purpose of simulating the two analytical reconstruction algorithms (FDK and π-line) in
In the studies involving ROI reconstruction, the object of interest is the cylindrical polystyrene with water as the background. Two studies were performed using moderate (
It was observed that the extrapolation technique in truncation correction with the FDK reconstruction could largely recover ROI image with an extremely small ROI. The various filter kernels in the FDK algorithms were derived from the ramp filter, which is formed by Sinc function. The feature of the filters is that their impulse responses decay fast toward both sides. Since the FDK reconstruction is involved with convolution of the filter and projections, the reconstructed images toward isocenter do not undergo significant distortion. If one properly extends the truncated projections by extrapolation, it mitigates the variation of reconstructed images caused by filtering truncated projections. The distortion in gray values becomes more obvious for reconstructed images far away from isocenter. Therefore, if a larger angle of FOV is applied, the distortion of the object in ROI will be lessened. This effect was clearly demonstrated in
The radiation dose calculation in computer simulations can be performed by Monte Carlo calculations in which photoelectric effect and scattering of photons are randomly chosen according to the probability distribution of the mass attenuation coefficient of a voxel in a simulated object. Validation of the Monte Carlo simulation has been described in detail elsewhere23-25. However, the dose calculation is still impractical, since the calculated dose distribution varies from object to object using the same CT system. The calculation of entrance dose in
The following references, to the extent that they provide exemplary procedural or other details supplementary to those set forth herein, are specifically incorporated herein in their entirety by explicit reference thereto:
All of the methods disclosed and claimed herein can be executed without undue experimentation in light of the present disclosure. While the methods of this invention and devices, imaging systems, and computer program products employing them have been described in terms of exemplary embodiments, it will be apparent to those of ordinary skill in the art that variations may be applied to the composition, methods and in the steps or in the sequence of steps of the method described herein without departing from the concept, spirit and scope of the invention as defined by the appended claims. Accordingly, the exclusive rights sought to be patented are as set forth in the following claims:
The present application claims priority from U.S. Provisional Application No. 61/113,655, filed Nov. 12, 2008; the entire contents of which is specifically incorporated herein by express reference thereto. Not Applicable.
Number | Date | Country | |
---|---|---|---|
61113655 | Nov 2008 | US |