CN103126704B - 图像数据校正设备和方法 - Google Patents

Abstract

本发明的实施例涉及成像系统和方法，其包括：获取图像检测单元阵列中待校正图像检测单元的相邻图像检测单元的图像数据；基于待校正图像检测单元的相邻图像检测单元的图像数据对所述待校正图像检测单元进行插值，得到所述待校正图像检测单元的估计图像数据；判断所述待校正图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点；基于所述判断结果选择不同的校正模式对所述估计图像数据进行校正。

Description

图像数据校正设备和方法

技术领域

本发明的实施例一般涉及成像系统和方法，具体涉及用于降低成像系统中图像检测单元阵列中故障图像检测单元对图像质量影响的设备和方法。

背景技术

医学诊断成像系统包括各种成像形式，诸如平面x射线、超声、磁共振(MR)、电子束x线断层照相术(EBT)、正电子发射计算机断层照相(PET)、单光子发射计算机断层照相(SPECT)、微计算机断层照相(micro CT)、宏计算机断层照相(macro CT)成像系统等。医学诊断成像系统通过将诸如患者的对象暴露到能量源(诸如通过患者的x射线)来生成对象的图像。所生成的图像可以用于许多目的。例如，可以检测对象的内部缺陷。此外，可以确定对准或内部结构中的变化。还可以表现对象中的流体流。而且，图像可以显示对象中部件的存在或不存在。从医学诊断成像得到的信息应用于许多领域，包括医疗和制造领域。

典型的成像系统采用图像检测单元阵列来检测对象，然后重构和显示所检测的图像。上述阵列包括多个图像检测单元行(row)。每个图像检测单元行包括多个图像检测单元(cell)，每个图像检测单元连接到不同的数据采集系统(DAS)通道(channel)。也就是说，DAS通道可以映射到图像检测单元。每个图像检测单元产生信号。大容量的阵列包括大数量的图像检测单元和DAS通道。由于图像检测单元和DAS通道的数量增加，图像检测单元、DAS通道故障、或DAS专用集成电路(ASIC)故障的可能性增加。此外，由于图像检测单元和DAS通道的数量增加，希望图像检测单元阵列和成像系统组件变得更加集成。

在图像检测单元阵列或DAS通道中任何一个图像检测单元中的问题(problem)可能引起重构图像中的伪影(artifact)。图像检测单元阵列和/或DAS通道中有问题的图像检测单元称为故障图像检测单元(malfunctioning cell)。故障图像检测单元可能以几种不同方式故障，诸如间歇地故障，引起低于其它图像检测单元信号准确性某一百分比的信号，产生弱于相邻图像检测单元信号某一百分比的信号，以及完全不起作用。由成像系统或成像系统中图像检测单元产生的图像中的任何错误或“伪影”可能导致医生、医学从业者或其它观察者基于错误信息采取行动。

增加图像检测单元阵列的容量覆盖允许用户对更大对象更快地进行成像，因为更大阵列的一次扫掠(sweep)对对象的更大部分进行成像。增加图像检测单元阵列的容量覆盖还更准确地对对象进行成像，因为在采用较少次数对对象进行成像时在成像期间消逝较少的时间。随着对更大容量覆盖的追求，图像检测单元阵列和DAS通道数量迅速增加。结果，增加故障图像检测单元的可能性。替换故障图像检测单元较大地增加了系统的成本。替换具有大量图像检测单元阵列和DAS通道的系统上的所有故障组件可能不经济。另外，故障组件中断了医院里的工作流程。因此，非常希望由于故障图像检测单元导致的图像质量恶化或图像伪影最小化的系统。

发明内容

根据本发明的实施例的一个方面，提供一种用于处理图像数据的方法，包括：图像数据获取步骤，获取图像检测单元阵列中待校正图像检测单元的相邻图像检测单元的图像数据；图像数据估计步骤，基于所述相邻图像检测单元的所述图像数据中的与所述待校正图像检测单元相同行的相同行相邻图像检测单元的图像数据对所述待校正图像检测单元进行插值，得到所述待校正图像检测单元的估计图像数据；以及图像数据修正步骤，基于所述相邻图像检测单元的图像数据的特性对所述估计图像数据进行修正。

根据本发明的实施例的另一个方面，提供一种用于处理图像数据的设备，包括：图像数据获取装置，获取图像检测单元阵列中待校正图像检测单元的相邻图像检测单元的图像数据；图像数据估计装置，基于所述相邻图像检测单元的所述图像数据中的与所述待校正图像检测单元相同行的相同行相邻图像检测单元的图像数据对所述待校正图像检测单元进行插值，得到所述待校正图像检测单元的估计图像数据；以及图像数据修正装置，基于所述相邻图像检测单元的图像数据的特性对所述估计图像数据进行修正。

进一步地，按照本发明的实施例的上述方面的方法和设备，还包括：判断所述待校正图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点；以及基于所述判断结果对所述估计图像数据进行校正。

更进一步地，在所述相邻图像检测单元的图像数据中存在图像数据的阶跃点时，按照第一校正模式对所述估计图像数据进行校正。

更进一步地，所述第一校正模式包括：沿与所述图像检测单元阵列的行不同的方向选择所述待校正图像检测单元的相邻图像检测单元，所述方向反映了所述相邻图像检测单元的图像数据中阶跃点的分布方向；计算所选相邻图像检测单元的图像数据的方向偏差Δ，以及基于所述方向偏差Δ来校正所述估计图像数据，得到方向性校正后的图像数据。

更进一步地，计算所述待校正图像检测单元的相邻图像检测单元的图像数据的方向偏差Δ包括：计算所选相邻图像检测单元的图像数据与所述待校正图像检测单元的估计图像数据的差值以及计算所述差值的平均值作为所述方向偏差Δ。

更进一步地，所述第一校正模式包括：对所述待校正图像检测单元进行单边外插，得到所述待校正图像检测单元的单边外插后的图像数据。

