CN102448374B

CN102448374B - 图像处理装置、x线摄影装置及图像处理方法

Info

Publication number: CN102448374B
Application number: CN201080022803.7A
Authority: CN
Inventors: 胜间田真弥; 铃木克己
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Ltd
Priority date: 2009-05-27
Filing date: 2010-05-17
Publication date: 2014-05-14
Anticipated expiration: 2030-05-17
Also published as: JP5670327B2; JPWO2010137482A1; CN102448374A; WO2010137482A1

Abstract

提供一种可抑制作为接合了多个X线图像的接合图像的长跨距图像的浓度和对比度的不均的图像处理装置，其中，图像处理装置6具有：图像存储器62，存储部分重复的多个X线图像(20、21、22)的数据；浓度对比度校正部64，进行使多个X线图像的浓度和对比度对应的处理；长跨距图像取得部65，接合通过浓度对比度校正部64对应了浓度和对比度的所有X线图像(20、21’、22”)，取得和被测体M的长跨距摄影区域对应的长跨距X线图像，浓度对比度校正部64具有：第1表格做成部642，其根据多个X线图像中相邻的二个X线图像(20、21或21、22)的重复部分(20a、21a、21b、22a)的图像信息，取得可使二个X线图像(20、21或21、22)的浓度和对比度对应的灰度变换表格(AB、BC)；第2表格做成部643，根据灰度变换表格(AB、BC)，取得可使从多个X线图像(20～22)选择的基准X线图像(20)和该基准X线图像以外的剩余X线图像(21、22)的浓度和对比度对应的第2灰度变换表格(BA、CB)；平均化处理部644，根据第2灰度变换表格(BA、CB)，进行使剩余X线图像(21、22)的浓度和对比度对应的处理。

Description

图像处理装置、X线摄影装置及图像处理方法

技术领域

本发明涉及一种图像处理装置、X线摄影装置及图像处理方法。

背景技术

使用具有部分重叠的多个X线图像取得作为接合图像的长跨距图像时，在使接合前的各图像的显示浓度对应后结合并取得长跨距图像，这一技术为世人所知。例如专利文献1公开的技术中，首先对接合前的各图像进行对数变换。接着根据对数变换后的各图像取得各图像的柱状图。接着检测出取得的各柱状图的重心位置的偏移，修正该偏移。并且结合偏移修正后的各图像，取得作为接合图像的长跨距图像。

专利文献1：特开2007-275228号公报

发明内容

但是，现有的方法是这样的技术：为了使应取得的接合图像整体的显示浓度一致，将接合前的各图像的浓度调整为规定值，而对于接合图像整体的对比度如何并未关注。

一般情况下，通过X线摄影获得的X线图像的对比度取决于摄影时使用的X线的线质(能量大小)。因此，在拍摄具有部分重叠的多个图像时，各摄影中使用改变了线质的X线时，虽可通过现有方法调整接合前的各图像的浓度，但浓度调整后的各图像之间，对比度有时会不同。各图像之间对比度不同时，对于某个骨骼部分，会产生在某个图像中可明确确认，但在其他图像中无法确认的问题。

本发明要解决的课题是，提供一种至少可抑制作为接合了多个X线图像的接合图像的长跨距图像的对比度的不均的图像处理装置、及包括该装置的X线摄影装置。

本发明通过以下解决方式解决上述课题。并且在下述解决方式中，对表示发明的实施方式的附图附加对应的标记进行说明，但该标记仅用于容易地理解发明，不用作限定发明。

发明涉及的图像处理装置(6)具有：存储单元(62)，存储部分重复的多个X线图像(20、21、22)的数据；长跨距图像取得部(65)，接合多个X线图像，取得和被测体(M)的长跨距摄影区域对应的长跨距X线图像。

除此之外还具有校正单元(64)，其进行使多个X线图像(20、21、22)的对比度(优选浓度和对比度)对应的处理，长跨距图像取得单元(65)接合通过校正单元(64)使对比度(优选浓度和对比度)对应的X线图像(20、21’、22”)，取得长跨距X线图像。校正单元(64)具有：表格取得单元(642)，根据多个X线图像中相邻的二个X线图像(20、21或21、22)的重复部分(20a、21a、21b、22a)的图像信息，取得可使二个X线图像(20、21或21、22)的对比度(优选浓度和对比度)对应的灰度变换表格(AB、BC)；平均化单元(644)，根据灰度变换表格(AB、BC)，进行使多个X线图像(20～22)的对比度(优选浓度和对比度)对应的处理。

