Classifications

• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T12/00—Tomographic reconstruction from projections
  • G06T12/20—Inverse problem, i.e. transformations from projection space into object space
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2211/00—Image generation
  • G06T2211/40—Computed tomography
  • G06T2211/416—Exact reconstruction
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2211/00—Image generation
  • G06T2211/40—Computed tomography
  • G06T2211/421—Filtered back projection [FBP]

Definitions

• the present disclosure relates to computed tomography (CT) scanning, and more particularly to a CT scanner and image reconstruction technique that employs a process that is capable of (1) processing CT data derived from any arbitrary beamline geometry, (2) reducing artifacts if the scanner has a large field of view, and (3) improving image quality while simultaneously enabling an efficient accelerated implementation on hardware.
• CT computed tomography
• Image reconstruction processes for helical CT scanning are classified into 2D and 3D processes.
• 2D processes [see Defrise, M., Noo, F., and Kudo, H., "Rebinning- based algorithms for helical cone-beam CT", Phys. Med. Biol. 46, 2001] estimate parallel tomographic datasets for tilted slices, perform 2D backprojection in each of the tilted slices and reconstruct the volume slice by slice.
• geometrical (cone- beam) artifacts increase with the cone-angle (angular extent of the detector in the direction orthogonal to the plane of the gantry rotation).
• 3D processes reconstruct voxels in the volume directly from projection data without conversion of the backprojection task to intermediate 2D problems. For a given cone angle, images reconstructed with 3D processes generate reduced cone-beam artifacts when compared with images reconstructed with 2D processes.
• filtering data along lines aligned with the source trajectory significantly reduces cone-beam artifacts.
• the tangential filtering is derived as an approximation to exact reconstruction processes (such as proposed in Katsevich, A., "An improved exact filtered backprojection algorithm for spiral computed tomography", Advances in Applied Mathematics, 32, 2004) and modifies the FDK process to be "quasi-exact”.
• the FDK processes in the literature are designed for flat or cylindrical detector arrays, which may not be applicable to modern CT scanning systems with unconventional beamline geometry (non-flat, non-cylindrical detector arrays of various geometric (sometimes arbitrary) shapes and with non- uniformly spaced detectors).
• the quasi-exact FDK process as described in [Kudo, H., Noo, F., and Defrise, M., "Cone-beam filtered-backprojection for truncated helical data", Phys. Med. Biol. 43, 1998] is only applicable to a flat detector geometry. Therefore, the use of non-conventional detector geometry requires special treatment.
• the FDK processes in the literature generate aliasing artifacts in scanning systems with large fields of view (FOV) due to the increase in the sampling interval of voxel locations as projected on the detector array with increasing distance from the center of the FOV.
• FOV fields of view
• FIG. 1 is a perspective view of a baggage scanning system including the X-ray source and detector array mounted on the rotating gantry, which can be adapted to incorporate the system and perform method described herein;
• FIG. 2 is a cross-sectional end view of the system of FIG. 1;
• FIG. 3 is a cross-sectional radial view of the system of FIG. 1;
• FIG. 4 is a simplified perspective view of helical scan geometry when a flat detector array is employed
• FIGs. 5(A) and 5(B) are simplified perspective views illustrating principles of 2D and 3D reconstruction techniques
• FIG. 6 is a simplified perspective view illustrating the application of the FDK process of creating image data from full rotation data
• FIG. 7 is a flow chart illustrating the basic steps of the FDK process
• FIGs. 8(A) and 8(B) are simplified illustrations showing the path of X-ray photons from the source to the detector array, and an example of a sampling region for a projection angle segment of 180° + 2 ⁇ m ;
• FIG. 9 is a flow chart illustrating the steps of a basic helical FDK process
• FIG. 10 is a simplified perspective drawing illustrating the use of tangential filtering to reduce cone beam artifacts in data acquired from a flat detector array;
• FIG. 11 is a flow chart showing the additional step utilized to provide tangential filtering as a part of a helical FDK process;
• FIG. 12 is a simplified view showing the relative differences in the number of samples required to resample the data using conventional flat and cylindrical detector arrays from the data acquired by an array having an arbitrary beamline geometry;
• FIG. 13 is a simplified perspective view showing how a scanner employing a non-flat detector array of an arbitrary beamline geometry can employ the helical FDK process by resampling data onto a virtual flat detector, according to the improvements described herein;
• FIG. 14 is a flow chart illustrating the steps of the FDK process using additional improvements in a scanner of the type employing a non-flat detector array of an arbitrary beamline geometry;
• FIG. 15 (A) shows a reconstruction coordinate system (CS) including an axial gantry CS and a detector array in a cylindrical CS, while FIG. 15 (B) shows detectors in a uw-plane of the detector cylindrical CS;
• CS reconstruction coordinate system
• FIG. 16 is a schematic illustrating resampling data from a detector array of an arbitrary beamline geometry onto curved lines of a virtual cylindrical detector array corresponding to titled lines of a virtual flat detector array;
• FIG. 17 is a graphical illustration of the filtering lines for an equiangular detector array visualized in cylindrical detector coordinates
• FIG. 18 is an illustration of an example of a reconstruction plane R, and helical source trajectory
• FIG. 19 is detector coordinates representation of a trajectory traced by a reconstruction point, illustrating the results of upsampling and sparse projection computation
• FIGs. 20(A) and 20(B) illustrate resampling approach on detector arrays respectively without and with row dithering
• FIG. 21 is a schematic view illustrating the problem of the increase in sampling interval of projected voxel locations with increasing distance from the center of the FOV;
• FIG. 22 is a schematic view illustrating the technique of upsampling of the projection trajectory;
• FIG. 23 is a schematic view illustrating an adaptive upsampling technique of reducing computation;
• FIGs. 24(A) and 24(B) are schematic views illustrating a backprojection technique using sparse projections.
• FDK processes in the literature are designed for flat or cylindrical detector arrays.
• the enhanced FDK process including tangential filtering is only designed for a flat detector array. If directly applied, such processes may not be applicable to modern CT scanning systems with unconventional beamline geometry.
• the terms "unconventional beamline geometry” and "arbitrary beamline geometry” each mean a beamline geometry in which the detector array represents a non- flat, non-cylindrical surface with possibly non-uniformly spaced detectors. Therefore, as described hereinafter, the data can be resampled into a flat or cylindrical detector geometry before reconstruction.
• FDK processes in the literature generate aliasing artifacts in scanners with large fields of view (FOV) due to the increase in sampling intervals of projected voxel locations on the detector array with increasing distance from the center of the FOV.
• FOV fields of view
• the sampling interval of projected voxel locations significantly exceeds the detector size, resulting in prominent aliasing artifacts.
• the present disclosure describes a CT scanner and scanning technique that employs a novel quasi-exact process derived from FDK, which handles arbitrary beamline geometry, reduces artifacts in scanners with large FOVs, and improves image quality while simultaneously enabling an efficient accelerated implementation on hardware such as graphics processing units (GPUs).
• GPUs graphics processing units
• the illustrated baggage scanning system indicated 100 includes a conveyor system 110 for continuously conveying baggage or luggage 112 in a direction indicated by arrow 114 through a central aperture of a CT scanning system 120 so that helical scans can be performed on the luggage.
• the conveyor system includes motor driven belts for supporting the baggage.
• Conveyer system 110 is illustrated as including a plurality of individual conveyor sections 122; however, other forms of conveyor systems may be used.
• the CT scanning system 120 includes an annular shaped rotating platform, or disk, 124 disposed within a gantry support 125 for rotation about a rotation axis 127 (shown in FIG. 3) that is preferably parallel to the direction of travel 114 of the baggage 112. Disk 124 is driven about rotation axis 127 by any suitable drive mechanism, such as a belt 116 and motor drive system 118, or other suitable drive mechanism.
• Rotating platform 124 defines a central aperture 126 through which conveyor system 110 transports the baggage 112.
• the rotation axis defines the Z-axis of the scanning system, while the X and F-axes (perpendicular to the Z-axis) are disposed in the center scanning plane (normal to the Z-axis).
• the system 120 includes an X-ray tube 128 and a detector array 130 which are disposed on diametrically opposite sides of the platform 124.
• the detector array 130 is preferably a two-dimensional array of any arbitrary geometry, as will be explained more fully hereinafter.
• the system 120 further includes a data acquisition system (DAS) 134 for receiving and processing signals generated by detector array 130, and an X-ray tube control system 136 for supplying power to, and otherwise controlling the operation of X-ray tube 128.
• DAS data acquisition system
• X-ray tube control system 136 for supplying power to, and otherwise controlling the operation of X-ray tube 128.
• the system 120 is also preferably provided with a computerized system (shown in FIG. 1 at 140) for processing the output of the data acquisition system 134 and for generating the necessary signals for operating and controlling the system 120.
• the computerized system 140 can also include a monitor for displaying information including generated images.
• System 120 also includes shields 138, which may be fabricated from
• the X-ray tube 128 includes at least one cathode and one anode for creating at least one separate focal spot from which an X-ray beam can be created and generated.
• the beam shown generally at 132 in FIGS. 1-3 passes through a three dimensional imaging field, through which conveying system 110 transports baggage 112. After passing through the baggage disposed in the imaging field, detector array 130 can receive each beam 132. The detector array then generates signals representative of the densities of exposed portions of baggage 112. The beams 132 therefore define a scanning volume of space.
• the data acquisition system 134 includes a processor subsystem within system 140 for carrying out the data processing described herein.
• the gantry can be represented in a coordinate system affixed to the imaged object 152.
• the beam from the X-ray source 150 is directed towards detector array 154 comprising rows of detectors 156.
• Each detector row 156 of the array 154 comprises detector channels 158. Intersection of a row and a channel defines one detector element (usually defined by the locus of the center of each detector).
• the x-ray source 150 traverses a helical trajectory 160 along the z-axis 162 relative to the imaged object 152 as shown.
• image reconstruction processes for helical CT scans can be classified into 2D and 3D processes.
• 3D processes produce fewer cone-beam artifacts when compared with the 2D processes.
• Relative image quality improvement achieved by 3D processes increases with scanner pitch and the number of detector rows.
• FIG. 5(A) A difference between 2D and 3D processes is illustrated in FIG. 5.
• approximated ID parallel projections 164 for a plurality of tilted planes 166 are extracted, indicated at 168. In doing so, one set of approximated parallel projections is generated for all voxels in the given tilted plane.
• the parallel projection data are then filtered.
• the 2D back-projection process is then used to reconstruct an image on tilted planes using a set of ID parallel data for a given tilted plane.
• 3D processes perform image reconstruction without conversion to intermediate 2D processes.
• One of the most widely used 3D image reconstruction processes is known as the Feldkamp Davis Kress (FDK) process.
• the FDK process is relatively easy to implement, stable, tractable, and is widely adopted by the industry.
• the FDK process includes filtering and 3D back-projection. First, the projection data are filtered along the detector rows. Second, 3D back-projection of the filtered data is performed. As shown in FIG. 5(B), in the 3D back-projection step, voxels 170 are individually and exactly projected directly in 3D space onto the detector array 172, filtered projection values from corresponding detector positions are accumulated to produce the density value for the reconstructed voxel 170.
• FIGs. 6 and 7 Aspects of the original FDK process for an axial scan are illustrated in FIGs. 6 and 7. As shown in FIGs. 6 and 7, the FDK process can be summarized as follows: The full rotation data 176 are weighted at 178 to account for different path lengths for different rays. Filtering is performed on data acquired along detector rows at 180, and then the filtered data are back projected at 182 voxels to create the backprojected image 184. During the backprojection step, the reconstructed voxel is projected onto the detector array for each view. Filtered projections are interpolated to obtain values where the projected voxel position falls between detectors as indicated at 186 in FIG. 6. The interpolated filtered projection values are then accumulated.
• the original FDK process has gained popularity due to its stability and ease of implementation. Numerous extensions of the process have been proposed to adapt the process to different scanning geometries and to improve image quality.
• FIG. 8(A) In the FDK process for the helical scan, in order to reconstruct a voxel, only ⁇ + ⁇ m (wherein ⁇ m is the detector array fan angle) worth of projection data (half- segment) is used, as illustrated in FIG. 8(A).
• the sampling region 194 of FIG. 8(B) for the half-segment set of data is represented in the space of detector angle and projection angle. If a given projection direction has two corresponding overlapping measurements (shown by way of example at 190 in FIG. 8(A), and in the shaded areas 192 in sampling region 194 of FIG. 8(B)), the data from rays projected in the two opposite directions must be weighted and summed to equalize contributions of all projection directions. Accordingly, as shown in FIG. 9 redundancy weighting was also introduced (shown at 196) prior to row filtering step 180 to take into account these redundant readings, which are weighted accordingly.
• Enhancements of the helical FDK process to reduce artifacts were recently proposed, including tilted plane reconstruction, (see, for example, Yan, M. and Zhang, C, "Tilted Plane Feldkamp Type Reconstruction Algorithm for Spiral Cone Beam CT", Med. Phys. 32 (11), 2005), and tangential filtering (see, for example, Sourbelle, K., and208, W.A., Generalization of Feldkamp reconstruction for Clinical Spiral Cone- Beam CT).
• the tilted plane reconstruction enhancement is as follows.
• the FDK process results in an exact reconstruction for a given plane when the source trajectory segment that is utilized lies in the plane. However, a constant speed helical source trajectory segment is not contained in any one plane. Inconsistency between the actual out-of-plane source trajectory and the in-plane assumption used in FDK process generates cone-beam artifacts. It was shown that by selecting planes of voxels to be reconstructed using closely fitted segments of the source trajectory, associated artifacts can be reduced.
• the tilted plane reconstruction approach will be used herein.
• Tangential filtering for flat detector arrays is generally illustrated in FIGs. 10 and 11.
• half-segment data 202 acquired from a flat detector array 200 (shown in FIG. 10) is filtered along tilted lines as shown at 204 so as to reduce artifacts.
• the half segment data 202 acquired from the flat detector array 200 is weighted at step 206 of the FDK process, and filtered at step 210 along the lines aligned with the helical source trajectory, followed by redundant data weighting at 208.
• the data are then back-projected at 212.
• FDK processes in the literature are usually designed for cylindrical detector arrays, which may not be applicable to modern CT scanners with unconventional (i.e., non-flat, non-cylindrical) beamline geometry.
• tangential filtering approach was proposed only for flat detector geometry. Therefore, in order to use the FDK process with tangential filtering enhancement in conjunction with unconventional detector geometry, the data must be resampled into a suitable detector geometry before reconstruction. This operation is performed while keeping resampled detector size as small as possible and reducing computation without compromising IQ.
• FDK processes in the literature generate aliasing artifacts in scanners with large fields of view (FOV) due to the increase in sampling interval of voxel locations projected on the detector array with increasing distanced from the center of the FOV.
• FOV fields of view
• scanners do not always include detector arrays that are made from one of two conventional geometries, i.e, flat and cylindrical, and in fact can physically be made in any one of a variety of shapes, particularly scanners used in security applications.
• a detector array of arbitrary geometry defined as having a detector array where data samples are not all equiangularly or equidistantly spaced
• resampling is taken at equiangular or equidistant positions, using for example a virtual array of one of two types of conventional geometry: flat detector arrays and cylindrical detector arrays.
• the actual detector array 220 can have any arbitrary non-flat geometry.
• the data acquired by the array 220 will be resampled at equiangular or equidistant spacings, as for example on a virtual flat detector array, or onto a virtual cylindrical detector array 226, in order to enable subsequent filtering step and to reduce artifacts.
• the data can be resampled onto tilted lines as shown at 228.
• the CT scanner employing the detector array with an arbitrary geometry is configured to carry out the exemplary steps of the process illustrated by the flow diagram of FIG. 14.
• the illustrated process receives half-segment data 232 and processes the data including the following steps: equiangular resampling with row dithering 234, fan-beam and cone-beam weighting 236, redundancy weighting 238, tangential filtering 240, adaptive upsampling 242, backprojection with sparse computation 244, and if necessary, as an added step Z-axis interpolation of the backprojected voxels 246.
• the detector data is resampled onto a virtual detector array (array 224 in FIG. 13).
• a virtual detector array array 224 in FIG. 13
• the detector data acquired with arbitrary detector geometry must be filtered along tilted along lines such as shown at 228.
• this filtering can be realized by resampling the data onto the tilted lines on a virtual flat detector array, such as shown at 224 (so as to provide equispatial samples).
• the data can be resampled onto curves in the virtual cylindrical detector array, which project onto required lines in the virtual flat detector, (so as to provide equiangular samples) all as shown by way of example in FIG. 16.
• equiangular resampling onto virtual cylindrical array
• the process is adapted for arbitrary beamline geometry, and required resolution is achieved with a minimum number of detector elements, while reducing computational cost.
• row dithering is used during step 234 (shown in FIG. 14).
• row dithering includes a shift in positions of the resampled detector rows for different views.
• resampled detectors for different detector angles have different bandwidths for different views, but on average have approximately the same effective bandwidth.
• resampled detectors do not exhibit bandwidth dependency upon the detector angle and the artifacts are eliminated.
• the detector sampling points using row dithering shown in FIG. 20(B) are compared to the sampling points of a prior art approach with no row dithering as shown in FIG. 20(A). The dithering approach in FIG.
• FIG. 20(B) is illustrated for simplicity using 2 views.
• the detector rows are shifted by half of the row width.
• the centers of the detectors are indicated by the intersection of the horizontal and vertical lines overlaid on the detector array.
• FIG. 20(A) high bandwidth is achieved where the sampling point lies at the center of a detector (such as detector 3 in FIG. 20(A)).
• data for detector 4 must be interpolated resulting in low bandwidth.
• the relative position of the detector rows and resampled curves is the same of all views.
• Each detector channel has a fixed bandwidth for all views, resulting in artifacts that add in phase during backprojection.
• the resampled detector data is weighted using a fan-beam and cone-beam weighting scheme, such as the one described in Wang, G. Lin, T., Cheng, P. and Schiozaki, D. M., "A General Cone-Beam Reconstruction Process", IEEE Trans. Med. Imaging, 12, 1993.
• the resampled and cone-beam weighted data is weighted to compensate for presence of redundant rays.
• one or two rays can pass through the voxel in a given direction.
• these two rays are called complementary rays.
• Complementary rays produce redundant data, which must be weighted and summed, to ensure that backprojection contributions from opposite sides of the segment transition smoothly into one another.
• a weighting scheme such as described in Parker, D.
• the 3D image is reconstructed by back-projecting the filtered data.
• Back-projection is carried out by projecting each voxel onto the resampled detector for each view in the half-segment while accumulating interpolated filtered projection values.
• FOVs fields of view
• the back-projection step results in aliasing artifacts that vary with location in the image. This is due to the fact that view angle interval between adjacent views strongly depends upon relative position of the reconstructed voxel and the X-ray source. Image locations close to the source show the largest aliasing artifacts (see FIG. 21).
• aliasing artifacts are corrected.
• the aliasing artifacts can be reduced by upsampling the voxel trajectory projected onto the detector. For example, as shown in FIG. 22, synthetic source positions 290 are created between each of two data views 292 and 294. The trajectory of the point projected onto the detector array is effectively upsampled. As shown at 296, the projection locations of the voxel corresponding to synthetic source positions are therefore generated. The backprojection sum is augmented with interpolated filtered projection values sampled at locations corresponding to synthetic source positions.
• the adaptive upsampling step 242 can be understood as performing filtering of the projections along the lines traced on the detector array surface by each reconstructed volume element, substantially eliminating aliasing artifacts.
• FIG 23 shows an example of adapting upsampling frequency as a function of image locations for large FOVs.
• Zones for two different views are shown in FIG. 23 by way of example, with zone 300 being closest to the source, zone 302 being of intermediate distance, and zone 304 farthest from the source. Voxels will fall into different zones for different views.
• the closest zone 300 will be the strong upsampling zone, for example, two additional locations 308 being provided between each location 310 (see FIG. 23(A)), while the intermediate zone 302 will be a moderate upsampling zone, for example, one additional location 312 being provided between each location 310 (see FIG. 23(B)), and zone 304 will be the no upsampling zone in which no additional locations are provided. Zones are thus rotated with the source S.
• a given voxel falls into different zones for different views. For any voxel, for any view, the resulting sampling interval is approximately the same. Zones with different upsampling are blended into one another in order to avoid discontinuities in the image.
• filtered projection data are backprojected onto the reconstructed image locations.
• An output driven process is utilized to reconstruct density values on tilted planes.
• a sparse computation approach can be implemented. Under this approach, projected locations for each voxel are computed exactly only for a set of equidistant key views. Projected locations for intermediate views are computed by interpolation from adjacent key views.
• FIG. 24 illustrates a simplified example of this backprojection step.
• the most computationally expensive part of back- projection is computing projected locations of voxels on the detector array.
• the projected locations are computed analytically only for key views. For example, the key views can be selected equispaced. Projected locations for intermediate views are determined by interpolation so that image quality is preserved.
• the projected locations of the voxel are computed exactly.
• the projected locations of the voxel are interpolated.
• the mathematically precise projected locations for intermediate views are shown at 332 along the true trajectory 334 of the projected voxel on the surface of the detector array.
• the locations 336 for the intermediate views are sampled along the approximating line 338, connecting projected locations for adjacent key views, shown at 330.
• CS Coordinate System
• ⁇ is the view angle and S be l t is the distance traveled by the belt in one rotation of the rotatable disk.
• ⁇ is the detector array tilt angle, given by the scanner design.
• the detector array tilt angle is the angle between the v-axis and the z-axis (detector is rotated around the s- axis of the axial gantry CS).
• R sd and R sc are source-to-detector and source-to-center distances respectively, given by the scanner design.
• detector channel coordinates in the ww-plane are given by the detector design.
• Detector channel u and w coordinates for the channel i are denoted u[i] and wfi] respectively.
• Detector angles y[i] for each channel ie [l..N ⁇ iet_DAs] are pre-computed and supplied to perform the operation of the scanner.
• the method is initialized by two preparation steps: (a) reconstruction geometry setup and (b) view segment selection.
• the reconstruction geometry setup step constructs the reconstruction image grid and computes the necessary parameters.
• the View segment selection step associates data view segments with reconstruction locations,
• the setup of reconstruction geometry includes specifying the reconstruction grid, resampled detector parameters, reconstruction planes, and other parameters utilized in subsequent steps of the method.
• One embodiment of the method includes the following:
• Image reconstruction grid is computed using image size parameters N x , N y , and N z , in a plane voxel size d pix , image center coordinates Xo, yo, and slice spacing d z .
• a set of reconstructed image voxels indices ⁇ is selected as follows:
• Detector elevation angles are computed.
• the elevation angle is defined for a given detector as the angle between the ray intersecting the detector and the plane containing the central row.
• the elevation angle is denoted as ⁇ (/,/).
• the elevation angle for each detector (rowy, channel /) is precomputed as follows:
• virtual filtering lines are defined in the virtual detector plane P. Virtual filtering lines are parallel and form filtering line angle, denoted as ⁇ , with t-axis. The angle ⁇ is equal to the tangent of the source trajectory:
• virtual filtering lines are optimal filtering lines required for tangential filtering.
• h f' is resampled detector height or the interval between filtering lines; and d rand ( ⁇ m °d Nj oop ) € [0..1] is the constant displacement for all filtering lines for the view m.
• Each filtering line in the virtual detector plane P forms a row of the resampled detector.
• N f channels for one row of the resampled detector.
• Detector channel angle ⁇ /'" of the channel / is given by: wherein wf' is the channel width in the resampled detector; and
• the resampled detectors are equiangular in the detector channels.
• Filtering lines computed using Eq. 10 are shown in FIG, 17 in detector coordinate system (detector angle ⁇ and row r).
• filtering lines computed using Eq. 10 are the tilted lines defined per equiangular resampling step with row dithering.
• Center detector for resampled detector s f n * is computed as follows:
• the half- fan angle ⁇ is computed as follows:
• Reconstruction planes with corresponding helical source trajectory segments are defined. Each reconstruction plane is defined by the center view of the associated helical trajectory segment. For each center view, the source lies in the corresponding reconstruction plane. The reconstruction plane approximately fits the helical source trajectory over the segment of 2N HS views.
• center views and corresponding reconstruction planes are defined as follows:
• Center views are selected with a view spacing of N vjn( [ ex views.
• the view index k m of the center view m is given by
• Reconstruction plane tilt angle y p (angle between tilted reconstruction plane and axial plane) is computed.
• y p defines two locations of the reconstruction plane intersection with the helical source trajectory, symmetric with respect to view v m , as illustrated in FIG. 18.
• the reconstruction plane R 1n for each center view m is defined as follows.
• the plane contains the source location P 1n corresponding to the view k m , given by (b)
• the plane normal vector n m is given by where ⁇ p is the reconstruction plane tilt angle computed in Eq. 15.
• the reconstruction plane Ri is illustrated in FIG. 18.
• set of reconstructed image voxels indices ⁇ m is selected as follows: where coordinates (x,y,z) for a voxel (i, j ) are computed as follows:
• the set of voxels defined by Eq. 19 defines the set of xy-grid voxels projected along the z- axis onto the tilted plane corresponding to the center view m.
• equiangular resampling with row dithering is employed, as illustrated in FIG. 20(B).
• the exemplary steps associated with equiangular resampling are as follows:
• elevation angle ⁇ is computed as where row coordinate rj(i) is computed using Eq. 10.
• Detector angle ⁇ f" computed using Eq. 9 is denoted as ⁇ 2.
• Real valued detector channel index c corresponding to the angle ⁇ is obtained using linear interpolation from the given set of detector angles ⁇ [ ⁇ ], i ⁇ [l..Ndet_DAs] as follows:
• Interpolated projection value P ⁇ (ij) for filtering line y and channel / is computed using bi-linear interpolation from the acquired detector data P m (.,.) as follows
• the process includes weighting and filtering steps. For each center view, resampled detector data for the views in the corresponding segment are weighted and filtered along rows. Weighting consists of half-scan redundant data weighting, and fan-beam cone-beam weighting. [0092] Fan-beam cone-beam weighting compensates for ray divergence in the fan-beam and in z-direction. Fan-beam weighting is given by cosine of the detector angle.
• Redundant data weighting is required because incomplete helical source trajectory segment and fan-beam geometry is used directly at the backprojection step. Depending on the location in the FOV, there is either one ray or two complementary rays passing through the reconstructed point. In case two rays emanating from the opposite sides of the helical source trajectory pass through the reconstructed point, corresponding samples must be weighted. Weighting ensures that two contributions add up to unity and that for any given reconstruction point, contributions from rays on opposite sides of the helical source trajectory transition smoothly into one another.
• the FDK process For each center view k m , the FDK process utilizes views from k m - N H s to k m + N HS to reconstruct the group of voxels associated with the center view k m . Each of the views is weighted and filtered. Weighting depends upon the view angle relative to the center view.
• ⁇ f' is the detector angle computed using Eq. 9 and W(k,k m i) is the
• Parker weight computed as follows. First, relative view angle ⁇ r is computed where ⁇ is half fan angle computed in Eq. 12. Parker weight W(k,k m i) is given by otherwise
• Row is zero-padded to length M p .
• Resulting row is denoted as Pfy ' ), i e [1..M p ].
• the backprojection step utilizes filtered views from fc m -N H s to fc m +N H s, to reconstruct the group of voxels ⁇ m associated with the center view k m . Selection of the group of voxels associated with each center view as described above. Each voxel in the group is projected into each of the filtered views. Interpolated filtered projection values are weighted and accumulated, yielding reconstructed voxel value.
• the backprojection step is the most computationally demanding step; therefore, in order to achieve required reconstruction rate, additional feature, sparse projection computation, is introduced.
• the concept of sparse projection computation is based on the fact that the projected location of a given reconstruction point onto the detector follow a smooth curve on the surface of the detector as a function of view angle. Therefore, it is sufficient to limit exact computation of projected coordinates of the point to only selected key views.
• the projected coordinates for intermediate views are computed by linearly interpolating projected coordinates computed for adjacent key views.
• Adaptive upsampling is integrated with sparse projection computation as follows. Projected coordinates of a reconstruction point are computed for the set of key views. Both projected coordinates of a point on intermediate views and projected coordinates corresponding to upsampled locations are computed by linear interpolation from coordinates computed for adjacent key view. The principle is illustrated in FIG. 19.
• the voxel (x, y, z) in the reconstruction coordinate system is transformed into the axial gantry CS using eq. 1.
• the voxel coordinates in the axial gantry CS are denoted as (s, t, p).
• Source-to-voxel distance L in the rf-plane is computed as
• Resampled detector real valued row coordinate is computed as
• Backprojection weight W k is computed as follows
• the backprojection step is carried out as follows. For a given center view k m , for each voxel ⁇ ' • > ' from the group m (coordinates ( x i Vi z ) computed in Eq. 4), the reconstructed value v 1 ⁇ ' Mn the tilted image is obtained as follows:
• detector indices kin and t kin and view index ⁇ -* are computed as
• the z-interpolation step performs interpolation in z-direction between tilted planes to form untilted image.
• the density value is computed as follows:
• the bottom adjacent tilted plane index ff ⁇ is found as follows m

Landscapes

• Physics & Mathematics (AREA)
• General Physics & Mathematics (AREA)
• Engineering & Computer Science (AREA)
• Theoretical Computer Science (AREA)
• Apparatus For Radiation Diagnosis (AREA)
• Analysing Materials By The Use Of Radiation (AREA)
• Algebra (AREA)
• Mathematical Analysis (AREA)
• Mathematical Optimization (AREA)
• Mathematical Physics (AREA)
• Pure & Applied Mathematics (AREA)

Priority Applications (4)

Applications Claiming Priority (1)

Publications (1)

Family

ID=41820542

Family Applications (1)

Country Status (4)

Cited By (2)

Families Citing this family (11)

Citations (2)

Family Cites Families (16)

• 2009
  • 2009-06-30 WO PCT/US2009/049138 patent/WO2011002442A1/en not_active Ceased
  • 2009-06-30 JP JP2012518516A patent/JP5507682B2/ja not_active Expired - Fee Related
  • 2009-06-30 DE DE112009005019.0T patent/DE112009005019B4/de active Active
  • 2009-06-30 US US13/318,653 patent/US9665953B2/en active Active

Patent Citations (2)

Non-Patent Citations (5)

Also Published As

Similar Documents

Legal Events