更进一步地，对所述待校正图像检测单元进行单边外插包括：计算所述待校正图像检测单元在所述图像检测单元阵列中的左侧权重w_left和右侧权重w_right，以及根据所述左侧权重w_left和右侧权重w_right来计算单边外插后的图像数据X：x＝w_left*x_left+w_right*x_right，其中，x_left和x_right分别是所述待校正图像检测单元左侧相邻图像检测单元和右侧相邻图像检测单元的图像数据。

更进一步地，按照以下方式计算所述左侧权重w_left和所述右侧权重w_right：

针对所述待校正图像检测单元，上一行左侧外插估计上一行的外插值x_{up_left}，上一行右侧外插估计上一行的外插值x_{up_right}，下一行左侧外插估计上一行的外插值x_{down_left}，下一行右侧外插估计上一行的外插值x_{down_right}，其中，所述外插值是采用相应行的相应侧的相邻图像检测单元的图像数据进行拉格朗日插值，以及

通过所述外插值x_{up_left}、x_{up_right}、x_{down_left}、x_{down_right}和所述待校正图像检测单元的上相邻图像检测单元和下相邻图像检测单元的图像数据x_up、x_down来求解以下方程而获得左外插权重w_{up_left}、w_{down_left}和右外插权重w_{up_right}、w_{down_right}：

x_up＝w_{up_left}*x_{up_left}+w_{up_right}*x_{up_right}

x_down＝w_{down_left}*x_{down_left}+w_{down_right}*x_{down_right}

以及

通过所获得的左外插权重w_{up_left}、w_{down_right}、和右外插权重w_{up_right}、w_{down_right}来计算所述左侧权重w_left和所述右侧权重w_right：

w_left＝a*w_{up_left}+(1-a)*w_{down_left}

w_right＝bw_{up_right}+(1-b)*w_{down_right}

其中，a和b是取值在0至1之间的常数。

更进一步地，将所述方向性校正后的图像数据作为所述待校正图像检测单元的校正后的图像数据。

更进一步地，采用所述单边外插后的图像数据作为所述待校正图像检测单元的校正后的图像数据；或者计算所述待校正图像检测单元的所选相邻图像检测单元的图像数据的方向偏差Δ，在所述方向偏差Δ大于0时，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最大者作为所述待校正图像检测单元的校正后的图像数据，否则，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最小者作为所述待校正图像检测单元的校正后的图像数据。

可选地，在所述相邻图像检测单元的图像数据中不存在图像数据的阶跃点时，按照第二校正模式对所述估计图像数据进行校正。

进一步地，所述第二校正模式包括：计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值，基于所述偏差的平均值对所述估计图像数据进行校正，得到所述待校正图像检测单元的校正后的图像数据。

进一步地，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

进一步地，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

进一步地，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：判断图像检测单元在所述图像检测单元阵列中列方向的尺寸与行方向的尺寸之比R是否大于或等于预定阈值V_threshold，如果R大于或等于所述预定阈值V_threshold，计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值，否则，计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

进一步地，其中所述所述预定阈值V_threshold是介于0.5至2之间的数值。

进一步地，获得所述待校正图像检测单元的估计图像数据包括：基于所述待校正图像检测单元所在行的相邻图像检测单元的图像数据对所述待校正图像检测单元进行拉格朗日插值，得到所述估计图像数据。

以下结合附图来描述本发明的实施例。

附图说明

在附图中：

图1是可以应用本发明的实施例的成像系统的示图。

图2是图1中所示的系统的示意框图。

图3是说明包含故障图像检测单元和无故障图像检测单元的图像检测单元阵列的示图。

图4是说明按照本发明的一个实施例的用于校正故障图像检测单元的方法的流程图。

图5a和5b是说明图像数据的阶跃点的示图。

图6是说明按照本发明的一个实施例的故障图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点的判断方法的流程图。

图7是说明图像数据的阶跃点的方向性的示图。

图8是用于说明单边外插的示图。

图9是说明按照本发明的一个优选实施例的故障图像检测单元的图像数据校正方法的流程图。

图10是说明按照本发明的一个实施例的用于校正故障图像检测单元的设备的方框图。

图11示出了按照本发明的实施例进行图像数据校正前后的图像比较。

具体实施方式

本发明的一些实施例提供用于对图像检测单元、DAS通道、和/或DAS ASIC中的错误进行纠错的系统和方法。

为了说明起见，以下详细描述涉及CT成像系统的一些实施例。本领域技术人员将意识到，本发明的实施例同样适合其它成像系统(诸如平面x射线、超声、磁共振(MR)、电子束x线断层照相术(EBT)、正电子发射计算机断层照相(PET)、单光子发射计算机断层照相(SPECT)、微计算机断层照相(micro CT)、宏计算机断层照相(macro CT)成像系统等)

图1是可以应用本发明的实施例的成像系统10的示图。本领域技术人员将意识到本发明适用于例如平面x射线、超声、磁共振(MR)、电子束X线断层照相术(EBT)、正电子发射计算机断层照相(PET)、单光子发射计算机断层照相(SPECT)、微计算机断层照相(micro CT)、宏计算机断层照相(macro CT)成像系统等。x射线系统10的下列讨论仅是一个这样的实现的示例并且从形态方面来说不意为限制性的。

参照图1，计算机断层摄影(CT)成像系统10包括机架12。机架12具有x射线管组件或x射线源组件14，其将x射线的锥形束朝图像检测单元组件或准直仪16投影在该机架12的对边上。现在参照图2，图像检测单元组件16由图像检测单元阵列18和数据采集系统(DAS)20形成。该图像检测单元阵列18感测通过医疗患者24的投影的x射线22，并且DAS 20将数据转换成数字信号用于后续处理。每个图像检测单元18产生模拟电信号，其代表碰撞x射线束的强度以及因此代表当它通过该患者24时的衰减束的强度。在扫描来采集x射线投影数据期间，机架12和安装在其上的部件绕着旋转中心26旋转。