发明涉及的X线摄影装置(100)具有：图像摄影单元(1，4)，取得部分重复的多个X线图像(20～22)的数据；图像处理单元(6)，对多个X线图像进行接合处理；显示单元(7)，显示图像处理单元(6)的处理结果。图像处理单元(6)具有上述发明涉及的图像处理装置(6)。显示单元(7)显示通过长跨距图像取得单元(65)取得的长跨距X线图像。

根据上述发明，可抑制作为接合图像的长跨距图像的对比度的不均。结果可取得易于诊断的长跨距X线图像，有助于提高诊断效率。

并且，现有技术中，被测体厚度较厚的部位包含于长跨距摄影区域时，因摄影条件的制约，对该厚度较厚的部位的摄影需要牺牲画质，但根据上述发明，可改变X线的线质来拍摄接合对象的X线图像，因此可避免牺牲画质，并且从被测体(患者)遭受辐射的角度而言也较为有利。

附图说明

图1是表示本发明的一个实施方式涉及的X线摄影装置的整体构成的框图。

图2(a)～图2(c)是表示作为长跨距X线图像的基础的接合对象的X线图像的一例的图。

图3是表示图1的图像处理部的动作顺序的流程图。

图4(a)是图2(a)所示的图像20的重复区域20a中的柱状图，图4(b)是图2(b)所示的图像21的重复区域21a中的柱状图。

图5(a)将累积了图4(a)的柱状图的累积柱状图比率化，图5(b)将累积了图4(b)的柱状图的累积柱状图比率化。

图6(a)及图6(b)是表示通过图3的步骤323生成的2个图像间灰度变换表格的一例的图。

图7是表示通过图3的步骤332生成的标准化灰度变换表格的一例的图。

图8是表示通过图3的步骤332生成的标准化灰度变换表格的一例的图。

图9是表示图3的步骤34的浓度对比度平均化处理的一例的图。

附图标记

100X线摄影装置、1X线管(图像摄影单元)、4X线检测器(图像摄影单元)、M被测体、6图像处理部(图像处理装置)、61图像取得部、62图像存储器(存储单元)、63接合位置设定部(接合位置设定单元)、64浓度对比度校正部(校正单元)、641对数变换处理部(对数变换处理单元)、642第1表格做成部(表格取得单元)、642a柱状图取得部(柱状图取得单元)、642b特征量取得部(特征量取得单元)、642c表格生成部、643第2表格做成部(第2表格取得单元)、643a基准图像选择部(基准图像选择单元)、643b表格生成部、644平均化处理部(平均化单元)、65长跨距图像生成部(长跨距图像取得单元)、7显示部(显示单元)、20～22X线图像、20a、21a、21b、22a重复区域、21’、22’灰度变换后的图像、22”灰度再变换后的图像、AB、BC 2个图像间灰度变换表格(灰度变换表格)、BA、CB反向搜索表格(第2灰度变换表格)

具体实施方式

以下参照附图说明上述发明的一个实施方式。

图1所示的本实施方式的X线摄影装置100例如医院等中用于诊断的装置，具有X线照射用的X线管1(图像摄影单元)、及作为透过X线图像检测用的二维X线检测器的平板型等的X线检测器4。该X线管1及检测器4彼此同步，可沿被测体M的体轴方向Z移动地构成。在使被测体M静止的状态下，使X线管1及检测器4向被测体M的体轴方向Z移动的同时反复进行摄影，从而可取得沿着被测体M的体轴的连续的、且部分重复的多个X线图像。

检测器4(图像摄影单元)在每次摄影时检测从X线管1照射、并通过被测体M的透过X线图像，将其检测结果变换为电气信号(X线检测信号)，输出到图像处理部6。

图像处理部6(图像处理装置)由CPU、ROM、RAM等构成，根据通过从检测器4输出的X线检测信号逐个取得的多个X线图像，进行用于显示的图像处理(灰度变换、接合)，生成和被测体M的长跨距摄影区域(例如从腹部到下肢为止的区域)所对应的一个长跨距X线图像。

显示部7(显示单元)接受图像处理部6的输出，显示长跨距X线图像。

此外也可如下构成：通过图像处理部6生成的长跨距X线图像在数字图像数据的状态下，保存到光磁盘装置等记录装置中，借助网络传送到外部设备。

接着根据以框图示出图像处理部6的构成的图1及图3的流程图，说明图像处理部6的动作(S1～S4)。

(S1)

首先，在图3所示的步骤(以下简称为“S”)1中，接合对象的多个X线图像的数据从X线检测器4提供到图像处理部6(接合对象的图像数据的供给)。提供到图像处理部6的图像数据是，从X线检测器4输出的X线检测信号借助A/D变换器(省略图示)变换为数字图像数据。

在本实施方式中，例如提供图2(a)～图2(c)所示的、沿着被测体M的体轴的连续的、且部分重复的3个X线图像20～22。并且，图像20和图像21的一部分重复，图像21和图像22的一部分重复。提供的图像数不限于3个，也可是2个或4个以上。

此外，在此设接合对象的X线图像的个数为N、对它们按照摄影顺序分配1到N为止的识别代码。在本实施方式中，接合对象的图像个数是3个，因此按照摄影顺序分配1到3的识别代码。识别代码是，图2(a)所示的图像20是1，图2(b)所示的图像21是2，图2(c)所示的图像22是3。

返回到图1及图3。提供到图像处理部6的图像20～22通过图像取得部61取得，暂时存储到图像存储器62(存储单元)中。

(S2)

接着在S2中，存储到存储器62的3个X线图像的数据读出到接合位置设定部63(接合位置设定单元)，在此设定相邻的图像各自的接合位置(接合位置的设定)。

在本实施方式中，设定作为上图像的图像20和作为中图像的图像21的接合位置，并且设定图像21和作为下图像的图像22的接合位置。接合位置的设定可根据用户的输入来确定，或者由公知的处理(例如特开2008-67916号公报等公开的接合部的判断处理)来决定。

此外，图2(a)及图2(b)中的标记20a、21a是指，在由位置设定部63设定的接合位置上接合了图像20、21的情况下的图像20、21的重复区域。图2(b)及图2(c)中的标记21b、22a也同样，是指在由位置设定部63设定的接合位置上接合了图像21、22时的图像21、22的重复区域。

(S3)

通过图像处理部6的处理获得的长跨距X线图像为了在整体上看上去象一个图像，需要平滑地连接接合部。为了平滑地连接，需要按照各图像调整灰度(浓度)并使显示浓度对应的作业。这是因为，如果获得的长跨距X线图像的接合边界线的内和外有浓度差，则难以进行影像判读，诊断能力可能下降。

但是，仅使显示浓度对应还不足够，调整接合前的各图像的对比度、与显示浓度一起也使显示对比度对应，这一点很重要。即使浓度相同但对比度不同时，对于某个骨骼部分，会产生在某个图像中可明确确认，但在其他图像中无法确认的可能性，变成难以诊断的长跨距X线图像。

因此，在本实施方式中，作为本发明的一个实施方式，为了抑制接合前的多个X线图像的显示浓度和对比度的不均，在图像处理部6上设置浓度对比度校正部64，在此进行通过S2设定了接合位置的图像间的浓度对比度的不均抑制的处理(浓度对比度校正处理)(S3)。此外，在本实施方式中，以下说明的校正浓度和对比度两者的方式，定位为校正对比度的本发明的一个方式。

本实施方式的浓度对比度校正处理基本上进行以下处理：以选定的基准图像为起点，使与之相邻的图像的显示浓度和对比度与基准图像的显示浓度和对比度对应。之后，将浓度和对比度对应的各图像在S2中设定的接合位置上接合，生成长跨距X线图像。

以下说明浓度对比度校正处理(S31～S34)的一例。

(S31)

首先，当设定了接合位置的各图像20～22提供到浓度对比度校正部64(浓度对比度校正单元)时，在S31中，通过对数变换处理部641(对数变换处理单元)对各图像20～22进行对数变换处理。

一般情况下，人体摄影中使用的X线的光子能量带中的、透过的X线量(入射到X线检测器的线量)，具有相对透过的物质的厚度和衰减系数对数性减少的特性。因此，通过进行对数变换处理，可抵消来自X线检测器的信息的对数性减少倾向，简单地置换为透过的物质的厚度和衰减系数的加算信息。这是进行对数变换的原因。

(S32)

