CN103327901A - X射线ct装置、物质确定方法及图像处理装置 - Google Patents
X射线ct装置、物质确定方法及图像处理装置 Download PDFInfo
- Publication number
- CN103327901A CN103327901A CN2012800062524A CN201280006252A CN103327901A CN 103327901 A CN103327901 A CN 103327901A CN 2012800062524 A CN2012800062524 A CN 2012800062524A CN 201280006252 A CN201280006252 A CN 201280006252A CN 103327901 A CN103327901 A CN 103327901A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- image
- ray
- projection
- primary standard
- 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
- 238000012545 processing Methods 0.000 title claims description 87
- 238000000034 method Methods 0.000 title claims description 80
- 239000012925 reference material Substances 0.000 claims abstract 8
- 239000000463 material Substances 0.000 claims description 340
- 239000000126 substance Substances 0.000 claims description 119
- 230000000875 corresponding effect Effects 0.000 claims description 64
- 238000009826 distribution Methods 0.000 claims description 45
- 230000009466 transformation Effects 0.000 claims description 26
- 238000010586 diagram Methods 0.000 claims description 23
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 230000002596 correlated effect Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 5
- 239000013077 target material Substances 0.000 abstract 1
- 230000008569 process Effects 0.000 description 43
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 39
- 239000002872 contrast media Substances 0.000 description 36
- 239000000203 mixture Substances 0.000 description 31
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 30
- LEHOTFFKMJEONL-UHFFFAOYSA-N Uric Acid Chemical compound N1C(=O)NC(=O)C2=C1NC(=O)N2 LEHOTFFKMJEONL-UHFFFAOYSA-N 0.000 description 20
- TVWHNULVHGKJHS-UHFFFAOYSA-N Uric acid Natural products N1C(=O)NC(=O)C2NC(=O)NC21 TVWHNULVHGKJHS-UHFFFAOYSA-N 0.000 description 20
- 229940116269 uric acid Drugs 0.000 description 20
- 230000009977 dual effect Effects 0.000 description 17
- 229910000019 calcium carbonate Inorganic materials 0.000 description 15
- 230000000694 effects Effects 0.000 description 13
- 210000004872 soft tissue Anatomy 0.000 description 12
- 210000000845 cartilage Anatomy 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 230000033228 biological regulation Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 8
- 230000009471 action Effects 0.000 description 7
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 6
- 230000005764 inhibitory process Effects 0.000 description 6
- 229910052740 iodine Inorganic materials 0.000 description 6
- 239000011630 iodine Substances 0.000 description 6
- 238000012937 correction Methods 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000003925 fat Substances 0.000 description 3
- 239000012467 final product Substances 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- 239000002904 solvent Substances 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- JYYOBHFYCIDXHH-UHFFFAOYSA-N carbonic acid;hydrate Chemical compound O.OC(O)=O JYYOBHFYCIDXHH-UHFFFAOYSA-N 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000009897 systematic effect 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
-
- 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]
-
- 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
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/20—Sources of radiation
- G01N2223/206—Sources of radiation sources operating at different energy levels
-
- 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
- G01N2223/6126—Specific applications or type of materials biological material tissue
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
实施方式的X射线CT装置根据扫描对象物所获得的投影数据,显示对象物内的图像,具有生成部、变换部、图像形成部和确定部。生成部通过能量不同的X射线分别扫描对象物,生成多个投影数据。变换部将多个投影数据变换为与多个基准物质对应的多个新投影数据。图像形成部通过分别重建由变换部变换的多个新投影数据,形成与多个基准物质对应的多个基准物质图像。确定部根据多个基准物质图像的像素的值的相关,确定对象物质。
Description
技术领域
本发明的实施方式涉及X射线CT装置、物质确定方法及图像处理装置。
背景技术
X射线计算机断层摄影(Computed Tomography:CT)是通过重建用X射线束扫描对象物获得的投影数据,形成表示对象物的信息的图像的技术。
作为X射线CT的应用,有确定对象物所包含的物质的类别的技术。该技术中,难以根据用基于单一的管电压的X射线束获得的图像来判别物质,因此,用基于不同的2个管电压的能量的不同的2个X射线束分别进行扫描的手法在近年备受关注。该手法称为“双能量CT(Dual Energy CT)”。
非专利文献1记载了应用不同的2个管电压来形成2个图像,用这些图像的CT值之比来确定物质的技术。另外,非专利文献1记载的全部事项作为构成本说明书的一部分而引用。
另外,专利文献1记载了通过将应用不同的2个管电压而获得的各投影数据所示的线衰减系数表现为2个基准物质(例如水和骨骼)的线衰减系数的线性组合,来形成表示各基准物质的分布的图像(基准物质图像)的技术。而且,专利文献1还记载了通过将这些基准物质图像进行组合,来形成有效原子序数图像、密度图像及单色X射线图像的手法。另外,专利文献1记载的全部事项作为构成本说明书的一部分而引用。
先行技术文献
专利文献
专利文献1:日本特开2009-261942号公报
非专利文献
非专利文献1:Johnson TR.et al.,“Material differentiation bydual energy CT:initial experience”,Eur Radiol(2007),17,1510-1517
发明内容
但是,如果考虑到CT值依赖于管电压,其结果,CT值之比依赖于管电压的组合,则以往的双能量CT中难高精度地判别CT值比相近的物质。
另外,以往的双能量CT中可以根据基准物质图像来掌握该基准物质的含有量是多还是少,但是难以确定该物质的类别。
本发明解决的课题是提供能够高精度地确定对象物所包含的物质的X射线CT装置、物质确定方法及图像处理装置。
实施方式的X射线CT装置根据扫描对象物而获得的投影数据显示对象物内的图像,该X射线CT装置具有生成部、变换部、图像形成部和确定部。生成部用能量不同的X射线分别扫描对象物,生成多个投影数据。变换部将多个投影数据变换为与多个基准物质对应的多个新的投影数据。图像形成部通过分别重建由变换部变换的多个新投影数据,形成与多个基准物质对应的多个基准物质图像。确定部根据多个基准物质图像的像素的值的相关,来确定对象物质。
附图说明
图1是表示第1实施方式的X射线CT装置的概略构成的一个例子的方框图。
图2是用于说明第1实施方式的X射线CT装置执行的处理的图。
图3是用于说明第1实施方式的X射线CT装置执行的处理的图。
图4是用于说明第1实施方式的X射线CT装置执行的处理的图。
图5是用于说明第1实施方式的X射线CT装置执行的处理的图。
图6是表示第1实施方式的X射线CT装置的动作的一个例子的流程图。
图7是表示第2实施方式的X射线CT装置的概略构成的一个例子的方框图。
图8是表示第2实施方式的X射线CT装置显示的图像的一个例子的概略图。
图9A表示是第2实施方式的X射线CT装置显示的图像的一个例子的概略图。
图9B是表示第2实施方式的X射线CT装置显示的图像的一个例子的概略图。
图9C是表示第2实施方式的X射线CT装置显示的图像的一个例子的概略图。
图10是表示第2实施方式的X射线CT装置显示的图像的一个例子的概略图。
图11是表示第3实施方式的X射线CT装置显示的图像的一个例子的方框图。
图12是用于说明第3实施方式的X射线CT装置执行的处理的概略图。
图13是用于说明第3实施方式的X射线CT装置执行的处理的概略图。
图14是表示第3实施方式的X射线CT装置的动作的一个例子的流程图。
图15A是用于说明X射线CT装置的变形例的概略图。
图15B是用于说明X射线CT装置的变形例的概略图。
图16是表示实施方式的图像处理装置的概略构成的一个例子的方框图。
(符号的说明)
1机架;2预处理部;3重建处理部;4数据处理部;41投影数据变换部;42图像形成部;43物质确定部;5控制部;6显示部;7操作部;8像素值比算出部;9物质确定部;10确定图像形成部;11确定控制部。
具体实施方式
参照图面说明实施方式的X射线CT装置、物质确定方法及图像处理装置。以下的实施方式中,以被检体(患者)作为对象物,但是对象物不限于此。
<第1实施方式>
[构成]
该实施方式的X射线CT装置的概略构成如图1所示。该X射线CT装置具备机架1、预处理部2、重建处理部3、数据处理部4。另外,图示虽然省略,该实施方式的X射线CT装置与一般的装置同样,设置有床、控制台、高电压发生装置等。
〔机架1〕
机架1用于通过X射线扫描被检体。机架1与一般的装置同样,设置有互相相对配置的X射线管及X射线检测器、使它们旋转的旋转机构、滑环机构、倾斜机构、数据收集部(Data Acqusition System:DAS)等。机架1使X射线管及X射线检测器旋转的同时,用X射线扫描被检体。由X射线检测器得到的检测数据通过数据收集部进行收集,发送到预处理部2。
特别地,机架1可以实施用能量不同的2个X射线束分别进行扫描的手法,即双能量CT。X射线的能量依赖于高电压发生装置向X射线管施加的电压(管电压)。作为双能量CT的手法,有Slow-kV Switching方式、Dual Source方式、Fast-kV Switching方式等。所谓Slow-kVSwitching方式是用第1管电压进行了扫描后,用第2管电压进行扫描的方式(双旋转方式)。Dual Source方式是采用具备2个X射线管的机架,向这些X射线管施加不同的管电压而进行扫描的方式(双管方式)。Fast-kV Switching方式是使X射线管及X射线检测器一边旋转,一边对每个视图切换管电压的方式(高速开关方式)。
〔预处理部2〕
预处理部2对从机架1送来的检测数据实施规定的预处理(图像重建处理前进行的处理)。作为该预处理,有计算数据的对数的处理、参考校正、水校正、射束硬化校正、体动校正等。由预处理部2生成的数据称为投影数据。由预处理部2生成的投影数据被发送到重建处理部3、数据处理部4。机架1及预处理部2“生成部”的一个例子而发挥作用。
〔重建处理部3〕
重建处理部3通过对预处理部2生成的投影数据实施重建处理,生成被检体的图像数据。重建处理是从投影数据对被检体的X射线吸收系数的分布进行逆运算的运算处理。作为该运算处理,有2维傅里叶变换法、卷积反投影法、扇束卷积反投影法等。
另外,重建处理部3也可以构成为对通过数据处理部4获得的投影数据实施重建处理而生成图像数据。该处理将随后叙述。
〔数据处理部〕
数据处理部4通过对由预处理部2生成的投影数据实施规定的数据处理,来确定被检体所包含的物质的类别。
在数据处理部4设置有投影数据变换部41、图像形成部42、物质确定部43。另外,投影数据变换部41作为“变换部”的一个例子、图像形成部42作为“图像形成部”的一个例子、物质确定部43作为“确定部”的一个例子而分别发挥作用。
(投影数据变换部41)
投影数据变换部41将通过双能量CT的手法获得的第1及第2投影数据变换为与预定的2个基准物质对应的2个投影数据。作为该处理的一个例子,投影数据变换部41采用专利文献1记载的手法,将第1及第2投影数据分别表现为与2个基准物质对应的预定的2个基准值、以及该2个投影数据所组成的线性组合。即,该线性组合中的2个系数成为目标的2个投影数据。该处理例的详细情况将随后叙述。
(图像形成部42)
图像形成部42通过将由投影数据变换部41获得的2个投影数据分别进行重建,形成与2个基准物质对应的2个基准物质图像。基准物质图像表示对象物内的线衰减系数的分布。线衰减系数表示入射X射线通过单一厚度的物质时衰减的能量的比例。由该重建处理获得的信息是2个基准值的系数和2个基准物质图像组成的线性组合图像。即,该线性组合图像中的2个系数成为目标的2个基准物质图像(重建图像)。
该重建处理以与重建处理部3同样的要领进行。另外,也可以取代图像形成部42,由重建处理部3进行该重建处理。该情况下,不需要图像形成部42。图像形成部42或重建处理部3执行的处理的例子的详细情况将随后叙述。
另外,2个基准值用于将投影数据的线衰减系数表示为线性组合。2个基准值也可以是2个基准物质的线衰减系数。另外,2个基准值也可以是任意的值。
(处理例)
作为上述的处理例,说明专利文献1记载的手法。由预处理部2生成的第1投影数据(高能量侧)用gH表示,第2投影数据(低能量侧)用gL表示。投影数据变换部41通过对这些投影数据gH、gL实施下式(1)所示的变换,生成2个投影数据L1、L2。
[式1]
其中:
D表示下式(2)的右边的2×2矩阵的矩阵式;
<μ>1,2 H,L表示专利文献1记载的能量平均化线衰减系数。
[式2]
图像形成部42通过分别重建式(2)所示2个投影数据L1、L2,形成与2个基准物质对应的2个基准物质图像(下式(3)的c1、c2)。任意的物质的线衰减系数μ如下式(3)所示,表现为采用了2个基准值μ1、μ2和2个基准物质图像c1、c2的线性组合(线性组合图像)。
[式3]
μ(E,x,y)=μ1(E)c1(x,y)+μ2(E)c2(x,y).......(3)
其中:
E表示X射线的能量;
μ1(E)表示能量E中的第1基准物质的线衰减系数(基准值);
μ2(E)表示能量E中的第2基准物质的线衰减系数(基准值);
c1(E,x,y)表示位于坐标(x,y)的像素中的第1基准物质的存在率;
c2(E,x,y)表示位于坐标(x,y)的像素中的第2基准物质的存在率。
基准物质图像(存在率)c1及c2是将任意的物质的线衰减系数μ表示为2个基准物质的线衰减系数μ1及μ2的函数时的系数,是表示该任意的物质与各基准物质的类似程度的指标。
以下,特别以第1基准物质为造影剂(碘浓度50[mgI/ml])、第2基准物质为水的情况为例进行说明。
(物质确定部43)
物质确定部43根据由图像形成部42形成的2个基准物质图像的相关,确定对象物质的类别。另外,所谓对象物质表示该实施方式中作为类别的确定处理的对象的物质。应用上述处理例的情况下,物质确定部43首先对于规定的对象物质,根据2个基准值μ1及μ2来决定预定的坐标系中的与2个基准物质图像c1及c2对应的坐标。而且,物质确定部43对于规定的多个物质,根据预定的上述坐标系中的多个坐标和与由前级的处理决定的2个基准物质图像c1及c2对应的坐标,来确定对象物质的类别。
对于上述坐标系,可将其任意设定。作为一个例子,可应用如图2所示的以2个存在率c1及c2为基的2维坐标系。该坐标系将作为第1基准物质的造影剂的存在率c1设为纵轴,将作为第2基准物质的水的存在率c2设为横轴。
对于该坐标系中的坐标,按纵轴的坐标、横轴的坐标的顺序记载。即,该坐标系中的坐标记载为(c1,c2)。作为纵轴中的基准的坐标P(1,0)、即纵轴方向中的基,应用造影剂的存在率c1与100%相当的矢量。另外,作为横轴中的基准的坐标Q(0,1)、即横轴方向中的基,应用水的存在率c2与100%相当的矢量。
物质确定部43将由图像形成部42获得的线性组合图像中的基准物质图像的组(c1,c2)作为该坐标系的坐标进行处理。
另外,本例中,以基准物质图像的组和坐标成为相同的表现的方式来设定坐标系(基),但是不限于此。例如,可以通过对c1、c2乘以系数而应用使坐标轴伸缩了的坐标系。另外,本例中,向纵轴分配c1,向横轴分配c2,但是也可以反过来。而且,也可以采用斜交坐标系等的其他坐标系来代替本例的正交坐标系。即,该实施方式中的坐标系只要可以表示线性组合图像中的2个基准物质图像的相关即可,其具体的形态是任意的。
供给由物质确定部43进行的处理的上述“多个物质”可以是任意的物质。另外,其个数也任意。图2所示例中,造影剂和水相当于多个物质。本例中,“多个物质”和基的生成中采用的物质相同,但是不限于此。作为其一个例子,将造影剂和水的2个物质作为“基准物质”,将造影剂、水、碳酸钙、脂肪及尿酸的5个物质作为“多个物质”而应用的情况将随后叙述。
与多个物质对应的坐标的取得方法也是任意的。作为其例子,有实际进行该物质的测定而求出坐标的方法、还考虑其他物质的线衰减系数而算出坐标的方法等。
例如,通过用X射线扫描该物质而生成投影数据,将基于此的各个像素的线衰减系数表现为线性组合,决定与该系数的组对应的坐标来进行前者的方法。
以下说明后者的方法。由于各物质的线衰减系数已知,因此通过将不同的2个能量对应的线衰减系数分别代入上述的式(3),获得以下的联立方程式(4a)、(4b)。
[式4]
μ(ELow)=μ1(ELow)c1+μ2(ELow)c2.......(4a)
μ(EHigh)=μ1(Ehigh)c1+μ2(EHigh)c2.....(4b)
其中:
ELow、EHigh表示2种X射线能量;
μ(ELow)表示低侧的能量ELow中的该物质的线衰减系数;
μ(EHigh)表示高侧的能量EHigh中的该物质的线衰减系数。
该联立方程式(4a)、(4b)中未知数是c1和c2的2个,所以通过对其进行求解,可以算出系数的组c1、c2。然后,根据该系数的组,获得目标的坐标。
物质确定部43根据对于这样取得的多个物质的多个坐标、以及对于对象物质由物质确定部43决定的坐标,来确定该对象物质的类别。作为该处理的一个例子,物质确定部43:(1)求出基于该多个坐标的坐标系内的区域,(2)根据对于该区域的对象物质的坐标的位置,来确定对象物质的类别。以下,说明采用连结多个坐标的图形作为坐标系内的区域的情况。这里,除了特别必要的情况,不必实际显示由(1)获得的区域(图形等)。以下的说明中,为了帮助理解实施方式,只不过是采用所显示的图形。以下说明的坐标系内的区域也同样。另外,坐标系内的区域相当于该坐标系中的坐标的集合,图形相当于位于该区域本身、其外周线的坐标的集合。
以下说明求出作为坐标系内的区域的图形的处理。图2所示例中,作为上述多个坐标,应用造影剂的坐标和水的坐标。物质确定部43求出连结这些坐标的线段L1作为目标的图形。
一般,线段的公式通过根据2个坐标来算出通过它们的直线的公式(即斜率和截距)并抽出该直线中以该2个坐标为两端的部分而获得。考虑3个以上的坐标而求出2个以上的线段时,对任意的2个坐标的组合进行同样的运算即可。
现在,线段L1上的坐标与造影剂的碘浓度对应。例如,线段L1上的坐标P1、P2、P3、P4、P5分别与碘浓度25[mgI/ml]、20[mgI/ml]、15[mgI/ml]、10[mgI/ml]、5[mgI/ml]对应。另外,如前所述,坐标P与50[mgI/ml]对应,坐标Q与0[mgI/ml](单纯的水)对应。即,在线段L1上,接近坐标P时碘浓度增加,接近坐标Q时碘浓度减少。
这样,通过采用水作为多个物质的一个,可以将物质的浓度表示为图形。例如,采用造影剂、水、碳酸钙、脂肪及尿酸的5个物质时的图形的例子如图3所示。
图3表示了造影剂的坐标P、水的坐标Q、碳酸钙的坐标R、脂肪的坐标S、尿酸的坐标T。另外,连结坐标P和坐标Q的线段L1相当于造影剂的碘浓度,连结坐标R和坐标Q的线段L2相当于碳酸钙的浓度,连结坐标S和坐标Q的线段L3相当于脂肪的浓度,连结坐标T和坐标Q的线段L4相当于尿酸的浓度。各线段L1~L4中,越接近坐标Q,浓度越低,越远离坐标Q,浓度越高。
另外,从图3可知,4个线段L1~L4仅仅相互在坐标Q交差。这表示除了相当于浓度0(单纯的水)的坐标Q以外,线段L1~L4相互分离,即与4个物质(造影剂、碳酸钙、脂肪及尿酸)对应的坐标在该坐标系中相互分离。换言之,如果该坐标系中的线段的位置不同,则与之对应的物质也不同。
通过利用这样的图形,物质确定部43根据由物质确定部43决定的对象物质的坐标和图形的位置关系,来确定该对象物质的类别。
以将图3所示线段L1~L4用作图形的情况为例进行说明,物质确定部43首先判断对象物质的坐标是否位于这些线段L1~L4中的一个之上。
在判断为对象物质的坐标位于线段Li(i=1~4之一)上的情况下,物质确定部43根据该线段Li中的该坐标的位置,来确定对象物质的类别。该类别不仅包含对象物质的物质名,也包含其浓度(即对象物质和水的成分比)。
另外,也可以构成为仅仅确定物质名。该情况下,仅仅确定对象物质的坐标所位于的线段Li即可。
另一方面,在判断为对象物质的坐标都不位于任一线段L1~L4上的情况下,物质确定部43获得该对象物质与相当于这些线段L1~L4的物质的任一个都不符合的结果。
对3个以上的基准物质获得坐标的情况下,作为坐标系内的区域,可以形成连结这些坐标而成的多角形。由该多角形规定的区域相当于这些基准物质的混合物。另外,所谓该混合物意味着将各基准物质的成分比设为0~100%,全部基准物质的成分比之和设为100%时定义的物质。因此,该混合物中也包含仅将3个以上的基准物质中的一个或2个作为成分的物质。仅将一个基准物质作为成分的物质相当于多角形的顶点,仅将2个基准物质为成分的物质相当于多角形的边。
作为采用3个以上的基准物质时的例子,说明将造影剂、水和碳酸钙的3个作为基准物质的情况。该情况下,如图4所示,获得以相当于造影剂、水、碳酸钙的坐标P、Q、R为顶点的多角形(三角形)A。相当于多角形A的边PQ(线段L1)及QR(线段L2)的区域分别如前所述,表示造影剂的浓度及碳酸钙的浓度。另外,相当于边PR(线段L5)的区域表示造影剂和碳酸钙的混合物中的造影剂和碳酸钙的成分比。边PR中,越接近顶点P,造影剂的成分比越增大,越接近顶点R,碳酸钙的成分比越增大。
另外,多角形A的内部区域,即从多角形A上的区域去除了线段L1、L2、L5的区域,相当于这3个基准物质全部的混合物。对于该内部区域中的坐标,也是越接近顶点P,造影剂的成分比越增大,越接近顶点R,碳酸钙的成分比越增大。而且,越接近相当于水的顶点Q,该混合物的浓度越减少。
物质确定部43判断所决定的对象物质的坐标是否位于多角形上。在判断为该坐标位于多角形的外部的情况下,物质确定部43判断为该对象物质不是该多角形对应的混合物。
另一方面,在判断为该坐标位于多角形上的情况下,物质确定部43判断为该对象物质是该混合物。而且,物质确定部43根据该对象物质的坐标的位置,求出构成该对象物质的3个以上的基准物质的成分比。
另外,以上的例中,未考虑X射线CT装置的测定误差,因此,假定相当于各物质的坐标被唯一决定。该测定误差起因于装置的机械误差(公差等)、噪音的混入等。以下,说明考虑测定误差时的处理。另外,在测定误差足够小等可以容许的情况下,应用上述的处理即可。
说明考虑测定误差时的一个例子。首先,通过对各种物质进行反复测定,取得对于各物质的坐标的分布。该分布信息例如是混入到由该X射线CT装置获得的投影数据的噪音的标准偏差信息。该分布信息存储在物质确定部43。
造影剂、水及碳酸钙的各自的坐标的分布的例子如图5所示。关于造影剂的坐标的分布范围用Pa表示,关于水的坐标的分布范围用Qa表示,关于碳酸钙的坐标的分布范围用Ra表示。
物质确定部43根据标准偏差信息,求出包含与该物质对应的图形的2维区域。以下说明该处理的一个例子。首先,物质确定部43求出连结分布范围Pa和分布范围Qa的2条线段L1a、L1b。作为线段L1a、L1b,例如,采用与分布范围Pa、Qa的双方相接,且相互不交差的线段。从而,获得由分布范围Pa、Qa及线段L1a、L1b包围的区域B。同样,获得由分布范围Ra、Qa及线段L2a、L2b包围的区域C。
区域B用作反映了测定误差的造影剂的浓度的坐标的分布范围。根据线段L1使区域B中不在线段L1上的坐标与浓度的值相对应。作为该相对应的例子,可以求出线段L1上的各位置中与线段L1正交的直线,相等地设定位于该直线上的坐标的浓度。对于表示反映了测定误差的碳酸钙的浓度的分布范围的区域C,也同样设定浓度。
物质确定部43根据所决定的对象物质的坐标和区域B、C的位置关系,来确定该对象物质的类别。例如,在对象物质的坐标位于区域B上的情况下,物质确定部43确定该对象物质是造影剂,而且,根据该坐标和设定在区域B上的坐标的浓度,来确定该对象物质的浓度。
[动作]
说明该实施方式的X射线CT装置的动作。该X射线CT装置的动作例子如图6所示。另外,由于装置各部的动作的详细情况已上述,这里仅进行简单说明。
(S1)
首先,用机架1进行基于双能量CT的摄影。预处理部2将由机架1收集的数据变换为投影数据。从而,生成能量不同的第1及第2投影数据。第1及第2投影数据发送到数据处理部4。
(S2)
投影数据变换部41通过将基于由步骤1生成的各投影数据的各个像素的线衰减系数表现为2个基准物质的线衰减系数的线性组合,进行基于2个基准物质的各投影数据的分解。通过该步骤2,第1及第2投影数据变换为与2个基准物质对应的2个投影数据。
(S3)
而且,图像形成部42通过重建由步骤2获得的线性组合,形成线性组合图像。从而,获得2个基准物质图像。
(S4)
物质确定部43求出由步骤3获得的2个基准物质图像的相关。上述的例子中,对于2个基准物质图像的各像素(x,y),获得规定的2维坐标系中的对应的坐标(c1(x,y),c2(x,y))。
(S5)
物质确定部43根据步骤4获得的相关,将对象物质分解为2个基准物质。上述的例子中,对象物质分解为造影剂和水等。
(S6)
物质确定部43根据步骤5的分解结果,来确定对象物质的类别。该确定结果在例如显示于未图示显示器。另外,该确定结果存储在X射线CT装置的存储装置、网络上的存储装置。
[作用·效果]
以下说明该实施方式的X射线CT装置的作用和效果。
该X射线CT装置通过机架1及预处理部2来生成投影数据。特别地说,通过双能量CT,即、通过用能量不同的第1及第2X射线分别扫描对象物,来生成第1及第2投影数据。
而且,该X射线CT装置具备投影数据变换部41、图像形成部42、物质确定部43。投影数据变换部41将第1及第2投影数据变换为与预定的2个基准物质对应的2个投影数据(新的投影数据)。图像形成部42通过分别重建2个投影数据,形成与2个基准物质对应的2个基准物质图像。物质确定部43根据2个基准物质图像的相关,来确定对象物质的类别。
作为上述的变换处理,投影数据变换部41也可以构成为将第1及第2投影数据分别表现为由与2个基准物质对应的预定的2个基准值和2个投影数据所组成的线性组合。
另外,图像形成部42也可以构成为通过重建该线性组合,形成由2个基准值和2个基准物质图像组成的线性组合图像。
另外,物质确定部43也可以构成为,对于对象物质,决定根据2个基准值预定的坐标系中的与2个基准物质图像对应的坐标,且,根据对于多个物质预先获得的坐标系中的多个坐标和所决定的坐标,来确定该对象物质的类别。
另外,对于第1及第2X射线的各自的能量,2个基准值可以是预定的2个基准物质的线衰减系数。而且,坐标系可以是由基于2个基准物质的线衰减系数的2个基展开的2维坐标系。
另外,物质确定部43也可以构成为,根据对于基于多个坐标的坐标系内的区域的对象物质的坐标的位置,来确定对象物质的类别。
另外,多个物质包含第1及第2物质时,物质确定部43也可以构成为,将连结坐标系中的第1物质的坐标和第2物质的坐标的线段作为该坐标系内的区域而求出,而且,对象物质的坐标位于该线段上时,根据该位置,求出第1物质和第2物质的成分比。
另外,多个物质包含3个以上的物质时,物质确定部43也可以构成为,将连结坐标系中的3个以上的物质的坐标的多角形作为上述坐标系内的区域而求出,对象物质的坐标位于该多角形上时,根据该位置,求出3个以上的物质的成分比。
这里,上述的例子中,考虑3个以上的物质的坐标时的上述坐标系内的区域不限于多角形。该情况下的上述坐标系内的区域一般可以是基于3个以上的物质的坐标的图形,即考虑这些坐标而形成的任意的图形。例如,可以将通过或包含3个以上的物质的坐标的图形用作上述坐标系内的区域。该例子中,图形的外周线不必是直线,另外,物质的坐标也不必处于该图形的外周线上。对象物质的坐标位于图形上的情况下,根据该位置,算出3个以上的物质的成分比。
另外,多个物质中的一个为水时,物质确定部43也可以构成为,将对第1及第2对象物质所决定的第1及第2坐标的每一个与水的坐标连结而成的2个线段作为坐标系内的区域而求出,根据这些线段的位置关系,来确定第1对象物质和第2对象物质是否为同种。
另外,物质确定部43也可以构成为,预先存储混入投影数据的噪音的标准偏差信息,根据该标准偏差信息,求出包含坐标系内的区域的2维区域,根据相对于该2维区域的对象物质的坐标的位置,来确定该对象物质的类别。
根据这样的X射线CT装置,可以将对象物质的特性表现为规定的坐标系中的坐标,根据该坐标和关于多个物质的多个坐标,来确定对象物质的类别。从而,即使是CT值比相近的物质,也可以通过参照坐标来进行更高精度的识别。
另外,通过应用采用了上述坐标系内的区域的构成,可以求出构成对象物质的物质的成分比。而且,通过参照上述坐标系内的区域,可以判别2个以上的对象物质分别是同种还是异种。
另外,通过应用考虑噪音的影响的构成,可以进行更高精度、更高准确率的物质的确定、成分比的确定、对象物质的分离等。
<第2实施方式>
该实施方式的X射线CT装置的概略构成如图7所示。该X射线CT装置包含通过第1实施方式说明了的构成。该X射线CT装置具备机架1、预处理部2、重建处理部3、数据处理部4、控制部5、显示部6、操作部7。数据处理部4设置有投影数据变换部41、图像形成部42、物质确定部43。与第1实施方式相同构成的部分只要没有特别说明,具有与第1实施方式同样的功能。
设为基准物质的一个是水。另外,物质确定部43如在第1实施方式中所说明那样,求出连结对象物质的坐标和水的坐标的线段(参照图2等)。该线段表示对象物质和水的混合物中的成分比,即对象物质的浓度。
数据处理部4通过对对象物(被检体等)的摄影区域的各位置应用第1实施方式的处理,求出摄影区域中的对象物质的分布。该分布信息被输入到重建处理部3。另外,分布信息包含表示摄影区域中的对象物质的存在位置的信息和表示浓度的信息。
接受了分布信息和投影数据的重建处理部3,通过对投影数据中的对象物质的分布区域实施重建处理,形成表示对象物质的分布(例如浓度的分布)的分布图像。作为该处理的一个例子,重建处理部3根据表示分布信息的对象物质的存在位置,来确定与该存在位置对应的像素。然后,重建处理部3根据投影数据,仅仅对所确定的像素重建图像。作为其他处理,重建处理部3根据投影数据来进行通常的重建处理,由此获得的图像中仅仅抽出与上述存在位置对应的像素。通过这些手法获得的重建图像是上述的分布图像。
控制部5接受表示对象物质的浓度的线段上的位置的指定,以与该指定位置相应的形态在显示部6显示分布图像。以下说明该处理。
在该X射线CT装置设置有用于变更对象物质的浓度的用户接口。该用户接口可以根据软件在显示部6显示,也可以是硬件。作为硬件,有拨号盘、滑动条等。采用软件时,例如,可以构成为使显示部6显示模仿拨号盘、滑动杆的图像,通过操作部7对其进行操作。
作为采用软件时的具体例,控制部5根据图2所示图形,使显示部6显示图8所示图像。该图像包含图2所示坐标系、线段图像D、滑块E。在坐标系的坐标轴附有表示物质名的字符串“造影剂”、“水”。线段图像D表示线段L1。滑块E设为可在线段图像D上移动(参照图中的两侧箭头)。滑块E通过例如操作部7的鼠标进行的拖动操作来移动。
用户通过移动滑块E到期望的位置,来进行造影剂的浓度的指定。更具体地说,线段图像D和线段L1预先相对应,控制部5指定与滑块E的位置对应的线段L1的位置,即造影剂的浓度。控制部5根据指定的浓度,变更分布图像的显示形态。作为显示形态的变更处理,有变更像素值(辉度、颜色等)的处理和变更图案的处理等。
存在多个分布图像时,可对各分布图像进行浓度的指定。该情况下,用户采用例如定位装置来指定期望的分布图像。控制部5在所指定的分布图像中与规定的浓度对应的位置显示滑块E。用户使滑块E移动到期望的位置。控制部5根据移动后的滑块E的位置,来变更该分布图像的显示形态。
作为浓度的变更处理,也可以采用输入数值的方法。该情况下,控制部5使显示部6显示用于输入对象物质的浓度信息的输入区域。用户采用操作部7(例如键盘),向输入区域输入期望的浓度的值。输入区域也可以是例如下拉菜单那样,可以从浓度的多个选择项选择期望的选择项。
浓度比初始值提高时的分布图像有时也被称作强调图像。另外,浓度比初始值降低时的分布图像有时也被称作抑制图像。
说明强调图像及抑制图像的例子。图9A表示作为60keV的单色X射线图像的分布图像(原图像)G0。另外,图9A(表示图像的以下的图也同样)中,实际的图像和显示辉度被反转。即,实际的图像中,对象物质的浓度越高,显示辉度越高,而图9A中,浓度越高,显示辉度越低。
单色X射线图像用下式(5)定义。
[式5]
其中:
CTnumber表示CT值;
μ表示式(5)所示对象物质的线衰减系数;
μwater表示水的线衰减系数。
在原图像G0描绘了各种浓度的造影剂的分布图像H01、H02、H03、H04、H05。分布图像H01、H02、H03、H04、H05分别相当于浓度25[mgI/ml]、20[mgI/ml]、15[mgI/ml]、10[mgI/ml]、5[mgI/ml]。
图9B表示基于原图像G0的强调图像G1。在强调图像G1描绘了分别与原图像G0中的分布图像H01、H02、H03、H04、H05对应的分布图像H11、H12、H13、H14、H15。各分布图像H1i(i=1~5)表示提高了分布图像H0i所示浓度的状态。
图9C表示基于原图像G0的抑制图像G2。在抑制图像G2描绘了分别与原图像G0中的分布图像H01、H02、H03、H04、H05对应的分布图像H21、H22、H23、H24、H25。各分布图像H2i(i=1~5)表示降低了分布图像H0i所示浓度的状态。
说明这样的强调图像G1及抑制图像G2的生成方法。首先,控制部5确定坐标位于表示造影剂的浓度的线段L1(参照图2)上的像素。接着,控制部5对于所确定的各像素,通过将式(3)所示造影剂的存在率c1设定成1000,将水的存在率c2设定成0,将其代入式(2a)、(2b),来生成强调图像G1。另外,控制部5通过将造影剂的存在率c1设定成0,将水的存在率c2设定成1000,将这些代入式(2a)、(2b),来生成抑制图像G2。控制部5使显示部6显示所生成的强调图像G1、抑制图像G2。
接着,说明浓度的算出方法的例子。如前述,坐标系中的线段表示物质和水的混合物,线段上的位置相当于浓度。图2所示情况下,表示造影剂的浓度的线段L1处于斜率为-1,截距为1的直线上,因此,造影剂的浓度可以下次式算出。
[式6]
图10示意性地示出了将抽出坐标位于表示造影剂的浓度的线段L1上的系数c1(x,y)、c2(x,y)并用式(6)算出的浓度进行彩色显示时的显示形态的一个例子。在分布图像G描绘了5个分布图像H1、H2、H3、H4、H5。各分布图像Hi(i=1~5)用与其浓度相应的颜色显示。浓度和显示颜色的关系由颜色条J表示。另外,对多个物质的浓度进行彩色显示时,也可以构成为对各个物质通过不同的颜色来显示分布图像。
说明该实施方式的X射线CT装置的作用及效果。
该实施方式的X射线CT装置除了第1实施方式的作用,还具有如下所示作用。另外,物质之一是溶剂(例如水)。物质确定部43求出连结对象物质的坐标和溶剂的坐标的线段。重建处理部3通过对投影数据中的对象物质的分布区域实施重建处理,形成表示对象物质的分布的分布图像。控制部5(显示图像生成部)接受线段上的位置的指定,生成与该指定位置相应的形态的分布图像(显示图像),在显示部6显示。
另外,通过控制部5使显示部6显示表示线段的线段图像,而且,接受用操作部7指定该线段图像上的位置的情况,控制部5可以构成为确定与该指定位置对应的线段上的位置,以与该确定位置相应的形态来显示分布图像。
另外,使显示部6显示用于输入物质的浓度信息的输入区域,而且,接受用操作部7向输入区域输入对象物质的浓度信息的情况,控制部5(浓度算出部、显示图像生成部)可以构成为确定与该浓度信息对应的线段上的位置,生成并显示与该确定位置相应的形态的分布图像。
另外,控制部5可构成为根据线段上的指定位置来变更分布图像的显示颜色。另外,控制部5可构成为根据由物质确定部43所确定的对象物质的类别来变更分布图像的显示颜色。而且,控制部5也可以构成为根据对象物质的类别和线段上的指定位置的双方来变更分布图像的显示颜色。
根据这样的X射线CT装置,除了第1实施方式的效果,还可以形成并显示强调图像、抑制图像。另外,由于可以在视觉上表现物质的浓度,所以可进行浓度的直观理解。
另外,上述的例子中,说明了求出连结对象物质的坐标和溶剂的坐标的线段而形成对象物质的分布图像的情况,但是也可以应用溶剂以外的物质而进行同样的处理。例如,通过求出连结前述的“多个物质”中的一个物质的坐标和对象物质的坐标的线段,而且,对投影数据中的对象物质的分布区域实施重建处理,来形成分布图像。该分布图像表示对象物质及一个物质的分布。即,获得表现2个物质的分布的单一的分布图像。
<第3实施方式>
第1实施方式说明的方法中,连结第1物质的坐标和水的坐标的第1线段与连结第2物质的坐标和水的坐标的第2线段相近时,难以判别这些物质。另一方面,非专利文献1记载的方法中,第1物质的CT值和第2物质的CT值相近时,难以判别这些物质。
因而,该实施方式中,通过兼用这2个方法,提供能够更可靠地进行物质判别的技术。另外,该实施方式中,利用专利文献1记载的技术,还说明形成有效原子序数图像、密度图像及单色X射线图像的处理。
图11中示出了该实施方式的X射线CT装置的概略构成。该X射线CT装置包含通过第1实施方式说明了的构成。另外,虽然图示省略,但是该X射线CT装置还可以包含第2实施方式中说明了的构成。
该X射线CT装置具备机架1、预处理部2、重建处理部3、数据处理部4、像素值比算出部8、物质确定部9、确定图像形成部10、确定控制部11。数据处理部4中设置有投影数据变换部41、图像形成部42、物质确定部43。另外,物质确定部9和确定控制部11分别作为“确定部”的一个例子而发挥作用。
与第1实施方式相同构成的部分只要没有特别说明,就具有与第1实施方式同样的功能。另外,第1实施方式说明的物质确定处理由数据处理部4执行,非专利文献1所述的物质确定处理由像素值比算出部8及物质确定部9执行。
重建处理部3通过对由双能量CT获得的第1及第2投影数据分别实施重建处理,形成第1及第2图像。
像素值比算出部8如非专利文献1所述,算出第1及第2图像中的像素值(CT值)之比。该处理在第1图像和第2图像之间进行像素的对应,算出所对应的2个像素的CT值的比。
物质确定部9通过采用非专利文献1所述的方法,根据由像素值比算出部8所算出的CT值比来确定对象物质的类别。
以下,作为具体例子,说明尿酸、软骨及软组织混合状况下确定尿酸的情况。对这些物质获得的CT值比的分布的例子如图12所示。本例中,设为用管电压80kV及135kV进行基于双能量CT的摄影。图12中,在x轴取管电压80kV的CT值,在y轴取管电压135kV的CT值。从而,通过原点和任意的坐标的直线的斜率相当于CT值比。
分布区域K1、K2、K3分别表示尿酸、软骨、软组织的分布。从分布区域K1、K2、K3的位置关系可知,尿酸和软骨的CT值比不同,可进行相互的判别,而尿酸和软组织中,CT值比相近,所以难以判别。另外,直线M1、M2分别表示尿酸的CT值比的最大值和最小值。
另一方面,图13表示通过第1实施方式说明了的方法所获得的尿酸、软骨骼及软组织的坐标的分布的例子。另外,纵轴及横轴分别根据造影剂及水来设定。分布区域N1、N2、N3分别表示尿酸、软骨、软组织的分布。直线L6是表示尿酸的坐标的分布状态的直线。另外,直线L6a、L6b分别表示与尿酸的分布区域的广度相应的直线的分布范围的上限及下限(参照第1实施方式)。对于软骨及软组织,虽然省略了直线的图示,但是从其形态可以容易地掌握直线的朝向。从分布区域N1、N2、N3的位置可知,本例中,尿酸和软骨中,由于获得相近的直线,因此难以进行物质相互的判别,而尿酸和软组织中,直线的朝向足够不同,可进行判别。
以上的考察下,说明兼用第1实施方式中说明了的方法和非专利文献1所述的方法的处理的一个例子(参照图14)。首先,与第1实施方式同样,通过双能量CT取得投影数据(S21)。
重建处理部3、像素值比算出部8及物质确定部9按前述的要领,获得图12所示CT值比的分布信息,即分布区域K1、K2、K3(S22)。另外,数据处理部4按前述的要领,获得图13所示坐标的分布信息,即分布区域N1、N2、N3(S23)。这些分布信息被发送到确定控制部11。以下,将CT值比的分布信息称为第1分布信息,将坐标的分布信息称为第2分布信息。
确定控制部11根据第1及第2分布信息,分别确定3个物质。作为其处理的例子,确定控制部11首先求出图12所示直线M1、M2和图13所示直线L6a、L6b。另外,直线M1、M2、直线L6a、L6b可以根据由本次的测定所获得的数据作成,也可以是事先获得是测定数据、规定的值(理论值、标准值等)。另外,确定控制部11对于各物质,也可以求出基于该CT值比的分布区域的直线,或者基于该坐标的分布区域的直线。这样,获得由直线M1、M2夹着的确定区域和由直线L6a、L6b夹着的确定区域(S24,S25)。
接着,确定控制部11判断软骨及软组织的各自的CT值比的分布区域K2、K3是否包含在由直线M1、M2夹着的确定区域(S26)。该处理可以是判断分布区域K2、K3本身是否包含在该确定区域,在上述求出直线的情况下,也可以判断该直线是否包含在该确定区域。
另外,为前者的情况下,可以判断分布区域K2、K3的整体是否包含在该确定区域,也可以判断分布区域K2、K3的规定比例是否包含在该确定区域。
另外,为后者的情况下,考虑即使因直线的长度的不同导致稍微的斜率的不同,直线也从该确定区域伸出的情况。该情况下,也可以预先决定直线的长度的范围,在该范围内直线包含在该确定区域时,判断为该直线包含在该确定区域。
接着,确定控制部11判断软骨及软组织的各自的坐标的分布区域N2、N3是否包含在由直线L6a、L6b夹着的确定区域(S27)。该处理可以与CT值比的情况同样地执行。
接着,确定控制部11根据CT值比的分布区域K2、K3是否包含在确定区域的判断结果和坐标的分布区域N2、N3是否包含在确定区域的判断结果,进行尿酸、软骨及软组织的判别。
作为其具体例子,确定控制部11判断为CT值比的分布区域包含在确定区域(S28:是),且,坐标的分布区域包含在确定区域(S29:是)时,判断为与这些分布区域对应的物质是尿酸(S30)。
另外,确定控制部11判断为CT值比的分布区域包含在确定区域(S28:是),且,坐标的分布区域不包含在确定区域(S29:否)时,判断为与这些分布区域对应的物质是软组织(S31)。
另外,确定控制部11判断为CT值比的分布区域不包含在确定区域(S28:否),且,坐标的分布区域包含在确定区域(S32:是)时,判断为与这些分布区域对应的物质是软骨(S33)。
另外,确定控制部11判断为CT值比的分布区域不包含在确定区域(S28:否),且,坐标的分布区域不包含在确定区域(S32:否)时,判断为与这些分布区域对应的物质不是尿酸,不是软骨,也不是软组织。以上,该处理例的说明结束。
接着,说明确定图像形成部10。确定图像形成部10形成有效原子序数图像、密度图像及单色X射线图像中的至少一个(也称为确定图像)。通过用专利文献1所述的方法组合根据第1实施方式中说明了的线性组合(式(3))中的系数c1及c2而形成的、表示各基准物质的分布的图像(基准物质图像)来执行该处理。
另外,有效原子序数图像是表示对象物中的有效原子序数的分布的图像。另外,密度图像是表示对象物质中的物质的密度的分布的图像。另外,单色X射线图像是模拟地再现通过具有单一能量的X射线扫描对象物质时的图像。能够形成这些图像的确定图像形成部10作为“形成部”的一个例子而发挥功能。
该实施方式的物质确定部43根据通过第1实施方式中说明了的处理所获得的多个对象物质的各自的确定结果和对这些对象物质分别由确定图像形成部10所形成的有效原子序数图像、密度图像及/或单色X射线图像,可以分别确定这些对象物质的类别。另外,第1实施方式中说明的处理根据对多个物质预先获得的坐标系中的多个坐标和对于对象物质由物质确定部43所决定的坐标,来确定对象物质的类别。
根据该实施方式的X射线CT装置,可以获得第1实施方式的效果,并且,根据与第1实施方式同样的坐标系中的各物质的坐标的分布状态和各物质的CT值比的分布状态的双方,可以确定物质的类别。从而,可以进一步提高对象物所包含的物质的确定精度。
另外,根据该实施方式的X射线CT装置,可以形成有效原子序数图像、密度图像及/或单色X射线图像,因此,可以在视觉上掌握物质的分布状态。另外,通过兼用采用了坐标系的物质确定处理和采用了这些图像的物质确定处理,可以进一步提高对象物所包含的物质的确定精度。
<变形例>
说明上述实施方式的变形例。
说明第1变形例。对象物质中包含未知的成分时,有时得不到图5所示的线段Li那样的直线的相关关系,难以适当决定区域B等。例如图15A所示,设为对于3个对象物质获得3个坐标V1、V2、V3。这些坐标V1、V2、V3不位于同一直线上。该情况下,使显示部显示该坐标系和坐标V1、V2、V3。用户参照显示的坐标系及坐标V1、V2、V3,用操作部在该坐标系内任意设定区域。此时,例如图15B所示,以包含坐标V1、V2、V3的方式来设定区域W。
区域W用作反映了由对象物的未知的成分等所引起的误差的造影剂的浓度的坐标的分布范围。对于区域W中不在连结坐标(c1,c2)=(1,0)和坐标(c1,c2)=(0,1)的线段(相当于图5的线段L1)上的坐标,根据该线段使浓度的值对应。该对应与例如图5所示情况同样地进行。物质确定部43通过与例如图5所示情况同样的手法,根据所决定的对象物质的坐标和区域W的位置关系,来确定该对象物质的类别。
根据该变形例,即使在对象物质中包含未知成分等的误差要因介入的情况下,也可以通过用户的判断高精度地进行物质的确定。
说明第2变形例。上述实施方式中,虽然应用了利用能量不同的2个X射线的双能量CT的手法,但是实施方式不限于此。
例如,能够采用能量不同的3个以上的X射线进行扫描的X射线CT装置中,生成部通过能量不同的N(3以上)的X射线分别扫描对象物,生成N个投影数据。变换部将这N个投影数据变换为与N个基准物质对应的N个新投影数据。图像形成部通过分别重建由变换部所变换的N个新投影数据,形成与N个基准物质对应的N个基准物质图像。确定部根据N个基准物质图像的像素的值的相关来确定对象物质。
该情况下在所执行的处理中,上述实施方式中,2个线性组合的基的个数成为N个,上述实施方式中,2维坐标系的维成为N维。即,在能量不同的X射线的数为2的情况和3以上的情况之间,没有实质的差异。因此,上述实施方式作为采用能量不同的多个(2以上)的X射线的情况,可以如下进行一般化。
实施方式的X射线CT装置根据扫描对象物所获得的投影数据,显示对象物内的图像,具有生成部、变换部、图像形成部和确定部。生成部用能量不同的X射线分别扫描对象物,生成多个投影数据。变换部将多个投影数据变换为与多个基准物质对应的多个新投影数据。图像形成部通过分别重建由变换部变换的多个新投影数据,形成与多个基准物质对应的多个基准物质图像。确定部根据多个基准物质图像的像素的值的相关,确定对象物质。
根据这样的X射线CT装置,与上述实施方式同样,可以高精度地确定对象物所包含的物质。而且,通过采用能量不同的3以上的X射线,与采用2个X射线的上述实施方式相比,可以更高精度进行物质确定。
也可以对上述实施方式说明了的任意的构成和处理进行同样的一般化处理。该情况下,可以获得上述实施方式同样的效果,并进一步提高物质确定的精度。
<物质确定方法>
实施方式的物质确定方法例如由上述实施方式的X射线CT装置来执行。实施方式的物质确定方法包含生成步骤、变换步骤、图像形成步骤、确定步骤。生成步骤中,用能量不同的X射线分别扫描对象物,生成多个投影数据。变换步骤中,将多个投影数据变换为与多个基准物质对应的多个新投影数据。图像形成步骤中,通过分别重建变换步骤中所变换的多个新投影数据,形成与多个基准物质对应的多个基准物质图像。确定步骤中,根据多个基准物质图像的像素的值的相关,来确定对象物质。
确定步骤也可以通过判定多个基准物质图像中的对应像素的值是否具有预先设定的相关,来确定对象物质。
变换步骤也可以采用包含与基准物质对应的X射线衰减系数的计算式,将多个投影数据变换为与多个基准物质对应的多个新投影数据。
还可以包含显示图像生成步骤,该步骤生成能够识别由确定步骤所确定的对象物质的区域的对象物内的图像。
还可以包含显示图像生成步骤,该步骤生成与其他区域的像素相比强调或抑制了由确定步骤所确定的对象物质的区域的像素的显示图像。
还可以包含浓度算出步骤,该步骤根据预先存储的计算式,从多个基准物质图像中的对应像素的值求出对象物质的浓度信息。
还可以包含显示图像生成步骤,该步骤生成表示浓度信息的显示图像。
也可以在图像形成步骤中,生成通过重建多个投影数据而获得的多个能量图像,且,在确定步骤中,根据多个能量图像的多个基准物质图像的像素的值的相关来确定对象物质。
还可以包含形成步骤,该步骤根据基准物质图像,形成表示对象物中的有效原子序数的分布的有效原子序数图像、表示对象物质中的物质的密度的分布的密度图像及模拟地再现用具有单一能量的X射线扫描对象物质的情况的单色X射线图像中的至少一个,且在确定步骤中,根据形成步骤中所形成的图像,来进行对象物质的确定。
根据这样的物质确定方法,与第1实施方式同样,即使是CT值比相近的物质,也可以通过参照坐标,进行高精度的识别。另外,可以求出构成对象物质的物质的成分比,或者可以判别2个以上的对象物质分别是同种还是异种。另外,通过如第1实施方式那样应用考虑了噪音影响的物质确定方法,可以进行更高精度、高准确率的物质的确定、成分比的确定、对象物质的分离等。
另外,根据实施方式的物质确定方法,与第2实施方式同样,除了上述效果,还可以形成并显示强调图像、抑制图像。另外,由于可以在视觉上表现物质的浓度,因此可以进行浓度的直观理解。
另外,根据实施方式的物质确定方法,与第3实施方式同样,除了上述效果,通过根据各物质的坐标的分布状态和各物质的CT值比的分布状态的双方来确定物质的类别,可以进一步提高对象物所包含的物质的确定精度。另外,由于可以形成有效原子序数图像、密度图像及/或单色X射线图像,因此可以在视觉上掌握物质的分布状态。另外,通过兼用采用了坐标系的物质确定处理和采用了这些图像的物质确定处理,可以进一步提高对象物所包含的物质的确定精度。
<图像处理装置>
说明实施方式的图像处理装置。实施方式的图像处理装置的构成例如图16所示。图像处理装置100经由院内LAN等的网络与图像保管装置300连接。图像保管装置300存储由X射线CT装置200形成的图像。
X射线CT装置200具有从例如第1实施方式的X射线CT装置(参照图1)除去了物质确定部43的构成。即,X射线CT装置200是根据扫描对象物所获得的投影数据来显示上述对象物内的图像的装置,具有以下功能:用能量不同的X射线分别扫描对象物而生成多个投影数据的生成功能;将多个投影数据变换为与多个基准物质对应的多个新投影数据的变换功能;通过分别重建由变换功能所变换的多个新投影数据,来形成与多个基准物质对应的多个基准物质图像的图像形成功能。X射线CT装置200将所形成的多个基准物质图像经由院内LAN等的网络送到图像保管装置300。图像保管装置300例如是图像保存通信系统(Picture Archiving and Communication Systems;PACS)。
图像处理装置100处理由X射线CT装置200形成而在图像保管装置300存储的图像。另外,也可以经由其他路由向图像处理装置100输入图像。作为其他路由的例,有存储图像的DVD等的记录介质。
图像处理装置100具有进行从外部取得图像的处理的图像取得部101。取得图像保管装置300存储的图像时,图像取得部101包含经由网络与图像保管装置300进行通信的通信部。而且,图像取得部101例如包含:用于选择取得的图像的用户接口(显示部、操作部等);生成用于取得选择的图像的信号,控制通信部,将该信号向图像保管装置300发送的通信控制部。经由记录介质取得图像时,图像取得部101包含可读取在该记录介质记录的信息的读取部(驱动器装置等)。
由图像取得部101所取得的图像在存储部102存储。存储部102构成为包含硬盘驱动器等的存储装置。该实施方式中,X射线CT装置200按上述的要领形成的图像在存储部102存储。该图像是,通过分别重建根据多个基准物质来变换用能量不同的X射线分别扫描对象物而生成的多个投影数据而获得的多个新投影数据,而形成的与多个基准物质对应的多个基准物质图像。
物质确定部103根据在存储部102存储的多个基准物质图像的像素的值的相关来确定对象物质。该处理与例如第1实施方式的物质确定部43同样地执行。物质确定部103作为“确定部”的一个例子而发挥作用。
图像处理装置100具有显示部时,使显示部显示物质确定部103的处理结果。另外,图像处理装置100为例如网络上的服务器时,图像处理装置100具有将物质确定部103的处理结果经由网络发送给规定的用户终端的功能(发送部)。另外,图像处理装置100也可以具有将物质确定部103的处理结果记录在记录介质的记录部(驱动器装置等)。
根据这样的图像处理装置100,能够高精度地确定对象物所包含的物质。另外,图像处理装置100也可以具有上述实施方式的X射线CT装置的任意的功能。应用以下任一构成时,可实现上述实施方式中说明的该构成对应的作用及效果。
例如,物质确定部103也可以通过判定多个基准物质图像中的对应像素的值是否具有预先设定的相关,来确定对象物质。
另外,图像处理装置100也可以具有显示图像生成部,该显示图像生成部生成可识别由物质确定部103所确定的对象物质的区域的对象物内的图像。
另外,图像处理装置100也可以具有显示图像生成部,该显示图像生成部生成与其他区域的像素相比强调或抑制了由物质确定部103所确定的对象物质的区域的像素的显示图像。
另外,图像处理装置100也可以具有浓度算出部,该浓度算出部根据预先存储的计算式,从多个基准物质图像中的对应像素的值求出对象物质的浓度信息。此外,图像处理装置100也可以具有生成表示浓度信息的显示图像的显示图像生成部。
另外,物质确定部103也可以根据多个能量图像的多个基准物质图像的像素的值的相关,来确定对象物质。
另外,图像处理装置100还可以构成为具有形成部,该形成部根据基准物质图像,形成表示对象物中的有效原子序数的分布的有效原子序数图像、表示对象物质中的物质的密度的分布的密度图像、以及模拟地再现用具有单一能量的X射线扫描对象物质的情况的单色X射线图像中的至少一个,且,物质确定部103根据由形成部形成的图像,进行对象物质的确定。
<程序及存储介质>
可以构成使X射线CT装置或其所包含的计算机执行上述实施方式说明的处理的程序。另外,使图像处理装置(计算机)执行上述实施方式中说明的处理的程序可以在例如DVD等的存储介质存储。另外,可以构筑经由因特网、LAN等的网络来发送这些程序之一的系统。
虽然说明了本发明的实施方式,但是这些实施方式只是例示,而不是意图限定发明的范围。这些新实施方式可以通过各种形态实施,在不脱离发明的要旨的范围,可以进行各种省略、置换、变更。这些实施方式及其变形包含于发明的范围和要旨,并且包含在权利要求的范围记载的发明及其等同的范围。
Claims (19)
1.一种X射线CT装置,根据对对象物进行扫描而获得的投影数据,显示上述对象物内的图像,其特征在于,具备:
生成部,用能量不同的X射线分别对对象物进行扫描,生成多个投影数据;
变换部,将上述多个投影数据变换为与多个基准物质对应的多个新投影数据;
图像形成部,通过分别重建由上述变换部所变换的多个新投影数据,形成与多个基准物质对应的多个基准物质图像;以及
确定部,根据上述多个基准物质图像的像素的值的相关,确定对象物质。
2.如权利要求1所述的X射线CT装置,其特征在于,
上述确定部通过判定上述多个基准物质图像中的对应的像素的值是否具有预先设定的相关,来确定上述对象物质。
3.如权利要求1或2所述的X射线CT装置,其特征在于,
上述变换部采用包含与上述基准物质对应的X射线衰减系数的计算式,将上述多个投影数据变换为与上述多个基准物质对应的多个新投影数据。
4.如权利要求1至3的任一项所述的X射线CT装置,其特征在于,具备:
显示图像生成部,生成能够识别由上述确定部确定的上述对象物质的区域的上述对象物内的图像。
5.如权利要求1至3的任一项所述的X射线CT装置,其特征在于,具备:
显示图像生成部,生成与其他区域的像素相比强调或抑制了由上述确定部确定的上述对象物质的区域的像素的显示图像。
6.如权利要求1至3的任一项所述的X射线CT装置,其特征在于,具备:
浓度算出部,根据预先存储的计算式,从上述多个基准物质图像中的对应的像素的值,求出上述对象物质的浓度信息。
7.如权利要求6所述的X射线CT装置,其特征在于,具备:
显示图像生成部,生成表示上述浓度信息的显示图像。
8.如权利要求1至7的任一项所述的X射线CT装置,其特征在于,
上述图像形成部生成重建上述多个投影数据而获得的多个能量图像,
上述确定部根据上述多个能量图像的上述多个基准物质图像的像素的值的相关来确定对象物质。
9.如权利要求1所述的X射线CT装置,其特征在于,还具备:
形成部,根据上述基准物质图像,形成表示上述对象物中的有效原子序数的分布的有效原子序数图像、表示上述对象物质中的物质的密度的分布的密度图像、以及模拟地再现通过具有单一能量的X射线对上述对象物质进行扫描的情况的单色X射线图像中的至少一个,
上述确定部根据由上述形成部形成的图像,进行上述对象物质的确定。
10.一种物质确定方法,其特征在于,包含:
生成步骤,通过能量不同的X射线分别对对象物进行扫描,生成多个投影数据;
变换步骤,将上述多个投影数据变换为与多个基准物质对应的多个新投影数据;
图像形成步骤,通过分别重建在上述变换步骤中变换了的多个新投影数据,形成与多个基准物质对应的多个基准物质图像;以及
确定步骤,根据上述多个基准物质图像的像素的值的相关,来确定对象物质。
11.如权利要求10的物质确定方法,其特征在于,
在上述确定步骤中,通过判定上述多个基准物质图像中的对应的像素的值是否具有预先设定的相关,来确定上述对象物质。
12.如权利要求10或11所述的物质确定方法,其特征在于,
在上述变换步骤中,采用包含与上述基准物质对应的X射线衰减系数的计算式,将上述多个投影数据变换为与上述多个基准物质对应的多个新投影数据。
13.如权利要求10至12的任一项所述的物质确定方法,其特征在于,包含:
显示图像生成步骤,生成能够识别由上述确定步骤确定的上述对象物质的区域的上述对象物内的图像。
14.如权利要求10至12的任一项所述的物质确定方法,其特征在于,包含:
显示图像生成步骤,生成与其他区域的像素相比强调或抑制了由上述确定步骤确定的上述对象物质的区域的像素的显示图像。
15.如权利要求10至12的任一项所述的物质确定方法,其特征在于,包含:
浓度算出步骤,根据预先存储的计算式,从上述多个基准物质图像中的对应的像素的值,求出上述对象物质的浓度信息。
16.如权利要求15所述的物质确定方法,其特征在于,包含:
显示图像生成步骤,生成表示上述浓度信息的显示图像。
17.如权利要求10至16的任一项所述的物质确定方法,其特征在于,
在上述图像形成步骤中,生成对上述多个投影数据进行重建而获得的多个能量图像,
在上述确定步骤中,根据上述多个能量图像的上述多个基准物质图像的像素的值的相关来确定对象物质。
18.如权利要求10所述的物质确定方法,其特征在于,还包含:
形成步骤,根据上述基准物质图像,形成表示上述对象物中的有效原子序数的分布的有效原子序数图像、表示上述对象物质中的物质的密度的分布的密度图像、以及模拟地再现通过具有单一能量的X射线对上述对象物质进行扫描的情况的单色X射线图像中的至少一个,
在上述确定步骤中,根据上述形成步骤中形成的图像,进行上述对象物质的确定。
19.一种图像处理装置,其特征在于,具备:
存储部,存储通过分别重建多个新投影数据而形成的与多个基准物质对应的多个基准物质图像,其中,根据多个基准物质对通过能量不同的X射线分别扫描对象物而生成的多个投影数据进行变换而获得上述多个新投影数据;以及
确定部,根据上述存储部存储的上述多个基准物质图像的像素的值的相关,来确定对象物质。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011-213047 | 2011-09-28 | ||
JP2011213047 | 2011-09-28 | ||
JP2012206419A JP6073616B2 (ja) | 2011-09-28 | 2012-09-20 | X線ct装置、画像処理装置及びプログラム |
JP2012-206419 | 2012-09-20 | ||
PCT/JP2012/074329 WO2013047403A1 (ja) | 2011-09-28 | 2012-09-24 | X線ct装置、物質特定方法、及び画像処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103327901A true CN103327901A (zh) | 2013-09-25 |
CN103327901B CN103327901B (zh) | 2016-01-27 |
Family
ID=47995436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280006252.4A Active CN103327901B (zh) | 2011-09-28 | 2012-09-24 | X射线ct装置、物质确定方法及图像处理装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9418451B2 (zh) |
JP (1) | JP6073616B2 (zh) |
CN (1) | CN103327901B (zh) |
WO (1) | WO2013047403A1 (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015051463A1 (en) * | 2013-10-09 | 2015-04-16 | Voti Inc. | Techniques for imaging a scanned object |
CN105188545A (zh) * | 2013-05-24 | 2015-12-23 | 株式会社日立医疗器械 | X射线ct装置以及处理方法 |
CN108135560A (zh) * | 2015-10-27 | 2018-06-08 | 株式会社日立制作所 | X射线ct数据处理装置以及搭载其的x射线ct装置 |
CN113361632A (zh) * | 2021-06-25 | 2021-09-07 | 西门子数字医疗科技(上海)有限公司 | 图像中确定生物组织类别方法、装置、计算机设备和介质 |
CN114813798A (zh) * | 2022-05-18 | 2022-07-29 | 中国工程物理研究院化工材料研究所 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
US11885752B2 (en) | 2021-06-30 | 2024-01-30 | Rapiscan Holdings, Inc. | Calibration method and device therefor |
US11977037B2 (en) | 2018-10-22 | 2024-05-07 | Rapiscan Holdings, Inc. | Insert for screening tray |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5798526B2 (ja) * | 2012-07-18 | 2015-10-21 | 株式会社リガク | X線分析装置、x線分析システム、x線分析方法およびx線分析プログラム |
US9978158B2 (en) * | 2013-08-30 | 2018-05-22 | Koninklijke Philips N.V. | Spectral projection data de-noising with anti-correlation filter |
JP6482934B2 (ja) | 2014-06-03 | 2019-03-13 | キヤノンメディカルシステムズ株式会社 | 画像処理装置、放射線検出装置および画像処理方法 |
JP6351164B2 (ja) * | 2014-06-12 | 2018-07-04 | 国立研究開発法人量子科学技術研究開発機構 | ビーム照射対象確認装置、ビーム照射対象確認プログラム、および阻止能比算出プログラム |
US9964499B2 (en) * | 2014-11-04 | 2018-05-08 | Toshiba Medical Systems Corporation | Method of, and apparatus for, material classification in multi-energy image data |
CN104574292B (zh) * | 2014-11-26 | 2018-06-26 | 沈阳东软医疗系统有限公司 | 一种ct图像的校正方法和装置 |
KR101725099B1 (ko) * | 2014-12-05 | 2017-04-26 | 삼성전자주식회사 | 컴퓨터 단층 촬영장치 및 그 제어방법 |
US9913622B2 (en) | 2015-02-23 | 2018-03-13 | Toshiba Medical Systems Corporation | X-ray CT apparatus and image processing device |
CN107427276B (zh) * | 2015-03-30 | 2021-01-26 | 株式会社日立制作所 | 图像生成装置、图像生成方法以及x射线ct装置 |
CN107924081B (zh) * | 2015-08-31 | 2021-12-10 | 富士胶片株式会社 | 液晶单元以及三维结构液晶单元 |
TW201736865A (zh) * | 2016-04-13 | 2017-10-16 | Nihon Medi-Physics Co Ltd | 來自核子醫學影像的生理累積之自動去除及ct影像之自動分段 |
EP3441005B1 (de) * | 2017-08-11 | 2022-04-06 | Siemens Healthcare GmbH | Analyse von läsionen mit hilfe der multi-energie-ct-bildgebung |
CN111201452B (zh) * | 2017-10-09 | 2023-09-12 | 皇家飞利浦有限公司 | 体积图像数据的材料选择性自适应混合 |
WO2020003744A1 (ja) * | 2018-06-27 | 2020-01-02 | キヤノン株式会社 | 放射線撮影装置、放射線撮影方法およびプログラム |
US10718724B2 (en) * | 2018-11-06 | 2020-07-21 | Jon M Frenn | Systems and methods of comparative computed tomography (CT) for qualification of commercial grade items |
JP2020185223A (ja) * | 2019-05-15 | 2020-11-19 | 国立大学法人北海道大学 | 単色ct画像作成方法、単色ct画像作成装置および単色ct画像作成プログラム |
KR20210126647A (ko) * | 2019-05-27 | 2021-10-20 | 다이아트렌드 코포레이션 | 데이터 처리 장치 및 데이터 처리 방법, 및, 그 장치 또는 방법을 탑재 또는 실시하는 x선 시스템 |
WO2021035477A1 (en) * | 2019-08-26 | 2021-03-04 | Guangxi Liugong Machinery Co., Ltd. | Electric excavator |
CN112730468B (zh) * | 2019-10-28 | 2022-07-01 | 同方威视技术股份有限公司 | 物品检测设备及其检测物品的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1509686A (zh) * | 2002-11-27 | 2004-07-07 | GEҽҩϵͳ����Ƽ���˾ | 定量组织脂肪含量的方法和装置 |
CN101028199A (zh) * | 2006-02-28 | 2007-09-05 | 西门子公司 | 借助多能量计算机断层造影确定体内物质浓度的方法和装置 |
JP2009153829A (ja) * | 2007-12-27 | 2009-07-16 | Ge Medical Systems Global Technology Co Llc | 画像処理装置、プログラムおよびx線ct装置 |
US20090262997A1 (en) * | 2008-04-21 | 2009-10-22 | Kabushiki Kaisha Toshiba | Method, apparatus, and computer-readable medium for pre-reconstruction decomposition and calibration in dual energy computed tomography |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5367443B2 (ja) * | 2009-04-28 | 2013-12-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US8705688B2 (en) * | 2010-10-13 | 2014-04-22 | Kabushiki Kaisha Toshiba | Medical image processing apparatus, X-ray computed tomography apparatus, and medical image processing method |
US20140376686A1 (en) * | 2011-02-18 | 2014-12-25 | Smiths Heimann Gmbh | System and method for multi-scanner x-ray inspection |
-
2012
- 2012-09-20 JP JP2012206419A patent/JP6073616B2/ja active Active
- 2012-09-24 US US13/994,346 patent/US9418451B2/en active Active
- 2012-09-24 WO PCT/JP2012/074329 patent/WO2013047403A1/ja active Application Filing
- 2012-09-24 CN CN201280006252.4A patent/CN103327901B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1509686A (zh) * | 2002-11-27 | 2004-07-07 | GEҽҩϵͳ����Ƽ���˾ | 定量组织脂肪含量的方法和装置 |
CN101028199A (zh) * | 2006-02-28 | 2007-09-05 | 西门子公司 | 借助多能量计算机断层造影确定体内物质浓度的方法和装置 |
JP2009153829A (ja) * | 2007-12-27 | 2009-07-16 | Ge Medical Systems Global Technology Co Llc | 画像処理装置、プログラムおよびx線ct装置 |
US20090262997A1 (en) * | 2008-04-21 | 2009-10-22 | Kabushiki Kaisha Toshiba | Method, apparatus, and computer-readable medium for pre-reconstruction decomposition and calibration in dual energy computed tomography |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105188545A (zh) * | 2013-05-24 | 2015-12-23 | 株式会社日立医疗器械 | X射线ct装置以及处理方法 |
CN105188545B (zh) * | 2013-05-24 | 2018-04-06 | 株式会社日立制作所 | X射线ct装置以及处理方法 |
WO2015051463A1 (en) * | 2013-10-09 | 2015-04-16 | Voti Inc. | Techniques for imaging a scanned object |
US10089956B2 (en) | 2013-10-09 | 2018-10-02 | Voti Inc. | Techniques for imaging a scanned object |
CN108135560A (zh) * | 2015-10-27 | 2018-06-08 | 株式会社日立制作所 | X射线ct数据处理装置以及搭载其的x射线ct装置 |
CN108135560B (zh) * | 2015-10-27 | 2021-07-30 | 株式会社日立制作所 | X射线ct数据处理装置以及搭载其的x射线ct装置 |
US11977037B2 (en) | 2018-10-22 | 2024-05-07 | Rapiscan Holdings, Inc. | Insert for screening tray |
CN113361632A (zh) * | 2021-06-25 | 2021-09-07 | 西门子数字医疗科技(上海)有限公司 | 图像中确定生物组织类别方法、装置、计算机设备和介质 |
US11885752B2 (en) | 2021-06-30 | 2024-01-30 | Rapiscan Holdings, Inc. | Calibration method and device therefor |
CN114813798A (zh) * | 2022-05-18 | 2022-07-29 | 中国工程物理研究院化工材料研究所 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
Also Published As
Publication number | Publication date |
---|---|
US9418451B2 (en) | 2016-08-16 |
WO2013047403A1 (ja) | 2013-04-04 |
US20130287260A1 (en) | 2013-10-31 |
CN103327901B (zh) | 2016-01-27 |
JP2013081765A (ja) | 2013-05-09 |
JP6073616B2 (ja) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103327901A (zh) | X射线ct装置、物质确定方法及图像处理装置 | |
US20170024113A1 (en) | Image handling and display in x-ray mammography and tomosynthesis | |
CN1981717B (zh) | 利用医学诊断成像工作站整合三维和二维监视器的方法和装置 | |
EP1685538B1 (en) | Device and method for generating a three-dimensional vascular model | |
CN102125439B (zh) | 图像处理装置、x射线计算机断层摄像装置以及图像处理方法 | |
Lyu et al. | Estimating dual-energy CT imaging from single-energy CT data with material decomposition convolutional neural network | |
US20080019581A1 (en) | Image Handling and display in X-ray mammography and tomosynthesis | |
US8422764B2 (en) | Method and apparatus of reconstructing 3D image from 2D images | |
CN101232845A (zh) | 磁共振成像装置 | |
CN109389655A (zh) | 时变数据的重建 | |
CN103153192A (zh) | X射线ct装置以及图像再构成方法 | |
CN104899827A (zh) | 基于固定分辨率条件下的离散Radon投影和Mojette投影转换方法 | |
CN106663319A (zh) | 谱图像数据的可视化 | |
CN113272869A (zh) | 医学成像中从定位片进行三维形状重构 | |
Guo et al. | Image reconstruction model for the exterior problem of computed tomography based on weighted directional total variation | |
CN103854297A (zh) | 利用紧框架学习的动态图像重建 | |
CA3057216A1 (en) | Method, device and system for simulating shadow images | |
Sidky et al. | First‐order convex feasibility algorithms for x‐ray CT | |
CN103455703A (zh) | 计算机系统、医用图像诊断装置以及图像显示方法 | |
CN103810712A (zh) | 一种能谱ct图像质量的评价方法 | |
Tolpadi et al. | K2S challenge: From undersampled K-space to automatic segmentation | |
Fessler | Image reconstruction: Algorithms and analysis | |
Li et al. | 3D coronary artery reconstruction by 2D motion compensation based on mutual information | |
Cooper et al. | Computer visualisation of the moving human lumbar spine | |
Anderst | Automated measurement of neural foramen cross-sectional area during in vivo functional movement |
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: 20160726 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. |