机架12的旋转和x射线源组件14的操作由CT系统10的控制机构28管理。控制机构28包括：x射线控制器30，其向x射线源组件14提供电力和定时信号；和机架马达控制器32，其控制机架12的转速和位置。图像重建器34从DAS 20接收采样的和数字化的x射线数据并且进行高速重建。重建的图像提供给计算机36(其将该图像存储在大容量存储装置38中)。计算机36还具有存储在其上的对应于电子束定位和磁场控制的软件，如在下文详细描述的。

计算机36还经由控制台40从操作者接收命令和扫描参数，该控制台40具有例如键盘、鼠标、语音激励的控制器或任何其他合适的输入设备等某个形式的操作者接口。关联的显示器42允许该操作者观察来自计算机36的重建的图像和其他数据。该操作者提供的命令和参数被计算机36用来向DAS 20、x射线控制器30和机架马达控制器32提供控制信号和信息。另外，计算机36操作台架马达控制器44，其控制电动台架46来安置患者24和机架12。特别地，台架46使患者24全部或部分移动通过图1的机架开口48。

图3是说明包含故障图像检测单元(或待校正图像检测单元)和无故障图像检测单元的图像检测单元阵列300的示图。图像检测单元阵列300具有图像检测单元行(图中水平方向延伸)和图像检测单元列(图中竖直方向延伸)。图像检测单元阵列300包括故障图像检测单元330，位于行310和列320的交叉点。故障图像检测单元330可能不会报告准确的图像数据。故障图像检测单元330周围的图像检测单元的图像数据可以用来估计故障图像检测单元330的图像数据。

图4是说明按照本发明的一个实施例的用于校正位于图像检测单元阵列中位于行310和列320的故障图像检测单元330的图像数据的方法400的流程图。

在步骤410，输入对应于图像检测单元阵列300中各图像检测单元的图像数据。

在步骤420，基于故障图像检测单元330的相邻图像检测单元进行插值，得到故障图像检测单元330的估计图像数据(Xbase)。按照本发明的一个实施例，基于与故障图像检测单元330同一行310的故障图像检测单元330的相邻图像检测单元的图像数据进行插值。在本发明的一个实施例中，可以基于与故障图像检测单元330同一行310的故障图像检测单元330的相邻图像检测单元的图像数据，采用拉格朗日插值法，来得到故障图像检测单元330的估计图像数据。

在本发明的实施例中，故障图像检测单元的相邻图像检测单元包括故障图像检测单元的最紧邻的相邻图像检测单元，或者进一步地还包括故障图像检测单元的次紧邻的相邻图像检测单元，或者进一步地还包括故障图像检测单元的次次紧邻的相邻图像检测单元，等等。

可以采用多种方法来定位故障图像检测单元。例如，故障图像检测单元330可以采用诸如以下几种方法的至少一种方法来定位：(1)将图像检测单元阵列中的所有图像检测单元暴露到x射线下，然后将每个图像检测单元的量度与图像检测单元的相邻图像检测单元进行比较以寻找信号变化(variations)，以及(2)寻找可以标识故障图像检测单元的相邻图像检测单元之间的平均信号(数据)中的改变(changes)(对于每个图像检测单元行的所有发射观察)。

第一种方式称为空气校准扫描(air-calibration scan)。在一个实施例中，每个成像系统经历一组快速校准处理(称为“空气校准”)。空气校准处理将图像检测单元直接暴露到x射线光子射束，而没有任何对象在射束中。通过检查从相邻图像检测单元340获得的不同的读数，故障图像检测单元330是可以很快地标识的。随着时间的过去，故障图像检测单元330可能在某些时候适当运行，而在另一些时候不适当运行。故障图像检测单元330可以标识成作为相邻图像检测单元340、350的读数以下的一组百分比的图像检测单元。上述百分比还可以按照图像检测单元的位置而变化。

用于标识故障图像检测单元的第二种方式是自始至终使用患者扫描数据。第二种方式的优点是：并非所有医院或成像系统的操作者每天进行空气校准。因此，如果差错在空气校准扫描之间引入图像检测单元，图像检测单元可能不被适当地标识。通过使用平均读数(例如对于每个图像检测单元行的所有发射观察)，比较相邻图像检测单元之间的平均信号(数据)可以标识故障图像检测单元330。图像检测单元的发射观察平均可以保持每个信号的读数接近另一个，除非图像检测单元故障。

空气校准方式和患者扫描方式都可以通过位于各个位置的分开的硬件或软件来实施，例如在图像检测单元阵列18、DAS 32、图像重建器34、计算机36等。图像检测单元阵列随时间恶化可以被检查到。进行每次检查以标识故障图像检测单元330。图像检测单元读数中的改变可以被存储。可以检查图像检测单元读数中的改变以查找在图像检测单元阵列18中显著改变。可以标识改变，从而标识故障图像检测单元330和潜在地要替换的图像检测单元阵列18。可以保持超过一组阈值的所有已检测的变化的日志，从而标识故障图像检测单元330。该日志可以通过位于各个位置的分开的硬件或软件来保持，例如在图像检测器阵列18、DAS 32、图像重建器34、计算机36，或者例如保持在用于成像系统10中先前存在的硬件中。

按照本发明的实施例，判断故障图像检测单元330的相邻图像检测单元的图像数据的特性，基于该图像数据的特性对故障图像检测单元330的估计图像数据(Xbase)进行校正。例如，基于故障图像检测单元330的相邻图像检测单元的图像数据的阶跃特性对故障图像检测单元330的估计图像数据(Xbase)进行校正。本发明的实施例不限于基于故障图像检测单元330的相邻图像检测单元的图像数据的阶跃特性对故障图像检测单元330的估计图像数据(Xbase)进行校正，也可以基于故障图像检测单元330的相邻图像检测单元的图像数据的其他特性对故障图像检测单元330的估计图像数据(Xbase)进行校正。

