Method of reconstructing images for spiral and non-spiral computer tomography
View Patent ↗General scheme processes and systems for constructing algorithms for reconstructing images of objects that have been scanned in a spiral or non-spiral fashions with detectors. Application of the scheme requires finding of a weight function, which would lead to the required reconstruction algorithm. This general scheme can use a C-arm scan with the closed x-ray source trajectory and gives a new, theoretically exact and efficient (i.e., with the convolution-based FBP structure) reconstruction algorithm. The invention can also utilize the algorithms disclosed in an earlier application U.S. patent application Ser. No. 10/143,160 filed May 10, 2002, entitled: Exact Filtered Back Projection (FBP) Algorithm For Spiral Computer Tomography, which claims the benefit of U.S. Provisional Application No. 60/312,827 filed Aug. 16, 2001, also fit into the general scheme.
1. A method of computing images derived from computer tomography detectors, comprising the steps of:
(a) collecting cone beam data from a detector during a scan of an object;
(b) preprocessing the cone beam data, wherein the step of preprocessing includes
(bi) computing the derivative of the said data along the source trajectory
∂
∂
q
D
f
(
y
(
q
)
,
Θ
)
q
=
s
,
where
D ƒ is the cone beam data,
q,s are variables along the source trajectory,
y(q) is the source position as a function of q,
Θ is a unit vector;
(c) identifying lines on a plane DP(s), where s is the parameter along the scan path and DP(s) is any plane intersecting the cone beam;
(e) convolving the said preprocessed data along said lines with filter 1/sin(γ), where γ is polar angle in the plane containing the current source position and a said line;
(f) back projecting said filtered data to form a precursor of said image; and
(g) repeating steps a, b, c, d, e, and f until an image of the object is completed.
2. A method of computing images derived from computer tomography detectors, comprising the steps of:
(a) collecting cone beam data from a detector during a scan of an object;
(b) identifying lines on a plane Π intersecting the cone beam, the said lines being filtering lines,
(c) preprocessing said data along said filtering lines, wherein the step (c) of preprocessing includes computing the derivative of D ƒ (y(s),Θ) with respect to Θ along a direction non-parallel to the plane determined by y(s) and a filtering line, the said plane being a filtering plane, here
s is parameter along the scan path, which determines point y(s) on the said path,
D ƒ (y,Θ) is the cone beam transform of ƒ corresponding to the x-ray source located at the point y and the direction Θ,
ƒ is a function describing the object being scanned;
(d) shift invariant filtering said preprocessed data along said filtering lines, wherein the step (d) of filtering includes
(i) convolving the data D ƒ (y(s),Θ) with kernel 1/sin(γ) within a filtering plane, where γ is polar angle in the plane,
(ii) convolving the derivative of D ƒ (y(s),Θ) with a kernel cot(γ) within a filtering plane, the derivative of D ƒ (y(s),Θ) is the derivative with respect to Θ along a direction non-parallel to the filtering plane,
(e) back projecting said filtered data to form a precursor of said image; and
(f) repeating steps a, b, c, d, and e until an image of the object is reconstructed.