接着在S32中，通过第1表格做成部642(表格取得单元)做成使相邻的X线图像20、21中的一个与另一个的图像浓度及对比度一致的2个图像间灰度变换表格AB(灰度变换表格)。同时在本实施方式中，也做成相邻的X线图像21、22的2个图像间灰度变换表格BC(灰度变换表格)。

以下说明用于做成表格AB(参照图6(a))的方法的一例(S321～S323)。此外，表格BC(参照图6(b))也可通过同样的步骤做成，因此省略其说明。

(S321)

首先在S321中，通过柱状图取得部642a(柱状图取得单元)取得柱状图，该柱状图表示进行了对数变换的图像20、21的重复区域20a、21a中的像素信号的强度(浓度、横轴)和出现频率(频率、纵轴)的对应关系。图4(a)及图4(b)表示在此导出的柱状图A、B的一例。

根据图4(a)及图4(b)的柱状图A、B可理解以下内容。首先，图像21的区域21a的图像浓度大于图像20的区域20a。这是因为，柱状图B的浓度范围(波形开始位置(最小值)到结束位置(最大值)的范围)和柱状图A相比靠右，即位于高浓度一侧。其次，图像20的区域20a的对比度大于图像21的区域21a。这是因为，柱状图A的浓度范围的幅度W1和柱状图B的浓度范围的幅度W2相比较大。由此可知，图像20和图像21相比，图像浓度较低，但对比度较大。浓度和对比度的校正是使图像20、21的柱状图一致的处理。

(S322)

返回到图1及图3。接着在S322中，通过特征量提取部642b(特征量提取单元)，从由取得部642a导出的各柱状图A、B提取特征量A’、B’。

作为提取的特征量A’、B’的一例包括柱状图的最小值、面积平均值、最大值、中间值、最频值中的任意一个或它们的组合。并且，求出累积了各柱状图的累积柱状图，使该累积柱状图比率化(例如参照图5(a)及图5(b))，可将与多个比率对应的浓度作为特征量A’、B’。其中，使累积柱状图比率化是指，以浓度值(设图像20的浓度为A1、A2……An、图像21的浓度为B1、B2……Bn)为横轴，以累积了该浓度值以下的浓度的出现频率(频度)的值的比率(图像20、21均为R1、R2……Rn)为纵轴来表示的图。

在本实施方式中，从各柱状图A、B提取特征量A’、B’时，优选组合多个上述特征量(例如最大值和最大值等)，提取该组合特征量。这样一来，有利于做成较高精度的灰度变换表格。

(S323)

接着在S323中，通过表格生成部642c，根据由提取部642b提取的各柱状图A、B的特征量A’、B’，生成使通过处理部641进行了对数变换的图像20、21中的一个与另一个的图像浓度及对比度一致的2个图像间灰度变换表格AB。

根据从2个图像20、21获得的特征量A’、B’，可知任意一个图像的特征量与任意另一图像的特征量对应，但不清楚特征量以外的像素值如何对应。因此，需要根据多个特征量来内插特征量之间的像素值，做成将一个图像的像素值可变换为另一个图像的像素值的灰度变换表格。

在本实施方式中，通过上述特征量做成表格AB时，导出多个特征量的逼近直线，根据该导出的逼近直线进行灰度变换表格的做成。由此做成图6(a)所示的、基于特征量的2个图像间灰度变换表格AB。逼近直线式的导出如仅是2点则可使用联立方程，当表格做成中使用的特征量个数不确定时，可使用最小乘方法等导出。在此使用逼近直线的理由是，在普通医用诊断中使用的X线的线质(管电压30kV～200kV)的范围内，当对象是人体组织时，光电效果和康普顿效应的和的减弱系数可通过指数函数逼近，因此在使用对数变换后的图像时，指数函数部分抵消，像素值与衰减系数和厚度的和成比例，因此可逼近一次式的直线(在X线摄影方法中称为指数函数法则)。

其中，在表格做成中，不是用一个逼近直线内插所有像素值，而是内插重复区域20a、21a的像素值的范围以外时，优选使表格以低像素区域、重复区域的像素区域、及高像素区域使用其他直线逼近来插值。具体而言，在低像素区域中，为进行表格的低像素区域的插值，仅使用接近低像素区域的特征量(接近最低值、平均值、柱状图累积率的低像素一侧的数个百分比等的特征量)，求出逼近直线。在重复区域的像素区域中，为进行插值，使用全部特征量、或除了靠近低像素区域和高像素区域的特征量的特征量，求出逼近直线。