在步骤430，判断故障图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点。图像数据的阶跃点特征是：图像检测单元阵列中某一图像检测单元行中的某一图像检测单元的左侧图像检测单元和右侧图像检测单元的有关图像数据相对平稳，但是左右侧的图像数据有明显的数值差异。在CT投影图像数据中，典型的图像数据的阶跃点出现在高密度物质的边缘。如果高密度物质的形状缓慢变化则投影图像数据也会呈现“缓慢连续”的变化，参见图5a所示的情形。如果高密度物质的形状有棱角，则很可能会在投影图像数据中出现图像数据的阶跃点，参见图5b所示的情形。当故障图像检测单元的相邻图像检测单元的图像数据存在图像数据的阶跃点时，意味着该故障图像检测单元的实际图像数据可能处于图像数据的阶跃点。本申请的发明人洞察到，对于图像数据的阶跃点，传统的插值方法很难恢复其真实的图像数据。低阶插值对遇到准确性不够的问题，而高阶插值又会引入插值过冲，都会影响插值的准确性。这种不准确在图像上体现为伪影，以条纹状伪影最为常见。伪影的引入会对临床产生误诊，这是医学图像上极力避免的。本发明的实施例的方法考虑图像数据的阶跃点因素，对故障图像检测单元330的估计图像数据采取不同的方式进行校正，以克服现有技术的方法的不足。

图6是说明按照本发明的一个实施例的故障图像检测单元的相邻图像检测单元的图像数据是否存在图像数据的阶跃点判断方法的示图。图6中，X1、X2、X3代表位于图像检测单元C同一行的左侧的相邻图像检测单元的图像数据，X4、X5、X6代表位于图像检测单元C同一行的右侧的相邻图像检测单元的图像数据，其中X3是图像检测单元C左侧与图像检测单元C最相邻的图像检测单元的图像数据，X2是图像检测单元C左侧与图像检测单元C次相邻的图像检测单元的图像数据，X1是在图像检测单元C左侧与图像检测单元C次次相邻的图像检测单元的图像数据，X4是图像检测单元C右侧与图像检测单元C最相邻的图像检测单元的图像数据，X5是故障图像检测单元C右侧与图像检测单元C次相邻的图像检测单元的图像数据，X6是图像检测单元C右侧与图像检测单元C次次相邻的图像检测单元的图像数据。计算图像检测单元C所在行的有关相邻图像检测单元的图像数据之间δ值：δ1＝|X2－X3|，δ2＝|X3－X4|，δ3＝|X4－X5|，如果δ2>3×max(δ1，δ3)并且δ2>0.8×(max(Xi)－min(Xi))，其中则认为图像检测单元C处于图像数据的阶跃点。运用以上判断图像数据的阶跃点的方法来判断故障图像检测单元的相邻图像检测单元的是否处于图像数据的阶跃点，本发明的判断图像数据的阶跃点方法不限于上述方法，本领域技术人员可以对以上方法进行修改、或者采用其它方法判断图像检测单元的图像数据是否处于图像数据的阶跃点。

当在步骤420中确定故障图像检测单元330的相邻图像检测单元的图像数据处于图像数据的阶跃点时，进入步骤450来执行第一校正方法，否则，进入步骤470来执行第二校正方法。

在步骤490，提供来自步骤450或来自步骤470的校正估计图像数据。

在本发明的一个实施例中，第一校正方法包括，判断图像检测单元阵列中故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在阵列中的是否存某种趋势，例如，表现出一定的方向性。

图7是说明图像数据的阶跃点的方向性的示图。可以通过检查图像检测单元阵列中故障图像检测单元的相邻图像检测单元的图像数据的阶跃点的分布来确定其是否存在方向性。图7a表示故障图形检测单元(图中中心点，即灰度最深的块)的左下和右上相邻图形检测单元是阶跃点(灰度最浅的块)，图7a中的虚线反映阶跃点的方向性。图7b表示故障图形检测单元(图中中心点，即灰度最深的块)的右上和下方相邻图形检测单元是阶跃点(灰度最浅的块)，图7b中的虚线反映阶跃点的方向性。图7c表示故障图形检测单元(图中中心点，即灰度最深的块)的上和下相邻图形检测单元是阶跃点(灰度最浅的块)，图7c中的连接阶跃点的虚线反映阶跃点的方向性。本发明的实施例不限于图7中所例示的情形，其它反映阶跃点方向性的情形也适用本发明。

如果确定故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在图像检测单元阵列中存在方向性，则以故障图像检测单元为中心，按照代表所述方向性的方向选择故障图像检测单元的相邻图像检测单元(例如图7中虚线连接的相邻图像检测单元)的图像数据，基于所选择的相邻图像检测单元的图像数据对步骤420中得到的故障图像检测单元的估计图像数据进行方向校正。按照本发明的一个实施例，计算所选相邻图像检测单元中的每一个图像检测单元(i)的图像数据与故障图像检测单元的估计图像数据的偏差(Δ_i)，计算所选相邻图像检测单元中的所有图像检测单元的图像数据与故障图像检测单元的估计图像数据的偏差(Δ_i)的平均偏差将上述平均偏差与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的方向性校正后的校正估计图像数据(X)。

可选地，如果确定故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在图像检测单元阵列中存在方向性，则针对故障图像检测单元进行单边外插校正。按照本发明的一个实施例的校正步骤如下。

(1)邻行双方向外插估计

参见图8a，邻行双方向外插估计包括故障图像检测单元的上一行左侧外插估计上一行的外插值x_{up_left}，上一行右侧外插估计上一行的外插值x_{up_right}，下一行左侧外插估计上一行的外插值x_{down_left}，下一行右侧外插估计上一行的外插值x_{down_right}。

