CN1223309C - 图像重建方法和x-射线计算机断层装置 - Google Patents
图像重建方法和x-射线计算机断层装置 Download PDFInfo
- Publication number
- CN1223309C CN1223309C CNB021426570A CN02142657A CN1223309C CN 1223309 C CN1223309 C CN 1223309C CN B021426570 A CNB021426570 A CN B021426570A CN 02142657 A CN02142657 A CN 02142657A CN 1223309 C CN1223309 C CN 1223309C
- Authority
- CN
- China
- Prior art keywords
- detector line
- pixel
- data
- detector
- ray
- 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
- 238000000034 method Methods 0.000 title claims description 53
- 238000003325 tomography Methods 0.000 title description 2
- 230000000007 visual effect Effects 0.000 claims description 31
- 238000010276 construction Methods 0.000 description 14
- 101100141323 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) RNR2 gene Proteins 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
Images
Classifications
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Algebra (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
为减小在应用多行检测器时由相对于旋转平面倾斜的扫描平面所引起的伪影,基于在z-轴方向上从扫描中心平面到图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)选择穿过像素点g(x,y)的X-射线所进入的第j检测器行,以及使用第j检测器行的数据确定该像素点g的像素值,这里在一定视角上的扫描中心平面定义为x-y平面,而垂直于所说的x-y平面的方向(即检测器行的方向)定义为z-轴方向。
Description
发明领域
本发明涉及一种图像重建方法和X-射线CT(计算机断层)装置,更具体地说,涉及这样的图像重建方法和X-射线CT装置:在应用多行检测器时它能够提供由相对于旋转平面倾斜的扫描平面所引起的伪影减少的图像。
背景技术
附图6所示为说明在日本专利申请公开No.H8-187241中描述的图像重建方法的解释性附图。
X-射线管11和X-射线检测器53绕旋转中心IC旋转以在不同的视角上采集数据。该扫描平面定义为x-y平面。
基本以如下的步骤产生图像:
(0)对在图像重建平面中的所有的像素点g(x,y)的像素值初始化为0;
(1)确定检测以视角β穿过在图像重建平面P中的像素点g(x,y)的X-射线的通道i;
(2)将通道i的数据加入到像素点g(x,y)的像素值中;
(3)在所要求的角度范围(例如360°)内对每个视角β重复(1)和(2)以获得像素点g(x,y)的像素值;以及
(4)对在图像重建平面P中的每个像素点g(x,y)重复(0)-(3)。
一旦已知通过像素点g(x,y)的X-射线与X-射线束的中心轴Bc所形成的角度γ就可以唯一地确定通道i。
通过下式可以计算角度γ:
这里,s=x·cosβ-ysinβ,和
t=x·sinβ-ycosβ
这里,x,y,β,s和t在附图6中的箭头所示的方向上都为正;从X-射线管11到旋转中心IC的距离由D表示;s表示从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离;以及t表示从像素点g朝X-射线束的中心轴Bc所引的垂线的长度。
附图7所示为在使用具有一个以上的检测器行的多行检测器13时穿过像素点g(x,y)的X-射线进入哪个检测器行的解释性附图。垂直于x-y平面的方向(即检测器行的方向)称为z-轴方向。应该注意的是附图7是从与X-射线束的中心轴Bc正交的方向看的视图。
多行检测器13在z-轴方向上、即在行方向上、越宽,扫描平面相对于旋转平面的倾斜越大。
通常,忽略这种倾斜,通过选择对应于图像重建平面P的z-轴位置的第m检测器行并使用第m检测器行的数据产生图像重建平面P的图像。
然而,在附图7中所示的穿过像素点g(x,y)的X-射线进入第j检测器行中,而没有进入第m检测器行中。因此,存在的问题是造成了伪影。
即存在这样的问题:在使用多行检测器时出现了由相对于旋转平面倾斜的扫描平面引起了伪影(锥角伪影)。
发明概述
因此本发明的一个目的是提供一种这样的图像重建方法和X-射线CT装置:在应用多行检测器时它能够提供由相对于旋转平面倾斜的扫描平面所引起的伪影减少的图像。
根据它的第一方面,本发明提供一种图像重建方法,在绕旋转中心旋转X-射线管和具有第一至第J(J≥2)检测器行的多行检测器中的至少一个的同时,该图像重建方法基于在不同的视角上所采集的数据重建指定的图像重建平面的图像,其特征在于包括:基于在z-轴方向上从扫描中心平面到图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)选择穿过像素点g(x,y)的X-射线进入的检测器行;以及使用所说的检测器行的数据确定所说的像素点g的像素值,这里在一定视角上的扫描中心平面定义为x-y平面,而垂直于所说的x-y平面的方向(即检测器行的方向)定义为z-轴方向。
在这种结构中,扫描中心平面定义为包含X-射线束的中心轴Bc并在通道方向上延伸的中心扫描平面。
根据第一方面的图像重建方法,基于在z-轴方向上从扫描中心平面到图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)确定穿过像素点g(x,y)的X-射线实际进入的检测器行;以及使用该检测器行的数据确定所说的像素点g的像素值。因此,在使用多行检测器时可减小相对于旋转平面倾斜的扫描平面所引起的伪影。
根据第二方面,本发明提供一种具有前述构造的图像重建方法,其特征在于选择通过如下方式确定的第j检测器行:
或者
这里
和
H=fdd·tanα,
这里从X-射线管到旋转中心的距离表示为D,从X-射线管到多行检测器的距离表示为fdd,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离表示为s,在z-轴方向上的一个检测器行的长度表示为Δd,以及rup{}表示为上舍入取整函数,rdwn{}表示下舍入取整函数。
根据第二方面的图像重建方法,可适当地确定穿过像素点g(x,y)的X-射线实际进入的第j检测器行。
此外,在J为奇数时以及在J为偶数并且H为正时使用
而在J为偶数并且H为负时使用
根据本发明的第三方面,本发明提供一种具有前述构造的图像重建方法,其特征在于选择第j检测器行和最近的检测器行,以及使用通过对所选择的检测器行的数据执行插值运算所获得的数据。
穿过像素点g(x,y)的X-射线具有一定范围,并且它不仅进入第j检测器行而且还进入最近的检测器行。
因此,根据本发明的第三方面的图像重建方法,使用通过组合经插值运算的检测器行的数据所获得的数据来确定像素点g的像素值。这进一步减小了伪影。
根据本发明的第四方面,本发明提供一种具有前述构造的图像重建方法,其特征在于选择第j检测器行和第(j-1)或第(j+1)检测器行,以及使用通过对所选择的检测器行的数据执行线性插值运算所获得的数据。
根据第四方面的图像重建方法,由于应用线性插值运算,因此插值运算处理可能是最简单的。
根据本发明的第五方面,本发明提供一种具有前述构造的图像重建方法,其特征在于在j>1时使j=J,而在j<1时使j=1。
根据第五方面的图像重建方法,如果穿过像素点g(x,y)的X-射线在行方向上行进到多行检测器第一或第二端之外,则选择在第一或第二端上的检测器行。因此,可以使用具有最小误差的可重建的数据。
根据本发明的第六方面,本发明提供一种具有前述构造的图像重建方法,其特征在于根据在螺旋扫描中的视角改变在z-轴方向上从扫描中心平面到图像重建平面的距离Δz。
在螺旋扫描中,床相对于台架移动,因此,扫描中心平面的位置根据视角在z-轴方向上移动。
因此,根据第六方面的图像重建方法,对于每个视角距离Δz逐渐地改变等于床移动的距离的量。因此,即使在螺旋扫描中也可以精确地确定穿过像素点g(x,y)的X-射线所进入的检测器行。
根据本发明的第七方面,本发明提供一种具有前述构造的图像重建方法,其特征在于对于每个视角将与{(D+S)2+t2}成反比的权重应用到所选择的检测器行的数据中,以及将加权的数据相加以确定像素点g的像素值,这里从X-射线管到旋转中心的距离表示为D,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离表示为s,以及从像素点g朝X-射线束的中心轴Bc所引的垂线的长度表示为t。
由X-射线管发出的X-射线穿过靠近X-射线管的像素点,此后经过远离X-射线管的像素点,并到达X-射线检测器。在这种情况下,靠近X-射线管的像素点比远离X-射线管更大地受到影响。
因此,根据本发明第七方面的图像重建方法,数据与从X-射线管到像素点的距离的平方成反比的权重相加。因此,产生了更加自然的图像。
根据本发明的第八方面,本发明提供一种具有前述构造的图像重建方法,其特征在于基于在图像重建平面中的像素点g的位置(x,y)和视角选择从其中获得用于确定所说的像素点的像素值的数据的通道。
根据第八方面的图像重建方法,可以选择穿过像素点g(x,y)的X-射线实际进入的通道。
根据本发明的第九方面,本发明提供一种X-射线CT装置,该X-射线CT装置包括X-射线管、具有第一至第J(J≥2)检测器行的多行检测器和数据采集装置,该数据采集装置在绕旋转中心旋转所说的X-射线管和所说的多行检测器之中的至少一个的同时在不同的视角上采集数据,其特征在于进一步包括:基于在z-轴方向上从扫描中心平面到指定的图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)选择检测器行的检测器行选择装置,这里在一定视角上的扫描中心平面定义为x-y平面,而垂直于所说的x-y平面的方向即检测器行的方向定义为z-轴方向;以及通过将所选择的检测器行的每个视角上的数据相加以确定所说的像素点g的像素值来重建图像重建平面的图像的图像产生装置。
根据第九方面的X-射线CT装置,适合于实施如在第一方面所描述的图像重建方法。
根据本发明的第十方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于所说的检测器行选择装置选择通过如下方式确定的第j检测器行:
或者
这里
和
H=fdd·tanα,
这里从X-射线管到旋转中心的距离表示为D,从X-射线管到多行检测器的距离表示为fdd,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离表示为s,在z-轴方向上的一个检测器行的长度表示为Δd,以及rup{}表示为上舍入取整函数,rdwn{}表示下舍入取整函数。
根据第十方面的X-射线CT装置,适合于实施如在第二方面所描述的图像重建方法。
根据本发明的第十一方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于所说的检测器行选择装置选择第j检测器行和最近的检测器行,以及所说的图像产生装置使用通过对所选择的检测器行的数据执行插值运算所获得的数据。
根据第十一方面的X-射线CT装置,适合于实施如在第三方面所描述的图像重建方法。
根据本发明的第十二方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于所说的检测器行选择装置选择第j检测器行和第(j-1)或第(j+1)检测器行,以及所说的图像产生装置使用通过对所选择的检测器行的数据执行线性插值运算所获得的数据。
根据第十二方面的X-射线CT装置,适合于实施如在第四方面所描述的图像重建方法。
根据本发明的第十三方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于在j>1时所说的检测器行选择装置使j=J,而在j<1时使j=1。
根据第十三方面的X-射线CT装置,适合于实施如在第五方面所描述的图像重建方法。
根据本发明的第十四方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于所说的检测器行选择装置根据在螺旋扫描中的视角改变在z-轴方向上从扫描中心平面到图像重建平面的距离Δz。
根据第十四方面的X-射线CT装置,适合于实施如在第六方面所描述的图像重建方法。
根据本发明的第十五方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于所说的图像产生装置对于每个视角将与{(D+S)2+t2}成反比的权重应用到所选择的检测器行的数据中,以及将经加权的数据相加以确定像素点g的像素值,这里从X-射线管到旋转中心的距离表示为D,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离表示为s,以及从像素点g朝X-射线束的中心轴Bc所引的垂线的长度表示为t。
根据第十五方面的X-射线CT装置,适合于实施如在第七方面所描述的图像重建方法。
根据本发明的第十六方面,本发明提供一种具有前述构造的X-射线CT装置,其特征在于包括通道选择装置以基于在图像重建平面中的像素点g的位置(x,y)和视角选择从其中获得用于确定所说的像素点的像素值的数据的通道。
根据第十六方面的X-射线CT装置,适合于实施如在第八方面所描述的图像重建方法。
根据本发明的图像重建方法和X-射线CT装置,在使用多行检测器时可以减小由相对于旋转平面倾斜的扫描平面所引起的伪影。
通过在附图中所示的本发明的优选实施例的描述,将会清楚本发明的进一步目的和优点。
附图概述
附图1所示为根据第一实施例的X-射线CT装置的方块图;
附图2所示为根据第一实施例的图像重建处理的解释性附图(从z-轴方向看)。
附图3所示为根据第一实施例的图像重建处理的解释性附图(从垂直于z-轴和X-射线束中心轴的方向看)。
附图4所示为根据第二实施例的图像重建处理的解释性附图(从垂直于z-轴和X-射线束中心轴的方向看)。
附图5所示为根据第四实施例的图像重建处理的解释性附图(从垂直于z-轴和X-射线束中心轴的方向看)。
附图6所示为常规的图像重建处理的解释性附图(从z-轴方向看)。
附图7所示为解释本发明的目的的附图(从垂直于z-轴和X-射线束中心轴的方向看)。
本发明的详细描述
下文根据附图更详细地描述本发明的几个实施例。
第一实施例
附图1所示为根据本发明的第一实施例的X-射线CT装置100的方块图。
该X-射线CT装置100包括操作台1、床装置8和扫描台架9。
操作台1包括用于接收由操作员输入的信息和指令的输入装置2、执行扫描处理和图像重建处理的中央处理装置3、与成像床8和扫描台架9进行控制信号等信号通信的控制接口4、采集由扫描台架9所获取的数据的数据采集缓冲器5、显示由该数据所重建的图像的CRT6和存储程序、数据和图像的存储装置7。
床装置8包括放置对象的托架8c和在z-和y-轴方向上移动托架8c的运动控制器8a。
y-轴代表垂直方向,而z-轴代表托架8c的纵向方向。此外,垂直于y-轴和z-轴的轴定义为x-轴。对象的身体轴沿着z-轴方向延伸。
扫描台架9包括X-射线控制器10、X-射线管11、准直器12、具有一个以上检测器行的多行检测器13、数据采集部分14和绕旋转中心IC旋转X-射线管11和多行检测器13的旋转控制器15。
附图2和3所示为根据本发明的第一实施例的图像重建方法的解释性附图。
附图2为z-轴方向的视图。X-射线管11和多行检测器13绕旋转中心IC旋转以从不同的视角采集数据。
附图3所示为垂直于z-轴和X-射线束的中心轴Bc的方向的视图。多行检测器13在z-轴方向上具有第一至第J(J≥2)检测器行。
扫描平面相对于旋转平面倾斜。包含X-射线束的中心轴Bc并在该通道的方向上延伸的中心扫描平面称为扫描中心平面。扫描中心平面与旋转平面一致,并且处于x-y平面中。
图像产生的基本过程如下:
(0)对在图像重建平面中的所有的像素点g(x,y)的像素值初始化为0;
(1)确定检测以视角β穿过在图像重建平面P中的像素点g(x,y)的X-射线的通道i,如附图2所示;
(2)确定用于检测穿过在图像重建平面P中的像素点g(x,y)的X-射线的第j检测器行,如附图3所示;
(3)计算在第j检测器行中的通道i的数据的权重W;
(4)将在第j检测器行中的通道i的数据乘以权重W并将经加权的数据加入到像素点g(x,y)的像素值中;
(5)在所要求的角度范围(例如360°)内对每个视角β重复(1)-(4)以获得像素点g(x,y)的像素值;以及
(6)对在图像重建平面P中的每个像素点g(x,y)重复(0)-(5)。
在(1)中,一旦已知通过像素点g(x,y)的X-射线与X-射线束的中心轴Bc所形成的角度γ就可以唯一地确定通道i,这可以从附图2中看出。
通过下式可以计算角度γ:
这里,s=x·cosβ-y·sinβ,和
t′=x·sinβ-y·cosβ
这里,x,y,β,s和t′在附图6中的箭头所示的方向上都是正的;从X-射线管11到旋转中心IC的距离由D表示;s表示从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离;以及t′表示从像素点g朝X-射线束的中心轴Bc所引的垂线在扫描中心平面上的投影的长度。
在(2)中,当已知通过像素点g(x,y)的X-射线与X-射线束的中心轴Bc所形成的角度α就可以唯一地确定第j检测器行,这可以从附图3中看出。
通过下式可以计算角度α:
这里,s和Δz在附图3中的箭头所示的方向上都是正的;以及Δz表示在z-轴方向上从扫描中心平面到图像重建平面的距离。
现在可以通过下式确定j:
或者
这里 H=fdd·tanα,
这里,从X-射线管11到多行检测器13的距离表示为fdd,在z-轴方向上的一个检测器行的长度表示为Δd,以及rup{}表示为上舍入取整函数,rdwn{}表示下舍入取整函数。
在J为奇数时以及在J为偶数并且H为正时使用
而在J为偶数并且H为负时使用
在螺旋扫描中,虽然图像重建平面P的位置并不改变,但是扫描中心平面的位置根据视角移动。因此,Δz应该随着视角变化。
根据本发明的第一实施例的X-射线CT装置100,确定穿过像素点g的X-射线实际进入的通道i和第j检测器行,并使用第j检测器行和通道i的数据以确定像素点g的像素值。因此,在使用多行检测器13时可以减小由相对于旋转平面倾斜的扫描平面所引起的伪影。
第二实施例
在第二实施例中,不仅使用第j检测器行的数据而且还使用最近的检测器行的数据,并且从通过对这些数据执行线性插值运算所获得的数据中重建图像。
如在附图4中所示,首先,与第一实施例一样,选择第j检测器。
如果通过
选择第j检测器,并且
则选择第j-1检测器行。然后,用d(j,i)表示第j检测器和通道i的数据和用d(j-1,i)表示第(j-1)检测器和通道i的数据,将通过下述线性插值运算获得数据d乘以权重W:
d=a·d(j-1,i)+b·d(j,i),
这里
和
b=1-a,
并将经加权的数据加入到像素点g的像素值。
此外,如果通过
选择第j检测器,并且 则选择第j+1检测器行。然后,用d(j,i)表示第j检测器和通道i的数据和用d(j+1,i)表示第(j+1)检测器和通道i的数据,将通过下述线性插值运算获得数据d乘以权重W:
d=a·d(j+1,i)+b·d(j,i),
这里
和
b=1-a,
并将经加权的数据加入到像素点g的像素值。
在另一方面,如果通过
选择第j检测器,并且 则选择第j-1检测器行。然后,通过d(j,i)表示第j检测器和通道i的数据和通过d(j-1,i)表示第(j-1)检测器和通道i的数据,将通过下述线性插值运算获得数据d乘以权重W:
d=a·d(j-1,i)+b·d(j,i),
这里
和
b=1-a,
并将经加权的数据加入到像素点g的像素值。
此外,如果通过
选择第j检测器,并且 则选择第j+1检测器行。然后,通过d(j,i)表示第j检测器和通道i的数据和通过d(j+1,i)表示第(j+1)检测器和通道i的数据,将通过下述线性插值运算获得数据d乘以权重W:
d=a·d(j+1,i)+b·d(j,i),
这里
和
b=1-a,
并将经加权的数据加入到像素点g的像素值。
根据第二实施例的X-射线CT装置,可以比第一实施例更多地减小伪影。
可取的是,在通道方向即对通道i的数据和通道(i-1)或(i+1)的数据也执行线性插值。
第三实施例
不使用线性插值运算,而是在检测器行方向和通道方向上可以执行二维汉宁插值(Hanning interpolation)。可替换的是,也可以执行三次插值或拉格朗日插值。虽然在这种插值中运算处理比较复杂,但是与线性插值运算相比可以获得良好的图像质量。
第四实施例
如附图5所示,如果j>J,则使j为j=J以选择第J检测器行。
如果j<1,则使j为j=1以选择第一检测器行。
根据第四方面的图像重建方法,即使穿过像素点g的X-射线行进到多行检测器13的第一或第二端之外,则仍然可以使用可能的最适合的数据。
在不脱离本发明的精神范围的前提下可以构造出许多不同的实施例。应该理解的是本发明并不限于在说明书中所描述的特定的实施例,而是以所附加的权利要求来限定。
Claims (16)
1.一种图像重建方法,在绕旋转中心旋转X-射线管和具有第一至第J检测器行的多行检测器中的至少一个的同时,该图像重建方法基于在不同的视角上所采集的数据重建指定的图像重建平面的图像,其中J≥2,该图像重建方法包括如下的步骤:
基于在z-轴方向上从扫描中心平面到图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)选择穿过像素点g(x,y)的X-射线所进入的检测器行;以及
使用所说的检测器行的数据确定所说的像素点g的像素值,这里在一定视角上的扫描中心平面定义为x-y平面,而垂直于所说的x-y平面的方向即检测器行的方向定义为z-轴方向。
2.根据权利要求1所述的图像重建方法,包括选择通过如下方式确定的第j检测器行的步骤:
或者
这里
,和
H=fdd·tanα,
这里,从X-射线管到旋转中心的距离表示为D,从X-射线管到多行检测器的距离表示为fdd,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转中心的距离表示为s,在z-轴方向上的一个检测器行的长度表示为Δd,以及rup{}表示为上舍入取整函数,rdwn{}表示下舍入取整函数。
3.根据权利要求2所述的图像重建方法,包括如下的步骤:选择与第j检测器行相邻的检测器行,以及使用通过对所选择的检测器行的数据执行插值运算所获得的数据。
4.根据权利要求2所述的图像重建方法,包括如下的步骤:选择第j-1或第j+1检测器行,以及使用通过对所选择的检测器行的数据执行线性插值运算所获得的数据,其中所选择的检测器行是第J检测器行和第j-1或第j+1检测器行。
5.根据权利要求2所述的图像重建方法,包括如下的步骤:在j>1时使j=J,而在j<1时使j=1。
6.根据权利要求1所述的图像重建方法,包括如下的步骤:根据在螺旋扫描中的视角改变在z-轴方向上从扫描中心平面到图像重建平面的距离Δz。
7.根据权利要求1所述的图像重建方法,包括如下的步骤:对于每个视角将与{(D+s)2+t2}成反比的权重应用到所选择的检测器行的数据中,以及将经加权的数据相加以确定像素点g的像素值,这里从X-射线管到旋转中心的距离表示为D,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转中心的距离表示为s;以及从像素点g朝X-射线束的中心轴Bc所引的垂线长度表示为t。
8.根据权利要求1所述的图像重建方法,包括如下的步骤:基于在图像重建平面中的像素点g的位置(x,y)和视角选择从其中获得用于确定所说的像素点g的像素值的数据的通道。
9.一种X-射线CT装置,该X-射线CT装置包括:
X-射线管;
具有第一至第J检测器行的多行检测器,其中J≥2;
数据采集装置,该数据采集装置在绕旋转中心旋转所说的X-射线管和所说的多行检测器之中的至少一个的同时在不同的视角上采集数据;
检测器行选择装置,该检测器行选择装置基于在z-轴方向上从扫描中心平面到所指定的图像重建平面的距离Δz和在该图像重建平面中的像素点g的位置(x,y)选择穿过像素点g(x,y)的X-射线所进入的检测器行,这里在一定视角上的扫描中心平面定义为x-y平面,而垂直于所说的x-y平面的方向即检测器行的方向定义为z-轴方向;以及
图像产生装置,该图像产生装置用于通过将所选择的检测器行在每个视角上的数据相加以确定所说的像素点g的像素值来重建图像重建平面的图像。
10.根据权利要求9所述的X-射线CT装置,其中所说的检测器行选择装置选择通过如下方式确定的第j检测器行:
或者
这里
,和
H=fdd·tanα,
这里从X-射线管到旋转中心的距离表示为D,从X-射线管到多行检测器的距离表示为fdd,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转的中心的距离表示为s,在z-轴方向上的一个检测器行的长度表示为Δd,以及rup{}表示为上舍入取整函数,rdwn{}表示下舍入取整函数。
11.根据权利要求10所述的X-射线CT装置,其中在j>1时所说的检测器行选择装置使j=J,而在j<1时使j=1。
12.根据权利要求10所述的X-射线CT装置,其中所说的检测器行选择装置选择与第j检测器行相邻的检测器行,以及所说的图像产生装置使用通过对所选择的检测器行的数据执行插值运算所获得的数据。
13.根据权利要求10所述的X-射线CT装置,其中所说的检测器行选择装置选择第j-1或第j+1检测器行,以及所说的图像产生装置使用通过对所选择的检测器行的数据执行线性插值运算所获得的数据。
14.根据权利要求9所述的X-射线CT装置,其中所说的检测器行选择装置根据在螺旋扫描中的视角改变在z-轴方向上从扫描中心平面到图像重建平面的距离Δz。
15.根据权利要求9所述的X-射线CT装置,其中所说的图像产生装置对于每个视角将与{(D+s)2+t2}成反比的权重应用到所选择的检测器行的数据中,以及将经加权的数据相加以确定像素点g的像素值,这里从X-射线管到旋转中心的距离表示为D,从像素点g朝X-射线束的中心轴Bc所引的垂线的垂足到旋转中心的距离表示为s;以及从像素点g朝X-射线束的中心轴Bc所引的垂线的长度表示为t。
16.根据权利要求9所述的X-射线CT装置,进一步包括通道选择装置,用于基于在图像重建平面中的像素点g的位置(x,y)和视角选择从其中获得用于确定所说的像素点g的像素值的数据的通道。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP281732/2001 | 2001-09-17 | ||
JP2001281732 | 2001-09-17 | ||
JP281732/01 | 2001-09-17 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1404797A CN1404797A (zh) | 2003-03-26 |
CN1223309C true CN1223309C (zh) | 2005-10-19 |
Family
ID=19105503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB021426570A Expired - Fee Related CN1223309C (zh) | 2001-09-17 | 2002-09-17 | 图像重建方法和x-射线计算机断层装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US7099503B2 (zh) |
CN (1) | CN1223309C (zh) |
DE (1) | DE10242910A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL156569A (en) * | 2003-06-22 | 2009-11-18 | Ultraspect Ltd | Improved single-plane plane emission simulation imaging method |
JP2005177203A (ja) * | 2003-12-22 | 2005-07-07 | Ge Medical Systems Global Technology Co Llc | 複数位置のct画像生成方法およびx線ct装置 |
JP2007044207A (ja) * | 2005-08-09 | 2007-02-22 | Ge Medical Systems Global Technology Co Llc | 放射線ct撮影方法およびx線ct装置 |
JP2009514636A (ja) * | 2005-11-09 | 2009-04-09 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | X線検出器において3dゴーストアーチファクトを低減させる方法 |
CN100435733C (zh) * | 2005-12-31 | 2008-11-26 | 清华大学 | X-ct扫描*** |
DE102011003240B4 (de) * | 2011-01-27 | 2017-02-02 | Siemens Healthcare Gmbh | Verfahren und Computersystem zur Reduktion von Artefakten in rekonstruierten CT-Bilddatensätzen |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416815A (en) | 1993-07-02 | 1995-05-16 | General Electric Company | Adaptive filter for reducing streaking artifacts in x-ray tomographic images |
US5515409A (en) | 1994-12-22 | 1996-05-07 | General Electric Company | Helical interpolative algorithm for image reconstruction in a CT system |
DE19528261C1 (de) | 1995-08-01 | 1996-10-31 | Siemens Ag | Verfahren zur Bildrekonstruktion bei einem Computertomographen |
US5727041A (en) | 1996-11-13 | 1998-03-10 | General Electric Company | Methods and apparatus for reducing partial volume image artifacts |
US5812628A (en) | 1996-12-12 | 1998-09-22 | General Electric Company | Methods and apparatus for detecting partial volume image artifacts |
US6115487A (en) | 1998-01-08 | 2000-09-05 | General Electric Company | Correction algorithm for bone-induced spectral artifacts in computed tomograph imaging |
US6035012A (en) | 1998-05-14 | 2000-03-07 | Gen Electric | Artifact correction for highly attenuating objects |
US6134292A (en) | 1998-07-13 | 2000-10-17 | General Electric Company | Methods and apparatus for reducing z-axis non-uniformity artifacts |
US6061419A (en) | 1998-08-25 | 2000-05-09 | General Electric Company | Methods and apparatus for noise compensation in an imaging system |
US6061420A (en) | 1998-08-25 | 2000-05-09 | General Electric Company | Methods and apparatus for graphical Rx in a multislice imaging system |
US6215841B1 (en) | 1998-09-29 | 2001-04-10 | General Electric Company | Methods and apparatus for 3D artifact reduction |
US6233308B1 (en) | 1999-03-19 | 2001-05-15 | General Electric Company | Methods and apparatus for artifact compensation with variable angular sampling |
US6421411B1 (en) | 2001-05-10 | 2002-07-16 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for helical image artifact reduction |
-
2002
- 2002-09-16 DE DE2002142910 patent/DE10242910A1/de not_active Ceased
- 2002-09-16 US US10/244,287 patent/US7099503B2/en not_active Expired - Fee Related
- 2002-09-17 CN CNB021426570A patent/CN1223309C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US20030063786A1 (en) | 2003-04-03 |
CN1404797A (zh) | 2003-03-26 |
US7099503B2 (en) | 2006-08-29 |
DE10242910A1 (de) | 2003-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1264480C (zh) | 三维反投影方法和一种x射线计算机层析成像装置 | |
CN1682658A (zh) | 放射层析成像装置和使用该装置的成像方法 | |
CN1775179A (zh) | X射线ct设备和x射线ct成像方法 | |
CN1247028C (zh) | 运动检测装置及搜索区域形状可变运动检测器 | |
CN1726870A (zh) | X射线ct图像处理方法和x射线ct设备 | |
CN1698540A (zh) | 图像处理方法、图像处理***以及x-射线ct*** | |
CN1159574C (zh) | Ct装置 | |
CN1573321A (zh) | X光照相设备 | |
CN1403058A (zh) | 图像产生方法和x-射线计算机断层成像装置 | |
CN1254216C (zh) | 放射线断层成像装置和方法 | |
CN1918605A (zh) | 运动假影补偿 | |
CN1474358A (zh) | 三维空间滤波设备和方法 | |
CN101079952A (zh) | 图像处理方法和图像处理设备 | |
CN1637380A (zh) | 宽度测量方法和表面特性测量机 | |
CN1518956A (zh) | 三维反向投影方法和x射线ct装置 | |
CN101032409A (zh) | 图像显示设备和x射线ct设备 | |
CN1678036A (zh) | 射线照相装置和射线检测信号处理方法 | |
CN1319495C (zh) | 一种锥束x射线ct***的定标模板 | |
CN1879562A (zh) | X射线ct图象重建方法和x射线ct*** | |
CN1959740A (zh) | 图像处理方法及装置、图像处理用程序及其存储介质 | |
CN1286431C (zh) | 射线束固化的后处理方法和x射线ct设备 | |
CN1572253A (zh) | 锥面射束重建设备以及计算机断层摄影设备 | |
CN1539377A (zh) | 辐射计算机断层扫描装置及断层图像产生的方法 | |
CN1275033C (zh) | 操纵产品并处理产品的射线检验图像的装置和方法及用途 | |
CN1596828A (zh) | 放射计算机断层摄像装置及断层图像数据生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20051019 Termination date: 20091019 |