CN100579453C - X射线测量装置 - Google Patents
X射线测量装置 Download PDFInfo
- Publication number
- CN100579453C CN100579453C CN200580046148A CN200580046148A CN100579453C CN 100579453 C CN100579453 C CN 100579453C CN 200580046148 A CN200580046148 A CN 200580046148A CN 200580046148 A CN200580046148 A CN 200580046148A CN 100579453 C CN100579453 C CN 100579453C
- Authority
- CN
- China
- Prior art keywords
- ray
- data
- filter
- projection
- thickness
- 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
- 238000005259 measurement Methods 0.000 claims abstract description 63
- 238000012937 correction Methods 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims description 34
- 230000008569 process Effects 0.000 claims description 27
- 238000007689 inspection Methods 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 description 11
- 238000010521 absorption reaction Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 229920006395 saturated elastomer Polymers 0.000 description 7
- 238000002372 labelling Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000003050 axon Anatomy 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- -1 pottery Substances 0.000 description 1
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 description 1
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. 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 for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
Abstract
本发明提供一种X射线测量装置,其对检查对象照射X射线,并检测与检查对象相关的测量数据,并且在X射线源与检查对象之间设置调整X射线透过量的过滤器,对于X射线源与检查对象,使相对位置变化,并对得到的测量数据进行运算处理,将测量数据对数变换从而得到投影数据,并得到与获得的投影数据对应的所述过滤器的X射线吸收量,对于获得的X射线吸收量使用规定的变换式计算出过滤器的厚度,并根据计算出的过滤器的厚度得到与获得的投影数据对应的修正系数,将得到的修正系数与所述投影数据相乘,重建运算乘以了所述修正系数的投影数据,从而得到三维图像。
Description
技术领域
本发明涉及可抑制检测器的饱和,并得到提高了值的均匀性的良好的三维图像的X射线测量装置。
背景技术
具有X射线源与二维X射线检测器相对而配置在C字形的支柱(以下称为C臂)的开放的两端上X射线测量装置。具有从顶棚吊挂C臂的构造或在地面上支撑C臂的构造。此外,具有X射线源与二维X射线检测器相对而配置在桶架上的X射线测量装置。在这些装置中,通过使C臂或桶架旋转,可使X射线和检测器对在被摄体周围旋转,同时进行X射线检测。此外,对由旋转测量得到的多个测量数据分别实施修正处理,从而得到用于三维重建的一组投影数据,并对得到的一组投影数据使用三维重建算法进行重建处理,从而可得到三维图像。这些三维测量记载于非专利文献1中。
非专利文献1:新医疗,2002年10月,Vol.29,No.10,pp.102-105
在以往的X射线测量装置中,如果测量厚的被摄体,则透过被摄体到达检测器的X射线量变少,测量的X射线图像的画质低劣。如果采用强化从X射线源照射的X射线量,或缩短X射线源与检测器的距离等方法,使到达检测器的X射线量增加,则能够提高被摄体区域的画质。但是,另一方面,在不存在被摄体的区域,或者被摄体薄的区域中,到达检测器的X射线量增大,在X射线量超过检测界限的检测元件中,产生饱和现象,无法显示正确的值。如果为了抑制检测器的饱和,而在X射线源与被摄体之间设置由金属等制成的过滤器,则由于过滤器的种类或厚度,入射到被摄体的X射线的能量分布不同。如果使用能量分布不同的X射线图像进行重建处理,则存在得到的三维图像不能表示正确的值的问题。
此外,在未设置过滤器的情况或过滤器的种类或厚度一样的情况下,入射到被摄体的能量分布一样,但由于被摄体的种类或厚度,在透过被摄体的过程中,X射线的能量分布不同。如果使用能量分布不同的X射线图像进行重建处理,则存在所得到的三维图像不能表示正确的值的问题。
此外,从过滤器或被摄体产生的散射X射线入射到检测器并混杂在数据中。如果使用散射X射线混杂的X射线图像进行重建处理,则存在得到的三维图像不能表示正确的值的问题。
发明内容
本发明的目的在于针对这些问题,提供一种可修正三维图像的值的不均匀性,得到良好的三维图像的X射线测量装置。
为了达到上述目的,提供一种X射线测量装置,其具有:X射线源,其产生向检查对象照射的X射线;X射线检测器,其以隔着检查对象的方式与该X射线源相对配置,并将检查对象的透过X射线作为测量数据进行检测;过滤器,其设置在X射线源与检查对象之间,对X射线的透过量进行调整;保持装置,其保持X射线源和X射线检测器;旋转装置,其使X射线源及X射线检测器在检测对象的周围旋转;以及控制处理装置,其对由X射线检测器检测的测量数据进行运算处理,测量数据是X射线检测器以通过旋转装置而旋转的、且相对于检查对象的多个角度来进行检测的;该X射线测量装置的特征在于,控制处理装置对测量数据进行对数变换而得到投影数据,并得到与得到的投影数据对应的过滤器的X射线吸收系数,对得到的X射线吸收系数使用规定的变换式来计算过滤器的厚度,并根据计算出的过滤器的厚度得到与得到的投影数据对应的修正系数,将得到的修正系数与投影数据相乘,重建运算乘以了修正系数的投影数据,从而得到三维图像。
根据本发明,能够得到消除了X射线的检测元件的饱和现象、由过滤器或被摄体导致的能量分布的变化、或者从过滤器或被摄体发生的散射X射线所产生的影响的三维图像。
附图说明
图1是表示实施例一的不均匀性修正处理的步骤的流程图。
图2是表示在实施例一的不均匀性修正处理中使用的变换表A的图。
图3是表示在实施例一的不均匀性修正处理中使用的变换表B的图。
图4是以剖面的形式表示适用本发明的X射线测量装置的一例的侧视图的概念图。
图5(A)-(C)是表示将过滤器410沿X射线源401的旋转面切开的剖面(左侧)和立体图(右侧)的三种例子的图,(A)是整体为同样厚度的过滤器410的剖面图和立体图,(B)是由凹面形成的区域和同样厚度的区域组合而成的过滤器410的剖面图和立体图,(C)是由凹面形成的区域、由凸面形成的区域、和同样厚度的区域组合而成的过滤器410的剖面图和立体图。
图6(A)-(C)是表示过滤器410的其他形状的剖面图及俯视图,(B)是从X射线的入射方向观察的俯视图,(A)是在A-A位置沿箭头方向观察的剖面图,(C)是在B-B位置沿箭头观察的剖面图。
图中:401-X射线源;402-检测器;403-支柱;404-旋转装置;405-卧台;406-控制处理装置;407-旋转轴;408-被摄体;409-存储装置;410-过滤器。
具体实施方式
图4是以剖面的形式表示适用本发明的X射线测量装置的一例的侧视图的概念图。X射线测量装置包括:X射线源401、检测器402、支柱403、旋转装置404、卧台405、控制处理装置406。X射线源401与检测器402设置在支柱403上。支柱403使用C字形的臂、“コ”字形的臂或桶架等。在图4中,表示C字形的臂。并考虑从顶棚吊挂支柱403的方式或在地面上支撑支柱403的方式。支柱403通过旋转装置401,以旋转轴407为中心,在横置于卧台405上的被摄体408的周围旋转。
图4表示了作为最普通的方式,旋转轴407及卧台405与地面平行的情况。通过使支柱403及卧台405移动,也可将旋转轴407相对于体轴倾斜地设定。此外,也可以是旋转轴407相对于地面垂直,支柱403在站立或坐在椅子上的被摄体408的周围旋转的方式。
图4表示了X射线源401与检测器402通过旋转装置404以旋转轴407为中心在被摄体408的周围旋转的情况。也可以是X射线源401与检测器402固定,被摄体408通过旋转装置404以旋转轴407为中心旋转的方式。或者也可以是X射线源401、检测器402与被摄体408以旋转轴407为中心旋转的方式。无论哪种情况,通过相对的旋转运动实现同样的旋转测量。
检测器402使用平面形X射线检测器、X射线图像增强器和CCD照相机的组合、成像板、CCD检测器、固体检测器等。作为平面形X射线检测器,具有将光电二极管与非结晶型硅TFT作为一对,将其配置在正方矩阵上,并且将其与荧光体直接组合的检测器等。
从X射线源401产生的X射线透过被摄体408,通过检测器402变换为与X射线强度对应的电信号,并作为测量数据输入控制处理装置406中。控制处理装置406控制X射线源401中X射线的发生、检测器402中数据的取得、旋转装置404中支柱403的旋转。由此,X射线测量装置可进行如下的旋转测量:旋转支柱403,同时进行X射线的产生与测量数据的取得。控制处理装置406对于测量数据可执行对数变换处理或重建运算处理等,并取得三维数据。在图4中,控制处理装置406安装在检测器402的外部,但也可以是安装在检测器402的内部的方式。该情况下,可高速处理。
本发明所述的X射线测量装置可在X射线源401与检测器402之间设置过滤器410。过滤器410由铝、铜、黄铜等金属、陶瓷、树脂等构成。最简便的也可以是向塑料或树脂等的壳体中注入液体而制成。过滤器的形状任意,后述的数据修正处理可对应任意的过滤器形状。
实施例一
以下,对本发明的第一实施例进行说明。
本发明的控制处理装置406执行对三维图像的值的不均匀性进行修整的处理。控制处理装置406在内部具有存储装置409,存储修正处理需要的表等。作为控制处理装置406的输入机构,虽未图示,但可考虑来自键盘的键输入、来自文件的读入、存储芯片的更换。控制处理装置406作为操作菜单具有输入修正处理的执行的有无的模式或者开关等。
以下,表示修正三维图像的值的不均匀性的处理的步骤。
(1)灵敏度数据取得:在未设置被摄体的状态下,照射X射线进行测量,并取得检测器的灵敏度数据。
(2)偏移修正处理:在未照射X射线的状态下进行测量,并取得检测器的偏移数据。设置被摄体进行旋转测量,取得与旋转对应的多个被摄体数据。从各被摄体数据中减去偏移数据,得到偏移修正后被摄体数据。
(3)灵敏度修正处理:从灵敏度数据中减去偏移数据,得到偏移修正后灵敏度数据。将偏移修正后被摄体数据除以偏移修正后灵敏度数据,得到灵敏度修正后被摄体数据。
(4)对数变换处理:对灵敏度修正后被摄体数据进行对数变换,并进行-1倍运算,得到投影数据。
(5)不均匀性修正处理:对投影数据乘以修正系数G,得到不均匀性修正后数据。
(6)三维数据取得:对所有的被摄体数据完成上述的处理后,进行重建运算处理,得到三维数据。
图1是表示实施例一的不均匀性修正处理的步骤的流程图。图2是表示在实施例一的不均匀性修正处理中使用的变换表A的图。图3是表示在实施例一的不均匀性修正处理中使用的变换表B的图。
将投影数据的横向的坐标设为u,纵向的坐标设为v。对于被摄体的投影数据As(u,v),指定位置v(步骤101)。对于被摄体的投影数据As(u,v),指定位置u(步骤102)。(※修正图1的u环的位置)对于坐标(u,v)检索变换表A,得到过滤器的X射线吸收量Af(u,v)(步骤103)。对于过滤器X射线吸收量Af(u,v),利用变换式A(后述)计算出等效过滤器厚度Ef(u,v)作为与过滤器X射线吸收量Af(u,v)相当的过滤器的厚度(步骤104)。
对于计算出的等效过滤器厚度Ef(u,v)及被摄体投影数据As(u,v)检索变换表B,得到修正系数G(Ef,As)(步骤105)。(※修正图3的G的标记。G的标记在上次的修正中将Ef作为下标,但因为不能使用下标,所以再次修正为“G(Ef,As)”的形式。在文中及图中,请将G的标记统一为“G(Ef,As)”。)在变换表B中,在没有期望的等效过滤器厚度Ef(u,v)及被摄体投影数据As(u,v)的情况下,由相近的值通过插值运算或者外插运算(外挿演算)求得修正系数G(Ef,As)。对于被摄体投影数据As(u,v),乘以修正系数G(Ef,As)(步骤106)。对于所有的坐标u执行处理(步骤107)。对于所有的坐标v执行处理(步骤108)。对于所有的被摄体投影数据进行处理(步骤109)。
说明图2所示的变换表A的求得方法。设置测量被摄体时使用的过滤器410,未设置被摄体,照射X射线,取得测量数据。对于测量数据,执行所述的偏移修正、灵敏度修正、对数变换处理,从而得到投影数据。将在坐标(u,v)上的投影数据的值设为过滤器的X射线吸收量Af(u,v)。使投影数据的横向的坐标u从1变化到最大值Nu,使纵向的坐标v从1变化到最大值Nv,求得各坐标(u,v)上的Af(u,v),存储在表中。
在这些处理中,如果使用将在相同条件下取得的投影数据进行加法平均后的投影数据,则能够使统计噪声减少,使精度提高。
在变换表A的坐标中取多点,通过插值运算可计算出任意坐标上的X射线吸收量。在此情况下,能够降低表的存储容量。
变换表A也可作为以坐标为变量而计算出X射线吸收量的变换式保持,在此情况下,能够降低表的存储容量。
说明变换式A的求得方法。设置规定的同样厚度Ef的过滤器410,未设置被摄体,照射X射线,取得测量数据。对于测量数据执行所述的偏移修正、灵敏度修正、对数变换处理,从而取得投影数据。在投影数据中计算出平均值,并将其设为过滤器410的X射线吸收量Af。作为过滤器410,准备厚度不同的多个过滤器,对于各个厚度的过滤器,取得投影数据,得到对于各个过滤器厚度Ef的X射线吸收量Af。对于X射线吸收量Af的过滤器厚度Ef以式(1)近似,并求得系数a0、a1、a2…。
[式1]
Ef=a0+a1·Af+a2·Af 2+a3·Af 3+… (1)
说明图3所示的变换表B的求得方法。设置规定的同样厚度Ef的过滤器410,且设置任意的被摄体,照射X射线,取得测量数据。此处,作为任意的被摄体,使用在圆柱形容器或椭圆柱容器中装入水的圆柱或丙烯圆柱等。对于测量数据执行所述的偏移修正、灵敏度修正、对数变换处理,从而取得投影数据As(u,v)。设定重建图像的值均匀的理想的被摄体的投影数据Ao(u,v)。在各坐标(u,v)上,将理想投影数据Ao(u,v)除以被摄体投影数据As(u,v),计算出修正系数G(Ef,As)。作为过滤器410,准备厚度不同的多个过滤器,对于各个厚度Ef的过滤器使用任意的被摄体,取得投影数据,并求得各个过滤器厚度Ef的对于被摄体投影数据As的修正系数G(Ef,As)。对于被摄体投影数据As的修正系数G(Ef,As)以式(2)近似,并求得系数b0、b1、b2…。
[式2]
G=b0+b1·As+b2·As 2+b3·As 3+… (2)
对于多个厚度的过滤器,使过滤器厚度Ef从0变化到最大值Emax,使被摄体投影数据As从0变化到最大值Smax,并使用近似式(2)求得G(Ef,As),并存储在表中。
在这些处理中,通过使用将在相同条件下取得的投影数据进行加法平均后的投影数据,能够使统计噪声减少,使修正系数的算出精度提高。
在变换表B中,将过滤器厚度及被摄体投影数据取多点,通过插值运算可计算出修正系数。在此情况下,能够降低表的存储容量。
变换表B也可作为以过滤器厚度及被摄体投影数据为变量而计算出修正系数的变换式保持,在此情况下,能够降低表的存储容量。
图5(A)-(C)是表示将过滤器410沿X射线源401的旋转面切开的剖面(左侧)和立体图(右侧)的三种例子的图,(A)是整体为同样厚度的过滤器410的剖面图和立体图,(B)是由凹面形成的区域和同样厚度的区域组合而成的过滤器410的剖面图和立体图,(C)是由凹面形成的区域、由凸面形成的区域、和同样厚度的区域组合而成的过滤器410的剖面图和立体图。在剖面图中表示箭头并标记X射线是为了表示X射线的入射方向。也可将X射线的入射方向设为箭头的相反方向。如图5(A)-(C)所示,如果过滤器形状为左右对称形,则相对于过滤器中心,过滤器的X射线吸收量左右对称,且所述的吸收量的计算运算量及表A的存储容量能够减半。如果将过滤器的形状设为沿剖面的深度方向同形,则过滤器的X射线吸收量在任一个剖面都相同,从而能够减少吸收量的计算运算量及表A的存储容量。
图6(A)-(C)是表示过滤器410的其他形状的剖面图及俯视图。图6(B)是从X射线的入射方向观察的俯视图,图6(A)是在A-A位置沿箭头方向观察的剖面图,图6(C)是在B-B位置沿箭头观察的剖面图。如图6(A)所示,在A-A位置由凹面形成的区域与由一定值形成的区域的组合形成过滤器410。此外,在B-B位置,过滤器的厚度在上部变化。
在从X射线入射方向观察过滤器的形状的情况下,与图6(B)同样,在A-A剖面中,与图5(C)同样地,由凹面形成的区域、由凸面形成的区域、和由同样厚度形成的区域的组合形成过滤器的形状,在B-B剖面中,在过滤器上部厚度变化并且由凹面形成的区域与由同样厚度形成的区域的边界可设为凸面形状。
如图5(A)所示的过滤器410在被摄体在测量区域(检测器402的检测区域)中比较平坦,平均地防止检测器402的饱和之类的测量的情况下是有用的,如图5(B)、(C)所示的过滤器410在被摄体的厚度在测量区域中两侧急剧减少的测量的情况下是有用的。此外,图6的例子在即使被摄体的厚度在测量区域的上部也急剧减少的测量的情况下是有用的。
具体地说,如图5(C)所示的过滤器410如下构成,即,过滤器的剖面形状在凹形的圆弧的附近具有凸形的圆弧,并在其附近具有直线,凹形的圆弧与凸形的圆弧的交点的切线具有相同斜率,且凸形的圆弧与直线的交点的切线具有相同斜率,由此能够得到相对于头部、腹部、下肢等的人体形状良好的三维图像。此外,根据图6的例子,过滤器的厚度沿旋转装置的旋转中心轴向变化,由此,能够得到相对于头部等具有在体轴方向上厚度不同的形状的检查对象良好的三维图像。
根据实施例一,在设置规定的同样厚度的过滤器并未设置检查对象的状态下,对取得的投影数据计算出平均值,并对多个厚度的过滤器取得平均值,对于过滤器厚度以多项式来近似平均值,并作成表示过滤器的X射线吸收量与过滤器厚度的关系的变换式,由此,能对具有任意的厚度的过滤器进行修正处理。
此外,根据实施例一,在设置规定的同样厚度的过滤器并设置了任意的被摄体的状态下得到投影数据,并求得重建图像的值均匀的理想的检查对象的投影数据,将理想投影数据除以检查对象的投影数据来计算出修正系数,对于多个厚度的过滤器取得修正系数,并取得表示过滤器的厚度、检查对象的投影数据和修正系数的关系的变换表,由此能对具有任意形状的检查对象进行修正处理。
此外,如果使上述的实施例一发展,在代替同样厚度的过滤器而设置具有任意形状的过滤器且未设置检查对象的状态下求得投影数据,并使投影数据的横向的坐标和纵向的坐标变化而求得各坐标的投影数据的值,作成表示检查对象的投影数据和过滤器的X射线吸收量的关系的变换表,由此,能对具有该形状的过滤器进行修正处理。
此外,如果使上述的实施例一发展,在代替同样厚度的过滤器而设置具有任意形状的过滤器且设置了任意的被摄体的状态下求得投影数据,并求得重建图像的值均匀的理想的检查对象的投影数据,并将理想投影数据除以检查对象的投影数据,计算出修正系数,相对于多个厚度的过滤器取得修正系数,作成表示过滤器的厚度、检查对象的投影数据和修正系数的关系的表,由此,能通过具有该形状的过滤器对具有任意形状的检查对象进行修正处理。
此外,根据实施例一,作为任意的被摄体,通过使用水圆柱或者水椭圆柱,能够简便且精度良好地模拟检查对象,并求得修正系数。
(其他实施例)
在上述实施例一的修正处理中,对逐一准备变换表及变换式的情况进行了说明。在其他的实施例中,与X射线管电压、散射X射线屏蔽栅极、准直仪等摄影条件对应,考虑预先准备多个变换表及变换式。在该情况下,控制处理装置406存储这多个变换表及变换式的组,并根据摄影条件进行选择,用于修正处理。通过使用与摄影条件对应的变换表及变换式,能够提高修正精度。此外,在其他实施例中,考虑从存储的多个变换表及变换式的组中,根据摄影条件选择多个组,通过其插值运算或者外插运算来计算出适合于任意摄影条件的修正系数,且用于修正处理。由此,能够对任意的摄影条件实现精度高的修正。
工业上的可利用性
根据本发明,能够提供一种X射线测量装置,其能够得到可消除X射线的检测元件的饱和现象、由过滤器或被摄体导致的能量分布的变化、或者从过滤器或被摄体发生的散射X射线所产生的影响,能抑制检测器的饱和,值的均匀性提高的良好的三维图像。
Claims (17)
1.一种X射线测量装置,具有:
X射线源,其产生向检查对象照射的X射线;
X射线检测器,其以隔着所述检查对象的方式与该X射线源相对配置,并将透过所述检查对象的X射线作为测量数据进行检测;
过滤器,其设置在所述X射线源与所述检查对象之间,对X射线的透过量进行调整;
保持装置,其保持所述X射线源和所述X射线检测器;
旋转装置,其使所述X射线源及所述X射线检测器在所述检测对象的周围旋转;以及
控制处理装置,其对由所述X射线检测器检测的测量数据进行运算处理,所述测量数据是所述X射线检测器以通过所述旋转装置而旋转的、且相对于所述检查对象的多个角度来进行检测的;
该X射线测量装置的特征在于,
所述控制处理装置对所述测量数据进行对数变换而得到投影数据,并得到与所述得到的投影数据对应的所述过滤器的X射线吸收量,对所述得到的X射线吸收量使用规定的变换式来计算所述过滤器的厚度,并根据所述计算出的过滤器的厚度得到与所述得到的投影数据对应的修正系数,将所述得到的修正系数与所述投影数据相乘,重建运算乘以了所述修正系数的投影数据,从而得到三维图像。
2.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置对于所述投影数据检索第一变换表,得到所述过滤器的X射线吸收量。
3.如权利要求2所述的X射线测量装置,其特征在于,
所述第一变换表如以下制成,即:在检查对象的测量前,在设置所述过滤器但不设置检查对象的状态下,得到测量数据,将该得到的测量数据进行对数变换从而得到投影数据,改变该得到的投影数据的横向坐标和纵向坐标,求得各坐标下的投影数据的值,得到检查对象的投影数据与过滤器的X射线吸收量的关系。
4.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置对于所述投影数据检索第二变换表得到所述修正系数。
5.如权利要求4所述的X射线测量装置,其特征在于,
所述第二变换表如以下制成,即:在设置规定的同样厚度的过滤器且设置了任意的被摄体的状态下,得到测量数据,将该得到的测量数据进行对数变换,从而得到被摄体的投影数据,求得重建成像的值均一的理想的检查对象的投影数据,将理想的检查对象的投影数据除以被摄体的投影数据,计算出修正系数,并对多个厚度的过滤器取得修正系数,从而求得过滤器的厚度、被摄体的投影数据、及修正系数的关系。
6.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的X射线吸收量时,使用将在同样条件下取得的投影数据相加平均后的投影数据。
7.如权利要求3所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的X射线吸收量时,在所述第一变换表的坐标中取多点,并通过插值运算计算出任意坐标下的X射线吸收量。
8.如权利要求3所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的X射线吸收量时,将所述第一变换表保持为以坐标为变量计算出X射线吸收量的变换式,并通过所述变换式计算出X射线吸收量。
9.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的厚度时,使用将在同样条件下取得的投影数据相加平均后的投影数据。
10.如权利要求5所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的厚度时,在所述第二变换表中将过滤器厚度及被摄体投影数据取多点,并通过插值运算计算出修正系数。
11.如权利要求5所述的X射线测量装置,其特征在于,
所述控制处理装置在获得所述过滤器的厚度时,将所述第二变换表保持为以过滤器厚度及被摄体投影数据为变量计算出修正系数的变换式,并通过所述变换式计算出修正系数。
12.如权利要求1所述的X射线测量装置,其特征在于,
所述过滤器的形状为左右对称形。
13.如权利要求1所述的X射线测量装置,其特征在于,
所述过滤器的剖面的深度方向的形状相同。
14.如权利要求1所述的X射线测量装置,其特征在于,
所述过滤器剖面的形状构成为:凹形的圆弧邻接凸形的圆弧,且凸形的圆弧邻接直线,凹形圆弧与凸形圆弧的交点切线的斜率相同,凸形圆弧与直线的交点切线的斜率相同。
15.如权利要求1所述的X射线测量装置,其特征在于,
所述过滤器的厚度沿所述旋转装置的旋转中心轴方向有所变化。
16.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置具备从水圆柱或水椭圆柱形状的任意被摄体的投影数据中得到的修正系数。
17.如权利要求1所述的X射线测量装置,其特征在于,
所述控制处理装置根据用于拍摄所述检查对象的摄影条件来从多个修正系数中选择修正系数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP019285/2005 | 2005-01-27 | ||
JP2005019285 | 2005-01-27 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101098660A CN101098660A (zh) | 2008-01-02 |
CN100579453C true CN100579453C (zh) | 2010-01-13 |
Family
ID=36740178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200580046148A Expired - Fee Related CN100579453C (zh) | 2005-01-27 | 2005-12-08 | X射线测量装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7515689B2 (zh) |
JP (1) | JP4599392B2 (zh) |
CN (1) | CN100579453C (zh) |
WO (1) | WO2006080144A1 (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10201321B4 (de) * | 2002-01-15 | 2011-02-24 | Siemens Ag | Computertomographie-Gerät und Verfahren mit aktiver Anpassung der Mess-Elektronik |
JP4512187B2 (ja) * | 2004-07-06 | 2010-07-28 | 株式会社日立メディコ | X線計測装置 |
EP2149284B1 (en) * | 2007-05-31 | 2014-08-13 | General Electric Company | Methods and systems to facilitate correcting gain fluctuations in image |
DE102009012631B4 (de) * | 2009-03-11 | 2011-07-28 | Bayer Schering Pharma Aktiengesellschaft, 13353 | Filter für einen Computertomographen sowie Computertomograph |
CN101587052B (zh) * | 2009-06-15 | 2011-03-30 | 浙江大学 | 基于x射线的密度、浓度和厚度测试装置及方法 |
WO2012093440A1 (ja) | 2011-01-07 | 2012-07-12 | 株式会社 東芝 | X線ct装置 |
US10667767B2 (en) * | 2014-05-02 | 2020-06-02 | General Electric Company | Systems and methods for selecting bowtie filter configuration |
JP6165695B2 (ja) * | 2014-09-24 | 2017-07-19 | 富士フイルム株式会社 | 放射線画像解析装置および方法並びにプログラム |
CN104475462B (zh) * | 2014-11-24 | 2016-08-24 | 北京金自天正智能控制股份有限公司 | 一种x射线测厚仪的在线校正装置及方法 |
CN105433973B (zh) * | 2015-12-30 | 2018-09-18 | 沈阳东软医疗系统有限公司 | Ct扫描设备、ct系统和控制过滤器组件的方法及装置 |
CN106955118B (zh) * | 2017-04-11 | 2020-04-17 | 沈阳开普医疗影像技术有限公司 | Ct形状过滤器缺陷的单能矫正方法 |
CN107179650A (zh) * | 2017-05-12 | 2017-09-19 | 太仓诚泽网络科技有限公司 | 一种基于射线的摄影方法 |
CN110353712B (zh) * | 2019-07-19 | 2023-07-21 | 深圳安科高技术股份有限公司 | 一种滤波器、调节ct系统射线束的方法及装置 |
JP7249666B2 (ja) | 2020-10-30 | 2023-03-31 | 株式会社リガク | 蛍光x線分析装置 |
KR102530973B1 (ko) * | 2021-06-25 | 2023-05-16 | (주)아이스트 | 2차 전지 x-선 검사를 위한 보정 데이터 취득 장치 및 방법 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE347859B (zh) * | 1970-11-30 | 1972-08-14 | Medinova Ab | |
US6418193B1 (en) * | 1999-11-01 | 2002-07-09 | General Electric Company | Imaging system including radiation filter for x-ray imaging |
DE10035984C1 (de) * | 2000-07-24 | 2002-01-31 | Siemens Ag | Röntgen-Computertomographieeinrichtung |
JP2002102217A (ja) * | 2000-09-28 | 2002-04-09 | Ge Medical Systems Global Technology Co Llc | X線ctシステム、ガントリ装置、コンソール端末及びその制御方法及び記憶媒体 |
DE10048775B4 (de) * | 2000-09-29 | 2006-02-02 | Siemens Ag | Röntgen-Computertomographieeinrichtung |
JP2003339686A (ja) * | 2002-05-29 | 2003-12-02 | Shimadzu Corp | X線撮影装置 |
JP4142482B2 (ja) * | 2003-04-04 | 2008-09-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US6990171B2 (en) * | 2003-10-27 | 2006-01-24 | General Electric Company | System and method of determining a user-defined region-of-interest of an imaging subject for x-ray flux management control |
-
2005
- 2005-12-08 JP JP2007500431A patent/JP4599392B2/ja not_active Expired - Fee Related
- 2005-12-08 WO PCT/JP2005/022522 patent/WO2006080144A1/ja active Application Filing
- 2005-12-08 US US11/722,802 patent/US7515689B2/en not_active Expired - Fee Related
- 2005-12-08 CN CN200580046148A patent/CN100579453C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US7515689B2 (en) | 2009-04-07 |
JP4599392B2 (ja) | 2010-12-15 |
CN101098660A (zh) | 2008-01-02 |
JPWO2006080144A1 (ja) | 2008-06-19 |
WO2006080144A1 (ja) | 2006-08-03 |
US20080043900A1 (en) | 2008-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100579453C (zh) | X射线测量装置 | |
Cho et al. | Accurate technique for complete geometric calibration of cone‐beam computed tomography systems | |
CN101726503B (zh) | 用于x射线相衬层析成像的系统和方法 | |
US7302031B2 (en) | Method and arrangement relating to X-ray imaging | |
CN101262819B (zh) | Ct散射的直接测量及校正 | |
CN101578535A (zh) | 能谱重建 | |
CN101416073A (zh) | 用于重建图像的双能量衰减数据的信噪比的动态优化 | |
CN103339497B (zh) | 探测值处理设备 | |
CN101772324A (zh) | X射线ct装置 | |
CN102272860B (zh) | X射线成像装置和x射线成像方法 | |
CN102451012A (zh) | 放射线照相装置和放射线照相系统 | |
Campbell et al. | A prototype fan‐beam optical CT scanner for 3D dosimetry | |
US6560308B1 (en) | Method and system for approximating missing data in cone beam x-ray CT reconstruction | |
CN102160085A (zh) | 包括针对散射辐射的校正单元的成像装置 | |
CN105424731B (zh) | 一种锥束ct的分辨率性能测量装置及标定方法 | |
US20140003573A1 (en) | X-Ray CT System for Measuring Three Dimensional Shapes and Measuring Method of Three Dimensional Shapes by X-Ray CT System | |
Ferrucci et al. | Measurement of the X-ray computed tomography instrument geometry by minimization of reprojection errors—Implementation on experimental data | |
Yang et al. | Design and characterization of high energy micro-CT with a laser-based X-ray source | |
Berggren et al. | Characterization of photon‐counting multislit breast tomosynthesis | |
Dekker et al. | Optical CT imaging of solid radiochromic dosimeters in mismatched refractive index solutions using a scanning laser and large area detector | |
CN104132950B (zh) | 基于原始投影信息的cl扫描装置投影旋转中心标定方法 | |
Rawlik et al. | Refraction beats attenuation in breast CT | |
CN103505231B (zh) | 调整计算机断层扫描系统的半扇形束角度的方法及其系统 | |
Campbell et al. | Radiation‐induced refraction artifacts in the optical CT readout of polymer gel dosimeters | |
Ghani et al. | An empirical method for geometric calibration of a photon counting detector-based cone beam CT system |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100113 Termination date: 20121208 |