x_{up_left}＝extra_lagrange(x_{up_left_1}，x_{up_left_2}，x_{up_left_3})；

x_{up_right}＝extra_lagrange(x_{up_right_1}，x_{up_right_2}，x_{up_right_3})；

x_{down_left}＝extra_lagrange(x_{down_left_1}，x_{down_left_2}，x_{down_left_3})；

x_{down_right}＝extra_lagrange(x_{down_right_1}，x_{down_right_2}，x_{down_right_3})；

式中extra_lagrange()表示针对括号内的参数进行拉格朗日计算。

(2)邻行双方向的权重分析

参见图8b，通过外插值和故障图像检测单元竖直方向点的图像数据，可以计算出左右外插点得权重。

即，通过将外插值x_{up_left}、x_{up_right}、x_{down_left}、x_{down_right}、和故障图像检测单元竖直方向的相邻图像检测单元的图像数据x_up、x_down代入下式来求解方程，

x_up＝w_{up_left}*x_{up_left}+w_{up_right}*x_{up_right}

x_down＝w_{down_left}*x_{down_left}+w_{down_right}*x_{down_right}

w_{up_left}+w_{up_right}＝1

w_{down_left}+w_{down_right}＝1

从而得到权重w_{up_left}、w_{down_left}、w_{up_right}、和w_{down_right}。

(3)本行双方向权重计算

根据下式计算本行双方向权重w_left和w_right。

w_left＝a*w_{up_left}+(1-a)*w_{down_left}

w_right＝bw_{up_right}+(1-b)*w_{down_right}

其中，常数a和b是0至1之间的数值。

可选地，根据下式计算本行双方向权重w_left和w_right。

w_left＝0.5*w_{up_left}+0.5*w_{down_left}

w_right＝0.5*w_{up_right}+0.5*w_{down_right}

(4)单边外插结果

单边外插的结果综合考虑两个方向的外插值，通过下式计算X。

x＝w_left*x_left+w_right*x_right

然后，用所求得的X作为故障图像检测单元的单边外插校正后的校正估计图像数据(X)。

按照本发明的一个实施例，如果确定故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在图像检测单元阵列中存在方向性，则针对故障图像检测单元同时进行方向性校正和单边外插校正。然后，判断方向性校正中算得的平均偏差是否大于零，如果平均偏差大于零，则取方向性校正后的校正估计图像数据和单边外插校正后的校正估计图像数据中最大者作为故障图像检测单元的校正估计图像数据(X)，否则，取方向性校正后的校正估计图像数据和单边外插校正后的校正估计图像数据中最小者作为故障图像检测单元的校正估计图像数据(X)。

在本发明的一个实施例中，第二校正方法包括，通过故障图像检测单元上下行相邻图像检测单元的图像数据进行delta校正。可选地，计算故障图像检测单元在图像检测单元阵列中竖直方向上的上下相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta(与故障图像检测单元的估计图像数据Xbase的差值)，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)。可选地，计算故障图像检测单元上、下、左、右、左上、右上、左下、右下的相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)。

按照本发明的一个实施例，第二校正方法还包括判断图像检测单元在图像检测单元阵列中竖直方向的尺寸与行方向的尺寸之比R是否大于或等于预定阈值V_threshold。如果R大于或等于预定阈值V_threshold，则计算故障图像检测单元上、下、左、右、左上、右上、左下、右下的相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)。如果R小于预定阈值V_threshold，计算故障图像检测单元在图像检测单元阵列中竖直方向上的上下相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)。预定阈值V_threshold可以在0.5至2数值范围内取值。

图9是说明按照本发明的一个优选的实施例的故障图像检测单元的图像数据校正方法的示图。在步骤910，接收图像检测单元阵列的图像数据，其中图像检测单元阵列中包含有故障图像检测单元或待校正的图像检测单元。在步骤920，基于故障图像检测单元在图像检测单元阵列所在行的相邻图像检测单元的图像数据对故障图像检测单元进行插值计算，得到故障图像检测单元的估计图像数据(Xbase)。在步骤930，判断故障图像检测单元的相邻图像检测单元的图像数据是否存在图像数据的阶跃点，如果存在图像数据的阶跃点，则进入步骤940，否则，进入步骤970。在步骤940，判断故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在图像检测单元阵列中是否具有方向性，如果具有方向性，进入步骤970，否则，进入步骤952和步骤953。在步骤952，计算上述方向性所表示的方向上处于故障图像检测单元的相邻图像检测单元的图像数据与故障图像检测单元的估计图像数据的差值(Δ_i)。在步骤954，对在步骤952中计算的所有差值求(Δ_i)平均，得到平均差值并将平均差值与故障图像检测单元的估计图像数据(Xbase)相加得到故障图像的方向性校正后的校正估计图像数据(Xin)。在步骤953，针对故障图像检测单元的相邻图像检测单元的图像数据来求解外插权重Wl和Wr。在步骤957，基于步骤953中得到的外插权重Wl和Wr，计算故障图像检测单元的单边外插值，得到故障图像检测单元的单边外插后的校正估计图像数据(Xextra)。在步骤960，判断上述平均差值是否大于0，如果大于0，进入步骤962，否则，进入步骤963。在步骤962，取方向性校正后的校正估计图像数据(Xin)和单边外插校正后的校正估计图像数据(Xextra)中最大者作为故障图像检测单元的校正估计图像数据(X)。在步骤963，取方向性校正后的校正估计图像数据(Xin)和单边外插校正后的校正估计图像数据(Xextra)中最小者作为故障图像检测单元的校正估计图像数据(X)。在步骤990，提供校正估计图像数据。

