CN106296595B - 一种平板探测器以及降低平板探测器图像残影的方法 - Google Patents
一种平板探测器以及降低平板探测器图像残影的方法 Download PDFInfo
- Publication number
- CN106296595B CN106296595B CN201510309168.5A CN201510309168A CN106296595B CN 106296595 B CN106296595 B CN 106296595B CN 201510309168 A CN201510309168 A CN 201510309168A CN 106296595 B CN106296595 B CN 106296595B
- Authority
- CN
- China
- Prior art keywords
- image
- target
- ghost
- point
- correction
- 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.)
- Active
Links
Landscapes
- Measurement Of Radiation (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种平板探测器以及降低平板探测器图像残影的方法，所述方法首先实现目标图像和对比图像的校正预处理工作，然后根据对比图像和目标图像之间的时间间隔、曝光窗口等变量，实现图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算，最后通过减法公式，将图像残影从目标图像中减去，以达到最终降低或者消除图像残影的目的。本发明通过在目标图像前插入多帧自动采集的对比图像的方式，利用校正后的对比图像具有相同性质的图像残影，对目标图像上含残影信息的各像素点的进行公式计算，实现在不增加现有硬件成本的基础上，降低或消除图像残影的目的，且该方法具备高度的集成性、灵活性和可观的市场应用价值。
Description
技术领域
本发明属于医疗设备及图像处理领域，特别是涉及一种平板探测器以及降低平板探测器图像残影的方法。
背景技术
X射线平板探测器主流为非晶硅和非晶硒，由于非晶硅和非晶硒平板探测器本身的物理特性：一定拍摄条件下，它们的曝光图像上都会存在着残影的问题，即第n次(n>1)成像时所获得的影像残留第1～n-1次被拍摄物体的影响，且这种影响程度根据第1～n-1次的顺序，呈现逐渐递增的效果。一旦出现图像残影，则需要很长时间进行自然衰减，才能完全消除图像残影的影响。
当前，大多研发平板探测器的企业或研究单位，主要还是通过硬件上的改进来降低图像残影，如改善玻璃的工艺制成。但是，这些技术手段繁琐、复杂且昂贵。此外，随着计算机处理技术不断应用于医疗器械领域，也有部分企业或研究机构从图像处理技术角度，设计理论算法，进行计算机图像后处理，如专利CN103377464A对含残影的图像进行时频变换，并进行图像去噪、图像特征提取等一系列相关算法的处理，这在一定程度上实现了图像残影的降低或消除。但是，这类理论算法，不具备很强的适用性，对于不同X射线剂量下的不同物体的残影，依然无法做到实时、有效合理地降低或者消除。而且，有时还会带入副作用，图像在降低或者消除残影的同时，往往会出现不同程度的图像失真，影响图像的完整性。
发明内容
鉴于以上所述现有技术的缺点，本发明的目的在于提供一种平板探测器以及降低平板探测器图像残影的方法，用于解决现有技术中平板探测器图像残影难以有效消除的问题。
为实现上述目的及其他相关目的，本发明提供一种降低平板探测器图像残影的方法，包括步骤：
1)采集第一目标图像，并记录采集时间点T1、曝光窗口的大小D1，采集原始暗场图像n帧，所述第一目标图像、原始暗场图像具有固定偏移量；
2)对第一目标图像进行固定偏移量校正，算法包括：
其中，n根据第一目标图像的固定偏移量校正方式以及第一帧原始暗场图像的残影量确定；Pi(x,y)是某幅图像上各个像素点，x为图像的宽度坐标，y为图像的高度坐标；
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x,y)是图像I经过固定偏移量校正后的亮场图；Avg_I是图像I的所有像素点累加后的平均值；
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，Dark_1为第一帧原始暗场图像，Light_1为第一目标图像，C为预设的固定常数；
3)采集第二目标图像，并记录采集时间点T2、曝光窗口的大小D2，采集原始暗场图像n帧，并根据步骤2)的固定偏移量校正的算法，获得Ghost_template2以及Target_image2；
4)根据图像残影的自然衰减规律，利用T1和T2之间的时间间隔以及D1和D2的变化，设置加权系数，作用到Ghost_template1，用以定量对比图像Ghost_template1到目标图像Target_image2之间的图像残影衰减量；
5)对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像。
作为本发明的降低平板探测器图像残影的方法的一种优选方案，所述固定偏移量为由于平板探测器中的电子电路、ROIC、Gate IC及像素本身存在的干扰信息造成原始暗场图像、及第一目标图像存在的本底噪声。
作为本发明的降低平板探测器图像残影的方法的一种优选方案，步骤2)中，还包括对原始暗场图像、及第一目标图像进行Gain校正及Defect校正的步骤，步骤3)中，还包括对原始暗场图像、及第二目标图像进行Gain校正及Defect校正的步骤。
作为本发明的降低平板探测器图像残影的方法的一种优选方案，步骤4)中，获得加权系数K的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1,T2)是一个自变量是T1、T2的函数，I0(x,y)是某一(T1,T2)点下经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
进一步地，步骤4)中，由于前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，需根据曝光窗口设置漏电流补偿系数Kd，拟合各个曝光窗口下的图像残影信息情况，确定其符合的公式：
其中，W(D1,D2)是自变量为D1、D2的函数，G(T1,T2)是某一(T1,T2)点下的理论计算值，I0(x,y)是某一(D1,D2)点下，对应某一(T1,T2)点的经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(D1,D2)点和(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
进一步地，步骤5)中，对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x,y)是目标图像Target_image2上的各个像素点灰度值，PG1(x,y)是对比图像Ghost_template1上的各个像素点灰度值，K1是微调系数。
本发明还提供一种平板探测器，包括：
图像采集模块，用于采集图像；
图像处理模块，用于实现目标图像和对比图像的图像校正工作，根据对比图像和目标图像之间的时间间隔以及曝光窗口大小，实现影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算，并通过减法公式，将图像残影从目标图像中减去，以终降低或者消除图像残影，获得最终的图像；
系统控制模块，用于实现图像采集模块、图像处理模块及图像显示模块的控制；
图像显示模块，用于显示最终的图像。
作为本发明的平板探测器的一种优选方案，所述图像处理模块中，目标图像和对比图像的图像校正工作包括：offset校正、Gain校正及Defect校正。
进一步地，所述offset校正包括：
依据采集的第一目标图像的时间点T1、曝光窗口的大小D1，以及采集的n帧原始暗场图像，采用算法：
其中，n根据第一目标图像的固定偏移量校正方式以及第一帧原始暗场图像的残影量确定；Pi(x,y)是某幅图像上各个像素点，x为图像的宽度坐标，y为图像的高度坐标；
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x,y)是图像I经过固定偏移量校正后的亮场图；Avg_I是图像I的所有像素点累加后的平均值；
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，Dark_1为第一帧原始暗场图像，Light_1为第一目标图像，C为预设的固定常数；
依据采集第二目标图像的时间点T2、曝光窗口的大小D2，以及采集的n帧原始暗场图像，依据上述的offset校正的算法，获得第二目标图像的Ghost_template2以及Target_image2。
进一步地，所述图像处理模块中，影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算包括：
获得加权系数K的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1,T2)是一个自变量是T1、T2的函数，I0(x,y)是某一(T1,T2)点下经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值；
由于前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，需根据曝光窗口设置漏电流补偿系数Kd，拟合各个曝光窗口下的图像残影信息情况，确定其符合的公式：
其中，W(D1,D2)是自变量为D1、D2的函数，G(T1,T2)是某一(T1,T2)点下的理论计算值，I0(x,y)是某一(D1,D2)点下，对应某一(T1,T2)点的经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(D1,D2)点和(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
进一步地，所述图像处理模块中，通过减法公式，将图像残影从目标图像中减去，以终降低或者消除图像残影，获得最终的图像的步骤包括：
对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x,y)是目标图像Target_image2上的各个像素点灰度值，PG1(x,y)是对比图像Ghost_template1上的各个像素点灰度值，K1是微调系数。
如上所述，本发明的平板探测器以及降低平板探测器图像残影的方法，具有以下有益效果：本发明解决的技术问题是在保证图像质量的同时，降低或者消除当前帧(第n次，且n>1)图像上残留的前几帧图像(第1～n-1次)的信息。通过在目标图像前插入多帧自动采集的对比图像的方式，利用校正后的对比图像具有相同性质的图像残影，对目标图像上含残影信息的各像素点的进行公式计算，实现在不增加现有硬件成本的基础上，降低或消除图像残影的目的，且该方法具备高度的集成性、灵活性和可观的市场应用价值。
附图说明
图1显示为本发明的降低平板探测器图像残影的方法的步骤流程示意图。
图2显示为本发明的平板探测器的结构框图。
元件标号说明
1 图像采集模块
11 闪烁体以及TFT面板
12 信号读出、放大、模数转换子模块
2 图像处理模块
21 预设图像基本校正子模块
22 图像Ghost校正子模块
3 系统控制模块
4 图像显示模块
S11～S15 步骤
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。
请参阅图1～图2。需要说明的是，本实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
如图1所示，本实施例提供一种降低平板探测器图像残影的方法，从成像角度来说，所谓的图像残影是指余辉增殖效应，是采集两张或两张以上亮场图时会出现的一种伪影现象。所述降低平板探测器图像残影的方法包括以下步骤：
步骤1)S11当平板探测器开机后预热一定时间后，根据预先设置好的曝光窗口(Exposure window)，进行第一次曝光，采集到第一目标图像，命名为Light_1，并记录采集时间点T1、曝光窗口(Exposure window)的大小D1；
当采集Light_1完成后，平板探测器会自动采集原始暗场图像n帧，由于平板探测器中的电子电路、ROIC、Gate IC及像素本身等存在干扰信息，会造成原始暗场图像、Light_1存在本底噪声，这是一种固定的偏移量(offset)。因此，需要对原始暗场图像、Light_1进行offset、gain、defect校正；
步骤2)S12，对原始暗场图像、Light_1的offset校正的具体算法如下：
其中，n是根据Light_1的offset校正方式以及Dark_1的残影量确定的；Pi(x,y)是某幅图像上各个像素点，x是图像宽度坐标，y是图像高度坐标。
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x,y)是经过offset校正后的亮场图I；Avg_I是图像I的所有像素点累加后的平均值。
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，C是一个预设的固定常数，通常是整数。
另外，需要说明的是，平板探测器在两次曝光之间有一个时间限定，需满足第二次曝光在Dark_n采集完成以后进行，即设置最短循环时间(Cycle time)。
步骤3)S13，进行第二次曝光采集到第二目标图像，命名为Light_2，并同时记录采集时间点T2、曝光窗口(Exposure window)的大小D2，并按照Light_1的处理方式处理Light_2，获得Target_image2；其中，所述第二目标图像与第一目标图像为相邻的两帧图像。
当采集Light_2完成后，平板探测器会自动采集原始暗场图像n帧，命名为Dark_1、Dark_2、……、Dark_n，并根据步骤2)的处理方式，获得Ghost_template2；
步骤4)S14，根据图像残影的自然衰减规律，利用T1和T2之间的时间间隔以及D1和D2的变化，设置加权系数，作用到Ghost_template1，用以定量对比图像到目标图像之间的图像残影衰减量；
具体地，需提前获得图像残影衰减系数K，本实施例根据拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1,T2)是一个自变量是T1、T2的函数；
其中，I0(x,y)是某一(T1,T2)点下经过offset校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值；
其中，I1(x,y)是对应某一(T1,T2)点的经过offset校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
另外，因前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，因此，需根据曝光窗口(Exposure window)设置漏电流补偿系数Kd，本发明根据拟合各个曝光窗口(Exposure window)下的图像残影信息情况，确定其符合的公式：
其中，W(D1,D2)是自变量为D1、D2的函数；
其中，G(T1,T2)是某一(T1,T2)点下的理论计算值；
其中，I0(x,y)是某一(D1,D2)点下，对应某一(T1,T2)点的经过offset校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值；
其中，I1(x,y)是对应某一(D1,D2)点和(T1,T2)点的经过offset校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值；
步骤5)S15，对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x,y)是目标图像Target_image2上的各个像素点灰度值，PG1(x,y)是对比图像Ghost_template1上的各个像素点灰度值；
其中，K、Kd是根据T1、T2、D1、D2实时计算G(T1,T2)、W(D1,D2)获得；
其中，K1是微调系数。
如图2所示，本实施例还提供一种平板探测器，在本实施例中，所述平板探测器为非晶硅X射线非晶硅平板探测器，其包括：图像采集模块1、图像处理模块2、系统控制模块3以及图像显示模块4。
如图2所示，所述图像采集模块1用于采集图像。
在本实施例中，所述图像采集模块1包括闪烁体以及TFT面板11，以及信号读出、放大、模数转换子模块12。
所述闪烁体用于将X射线转换为可见光，所述TFT面板用于感应所述可见光，并将该可见光转换为电信号。
所述信号读出、放大、模数转换子模块12用于对TFT面板输出的电信号进行读出、放大后，进行模数转换，最后输出数字信号。
如图2所示，所述图像处理模块2用于实现目标图像和对比图像的图像校正工作，根据对比图像和目标图像之间的时间间隔以及曝光窗口大小，实现影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算，并通过减法公式，将图像残影从目标图像中减去，以终降低或者消除图像残影，获得最终的图像。
在本实施例中，所述图像处理模块2中，目标图像和对比图像的图像校正工作包括：offset校正、Gain校正及Defect校正。
具体地，所述offset校正包括：
依据采集的第一目标图像的时间点T1、曝光窗口的大小D1，以及采集的n帧原始暗场图像，采用算法：
其中，n根据第一目标图像的固定偏移量校正方式以及第一帧原始暗场图像的残影量确定；Pi(x,y)是某幅图像上各个像素点，x为图像的宽度坐标，y为图像的高度坐标；
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x,y)是图像I经过固定偏移量校正后的亮场图；Avg_I是图像I的所有像素点累加后的平均值；
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，Dark_1为第一帧原始暗场图像，Light_1为第一目标图像，C为预设的固定常数；
依据采集第二目标图像的时间点T2、曝光窗口的大小D2，以及采集的n帧原始暗场图像，依据上述的offset校正的算法，获得第二目标图像的Ghost_template2以及Target_image2，其中，所述第二目标图像与第一目标图像为相邻的两帧图像。
进一步地，所述图像处理模块2中，影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算包括：
获得加权系数K的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1,T2)是一个自变量是T1、T2的函数，I0(x,y)是某一(T1,T2)点下经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值；
由于前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，需根据曝光窗口设置漏电流补偿系数Kd，拟合各个曝光窗口下的图像残影信息情况，确定其符合的公式：
其中，W(D1,D2)是自变量为D1、D2的函数，G(T1,T2)是某一(T1,T2)点下的理论计算值，I0(x,y)是某一(D1,D2)点下，对应某一(T1,T2)点的经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x,y)的所有像素点累加后的平均值，I1(x,y)是对应某一(D1,D2)点和(T1,T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
进一步地，所述图像处理模块2中，通过减法公式，将图像残影从目标图像中减去，以降低或者消除图像残影，获得最终的图像的步骤包括：
对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x,y)是目标图像Target_image2上的各个像素点灰度值，PG1(x,y)是对比图像Ghost_template1上的各个像素点灰度值，K1是微调系数。
在本实施示例中，所述图像处理模块2包括预设图像基本校正子模块21、以及图像Ghost校正子模块22。
所述预设图像基本校正子模块21由ARM处理器实现，其用于实现目标图像和对比图像的校正预处理工作(包括offset校正、Gain校正及Defect校正)，然后根据对比图像和目标图像之间的时间间隔、曝光窗口(Exposure window)等变量，实现影响目标图像残影信息量的影响因子(图像残影衰减系数K以及玻璃漏电流补偿系数Kd)的函数计算。
所述图像Ghost校正子模块22则用于通过减法公式，将图像残影从目标图像中减去，以达到最终降低或者消除图像残影的目的。
如图2所示，在本实施例中，所述系统控制模块3包括显示系统、电源系统、通信控制系统。所述系统控制模块3用于实现图像采集模块1、图像处理模块2及图像显示模块4的控制。整个采集处理过程由系统控制模块3控制完成，例如，系统控制模块3控制图像采集模块1在目标图像之前动态的采集多帧对比图像。
图像显示模块4，用于显示最终的图像。
最终处理后图像需在显示模块进行呈现，通过图像处理模块2校正处理后的图像最终由图像显示模块4显示出来。
本发明的技术关键点是根据相邻两帧校正后图像上的残影形状一致的特点以及残影信息量具有自然衰减的规律的特征，在目标图像之前插入多帧的对比图像，将对比图像进行相应地等量转换后对目标图像进行减法运算，实现图像残影的最大化降低或者消除。
综上所述，本发明的平板探测器以及降低平板探测器图像残影的方法，具有以下有益效果：本发明解决的技术问题是在保证图像质量的同时，降低或者消除当前帧(第n次，且n>1)图像上残留的前几帧图像(第1～n-1次)的信息。通过在目标图像前插入多帧自动采集的对比图像的方式，利用校正后的对比图像具有相同性质的图像残影，对目标图像上含残影信息的各像素点的进行公式计算，实现在不增加现有硬件成本的基础上，降低或消除图像残影的目的，且该方法具备高度的集成性、灵活性和可观的市场应用价值。
本发明在确保X射线平板探测器功能的前提下，不改动现有的硬件设备，对图像数据进行计算机处理，不但实现降低或者消除图像残影的目的，而且增强了显示系统、通信控制系统与X射线平板探测器的融合度、集成性。在后续产品升级或者维护中，并不需要技术人员的现场服务，只需发布优化升级的软件包即可。
所以，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。根据本发明的技术方案或者构思进行相应的改动或者替换，比如对目标图像和对比图像的预处理方式(校正方式)的改动，对加权系数的产生采用类似的数学公式进行计算(如指数函数、高斯函数、高次非线性函数等)以及对连续三帧或者以上曝光图像基于本发明方式进行重复迭代计算，都应该属于本发明所附有的权利要求和保护范围。
任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
Claims (10)
1.一种降低平板探测器图像残影的方法，其特征在于，包括步骤：
1)采集第一目标图像，并记录采集时间点T1、曝光窗口的大小D1，采集原始暗场图像n帧，所述第一目标图像、原始暗场图像具有固定偏移量；
2)对第一目标图像进行固定偏移量校正，算法包括：
其中，n根据第一目标图像的固定偏移量校正方式以及第一帧原始暗场图像的残影量确定；Pi(x，y)是某幅图像上各个像素点，x为图像的宽度坐标，y为图像的高度坐标；
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x，y)是图像I经过固定偏移量校正后的亮场图；Avg_I是图像I的所有像素点累加后的平均值；
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，Dark_1为第一帧原始暗场图像，Light_1为第一目标图像，C为预设的固定常数；
3)采集第二目标图像，并记录采集时间点T2、曝光窗口的大小D2，采集原始暗场图像n帧，并根据步骤2)的固定偏移量校正的算法，获得Ghost_template2以及Target_image2；
4)根据图像残影的自然衰减规律，利用T1和T2之间的时间间隔以及D1和D2的变化，设置加权系数，作用到Ghost_template1，用以定量对比图像Ghost_template1到目标图像Target_image2之间的图像残影衰减量；
5)对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像。
2.根据权利要求1所述的降低平板探测器图像残影的方法，其特征在于：所述固定偏移量为由于平板探测器中的电子电路、ROIC、Gate IC及像素本身存在的干扰信息造成原始暗场图像、及第一目标图像存在的本底噪声。
3.根据权利要求1所述的降低平板探测器图像残影的方法，其特征在于：步骤2)中，还包括对原始暗场图像、及第一目标图像进行Gain校正及Defect校正的步骤，步骤3)中，还包括对原始暗场图像、及第二目标图像进行Gain校正及Defect校正的步骤。
4.根据权利要求1所述的降低平板探测器图像残影的方法，其特征在于：步骤4)中，获得加权系数K的方法为拟合各个(T1，T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1，T2)是一个自变量是T1、T2的函数，I0(x，y)是某一(T1，T2)点下经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x，y)的所有像素点累加后的平均值，I1(x，y)是对应某一(T1，T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
5.根据权利要求4所述的降低平板探测器图像残影的方法，其特征在于：步骤4)中，由于前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，需根据曝光窗口设置漏电流补偿系数Kd，拟合各个曝光窗口下的图像残影信息情况，确定其符合的公式：
其中，W(D1，D2)是自变量为D1、D2的函数，G(T1，T2)是某一(T1，T2)点下的理论计算值，I0(x，y)是某一(D1，D2)点下，对应某一(T1，T2)点的经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x，y)的所有像素点累加后的平均值，I1(x，y)是对应某一(D1，D2)点和(T1，T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
6.根据权利要求5所述的降低平板探测器图像残影的方法，其特征在于：步骤5)中，对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x，y)是目标图像Target_image2上的各个像素点灰度值，PG1(x，y)是对比图像Ghost_template1上的各个像素点灰度值，K1是微调系数。
7.一种平板探测器，其特征在于，包括：
图像采集模块，用于采集图像；
图像处理模块，用于实现目标图像和对比图像的图像校正工作，根据对比图像和目标图像之间的时间间隔以及曝光窗口大小，实现影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算，并通过减法公式，将图像残影从目标图像中减去，以终降低或者消除图像残影，获得最终的图像；
所述图像处理模块中，目标图像和对比图像的图像校正工作包括offset校正，所述offset校正包括：
依据采集的第一目标图像的时间点T1、曝光窗口的大小D1，以及采集的n帧原始暗场图像，采用算法：
其中，n根据第一目标图像的固定偏移量校正方式以及第一帧原始暗场图像的残影量确定；Pi(x，y)是某幅图像上各个像素点，x为图像的宽度坐标，y为图像的高度坐标；
Map_2(x，y)＝(I(x，y)-Avg_I)/Avg_I
其中，I(x，y)是图像I经过固定偏移量校正后的亮场图；Avg_I是图像I的所有像素点累加后的平均值；
Ghost_template1＝(Dark_1-Map_1+C)/(1+Map_2)
Target_image1＝(Light_1-Map_1+C)/(1+Map_2)
其中，Dark_1为第一帧原始暗场图像，Light_1为第一目标图像，C为预设的固定常数；
依据采集第二目标图像的时间点T2、曝光窗口的大小D2，以及采集的n帧原始暗场图像，依据上述的offset校正的算法，获得第二目标图像的Ghost_template2以及Target_image2；
系统控制模块，用于实现图像采集模块、图像处理模块及图像显示模块的控制；
图像显示模块，用于显示最终的图像。
8.根据权利要求7所述的平板探测器，其特征在于：所述图像处理模块中，目标图像和对比图像的图像校正工作包括：offset校正、Gain校正及Defect校正。
9.根据权利要求7所述的平板探测器，其特征在于：所述图像处理模块中，影响目标图像残影信息量的图像残影衰减系数K以及玻璃漏电流补偿系数Kd的函数计算包括：
获得加权系数K的方法为拟合各个(T1，T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
其中，G(T1，T2)是一个自变量是T1、T2的函数，I0(x，y)是某一(T1，T2)点下经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是10(x，y)的所有像素点累加后的平均值，I1(x，y)是对应某一(T1，T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值；
由于前后二次曝光窗口的差异会引起玻璃漏电流的变化，从而改变图像残影的大小，需根据曝光窗口设置漏电流补偿系数Kd，拟合各个曝光窗口下的图像残影信息情况，确定其符合的公式：
其中，W(D1，D2)是自变量为D1、D2的函数，G(T1，T2)是某一(T1，T2)点下的理论计算值，I0(x，y)是某一(D1，D2)点下，对应某一(T1，T2)点的经过固定偏移量校正后的目标图像的各个像素点的灰度值，Avg_I0是I0(x，y)的所有像素点累加后的平均值，I1(x，y)是对应某一(D1，D2)点和(T1，T2)点的经过固定偏移量校正后的对比图像的各个像素点的灰度值，Avg_I1是I1的所有像素点累加后的平均值。
10.根据权利要求9所述的平板探测器，其特征在于：所述图像处理模块中，通过减法公式，将图像残影从目标图像中减去，以终降低或者消除图像残影，获得最终的图像的步骤包括：
对目标图像Target_image2和对比图像Ghost_template1进行对应像素点的减法运算，获得最终校正图像：
Corrected_image＝PT2(x，y)-K1*K*Kd*(PG1(x，y)-C)
其中，PT2(x，y)是目标图像Target_image2上的各个像素点灰度值，PG1(x，y)是对比图像Ghost_templatel上的各个像素点灰度值，K1是微调系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510309168.5A CN106296595B (zh) | 2015-06-08 | 2015-06-08 | 一种平板探测器以及降低平板探测器图像残影的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510309168.5A CN106296595B (zh) | 2015-06-08 | 2015-06-08 | 一种平板探测器以及降低平板探测器图像残影的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106296595A CN106296595A (zh) | 2017-01-04 |
| CN106296595B true CN106296595B (zh) | 2019-01-01 |
Family
ID=57659770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510309168.5A Active CN106296595B (zh) | 2015-06-08 | 2015-06-08 | 一种平板探测器以及降低平板探测器图像残影的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106296595B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171765B (zh) * | 2017-12-08 | 2021-06-18 | 上海奕瑞光电子科技股份有限公司 | 平板探测器残影校正方法及其校正装置 |
| CN108416817B (zh) * | 2018-02-01 | 2022-04-19 | 上海奕瑞光电子科技股份有限公司 | 一种残影校正系数自动获得方法 |
| CN111445397B (zh) * | 2020-03-09 | 2023-07-04 | 东软医疗系统股份有限公司 | 平板探测器残影校正方法、装置、存储介质及医疗设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102426376A (zh) * | 2011-08-11 | 2012-04-25 | 西北工业大学 | 一种平板探测器的监控与校正方法 |
| CN102567954A (zh) * | 2010-12-30 | 2012-07-11 | 深圳迈瑞生物医疗电子股份有限公司 | 平板探测器噪声抑制方法及装置 |
| CN103091968A (zh) * | 2013-01-17 | 2013-05-08 | 浙江康源医疗器械有限公司 | 具有自动校正功能的x-射线成像方法 |
| CN103377464A (zh) * | 2012-04-28 | 2013-10-30 | 北京国药恒瑞美联信息技术有限公司 | 一种消除残影的图像处理方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5405093B2 (ja) * | 2008-12-05 | 2014-02-05 | 富士フイルム株式会社 | 画像処理装置及び画像処理方法 |
-
2015
- 2015-06-08 CN CN201510309168.5A patent/CN106296595B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102567954A (zh) * | 2010-12-30 | 2012-07-11 | 深圳迈瑞生物医疗电子股份有限公司 | 平板探测器噪声抑制方法及装置 |
| CN102426376A (zh) * | 2011-08-11 | 2012-04-25 | 西北工业大学 | 一种平板探测器的监控与校正方法 |
| CN103377464A (zh) * | 2012-04-28 | 2013-10-30 | 北京国药恒瑞美联信息技术有限公司 | 一种消除残影的图像处理方法及系统 |
| CN103091968A (zh) * | 2013-01-17 | 2013-05-08 | 浙江康源医疗器械有限公司 | 具有自动校正功能的x-射线成像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106296595A (zh) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106097282B (zh) | 消除平板探测器图像残影的方法及平板探测器 | |
| Newberry | Signal-to-noise considerations for sky-subtracted CCD data | |
| CN101739672A (zh) | 一种基于分区域插值的直方图均衡方法和装置 | |
| CN106296595B (zh) | 一种平板探测器以及降低平板探测器图像残影的方法 | |
| US9571743B2 (en) | Dynamic exposure adjusting method and electronic apparatus using the same | |
| CN105100632B (zh) | 成像设备自动曝光的调整方法及装置、成像设备 | |
| CN102156962B (zh) | 信息处理设备及其处理方法 | |
| CN108172659B (zh) | 平板探测器及其残影数据表的生成方法、残影补偿校正方法 | |
| CN108833775B (zh) | 一种抗运动鬼影的hdr方法、装置及便携式终端 | |
| JP2013255606A5 (zh) | ||
| CN110211056A (zh) | 基于局部中值直方图的自适应红外图像去条纹算法 | |
| RU2441281C1 (ru) | Способ оценки шума цифровых рентгенограмм | |
| Saunders Jr et al. | Resolution and noise measurements of five CRT and LCD medical displays | |
| CN102113891B (zh) | 放射线摄像装置 | |
| CN112288646A (zh) | 堆叠降噪方法、装置、电子设备及存储介质 | |
| CN108055486A (zh) | 一种亮度校正方法及装置 | |
| CN108171765B (zh) | 平板探测器残影校正方法及其校正装置 | |
| TWI661812B (zh) | 造影系統與造影方法 | |
| CN106252371B (zh) | 一种x射线图像传感器及校正图像干扰的方法 | |
| CN106526685A (zh) | 一种行李物品的安检方法和系统 | |
| CN107610072A (zh) | 一种基于梯度引导滤波的微光视频图像自适应降噪方法 | |
| JP2011151549A5 (zh) | ||
| JP2018201685A (ja) | 放射線撮影装置、放射線撮影方法およびプログラム | |
| CN106204505A (zh) | 基于红外图像增强直方图映射表的视频去闪动方法 | |
| CN106206637A (zh) | 一种x射线图像传感器及校正图像干扰的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 201201 Shanghai City, Pudong New Area Ruiqinglu 590 No. 9 Building 2 Room 202 Applicant after: Shanghai Yi Ruiguang electronic Polytron Technologies Inc Address before: 201201 Shanghai City, Pudong New Area Zhangjiang hi tech Industrial District Ruiqinglu No. 590 9 2 storey 202 room Applicant before: Shanghai Yirui Optoelectronics Technology Co., Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |