CN103679768B - 以加权反投影重建ct图像数据的方法、计算单元和系统 - Google Patents
以加权反投影重建ct图像数据的方法、计算单元和系统 Download PDFInfo
- Publication number
- CN103679768B CN103679768B CN201310435145.XA CN201310435145A CN103679768B CN 103679768 B CN103679768 B CN 103679768B CN 201310435145 A CN201310435145 A CN 201310435145A CN 103679768 B CN103679768 B CN 103679768B
- Authority
- CN
- China
- Prior art keywords
- detector
- weighting function
- projection data
- row
- weighting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Pulmonology (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
Abstract
本发明涉及一种用于通过加权滤波反投影重建由多个图像体素(Vx,y,z)组成的CT图像数据的方法,其中,为了对滤波反投影中的投影数据(Pconv)进行加权,使用加权函数所述加权函数将在行边缘附近击中检测器行的中央射线的投影数据更小地加权,并且将中央地击中检测器行的中央射线的投影数据更强地加权。
Description
技术领域
本发明涉及一种用于通过加权滤波反投影(=weighted filtered back-projection=WFBP)重建由多个图像体素组成的CT图像数据的方法,其中从至少一个围绕检查对象旋转的、其检测器行在圆周方向上垂直于CT系统的系统轴延伸的多行检测器接收检测器数据,从相邻射线的检测器数据插值投影数据;其中对于每个预先给出的投影角,与每个图像体素对应一个通过该图像体素的中央射线的投影数据;并且其中为每个投影数据分配与检测器相关的检测器行坐标,所述坐标对应于检测器上的中央射线的位置,并且根据为每个投影数据分配的检测器行坐标进行投影数据的加权。此外,本发明还涉及一种用于执行该重建方法的计算单元和CT系统。
背景技术
这样的重建方法一般作为WFBP重建公知。与此相关地,例如参见Stierstorfer etal.,Phys.Med.Biol.49pp.2209-2218,2004,其中公开了通过加权滤波反投影的这样的重建。
这样的WFBP方法在大多数情况下由四个步骤构成,其中出于实际原因通常使用锥形射线检测器数据到平行数据的转换,以便在计算技术上简化实际的重建。
在第一步骤中利用CT系统测量检测器原始数据,其中,在此使用的离散坐标与连续坐标具有以下关系:
β=β0+kβΔβ,
q=q0+kqΔq,
α=α0+kαΔα,
其中β描述了测量射线的扇形角,q描述了检测器通道,并且α描述了焦点的角度位置,并且β0,q0,α0分别描述了测量的开始并且Δβ,Δq,Δα分别描述了测量步骤的增量。
在WFBP方法的第二步骤中将存在的锥形射线检测器数据Praw[kβ,kq,kα]重整为重整的平行投影数据Preb[kp,kq,kθ]。在此从离散的锥形射线坐标β,q,α过渡到离散的平行射线坐标p,q,θ,其中p是通道索引,q是检测器行坐标,并且θ是平行投影角。
在第三步骤中利用斜坡滤波器将重整的平行投影数据Preb[kp,kq,kθ]卷积,并且形成卷积的投影数据Pconv(p,q,θ),其中,在此不改变数据的几何特征。
在第四并且是最后的步骤中将卷积的投影数据反投影到待重建的体积。这在数学上可以通过公式
来描述,其中Vx,y,z对应于CT图像的笛卡尔坐标x、y和z上的重建的图像体素的值,Pconv(p,q,θ+kπ)是在平行坐标上的卷积的并进行了平行重整的投影数据,其中k是整数,其反映了在扫描时检测器的半周旋转数量,检测器行坐标q描述了在点x,y,z处的图像体素Vx,y,z在系统轴方向上投影到检测器上的中央射线,具有关于检测器宽度标准化的、在-1和+1之间的值,并且W(q)是取决于检测器行坐标的加权函数,其为了降低伪影将检测器边缘侧的行向着边缘的方向降低地加权。
原则上在这样的方法中存在问题,即,图像体素到扫描的检测器上的投影仅很少击中检测器行的中心。由此必要的是,首先从相邻的检测器数据中,必要时在事先重整为平行数据之后,通过插值计算投影数据,利用所述投影数据进行WFBP重建。在特别简单的变形方案中,在此仅使用来自于直接相邻的检测器行的检测器数据。
在此的缺陷是,通过从检测器数据插值投影数据,形成的图像的空间分辨率在z方向上变差。
发明内容
由此,本发明要解决的技术问题是,找到一种可能性来改善由于检测器数据的插值而变差的投影数据的空间分辨率,并由此改善从中重建的图像数据。
发明人已经认识到,可以通过如下改善重建的图像数据的分辨率,即,对为了WFBP重建而使用的、与相对的检测器行坐标对应的投影数据进行加权,该加权将其检测器坐标更靠近检测器行的中心的投影数据比其检测器坐标远离检测器行的中心的投影数据更强地加权。
相应地得到取决于检测器行坐标的加权函数,其交替地关于检测器行延伸,并且在此,其检测器行坐标靠近行中心的投影数据被更强地加权,而其检测器行坐标远离行中心的投影数据被更弱地加权。该加权函数优选也可以与将在z方向上靠近边缘的检测器行比位于中央的检测器行更小地加权的公知的加权函数叠加,从而按照趋势,边缘附近的检测器行被更小地加权,而在行宽内部,加权的最大值在行中央处。也就是公知的加权函数通过随着行宽交替的加权叠加。
该按照本发明的过程的一个主要优点在于,已经存在的、硬件方面实现的反投影可以无需改变就用于FBP重建,因为仅必须对传输到那里的投影数据进行按照本发明的加权。
相应地,发明人建议一种用于改善CT图像数据的重建的方法,其中CT图像数据由多个图像体素组成,并且通过加权滤波反投影(=weighted filtered backprojection=WFBP)进行重建,具有以下方法步骤:
-从至少一个围绕检查对象旋转的、其检测器行在圆周方向上垂直于CT系统的系统轴延伸的多行检测器接收检测器数据,
-从相邻射线的检测器数据计算卷积的投影数据,其中对于每个预先给出的投影角,与每个图像体素分别对应一个通过该图像体素的中央射线的投影数据,并且其中为每个投影数据分配与检测器相关的检测器行坐标,所述坐标对应于检测器上的中央射线的位置,
-根据为每个投影数据分配的检测器行坐标进行投影数据的加权。
关于相对检测器行坐标的确定,例如参见已经引用的文献Stierstorfer et al.,Phys.Med.Biol.49pp.2209-2218,2004,其中在那里在等式(7)和(14)中描述了在此使用的相对检测器行坐标的计算。
该重建方法的按照本发明的改进在于,为了在滤波反投影中对投影数据的加权,使用加权函数,所述加权函数将在行边缘附近击中检测器行的中央射线的投影数据更小地加权,并且将中央地击中检测器行的中央射线的投影数据更强地加权。
由此有针对性地将更好地击中检测器行的中心的投影数据比其检测器行坐标位于两个检测器行之间的那些投影数据更高地评价。也就是进行投影数据的取决于位置的加权,在该加权中,在检测器行宽的区域中的距离导致关于使用的权重的强烈变化。
虽然最后描述的方法原则上针对FBP重建的改进,但是在特别优选的实施方式变形中也可以改进WFBP重建方法,方法是,附加地这样构造加权函数,使得通过加权函数对投影数据的加权当逐行地来考虑时,将具有检测器边缘附近的检测器行坐标的投影数据总体上更小地加权,而将具有检测器中央的检测器行坐标的投影数据更强地加权。
尽管建议的按照本发明的方法在使用锥形射线重建的条件下工作,但是主要建议,在WFBP重建之前进行重整,据此然后所有使用的投影按照平行射线几何特征呈现,由此按照FBP方法更简单地进行实际的重建。
该按照本发明的技术相对于已知的使用共轭射线对的插值相比的优点在于,不仅自动地使用两个相邻的射线的信息,而且使用所观察的体素附近的任意数量的射线的信息。由此不仅将冗余的信息用于改善分辨率,而且也用于改善信噪比。
具体来说,为了按照WFBP方法重建图像体素使用以下公式:
其中:
-Vx,y,z对应于重建的体素在笛卡尔坐标x,y和z处的值,
-Pconv(p,q,θ+kπ)对应于重整并卷积的投影数据,
-θ对应于在0和180°之间的平行投影角,
-k对应于反映了在扫描时检测器的半周旋转的数量的整数,
-q对应于在点x,y,z处的体素Vx,y,z在系统轴方向上投影到检测器上的中央射线的检测器行坐标,具有关于检测器宽度标准化的、在-1和+1之间的值,和-对应于在WFBP重建中的取决于检测器行坐标的加权函数。
此外发明人建议,如下修改本身公知的WFBP重建方法(在所述WFBP重建方法中已经进行了全体检测器行的加权,该加权将边缘行比位于中央的行更少地加权),使得
-为了加权在滤波反投影中的投影数据使用加权函数,所述加权函数是第一加权函数和第二加权函数的乘积,其中:
-第一加权函数将位于检测器中在系统轴方向上靠近边缘的检测器行的投影数据向着边沿方向比位于里面的检测器行的投影数据越来越少地加权,和
-第二加权函数具有在系统轴方向上关于检测器行的多次交替的值曲线。
由此可以非常简单地在已经实施了WFBP方法的并且为此具有优选由硬件实现的WFBP投影器的现有CT系统上实现该方法,因为仅需将按照本发明的小结构的行加权,也就是关于行的梳状的加权,与已经存在的对检测器全面进行的、具有检测器边缘侧降低的权重的加权相叠加。
在此特别有利的是,第二加权函数的交替的值曲线分别在检测器行的中心具有最大值,其中该最大值可以形成曲线平坦区,和/或第二加权函数的交替的值曲线分别在检测器行之间具有最小值。也就是例如对于第二加权函数可以使用正弦函数,其最大值和最小值分别与行中心和行边缘同步。
在此可以按照公式计算加权函数的叠加,其中成立:
W(q)=第一加权函数,
D(q)=第二加权函数,q=系统轴方向上的检测器行的坐标。
此外第一加权函数可以具有值曲线,对于所述值曲线成立:
对于|q|<Q,W(q)=1
对于Q≤|q|<1,
对于q的所有其他值,W(q)=0,
其中q是系统轴方向上的检测器行的坐标,并且Q是在0和1之间的预先给出的恒定值,其确定了加权函数的曲线平坦区宽度。
第二加权函数的值曲线的示例计算可以通过以下公式进行:
并且:
其中:
Nq对应于在使用的检测器中的检测器行的数量,
参数Wm可以取在0和1之间的值,对应于对于D(q)的最小值,并且
用于定义峰值宽度的参数rw可以取在0和0.5之间的值。
通过常数rw的变化例如可以从事先确定的检测器数据中计算具有不同的大的峰值宽度的不同图像数据。在此,在更大的峰值宽度的情况下产生更好的剂量使用,并且随着峰值宽度减小,对强加权的投影数据产生空间上更精确的选择。
为了按照WFBP方法重建体素,此外还建议使用以下公式:
其中:
-Vx,y,z对应于笛卡尔坐标x,y和z处的重建的体素的值,
-Pconv(p,q,θ+kπ)对应于重整并卷积的投影数据,
-θ对应于在0和180°之间的平行投影角,
-k对应于反映了在扫描时检测器的半周旋转的数量的整数,
-q对应于在系统轴方向上在点x,y,z处的体素Vx,y,z的投影到检测器上的检测器行坐标,具有关于检测器宽度标准化的、在-1和+1之间的值,和
-W(q)对应于第一加权函数,
-D(q)对应于第二加权函数。
在此对于每个投影,对于WFBP重建的加权可以独立于来自于检测器的其他旋转的其他投影来进行。
特别有利地可以将按照本发明的方法应用于来源于CT螺旋扫描的检测器数据,特别是当使用产生重叠的扫描区域的进给速度时。
但是替换地,按照本发明的方法也可以结合来源于CT圆扫描的检测器数据进行。在此,圆扫描也可以划过大于360°,或者利用多个重叠的扫描区域进行。
此外指出,也可以结合迭代的重建方法应用按照本发明的方法。
相应地,发明人还建议一种用于迭代重建由多个图像体素组成的CT图像数据的方法,其中:
-在迭代中对卷积的投影数据进行至少一次加权反投影,
-从相邻射线的检测器数据计算卷积的投影数据,并且
-对于每个预先给出的投影角,与每个图像体素分别对应一个通过该图像体素的中央射线的投影数据,
-并且为每个投影数据分配与检测器相关的检测器行坐标,所述坐标对应于检测器上的中央射线的位置。
该方法的按照本发明的改进通过如下实现,即,为了投影数据的加权,使用加权函数,所述加权函数将在行边缘附近击中检测器行的中央射线的投影数据更小地加权,并且将中央地击中检测器行的中央射线的投影数据更强地加权。
有利地,可以通过加权函数这样进行对投影数据的加权,使得当逐行地来考虑时,具有检测器边缘附近的检测器行坐标的投影数据总体上被更小地加权,而具有检测器中央的检测器行坐标的投影数据被更强地加权。
此外还建议:
-为了对卷积的投影数据加权,使用加权函数,所述加权函数是第一加权函数和第二加权函数的乘积,其中:
-第一加权函数将位于检测器中在系统轴方向上靠近边缘的检测器行的投影数据向着边沿方向比位于里面的检测器行的投影数据越来越少地加权,和
-第二加权函数具有在系统轴方向上关于检测器行的多次交替的值曲线。
关于第二加权函数的交替的值曲线的构造建议,其分别在检测器行的中心具有最大值,或替换地,第二加权函数的交替的值曲线分别在检测器行的中心比检测器行的边缘具有更大的值。
除了按照本发明的方法之外,属于本发明范围的还有一种用于执行WFBP重建的计算单元,具有至少一个处理器、用于计算机程序的存储器和输出单元,其中在程序存储器中存储了计算机程序,其执行按照本发明的方法的方法步骤。
此外,计算单元还可以具有至少一个可手动操作的调节器,通过所述调节器可以在每次图像计算之前调整至少一个参数,所述参数定义了至少一个加权函数的形式。由此,系统的用户可以利用单独选择的参数多次进行图像计算,并且通过视觉上观察所得到的图像结果,这样找到在剂量利用和使用的投影数据的宽度之间的最佳折衷。
属于本发明的还有一种具有上面描述的计算单元的CT系统。
附图说明
以下借助附图详细解释本发明,其中仅示出对于理解本发明所需的特征。使用了以下附图标记和变量:1:CT系统;2:第一X射线管;2.1:焦点;3:第一多行检测器;4:第二X射线管;5:第二多行检测器;6:机架壳体;6.1:机架;7:患者;8:患者卧榻;9:系统轴;10:计算单元;q:检测器行坐标;M:中线;Prg1–Prgn:计算机程序;S:射线束;W(q):第一加权函数;D(q):第二加权函数;:新的加权函数;WQ1(q)-WQ4(q):第一加权函数的曲线。各个图示出:
图1示出具有用于执行按照本发明的方法的计算单元的CT系统;
图2示出机架中的多行检测器的示意图;
图3示出8行检测器的俯视图;
图4示出分别对应检测器的半周旋转位置和焦点在k位置处的投影的射线的示意图;
图5示出关于多行检测器的行的第一加权函数W(q)的示例曲线,其对于位于边缘处的行具有减小的权重;
图6示出关于多行检测器的行的第二加权函数D(q)的示例曲线,具有行中线处的最大值和行边缘处的最小值的关于行交替的加权;
图7示出图5和6的第一和第二加权函数通过相乘的叠加;
图8示出图5的第一加权函数与具有更窄的峰值的第二加权函数的叠加;
图9示出对具有不同直径的多个圆柱体结构的模体的按照现有技术的重建结果和按照本发明的重建结果的对比显示。
具体实施方式
为了详细解释本发明在图1中示出具有机架壳体6的CT系统1,其中在机架上具有多行检测器3以及对置的X射线管2。从X射线管2出发或者更确切来说从在那里产生的焦点出发,在机架的旋转期间,射线束指向对置的检测器,并且借助患者卧榻8将患者7顺序地或优选连续地对于螺旋扫描在系统轴9的方向上移动通过射线束的测量场,并且测量射线的衰减。可选地,为了扫描还可以同时使用第二射线束,所述第二射线束从第二X射线管4指向角度偏移地布置在机架上的第二多行检测器5。如果在两个管-检测器-系统中使用相同的X射线谱,则可以将所确定的检测器数据共同地用于图像数据的重建。
CT系统1通过计算单元10控制,在所述计算单元中存储了计算机程序Prg1–Prgn,所述计算机程序在运行中也执行按照本发明的重建方法等。
图2此时示出了从具有机架6.1的CT系统的示意性三维截图,所述机架布置在前面示出的机架壳体中。示出了具有在X射线管中布置的焦点2.1的辐射器-检测器-系统,从所述焦点出发,射线束S指向具有检测器行的多行检测器3,所述检测器行沿着圆周方向对应于画出的坐标延伸。检测器行的位置在此通过检测器行坐标q确定,其通常平行于CT系统的系统轴取向。当机架不倾斜时,系统轴同样对应于画出的CT系统的笛卡尔坐标系的z轴。
如果在从焦点俯视观察按照平面图所示的图1和2的多行检测器3,则得到图3中的图示。在此示出了具有Nq=8个检测器行和24个检测器通道的示例性多行检测器3。在沿着轴延伸的检测器行中,虚线地示出了检测器行的中线M。此外以在-1和+1之间的值标出了优选关于检测器宽度标准化的检测器行坐标q。中线M相应地具有检测器行坐标q∈{-1+1/Nq,-1+3/Nq,-1+5/Nq,....,1-1/Nq,},而在检测器行之间的边界取检测器行坐标q∈{-1,-1+2/Nq,-1+4/Nq,....,1}。
当通过多行检测器进行螺旋形扫描时,如果观察沿着z轴的片段,则得到在图4中示出的、在位置k处的焦点和对置的检测器之间的射线的图像。
相对的检测器行坐标如示例性在已经引用的文献Stierstorfer et al.,Phys.Med.Biol.49pp.2209-2218,2004中示出的那样计算。那里也给出了用于计算在此使用的相对检测器行坐标的等式(7)和(14)。
为了基于已经存在的、对检测器行进行检测器边缘附近降低加权的第一加权函数来描述本发明,可以使用如在图5中示出的第一加权函数W(q)。该图示出了第一加权函数W(q)的曲线,对于其成立:
对于|q|<Q,W(q)=1
对于Q≤|q|<1,
对于q的所有其他值,W(q)=0,
其中q是系统轴方向上的检测器行的坐标,并且Q是在0和1之间的预先给出的恒定值,其确定了加权函数的曲线平坦区宽度。
示出了具有对于参数Q的不同值的五条曲线WQ1(q)至WQ5(q)。
按照本发明可以将第一加权函数与第二加权函数叠加,后者用于保证,从存在的投影数据的总和中,将FBP重建中的以下投影数据加权,所述投影数据具有与理论上的中央射线的最佳空间一致性,所述中央射线通过待重建的图像体素确定。相应地得到关于确定的投影数据的检测器行坐标的加权函数,其交替地分别在检测器行的中线处和附近具有最大值,并且在从检测器行到检测器行的过渡处具有最小值。在图6中示例性示出具有相对大的峰值宽度的这样的第二加权函数D(q)的曲线。
第二加权函数D(q)在此遵循值曲线,对于所述值曲线成立:
并且:
其中:
Nq对应于在使用的检测器中的检测器行的数量,
参数Wm可以取在0和1之间的值,对应于对于D(q)的最小值,并且
用于定义峰值宽度的参数rw可以取在0和0.5之间的值。
在图6具体示出的情况下,选择了参数rw=0.25和Wm=0.1,其中假定是具有Nq=16的16行多行检测器。
按照本发明可以将第一和第二加权函数以相乘W(q)*D(q)的形式叠加,从而总体上得到具有加权的边缘侧低权重的梳子形(kammartiger)的值曲线,如在图7中示出的。示例性地在图8中基于相同的第一加权函数和具有窄的峰值的第二加权函数,对应于改变的参数rw=0.1,示出总加权W(q)*D(q)的另一个曲线,其导致空间上更精确的图像数据。
最后在图9中示出了仅采用第一加权函数的WFBP方法的重建结果(左侧A列)和利用将第一加权函数附加地与按照第二加权函数的梳状加权(Kammgewichtung)叠加的WFBP方法的重建结果(右侧B列)的对比图示。该图示在此示出具有在0.5mm和3.0mm之间的不同直径的多个圆柱体结构的模体的扫描。从对比示出的重建结果中可以清楚看出,在按照本发明的重建中明显改善了空间分辨率和细节显示。
总体上本发明描述了一种用于通过加权滤波反投影来重建由多个图像体素组成的CT图像数据的方法、计算单元和CT系统,其中,为了加权在滤波反投影中的投影数据,使用加权函数,所述加权函数将在行边缘附近击中检测器行的中央射线的投影数据更小地加权,并且将中央地击中检测器行的中央射线的投影数据更强地加权。
尽管详细地通过优选实施例详细示出并描述了本发明,但是本发明不限于公开的例子,并且专业人员可以从中导出其他变形,而不脱离本发明的保护范围。
Claims (24)
1.一种用于通过加权滤波反投影重建由多个图像体素(Vx,y,z)组成的CT图像数据的方法,具有以下方法步骤:
-从至少一个围绕检查对象旋转的、其检测器行在圆周方向上垂直于CT系统(1)的系统轴(9)延伸的多行检测器(3,5)接收检测器数据,
-从相邻射线的检测器数据计算卷积的投影数据(Pconv),其中对于每个预先给出的投影角,与每个图像体素(Vx,y,z)分别对应一个通过该图像体素的中央射线的投影数据,并且其中为每个投影数据分配与所述检测器有关的检测器行坐标(q),所述坐标对应于检测器上的中央射线的位置,
-根据为每个投影数据所分配的检测器行坐标(q)进行所述投影数据(Pconv)的加权,
其特征在于,为了对所述滤波反投影中的所述投影数据(Pconv)进行加权,使用加权函数所述加权函数将在行边缘附近击中检测器行的中央射线的投影数据(Pconv)更小地加权,并且将中央地击中检测器行的中央射线的投影数据更强地加权,
其中,
-为了加权在所述滤波反投影中的卷积的投影数据(Pconv),使用加权函数所述加权函数是第一加权函数和第二加权函数的乘积,其中:
-所述第一加权函数(W(q))将位于检测器中在系统轴方向上靠近边缘的检测器行的投影数据向着边沿方向比位于里面的检测器行的投影数据越来越少地加权,和
-所述第二加权函数(D(q))具有在系统轴(9)方向上关于检测器行的多次交替的值曲线。
2.根据上述权利要求1所述的方法,其特征在于,通过所述加权函数对所述投影数据的加权当逐行地来考虑时,将具有检测器边缘附近的检测器行坐标的投影数据总体上更小地加权,而将具有检测器中央的检测器行坐标的投影数据更强地加权。
3.根据上述权利要求1至2中任一项所述的方法,其特征在于,在WFBP重建之前进行重整。
4.根据权利要求1或2所述的方法,其特征在于,对于每个投影,对于WFBP重建的加权能够独立于来自于检测器(3,5)的其他半周旋转的其他投影来进行。
5.根据权利要求1或2所述的方法,其特征在于,为了按照WFBP方法重建图像体素(Vx,y,z),使用以下用于反投影的公式:
其中:
-Vx,y,z对应于笛卡尔坐标x,y和z处的重建的体素的值,
-Pconv(p,q,θ+kπ)对应于重整并卷积的投影数据,
-θ对应于在0和180°之间的平行投影角,
-k对应于反映了在扫描时检测器的半周旋转的数量的整数,
-q对应于在系统轴方向上在点x,y,z处的体素Vx,y,z的投影到检测器上的中央射线的检测器行坐标,具有关于检测器宽度标准化的、在-1和+1之间的值,和
-对应于在WFBP重建中的取决于检测器行坐标的加权函数。
6.根据上述权利要求1所述的方法,其特征在于,所述第二加权函数(D(q))的交替的值曲线分别在检测器行的中心具有最大值。
7.根据上述权利要求6所述的方法,其特征在于,所述最大值是曲线平坦区。
8.根据权利要求1所述的方法,其特征在于,所述第二加权函数(D(q))的交替的值曲线分别在检测器行之间具有最小值。
9.根据权利要求1所述的方法,其特征在于,按照公式计算所述加权函数的叠加,其中成立:
W(q)=第一加权函数,
D(q)=第二加权函数,
q=系统轴方向上的检测器行的坐标。
10.根据权利要求1所述的方法,其特征在于,所述第一加权函数(W(q))具有值曲线,对于所述值曲线成立:
对于|q|<Q,W(q)=1
对于Q≤|q|<1,
对于q的所有其他值,W(q)=0,
其中q是系统轴方向上的检测器行的坐标,并且Q是在0和1之间的预先给出的恒定值,其确定了加权函数的曲线平坦区宽度。
11.根据权利要求1所述的方法,其特征在于,为了按照所述用于反投影的WFBP方法重建体素(Vx,y,z),使用以下公式:
其中:
-Vx,y,z对应于笛卡尔坐标x,y和z处的重建的体素的值,
-Pconv(p,q,θ+kπ)对应于重整并卷积的投影数据,
-θ对应于在0和180°之间的平行投影角,
-k对应于反映了在扫描时检测器的半周旋转的数量的整数,
-q对应于在系统轴方向上在点x,y,z处的体素Vx,y,z的投影到检测器上的中央射线的检测器行坐标,具有关于检测器宽度标准化的、在-1和+1之间的值,和
-W(q)对应于所述第一加权函数,
-D(q)对应于所述第二加权函数。
12.根据权利要求1所述的方法,其特征在于,对于每个投影,对于WFBP重建的加权能够独立于来自于检测器的其他半周旋转的其他投影来进行。
13.根据权利要求1或2所述的方法,其特征在于,所述检测器数据来源于CT螺旋扫描。
14.根据权利要求13所述的方法,其特征在于,使用了产生重叠的扫描区域的进给速度。
15.根据权利要求1或2所述的方法,其特征在于,所述检测器数据来源于CT圆扫描。
16.根据上述权利要求15所述的方法,其特征在于,这样执行所述CT圆扫描,使得产生重叠的扫描区域。
17.根据权利要求1或2所述的方法,其特征在于,所述方法在迭代的重建方法中被应用。
18.一种用于迭代重建由多个图像体素(Vx,y,z)组成的CT图像数据的方法,其中:
-在迭代中对卷积的投影数据(Pconv)进行至少一次加权反投影,
-从相邻射线的检测器数据计算所述卷积的投影数据(Pconv),并且
-对于每个预先给出的投影角,与每个图像体素(Vx,y,z)分别对应一个通过该图像体素的中央射线的投影数据,
-并且为每个投影数据分配与检测器有关的检测器行坐标(q),所述坐标对应于检测器上的中央射线的位置,
其特征在于,为了所述投影数据(Pconv)的加权,使用加权函数所述加权函数将在行边缘附近击中检测器行的中央射线的所述投影数据(Pconv)更小地加权,并且将中央地击中检测器行的中央射线的所述投影数据更强地加权,
其中,
-为了对卷积的投影数据(Pconv)加权,使用加权函数所述加权函数是第一加权函数和第二加权函数的乘积,其中:
-所述第一加权函数(W(q))将位于检测器中在系统轴方向上靠近边缘的检测器行的投影数据向着边沿方向比位于里面的检测器行的投影数据越来越少地加权,和
-所述第二加权函数(D(q))具有在系统轴(9)方向上关于检测器行的多次交替的值曲线。
19.根据上述权利要求18所述的方法,其特征在于,通过所述加权函数对所述投影数据的加权,使得当逐行地来考虑时具有检测器边缘附近的检测器行坐标的投影数据总体上被更小地加权,而具有检测器中央的检测器行坐标的投影数据被更强地加权。
20.根据上述权利要求18所述的方法,其特征在于,所述第二加权函数的交替的值曲线分别在检测器行的中心具有最大值。
21.根据上述权利要求18所述的方法,其特征在于,所述第二加权函数(D(q))的交替的值曲线分别在检测器行的中心比检测器行的边缘具有更大的值。
22.一种用于执行WFBP重建的计算单元(10),具有至少一个处理器、用于计算机程序的存储器和输出单元,其特征在于,在程序存储器中存储了计算机程序(Prg1–Prgn),其执行按照上述权利要求1至21中任一项所述的方法的方法步骤。
23.根据上述权利要求22所述的计算单元(10),其特征在于,具有至少一个可手动操作的调节器,通过所述调节器能够在每次图像计算之前调整至少一个参数,所述参数定义了至少一个加权函数的形式。
24.一种具有按照上述权利要求23所述的计算单元(10)的CT系统(1)。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102012217163.0A DE102012217163B4 (de) | 2012-09-24 | 2012-09-24 | Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren |
DE102012217163.0 | 2012-09-24 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103679768A CN103679768A (zh) | 2014-03-26 |
CN103679768B true CN103679768B (zh) | 2017-04-12 |
Family
ID=50317215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310435145.XA Expired - Fee Related CN103679768B (zh) | 2012-09-24 | 2013-09-23 | 以加权反投影重建ct图像数据的方法、计算单元和系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9299140B2 (zh) |
CN (1) | CN103679768B (zh) |
DE (1) | DE102012217163B4 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102012217163B4 (de) * | 2012-09-24 | 2022-06-02 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren |
EP3286736B1 (en) * | 2015-02-03 | 2018-09-05 | Koninklijke Philips N.V. | Image reconstruction system, method, and computer program |
CN107341836B (zh) * | 2016-05-03 | 2021-03-09 | 上海联影医疗科技股份有限公司 | 一种ct螺旋扫描图像重建方法及装置 |
US10198812B2 (en) | 2016-09-30 | 2019-02-05 | Toshiba Medical Systems Corporation | Data fidelity weight design for iterative reconstruction |
CN109300167B (zh) * | 2017-07-25 | 2023-10-10 | 清华大学 | 重建ct图像的方法和设备以及存储介质 |
US10463332B2 (en) * | 2018-02-02 | 2019-11-05 | Siemens Healthcare Gmbh | Spectral topogram-based determination of imaging parameters |
CN111757703B (zh) * | 2018-09-26 | 2023-04-07 | 上海联影医疗科技股份有限公司 | 用于图像生成的系统和方法 |
CN109949411B (zh) * | 2019-03-22 | 2022-12-27 | 电子科技大学 | 一种基于三维加权滤波反投影和统计迭代的图像重建方法 |
WO2025007524A1 (zh) * | 2023-07-06 | 2025-01-09 | 有方(合肥)医疗科技有限公司 | Ct成像方法、装置及cbct设备 |
CN116543071B (zh) * | 2023-07-06 | 2023-09-19 | 有方(合肥)医疗科技有限公司 | 大视野ct成像方法及装置 |
CN116543088B (zh) * | 2023-07-07 | 2023-09-19 | 有方(合肥)医疗科技有限公司 | Cbct图像重建方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1489977A (zh) * | 2002-09-23 | 2004-04-21 | 在计算机断层造影中设置图像的方法及实施该方法的设备 | |
CN1917811A (zh) * | 2004-02-16 | 2007-02-21 | 株式会社日立医药 | 断层摄影图象的再构成方法及断层摄影装置 |
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6055291A (en) * | 1998-08-25 | 2000-04-25 | General Electric Company | Systems, methods and apparatus for reconstructing images |
US6459754B1 (en) * | 1999-10-27 | 2002-10-01 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for cone beam multislice CT correction |
US6678346B2 (en) * | 2001-10-11 | 2004-01-13 | Ge Medical Systems Global Technology Company Llc | Cone-beam CT scanner with image reconstruction using multiple sub-images |
US6522714B1 (en) * | 2001-11-09 | 2003-02-18 | Ge Medical Systems Global Technologies Company Llc | Row-wise full helical view weighting method and apparatus for CT scanners |
US7215734B2 (en) | 2004-06-30 | 2007-05-08 | General Electric Company | Method and system for three-dimensional reconstruction of images |
US7672423B2 (en) * | 2005-10-20 | 2010-03-02 | Koninklijke Philips Electronics N.V. | Short scan cardiac CT on a quasi axial trajectory |
CN100565336C (zh) * | 2005-11-21 | 2009-12-02 | 清华大学 | 成像系统 |
DE102007021023A1 (de) * | 2007-05-04 | 2008-11-13 | Siemens Ag | Verfahren zur Bilderstellung für die Spiral-CT mit veränderlichem Pitch und CT-Gerät zur Durchführung des Verfahrens |
US7792238B2 (en) * | 2008-02-18 | 2010-09-07 | General Electric Company | Method and system for reconstructing cone-beam projection data with reduced artifacts |
WO2011002442A1 (en) * | 2009-06-30 | 2011-01-06 | Analogic Corporation | Efficient quasi-exact 3d image reconstruction algorithm for ct scanners |
JP6214128B2 (ja) * | 2010-11-22 | 2017-10-18 | キヤノン株式会社 | 画像処理装置、画像処理方法、及び記憶媒体 |
US8571291B2 (en) * | 2011-10-19 | 2013-10-29 | Kabushiki Kaisha Toshiba | Combination weight applied to iterative reconstruction in image reconstruction |
DE102012217163B4 (de) * | 2012-09-24 | 2022-06-02 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren |
-
2012
- 2012-09-24 DE DE102012217163.0A patent/DE102012217163B4/de not_active Expired - Fee Related
-
2013
- 2013-08-22 US US13/973,031 patent/US9299140B2/en not_active Expired - Fee Related
- 2013-09-23 CN CN201310435145.XA patent/CN103679768B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1489977A (zh) * | 2002-09-23 | 2004-04-21 | 在计算机断层造影中设置图像的方法及实施该方法的设备 | |
CN1917811A (zh) * | 2004-02-16 | 2007-02-21 | 株式会社日立医药 | 断层摄影图象的再构成方法及断层摄影装置 |
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
DE102012217163A1 (de) | 2014-06-12 |
US20140086466A1 (en) | 2014-03-27 |
CN103679768A (zh) | 2014-03-26 |
DE102012217163B4 (de) | 2022-06-02 |
US9299140B2 (en) | 2016-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103679768B (zh) | 以加权反投影重建ct图像数据的方法、计算单元和系统 | |
US5825842A (en) | X-ray computed tomographic imaging device and x-ray computed tomographic method | |
US8724889B2 (en) | Method and apparatus for CT image reconstruction | |
US8270561B2 (en) | Motion weighting in computed tomography (CT) with cone angle | |
US7403587B2 (en) | Computer tomography method using a cone-shaped bundle of rays | |
US6130930A (en) | Exact region of interest cone beam imaging without circle scans | |
US8116426B2 (en) | Computed tomography device and method using circular-pixel position-adaptive interpolation | |
US7245755B1 (en) | Algorithm for image reconstruction and image noise analysis in computed tomography | |
EP1308898A2 (en) | System and method for image reconstruction | |
US9805448B2 (en) | Method and apparatus for computed tomography using asymetric filter for volume half reconstruction | |
EP1800264B1 (en) | Image reconstruction with voxel dependent interpolation | |
US6084937A (en) | Adaptive mask boundary correction in a cone beam imaging system | |
EP1828985A1 (en) | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data | |
JP2003144429A (ja) | コーンビームctの再構成において欠落データを近似するための方法及びシステム、x線ct装置 | |
US6324245B1 (en) | Method and apparatus for identifying and correcting image inaccuracies caused by simplified processing of masked cone beam projection data | |
US6490335B1 (en) | Helical segment image reconstruction | |
US7430270B2 (en) | Computerized tomography method with helical relative movement and conical beam | |
JP3277167B2 (ja) | X線コンピュータ断層撮影装置 | |
JP3490699B2 (ja) | X線コンピュータ断層撮影装置 | |
EP1769460B1 (en) | A computed tomography method for the reconstruction of object images from real and fictitious measured values | |
JP2004089720A (ja) | X線コンピュータ断層撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220128 Address after: Erlangen Patentee after: Siemens Healthineers AG Address before: Munich, Germany Patentee before: SIEMENS AG |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170412 |