图10是说明按照本发明的一个实施例的用于校正图像检测阵列中的故障图像检测单元330的设备的方框图。该设备包括：图像数据输入装置1010，用于输入对应于图像检测单元阵列300中各图像检测单元的图像数据；图像数据估计装置1020，用于针对故障图像检测单元基于该故障图像检测单元330的相邻图像检测单元进行插值，得到故障图像检测单元330的估计图像数据(Xbase)；图像数据阶跃点判断装置1030；用于判断故障图像检测单元的相邻图像检测单元的图像数据是否存在图像数据的阶跃点；第一校正装置1050，用于在图像数据阶跃点判断装置在判断故障图像检测单元的相邻图像检测单元的图像数据存在图像数据的阶跃点时，采用第一校正模式对故障图像检测单元的估计图像数据进行校正；第二校正装置1070，用于在图像数据阶跃点判断装置在判断故障图像检测单元的相邻图像检测单元的图像数据不存在图像数据的阶跃点时，采用第二校正模式对故障图像检测单元的估计图像数据进行校正；图像数据输出装置1090，用于输出故障图像检测单元的校正估计图像数据。

图像数据阶跃点判断装置1030可以进一步地包括：第一判断装置，用于判断故障图像检测单元的相邻图像检测单元的图像数据是否存在图像数据的阶跃点；和第二判断装置，用于在判断故障图像检测单元的相邻图像检测单元的图像数据存在图像数据的阶跃点时判断图像数据的阶跃点在图像检测单元阵列中的分布是否存在方向性。

第一校正装置1050可以进一步包括：方向性校正装置，用于在确定故障图像检测单元的相邻图像检测单元的图像数据的阶跃点在图像检测单元阵列中存在方向性时，以故障图像检测单元为中心，按照代表所述方向性的方向选择故障图像检测单元的相邻图像检测单元的图像数据，基于所选择的相邻图像检测单元的图像数据对故障图像检测单元的估计图像数据进行方向性校正，得到故障图像检测单元的方向性校正后的校正估计图像数据；单边外插校正装置，用于针对故障图像检测单元进行单边外插，以得到故障图像检测单元的单边外插校正后的校正估计图像数据；方向性－单边外插校正选择装置，用于选择是采用故障图像检测单元的方向性校正的校正估计图像数据还是选择故障图像检测单元的单边外插校正后的校正估计图像数据。

第二校正装置可以进一步包括判断装置，用于判断图像检测单元在图像检测单元阵列中竖直方向的尺寸与行方向的尺寸之比R是否大于或等于预定阈值V_threshold。当判断装置判断R大于或等于预定阈值V_threshold，第二校正装置计算故障图像检测单元上、下、左、右、左上、右上、左下、右下的相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)，否则，计算故障图像检测单元在图像检测单元阵列中竖直方向上的上下相邻图像检测单元中每一个相邻图像检测单元的图像数据的delta，对所计算的所有这些delta计算平均值，将该平均值与故障图像检测单元的估计图像数据(Xbase)相加，得到故障图像检测单元的校正估计图像数据(X)。

本发明的实施例的用于校正故障图像检测单元330的设备1000的各个组成部分可以整体地布置在例如图2所示系统的一个组成部分中(例如，位于DAS 20、或者位于图像重建器34等)，也可以分布到图2所示系统的不同组成部分(例如，分布于DAS 20和图像重建器34等)。本发明的实施例的用于校正故障图像检测单元330的设备1000的各个组成部分不仅可以是硬件，还可以是固件、或者硬件与固件的结合，或者硬件、固件与软件的结合。

图11示出了按照本发明的实施例进行图像数据校正前后的图像比较。右边的图像是按照本发明的实施例进行图像数据校正之前的图像(采用现有校正技术进行校正的图像结果)，图中存在“伪影”，如图中箭头所指示的条纹；左边的图像是按照本发明的实施例进行图像数据校正之后的图像，其中没有象右边图像中那样的“伪影”。

以上通过具体的实施例对本发明进行了说明，但本发明并不限于这些具体的实施例。本领域技术人员将领会，还可以对本发明进行各种修改、等同替换、变化等。例如将上述实施例中的一个步骤或模块分为两个或更多个步骤或模块来实现，或者相反，将上述实施例中的两个或更多个步骤或模块或装置的功能放在一个步骤或模块中来实现。但是，这些变换只要未背离本发明的精神，都应在本发明的保护范围之内。另外，本申请说明书和权利要求书所使用的一些术语并不是限制，而仅仅是为了便于描述。

Claims (35)

1.一种用于处理图像数据的方法，所述方法包括：

图像数据获取步骤，获取图像检测单元阵列中待校正图像检测单元的相邻图像检测单元的图像数据；

图像数据估计步骤，基于所述相邻图像检测单元的所述图像数据中的与所述待校正图像检测单元相同行的相邻图像检测单元的图像数据对所述待校正图像检测单元进行插值，得到所述待校正图像检测单元的估计图像数据；以及

图像数据修正步骤，基于所述相邻图像检测单元的图像数据的特性对所述估计图像数据进行修正。

2.如权利要求1所述的方法，其中，所述图像数据修正步骤包括：

阶跃点判断步骤，判断所述待校正图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点；以及

图像数据校正步骤，基于上述判断结果选择不同的校正模式对所述估计图像数据进行校正。

3.如权利要求2所述的方法，其中，所述图像数据校正步骤还包括：

在所述相邻图像检测单元的图像数据中存在图像数据的阶跃点时，按照第一校正模式对所述估计图像数据进行校正，否则，按照第二校正模式对所述估计图像数据进行校正。

4.如权利要求3所述的方法，其中，所述第一校正模式包括：

沿这样的方向选择所述待校正图像检测单元的相邻图像检测单元，所述方向与所述图像检测单元阵列的行不同，反映了所述相邻图像检测单元的图像数据中阶跃点的分布方向；

计算所选相邻图像检测单元的图像数据的方向偏差Δ，以及

基于所述方向偏差Δ来校正所述估计图像数据，得到方向性校正后的图像数据。

5.如权利要求4所述的方法，其中，计算所述待校正图像检测单元的相邻图像检测单元的图像数据的方向偏差Δ包括：

