CN101065065A - X射线ct装置和图像处理装置 - Google Patents
X射线ct装置和图像处理装置 Download PDFInfo
- Publication number
- CN101065065A CN101065065A CN 200580040425 CN200580040425A CN101065065A CN 101065065 A CN101065065 A CN 101065065A CN 200580040425 CN200580040425 CN 200580040425 CN 200580040425 A CN200580040425 A CN 200580040425A CN 101065065 A CN101065065 A CN 101065065A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- image
- ray
- unit
- value
- 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
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明的X射线CT装置具备:从图像中抽出X射线衰减系数比较高的高对比度区域的单元(20);根据抽出的高对比度区域的位置和装置固有的点像强度分布函数,产生与高对比度区域有关的模糊图像的单元(21,22);为了产生与X射线衰减系数比较低的低对比度区域有关的低对比度图像,从图像减去模糊图像的单元(23)。
Description
技术领域
本发明涉及根据对被检体照射X射线得到的投影数据,生成被检体内部的图像的X射线CT装置和图像处理装置。
背景技术
X射线CT装置的进步很显著,与希望更高精细(高分辨率)并且在广范围内进行摄影的来自医疗现场的强烈愿望相应地,近年来,开发出了多切片X射线CT装置,它已经相当地普及。该多切片X射线CT装置具备:照射向切片方向(卧台的长度方向)扩展而具有宽度的扇形射束X射线的X射线源;具有在切片方向排列了多列(4列、8列、16列等)的检测元件列的构造的2维检测器,是通过多扫描或螺旋扫描使其动作的扫描器。由此,与单切片X射线CT装置相比,能够高精度并且在短时间中得到被检体的广范围内的体数据(volume data)。
通过这样的X射线CT装置生成的CT图像中的血管在近旁有高对比度物质(具有高X射线衰减系数的物质)的情况下,受到其模糊的影响,血管内腔的样子不鲜明。现在,作为不受高对比度区域的影响而改善关注区域的画质的方法,已知以下的方法:在将高对比度区域的CT值设置为与周围的CT值一致的适当的值后进行图像处理,然后,恢复高对比度区域的值。但是,在该方法中,无法良好地除去高对比度区域周围的模糊(参考特开平10-40372号公报)。
发明内容
本发明的目的在于:在显示通过X射线CT装置得到的被检体内部的X射线吸收分布图像的情况下,降低在高X射线衰减系数的物体周围产生的模糊。
本发明的第一方面提供一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,具备:构成为从上述图像中抽出X射线衰减系数比较高的高对比度区域的单元;构成为根据上述抽出的高对比度区域的位置和装置固有的点像强度分布函数,产生与上述高对比度区域有关的模糊图像的单元;构成为为了产生与X射线衰减系数比较低的低对比度区域有关的低对比度图像,从上述图像减去上述模糊图像的单元。
本发明的第二方面提供一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,具备:构成为为了产生与X射线衰减系数比较低的区域有关的低对比度图像,从上述图像减去与包含在上述图像中的X射线衰减系数比较高的区域有关的模糊图像的单元;与CT值对应地将上述低对比度图像分类为多个区域的单元;将上述低对比度图像的像素值置换为上述分类了的多个区域各自固有的值的单元。
本发明的第三方面提供一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,具备:从上述图像中,根据多个阈值,顺序地抽出所包含的大小不同的多个区域候选的单元;计算上述抽出的多个区域候选各自的重心的单元;根据上述重心之间的距离,从上述抽出的多个区域候选中选择特定的区域的单元。
附图说明
图1是本发明的X射线CT装置的结构图。
图2是本发明的实施例1的图像处理部件的结构图。
图3是传递函数的说明图。
图4是本发明的实施例2的图像处理部件的结构图。
图5是本发明的实施例2的分类处理部件的结构图。
图6是表示通过分类处理分割的3个区域的图。
图7是表示图5的区域决定部件的区域决定处理步骤的流程图。
图8是图7的S13的补足图。
图9是表示在图7的S14中抽出的区域R1和在S15中计算的重心位置B1的图。
图10是表示在图7的S14中抽出的区域R2和在S15中计算的重心位置B2的图。
图11是表示在图7的S14中抽出的区域R3和在S15中计算的重心位置B3的图。
具体实施方式
(实施例1)
以下,说明本发明的实施例1。
图1是实施例1的X射线CT装置的结构图。X射线CT装置1具备:构成为收集与被检体有关的投影数据的架台2;用于载置被检体P并使其移动的卧台3;进行用于X射线CT装置的操作的输入和图像显示的操作控制台。
架台1具有X射线管5和X射线检测器6。X射线管5和X射线检测器6被安装在被架台驱动装置7旋转驱动的环状的旋转框架8上。卧台3具备载置被检体的顶板8、用于使顶板8移动的顶板驱动装置9。旋转框架8在其中央部分具有开口。将载置在顶板8上的被检体P插入到该开口中。另外,用Z轴(切片方向轴)定义旋转框架8的旋转中心轴,用XY的垂直2轴定义与Z轴垂直的平面。
在X射线管5的阴极阳极之间,从高电压产生器10施加管电压。从高电压产生器10向X射线管5的灯丝供给灯丝电流。通过管电压的施加和灯丝电流的供给,产生X射线。对于X射线检测器6,采用1维阵列型检测器或2维阵列型检测器(也称为多切片型检测器)。X射线检测元件例如具有0.5mm×0.5mm的正方形的受光面。例如在信道方向上排列916个X射线检测元件。在切片方向上例如并列设置40个该列,则是2维阵列型检测器。由单一列构成的是1维阵列型检测器。
数据收集装置11一般称为DAS(data acquisition system)。数据收集装置11将从检测器6对每个信道输出的信号变换为电压信号并放大,进而变换为数字信号。将该数据(原始数据)供给架台外部的操作控制台4。
操作控制台4的前处理部件12对从数据收集装置11输出的数据(原始数据)实施灵敏度修正等修正处理,输出投影数据。将该投影数据发送到重构处理部件13。重构处理部件13例如根据通过螺旋扫描、使用了锥形射束X射线的体扫描、或同时使用它们而收集到的投影数据,重构图像数据,并存储到图像存储部件14中。图像处理部件15根据存储在图像存储部件14中的图像数据,生成显示图像。根据操作者向输入装置16的输入,进行用于显示图像的条件设置、关注区域的设置等。将在后面详细说明图像处理部件。显示装置17显示在图像处理部件15中生成的图像。另外,操作控制台4的扫描控制部件18根据操作者的输入,控制高压产生装置10、架台驱动装置7、数据收集装置11、顶板驱动装置9,使得进行螺旋扫描等扫描。另外,操作控制台4可以用专用硬件构成,也可以使用计算机用软件实现同样的功能。
图2是说明图1的图像处理部件15的结构的图。另外,图2的(A)是作为图像处理对象的一个例子的血管断面构造的模式图。包括血管壁P1、支架P2、造影剂(血流部分)P3、脂肪P4。以下,以对这样的实体的图像进行了处理的情况为例子进行说明。
图像处理部件15具备关注区域设置图像抽出部件19、高对比度抽出部件20、卷积部件21、CT值比例缩放(scaling)部件22、减法器23、加法器24。
关注区域设置图像抽出部件19根据向输入装置16的输入,设置关注区域,从图像存储部件14抽出该关注区域的图像数据并输出。图2的(B)表示由关注区域设置图像抽出部件19抽出的关注区域的图像。
高对比度抽出部件20从由抽出部件19抽出的关注区域的图像中,抽出X射线衰减系数非常高的高对比度区域,即具有非常高的CT值的像素群。作为高对比度物质,有为了治疗而埋入被检体内的支架等金属性器具、石灰化了的钙等。另外,作为后述的除去了支架等的低对比度图像中的高对比度物质,不是支架等,而是主要将造影剂为对象。选择性地使用具有支架和血管壁之间的CT值的造影剂。
作为具体的抽出处理,通过求出关注像素的值与该关注像素周围的多个像素的值的平均值之间的差分,来求出输入的图像的微分图像,通过使用与支架的边缘对应的阈值对该微分图像实施阈值处理,得到二值化图像。该二值化图像为在可能存在高对比度物质的部分和除此以外的部分中像素值不同的图像。即,该二值化图像是包含可能存在高对比度物质的位置的信息的图像(以下称为高对比度位置图像)。图2的(C)表示该高对比度位置图像。在其微分、二值化处理中,能够以比较少的计算处理容易地得到高对比度位置图像。
卷积部件21针对高对比度位置图像,卷积计算该X射线计算机断层摄影装置固有的点像强度分布函数(PSF:point spread function)。PSF作为定义装置固有的模糊特性的函数,也被称为传递函数或模糊函数。PSF对比检测器6的检测器间距(分辨率限界)的例如0.5mm更微小的例如具有其1/10的直径(0.05mm)的线模拟像进行扫描,得到根据该投影数据通过重构得到的与线模拟像有关的2维图像(模糊图像)的数据。将该PSF与高对比度位置图像进行卷积,由此得到根据单独地对高对比度对象进行扫描取得的数据而重构出的图像,即包含高对比度对象和在其周围出现的模糊的图像。
CT值比例缩放部件22依照作为CT值并不特别的值,即在此为作为造影剂的标准CT值,对从卷积部件21输出的图像(模糊图像)进行标准化(比例缩放)。该标准化处理是用于后述的减法处理的像素值水平的处理。图2的(D)表示从CT值比例缩放部件22输出的图像(以下称为高对比度图像)。
减法部件23进行减法处理,求出关注区域的图像与高对比度图像的差分图像。由此,从关注区域的图像得到降低了高对比度物质及其周围的模糊成分的低对比度图像。该低对比度图像是表示低衰减系数的物质的图像。图2的(E)表示该低对比度图像。
加法部件24通过将高对比度位置图像加上该低对比度图像,来加上高对比度物质的位置信息。图2的(F)表示该输出图像。该输出图像是将高对比度物质的位置信息加到从原来的关注区域的图像减少了高对比度物质造成的模糊(伪像)的图像上的图像,经由后级的处理而显示在显示装置17上。在该输出图像中,能够良好地观察在原来的图像中由于模糊的影响而难以识别的高对比度物质周围的脂肪。
根据这样的实施例1,能够降低在高对比度物质的周围产生的模糊成分,因此能够良好地观察高对比度物质周围的物质。特别在心脏的冠状动脉那样的具有3mm到5mm的直径的毛细血管中有支架或石灰化的情况下,能够良好地观察其周围的血管壁、脂肪的积蓄状态、造影剂的状态等。
也可以对该实施例1进行各种变形而实施。例如,在上述实施例中,加上了高对比度位置图像,但也可以不相加,而只进行低对比度图像的显示。另外,还可以重叠表示高对比度物质的位置的颜色不同的标志。
进而,在上述实施例中,说明了X射线CT装置,但也可以在根据从X射线CT装置输出的投影数据、CT图像等进行图像的显示的医用图像处理装置中实施。进而,在上述实施例中,说明了2维图像的处理,但也可以通过对多个2维图像进行上述的处理,来进行3维图像的处理。
进而,在上述实施例中,说明了图像的像素值是CT值的情况,但如果是表示X射线衰减系数,则也可以是CT值以外的值。另外,在上述实施例中,通过微分、二值化处理求出了高对比度位置图像,但也可以使用以下的方法:对图像逆卷积表示高对比度物质造成的模糊的传递函数,求出高对比度位置图像的方法;预先存储支架的位置、形状的信息,根据该信息求出位置信息的方法等。
(实施例2)
以下,说明本发明的实施例2。另外,省略对与实施例1相同的部分的说明。
实施例2通过分类(clustering)处理,自动地进行分类使得容易识别图像中的类似物质。已知在理想的条件下,血小板或脂肪表示-100到50的CT值,血管壁表示50到129的CT值,包含在血流中的造影剂表示130-350的CT值。但是,在实际的检查中,受到被检体的体格、射束硬化、重构函数、物质的大小、关注区域以外的物质的状态、造影剂的浓度等的影响,CT有变动。因此,在进行区域抽出(区域分割)的方法中,无法对应这样的CT值的变动,无法良好地进行分类。在本实施例中,从低对比度图像中,高精度地抽出主要与血小板对应的区域A、主要与血管壁对应的区域B、主要与造影剂对应的区域C。
图4是实施例2的X射线CT装置1的结构图。在图4中,在减法部件23与加法部件24之间具有分类处理部件25。该分类处理部件25根据CT值求出规定的多个区域,并进行CT值的变换处理,使得用同一亮度或颜色显示各区域。另外,在本实施例中,说明了将CT值分为3个区域的情况,但该区域的个数也可以是其他数目。另外,也可以与操作者的输入对应地变更区域数目。
图5是分类处理部件25的结构图。分类处理部件25具备直方图生成部件26、区域决定部件27、CT值变换部件28。直方图生成部件26求出所输入的低对比度图像的直方图。直方图表示每个CT值的像素出现的频度分布。在此,CT值分布在-100到537之间的范围内。
区域决定部件27如图6所示那样,从低对比度图像中例如抽出3个区域A、B、C。区域A主要对应于血小板。区域B主要对应于血管壁。区域C主要对应于造影剂。
图7表示由区域决定部件27从低对比度图像中抽出3个区域A、B、C的处理的步骤。首先,从低对比度图像中抽出血管壁区域B(S11)。从由操作者指定的血管的大致中心放射状地设置多个跟踪线。沿着各跟踪线跟踪CT值。将CT值的变动比较大的位置,即表示超过规定值的微分值的2个位置确定为血管壁的内点和外点。接着,从低对比度图像中抽出血小板区域A(S12)。作为血小板区域A,抽出位于血管壁的内侧的具有比血管壁低的CT值的像素群。
接着,抽出造影剂的区域C(S13~S19)。首先,如图8所示例的那样,通过k-means法,将直方图上的CT值所分布的范围(在图6中,为-100到537的范围)分割为N个段SEG1~SEGN(S13)。设CT值最高的段为SEG1。将N设置为超过最终分类的区域数n(在此n=3)的整数。理想的是将N设置为最终分类的区域数n的2倍或3倍。m是处理变量。将m初始化为1。
首先,如图9所示那样,从低对比度图像中抽出与段SEG1对应的区域候选R1(S14)。将段SEG1的最低值作为阈值,从低对比度图像中抽出区域R1。计算抽出了的区域候选R1的重心位置B1(S15)。在S17中,将变量m加1(S17)。另外,在S18中,作为用于防止扩散的停止条件,判断“CTmin seg m<CTmax/2”CTmin seg m是段m的最小值,CTmax是低对比度图像的最大CT值。即,在段m的最小值小于低对比度图像的最大CT值的1/2时,停止区域放大处理。在S18中的No时,执行S19。在S18中的Yes时,返回到S14。
如图10所示那样,从低对比度图像中抽出与段SEG1和段SEG2对应的区域候选R2(S14)。实际上,将段SEG2的最低值作为阈值,从低对比度图像中抽出区域候选R2。计算抽出了的区域候选R2的重心位置B2(S15)。计算该区域候选R2的重心位置B2与一个之前的区域候选R1的重心位置B1之间的距离D2,并与规定的阈值Th进行比较(S16)。
直到重心位置的变位量Dm超过阈值Th为止,循环进行S14~S18的处理。即,随着变量m的递增,区域候选R一点一点地扩大。随着区域候选R的扩大,重心B变位。在抽出对象相同时,重心B的变位量比较少。在进行扩大使得抽出对象包含其他物质时,重心B的变位量比较大。即,在抽出对象是造影剂时,即使扩大区域候选R,重心位置也大致不变化。但是,在向抽出对象中加入了血管壁等时,区域候选R飞跃地变大,重心位置也有很大变化。搜索重心位置有很大变化的段。
在m=3时,如图11所示,在区域候选R3的重心位置B3与一个之前的区域候选R2的重心位置B2之间的距离D3超过了阈值Th时(S16),将一个之前的区域候选R2确定为造影剂的区域C(S19)。即,将重心的变位量收敛为小于阈值Th的最大的区域候选R2选择为特定的区域C。
这样,通过着眼于区域扩大造成的重心位置的变位,能够几乎不受到CT值变动因素的影响,高精度地抽出造影剂的区域C。另外,对于血管壁的抽出也一样,通过着眼于CT值变化,能够几乎不受到CT值变动因素的影响,高精度地抽出血管壁的区域B。进而,对于血小板的抽出也一样,通过设置血管壁内侧这样的条件,能够几乎不受到CT值变动因素的影响,高精度地抽出血小板的区域A。
另外,区域决定部件27也可以使用k-mean法,将该直方图分类为3个区域。通过k-mean法,求出用于求出各区域的阈值Th,使得各区域内的CT值的分散值相等。通过决定各区域使得各区域内的CT值的分散值相等,从而容易地使类似性高的物质包含在同一区域内。
图6表示通过本实施例进行分区的例子,将-100到537的CT值分类为3个区域。具有区域A的范围的CT值的像素主要与脂肪的部分对应。具有区域B的范围的CT值的像素主要与血管壁的部分对应。具有区域C的范围的CT值的像素主要与血小板的部分对应。
CT值变换部件28将由区域决定部件27抽出或分类的3个区域分别包含的像素的CT值置换为预先设置的每个区域的固有值。由此,在显示图像上用同一亮度或颜色显示同一区域部分,以能够掌握各区域的形式进行显示。
根据这样的本实施例,由于与CT值对应地显示分区了的图像,能够提供可以良好地掌握被检体内部的组织、脂肪、造影剂等不同物质的图像。另外,根据直方图自动地求出用于决定区域的阈值,因此即使在CT值有变动的情况下,也能够提供良好地分类了的图像。
另外,在该实施例中,在降低了在高对比度物质周围产生的模糊成分(伪像成分)后,将低对比度图像作为对象,进行分类处理,因此能够良好地观察高对比度物质周围的物质。特别在心脏的冠状动脉那样的具有3mm到5mm的直径的毛细血管中有支架或石灰化的情况下,能够良好地观察其周围的血管壁、脂肪的积蓄状态、造影剂的状态等。
另外,由于针对不包含高对比度物质的图像进行分类处理,然后加入高对比度物质的信息,所以能够良好地对具有低CT值的部分进行分区。
另外,本发明并不只限于上述实施例,在实施阶段,在不脱离其宗旨的范围内,能够对构成要素进行变形而具体化。另外,通过适当地组合上述实施例所揭示的多个构成要素,能够形成各种发明。例如,可以从实施例所示的全部构成要素中删除若干个构成要素。进而,也可以适当地组合不同的实施例中的构成要素。
例如,在上述实施例中,使用k-mean法求出阈值,但也可以使用群平均法、Ward法等其他分类方法求出阈值。另外,在上述实施例中,针对除去了高对比度物质及其周围的模糊的图像进行分类处理,但也可以不是这样的处理,而针对关注区域的原始图像,直接进行分类处理,生成显示图像。
进而,在上述实施例中,对高对比度位置图像进行相加,但也可以不进行相加,而只显示低对比度图像。另外,也可以重叠显示表示高对比度物质的位置的颜色不同的标记。
进而,在上述实施例中,说明了X射线CT装置,但也可以在根据从X射线CT装置输出的投影数据CT图像等进行图像显示的医用图像处理装置中实施。进而,在上述实施例中,说明了2维图像的处理,但也可以通过对多个2维图像进行上述处理,来进行3维图像的处理。进而,在上述实施例中,说明了图像的像素值是CT值的情况,但如果是表示X射线衰减系数的值,则也可以是CT值以外的值。
根据本发明,在显示通过X射线CT装置得到的被检体内部的X射线吸收分布图像的情况下,能够降低在高X射线衰减系数的物体周围产生的模糊。
Claims (23)
1.一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,其特征在于包括:
构成为从上述图像中抽出X射线衰减系数比较高的高对比度区域的单元;
构成为根据上述抽出的高对比度区域的位置和装置固有的点像强度分布函数,产生与上述高对比度区域有关的模糊图像的单元;
构成为为了产生与X射线衰减系数比较低的低对比度区域有关的低对比度图像,从上述图像减去上述模糊图像的单元。
2.根据权利要求1所述的X射线CT装置,其特征在于:
上述抽出高对比度区域的单元根据上述图像产生微分图像,通过对上述微分图像进行阈值处理,而抽出上述高对比度区域。
3.根据权利要求1所述的X射线CT装置,其特征在于:
上述产生模糊图像的单元具备:
对上述抽出的高对比度区域的位置图像,卷积上述点像强度分布函数的单元;
依照上述高对比度区域的标准CT值,对上述进行了卷积的高对比度区域的位置图像进行标准化的单元。
4.根据权利要求1所述的X射线CT装置,其特征在于:
上述高对比度区域是为了治疗而埋入到被检体内的金属性的器具或石灰化了的钙。
5.根据权利要求1所述的X射线CT装置,其特征在于还包括:
将上述抽出的高对比度区域的位置图像与上述低对比度图像相加的单元。
6.根据权利要求1所述的X射线CT装置,其特征在于还包括:
与CT值对应地将上述低对比度图像分类为多个区域的单元;
将上述低对比度图像的像素值置换为上述分类了的多个区域各自固有的值的单元。
7.根据权利要求6所述的X射线CT装置,其特征在于还包括:
将上述抽出的高对比度区域的位置图像与上述置换后的低对比度图像相加的单元。
8.根据权利要求6所述的X射线CT装置,其特征在于:
上述分类了的多个区域包括主要是脂肪的区域、主要是血管壁的区域、主要是造影剂的区域。
9.根据权利要求6所述的X射线CT装置,其特征在于:
上述分类单元具备:
根据多个阈值,从上述低对比度图像中抽出大小不同的多个区域的单元;
计算上述抽出的多个区域各自的重心的单元;
根据上述重心之间的距离,从上述抽出的多个区域中选择与抽出对象物对应的特定的区域的单元。
10.根据权利要求9所述的X射线CT装置,其特征在于:
根据与上述图像有关的每个CT值的像素出现的频度分布,决定上述多个阈值。
11.一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,其特征在于包括:
构成为为了产生与X射线衰减系数比较低的区域有关的低对比度图像,从上述图像减去与包含在上述图像中的X射线衰减系数比较高的区域有关的模糊图像的单元;
与CT值对应地将上述低对比度图像分类为多个区域的单元;
将上述低对比度图像的像素值置换为上述分类了的多个区域各自固有的值的单元。
12.根据权利要求11所述的X射线CT装置,其特征在于:
上述分类单元将上述低对比度图像分类为造影剂区域、血管壁区域和血小板区域。
13.根据权利要求12所述的X射线CT装置,其特征在于:
上述分类单元为了确定上述造影剂区域,具备:
根据多个阈值,从上述低对比度图像中顺序抽出内部包含的大小不同的多个区域候选的单元;
计算上述抽出的多个区域候选各自的重心的单元;
根据上述重心之间的距离,从上述抽出的多个区域候选中选择上述造影剂区域的单元。
14.根据权利要求13所述的X射线CT装置,其特征在于:
根据与上述图像有关的每个CT值的像素出现的频度分布,决定上述多个阈值。
15.根据权利要求13所述的X射线CT装置,其特征在于:
上述选择单元选择上述重心的变位量收敛为小于规定量的最大的区域作为特定的区域。
16.根据权利要求13所述的X射线CT装置,其特征在于:
上述分类单元为了确定上述血管壁区域,从上述低对比度图像的血管大致中心放射状地跟踪CT值,将上述CT值的变化超过规定值的2个位置确定为血管的内壁和外壁。
17.根据权利要求16所述的X射线CT装置,其特征在于:
上述分类单元为了确定上述血小板区域,抽出上述血管壁内侧的、具有比上述血管壁低的CT值的像素群。
18.一种X射线CT装置,是在收集与被检体有关的投影数据,根据上述收集到的投影数据而重构上述被检体的内部的图像的X射线CT装置中,其特征在于包括:
从上述图像中,根据多个阈值,顺序地抽出内部包含的大小不同的多个区域候选的单元;
计算上述抽出的多个区域候选各自的重心的单元;
根据上述重心之间的距离,从上述抽出的多个区域候选中选择特定的区域的单元。
19.根据权利要求18所述的X射线CT装置,其特征在于:
根据与上述图像有关的每个CT值的像素出现的频度分布,决定上述多个阈值。
20.根据权利要求18所述的X射线CT装置,其特征在于:
上述选择单元选择上述重心的变位量收敛为小于规定量的最大的区域作为特定的区域。
21.一种图像处理装置,其特征在于包括:
构成为从X射线图像中抽出X射线衰减系数比较高的高对比度区域的单元;
构成为根据上述抽出的高对比度区域的位置和装置固有的点像强度分布函数,产生与上述高对比度区域有关的模糊图像的单元;
构成为为了产生与X射线衰减系数比较低的低对比度区域有关的低对比度图像,从上述图像减去上述模糊图像的单元
22.一种图像处理装置,其特征在于包括:
构成为从X射线图像减去与包含在上述X射线图像中的X射线衰减系数比较高的区域有关的模糊图像,产生与X射线衰减系数比较低的区域有关的低对比度图像的单元;
与像素值对应地将上述低对比度图像分类为多个区域的单元;
将上述低对比度图像的像素值置换为上述分类了的多个区域各自固有的值的单元。
23.一种图像处理装置,其特征在于包括:
根据多个阈值,从X射线图像中顺序抽出内部包含的大小不同的多个区域候选的单元;
计算上述抽出的多个区域候选各自的重心的单元;
根据上述重心之间的距离,从上述抽出的多个区域候选中选择特定的区域的单元。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP342273/2004 | 2004-11-26 | ||
JP342272/2004 | 2004-11-26 | ||
JP2004342272 | 2004-11-26 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101065065A true CN101065065A (zh) | 2007-10-31 |
CN100589760C CN100589760C (zh) | 2010-02-17 |
Family
ID=38965623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200580040425A Expired - Fee Related CN100589760C (zh) | 2004-11-26 | 2005-11-24 | X射线ct装置和图像处理装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100589760C (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102688060A (zh) * | 2011-03-23 | 2012-09-26 | 株式会社东芝 | 图像处理装置、x 线ct 装置以及图像处理方法 |
CN101727661B (zh) * | 2008-10-10 | 2012-10-03 | 株式会社东芝 | 医用图像处理装置 |
CN104517263A (zh) * | 2013-09-30 | 2015-04-15 | Ge医疗系统环球技术有限公司 | 减少计算机断层扫描图像重构中伪像的方法和装置 |
CN104580879A (zh) * | 2013-10-09 | 2015-04-29 | 佳能株式会社 | 图像处理设备、图像拾取设备以及图像处理方法 |
CN104602606A (zh) * | 2012-08-31 | 2015-05-06 | 株式会社岛津制作所 | 放射线断层图像生成装置、放射线断层摄影装置以及放射线断层图像生成方法 |
CN105208937A (zh) * | 2013-05-10 | 2015-12-30 | 株式会社岛津制作所 | 图像处理装置 |
CN107106109A (zh) * | 2014-11-06 | 2017-08-29 | 皇家飞利浦有限公司 | 计算机断层扫描系统 |
CN107392893A (zh) * | 2017-06-30 | 2017-11-24 | 上海联影医疗科技有限公司 | 组织密度分析方法及系统 |
WO2019000418A1 (zh) * | 2017-06-30 | 2019-01-03 | 上海联影医疗科技有限公司 | 组织密度分析方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6266388B1 (en) * | 1999-09-07 | 2001-07-24 | General Electric Company | Methods and apparatus for two-pass cone beam image reconstruction |
CN1261910C (zh) * | 2003-12-08 | 2006-06-28 | 西安理工大学 | 彩色多窗ct图像的自动生成方法 |
-
2005
- 2005-11-24 CN CN200580040425A patent/CN100589760C/zh not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727661B (zh) * | 2008-10-10 | 2012-10-03 | 株式会社东芝 | 医用图像处理装置 |
CN102688060B (zh) * | 2011-03-23 | 2015-04-01 | 株式会社东芝 | 图像处理装置、x 线ct 装置以及图像处理方法 |
CN102688060A (zh) * | 2011-03-23 | 2012-09-26 | 株式会社东芝 | 图像处理装置、x 线ct 装置以及图像处理方法 |
CN104602606A (zh) * | 2012-08-31 | 2015-05-06 | 株式会社岛津制作所 | 放射线断层图像生成装置、放射线断层摄影装置以及放射线断层图像生成方法 |
CN104602606B (zh) * | 2012-08-31 | 2017-06-23 | 株式会社岛津制作所 | 放射线断层图像生成装置以及放射线断层图像生成方法 |
CN105208937B (zh) * | 2013-05-10 | 2018-04-03 | 株式会社岛津制作所 | 图像处理装置 |
CN105208937A (zh) * | 2013-05-10 | 2015-12-30 | 株式会社岛津制作所 | 图像处理装置 |
CN104517263A (zh) * | 2013-09-30 | 2015-04-15 | Ge医疗系统环球技术有限公司 | 减少计算机断层扫描图像重构中伪像的方法和装置 |
CN104517263B (zh) * | 2013-09-30 | 2019-06-14 | Ge医疗系统环球技术有限公司 | 减少计算机断层扫描图像重构中伪像的方法和装置 |
CN104580879A (zh) * | 2013-10-09 | 2015-04-29 | 佳能株式会社 | 图像处理设备、图像拾取设备以及图像处理方法 |
CN104580879B (zh) * | 2013-10-09 | 2018-01-12 | 佳能株式会社 | 图像处理设备、图像拾取设备以及图像处理方法 |
US9747672B2 (en) | 2013-10-09 | 2017-08-29 | Canon Kabushiki Kaisha | Image processing apparatus, image pickup apparatus, image processing method, and non-transitory computer-readable storage medium |
CN107106109A (zh) * | 2014-11-06 | 2017-08-29 | 皇家飞利浦有限公司 | 计算机断层扫描系统 |
CN107106109B (zh) * | 2014-11-06 | 2020-11-06 | 皇家飞利浦有限公司 | 计算机断层扫描系统 |
CN107392893A (zh) * | 2017-06-30 | 2017-11-24 | 上海联影医疗科技有限公司 | 组织密度分析方法及系统 |
WO2019000418A1 (zh) * | 2017-06-30 | 2019-01-03 | 上海联影医疗科技有限公司 | 组织密度分析方法及系统 |
US11664112B2 (en) | 2017-06-30 | 2023-05-30 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for tissue density analysis |
Also Published As
Publication number | Publication date |
---|---|
CN100589760C (zh) | 2010-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101065065A (zh) | X射线ct装置和图像处理装置 | |
JP5159937B2 (ja) | X線ct装置及び画像処理装置 | |
US7835562B2 (en) | Methods and apparatus for noise reduction filtering of images | |
CN100350434C (zh) | 图像数据区域抽取系统以及图像数据区域抽取方法 | |
JP4751495B2 (ja) | 構造画像成分を識別する方法及び装置 | |
US7916912B2 (en) | Efficient border extraction of image feature | |
EP1048008B1 (en) | Methods and apparatus for calcification scoring | |
EP2831627B1 (en) | Projection data binning and image generation in photon counting imaging modality | |
EP2891455B1 (en) | Radiation tomographic image-generating apparatus and radiation tomographic image-generating method | |
US8045781B2 (en) | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus | |
US9659390B2 (en) | Tomosynthesis reconstruction with rib suppression | |
KR101493375B1 (ko) | 화상처리장치, 화상처리방법, 및 컴퓨터 판독가능한 기억매체 | |
WO2005008586A2 (en) | Metal artifact correction in computed tomography | |
CN1188949A (zh) | 通过消除主结构增强数据显象能力 | |
CN1947151A (zh) | 用于在图像中使用发散梯度场响应进行基于滑降的对象分割的系统和方法 | |
JP5804340B2 (ja) | 放射線画像領域抽出装置、放射線画像領域抽出プログラム、放射線撮影装置および放射線画像領域抽出方法 | |
CN1706344A (zh) | X射线ct装置和心肌灌注图像产生系统 | |
JP5329103B2 (ja) | 画像処理装置及びx線ct装置 | |
JP5214110B2 (ja) | X線ct装置 | |
JP2004222776A (ja) | 異常陰影候補検出装置 | |
WO2009004571A1 (en) | Method and apparatus for image reconstruction | |
KR102422613B1 (ko) | 인공지능 기반 실시간 자동 엑스레이 폐결절 판독 방법 및 시스템 | |
CN1305298C (zh) | 图像重构方法和x射线ct装置 | |
Koonsanit et al. | Metal Artifact Removal on Dental CT Scanned Images by Using Multi-Layer Entropic Thresholding and Label Filtering Techniques for 3-D Visualization of CT Images |
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: 20160706 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. |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100217 Termination date: 20181124 |
|
CF01 | Termination of patent right due to non-payment of annual fee |