在高像素区域中，为进行插值，使用接近高像素区域的特征量(最高值、平均值、柱状图累积率的高像素值一侧的数个百分比等的特征量)，求出逼近直线。并且，将各区域作为以各区域使用的所求出的逼近直线来内插并连接的表格做成。此时，为使各区域的连接部分的表格的倾斜不剧烈变化，优选使用平滑化处理(移动平均等)。

(S33)

接着在S33中，通过第2表格做成部643(第2表格取得单元)，对作为浓度对比度校正处理的基准的特定的X线图像，做成与其他图像的图像浓度和对比度一致的标准化灰度变换表格ABC(第2灰度变换表格)。

以下说明做成标准化灰度变换表格ABC的方法的一例(S331、S332)。

(S331)

首先，在S331中，通过基准图像选择部643a(基准图像选择单元)，选择作为浓度对比度校正处理的起点的基准图像。选择的基准图像只要是已经存在的图像20、21、22中的任意一个即可，或者也可选择对图像20、21、22进行了浓度对比度校正的其他图像。基准图像的选择基准没有特别限定，例如可根据通过做成部642做成的2个图像间灰度变换表格AB(参照图6(a)、BC(参照图6(b))，将具有最小值或最大值的图像选择为基准图像，或者将和基准图像对应时的变换表格的倾斜的平均最小(对比度高)的图像选择为基准图像。以下示例将图像20作为基准图像选择的情况。

(S332)

接着在S332中，通过表格生成部643b，生成使通过选择部643a选择的作为基准图像的图像20、与作为其他图像的图像21、22的图像浓度和对比度一致的标准化灰度变换表格。

在本实施方式中，作为第1阶段，首先做成相对图像20的图像21的标准化灰度变换表格。因图像20是基准图像，所以可做成通过S32做成的2个图像间灰度变换表格AB(参照图6(a))的反向搜索表格BA(参照图7)，并使用它。

作为第2阶段，做成相对图像20的图像22的标准化灰度表格。做成通过S32做成的2个图像间灰度变换表格BC(参照图6(b))的反向搜索表格CB(参照图8)，做成它与上述反向搜索表格BA的合成反向搜索表格CBA(省略图示)，并使用它。

(S34)

接着在S34中，通过平均化处理部644(平均化单元)，根据做成的标准化灰度变换表格，进行作为基准图像的图像20以外的图像21、22的平均化处理。

例如如图9所示，首先将图像21的全部像素值变换为基于通过S332做成的反向搜索表格BA的灰度，作为图像21’。由此使图像20和与其下游侧相邻的灰度变换后的图像21’的浓度及对比度一致(图像21→图像21’)。接着，将图像22的全部像素值变换为基于通过S332做成的合成反向搜索表格CBA的灰度，作为图像22”。这样一来，使灰度变换后的图像21’、及与其下游侧相邻的灰度再变换后的图像22”的浓度及对比度一致(图像22→图像22”)。

此外，图像22→图像22”的灰度变换也可通过2个阶段进行。

例如，首先将图像22的全部像素值变换为基于通过S332做成的反向搜索表格CB的灰度，作为图像22’。这样一来，使灰度变换前的图像21、及与其下游侧相邻的灰度变换后的图像22’的浓度及对比度一致(图像22→图像22’)。接着，将图像22’的全部像素值变换为基于通过S332做成的反向搜索表格BA的灰度，作为图像22”(图像22’→图像22”)。

经过以上S31～S34，可获得浓度和对比度对应的图像20、21’、22”(参照图9)。

(S4)

接着在S4中，通过长跨距图像生成部65(长跨距图像取得单元)，在通过S2设定的接合位置上接合浓度和对比度对应的图像20、21’、22”，做成作为一个图像的长跨距图像。

各图像的接合如公知方法所示，可在对重复区域20a、21a、21b、22a的任意区域加权并校正后接合，或者可使用相邻的图像20、21’或21’、22”的一个图像。

此外，一般情况下，做成的长跨距图像输出到图像输出装置，但该输出装置可以是监视器、胶片打印机，也可作为数据保存到服务器、记录介质。

在本实施方式中，将做成的长跨距图像输出到显示部7并显示，结束流程。

根据本实施方式，在接合多个X线图像20～22前，使各图像的显示浓度和对比度对应，因此接合后的长跨距X线图像不会产生浓度及对比度的不均。其结果是，可取得易诊断的长跨距X线图像。

以上说明的实施方式用于容易地理解上述发明，不用于限定上述发明。因此，上述实施方式中公开的各要素也包括属于上述发明的技术范围内的所有设计变更、均等物。

Claims

1.一种图像处理装置，具有：存储单元，存储部分重复的多个X线图像的数据；长跨距图像取得单元，接合上述多个X线图像，取得和被测体的长跨距摄影区域对应的长跨距X线图像，其特征在于，

具有校正单元，其进行使上述多个X线图像的对比度对应的处理，上述校正单元具有：

表格取得单元，其具有：对数变换处理单元，对上述多个X线图像中相邻的二个X线图像的重复部分进行对数变换处理；柱状图取得单元，从经对数变换的二个X线图像的重复部分的图像，取得表示像素浓度和出现率的对应关系的柱状图；特征量提取单元，从上述柱状图提取至少一个特征量，上述表格取得单元根据上述特征量，生成使上述二个X线图像的对比度及浓度对应时使用的灰度变换表格；

平均化单元，根据上述灰度变换表格，进行使上述多个X线图像的对比度及浓度对应的处理，

上述长跨距图像取得单元接合通过上述校正单元使对比度及浓度对应的X线图像，取得上述长跨距X线图像。

2.根据权利要求1所述的图像处理装置，其特征在于，

上述特征量提取单元提取上述柱状图的最小值、面积平均值、中间值、最大值、最频值中的任意一个以上，作为上述特征量。

3.根据权利要求1所述的图像处理装置，其特征在于，

上述特征量由上述柱状图的累积值的比率提取。

4.根据权利要求1所述的图像处理装置，其特征在于，

上述表格取得单元在取得上述灰度变换表格时，对上述重复部分以外的图像内的低像素区域和高像素区域，使用利用最小二乘法算出的逼近直线进行线性插值。

5.根据权利要求1所述的图像处理装置，其特征在于，

上述校正单元具有基准X线图像选择单元，从上述多个X线图像选择基准X线图像，上述平均化单元进行使上述基准X线图像以外的图像的对比度与上述基准X线图像的对比度对应的处理。

6.根据权利要求5所述的图像处理装置，其特征在于，

上述校正单元具有第2表格取得单元，其根据上述灰度变换表格，取得可使上述基准X线图像和上述基准X线图像以外的图像的对比度对应的第2灰度变换表格，上述平均化单元根据上述第2灰度变换表格，进行使上述基准X线图像以外的图像的对比度和上述基准X线图像的对比度对应的处理。

7.一种X线摄影装置，具有：图像摄影单元，取得部分重复的多个X线图像的数据；图像处理单元，对上述多个X线图像进行接合处理；显示单元，显示上述图像处理单元的处理结果，其特征在于，

上述图像处理单元具有：

存储单元，存储上述多个X线图像的数据；

校正单元，进行使上述多个X线图像的对比度对应的处理；

长跨距图像取得单元，接合通过上述校正单元使对比度对应的X线图像，取得和被测体的长跨距摄影区域对应的长跨距X线图像，

上述校正单元具有：

上述显示单元显示通过上述长跨距图像取得单元取得的上述长跨距X线图像。

8.一种图像处理方法，具有：存储步骤，存储部分重复的多个X线图像的数据；长跨距图像取得步骤，接合上述多个X线图像，取得和被测体的长跨距摄影区域对应的长跨距X线图像，其特征在于，

具有校正步骤，其进行使上述多个X线图像的对比度对应的处理，上述校正步骤具有：

表格取得步骤，其具有：对数变换处理步骤，对上述多个X线图像中相邻的二个X线图像的重复部分进行对数变换处理；柱状图取得步骤，从经对数变换的二个X线图像的重复部分的图像，取得表示像素浓度和出现率的对应关系的柱状图；特征量提取步骤，从上述柱状图提取至少一个特征量，上述表格取得步骤根据上述特征量，生成使上述二个X线图像的对比度及浓度对应时使用的灰度变换表格；

平均化步骤，根据上述灰度变换表格，进行使上述多个X线图像的对比度及浓度对应的处理，

上述长跨距图像取得步骤接合通过上述校正步骤使对比度及浓度对应的X线图像，取得上述长跨距X线图像。