CN100528088C - 用圆-线快速扫描算法重构计算机断层分析图像的方法 - Google Patents
用圆-线快速扫描算法重构计算机断层分析图像的方法 Download PDFInfo
- Publication number
- CN100528088C CN100528088C CN200610009268.7A CN200610009268A CN100528088C CN 100528088 C CN100528088 C CN 100528088C CN 200610009268 A CN200610009268 A CN 200610009268A CN 100528088 C CN100528088 C CN 100528088C
- Authority
- CN
- China
- Prior art keywords
- ideal
- actual
- circle
- track
- incomplete circle
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 230000005855 radiation Effects 0.000 claims description 10
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 238000011017 operating method Methods 0.000 claims 1
- 238000001914 filtration Methods 0.000 description 11
- 238000002591 computed tomography Methods 0.000 description 8
- 239000013598 vector Substances 0.000 description 6
- 239000013256 coordination polymer Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 208000030984 MIRAGE syndrome Diseases 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- TVLSRXXIMLFWEO-UHFFFAOYSA-N prochloraz Chemical compound C1=CN=CN1C(=O)N(CCC)CCOC1=C(Cl)C=C(Cl)C=C1Cl TVLSRXXIMLFWEO-UHFFFAOYSA-N 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 206010019133 Hangover Diseases 0.000 description 1
- 208000032843 Hemorrhage Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000007408 cone-beam computed tomography Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000001361 intraarterial administration Methods 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010603 microCT Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036651 mood Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000002792 vascular Effects 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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; 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Pure & Applied Mathematics (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)
- Image Analysis (AREA)
Abstract
在根据从检查对象获取的数据重构CT图像的方法中,采用基于理想快速扫描的包括圆和线的轨道的重构算法。为了使该重构算法适于“真实世界”扫描轨道,用C臂CT设备获取数据,其中沿着实际快速扫描的包括圆和线的轨道移动焦点。对于实际轨道中的焦点的每一位置,以电子方式产生投影矩阵,并利用该投影矩阵而使具有理想轨道的重构算法适于实际轨道。
Description
技术领域
本发明针对一种重构CT图像的方法,并更具体地,针对一种利用用于快速扫描的结合圆和各种线的重构算法重构CT图像的方法。
背景技术
在临床干涉环境中已确立了基于在C臂系统中获取的锥面光束投影的高对比度对象的断层分析成像,例如如M.Grass,R.Koppe,E.Klotz,R.Proksa,M.Kuhn,H.Aerts,J.O.de Beck,和R.Kempkers,“Three-dimensional reconstructionof high-contrast objects using C-arm image intensifier projection data”Comp.Med.Imag.and Graphics 23,pp.311-321,1999所述。特别是,在神经系放射学中,复杂维管树的3D表现具有规划或验证治疗的高度临床价值。由于造影剂的侵入式动脉注射,该维管树具有比周围组织(例如骨头)高很多的对比度。由此,该过程对失真和图像假象相对不敏感。近期例如通过使用平板检测器对数据获取的改善将临床应用移动到低对比度对象成像。例如,在同一C臂系统上对中风的诊断和治疗是非常期望的目标。这将需要在治疗局部缺血之前排除脑实体的溢血。根据当前的临床规则,这通过本地的计算机断层分析(CT)扫描而完成。软组织成像需要精确的数据获取和处理,如J.Wiegert,M.Bertram,D.Schaefer,N.Conrads,N.Noordhoek,K.de Jong,T.Aach,和G.Rose,“Soft tissuecontrast resolution within the head of human cadaver by means of flat detector basedcone-beam CT”in Proc.SPIE 5368,pp.330-337,2004所述。严重的局限是由传统快速扫描圆形源轨道获取的投影数据的不完全性。由该不完全性导致的锥面假象作为拖尾和阴影假象而存在,并可叠加特别重要的低对比度细节。
在以下文献中可发现对满足Tuy完全性条件(见H.K.Tuy,“An inversionformula for cone-beam reconstruction”in SIAM J.Appl.Math.1983)的源轨道的很多调查:鞍状轨道(J.Pack,F.Noo,和H.Kudo,“Investigation of saddle trajectoriesfor cardiac CT imaging in cone-beam geometry”Phys.Med.Biol.49,pp.2317-2336,2004)、为C臂装置优化的非平面、非封闭的轨道的选择(H.Schomberg,“Complete source trajectories for C-arm systems and a method for coping withtruncated cone-beam projections,”in Proc.Meeting on Fully 3-D ImageReconstruction in Radiology and Nucl.Med.,2001)、为CT台架优化的圆形和弧形轨道(R.Ning,X.Tang,D.Conover,和R.Yu,“Flat panel detector-based conebeam computed tomography with a circle-plus-two-arcs data acquisition orbit:Preliminary phantom study,”Med.Phys.30,pp.1694-1705,2003)、圆形和线形轨道(G.L.Zeng和G.T.Gullberg,“A cone-beam tomography algorithm for orthogonalcircle-and-line orbit,”Phys.Med.Biol.37,pp.563-577,1992,以及R.Johnson,H.Hu,S.Haworth,P.Cho,C.Dawson,和J.Linehan,“Feldkamp and circle-and-linecone-beam reconstruction for 3D micro-CT of vascular networks,”Phys.Med.Biol.43,pp.929-940,1998以及H.Kudo和T.Saito,“Fast and stable cone-beamfiltered back-projection method for non-planar orbits,”Phys.Med.Biol.43,pp.747-760,1998)等等。
发明内容
本发明的目的是提供一种用于采用利用快速扫描的包括圆和线的轨道的算法来重构CT图像的方法,其可在现有C臂装置中容易地实现,而无需任何硬件修改。该线形扫描可被看作是传统快速扫描圆圈路径的扩充。该方案应该在理论上是精确的,具有有效的、偏移不变(shift-invariant)的滤波后方投影(FBP)结构,并解决长对象问题。该算法在处理各种圆和线的配置的情况下应该是灵活的。该重构方法仅需要扫描轨道的理论上的最小长度。
这些目的根据本发明被实现为一种用于根据用C臂设备从对象获取的数据重构该对象的CT图像的方法,该C臂设备具有安装在C臂上的x射线源和辐射检测器,所述x射线源具有从中以锥面光束发射x射线的焦点,所述操作通过以下步骤而完成:沿着聚焦轨道而围绕对象旋转所述x射线源的焦点,并用所述辐射检测器检测由所述对象衰减的辐射。根据本发明,操作所述C臂以沿着包括实际不完全圆圈和与所述实际不完全圆圈的一端连接的实际直线段的实际聚焦轨道而移动所述x射线源的所述焦点,并对于实际聚焦轨道中的每一焦点位置而检测由对象衰减的辐射。对于所述实际聚焦轨道中的所述焦点的每一位置,以电子方式计算投影矩阵,该矩阵对于该焦点位置而描述该辐射检测器上的对象的透视锥面光束投影。利用基于包括理想不完全圆圈和与所述理想不完全圆圈的一端连接的理想直线段的理想聚焦轨道的已知重构算法而重构该对象的图像,其中利用所述投影矩阵而使所述已知重构算法中的所述理想轨道适于所述C臂设备的所述实际聚焦轨道。
该创造性重构算法基于使用理想源轨道的重构算法,参见A.Katsevich,“Image reconstruction for the circle and line trajectory,”Phys.Med.Bial.49,pp.5059-5072,2004(根据其给出以下对已知反演算法和图1-4的讨论)以及A.Katsevich,“A general scheme for constructing inversion algorithms for cone beamCT,”International Journal of Mathematics and Mathematical Sciences 21,pp.1305-1321,2003。然而,C臂装置呈现了不得不考虑的某些机械不稳定性。幸运的是,与理想源路径的几何偏离几乎是可再现的,并通过几何校准处理来说明。非理想源轨道的投影几何结构在投影矩阵的框架中进行了方便的描述。用于描述投影几何结构的投影矩阵的一般用途在R.Hartley and A.Zisserman,Multiple View Geometry in Computer Vision,Cambridge University Press,2000中进行了描述。通过直接使用投影矩阵而精确执行后方投影步骤。滤波步骤需要更精细的适应策略。本创造性方法是使重构算法适于非理想采样图案的简单但鲁棒的方案,因为这些非理想采样图案出现在利用真实世界C臂装置的成像中。
附图说明
图1图示了用于解释该创造性方法的基本圆形和线形轨道。
图2图示了当源在线上时的检测器平面上的投影。
图3图示了当源在圆上时的检测器平面上的投影。
图4a和4b分别图示了可利用该创造性方法处理的轨道的例子。
图5图示了根据该创造性方法的理想圆形和线形轨道在真实世界(实际)焦点位置集合中的拟合。
具体实施方式
在以下讨论中,将讨论用于理想快速扫描(不完全的)圆形和线形轨道的公知反演算法。随后讨论根据该创造性方法使得算法适应偏离理想轨道的源路径。最后为论证该创造性重构算法的可利用性的实验的讨论。
用于理想快速扫描圆和线形轨道的公知反演算法
设想源轨道包括不完全的圆圈C和在C的一个端点与C连接的线段L(见图1)。首先,我们假设C充分接近完整的圆圈,并且L充分长。假设y0是它们的交点。假设检测器阵列DP(s)是扁平的,包括x3轴(C的轴),并且与连接源y(s)和x3轴的最短线段垂直。
图1图示了该圆形和线形轨道。
使用以下符号表示法。S2是IR3中的单位球面,并且
II(x,ξ):={z∈IR3:(z-x)·ξ=0}.
假设在源轨道附近,f是平滑的、紧密支持的、并总是等于零。
假设 和 分别是线和圆圈的参数化。假设圆圈的半径为R,而圆心为原点。假设U是开放集,使得
挑选重构点x∈U,并设想平面∏(x)通过x和L。∏(x)与C相交两点。一个点是y0,另一个被标为yC(x)。假设L1∏(x)是包括x并且将yC(x)连接到L的线段(见图1)。然后,YL(x)∈L表示L1∏(x)的另一端点。公知过程确定两个参数区间。第一个 对应于y0和YL(x)之间的L部分。第二个 对应于y0和YC(x)之间的C部分。由L1∏(x)限制范围的C∪L部分表示为Λ1∏(x)。可以容易地看出Λ1∏(x)在Tuy的意义上是完全的。
设想通过x的平面与Λ1∏(x)相交。忽略与轨道相切的平面,可以存在一个或三个交点(IP)。此外,可能最多有一个IP属于L。考虑到这一点,表1中的数据定义了最高为零测度集合的加权函数n。
n的作用是双重的。首先,不得不通过向Radon平面和源轨道之间的IP分配权重而处理锥面光束数据中的冗余。其次,n的正确选择得出Katsevich一般反演公式的框架中的有效偏移不变卷积后方投影算法。可如下描述表1所述的函数n。如果存在一个IP,则给定权重1。如果存在三个IP,则圆圈上的两个IP的每个都具有权重1,而线段上的IP具有权重-1。可以容易地看出,n被归一化为:对于几乎所有α∈S2,∑jn(sj,x,α)=1。这里,在所有交点上求和,y(sj)∈II(x,α)∩Λ1π(x)。
表1:加权函数n(s,x,α)的定义
情况 | n |
1个IP,s<sub>1</sub>∈I<sub>1</sub>(x) | n(s<sub>1</sub>,x,α)=1 |
1个IP,s<sub>1</sub>∈I<sub>2</sub>(x) | n(s<sub>1</sub>,x,α)=1 |
3个IP,s<sub>1</sub>∈I<sub>1</sub>(x)s<sub>2</sub>,s<sub>3</sub>∈I<sub>2</sub>(x) | n(s<sub>1</sub>,x,α)=-1n(s<sub>k</sub>,x,α)=1,k=2,3 |
表示
其中θ是与β(s,x)正交的平面中的极角。根据Katsevich描述的一般方案,必须在θ中找到Φ(s,x,θ)的跳变。通过研究两种情况:s∈I1(x)和s∈I2(x)下的跳变并利用该一般方案,可获得以下反演算法。挑选s∈I1(x)(即y(s)在线上)。找到通过x和y(s)的平面,其在某yt(s,x),s∈I2(x)处与C相切。假设u1(s,x)是与该平面正交的单位矢量:
现在挑选s∈I2(x)(即y(s)在圆上)并定义
通过构造,u2(s,x)是与包括x和y(s)并在y(s)处与C相切的平面正交的单位矢量。利用(3)和(4),我们获得用于 的以下重构公式:
其中
Θk(s,x,γ):=cosγβ(s,x)+sinγek(s,x),ek(s,x):=β(s,x)xuk(s,x). (6)
并且δk被定义为:
假设例如L以这样的方式被参数化:使得源随着s的增加而沿着L向下移动。然后,δ1(s,x)=1,s∈I1(x)。如果源随着s的增加而沿着L向上移动,则δ1(s,x)=1,s∈I1(x)。
图2图示了当源在线上时到检测器平面上的投影。
现在设想该算法的计算结构。挑选y(s)∈L。对于点x∈U,必须得到si∈I2(X)。这确定了检测器上的滤波线,其在与相切。这里,和分别是C和y(st)到检测器平面上的投影。容易看出在y(st)左边的投影在该线上的所有其他x∈U都将该线共享为它们的滤波线。因此,我们可首先沿着与相切的检测器上的线执行滤波(见图2的族L1),并然后执行后方投影。St值的范围 取决于感兴趣的区域(ROI)并如图2所示。可以容易地看出,该滤波是偏移不变的,并包括用1/sinγ卷积
图3图示了当源在圆上时到检测器平面上的投影。
如果y(s)∈C,则必须沿着与平行的检测器上的线执行滤波。得到的族用图3中的L2表示。从L2中挑选任一线。这示出了其投影属于该线并出现在的右边的所有x都将该线共享为它们的滤波线。和前面一样,可首先沿着这些线执行滤波(即用1/sinγ卷积),并随后进行后方投影。由此,得到的算法是基于卷积的FBP类型。
该算法的一些性质如下。根据L1π(x)的构造,随着x3→0,yL(x)→y0。在极限点x3=0,yL(x)=y0,使得(5)中L上的积分消失,而C上的积分变为超快速扫描扇形光束重构公式。
给定特定C和L,可确定可由该算法精确重构的f的支持的部分。这是由以下三个表面限制范围的空间:C平面、由L和不在L上的C的端点限定的平面、以及连接C的点和不在C上的L的端点的线的锥形面。该空间将由U(C,L)表示。然而,应该注意对象f可延伸到U(C,L)外面,只要其远离源轨道C∪L即可。
包括不完全圆圈和线段的轨道可用作用于构造其他轨道的构件块。例如,可设想这样的不完全圆圈C:线段连接到C的每一端点。这些线段可以在C的对侧(见图4)、或在C的同一侧(见图4b)。通过向每一圆圈+线段的子集应用(5)并然后将结果相加(如果必要的话),而从(5)获得这些轨道的反演算法。事实上,假设这些线段在C的对侧。然后,利用轨道C∪L重构z≥0的半空间(half-space)中的空间U(C,L),并利用C∪L′重构z≤0的空间U(C,L’)。在该情况下,不需要求和。如果L和L’在C的同一侧,则仅在z≥0的半空间进行该重构。在该情况下,利用C∪L和C∪L′对空间U(C,L)∩U(C,L’)中的每一三维像素(voxel)进行两次重构,从而使用求和。这并不意味着重构时间是两倍的长度。首先,扫描的线部分L和L’每个仅被使用一次。其次,仅圆圈部分C被使用两次。这并不导致计算时间的任何增加,因为两种情况下的滤波和后方投影是一样的。结果,简单的后滤波权重解决了对任何给定三维像素的多重贡献的问题。
现在考虑整体检测器需求。假设L和L’在C的对侧,并且重构空间是 由此,圆形扫描需要具有以下尺寸的矩形检测器:
这里d1和d2是检测器上水平和垂直轴。Katsevich已示出了线形扫描需要具有以下尺寸的检测器:
由此,与传统Feldkamp型圆形重构相比,线形扫描的添加使得检测器高度仅增加了因子1/(1-(r/R)2)。
已知算法对非理想源轨道的适应
上述精确的已知重构算法假设了理想的获取几何结构。然而,利用C臂装置的数据获取从未达到这些理想几何假设。获取系统的移动受到机械现象的影响,例如重力和惯性,导致不同的非理想类型的聚焦轨道-在重构方案中不得不考虑的现象。
真实世界C臂装置的非理想获取几何结构在该创造性方法中由齐次(homogenous)投影矩阵的序列Ps∈IR3×4表示。对于每个源位置s,并由此对于每个测量的投影图像,矩阵Ps完整地描述了对象的透视锥面光束投影。更准确地,该矩阵定义了对象的每一三维像素xh和对应检测器图像点的坐标wh之间的关系
wh=Psxh (10)
其中具有笛卡尔坐标向量x的三维像素由齐次向量xh-(b·xT,b)(其中b∈IR\0)表示,并且笛卡尔检测器位置w=(upix,vpix)T的模拟表示是wh=(c·wT,c)T(其中c∈IR\0)。upix和vpix是分别沿着与检测器DP(s)的像素网格的行或列方向相符的两条正交轴的向量eu,s和ev,s测量的图像点的坐标值。
因为获取几何结构的偏差在不同C臂装置之间变化、而在同一装置上对于相继扫描基本维持恒定,所以基本的任务在于确定用于给定C臂系统的单独有效的序列Ps。通过以下自动化过程而完成该几何校准:该自动化过程涉及具有精确定义的结构的校准幻象和对于给定s计算有效矩阵Ps的合适的校准算法。
可以在投影参数的全集中分解矩阵Ps。特别是,非本征参数是感兴趣的,因为它们包括所涉及的检测器和聚焦实体的位置和方位。该创造性方法从每一Ps计算像素坐标系统的轴的焦点位置和方向向量。由此可组成获取轨道,并且当使用从真实C臂装置下载的矩阵时,可以确定该获取几何结构与重构方案假设的理想圆形和线形几何结构的偏差。为了方便而引入了矩阵Ms∈IR3×3,其包括Ps的前三列。注意所有Ms是可逆的。
焦点位置y(s)被计算为:
投影矩阵仅按照缩放比例定义了检测器。为了具有关于C臂获取系统的精确结构的知识,需要焦点-检测器距离或检测器像素间隔的规格。然而,检测器像素坐标系统的两个轴的方向普遍有效。eu,s平行于向量((0,0,1)·Ms)T×((1,0,0)·Ms)T,而ev,s指向((0,0,1)·Ms)T×((0,1,0)·Ms)T方向。此外,光轴和检测器平面的交点的检测器坐标w0,s被计算为w0,s=Ms·(0,0,1)Ms。结果表明真实路径可比理想圆圈沿径向变化最多2%。此外,焦点位置没有定位在平面内,而是沿纵向变化。当C臂运动时,出现焦点和检测器的相对移动。倾斜和旋转偏差不是非常突出,但是在获取工作期间的检测器的平面内平移运动是显著的。为理想圆形和线形轨道导出的已知重构算法不得不适应于真实C臂装置的非理想源路径。
根据本发明,对于圆形和线形重构方案的应用,应用了以下一般策略。投影矩阵Ps精确描述了对于每一s的对象与其锥面光束投影图像之间的关系。由此,通过直接使用所述投影矩阵,后方投影步骤中的非理想获取几何结构的精确考虑是可能的。对于滤波步骤,非理想轨道的情况被近似转移到理想情况中。所假设的包括部分圆圈和正交连接的线段的理想轨道被拟合到真实世界焦点位置的集合中。所拟合的圆形路径再次作为抛物线而投影在检测器上,并且与上述相同的方案可用来将滤波方向确定为所发生抛物线的切线。
图5图示了理想圆形和线形轨道yfitted(s)在真实世界焦点位置集合y(s)中的拟合。有点的源位置被定位在圆面CP之下。
任何合适的成本函数可用于将理想轨道yfitted(s)拟合到路径y(s)。下面将描述最小方差拟合。如图5所示的理想轨道yfitted(s)到路径y(s)的最小方差拟合对应于总估计误差的最小化
并在三步骤方案中完成。
首先,执行平面到圆圈的焦点位置的最小方差代数拟合,随后执行路径y(s)在确定的圆面CP上的正交投影。在CP上,利用2D代数最小方差估计方法而将部分圆圈拟合到所投影的焦点位置,并然后最优地代表圆形扫描。最后,所述线段被确定为与圆面垂直,并连接到该圆圈段的一端。所拟合的轨道可通过圆面CP、圆心Xcenter、圆圈段的半径R、以及线段的长度和位置来描述。
应记住,空间坐标系统的位置和方位由投影矩阵Ps给出。执行坐标系统的归一化。这通过将每一Ps从右侧乘以与s无关的变换矩阵Tv∈IR4×4而完成。该归一化包括两个操作:平移Tv,t以将原点定位在Xcenter,以及旋转Tv,t以使线方向与x3轴平行。由此,Tv=Tv,r·Tv,t。在坐标变换之后,圆面CP等于x3=0的面,并且轨道以转轴为中心,所述线指向正x3方向。
此外,不得不处理焦点和检测器的相对位置的改变。改编该检测器坐标系统,使得拟合的源轨道被投影到与理想情况下相同的位置上的检测器上。然后,在无需任何进一步修正的情况下,可以应用已知反演算法的滤波指令。关于所使用的C臂硬件,校正与理想几何结构的情况有最显著的偏离的像素坐标系统的平面内平移运动就足够了。然而,可类似地处理任何其他几何偏离。对于每一s,确定平移矩阵Td,s∈IR3×3并从左边乘以Ps,使得空间坐标原点总是投影在相同检测器坐标上。由此,的坐标变为独立于检测器区域的当前错位。
最后,修正的投影矩阵
描述对于所述空间和检测器的涉及归一化坐标系统的锥面光束投影。应注意,这些投影矩阵仍然代表非理想获取几何结构。
实验结果
研究中评估了与真实世界C臂扫描对应的非理想获取几何结构,以验证本创造性适应方法。此外,注意到几何偏离对得到的图像质量的影响。
该实验的感兴趣的对象包括用数学方法定义的几何对象(如椭球或立方体),并且用嵌入的低对比度对象以及高对比度结构仿真包括均匀区域的人头的基本解剖结构。因为该组成,所以该所谓的数学头部幻象(参照例如“http://www.imp.uni-erlargen.de/phantoms/head/head.html”,头部幻象描述)对重构方案非常苛刻,从而显著地揭示任何类型的假象。
通过将对象沿着轴偏移4cm,而仿真两种扫描。第一种涉及理想几何结构,而第二种包括从真实世界C臂装置下载并由此代表相关几何偏离的一连串投影矩阵。表2中列出了重构和仿真参数。
在FDK重构图像中出现了严重的锥面假象,从而与圆形和线形方法得到的高质量图像数据形成对比。理想和非理想轨道的仿真扫描开始于不同的角位。由此,假象的方位也不同。关于非理想获取几何结构,检测两种重构方案中的附加假象。一些低密度的条纹状假象出现在高对比度骨头结构附近。该假象是由于角度采样中的一些轻微的不规则性以及滤波步骤中的一些残余的不精确性导致的。特别是,线形和圆形扫描贡献的匹配在更苛刻的情况下可能是关键性的。无论如何,该实验表明:该创造性改编用于考虑真实世界C臂装置的几何失真是足够的。
表2:重构和仿真参数
尽管本领域普通技术人员可提出修改和改变,但是发明人的意图是在这里许可的专利中,实施合理和正确地落入其对本领域的贡献的范围内的所有改变和修改。
Claims (5)
1.一种用于根据从对象获取的数据重构该对象的CT图像的方法,包括步骤:
操作C臂x射线设备,该设备具有其上安装有x射线源和辐射检测器的可旋转C臂,所述x射线源具有从中以锥面光束发射x射线的焦点,所述操作步骤通过以下步骤而完成:沿着包括实际不完全圆圈和与所述实际不完全圆圈的一端连接的实际直线段的实际聚焦轨道而围绕对象旋转所述焦点,并在所述辐射检测器上、对于所述实际聚焦轨道中的每一焦点位置、而检测从所述焦点发射并由所述对象衰减的辐射;
对于所述实际聚焦轨道中的所述焦点的每一位置,以电子方式计算投影矩阵,该矩阵对于该焦点位置而描述该辐射检测器上的对象的透视锥面光束投影;和
利用基于包括理想不完全圆圈和与所述理想不完全圆圈的一端连接的理想直线段的理想轨道的重构算法、并利用所述投影矩阵而使所述重构算法中的所述理想轨道适于所述C臂设备的所述实际聚焦轨道,而重构该对象的图像。
2.根据权利要求1的方法,其中所述对象包括多个三维像素,并且其中所述辐射检测器包括多个检测器图像点,该方法包括:以电子方式计算每一投影矩阵,以定义所述对象的每一三维像素和对应检测器图像点的坐标之间的关系。
3.根据权利要求1的方法,包括:操作所述C臂设备以沿着包括实际不完全圆圈和与所述实际不完全圆圈的一端连接的一根实际直线段的实际聚焦轨道而移动所述x射线源的所述焦点,并在所述重构算法中采用包括理想不完全圆圈和与所述理想不完全圆圈的一端连接的一根理想直线段的理想轨道。
4.根据权利要求1的方法,包括:操作所述C臂设备以沿着包括具有相对端的实际不完全圆圈和分别与所述实际不完全圆圈的所述两端连接的两根实际直线段的实际聚焦轨道而移动所述x射线源的所述焦点,所述两根实际直线段被安排在所述实际不完全圆圈的同一侧,并在所述重构算法中采用包括具有相对端的理想不完全圆圈和分别与所述不完全圆圈的所述两端连接的两根理想直线段的理想轨道,所述两根理想直线段被安排在所述理想不完全圆圈的同一侧。
5.根据权利要求1的方法,包括:操作所述C臂设备以沿着包括具有相对端的实际不完全圆圈和分别与所述实际不完全圆圈的所述两端连接的两根实际直线段的实际聚焦轨道而移动所述x射线源的所述焦点,所述两根实际直线段被安排在所述实际不完全圆圈的对侧,并在所述重构算法中采用包括具有相对端的理想不完全圆圈和分别与所述不完全圆圈的所述两端连接的两根理想直线段的理想轨道,所述两根理想直线段被安排在所述理想不完全圆圈的对侧。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/057,978 | 2005-02-15 | ||
US11/057,978 US7359477B2 (en) | 2005-02-15 | 2005-02-15 | Method for reconstructing a CT image using an algorithm for a short-scan circle combined with various lines |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1839758A CN1839758A (zh) | 2006-10-04 |
CN100528088C true CN100528088C (zh) | 2009-08-19 |
Family
ID=36815607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200610009268.7A Expired - Fee Related CN100528088C (zh) | 2005-02-15 | 2006-02-15 | 用圆-线快速扫描算法重构计算机断层分析图像的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US7359477B2 (zh) |
CN (1) | CN100528088C (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8086010B2 (en) * | 2006-06-30 | 2011-12-27 | Kabushiki Kaisha Toshiba | Medical image diagnosis apparatus and the control method thereof |
JP5539719B2 (ja) * | 2006-08-31 | 2014-07-02 | コーニンクレッカ フィリップス エヌ ヴェ | 画像形成システム |
DE102006044661B4 (de) * | 2006-09-21 | 2008-07-24 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Ermittlungsverfahren für endgültige Projektionsmatrizen |
EP2074383B1 (en) * | 2006-09-25 | 2016-05-11 | Mazor Robotics Ltd. | C-arm computerized tomography |
US20090213985A1 (en) * | 2007-10-31 | 2009-08-27 | Kabushiki Kaisha Toshiba | System and method for tomographic reconstruction utilizing circular trajectory and scanogram to reduce artifacts |
JP2011502683A (ja) * | 2007-11-19 | 2011-01-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 関心領域の画像を決定するイメージング装置、イメージング方法及びコンピュータプログラム |
US8401144B2 (en) * | 2008-08-07 | 2013-03-19 | Koninklijke Philips Electronics N.V. | Method and apparatus for correcting artifacts in circular CT scans |
US20100202583A1 (en) * | 2009-02-03 | 2010-08-12 | Ge Wang | Systems and Methods for Exact or Approximate Cardiac Computed Tomography |
DE102009020400B4 (de) * | 2009-05-08 | 2016-04-21 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Bildbestimmung aus beim Durchlaufen einer Trajektorie aufgenommenen Röntgenprojektionen |
DE102009031165A1 (de) * | 2009-06-30 | 2011-01-05 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Aufnahme von Röntgenbildern für eine dreidimensionale Bildrekonstruktion |
BR112012015836A2 (pt) * | 2010-01-13 | 2016-06-14 | Univ Australian | processo e sistema de formação de imagem por tomografia computadorizada, e, meio de armazenagem legível por computador |
CN102682199B (zh) * | 2012-04-28 | 2015-08-26 | 昆明理工大学 | 基于二维索引计算投影系数构造系统矩阵的简易方法 |
DE102012207910B4 (de) * | 2012-05-11 | 2021-09-23 | Siemens Healthcare Gmbh | Erweiterung der Trajektorie zur 3D Rekonstruktion |
US10702226B2 (en) | 2015-08-06 | 2020-07-07 | Covidien Lp | System and method for local three dimensional volume reconstruction using a standard fluoroscope |
US10777315B2 (en) | 2015-10-13 | 2020-09-15 | Mazor Robotics Ltd. | Global spinal alignment method |
CN106228584B (zh) * | 2016-07-20 | 2019-08-13 | 中国人民解放军信息工程大学 | 锥束ct圆加直线轨迹反投影滤波重建方法 |
US10699448B2 (en) | 2017-06-29 | 2020-06-30 | Covidien Lp | System and method for identifying, marking and navigating to a target using real time two dimensional fluoroscopic data |
WO2019075074A1 (en) | 2017-10-10 | 2019-04-18 | Covidien Lp | SYSTEM AND METHOD FOR IDENTIFICATION AND MARKING OF A TARGET IN A THREE-DIMENSIONAL FLUOROSCOPIC RECONSTRUCTION |
US10905498B2 (en) | 2018-02-08 | 2021-02-02 | Covidien Lp | System and method for catheter detection in fluoroscopic images and updating displayed position of catheter |
CN110974278B (zh) * | 2019-12-21 | 2022-02-11 | 电子科技大学 | 一种dsa锥束精确滤波反投影断层成像系统及成像方法 |
CN111707687A (zh) * | 2020-03-26 | 2020-09-25 | 国网河南省电力公司电力科学研究院 | 一种直线轨迹扫描ct检测方法及装置 |
CN113364315B (zh) * | 2021-05-08 | 2022-04-12 | 南昌大学 | 一种基于海鸥优化算法的三相逆变器控制方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2700909B1 (fr) * | 1993-01-27 | 1995-03-17 | Gen Electric Cgr | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. |
US5963612A (en) * | 1997-12-31 | 1999-10-05 | Siemens Corporation Research, Inc. | Apparatus for C-arm calibration for 3D reconstruction in an imaging system utilizing planar transformation |
US6379041B1 (en) * | 1998-11-02 | 2002-04-30 | Siemens Aktiengesellschaft | X-ray apparatus for producing a 3D image from a set of 2D projections |
US7197105B2 (en) * | 2001-08-16 | 2007-03-27 | Research Foundation Of The University Of Central Florida, Inc. | Efficient image reconstruction algorithm for the circle and line cone beam computed tomography |
JP4423959B2 (ja) * | 2003-12-16 | 2010-03-03 | 株式会社島津製作所 | 断層撮影装置 |
-
2005
- 2005-02-15 US US11/057,978 patent/US7359477B2/en not_active Expired - Fee Related
-
2006
- 2006-02-15 CN CN200610009268.7A patent/CN100528088C/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
Image reconstruction for the circle and line trajectory. A. Katsevich.Physics in Medicine and Biology,Vol.49 No.22. 2004 |
Image reconstruction for the circle and line trajectory. A. Katsevich.Physics in Medicine and Biology,Vol.49 No.22. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
US20060182216A1 (en) | 2006-08-17 |
CN1839758A (zh) | 2006-10-04 |
US7359477B2 (en) | 2008-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100528088C (zh) | 用圆-线快速扫描算法重构计算机断层分析图像的方法 | |
Navab et al. | 3D reconstruction from projection matrices in a C-arm based 3D-angiography system | |
US9013471B2 (en) | 3D X-ray imaging of coronary vessels with ECG gating and motion correction | |
US7734009B2 (en) | Angiographic x-ray diagnostic device for rotation angiography | |
JP4415762B2 (ja) | 断層撮影装置 | |
EP0520778B1 (en) | Tomographic image reconstruction using cross-plane rays | |
JP4606414B2 (ja) | 円錐形状光線束を用いるコンピュータ断層撮影方法 | |
CN102525501A (zh) | 图像处理设备和图像处理方法 | |
JPH06214031A (ja) | イメージ表示装置およびイメージ表示方法 | |
US8213565B2 (en) | Method for correcting truncated projection data | |
US20110019791A1 (en) | Selection of optimal views for computed tomography reconstruction | |
Rit et al. | Filtered‐backprojection reconstruction for a cone‐beam computed tomography scanner with independent source and detector rotations | |
JP2008541982A (ja) | コーンビームct用高速再構成アルゴリズム | |
Lee et al. | A Grangeat‐type half‐scan algorithm for cone‐beam CT | |
US7187747B2 (en) | Computerized tomography method with helical relative movement and conical beam | |
JP2001238877A (ja) | X線画像を形成するための方法及び装置 | |
JP2010501270A (ja) | 2つの傾斜円についてのコンピュータ断層撮影再構成 | |
Lee et al. | Artifacts associated with implementation of the Grangeat formula | |
JP6437163B1 (ja) | コンピュータ断層撮影画像生成装置 | |
Zhong et al. | A dual‐view digital tomosynthesis imaging technique for improved chest imaging | |
Hoppe et al. | Accurate image reconstruction using real C-arm data from a Circle-plus-arc trajectory | |
Yang et al. | FDK Half‐Scan with a Heuristic Weighting Scheme on a Flat Panel Detector‐Based Cone Beam CT (FDKHSCW) | |
Dennerlein et al. | Exact and efficient cone-beam reconstruction algorithm for a short-scan circle combined with various lines | |
US10881369B2 (en) | Systems and methods for enhanced radiography | |
Zheng et al. | Position coordinates-based iterative reconstruction for robotic CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
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: 20090819 Termination date: 20200215 |