计算所选相邻图像检测单元的图像数据与所述待校正图像检测单元的估计图像数据的差值Δ_i，以及

计算所述差值Δ_i的平均值作为所述方向偏差Δ。

6.如权利要求3-5中任意一项权利要求所述的方法，其中，所述第一校正模式包括：

对所述待校正图像检测单元进行单边外插，得到所述待校正图像检测单元的单边外插后的图像数据。

7.如权利要求6所述的方法，其中，对所述待校正图像检测单元进行单边外插包括：

计算所述待校正图像检测单元在所述图像检测单元阵列中的左侧权重w_left和右侧权重w_right，以及

根据所述左侧权重w_left和右侧权重w_right来计算单边外插后的图像数据X：

X＝W_left*X_left+W_right*X_right，

其中，X_left和X_right分别是所述待校正图像检测单元左侧相邻图像检测单元和右侧相邻图像检测单元的图像数据。

8.如权利要求7所述的方法，其中，按照以下方式计算所述左侧权重w_left和所述右侧权重w_right：

针对所述待校正图像检测单元，上一行左侧外插估计上一行的外插值x_{up_left}，上一行右侧外插估计上一行的外插值x_{up_right}，下一行左侧外插估计上一行的外插值x_{down_left}，下一行右侧外插估计上一行的外插值x_{down_right}，其中，所述外插值是采用相应行的相应侧的相邻图像检测单元的图像数据进行拉格朗日插值，以及

通过所述外插值x_{up_left、}x_{up_right、}x_{down_left}、x_{down_right}和所述待校正图像检测单元的上相邻图像检测单元和下相邻图像检测单元的图像数据x_up、x_down来求解以下方程而获得左外插权重w′_{up_left}、W_{down_left}和右外插权重w_{up_right}、w_{down_righht}：

X_up＝W_{up_left}*X_{up_left}+W_{up_right}*X_{up_right}

X_down＝W_{down_left}*X_{down_left}+W_{down_right}*X_{down_right}

W_{up_left}+W_{up_right}＝1，以及

W_{down_left}+W_{down_right}＝1

通过所获得的左外插权重w_{up_left}、W_{down_left}、和右外插权重w_{up_right}、w_{down_right}来计算所述左侧权重w_left和所述右侧权重w_right：

W_left＝a*w_{up_left}+(1-a)*W_{down_left}

W_right＝bw_{up_right}+(1-b)*w_{down_right}

其中，a和b是取值在0至1之间的常数。

9.如权利要求4或5所述的方法，其中，将所述方向性校正后的图像数据作为所述待校正图像检测单元的校正后的图像数据。

10.如权利要求6所述的方法，其中，

采用所述单边外插后的图像数据作为所述待校正图像检测单元的校正后的图像数据；或者

计算所述待校正图像检测单元的所选相邻图像检测单元的图像数据的方向偏差Δ，在所述方向偏差Δ大于0时，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最大者作为所述待校正图像检测单元的校正后的图像数据，否则，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最小者作为所述待校正图像检测单元的校正后的图像数据。

11.如权利要求3所述的方法，其中，所述第二校正模式包括：

计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值，基于所述偏差的平均值对所述估计图像数据进行校正，得到所述待校正图像检测单元的校正后的图像数据。

12.如权利要求11所述的方法，其中，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：

计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

13.如权利要求11所述的方法，其中，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：

计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

14.如权利要求11所述的方法，其中，所述计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值包括：

判断图像检测单元在所述图像检测单元阵列中列方向的尺寸与行方向的尺寸之比R是否大于或等于预定阈值V_threshold，如果R大于或等于所述预定阈值V_threshold，计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值，否则，计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值，作为所述偏差的平均值。

15.如权利要求14所述的方法，其中所述所述预定阈值V_threshold是介于0.5至2之间的数值。

16.如权利要求1所述的方法，其中，所述图像数据估计步骤还包括：

基于所述待校正图像检测单元所在行的相邻图像检测单元的图像数据对所述待校正图像检测单元进行拉格朗日插值，得到所述估计图像数据。

17.一种用于处理图像数据的设备，所述设备包括：

图像数据获取装置，获取图像检测单元阵列中待校正图像检测单元的相邻图像检测单元的图像数据；

图像数据估计装置，基于所述相邻图像检测单元的所述图像数据中的与所述待校正图像检测单元相同行的相邻图像检测单元的图像数据对所述待校正图像检测单元进行插值，得到所述待校正图像检测单元的估计图像数据；以及

图像数据修正装置，基于所述相邻图像检测单元的图像数据的特性对所述估计图像数据进行修正。

18.如权利要求17所述的设备，其中，所述图像数据修正装置包括：

阶跃点判断装置，判断所述待校正图像检测单元的相邻图像检测单元的图像数据中是否存在图像数据的阶跃点；以及

图像数据校正装置，基于上述判断结果选择不同的校正模式对所述估计图像数据进行校正。

19.如权利要求18所述的设备，其中，所述图像数据校正装置还包括：

第一模式校正装置，用于在所述相邻图像检测单元的图像数据中存在图像数据的阶跃点时，按照第一校正模式对所述估计图像数据进行校正，以及

第二模式校正装置，用于在所述相邻图像检测单元的图像数据中不存在图像数据的阶跃点时，按照第二校正模式对所述估计图像数据进行校正。

20.如权利要求19所述的设备，其中，所述第一模式校正装置还包括：

方向选择装置，用于沿这样的方向选择所述待校正图像检测单元的相邻图像检测单元，所述方向与所述图像检测单元阵列的行不同，反映了所述相邻图像检测单元的图像数据中阶跃点的分布方向；

方向偏差计算装置，用于计算所选相邻图像检测单元的图像数据的方向偏差Δ，以及

方向性校正装置，用于基于所述方向偏差Δ来校正所述估计图像数据，得到方向性校正后的图像数据。

