专利名称：边缘评价方法、边缘检测方法、图像修正方法及图像处理系统的制作方法
技术领域：
0001
本发明涉及边缘评价方法、边缘检测方法、图像修正方法及图像处理 系统，特别是涉及评价构成图像数据的像素为边缘部的可能性的技术。
背景技术：
0002
以往，对从x射线拍摄装置或非破坏检査装置等得到的图像数据进行
边缘检测。该处理是在为了上下亮度差的假象回避等之类的图像质量提 高、或自动检测组织/病变部等中多采用的重要步骤。
0003
最基本的方法是取与相邻像素的像素值之差，且若该差大则将该像素
作为边缘部/组织边界部。另外，也存在对图像数据施加高通滤波(HPF:
High Pass Filter)来强调变动，且若该值大则将该像素作为边缘部/组织边 界部的方法(例如，参照专利文献l)。
专利文献l:日本特开平7-134418号公报0004
可是，上述的方法存在在像素值由噪声等影响而不稳定的情况下边缘 检测的精度降低的问题。对此，例如存在导入"轮廓线是封闭曲线"的假 设来精度维持/提高的方法(Snake法等)。但是，却或者存在需要手动操 作的预处理或者存在在从假设偏离的情况中精度降低的问题
发明内容
0005
本发明鉴于上述问题进行的，其目的在于提供一种能够高精度评价像 素为边缘部的可能性的边缘评价方法、边缘检测方法、图像修正方法及图像处理系统。
0006
本发明为了达成上述目的，选择以下的构成。
艮P，本发明的边缘评价方法其特征在于，在对构成图像数据的像素为
边缘部的可能性进行评价的边缘评价方法中，该边缘评价方法具有对各 像素计算相对于相邻像素的像素值的变动量的变动量计算过程、对任意一 个像素基于由该一个像素及其周边的像素构成的周边像素组的各变动量 来计算对所述一个像素为边缘部的可能性进行数值化的边缘程度的边缘 程度计算过程。
0007
根据本发明的边缘评价方法，由于边缘程度计算过程基于一个像素及 其周边的像素的多个变动量来计算一个像素的边缘程度，因此能够高精度 评价一个像素为边缘部的可能性。0008
另外，本发明的图像处理系统其特征在于，在对构成图像数据的像素
为边缘部的可能性进行评价的图像处理系统中，该图像处理系统具有对 各像素计算相对于相邻像素的像素值的变动量的变动量计算单元和对任 意一个像素基于由该一个像素及其周边的像素构成的周边像素组的各变 动量来计算对所述一个像素为边缘部的可能性进行数值化的边缘程度的 边缘程度计算单元。
0009
根据本发明的图像处理系统，由于边缘程度计算单元基于一个像素及 其周边的像素的多个变动量计算一个像素的边缘程度，因此能够高精度评 价一个像素为边缘部的可能性。0010
在这些发明的边缘评价方法及图像处理系统中，在将所述周边像素组 的各变动量置换为绝对值后进行合计的值作为变动幅度总和，且将直接合 计所述周边像素组的各变动量后的值的绝对值作为变动量增减幅度时，所 述边缘程度优选变动量增减幅度相对于所述变动幅度总和的比例或与该
比例相应的值。边缘程度表示像素的周边像素组的变动量与增加/减少的任一越相同则一个像素为边缘部的可能性越变高，且周边像素组的变动量 的增加/减少越偏差则为边缘部的可能性越变低的倾向。由于边缘程度具有 这种倾向，因此通过噪声影响等能够防止边缘程度的精度降低。
0011
另外，在这些发明的边缘评价方法及图像处理系统中，与一个像素对 应的周边像素组优选以一个像素为大致中心的区域内的像素集合。能够高 精度地求出像素的边缘程度。
0012
在本发明的边缘评价方法及图像处理系统中，图像数据优选二维图像 数据或三维图像数据。即使是二维图像数据或三维图像数据也能适用。0013
另外，在本发明的边缘评价方法中，所述相邻像素作为与各像素相邻 的方向不同的多个像素，所述变动量计算过程对各像素按每个所述相邻像 素相邻方向来分别计算所述变动量，所述边缘程度计算过程优选对一个像 素按每个与各方向对应的变动量来分别计算所述边缘程度。像素值变动的 斜率即使是任意方向，也能适当地计算边缘程度。0014
且有，在本发明的边缘评价方法中，所述边缘评价计算过程优选还对 各方向的所述边缘程度进行合成而取得合成边缘程度。即使像素值变动的 斜率是任意方向，也能适当计算合成边缘程度。0015
在本发明的图像处理系统中，所述相邻像素作为相对于各像素而相邻 的方向不同的多个像素，所述变动量计算单元对各像素按每个所述相邻像 素相邻的方向分别计算所述变动量，所述边缘程度计算单元优选对一个像 素按每个与各方向对应的变动量来分别计算所述边缘程度。像素值变动的 斜率即使是任意方向，也能适当计算合成边缘程度。0016
在本发明的图像处理系统中，所述边缘程度计算过程优选还对各方向 的所述边缘程度进行合成而取得合成边缘程度。像素值变动的斜率即使是 任意方向，也能适当计算合成边缘程度。0017
在本发明的图像处理系统中，所述系统优选还具备边缘部检测单元， 该边缘部检测单元根据将所述边缘程度与阈值比较后的结果来判断所述 一个像素为边缘部。能够高精度地检测为边缘部的像素。0018
在本发明的图像处理系统中，所述系统优选具备图像修正单元，该图 像修正单元根据所述边缘程度对所述一个像素的像素值进行修正。根据边 缘程度能够适当修正图像数据。0019
在本发明的图像处理系统中，所述边缘程度计算单元优选对构成所述 图像数据的各像素计算所述边缘程度。对图像数据的全部像素能够计算边 缘程度。
0020
在本发明的图像处理系统中，所述系统还具备指定作为所述边缘程度 的计算对象的像素的输入单元，所述边缘程度计算单元优选由所述输入单 元指定的像素来计算所述边缘程度。能够适当指定计算图像数据的边缘程 度的像素。0021
在本发明的图像处理系统中，所述系统优选还具有输出所述边缘程度 的输出单元。能够将计算的边缘程度适当地输出到外部。0022
在本发明的图像处理系统中，优选构成为可从外部输入图像数据。能 够对各种图像数据进行计算边缘程度的处理。0023
另外，本发明的边缘检测方法其特征在于，在检测构成图像数据的像
素为边缘部的边缘检测方法中，该边缘检测方法具有对各像素计算相对 于相邻像素的像素值的变动量的变动量计算过程、对任意一个像素基于由 该一个像素及其周边的像素构成的周边像素组的各变动量来计算对所述 一个像素为边缘部的可能性进行数值化的边缘程度的边缘程度计算过程 和根据将所述边缘程度与阈值比较后的结果来判断所述一个像素为边缘部的边缘检测过程。
0024
根据本发明的边缘检测方法，由于边缘程度计算过程基于一个像素及 其周边的像素的多个变动量来计算一个像素的边缘程度，因此能够高精度 地评价一个像素为边缘部的可能性。由此，能够高精度地检测为边缘部的 像素。
0025
另外，本发明的图像修正方法其特征在于，在对构成图像数据的像素
的像素值进行修正的图像修正方法中，该图像修正方法具有对各像素计 算相对于相邻像素的像素值的变动量的变动量计算过程、对任意一个像素 基于由该一个像素及其周边的像素构成的周边像素组的各变动量来计算 对所述一个像素为边缘部的可能性进行数值化的边缘程度的边缘程度计 算过程和根据所述边缘程度对所述一个像素的像素值进行修正的图像修 正过程。
0026
根据本发明的图像修正方法，由于边缘程度计算过程基于一个像素及 其周边的像素的多个变动量来计算一个像素的边缘程度，因此能够高精度 地评价一个像素为边缘部的可能性。由此，根据边缘程度能够适当生成修 正的图像数据。 (发明效果)0027
根据本发明的边缘评价方法、边缘检测方法、图像修正方法及图像处 理系统，由于边缘程度计算过程基于一个像素及其周边的像素的多个变动 量来计算一个像素的边缘程度，因此能够高精度地评价一个像素为边缘部 的可能性。
0028
图1是表示实施例1相关的X射线拍摄装置的大致构成的框图。 图2是FPD主要部分的垂直剖视图。图3是FPD平面图。
图4是表示X射线拍摄装置的动作步骤的流程图。 图5是表示输入图像数据的一个例子的示意图。 图6是周边像素组的示意图。
图7是表示实施例2相关的X射线拍摄装置的大致构成的框图。 图8是表示X射线拍摄装置的动作步骤的流程图。 图9是表示实施例3相关的X射线拍摄装置的大致构成的框图。 图IO是表示X射线拍摄装置的动作步骤的流程图。 图ll是周边像素组的示意图。0029
图中
13—图像处理部，
15 —输入部，
17—输出部，
21 —变动量计算部，
23 —边缘程度计算部，
25 —边缘检测部，
26—图像修正部，
M—被摄体。
具体实施方式
实施例0030
以下，参照附图对本发明的实施例l进行说明。
图1是表示实施例1相关的x射线拍摄装置的大致构成的框图。
0031
桌面1是由X射线透过材料等构成，载置了作为检査对象的被摄体M。 X射线管3和平板型X射线检测器(以下，单单作为"FPD"进行适当记 载)5夹持桌面1而对置配置。X射线管3将X射线照射到被摄体M上。 FPD5检测透过了被摄体M后的X射线并输出与X射线强度相应的检测信号。X射线管3和FPD5与摄影控制部7连接，在该摄影控制部7上连 接了省略图示的操作部。操作部接受来自摄影者的指示(例如，有关X射 线管3的X射线照射强度或照射时间或FPD5中的放大度等指示)而向摄 影控制部7输出。摄影控制部7基于该操作者的指示从X射线管3照射X 射线并由FPD5中检测X射线。0032
收集部11收集FPD5输出的检测信号，并将与收集到的检测信号相应 的输入图像数据输出到图像处理部13。图像处理部13具备变动量计算部 21、边缘程度计算部23及边缘检测部25，并对输入图像数据进行图像处 理。
0033
变动量计算部21对构成输入图像数据的各像素进行求出相对于相邻 像素的像素值的变动量的运算处理。变动量有增加和减少的任意方向，对 应其取正或负的值。0034
边缘程度计算部23基于一个像素及周边的像素(以下，将这些像素 统称为周边像素组)的各变动量而进行求出对该像素为边缘部的可能性进 行评价(数值化)了的边缘程度的运算处理。0035
更具体地说，边缘程度计算部23在将周边像素组的各变动量置换为 绝对值后将合计的值作为"变动幅度总和"且将直接合计了周边像素组的 各变动量的值的绝对值作为"变动量增减幅度"时，将"变动量增减幅度" 相对于"变动幅度总和"的比例或将与其相应的值作为该像素的边缘程 度。
0036
如上述，因为边缘程度是像素为边缘部的可能性的指标，所以边缘程 度优选利用从0到1的小数或百分比(％)来表示。且有，所谓"像素为 边缘部"是指与被摄体M整体或部分的轮廓/边界对应的像素。0037
边缘检测部25根据边缘程度与阈值比较后的结果来判断像素为边缘部。
0038
该图像处理部13是利用读取规定程序来执行的中央运算处理装置 (CPU)或存储各种信息的RAM(Random-Access Memory)或固定磁盘等 存储介质等来实现的。0039
另外，在图像处理部13连接了输入部15和输出部17。输入部15由 键盘或鼠标等构成，输出部17由监视器等构成。输入部15接受来自操作 者的指示并输出到图像处理部13。作为指示，例如存在阈值设定相关的指 示、边缘程度的计算或指定成为边缘检测对象的像素的指示、各运算处理 结果的输出相关的指示等。输出部17根据操作者指示适当输出图像处理 部13的运算处理结果(变动量、边缘程度、边缘部的像素等)。且有， 图像处理部13、输入部15及输出部17相当于本发明的图像处理系统。0040
参照图2、图3对FPD5的构成进行说明。图2是FPD主要部分的垂 直剖视图，图3是FPD的平面图。0041
FPD5具有施加电极3K X射线感应型的半导体膜33、载波收集电 极35和有源矩阵基板37，将他们从X射线入射侧起依次进行层叠。另外， FPD5还具备栅极驱动器41和放大器43等。0042
载波收集电极35俯视二维矩阵状地分离形成。在有源矩阵基板37上 分离形成了按每个载波收集电极35蓄积电荷信息的电容器Ca和作为取出 电荷信息的开关元件的薄膜晶体管(Thin Film Transistors) Tr。载波收集 电极35和电容器Ca与薄膜晶体管Tr的源极连接。这些一组的载波收集 电极35、电容器Ca和薄膜晶体管Tr构成一个检测元件d。检测元件d在 FPD5的检测面上矩阵状地排列了多个(例如，在纵30cmX横30cm程度 的宽度上为4096行X4096列)。各检测元件d由于与上述的像素存在对 应关系，因此输入到图像处理部13中的输入图像数据分辨率高。例如， 一个像素的大小为150pm见方。0043
在有源矩阵基板37上按每行检测元件d铺设形成了栅极总线45，按 每列检测元件d铺设形成了数据总线47。各栅极总线45分别将各行的薄 膜晶体管Tr的栅极汇总而与栅极驱动器41电连接。另外，各数据总线47 将各列的薄膜晶体管Tr的漏极与放大器43电连接。在放大器43的输出 侧设有省略图示的A/D变换器等。0044
并且，若以向施加电极31施加偏置电压的状态向FPD5射入X射线， 则在半导体膜33产生电荷，该电荷经由各载波收集电极35而被蓄积在电 容器Ca中。栅极总线45发送来自栅极驱动器41的扫描信号并给予薄膜 晶体管Tr的栅极。由此，经由转移到导通状态的薄膜晶体管Tr而将被蓄 积在电容器Ca中的电荷信息读取到数据总线47。通过各数据总线47而 读取的电荷信息分别利用放大器43进行放大。其后，利用A/D变换器进 行数字化并作为检测信号输出。
0045
接着，对实施例1相关的X射线拍摄装置的动作，以图像处理部13 中的处理为中心进行说明。图4是表示X射线拍摄装置的动作步骤的流程 图。另外，以下，对将全部像素作为对象而进行边缘检测的情况进行说明。0046
<步骤81>进行X射线摄影
在摄影控制部7的控制下，X射线管3向被摄体M照射X射线。另 外，FPD5检测透过了被摄体M的X射线，并将得到的检测信号输出到收 集部11。收集部11将与收集的检测信号相应的输入图像数据输出到图像 处理部13中。0047
<步骤82>计算变动量(变动量计算过程)
变动量计算部21对各像素进行求出变动量的运算处理。参照图5进 行具体说明。图5是表示输入图像数据的一个例子的示意图。如图示，将 输入图像数据作为在正交的纵横2个方向上矩阵状地排列像素的二维图像 数据。为了说明的方便，将横方向作为X轴，将纵方向作为Y轴，并通过横方向的行号码x和纵方向的列号码y来识别各像素。另外，将各像素 的像素值表示为A (x,y)。且有，在图5中，只用"A"表示，省略了矩 阵号码(x,y)。另外，在以下的说明中，也适当地省略行列号码(x，y) 而进行记载。0048
变动量计算部21对各像素计算相对于横方向的相邻像素的像素值A 的变动量SX。例如，将变动量SX(x,y)作为右边的相邻像素的像素值A( x + 1， y)与左边的相邻像素的像素值A — l， y)之差除以2后的 值，对各像素计算变动量SX。此时的变动量SX与像素值A的关系式表 示在公式1中。
0049
公式1
<formula>formula see original document page 15</formula>0050
计算出的各像素的变动量SX被输出到边缘程度计算部23中。0051
〈步骤S3〉计算边缘程度(边缘程度计算过程)
边缘程度计算部23进行基于与各像素相应的周边像素组的变动量SX 来分别求出各像素的边缘程度的运算处理。作为周边像素组，例示了将该 像素作为中心的9行9列、共81个像素。图6是表示该情况下的周边像 素组的示意图。如图示，各像素与图5同样沿着X轴及Y轴矩阵状地排 列，为了方便将各像素的变动量作为"SX"。并且，像素(x，y)的周边 像素组g(x,y)是图示的81个像素的集合。0052
边缘程度计算部23将边缘程度RX作为"变动量增减幅度"相对于 上述的"变动幅度总和"的比例，并对各像素计算边缘程度RX。此时的 边缘程度RX和变动量SX的关系式表示在公式2中。0053
公式2<formula>formula see original document page 16</formula>
0054
由此计算出的边缘程度RX取O以上1以下的值，越接近于1像素为 边缘部的概率越高。各像素的变动量RX被输出到边缘检测部25中。0055
<步骤84>检测边缘(边缘检测过程)
边缘检测部25将各像素的边缘程度RX与阈值进行比较。作为阈值 例示了0.5等。并且，仅在边缘程度RX为阈值以上的情况下，判断为像 素为边缘部。0056
并且，根据操作者的指示，图像处理部13适当地将变动量计算部21、 边缘程度计算部23及边缘检测部25的各运算处理的结果(例如，变动量 SX的分布、边缘程度RX的分布或判断为边缘部的像素的分布等)输出 到输出部17中。且有，边缘检测部25的处理结果构成对输入图像数据进 行边缘强调的2值化图像数据。0057
由此，根据实施例1相关的X射线拍摄装置具备的图像处理系统，因 为边缘程度计算部23基于周边像素组g的多个变动量SX算出边缘程度 RX，所以能够精度高地评价像素为边缘部。
0058
具体地说，边缘程度RX具有如下倾向周边像素组g的变动量SX 的方向与增加或减少越一致则越接近1，且周边像素组g的变动量SX的 方向越有偏差则越接近0。其中，变动量SX越小则通过噪声影响而变动 量SX的方向越容易反转。由此，成为如下的倾向周边像素组g的各变 动量SX越小则边缘程度RX越容易降低。换言之，求取边缘程度RX的 像素也只是周边像素组g的一个像素，虽说该像素接受了噪声，但边缘程 度RX值不会受到大影响。0059
另一方面，在各变动量sx大到超过噪声影响的程度的情况下，不必
考虑由噪声影响而引起的变动量SX的方向反转。此时，边缘程度RX成 为如下倾向周边像素组g的各变动量SX的方向(增加或减少)相同的 程度越高则边缘程度RX越高，虽说各像素接受了噪声影响但边缘程度RX 的值却没有降低。0060
由此，因为基于周边像素组g的各变动量sx而不是基于一个像素的 变动量SX来计算边缘程度RX，所以即使该像素接受了噪声也能够高精 度地评价是否是边缘部。
0061
特别是，根据如FPD5利用高分辨率收集的输入图像数据，与边缘部 相比像素尺寸小且边缘部的幅度涉及多个像素的情况多。由此，相对于边 缘部，像素的大小越小则边缘程度RX的精度越提高。其结果，在其后的 处理中，能够避免假象的发生。另外，能够提高自动诊断的可能性。0062
另外，通过将这种边缘程度RX与阈值进行比较的边缘检测部25，能 够适当地检测为边缘部的像素。
实施例0063
以下，参照附图对本发明的实施例2进行说明。图7是表示实施例2 相关的X射线拍摄装置的大致构成的框图。且有，对与实施例1相同的构 成附加同一符号并省略详细说明。
0064
实施例2中的图像处理部13具备基于边缘程度进行修正输入图像数 据的运算处理的图像修正部26。在本实施例中，图像修正部26生成将各 像素的像素值置换为与该边缘程度相应的灰度值而构成的修正图像。0065
接着，对实施例2相关的X射线拍摄装置的动作以图像处理部13中的处理为中心进行说明。图8是表示X射线拍摄装置的动作步骤的流程图。
且有，对在实施例1中说明的处理进行简略说明。
0066
〈步骤S1〉进行x射线摄影
通过摄影控制部7的控制进行X射线摄影。收集部11将输入图像数 据输出到图像处理部13中。0067
〈步骤S2〉计算变动量(变动量计算过程)
变动量计算部21对各像素计算变动量SX并且进行求出相对于纵方向 的相邻像素的像素值A的变动量SY的运算处理。以下，适当地区分横方 向的变动量SX、纵方向的变动量SY并进行称呼。0068
作为纵方向的变动量SY(x，y)，例示了将与下侧相邻的相邻像素的像 素值A (x，y+l)和与上侧相邻的相邻像素的像素值A (x,y-l)之差用2 除后的值。此时的变动量SY与像素A的关系式在公式3中示出。
0069
公式3<formula>formula see original document page 18</formula>
对各像素，按每个相邻像素相邻的方向，将计算出的变动量SX、 SY 被输出到边缘程度计算部23中。0071
<步骤83>计算边缘程度(边缘程度计算过程)
边缘程度计算部23与实施例1同样，基于与各像素对应的周边像素 组g的变动量SX来分别计算各像素的边缘程度RX。还基于与各像素相 应的周边像素组g的纵方向的变动量SY来分别计算边缘程度RY。以下， 将这些适当地区别为横方向的边缘程度RX、纵方向的边缘程度RY并进 行称呼。
0072且有，作为纵方向的边缘程度RY，例示了纵方向的"变动量增减幅 度"相对于纵方向的"变动幅度总和"的比例。此时的边缘程度RY与变 动量SY的关系式在公式4中表示。
0073<formula>formula see original document page 19</formula>0074
边缘程度计算部23还对各像素计算合成了横方向的边缘程度RX和 纵方向的边缘程度RY后的合成边缘程度R。例如，作为合成边缘程度R， 例示了边缘程度RX的平方与边缘程度RY的平方和用2除后的值的平方 根。此时的合成边缘程度R与横方向、纵方向的边缘程度RX、 RY的关 系在公式5中示出。且有，合成边缘程度R取0以上1以下的值。
0075<formula>formula see original document page 19</formula>0076
计算出的各像素的合成边缘程度R被输出到图像修正部26中。0077
<步骤86>生成修正图像数据(图像修正过程)
图像修正部26通过将从合成边缘程度R的最小值到最大值放大变换 为从可分别显示在输出部17中的最小灰度值到最大灰度值，而对各像素 进行取得与合成边缘程度对应的灰度值的处理。并且，进行将各像素的像 素值A变换为与该合成边缘程度R对应的灰度值的处理。由此，修正后
的各像素构成修正图像数据。0078
图像处理部13根据操作者的指示将变动量计算部21、边缘程度计算部23及图像修正部26的各运算处理结果输出到输出部17中。0079
由此，根据实施例2相关的X射线拍摄装置具备的图像处理系统，变 动量计算部21对各像素按相邻像素相邻的方向，来个别地计算变动量RX、 RY，边缘程度计算部23对各像素分别基于各方向的变动量RX、 RY来个 别地计算边缘程度RX、 RY。进而，计算出利用这些边缘程度RX、 RY来 合成的合成边缘程度R。根据这种边缘程度RX、 RY或合成边缘程度R， 即使在除像素值A的斜率在纵方向或横方向的情况以外的任意方向，也能 适当地评价该像素为边缘部的可能性。
0080
另外，因为具备根据合成边缘程度R来修正对应的像素的像素值A 的图像修正部26，所以能够生成对被精度良好评价为边缘部的像素进行强 调的图像数据(修正图像数据)。
实施例0081
以下，参照附图对本发明的实施例3进行说明。图9是表示实施例3 相关的X射线拍摄装置的大致构成的框图。且有，对与实施例1同样的构 成附加同一符号并省略详细说明。0082
实施例3相关的X射线拍摄装置进行断层摄影。本装置具备可旋转地 支撑X射线管3和FPD5的C字状臂9。摄影控制部7驱动C字状臂9而 从各种角度对被摄体M进行摄影。0083
另外，图像处理部13具备变动量计算部21、边缘程度计算部23、 再构成处理部27和存储部28。从收集部11输出的输入图像数据被输入到 再构成处理部27中。再构成处理部27基于输入图像数据进行再构成运算 处理，生成三维图像数据。生成的三维图像数据被存储到存储部28中。0084
接着，对实施例3相关的X射线拍摄装置的动作以图像处理部13中的处理为中心进行说明。图io是表示x射线拍摄装置的动作步骤的流程
图。且有，对在实施例1中说明的处理进行简略说明。
0085
〈步骤S1〉进行x射线摄影
由摄影控制部7的控制进行X射线摄影。收集部11将输入图像数据 输出到图像处理部13中。0086
〈步骤S6〉生成三维图像数据
再构成处理部27基于输入图像数据进行再构成运算处理，生成三维 图像数据。作为再构成运算处理，例如举出了利用适当的再构成函数进行 巻积积分并且对巻积积分结果进行逆投影的处理。生成的三维图像数据被 存储在存储部28中。0087
<步骤82>计算变动量(变动量计算过程)
变动量计算部21从存储部28中读取三维图像数据。读取的三维图像 数据作为在正交的3方向上配列了像素的三维图像数据，为了说明的方便 上，将这些3方向分别称为X轴、Y轴和Z轴。0088
变动量计算部21对各像素计算X轴方向的变动量SX。例如，作为变 动量SX(x,y,z)，例示了将像素值A ( x + 1 ， y ， z )与像素值A ( x — 1， y， z)之差用2除后的值。同样地，对于各像素也计算Y轴方向的 变动量SY和Z轴方向的变动量SZ。并且，对各像素分别计算出的变动量 SX、 SY、 SZ被输出到边缘程度计算部23中。0089
〈步骤S3〉计算边缘程度(边缘程度计算步骤)
边缘程度计算部23进行基于与各像素相应的周边像素组g的变动量 SX来求出各像素的边缘程度RX的运算处理。其中，作为周边像素组g， 例示了以该像素为中心的X、 Y、 Z轴方向的9行的公共部分共729个像 素。另外，将边缘程度RX与变动量SX的关系式在公式6,中例示。0090<formula>formula see original document page 22</formula>0091
边缘程度计算部23同样基于各像素的周边像素组g的变动量SY、 SZ 计算各像素的边缘程度RY、 RZ。并且，合成边缘程度RX、 RY、 RZ而 取得合成边缘程度R。在公式7中例示了合成边缘程度R和各边缘程度 RX、 RY、 RZ的关系式。
0092
<formula>formula see original document page 22</formula>
0093
图像处理部13根据操作者指示将再构成处理部27、变动量计算部21、 边缘程度计算部23的各运算处理的结果输出到输出部17中。0094
由此，根据实施例3相关的X射线拍摄装置具备的图像处理系统，通 过具备变动量计算部21和边缘程度计算部23，对三维图像数据也能够适 当地评价该像素为边缘部的可能性。
0095
本发明并不限定于上述实施方式，也可如下进行变形实施。0096
(1)在上述的各实施例中，虽然对各像素计算边缘程度，但并不限 定于此。例如，在操作者通过操作输入部15来指定边缘程度的计算对象 的像素(像素数为单个或多个中的之一即可)时，边缘程度计算部23也 可以按照对指定的像素计算边缘程度的方式进行变更。另外，变动量计算 部21也可以按照对与指定的像素对应的周边像素计算变动量的方式进行变更。
0097
(2) 在上述的各实施例中，虽然作为沿着与输入图像数据、三维图 像数据都正交的2轴或3轴排列像素进行了说明，但是并不限定于此，也 能够使用任意像素排列的图像数据。
0098
(3) 在上述的各实施例中，作为周边像素组g虽然在图6中进行了 例示，但是并不限定于此。例如，如图10所示，对于任意一个像素(x,y) 也可以将位于以该像素为大致中心的椭圆区域Q内的像素作为周边像素 组g(x,y)。由此，可以使从一个位置像素直到各方向最远的像素为止的 距离大致相等。另外，关于各周边像素组的像素数，只要是2以上即可适 当进行变更。
0099
(4) 在上述的各实施例中，虽然具体例示了变动量，但是并不限定 于此。例如，也可以将变动量SX (x,y)作为A(x+l，y)与A(x,y)之差。或 者，也可以将变动量SX (x,y)作为A(x,y)与A(x-l,y)之差。
0100
(5) 在上述的各实施例中，虽然作为变动量的方向，具体例示了横 方向、纵方向或X轴、Y轴、Z轴，但是并不限定于此，也可以变更为对 于倾斜方向等任意方向计算变动量。
0101
(6) 在上述的各实施例中，边缘程度虽然作为"变动量增减幅度" 相对于"变动幅度总和"的比例，但是并不限定于此。例如，为了用百分 比显示，也可以将在上述的比例上乘以IOO后的值作为边缘程度。或者， 也可以将上述的比例的平方作为边缘程度RX'。此时的边缘程度RX'和变 动量SX的关系式在公式8中表示。
010公式8<formula>formula see original document page 24</formula>0103
另外，此时在与同样求出的边缘程度RX'进行合成而取得合成边缘程 度R时，虽然也可以利用上述的公式5，但是也可以利用以下示出的公式 9。
0104
公式9<formula>formula see original document page 24</formula>
(7)在上述的实施例2中，虽然图像修正部26根据合成边缘程度R 进行修正，但是也可以根据边缘程度RX、 RY进行修正。
(8)在上述的各实施例中，输入图像数据是与FPD5的检测结果相应 的图像数据，但是并不限定于此，也可以变更为通过适当的检测单元或拍 摄单元等得到图像数据。另外，输入图像数据是由X射线摄影得到的，但 是并不限定于此，也可以变更为其他的放射线图像数据或光学图像数据。 另外，虽然未特定X射线拍摄装置的用途，但是例如，也可以是用于医用 领域或非破坏检测、RI (Radioisotope)检查等产业领域等中的X射线拍 摄装置。
0105
0106
权利要求
1.一种边缘评价方法，对构成图像数据的像素为边缘部的可能性进行评价，该边缘评价方法，具有变动量计算过程，对各像素计算相对于相邻像素的像素值的变动量；和边缘程度计算过程，对任意一个像素，基于由该一个像素及其周边的像素构成的周边像素组的各变动量，来计算对所述一个像素为边缘部的可能性进行数值化的边缘程度。
2. 根据权利要求1所述的边缘评价方法，其特征在于， 在将所述周边像素组的各变动量置换为绝对值后进行合计的值作为变动幅度总和，且在将直接合计所述周边像素组的各变动量后的值的绝对值作为变 动量增减幅度时，所述边缘程度是所述变动量增减幅度相对于变动幅度总和的比例或 与该比例相应的值。
3. 根据权利要求1或2所述的边缘评价方法，其特征在于， 与一个像素对应的周边像素组是以一个像素为大致中心的区域内的像素的集合。
4. 根据权利要求1 3中任意一项所述的边缘评价方法，其特征在于， 所述相邻像素为相对于各像素，相邻的方向不同的多个像素， 所述变动量计算过程对各像素按所述相邻像素相邻的方向，个别地计算所述变动量，所述边缘程度计算过程对一个像素按与各方向对应的变动量，个别地 计算所述边缘程度。
5. 根据权利要求4所述的边缘评价方法，其特征在于， 所述边缘程度计算过程还对各方向的所述边缘程度进行合成而取得合成边缘程度。
6. 根据权利要求1 5中任意一项所述的边缘评价方法，其特征在于，图像数据是二维图像数据或三维图像数据。
7. —种边缘检测方法，对构成图像数据的像素为边缘部进行检测， 该边缘检测方法，具有变动量计算过程，对各像素计算相对于相邻像素的像素值的变动量； 边缘程度计算过程，对任意一个像素，基于由该一个像素及其周边的像素构成的周边像素组的各变动量，来计算对所述一个像素为边缘部的可能性进行数值化的边缘程度；和边缘检测过程，根据将所述边缘程度与阈值进行比较后的结果来判断所述一个像素为边缘部。
8. —种图像修正方法，对构成图像数据的像素的像素值进行修正， 该图像修正方法，具有变动量计算过程，对各像素计算相对于相邻像素的像素值的变动量； 边缘程度计算过程，对任意一个像素，基于由该一个像素及其周边的像素构成的周边像素组的各变动量，来计算对所述一个像素值为边缘部的可能性进行数值化的边缘程度；和图像修正过程，根据所述边缘程度对所述一个像素的像素值进行修正。
9. 一种图像处理系统，其对构成图像数据的像素为边缘部的可能性进 行评价，该图像处理系统，具备变动量计算单元，其对各像素计算相对于相邻像素的像素值的变动量；禾口边缘程度计算单元，其对任意一个像素，基于由该一个像素及其周边 的像素构成的周边像素组的各变动量来计算对所述一个像素为边缘部的 可能性进行数值化的边缘程度。
10. 根据权利要求9所述的图像处理系统，其特征在于， 在将所述周边像素组的各变动量置换为绝对值后合计的值作为变动幅度总和，且在将直接合计所述周边像素组的各变动量后的值的绝对值作为变 动量增减幅度时，所述边缘程度是变动量增减幅度相对于所述变动幅度总和的比例或 与该比例相应的值。
11. 根据权利要求9或10所述的图像处理系统，其特征在于， 与一个像素对应的周边像素组是以一个像素为大致中心的区域内的像素的集合。
12. 根据权利要求9 11中任意一项所述的图像处理系统，其特征在于，所述相邻像素作相对于各像素相邻的方向不同的多个像素， 所述变动量计算单元对各像素按每个所述相邻像素相邻的方向，来分别计算所述变动量，所述边缘程度计算单元对一个像素按与各方向对应的变动量，来分别计算所述边缘程度。
13. 根据权利要求12所述的图像处理系统，其特征在于， 所述边缘程度计算过程还对各方向的所述边缘程度进行合成而取得合成边缘程度。
14. 根据权利要求9 13中任意一项所述的图像处理系统，其特征在于，图像数据是二维图像数据或三维图像数据。
15. 根据权利要求9 14中任意一项所述的图像处理系统，其特征在于，所述系统还具备边缘检测单元，该边缘检测单元根据将所述边缘程度 与阈值进行比较后的结果来判断所述一个像素为边缘部。
16. 根据权利要求9 15中任意一项所述的图像处理系统，其特征在于，所述系统还具备图像修正单元，该图像修正单元根据所述边缘程度对 所述一个像素的像素值进行修正。
17. 根据权利要求9 16中任意一项所述的图像处理系统，其特征在于，所述边缘程度计算单元对构成所述图像数据的各像素计算所述边缘 程度。
18. 根据权利要求9 16中任意一项所述的图像处理系统，其特征在于，所述系统还具备输入单元，该输入单元对作为所述边缘程度的计算对 象的像素进行指定，所述边缘程度计算单元对由所述输入单元指定的像素计算所述边缘 程度。
19. 根据权利要求9 18中任意一项所述的图像处理系统，其特征在于，所述系统还具备输出单元，该输出单元输出所述边缘程度。
20. 根据权利要求9 19中任意一项所述的图像处理系统，其特征在于，该图像处理系统构成为能从外部输入图像数据。
全文摘要
本发明的图像处理系统具有对各像素计算相对于相邻像素的变动量的像素值的变动量计算部和对任意一个像素基于由该一个像素及其周边的像素构成的周边像素组的各变动量来计算对一个像素为边缘部的可能性进行数值化的边缘程度的边缘程度计算部。根据由此构成的图像处理系统，即使像素接受了噪声也能够高精度评价是否为边缘部。
文档编号G06T7/60GK101617340SQ20078005192
公开日2009年12月30日 申请日期2007年3月6日 优先权日2007年3月6日
发明者冈村升一 申请人:株式会社岛津制作所