CN1765323A - 半径内图像依存检测器列滤波方法、x射线ct装置和程序 - Google Patents
半径内图像依存检测器列滤波方法、x射线ct装置和程序 Download PDFInfo
- Publication number
- CN1765323A CN1765323A CN200510054839.4A CN200510054839A CN1765323A CN 1765323 A CN1765323 A CN 1765323A CN 200510054839 A CN200510054839 A CN 200510054839A CN 1765323 A CN1765323 A CN 1765323A
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- mentioned
- filtering
- 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.)
- Granted
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims description 21
- 230000001419 dependent effect Effects 0.000 title description 3
- 238000002591 computed tomography Methods 0.000 claims description 35
- BEWYHVAWEKZDPP-UHFFFAOYSA-N bornane Chemical compound C1CC2(C)CCC1C2(C)C BEWYHVAWEKZDPP-UHFFFAOYSA-N 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- XETQTCAMTVHYPO-UHFFFAOYSA-N Isocamphan von ungewisser Konfiguration Natural products C1CC2C(C)(C)C(C)C1C2 XETQTCAMTVHYPO-UHFFFAOYSA-N 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 229930006742 bornane Natural products 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000002596 correlated effect Effects 0.000 claims 1
- 230000006866 deterioration Effects 0.000 abstract 1
- 230000002265 prevention Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 19
- 238000002224 dissection Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000003860 storage Methods 0.000 description 5
- 238000013480 data collection Methods 0.000 description 4
- 238000002203 pretreatment Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 229910052709 silver Inorganic materials 0.000 description 3
- 239000004332 silver Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004171 remote diagnosis Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- 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/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
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4021—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pulmonology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种根据计算机断层摄影扫描取得数据的方法,从在CT***中具有多个检测器列的检测器(103)中取得投影数据,在检测器(103)的列方向上对该投影数据进行滤波,在列方向滤波单元(113)中取得风车假象被削减后的已滤波数据,从已滤波数据中用再构成推荐值(114)来再构成图像数据。能够在防止z轴析象度降低的同时,实现抑制风车假象和改善图像质量。
Description
技术领域
本发明涉及多层(multi-slice)X射线计算机断层摄像***的滤波方法、X射线CT装置和程序。
背景技术
多层X射线计算机断层摄影(CT)***在1998年开发并进入医疗市场。现在,层数大致是从2到40的范围,但有望增加到64或者256(参照非专利文献1、非专利文献2以及非专利文献3。各自的内容通过参照在本说明书中一体化)。使用了CT***用螺旋扫描的圆锥射束图像再构成算法使用了导入柔性聚焦轨道在费罗德莰普(Feld kamp)算法中适用加权功能的、一般被加权的费罗德莰普再构成算法(参照专利文献1及专利文献2以及从非专利文献4到非专利文献10。各自的内容通过参照在本说明书中一体化)。
在其他的再构成算法中,有再构成未被排列在单一轴上的多个层(multi-slice)的、高度单层重束(ASSR)和适合型多平面再构成(AMPR)等的准圆锥射束算法。层的轴与旋转轨道一致,z滤波它可以得到垂直层(slice)(参照专利文献3以及从非专利文献1到非专利文献13。它们全部的内容通过参照在本说明书中一体化)。
通过使用圆锥平行扇型射束重束(rebinning)(例如,参照从非专利文献14到非专利文献17。这些全部内容通过参照在本说明书中一体化)以及扇形射束算法(参照非专利文献18,其内容通过参照在本说明书中一体化),具有同样混淆现象(aliasing)问题的其他算法有数个。本发明并未被再构成方案的选择所限定。
在圆锥射束或者准圆锥射束算法中共同的问题之一是在z轴(检测器列)方向上的采样间隔不充分。这些算法违反乃奎斯特定理(在1个检测器单元开口内需要2个试料),起因于高频成分,在再构成的图像内产生混淆(aliasing)假象(artifact)(参照非专利文献19。其内容通过参照在本说明书中一体化)。这些混淆假象作为风车假象(windmillartifact)为人们所知。目前,有克服该问题的2种方法。一种方法是改***件的方法,另一种方法是改***件。
硬件的变更在z轴方向上适用飞点聚焦技术(xy平面,在检测器通道方向),把2组投影与在z中的N个试料组合设置成在z中的2个N个试料的1组(参照非专利文献20。该内容通过参照在本说明书中一体化)。但是,其潜在要解决对策存在以下问题。(1)数据大小是现在使用的数据大小的2倍。(2)扫描器价格更高。(3)需要新的图像再构成算法。
在采用软件(software)的解决对策中,有削减z析象度(z中的高频)的不适合型z滤波,和适合型的对象依存z滤波。在采用不适合型滤波的解决对策中,有在具有固定核心的投影数据域中的z滤波(参照专利文献4。其内容通过参照和本说明书一体化),在图像域中的z滤波,以及在具有可变尺寸的核心的投影数据域中的z滤波(参照专利文献5。其内容通过参照在本说明书中一体化)。
这些解决对策的一个缺点是同样失去图像整体的z析象度。简单地说,z析象度到处丢失。肋骨、头盖骨以及脊柱等,风车(windmill)假象(artifact)原因的几乎全部位于图像的周边区域上,图像的中央区域需要高空间z析象度。
在适合型z滤波中,有在投影数据域中的滤波(参照J.Hsieh,非专利文献21。其内容通过参照在本说明书中一体化)以及图像域中的滤波。这些方法都要根据z中的梯度等的对象依存索引变更z滤波的核心,需要大量处理,而这是不理想的。
[专利文献1]荒馆博,南部恭二郎,“CT装置”日本专利第2,825,352号
[专利文献2]K.Taguchi,“X-ray Computed TomographyApparatus,”U.S.Patent No.5,825,842(1995年申请)
[专利文献3]町田好男,“计算机断层摄影装置”,日本专利申请公开第08-187240号
[专利文献4]K.Taguchi,U.S.Pat.No.5,825,824(1998)
[专利文献5]I.Zmora,U.S.Pat.No.6,560,308
[非专利文献1]Y.Saito,H.Aradate,H.Miyazaki,K.Igarashi,and H.Ide,“Development of a Large Area 2-dimensional Detector for Real-Time3-dimensional CT(4D-CT),”Radiology vol.217(P),405(2000)
[非专利文献2]Y.Saito,H.Aradate,H.Miyazaki,K.Igarashi,and H.Ide,“Large Area Two-Dimensional Detector System for Real-TimeThree-Dimensional CT(4D-CT),”Proc.Of SPIE Med.Imag.Conf,4320,775-782(2001)
[非专利文献3]http://www3.toshiba.co.jp/medical/4d-ct/
[非专利文献4]L.A.Feldkamp,L.C.Davis,and J.W.Kress,“Practical Cone-Beam Algorithm,”L.Opt.Soc.Am.A,6,612-19(1984)
[非专利文献5]L.G.Zeng and G.T.Gullberg,“Short-scan ConeBeam Algorithm for Circular and Noncircular Detector Orbit,”Proc.OfSPIE Med.Imag.Conf,1233,453-463(1990)
[非专利文献6]H.Kudo and T.Saito,“Three-DimensionalHelical-Scan Computer Tomography Using Cone-Beam Projections,”J.Electron.Information Commun.Soc.Japan,J74-D-II,1108-1114(1991)
[非专利文献7]G.Wang,T.H.Lin,P.C.Cheng,D.M.Shinozaki,“AGeneral Cone-Beam Reconstruction Algorithm,”IEEE Trans.Med.Imaging,12,486-496(1993)
[非专利文献8]K.Wiesent,K.Barth,N.Novab,et al.,“Enhanced3-D-Reconstruction Algorithm for C-arm Systems Suitable forInterventional Procedures,”IEEE Trans.Med.Imaging,19,391-403(2000)
[非专利文献9]M.D.Silver,K.Taguchi,and K.S.Han,“Field-of-View Dependent Helical Pitch in Multi-Slice CT,”Proc.OfSPIE Med.Imag.Conf.,4320,839-850(2001)
[非专利文献10]M.D.Silver,K.Taguchi,and I.A.Hein,“A SimpleAlgorithm for Increased Helical Pitch in Cone-Beam CT,”The SixthInternational Meeting on Fully Three-Dimensional Image Reconstructionin Radiology and Nuclear Medicine,70-73(2001)
[非专利文献11]M.Kachelriess,S.Schaller,W.A.Kalender,“Advanced Single-slice Rebinning in Cone-Beam Sprial CT,”MedicalPhysics vol.27,pp.754-772(2000)
[非专利文献12]S.Schaller,K.Stierstorfer,H.Bruder,M.Kachelriess,and T.Flohr,“Novel Approximate Approach forHigh-Quality Image Reconstruction in Helical Cone Beam CT atArbitrary Pitch,”Proc of SPIE Vol.4322,pp.113-127(2001)
[非专利文献13]T.Flohr,K.Stierstorfer,H.Bruder,J.Simon,A.Polacin,and S.Schaller,“Image Reconstruction and Image QualityEvaluation for a 16-slice CT scanner,”Medical Physics vol.30,pp.832-845(2003)
[专利文献14]H.Tuy,“3D Image Reconstruction for Helical PartialCone-beam Scanner,”proc of Fully 3D 1999,pp.7-10
[非专利文献15]H.Turbell,et al.,“Three-dimensional ImageReconstruction in Circular and Helical Computed Tomography,”Licentiate thesis No.760,Linkoping Univ,ISBN 91-7219-463-4,1999
[非专利文献16]H.Turbell,et al.,“An Improved PI-method forReconstruction from Helical Cone-beam Projections,”Conf record ofIEEE MIC 1999
[非专利文献17]R.Manzke et al.,“Extended CardiacReconstruction(ECR):A Helical Cardiac Cone-beam ReconstructionMethod,”proc of Fully 3D 2003,Mo-PM2-4
[非专利文献18]K.Taguchi and H.Anno,“High TemporalResolution for Multi-slice Helical CT,”Medical Physice vol.27,May 2000
[非专利文献19]M.Silver,K.Taguchi,I.Hein,B.Chiang,M.Kazama,I.Mori,“Windmill Artifact in Multi-Slice Helical CT,”Proc ofSPIE Vol.5032,pp.1918-1927(2003)
[非专利文献20]T.Flohr.H.K.Bruder,K.Stierstorfer,S.Schaller,“Evaluation of Approaches to Reduce Spiral Artifacts in Multi-SliceSpiral CT,”RSNA 2003 program,pp567
[非专利文献21]J.Hsieh,“Adaptive Interpolation Approach forMulti-slice Helical CT Reconstruction,”Proc of SPIE Vol.5032,pp.1876-1833(2003)
发明内容
本发明的目的在于:在防止z轴析象度降低的同时实现抑制假象(artifact)和改善图像质量。
本发明的第1非限定形态提供取得包含从CT***中的至少2个检测器列取得投影数据的步骤,至少在2个检测器列的方向上对该投影数据进行滤波,取得削减了风车假象的已滤波数据的步骤,以及从已滤波数据再构成图像数据的步骤的、从计算机断层摄影(CT)扫描中取得数据的方法。
本发明的另一形态提供一种X射线CT装置,包含,螺旋扫描装置,其构成是在台架和床台中的至少一方沿着该床台的轴方向移动期间,收集投影数据,具备构成为发生X射线而构成的X射线源,和沿着轴方向被配置成至少2个检测器列,具有为生成投影数据而构成的检测器部件的检测器;处理器,具备至少在2个检测器列的方向上对投影数据进行滤波,为了取得削减了风车假象的已滤波数据而构成的滤波装置,和为再构成已滤波数据而构成的再构成装置。
本发明的另一形态提供一致X射线CT装置,包含:螺旋扫描装置,其构成是在台架和床台中的至少一方沿着该床台的轴方向移动期间,收集投影数据,包含:具备为发生X射线而构成的X射线源,和沿着轴方向被配置成至少2个检测器列,具有为生成投影数据而构成的检测器部件的检测器;处理器,具备至少在2个检测器列的方向上对投影数据进行滤波,来取得削减了风车(windmill)假象的已滤波数据的装置,和为再构成已滤波数据而构成的再构成装置。
本发明的另一形态提供让计算机执行以下步骤的程序:从CT***中的至少2个检测器列中取得投影数据的步骤;至少在2个检测器列的方向上滤波其投影数据,取得风车(windmill)假象被削减后的滤波后的步骤;以及再构成滤波后的数据的步骤。
本发明的另一形态提供计算机断层摄影中的图像数据的取得方法,即,从CT***中的至少2个检测器列中取得投影数据,在上述至少2个检测器列的方向上对上述投影数据进行加权并加法计算,为了使旋转中心的图像清晰度高,离开上述旋转中心越远则图像清晰度越低,而使上述加权依存通道变化,从经上述加权加算的投影数据中再构成图像数据。
本发明的另一形态提供具备以下装置的X射线计算机断层摄影装置:从CT***中的至少2个检测器列中取得投影数据的装置;在上述至少2个检测器列的列方向上对上述投影数据进行加权并加法计算的装置;为了使旋转中心的图像清晰度高,离开上述旋转中心越远则图像清晰度越低,而使上述加权依存通道变化的装置;从经上述加权加算的投影数据中再构成图像数据的装置。
如果可以参照附图并参照以下详细说明进一步理解本发明,则可以更容易完整理解本发明以及随之的许多优点。
如果采用本发明则在防止z轴析象度降低的同时,可以实现抑制假象和改善图像质量。
附图说明
图1(a)是表示N=4的多层检测器的图,图1(b)是表示几何学配置的图。
图2(a)是表示σ(r2D)的例子的图,图2(b)是表示高斯分布w(k)的图。
图3(a)是距离定义的图,图3(b)是距离定义的图,图3(c)是距离定义的图。
图4是表示WGn(0)的例子的图。
图5(a)是表示在RF软件以及RF清晰算法两方面中使用的函数Wz L例子的图,图5(b)是在RF软件算法中使用的函数Wz r例子的图。
图6是表示在检测器中投影的象素的图。
图7是z析象度指数的图。
图8是本实施方式的X射线计算机断层摄影装置的构成图。
具体实施方式
图8用方框图表示本实施方式的X射线计算机断层摄影装置(X射线CT装置)的构成。旋转扫描装置具备台架100、在该台架100的摄影区域内***被检测体的床台。台架100具有旋转环102,在该旋转环102上圆锥射束型X射线管101和多层形X射线检测器103相对配置。圆锥形射束X射线管101的构成是接收从高电压发生器109经由套环106周期发生的高电压脉冲,以四棱锥形发出X射线。多层形X射线检测器103用电离箱形检测器箱或者半导体检测器构成,如果是可以同时检测多层(多列)投影数据那样的电离箱形X射线检测器,则以圆锥射束形X射线管101的焦点(圆锥射束的顶点)为中心,被形成圆弧形状的多通道型的电离箱形X射线检测器沿着和旋转环102的旋转轴大致平行的方向排列多个,另外如果是半导体X射线检测器,则以圆锥射束的顶点为中心在圆弧状上排列多个X射线检测原件,而且,该X射线检测原件列沿着和旋转环102的旋转轴大致平行的方向排列多个。而且,把与旋转轴方向平行的方向称为“分层方向”或者检测器列方向,另外把电离箱形X射线检测器和X射线检测原件列的圆弧的方向称为“通道方向”。在多层X射线检测器103上,连接一般被称为DAS(data acquisition system)的数据收集电路104。在该数据收集电路104中,对每个通道上设置把多层形X射线检测器103的各通道的电流信号变换为电压的I-V变换器;把该电压信号与X射线的辐射周期同步地周期积分的积分器;放大该积分器的输出信号的放大器;把该前置放大器的输出信号变换为数字信号的模拟数字变换器。在前处理单元105中,对于在该数据收集电路104中检测出的数据(称为原始数据),补正通道间的灵敏度不均匀,另外执行补正X射线强吸收体,主要是金属部分产生的极端的信号强度降低或者信号脱落等的前处理。
列方向滤波单元113的列方向滤波典型地包含在前处理中,为了方便分开说明。列方向滤波单元113把在前处理单元105中接收前处理的原始数据在列方向上滤波。例如,在把相邻3列作为对象时,加权加算通道号码相同的3通道的原始数据。各列的加权w在从加权发生单元115向列方向滤波单元113中发生。在本实施方式中,加权不是一定的、而是沿着通道方向变化。在中央通道中,给予中央列的原始数据的加权高,给予两侧列的原始数据的加权低。从中央通道向周边通道,给予中央列的原始数据的加权逐渐降低,给予两侧列的原始数据的加权逐渐提高。具体的说明在以后叙述。
再构成装置114根据在前处理单元9中接收前处理,在列方向滤波单元113中接收列方向滤波的数据(为了和原始数据区别,称为投影数据)再构成图像,把该图像显示在显示器116上。
扫描控制器110为了执行投影数据的收集动作,即扫描,总控与扫描器有关的制控制器,使得在以定速稳定地使旋转环102旋转的同时,从高电压发生装置109中以一定周期发生高电压脉冲,接着与该高电压脉冲同步地用数据收集电路104周期收集投影数据等。
说明列方向滤波单元113的列方向滤波的概要。
例如在列方向上加权加算附近3列的原始数据。加权(也称为系数)从检测元件(通道)对应的象素的旋转中心开始的依存于半径r变化。即,加权和观测仪以及列一同还依存于通道方向变化。在旋转中心中因滤波特性清晰度高,即中央列的数据影响高,随着离开旋转中心,清晰度降低。如果把加权组合表示为(相对一方列的原始数据的加权,相对中央列的原始数据的加权,相对另一方列的原始数据的加权),则在旋转中心(中心通道)上,用(0,1,0)给予,在最终通道上用(0.3,0.4,0.3)给予。其间从(0,1,0)向(0.3,0.4,0.3)逐渐变化。
在焦点位置接近分层面时清晰度高(清晰),即单纯依存于通道(用γ计算半径r)。与加权相应的滤波特性的变移点如图5(a)示例的那样,用z轴方向的距离z0(=2mm)定义。
3列(3点)滤波的中心加权在列方向上对称。这是因为为了抑止射线和层面不交叉的状况的影响而限幅和焦点-交叉象素的距离L的缘故。检测器的上下某一单边数据因为在关注图像的再构成中不使用所以有问题。在该限幅使用的列接近检测器中心列时也有平滑效果。
在其他例子中也假定加权加算附近3列。加权依存于检测原件对应的象素从旋转中心开始的半径r(参照图3(a))。因而加权依存于观测仪、列、以及通道。滤波特性随着从旋转中心离开,从高的清晰度(0,1,0)向低的清晰度(0.3,0.4,0.3)逐渐变化。在焦点位置接近层面时是低的清晰度,即在基准值r0下相对全部通道大致统一为一定值。
这样,核心的清晰度随着距离r2D(例如,从同中心向xy平面上的直至总和值投影的投影距离)的减少而增大,同中心附近的象素具有比周边还高的z析象度。因此,在周边区域上的风车(windmill)假象的原因通过适用与平滑核心对应的检测器通道可以抑止。
因此,可以在图像再构成前取得在检测器列方向上的滤波数据。滤波器的核心可以定义为光线角度的函数。即,可以根据从同中心向xy平面上直至总和值投影的投影距离调整核心。从下式(1)求所希望的函数。
[式1]
在式(1)中,Pin是投影数据,v表示与投影角度β对应的投影数指数,ch是与光线角度γ对应的检测器通道指数,row是与圆锥角度α对应的检测器列指数,WGn是z滤波器的系数,2K-1表示核心的大小。
从下式(2)到(4)表示核心。
[式2]
在从式(2)到(4)中,R表示从焦点到象素中心的距离,γm表示最大光线角度,Nch表示1列的检测器通道的数,Cch是与γ=0对应的检测器通道指数,r2D表示从象素中心向xy平面的直至总和值投影的投影距离,σ(r2D)定义核心的清晰度。图2(a)表示σ(r2D)的例子,图2B表示使用不同的σ的w(k)。
通过该方法,随着r2D减少核心的清晰度增大,得到同中心附近的象素比周边高的z析象度。因此,在周边区域上风车(windmill)假象通过适用与平滑核心对应的检测器通道抑止(例如,通过用平滑化滤波)。
接着,滤波后的数据用于图像再构成。图像再构成技术包含滤波补正反投影和高速傅立叶变换(FFT),但并不限于这些。如果采用本实施方式,则γ滤波在适用平滑核心前进行。γ滤波包含沿着灯滤波等的图像再构成用的检测器通道适用旋转(convolution)工序的步骤,但并不限于此。
图3A表示总和值投影、旋转轴(z轴),以及再构成的图像(反投影数据/总和值投影的三维象素(voxel)/象素)的侧面图。图3B是从同中心向xy平面上的直至总和投影的投影距离。表示2次元距离r2D。图3c表示从同中心到再构成的xy平面(特定的z位置)与总和值投影一致的总和值投影上的三维象素的距离的3维距离r3D。
作为非限定的代替策略,式5(以下叙述)表示固定核心(2K+1)的大小的情况。
[式3]
wGn(k,ch)|k≠0=(1-wGn(0,ch))/(2K) (5)WGn(0)的例子表示在图4中。
作为非限定的代替对策,在不是空间域而是频率域中,在检测器列方向上可以适用滤波工序。因此,可以把FFT对每个通道适用在检测器列方向上。接着,在用FFT变换的数据上,可以乘表示z滤波核心的频率。最后,可以在已滤波数据上适用反FFT。
而且,采用本实施方式通过把γ滤波和列滤波组合,还可以适用2维FFT。在检测器数据中适用2维FFT后,可以一边把平滑滤波适用在列方向上,一边在γ方向上适用灯滤波。
作为非限定的代替实施方式,在图像再构成前在检测器列方向上可以滤波数据。滤波器的核心可以根据从同中心到与检测单元对应的象素的距离变更。下式6表示滤波式。
[式4]
从下式(7)到(16)说明核心。
[式5]
z=-CS·β/2π=-CS·(v-vC)/NvRev (16)
在从式(7)到(16)中,rowc表示在α=0中的检测器列指数,d是至同中心的投影检测单元的高度,vc是β=0(例如,在再构成的平面上有焦点时)中的投影数指数,CS是每旋转1圈的进给(螺旋螺距),NvRev是每转1圈的投影数。
在β小时也可以有2个非限定的变形。即,是RF清晰和RF软件,在此,“RF”表示“半径依存滤波方案”。使用RF清晰算法,把wz r固定在1.0。如图5A所示,wz L是z(β)的函数。式(17)以及(18)表示RF清晰算法。
[式6]
RF软件算法如图5A、5B所示,是Wz L以及wz r都是z(β)的函数的情况。具体地说,图5A表示在RF软件以及RF清晰算法两方面使用的函数wz L的例子。图5B表示在RF软件算法中使用的函数wz r的例子。下式(19)以及(20)表示RF软件算法。
[式7]
(20)
wGn→wGn(γ=γ(r0))(fixed for all)
图6表示检测器单元和再构成的象素的关系。通过图6所示的关系,发明者想到了r3D的概念。在图6中,象素对于β[-π,π]被投影在检测器上。横轴表示检测器通道,纵轴表示检测器列。如图6所示,投影区域和β一同变化。图7表示使用了RF清晰算法以及RF软件算法的z析象度指数。
如计算机技术业内人员明白那样,本实施方式的全部实施方式可以根据本实施方式的指点使用经编程的以往的通用计算机或者微处理器很顺利地实施。如软件技术的从业人员明白那样,适合的软件根据本揭示的指点具有一般技能的编程人员容易制作。
如相互参照的在美国专利第6,236,051中揭示的那样,计算机可以实施本实施方式的方法。计算机机箱收纳CPU、存储器(例如,DRAM,ROM,EPROM,EEPROM,SRAM,SDRAM,以及闪光RAM),以及其他任意的专用逻辑器件(例如,ASICS)或者包含可以构成的逻辑器件(例如,GAL以及可以再编程FPGA)的母板。计算机还包含多个输入装置(例如,键盘和鼠标)和监视器控制用的显示卡。而且,计算机用适合的设备总线(例如,SCSI总线、扩展IDE总线,或者超DMA总线)连接。可以包含软盘驱动器、其他的可取下存储器件(例如,小型盘、带,以及可以取下的高磁性介质),以及硬盘或者其他固定高密度介质驱动器。计算机可以与同样的设备总线或者其他设备总线连接。还可以包含小型盘阅读器、小型盘阅读器/记录器单元,或者小型盘自动播放器。
在与本实施方式相关的计算机可以读取的存储介质的例子中,有小型盘、硬盘、软盘、带、光磁盘、PROM(例如,EPROM,EEPROM,闪光EPROM),DRAM,SRAM,SDRAM等。本实施方式具备存储在这些计算机可以读取的存储介质之一或者它们的组合中,控制计算机的硬件并且计算机可以和人对话的软件。该软件包含设备驱动、操作***,以及开发工具等的用户应用程序,但并不限于此。在本实施方式的计算机程序产品中,有存储有如果由计算机执行则让计算机进行本实施方式的方法的计算机程序命令(例如,计算机代码设备)的、任意的计算机可以读取的存储介质。本实施方式的计算机代码设备可以是能任意解释或者可以执行的代码机构,它包含脚本、翻译机、动态连接程序库、Java类别(class),以及可以完全执行的程序,但并不限于此。而且,本实施方式的处理部分为了更高的性能、可靠性,以及/或者成本可以分散(例如,(1)在多个CPU之间,或者(2)在至少1个CPU和至少1个可以构成的逻辑器件之间)。例如,用第1计算机选择轮廓或者图像,为了远程诊断可以发送到第2计算机。
本实施方式为了补正图像对比度、不规则性的程度、纹理特征等,还可以用追加的滤波技术和工具补充。
如专业人员容易明白那样,本实施方式还可以通过应用程序特定集成电路的制作,或者通过相互连接以往构成电路的适合的网络执行。
本实施方式的图像数据信号源可以是X射线摄影装置、CT装置,以及MRI装置等任意的适合的图像取得装置。另外,已取得的数据如果还不是数字形式,则可以进行数字化,或者,取得并处理的图像数据的信号源可以是存储用图像取得装置制成的数据的存储器。存储器可以是本地存储器或者远程存储器,这种情况下,使用PACS(图像存档计算机***)等的数据通信网络,为了采用本实施方式进行处理可以访问图像数据。
而且,本实施方式并不局限于上述实施方式,在实施阶段中在不脱离其主旨的范围内可以变形构成要素来具体化。另外,可以通过被揭示在上述实施方式中的多个构成要素的适合的组合来形成各种发明。例如,可以从实施方式所示的全部构成要素中删除几个构成要素。而且,也可以适当组合涉及不同实施方式的构成要素。
Claims (22)
1、一种根据计算机断层摄影扫描取得数据的方法,包含:
从CT***中的至少2个检测器列中取得投影数据的步骤;
在上述至少2个检测器列的方向上对上述投影数据进行滤波,来取得削减了风车假象的已滤波数据的步骤;和
根据上述已滤波数据来再构成图像数据的步骤。
2、如权利要求1所述的根据计算机断层摄影扫描取得数据的方法,上述滤波与从光线角度和同中心到检测器单元的距离中的至少一方相关联来进行。
3、如权利要求1所述的根据计算机断层摄影扫描取得数据的方法,上述滤波与从同中心到总和值投影与xy平面一致的该总和值投影上的三维象素的距离相关来进行。
4、如权利要求1所述的根据计算机断层摄影扫描取得数据的方法,上述再构成包含反投影。
5、如权利要求4所述的根据计算机断层摄影扫描取得数据的方法,上述反投影包含:适用费罗德莰普再构成、适用高度单层重束、FFT、适用反FFT、以及适合型多平面再构成中的至少1种。
6、一种X射线CT装置,包括:
螺旋扫描装置,为了在台架和床台中的至少一方沿着该床台的轴方向移动期间收集投影数据而构成,它包含:为发生X射线而构成的X射线源;和沿着轴方向被配置成至少2个检测器列,具有为生成投影数据而构成的检测器部件的检测器;
处理器,包含至少在2个检测器列的方向上对上述投影数据进行滤波,为了取得削减了风车假象的已滤波数据而构成的滤波装置;和
为再构成上述已滤波数据而构成的再构成装置。
7、如权利要求6所述的X射线CT装置,上述滤波装置的构成为:根据从光线角度和同中心到检测器单元的距离中的至少一方来对上述投影数据进行滤波。
8、如权利要求6所述的X射线CT装置,上述滤波与从同中心到总和值投影与xy平面一致的该总和值投影上的三维象素的距离相关联来进行。
9、如权利要求6所述的X射线CT装置,上述再构成装置包含反投影装置。
10、如权利要求9所述的X射线CT装置,上述反投影装置的构成是:通过适用费罗德莰普再构成、高度单层重束、FFT、适用反FFT、以及适合型多平面再构成中的至少1种,来对上述已滤波数据进行反投影。
11、一种X射线CT装置,包含:
螺旋扫描装置,为了在台架和床台中的至少一方沿着该床台的轴方向移动期间收集投影数据而构成,该螺旋扫描装置包含:
为发生X射线而构成的X射线源;和
沿着轴方向被配置成至少2个检测器列,具有为生成投影数据而构成的检测器部件的检测器;
处理器,它包含至少在2个检测器列的方向上对投影数据进行滤波,来取得削减了风车假象的已滤波数据的装置;和
为再构成上述已滤波数据而构成的再构成装置。
12、如权利要求11所述的X射线CT装置,上述滤波装置的构成为:根据从光线角度和同中心到检测单元的距离中的至少一方来对上述投影数据进行滤波。
13、如权利要求11所述的X射线CT装置,上述滤波与从同中心到总和值投影与xy平面一致的该总和值投影上的三维象素的距离相关联来进行。
14、如权利要求11所述的X射线CT装置,上述再构成装置包含反投影装置。
15、如权利要求14所述的X射线CT装置,上述反投影装置的构成是:通过适用费罗德莰普再构成、高度单层重束、FFT、适用反FFT、以及适合型多平面再构成中的至少1种,来对上述已滤波数据进行反投影。
16、一种程序,使计算机执行以下步骤:
从CT***中的至少2个检测器列中取得投影数据;
在上述至少2个检测器列的方向上对上述投影数据进行滤波,来取得削减风车假象的已滤波数据;
再构成上述已滤波数据。
17、如权利要求16所述的程序,上述滤波与从光线角度和同中心到检测器单元的距离中的至少一方相关联来进行。
18、如权利要求16所述的程序,上述滤波与从同中心到总和值投影与xy平面一致的该总和值投影上的三维象素的距离相关联来进行。
19、如权利要求16所述的程序,上述再构成包含反投影。
20、如权利要求19所述的程序,上述反投影包含适用费罗德莰普再构成、适用高度单层重束、FFT、适用反FFT、以及适合型多平面再构成中的至少1种。
21、一种计算机断层摄影中的图像数据取得方法,
从CT***中的至少2个检测器列取得投影数据;
在上述至少2个检测器列的列方向上对上述投影数据进行加权并加法计算;
为了使旋转中心的图像清晰度高,离开上述旋转中心越远则图像清晰度越低,而使上述加权依存通道而变化,
根据上述加权并加法计算后的投影数据来再构成图像数据。
22、一种X射线计算机断层摄影装置,包括:
从CT***中的至少2个检测器列取得投影数据的装置;
在上述至少2个检测器列的列方向上对上述投影数据进行加权并加法计算的装置;
为了使旋转中心的图像清晰度高,离开上述旋转中心越远则图像清晰度越低,而使上述加权依存通道而变化的装置;和
根据加权并加法计算后的投影数据来再构成图像数据的装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/812,187 US20050226365A1 (en) | 2004-03-30 | 2004-03-30 | Radius-in-image dependent detector row filtering for windmill artifact reduction |
US10/812,187 | 2004-03-30 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1765323A true CN1765323A (zh) | 2006-05-03 |
CN100469317C CN100469317C (zh) | 2009-03-18 |
Family
ID=35060531
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200510054839.4A Active CN100469317C (zh) | 2004-03-30 | 2005-03-17 | 半径内图像依存检测器列滤波方法、x射线ct装置和程序 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20050226365A1 (zh) |
JP (2) | JP2005279282A (zh) |
CN (1) | CN100469317C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103228215A (zh) * | 2011-09-08 | 2013-07-31 | 株式会社东芝 | 多层ct装置及数据预处理方法 |
CN109887051A (zh) * | 2019-02-28 | 2019-06-14 | 沈阳开普医疗影像技术有限公司 | Ct图像重建反投影过程中的线性插值方法 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6252965B1 (en) * | 1996-09-19 | 2001-06-26 | Terry D. Beard | Multichannel spectral mapping audio apparatus and method |
US20030171665A1 (en) * | 2002-03-05 | 2003-09-11 | Jiang Hsieh | Image space correction for multi-slice helical reconstruction |
US7760924B2 (en) * | 2002-11-27 | 2010-07-20 | Hologic, Inc. | System and method for generating a 2D image from a tomosynthesis data set |
JP4611168B2 (ja) * | 2005-10-07 | 2011-01-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像再構成方法、およびx線ct装置 |
JP4509971B2 (ja) * | 2006-06-09 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP5007982B2 (ja) * | 2006-06-22 | 2012-08-22 | 国立大学法人東北大学 | X線ct装置及び同装置の画像再構成方法並びに画像再構成プログラム |
JP4350738B2 (ja) * | 2006-10-27 | 2009-10-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線断層撮影装置およびアーチファクトの低減方法 |
JP4414420B2 (ja) * | 2006-10-27 | 2010-02-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線断層撮影装置およびアーチファクトの低減方法 |
US8116426B2 (en) * | 2008-11-11 | 2012-02-14 | Kabushiki Kaisha Toshiba | Computed tomography device and method using circular-pixel position-adaptive interpolation |
JP5711241B2 (ja) | 2009-10-06 | 2015-04-30 | コーニンクレッカ フィリップス エヌ ヴェ | コーンビームct画像におけるアーチファクト低減のための方法 |
WO2014165455A2 (en) * | 2013-04-04 | 2014-10-09 | Illinois Tool Works Inc. | Helical computed tomography |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11511375A (ja) * | 1996-03-13 | 1999-10-05 | アナロジック コーポレーション | コンピュータ断層撮影用リング抑圧フィルタ |
JPH1021372A (ja) * | 1996-07-05 | 1998-01-23 | Toshiba Corp | X線ct装置 |
US5708690A (en) * | 1996-10-11 | 1998-01-13 | General Electric Company | Methods and apparatus for helical image reconstruction in a computed tomography fluoro system |
US6295331B1 (en) * | 1999-07-12 | 2001-09-25 | General Electric Company | Methods and apparatus for noise compensation in imaging systems |
JP4025530B2 (ja) * | 2001-10-24 | 2007-12-19 | 株式会社日立メディコ | X線ct装置 |
JP2003263420A (ja) * | 2002-03-12 | 2003-09-19 | Hitachi Medical Corp | 医療機器を用いた並列演算処理システム |
US6587537B1 (en) * | 2002-04-01 | 2003-07-01 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for multi-slice image reconstruction |
-
2004
- 2004-03-30 US US10/812,187 patent/US20050226365A1/en not_active Abandoned
-
2005
- 2005-03-17 CN CN200510054839.4A patent/CN100469317C/zh active Active
- 2005-03-30 JP JP2005099446A patent/JP2005279282A/ja active Pending
-
2013
- 2013-10-30 JP JP2013226074A patent/JP5661897B2/ja active Active
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103228215A (zh) * | 2011-09-08 | 2013-07-31 | 株式会社东芝 | 多层ct装置及数据预处理方法 |
CN103228215B (zh) * | 2011-09-08 | 2015-12-23 | 株式会社东芝 | 多层ct装置及数据预处理方法 |
CN109887051A (zh) * | 2019-02-28 | 2019-06-14 | 沈阳开普医疗影像技术有限公司 | Ct图像重建反投影过程中的线性插值方法 |
CN109887051B (zh) * | 2019-02-28 | 2023-07-25 | 沈阳开普医疗影像技术有限公司 | Ct图像重建反投影过程中的线性插值方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5661897B2 (ja) | 2015-01-28 |
US20050226365A1 (en) | 2005-10-13 |
JP2005279282A (ja) | 2005-10-13 |
CN100469317C (zh) | 2009-03-18 |
JP2014054553A (ja) | 2014-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1765323A (zh) | 半径内图像依存检测器列滤波方法、x射线ct装置和程序 | |
EP1643446B1 (en) | Image reconstruction method and x-ray computed tomography apparatus | |
US8094910B2 (en) | Method of reconstructing an image function from Radon data | |
CN1292712C (zh) | 多模态成像方法和装置 | |
US8284892B2 (en) | System and method for image reconstruction | |
US6574297B2 (en) | System and method for image reconstruction in a cone beam imaging system | |
JP5142664B2 (ja) | X線コンピュータ断層撮影装置 | |
CN1781454A (zh) | X-射线ct设备 | |
JP4612347B2 (ja) | コーンビーム再構成装置およびコンピュータ断層撮影装置 | |
CN1251975A (zh) | 章动断层面ct图象重建 | |
JP4612294B2 (ja) | X線コンピュータ断層装置、x線コンピュータ断層装置制御方法、及びx線コンピュータ断層撮影プログラム | |
US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
CN1830392A (zh) | X射线计算机断层摄影装置、图象处理装置和图象处理方法 | |
CN1489976A (zh) | 运动检查对象的计算机断层造影中设置图像的方法及设备 | |
US6459756B1 (en) | System and method for providing normalization correction for image reconstruction in a reduced pitch spiral scan cone beam computed tomography imaging system | |
JP2008529637A (ja) | 螺旋相対運動及び円錐ビーム束を用いたコンピュータ断層撮影方法 | |
US7778387B2 (en) | Reconstruction method for helical cone-beam CT | |
CN108242066B (zh) | Ct图像的空间分辨率增强装置和方法以及ct成像*** | |
US7539281B2 (en) | Methods and apparatus for reconstruction in helical cone beam volumetric CT | |
Xie et al. | Optimization of data acquisition in axial CT under the framework of sampling on lattice for suppression of aliasing artifacts with algorithmic detector interlacing | |
EP1099961A2 (en) | Image processing for medical imaging | |
KR101379316B1 (ko) | 다중 팬빔 나선형 ct 시스템 및 이를 이용한 엑스선 단층촬영방법 | |
CN101077302A (zh) | 计算机断层摄影装置 | |
Tang et al. | Extending three‐dimensional weighted cone beam filtered backprojection (CB‐FBP) algorithm for image reconstruction in volumetric CT at low helical pitches | |
CN1305298C (zh) | 图像重构方法和x射线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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160804 Address after: Japan Tochigi Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |