CN1957847A - 再现对象的断层造影图像的方法和断层造影设备 - Google Patents
再现对象的断层造影图像的方法和断层造影设备 Download PDFInfo
- Publication number
- CN1957847A CN1957847A CN200610142407.3A CN200610142407A CN1957847A CN 1957847 A CN1957847 A CN 1957847A CN 200610142407 A CN200610142407 A CN 200610142407A CN 1957847 A CN1957847 A CN 1957847A
- Authority
- CN
- China
- Prior art keywords
- projection
- wave filter
- described method
- adopt
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000003325 tomography Methods 0.000 title abstract description 4
- 230000005855 radiation Effects 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000002591 computed tomography Methods 0.000 claims description 22
- 238000003384 imaging method Methods 0.000 claims description 13
- 238000002583 angiography Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 5
- 239000004744 fabric Substances 0.000 claims description 3
- 238000005481 NMR spectroscopy Methods 0.000 claims description 2
- 238000013500 data storage Methods 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims description 2
- 230000003278 mimic effect Effects 0.000 claims 1
- 239000011159 matrix material Substances 0.000 description 11
- 210000000481 breast Anatomy 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000002601 radiography Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- RKTYLMNFRDHKIL-UHFFFAOYSA-N copper;5,10,15,20-tetraphenylporphyrin-22,24-diide Chemical compound [Cu+2].C1=CC(C(=C2C=CC([N-]2)=C(C=2C=CC=CC=2)C=2C=CC(N=2)=C(C=2C=CC=CC=2)C2=CC=C3[N-]2)C=2C=CC=CC=2)=NC1=C3C1=CC=CC=C1 RKTYLMNFRDHKIL-UHFFFAOYSA-N 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. 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/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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/612—Specific applications or type of materials biological material
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种用于从运动的辐射源(11)穿过对象(18)而落在检测器(13)上的投影数据中再现对象(18)的断层造影图像的方法,其中在再现时对投影数据进行滤波和反投影,在辐射源(11)、检测器(13)和位于待扫描对象(18)位置处的测试对象(58)具有至少一种相同的空间布局的条件下通过测试投影和迭代分析再现技术确定滤波器,该滤波器在该给定布局下向该断层造影图像给出测试对象的投影数据的最佳滤波和反投影,在该给定布局下扫描位于该测试对象位置处的对象(18)并获得投影数据,利用该投影数据和根据步骤1.1确定的滤波器再现该断层造影图像。此外本发明还涉及一种用于执行该方法的断层造影设备。
Description
技术领域
本发明涉及一种用于从运动的辐射源落在检测器上的投影数据中再现对象的断层造影图像的方法,在该再现中对投影数据进行滤波和反投影。此外本发明还涉及一种采用不同的辐射来获取投影的断层造影设备。
背景技术
计算机断层造影(CT)为医学和检测技术提供了一种检查和测量方法,借助这种方法可以检查患者或受检对象的内部结构,而不必在此过程中介入患者体内或损坏受检对象。在此对所要检查的对象从不同角度拍摄多个投影,从这些投影中可以计算出该对象的三维描述。
普遍公知通过所谓的滤波反投影(Filtered Back Projection,FBP)来解决该问题,例如参考以下文献[Bu04]Buzug:Einführung in die Computertomographie.1.Auflage 2004.Springer-Verlag,ISBN3-540-20808-9和[KS84]Kak,Slaney:Principles of Computerized Tomographic Imaging.1987,IEEE Press.ISBN0-87942-198-3。FBP是一种高性能计算方法,其中对测量的投影进行滤波并反投影到该图像上。在该方法中图像质量取决于所采用的滤波器或卷积核。对于简单的扫描几何形状来说可以分析准确地给定该滤波器或卷积核。主要是其中以均匀的角度步长拍摄很多投影的圆形轨道。打破这一假设的更为复杂的拍摄几何形状在尝试分析地确定滤波器时会产生问题。对此的一个例子就是断层造影图像的合成,其中在最普通的情况下在一个空闲轨道上只能从有限的角度范围获得很少的投影。
对于这种再现问题适用迭代方法,如代数再现技术(ART)。对该技术参考以下文献[Bu04]Buzug:Einführung in die Computertomographie.1.Auflage 2004.Springer.ISBN3-540-20808-9;[KS84]Kak,Slaney:Principles of ComputerizedTomographic Imaging.1987,IEEE Press.ISBN 0-87942-198-3;[WZM04]T.Wu,J.Zhang,R.Moore,E.Rafferty,D.Kopans,W.Meleis,D.Kaeli:Digital TomosynthesisMammography Using a Parallel Maximum Likelihood Reconstruction Method,Medical Imaging 2004:Physics of Medical Imaging,Proceedings of SPIE VOL.,5368(2004)1-11。有利的是在这种ART中迭代方法不需要在FBP中需要的滤波器。但由于其迭代特性计算持续时间明显更长并因此在实践中通常无法应用。ART的另一个缺点在于与FBP不同无法利用该方法再现对象的部分区域。
还要补充说明专利申请US2005/0058240A1,其中发明人采用极度简化的启发式假设来分析计算再现滤波器。因此该方法局限于少数简单的拍摄几何形状和断层造影图像的合成。
因此存在要找到一种有效方法来从投影数据中再现对象的断层造影图像的问题,该方法一方面不对再现的必要计算效率提出太大的要求,而另一方面可以用于任意拍摄几何形状以及辐射源、检测器与对象之间在测量期间的相对运动。
发明内容
因此本发明要解决的技术问题在于提供一种方法来克服上述问题。
发明人已经认识到,在数字成像中、准确地说是在从一维或二维投影照片或来自一维投影的二维图像中获得三维图像和立体数据的过程中,可以从迭代方法、优选为ART中获得用于任意扫描几何形状的滤波器,然后将该滤波器用于根据FBP方法的再现。
为此在辐射源、检测器和作为待扫描对象的测试对象具有特定的空间布局的条件下借助对测量对象的投影用迭代再现技术对每个投影确定一个滤波器,该滤波器在给定的布局中将投影数据最佳地反投影到测试对象上。然后将这样确定的滤波器用于通过对检查对象的投影数据进行滤波和反投影来进行再现,该检查对象在给定布局下在位于测试对象的位置被扫描。
在此虽然正确滤波器的计算强度很大,但对于给定的扫描几何形状来说只有一次。采用该给定滤波器的FBP的计算开销与经典FBP没有什么不同。
与经典FBP相比,该方法提供了为任意扫描形体产生和应用与问题相匹配的滤波器的可能性。与ART相比显著提高了速度。此外该方法还提供了选择性地再现对象的部分区域(感兴趣区,ROI)的可能性,这在经典FBP中已公知。
原则上该方法既可以用于断层造影图像又可以用于断层造影合成,并且对所有的成像断层造影方法来说都可以与射线类型无关地使用。
根据这一认识,本发明人建议,改善本身公知的用于从运动的辐射源落在检测器上的投影数据中再现对象的断层造影图像的方法,其中为了再现对投影数据进行滤波和反投影,使得在辐射源、检测器和位于待扫描对象位置处的测试对象具有至少一种相同的空间布局的条件下通过测试投影和迭代分析再现技术确定滤波器,该滤波器在给定布局下向该断层造影图像给出测试对象的投影数据的最佳滤波和反投影,在给定布局下扫描位于测试对象位置处的对象并获得投影数据,利用该投影数据和确定的滤波器再现该断层造影图像。
按照本发明,可以再现对象的至少一幅二维截面图像或三维立体图像。此外为了确定滤波器和再现而采用一维投影或二维投影。
就拍摄几何形状而言不存在什么限制,从而可以例如以扇形辐射几何形状或锥形辐射几何形状来拍摄投影。通常辐射源和检测器在圆形轨道或螺旋轨道上相对于被扫描对象运动。在该方法中可以自由选择这些轨道。
还可以用测量的数据进行重整(Rebinning),从而用于再现的投影数据例如以平行射线几何形状出现。
本发明方法的另一个特征是,为了采集对象可以在<180度的角度范围内进行扫描。
此外在该方法中可以比较大的角度间隔来扫描对象,例如投影之间的角度间隔可以至少为2°。还可以各测量的投影之间的可变步长来扫描。
有利的,尤其是当辐射源的路径以非圆形变化时,为每个投影角确定自己的滤波器并用于再现。相应的还可以为辐射源和/或检测器的每个投影位置确定自己的滤波器并用于再现。该滤波器在最普通的情况下还可以与位置有关。
有利的,为了确定要采用的滤波器可以对至少一个测试对象进行测量,或者可以计算对模拟测试对象的投影。
测试对象应当包含尽可能多的位置频率。因此有利的是,例如采用一个或多个导线或小球体的设置。作为模拟测试对象合适的除了上述测试对象之外例如还有具有高斯形密度分布的噪声、球或棒。
为了避免在采用的各计算***中有太大的存储需要,有利的是从原始迭代确定的滤波器中通过其对位置和/或投影角的平均值来计算较少数量的新的平均滤波器。。如果采用这种数量减少的滤波器或者在开始时就只存在比在实际扫描期间测量投影的阶跃大的滤波器,则用于一个特定位置和/或一个特定投影的滤波器可以通过在属于其它位置和/或投影的滤波器之间进行插值来计算。
根据本发明的基本思想和上述方法,该方法可用于各种类型的断层造影设备、但尤其是用于从X射线成像、从磁共振成像、从超声波成像或从光学成像中获得投影的断层造影设备。此外对于这种断层造影设备发明人建议具有数据存储器来长期存储预先确定的滤波器。
尤其是应用于大型综合医院时,有利的是具有将获得的投影数据和存储的滤波器传送给其它分立图像计算机的程序。
附图说明
下面借助附图中的优选实施例详细描述本发明,尤其是还作为滤波器计算的数学基础。其中仅示出理解本发明所需的特征。在此采用以下附图标记:11:位于第一位置的X射线源;11’:位于另一位置的X射线源;12:第一投影的X射线束;12’:另一投影的X射线束;13:位于第一位置的检测器;13’:位于另一位置的检测器;14:再现场;15:分析计算机;16:显示单元;17:滤波器的存储器;18:对象/患者;21:测量的投影;22:滤波/卷积;23:滤波后的投影;24:滤波器;25:反投影:26:图像/立体数据;31:测量的投影;32:起始图像;33:第n次迭代之后变化的投影;34:第n次迭代之后再现的图像;35:所计算的投影的计算(投影器);36:确定计算的投影和测量的投影之差;37:形成差和变化后投影之和;38:变化后的投影的反投影(反投影器);41:测量的投影;42:迭代变化的投影;43:算法;44:滤波器;51:X射线源;53:检测器;54:再现区域;58:导线模型;61:X射线源;62:X射线束;63:检测器;65:分析计算机;66:显示单元;67:滤波器的存储器;68:胸部;69:压缩板。
图中具体示出:
图1示出具有X射线源的典型CT设置;
图2示出FBP计算规则的流程图;
图3示出迭代再现的流程图;
图4示出滤波器计算的示意性描述;
图5示出***X射线造影***的拍摄几何;
图6示出***X射线造影***。
具体实施方式
图1示出公知的典型CT设置,其具有位于第一位置的X射线源11,其发射出用于进行第一投影的X射线束,该射线束在穿透位于再现场14中的待检查对象、在此即患者18之后被位于该第一位置的检测器13检测到。检测器的数据传送到进行再现的分析计算机15中,然后显示在显示单元16上。X射线源11在此按照理想方式在圆形轨道上运动,其中从不同角度拍摄大量投影。在图1中还示出位于另一角度位置的X射线源11′,其中发射出用于进行另一投影的X射线束12′,然后该射线束12’被位于另一位置的检测器13’检测到。
这种CT设置的标准再现方法是滤波的反投影(FBP)。大多数情况下在计算之前还要对射线进行分类(重整),从而给出一组用平行射线以等间距的投影角度拍摄的投影。在最简单的情况下是平行的射线和圆形、等间距扫描的扫描轨道,此时FBP的计算规则分为两步,其在图2中示意性示出:
1.用对所有投影都相同的、频率特征为|w|的滤波器24对投影数据21进行卷积22,其中|w|是投影的傅立叶变换的频率的绝对值。
2.将滤波后的投影23反投影到已完成的图像或立体数据26上。
如果射线从焦点以锥形射出,则在卷积之前需要对数据进行依据几何形状的加权,或者需要重整为平行射线(rebinning)。此外还可能需要对滤波器进行依据几何形状的修改。该算法的详细描述在[KS84]Kak,Slaney:Principles ofComputerized Tomographic Imaging.1987,IEEE Press.ISBN 0-87942-198-3或[Bu04]Buzug:Einführung in die Computertomographie.1.Auflage 2004.Springer-Verlag.ISBN3-540-20808-9中。对普通的扫描几何形状不再能分析地给定具有代表性消耗的滤波器。
数学上FBP可以用下面的等式描述:
等式(1)X=RWY
其中X表示待再现对象的所要查找的向量,R表示反投影矩阵,Y表示测量的投影数据的向量。矩阵W包括滤波和加权的组合,并在下面简称为滤波器。
本发明的确定用于后面的FBP的最佳滤波器的迭代ART方法基于以下原理:将测量的投影与从已再现的对象中计算的投影进行比较,然后将误差用于校正对象的图像。在此第n次迭代Xn中的图像利用更新等式计算:
等式(2)Xn=Xn-1+RV(Y-PXn-1)
在该迭代开始时存在一幅合适的起始图像X0,如零图像。P在此表示***矩阵,利用该矩阵在已知扫描几何体的条件下从已扫描的对象图像中计算出投影。V是条件矩阵,利用它可以影响收敛速度。在最简单的情况下V是具有相同值、如1的对角矩阵。
在文献中用过的ART描述可以如下变换,其中Xn-1表示为“校正后数据”Yn-1的反投影:
等式(3)Xn-1=RYn-1,
从而等式(2)可改写如下:
Xn=RYn-1+RV(Y-PRYn-1)
等式(4)=R(Yn-1+V(Y-PRYn-1))
=R((1-VPR)Yn-1+VY)
因此对于Yn来说下式成立:
等式(5)Yn=(1-VPR)Yn-1+VY=KYn-1+VY
其中(6)K=(1-VPR)。
等式(5)中的递归表达可以展开为更清楚的表达
等式(7)
由此表明从Y中通过矩阵运算得到了Yn。
如果将等式(3)和等式(7)与等式(1)比较,则Un必须对应于等式(1)的FBP滤波器矩阵W,该矩阵在FBP中对数据进行滤波和加权。由此,通过对事先已知对象的迭代分析再现可以为特定的投影找到可以为该投影在给定扫描几何形状的条件下用于FBP的最佳滤波器。
迭代再现的过程在图3中示出。从起始图像32中借助投影器35获得计算的投影。然后建立计算的投影和测量的投影31之间的差值36,并在框37中与变化后的投影33相加。接着将该和一方面又存储到33中,另一方面反投影到38中。在第二次迭代时从这样获得的图像34中确定计算的投影。该迭代将一直进行下去,直到达到收敛或算法中断为止。该过程必须被这样初始化,开始时在33中存储投影,通过对该投影的简单反投影形成起始图像32。在零图像的情况下33只包含值0。
由于其结构这样获得的滤波器与要解决的再现问题相匹配。通常加权和滤波随位置变化,大多数情况下取决于正在观察的投影。但在很多情况下滤波器只随着位置和投影缓慢变化,从而通常足以通过对合适的位置区域或投影区域取平均来减少滤波器的数量。还可以在对应于不同位置区域或投影区域的滤波器之间进行插值。
迭代再现的经验显示,属于中间步骤的滤波器Un也能具有有利的再现特性。
由于投影矩阵和反投影矩阵的复杂性并首先由于它们的大小,多数情况下不能分析地或直接地计算U或Un。
根据本发明,滤波器的计算对预先给定的拍摄几何形状来说是通过确定与位置和投影有关的传输函数Un来进行的。
为了计算滤波器必须首先为给定的扫描几何形状确定投影Y。在一个应用中可以通过测量合适测量对象如CT设备中的细导线的投影来进行。在另一个应用中可以通过在期望的扫描几何形状中模拟相同的对象来获得投影。
利用这样获得的投影根据等式(5)迭代地再现要成像的对象。在此当达到期望的图像清晰度或期望的信噪比时中断迭代。在该迭代步骤之后可以借助普遍公知的方法通过比较Yn和Y来确定整体上与位置有关的滤波器Un。对此可参阅文献[OS75]Oppenheim,Alan V.Schafer,Ronald:Digital Signal Processing,Prentice Hall,1975,ISBN 0132146355。
下面给出确定Un的几个例子。如果测试对象例如是一个细导线,则可以将该导线***理论地表示为脉冲。然后可以容易地从校正后的投影Yn中获得与位置有关的传输函数Un。通过用多个相互错位的导线重复测量,可以与位置和投影有关地确定Un。
对于普通的测试对象例如可以借助短时傅立叶变换确定Un。
在另一个实施方式中投影还可以包含模拟或测量的噪声。这样例如可以借助局部自动相关函数来确定Un。为此可参阅文献[OS75]Oppenheim,Alan V.Schafer,Ronald:Digital Signal Processing,Prentice Hall,1975,ISBN 0132146355。
滤波器计算的示意描述在图4中示出。在此从测量的投影41及其由此迭代改变的投影42中利用与投影匹配的算法43计算滤波器44。
如上所述,本发明的方法例如还是用于在***X射线造影中的断层造影图像合成。例如,为此在图5中用示意性的三维视图示出在四方形的再现区域54内的X辐射源51、平面检测器53和由3个平行设置的导线组成的导线模型58的拍摄几何形状。利用这样的布局在一次测试测量中迭代地确定与相应投影对应的滤波器。
如果一次性确定滤波器或一组滤波器,则可以在相同的几何布局下扫描实际要被扫描的对象,并通过该滤波器借助非常快速的FBP来计算断层造影数据。
图6以截面示意图示出这种扫描情况。在此女性胸部68定位在平面监测器63上,并通过未详细示出的再现区域内的压缩板69来保持。X射线管61在用箭头示出的圆形轨道上围绕中心M运动,从而通过以该方式旋转的射线扇形62扫描胸部68,并通过具有多个在此未清楚示出的检测器元件的、固定的平面检测器确定多个投影角下的吸收值。根据事先拍摄的投影,在计算单元65中采用事先存放在存储器67中的滤波器,并利用该滤波器执行FBP。为此所需的程序Prgx同样可以存放在存储器67中并在需要时调用。再现后的断层造影数据可以基于高性能的方法在该扫描之后非常快速地直接显示在显示屏66上。
本发明的方法可用于图1的经典CT及其它采用多个检测器和辐射源来进行拍摄的设置中,例如在文献[KKK05]Kachelriess M.,Knaup M.,Kalender W.A.:Phase-Correlated Imaging from Multithreaded Spiral Cone-Beam CT Scans of theHeart,International Meeting on Fully Three Dimensional Image Reconstruction,SaltLake City;Utah,USA,July 6-9,2005;Proceedings of Fully 3D pp.159-162中所述。
应当理解,本发明的上述特征在不脱离本发明的范围的情况下不仅可以用于各种给出的组合,还能用于其它组合或单独使用。
文献:
Bu04Buzug:“Einführung in die Computertomographie”,1.Auflage 2004.Springer-Verlag.ISBN3-540-20808-9
KS84 Kak,Slaney:“Principles of Computerized Tomographic Imaging”,1987,IEEE Press.ISBN 0-87942-198-3
OS75 Oppenheim,Alan V.Schafer,Ronald:“Digital Signal Processing”,Prentice Hall,1975,ISBN 0132146355
WZM04:T.Wu,J.Zhang,R.Moore,E.Rafferty,D.Kopans,W.Meleis,D.Kaeli:“Digital Tomosynthesis Mammography Using a Parallel Maximum LikelihoodReconstruction Method”,Medical Imaging 2004:Physics of Medical Imaging,Proceedings of SPIE VOL.,5368(2004)1-11
KKK05:Kachelriess M.,Knaup M.,Kalender W.A.:“Phase-CorrelatedImaging from Multithreaded Spiral Cone-Beam CT Scans of the Heart”,International Meeting on Fully Three Dimensional Image Reconstruction,Salt LakeCity;Utah,USA,July 6-9,2005;Proceedings of Fully3D pp.159-162
Claims (32)
1.一种用于从运动的辐射源(11)穿过对象(18)而落在检测器(13)上的投影数据中再现该对象(18)的断层造影图像的方法,其中,在该再现中对投影数据进行滤波和反投影,其特征在于,
1.1.在辐射源(11)、检测器(13)和位于待扫描对象(18)位置处的测试对象(58)的具有至少一种相同的空间布局的条件下通过测试投影和迭代再现技术确定滤波器,该滤波器在该给定的布局下向该断层造影图像给出测试对象的投影数据的最佳滤波和反投影,
1.2.在该给定布局下扫描位于该测试对象位置处的对象(18)并获得投影数据,
1.3.利用该投影数据和根据步骤1.1确定的滤波器再现该断层造影图像。
2.根据权利要求1所述的方法,其特征在于,再现对象(18)的至少一幅二维截面图像。
3.根据权利要求1所述的方法,其特征在于,再现对象(18)的至少一幅三维立体图像。
4.根据权利要求1至3中任一项所述的方法,其特征在于,为了确定滤波器和为了再现而采用一维投影。
5.根据权利要求1至3中任一项所述的方法,其特征在于,为了确定滤波器和为了再现而采用二维投影。
6.根据权利要求1至5中任一项所述的方法,其特征在于,以锥形辐射几何形状拍摄投影。
7.根据权利要求1至5中任一项所述的方法,其特征在于,以扇形辐射几何形状拍摄投影。
8.根据权利要求1至5中任一项所述的方法,其特征在于,以平行射线几何形状拍摄投影。
9.根据权利要求1至8中任一项所述的方法,其特征在于,用测量的数据进行重整。
10.根据权利要求1至9中任一项所述的方法,其特征在于,根据步骤1.1确定的滤波器与位置有关。
11.根据权利要求1至10中任一项所述的方法,其特征在于,在<180度的角度范围内对对象(18)进行扫描。
12.根据权利要求1至11中任一项所述的方法,其特征在于,以投影之间至少2°的角度间隔来扫描对象(18)。
13.根据权利要求1至12中任一项所述的方法,其特征在于,以各测量的投影之间的可变步长来进行扫描。
14.根据权利要求1至13中任一项所述的方法,其特征在于,为每个投影角确定自己的滤波器并用于再现。
15.根据权利要求1至14中任一项所述的方法,其特征在于,为辐射源和/或检测器(13)的每个投影位置确定自己的滤波器并用于再现。
16.根据权利要求1至15中任一项所述的方法,其特征在于,为了确定要采用的滤波器对至少一个测试对象进行测量。
17.根据权利要求1至15中任一项所述的方法,其特征在于,为了确定要采用的滤波器计算模拟的测试对象的投影。
18.根据权利要求1至17中任一项所述的方法,其特征在于,为了确定要采用的滤波器采用导线(58)作为测试对象。
19.根据权利要求1至17中任一项所述的方法,其特征在于,为了确定要采用的滤波器采用小球体的设置作为测试对象。
20.根据权利要求1至17中任一项所述的方法,其特征在于,为了确定要采用的滤波器采用噪声图像作为测试对象。
21.根据权利要求1至17中任一项所述的方法,其特征在于,为了确定要采用的滤波器采用噪声作为给定投影。
22.根据权利要求1至21中任一项所述的方法,其特征在于,首先确定要采用的滤波器,然后存储并在以后才使用。
23.根据权利要求1至22中任一项所述的方法,其特征在于,从原始迭代确定的滤波器中通过其对位置和/或投影角的平均值来计算较少数量的新的平均滤波器。
24.根据权利要求1至23中任一项所述的方法,其特征在于,用于一个特定位置和/或一个特定投影的滤波器通过在属于其它位置和/或投影的滤波器之间进行插值来计算。
25.根据权利要求1至24中任一项所述的方法,其特征在于,采用多个检测器和/或辐射源。
26.一种用于从X射线成像中获得投影的断层造影设备,其特征在于,具有执行按照权利要求1至26中任一项的方法步骤的程序(Prgx)并在运行时执行该程序。
27.一种用于从磁共振成像中获得投影的断层造影设备,其特征在于,具有执行按照权利要求1至26中任一项的方法步骤的程序(Prgx)并在运行时执行该程序。
28.一种用于从超声波成像中获得投影的断层造影设备,其特征在于,具有执行按照权利要求1至26中任一项的方法步骤的程序(Prgx)并在运行时执行该程序。
29.一种用于从光学成像中获得投影的断层造影设备,具有执行按照权利要求1至26中任一项的方法步骤的程序(Prgx)并在运行时执行该程序。
30.根据权利要求26至29中任一项所述的断层造影设备,其特征在于,具有数据存储器(17),在其中长期存储预先确定的滤波器。
31.根据权利要求26至30中任一项所述的断层造影设备,其特征在于,具有将获得的投影数据和存储的滤波器传送给分立图像计算机的程序(Prgx)。
32.根据权利要求26至31中任一项所述的断层造影设备,其特征在于,该设备具有多个检测器和/或辐射源。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102005050917A DE102005050917A1 (de) | 2005-10-24 | 2005-10-24 | Verfahren und Tomographiegerät zur Rekonstruktion einer tomographischen Darstellung eines Objektes |
DE102005050917.7 | 2005-10-24 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1957847A true CN1957847A (zh) | 2007-05-09 |
CN1957847B CN1957847B (zh) | 2012-02-01 |
Family
ID=37905314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200610142407.3A Expired - Fee Related CN1957847B (zh) | 2005-10-24 | 2006-10-24 | 再现对象的断层造影图像的方法和断层造影设备 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7447295B2 (zh) |
JP (1) | JP5601748B2 (zh) |
CN (1) | CN1957847B (zh) |
DE (1) | DE102005050917A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103052970A (zh) * | 2010-08-03 | 2013-04-17 | 皇家飞利浦电子股份有限公司 | 从成像数据中移除对象支撑 |
CN109425843A (zh) * | 2017-08-25 | 2019-03-05 | 西门子保健有限责任公司 | 用于磁共振成像的方法和*** |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4617318B2 (ja) * | 2003-12-16 | 2011-01-26 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | フィルタリング後方投影を用いる画像化方法 |
DE102005051620A1 (de) * | 2005-10-27 | 2007-05-03 | Siemens Ag | Verfahren zur Rekonstruktion einer tomographischen Darstellung eines Objektes |
JP5011859B2 (ja) * | 2006-07-13 | 2012-08-29 | 株式会社島津製作所 | 放射線断層撮像装置 |
DE102006036327A1 (de) | 2006-08-03 | 2008-02-14 | Siemens Ag | Verfahren zum Bereitstellen von 3D-Bilddaten und System zum Aufnehmen von Röntgenbildern |
DE102008003173B4 (de) * | 2008-01-04 | 2016-05-19 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung für die Computertomographie zur Be- |
DE102008028387B4 (de) | 2008-06-13 | 2018-05-17 | Siemens Healthcare Gmbh | Tomographisches Bildrekonstruktionsverfahren zum Erzeugen eines Bildes von einem Untersuchungsobjekt und nach diesem Verfahren arbeitende bildgebende Einrichtung |
DE102009014726A1 (de) * | 2009-03-25 | 2010-10-07 | Siemens Aktiengesellschaft | Verfahren und Bildrekonstruktionseinrichtung zur Rekonstruktion von Bilddaten |
DE102009032060B4 (de) * | 2009-07-07 | 2016-12-01 | Siemens Healthcare Gmbh | Verfahren und Einrichtung zum Ermitteln der räumlichen Position und/oder Orientierung eines Objektes bei der Durchführung einer bildunterstützten medizinischen Maßnahme |
FR2953964B1 (fr) | 2009-12-15 | 2012-05-04 | Gen Electric | Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections |
DE102010019016B4 (de) * | 2010-05-03 | 2017-03-02 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von Bilddaten eines bewegten Untersuchungsobjektes aus Messdaten nebst zugehöriger Gegenstände |
GB201112359D0 (en) | 2011-07-19 | 2011-08-31 | Univ Antwerpen | Filter for tomographic reconstruction |
NL2007151C2 (en) * | 2011-07-20 | 2013-01-22 | Nucletron Bv | A gamma source tracking system. |
BR112014001596A2 (pt) | 2011-07-28 | 2017-02-21 | Koninklijke Philips Nv | aparelho de geração de imagem para gerar uma imagem de um objeto, método de geração de imagem para gerar uma imagem de um objeto e programa de computador para gerar uma imagem de um objeto |
JP6118324B2 (ja) | 2011-09-12 | 2017-04-19 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 制限角度トモグラフィーにおけるフィルターバックプロジェクションのための画像再構成方法 |
US8948337B2 (en) * | 2013-03-11 | 2015-02-03 | General Electric Company | Computed tomography image reconstruction |
US20160026894A1 (en) * | 2014-07-28 | 2016-01-28 | Daniel Nagase | Ultrasound Computed Tomography |
JP6486100B2 (ja) * | 2014-12-22 | 2019-03-20 | キヤノン株式会社 | 画像処理装置、画像処理方法、およびプログラム |
FR3040618B1 (fr) * | 2015-09-08 | 2022-07-29 | Orion France | Objet-test de controle qualite des images de mammographie numerique par tomosynthese |
US10893842B2 (en) | 2018-02-08 | 2021-01-19 | Covidien Lp | System and method for pose estimation of an imaging device and for determining the location of a medical device with respect to a target |
US10743822B2 (en) * | 2018-06-29 | 2020-08-18 | Carestream Health, Inc. | Fiducial marker for geometric calibration of bed-side mobile tomosynthesis system |
CN112085811B (zh) * | 2020-09-23 | 2021-04-23 | 赛诺威盛科技(北京)有限公司 | Ct局部重建的方法及装置 |
Family Cites Families (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4494141A (en) * | 1982-09-30 | 1985-01-15 | General Electric Company | Dual purpose array processor for a CT scanner |
US4649561A (en) * | 1983-11-28 | 1987-03-10 | Ben Arnold | Test phantom and method of use of same |
US4729100A (en) | 1984-08-28 | 1988-03-01 | Kabushiki Kaisha Toshiba | CT System which convolutes projection data with a frequency varying filter function |
JPS6156946A (ja) * | 1984-08-28 | 1986-03-22 | Toshiba Corp | 放射線断層検査装置 |
US5105451A (en) * | 1988-12-07 | 1992-04-14 | Eastman Kodak Company | Electrographic process utilizing fluorescent toner and filtered detector for generating an electrical image signal |
US5384861A (en) * | 1991-06-24 | 1995-01-24 | Picker International, Inc. | Multi-parameter image display with real time interpolation |
JP3409220B2 (ja) * | 1994-09-20 | 2003-05-26 | コニカ株式会社 | 画像読み取り装置 |
JPH11511375A (ja) * | 1996-03-13 | 1999-10-05 | アナロジック コーポレーション | コンピュータ断層撮影用リング抑圧フィルタ |
JP3373720B2 (ja) * | 1996-03-25 | 2003-02-04 | 株式会社日立メディコ | X線断層撮影装置 |
JP3657689B2 (ja) * | 1996-04-05 | 2005-06-08 | 株式会社日立メディコ | X線ct装置 |
JPH1021372A (ja) * | 1996-07-05 | 1998-01-23 | Toshiba Corp | X線ct装置 |
JP3673041B2 (ja) | 1996-11-21 | 2005-07-20 | ジーイー横河メディカルシステム株式会社 | Ctイメージング方法およびx線ct装置 |
US5909476A (en) * | 1997-09-22 | 1999-06-01 | University Of Iowa Research Foundation | Iterative process for reconstructing cone-beam tomographic images |
US5907593A (en) * | 1997-11-26 | 1999-05-25 | General Electric Company | Image reconstruction in a CT fluoroscopy system |
DE19845133A1 (de) * | 1998-10-01 | 2000-04-06 | Philips Corp Intellectual Pty | Computertomographie-Verfahren mit kegelförmigem Strahlenbündel |
US6307910B1 (en) * | 2000-01-07 | 2001-10-23 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for reduced radiation coronary computed tomography imaging |
US6674834B1 (en) * | 2000-03-31 | 2004-01-06 | Ge Medical Systems Global Technology Company, Llc | Phantom and method for evaluating calcium scoring |
JP3707347B2 (ja) * | 2000-04-07 | 2005-10-19 | 株式会社島津製作所 | X線ct装置の画像処理方法及びx線ct装置並びにx線ct撮影用記録媒体 |
DE10037491A1 (de) | 2000-08-01 | 2002-02-14 | Stryker Leibinger Gmbh & Co Kg | Verfahren zum dreidimensionalen Visualisieren von Strukturen im Körperinneren |
US6680995B2 (en) * | 2001-10-31 | 2004-01-20 | Ge Medical Systems Global Technology Co., Llc | Method and apparatus of determining and displaying a helical artifact index |
US7006594B2 (en) * | 2002-02-25 | 2006-02-28 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for reconstruction calibration of detector position and source motion based on a multi-pin phantom |
US6707878B2 (en) * | 2002-04-15 | 2004-03-16 | General Electric Company | Generalized filtered back-projection reconstruction in digital tomosynthesis |
US6529575B1 (en) * | 2002-04-29 | 2003-03-04 | Ge Medical Systems Global Technology Company, Llc | Adaptive projection filtering scheme for noise reduction |
JP3950782B2 (ja) * | 2002-05-22 | 2007-08-01 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 3次元逆投影方法およびx線ct装置 |
DE10244181A1 (de) | 2002-09-23 | 2004-04-01 | Siemens Ag | Verfahren zur Bilderstellung in der Computertomographie und CT-Gerät zur Durchführung des Verfahrens |
JP4142482B2 (ja) * | 2003-04-04 | 2008-09-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP4439202B2 (ja) * | 2003-05-09 | 2010-03-24 | 株式会社東芝 | X線コンピュータ断層撮影装置及び画像ノイズシミュレーション装置 |
JP4495926B2 (ja) * | 2003-07-01 | 2010-07-07 | 株式会社東芝 | X線立体再構成処理装置、x線撮影装置、x線立体再構成処理方法及びx線立体撮影補助具 |
US6987829B2 (en) * | 2003-09-16 | 2006-01-17 | General Electric Company | Non-iterative algebraic reconstruction technique for tomosynthesis |
JP2005095329A (ja) * | 2003-09-24 | 2005-04-14 | Toshiba Corp | 超解像処理装置及び医用画像診断装置 |
US7430269B2 (en) * | 2003-10-14 | 2008-09-30 | Koninklijke Philips Electronics N.V. | Iterative data reconstruction |
US6973157B2 (en) * | 2003-12-23 | 2005-12-06 | General Electric Company | Method and apparatus for weighted backprojection reconstruction in 3D X-ray imaging |
US7042974B2 (en) * | 2004-05-20 | 2006-05-09 | Eastman Kodak Company | Method for reconstructing tomographic images |
US8774355B2 (en) * | 2004-06-30 | 2014-07-08 | General Electric Company | Method and apparatus for direct reconstruction in tomosynthesis imaging |
-
2005
- 2005-10-24 DE DE102005050917A patent/DE102005050917A1/de not_active Withdrawn
-
2006
- 2006-10-23 US US11/584,544 patent/US7447295B2/en not_active Expired - Fee Related
- 2006-10-23 JP JP2006287630A patent/JP5601748B2/ja not_active Expired - Fee Related
- 2006-10-24 CN CN200610142407.3A patent/CN1957847B/zh not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103052970A (zh) * | 2010-08-03 | 2013-04-17 | 皇家飞利浦电子股份有限公司 | 从成像数据中移除对象支撑 |
CN103052970B (zh) * | 2010-08-03 | 2015-11-25 | 皇家飞利浦电子股份有限公司 | 从成像数据中移除对象支撑 |
CN109425843A (zh) * | 2017-08-25 | 2019-03-05 | 西门子保健有限责任公司 | 用于磁共振成像的方法和*** |
Also Published As
Publication number | Publication date |
---|---|
DE102005050917A1 (de) | 2007-04-26 |
US20070093711A1 (en) | 2007-04-26 |
JP5601748B2 (ja) | 2014-10-08 |
US7447295B2 (en) | 2008-11-04 |
CN1957847B (zh) | 2012-02-01 |
JP2007117736A (ja) | 2007-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1957847A (zh) | 再现对象的断层造影图像的方法和断层造影设备 | |
JP6564254B2 (ja) | X線ct装置及びその制御方法 | |
JP5019193B2 (ja) | スキャン対象に関するct画像内の再構成点における画像データ値を決定する再構成方法及びx線コンピュータ断層撮影装置 | |
O’Brien et al. | Recent advances in X-ray cone-beam computed laminography | |
JP6824633B2 (ja) | X線コンピュータ断層撮影装置、逐次近似再構成方法および医用画像処理装置 | |
US7706497B2 (en) | Methods and apparatus for noise estimation for multi-resolution anisotropic diffusion filtering | |
US20070116175A1 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
US11670017B2 (en) | Systems and methods for reprojection and backprojection via homographic resampling transform | |
KR20100133950A (ko) | 동적인 제약들에 따른 오브젝트 주변의 사용을 통한 단층 촬영에 있어서의 양 감소 및 이미지 강화 | |
Wu et al. | Spatial-temporal total variation regularization (STTVR) for 4D-CT reconstruction | |
WO2010135901A1 (zh) | 伪双能欠采样物质识别***和方法 | |
Melot et al. | Some proximal methods for Poisson intensity CBCT and PET | |
US9619902B2 (en) | Filter for tomographic reconstructions | |
JP2016152916A (ja) | X線コンピュータ断層撮像装置及び医用画像処理装置 | |
US20090086888A1 (en) | X-ray ct apparatus and image reconstruction method | |
Li et al. | Sparse CT reconstruction based on multi-direction anisotropic total variation (MDATV) | |
CN1742296A (zh) | 锥形光束计算机断层成像的立体像素驱动的螺旋重建方法 | |
CN115777114A (zh) | 针对ct图像去噪的3d-cnn处理 | |
Park et al. | A fully GPU-based ray-driven backprojector via a ray-culling scheme with voxel-level parallelization for cone-beam CT reconstruction | |
JP5214110B2 (ja) | X線ct装置 | |
KR20170087320A (ko) | 단층 영상 생성 장치 및 그에 따른 단층 영상 복원 방법 | |
Hou et al. | Analysis of compressed sensing based CT reconstruction with low radiation | |
CN103886627B (zh) | 带有保持边缘的滤波的ct图像重建 | |
Kong et al. | Analytic reconstruction approach for parallel translational computed tomography | |
WO2022165117A1 (en) | Simultaneous ct scanning and reconstruction |
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 |
Granted publication date: 20120201 Termination date: 20141024 |
|
EXPY | Termination of patent right or utility model |