21.如权利要求20所述的设备，其中，所述方向偏差计算装置还包括：

用于计算所选相邻图像检测单元的图像数据与所述待校正图像检测单元的估计图像数据的差值Δ_i的装置，以及

用于计算所述差值Δ_i的平均值并所述该平均值作为所述方向偏差Δ的装置。

22.如权利要求19-21中任意一项权利要所述的设备，其中，所述第一模式校正装置还包括：

用于对所述待校正图像检测单元进行单边外插以得到所述待校正图像检测单元的单边外插后的图像数据的装置。

23.如权利要求22所述的设备，其中，所述用于对所述待校正图像检测单元进行单边外插以得到所述待校正图像检测单元的单边外插后的图像数据的装置还包括：

用于计算所述待校正图像检测单元在所述图像检测单元阵列中的左侧权重w_left和右侧权重w_right的装置，以及

用于根据所述左侧权重w_left和右侧权重w_right来计算单边外插后的图像数据X的装置，其中，

X＝W_left*X_left+W_right*X_right，

其中，X_left和X_right分别是所述待校正图像检测单元左侧相邻图像检测单元和右侧相邻图像检测单元的图像数据。

24.如权利要求23所述的设备，其中，所述用于计算所述待校正图像检测单元在所述图像检测单元阵列中的左侧权重w_left和右侧权重w_right的装置：

针对所述待校正图像检测单元，上一行左侧外插估计上一行的外插值x_{up_left}，上一行右侧外插估计上一行的外插值x_{up_right}，下一行左侧外插估计上一行的外插值x_{down_left}，下一行右侧外插估计上一行的外插值x_{down_right}，其中，所述外插值是采用相应行的相应侧的相邻图像检测单元的图像数据进行拉格朗日插值，以及

通过所述外插值x_{up_left}、x_{up_right}、x_{down_left}、x_{down_right}和所述待校正图像检测单元的上相邻图像检测单元和下相邻图像检测单元的图像数据x_up、x_down来求解以下方程而获得左外插权重w_{up_left}、W_{down_left}和右外插权重w_{up_right}、w_{down_right}：

X_up＝W_{up_left}*X_{up_left}+W_{up_right}*X_{up_right}

X_down＝W_{down_left}*X_{down_left}+W_{down_right}*X_{down_right}

W_{up_left}+W_{up_right}＝1，以及

W_{down_left}+W_{down_right}＝1

通过所获得的左外插权重w_{up_left}、W_{down_left}、和右外插权重w_{up_right}、w_{down_right}来计算所述左侧权重w_left和所述右侧权重w_right：

W_left＝a*w_{up_left}+(1-a)*W_{down_left}

W_right＝bw_{up_right}+(1-b)*w_{down_right}

其中，a和b是取值在0至1之间的常数。

25.如权利要求20或21所述的设备，其中，所述图像数据校正装置将所述方向性校正后的图像数据作为所述待校正图像检测单元的校正后的图像数据。

26.如权利要求22所述的设备，其中，所述图像数据校正装置：

采用所述单边外插后的图像数据作为所述待校正图像检测单元的校正后的图像数据；或者

计算所述待校正图像检测单元的所选相邻图像检测单元的图像数据的方向偏差Δ，在所述方向偏差Δ大于0时，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最大者作为所述待校正图像检测单元的校正后的图像数据，否则，选择所述方向性校正后的图像数据和所述单边外插后的图像数据中的最小者作为所述待校正图像检测单元的校正后的图像数据。

27.如权利要求19所述的设备，其中，所述第二模式校正装置还包括：

用于计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值的装置，以及

用于基于所述偏差的平均值来校正所述估计图像数据以得到所述待校正图像检测单元的校正后的图像数据的装置。

28.如权利要求27所述的设备，其中，所述用于计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值的装置包括：

用于计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值并将所述算术差值的平均值作为所述偏差的平均值的装置。

29.如权利要求27所述的设备，其中，所述用于计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值的装置包括：

用于计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值并将所述算术差值的平均值作为所述偏差的平均值的装置。

30.如权利要求27所述的设备，其中，所述用于计算所述待校正图像检测单元的相邻图像检测单元的图像数据与所述估计图像数据的偏差的平均值的装置还包括：

阈值判断装置，用于判断图像检测单元在所述图像检测单元阵列中列方向的尺寸与行方向的尺寸之比R是否大于或等于预定阈值V_threshold，

用于在阈值判断装置判断R大于或等于所述预定阈值V_threshold的情况下计算所述待校正图像检测单元的上、下、左、右、上左、上右、下左、下右相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值并将所述算术差值的平均值作为所述偏差的平均值的装置，以及

用于在阈值判断装置判断R小于所述预定阈值V_threshold的情况下计算所述待校正图像检测单元的上、下相邻图像检测单元的图像数据与所述估计图像数据的算术差值的平均值并将所述算术差值的平均值作为所述偏差的平均值的装置。

31.如权利要求30所述的设备，其中，所述所述预定阈值V_threshold是介于0.5至2之间的数值。

32.如权利要求17所述的设备，其中，所述图像数据估计装置还包括：

用于基于所述待校正图像检测单元所在行的相邻图像检测单元的图像数据对所述待校正图像检测单元进行拉格朗日插值以得到所述估计图像数据的装置。

33.一种用于处理图像数据的产品，包括用于执行权利要求1-16中任何一项中的步骤的单元。

34.一种成像设备，其特征在于包括如权利要求17-32中任何一项所述的用于处理图像数据的设备。

35.如权利要求34所述的成像设备，其中，所述成像设备选自CT成像设备、平面x射线成像设备、超声成像设备、磁共振(MR)成像设备、电子束x线断层照相术(EBT)成像设备、正电子发射计算机断层照相(PET)成像设备、单光子发射计算机断层照相(SPECT)成像设备、微计算机断层照相(micro CT)成像设备、宏计算机断层照相(macro CT)成像设备之一。