CN106097282B - 消除平板探测器图像残影的方法及平板探测器 - Google Patents
消除平板探测器图像残影的方法及平板探测器 Download PDFInfo
- Publication number
- CN106097282B CN106097282B CN201610513009.1A CN201610513009A CN106097282B CN 106097282 B CN106097282 B CN 106097282B CN 201610513009 A CN201610513009 A CN 201610513009A CN 106097282 B CN106097282 B CN 106097282B
- Authority
- CN
- China
- Prior art keywords
- image
- original
- flat panel
- acquisition
- panel detector
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 206010047571 Visual impairment Diseases 0.000 claims abstract description 73
- 238000003702 image correction Methods 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims description 29
- 238000003705 background correction Methods 0.000 claims description 27
- 238000012545 processing Methods 0.000 claims description 26
- 230000005611 electricity Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 abstract description 4
- 230000008030 elimination Effects 0.000 description 16
- 238000003379 elimination reaction Methods 0.000 description 16
- 230000003068 static effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 3
- 229910021417 amorphous silicon Inorganic materials 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 229910052711 selenium Inorganic materials 0.000 description 3
- 239000011669 selenium Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- XQPRBTXUXXVTKB-UHFFFAOYSA-M caesium iodide Chemical compound [I-].[Cs+] XQPRBTXUXXVTKB-UHFFFAOYSA-M 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229940075613 gadolinium oxide Drugs 0.000 description 1
- 229910001938 gadolinium oxide Inorganic materials 0.000 description 1
- CMIHHWBVHJVIGI-UHFFFAOYSA-N gadolinium(iii) oxide Chemical compound [O-2].[O-2].[O-2].[Gd+3].[Gd+3] CMIHHWBVHJVIGI-UHFFFAOYSA-N 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供一种平板探测器以及消除平板探测器图像残影的方法，方法包括：1）使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1；2）使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2；3）根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数；4）利用步骤3）中加权系数对步骤1）采集的原始暗场图像进行补偿后通过差分运算对步骤2）中采集的原始亮场图像进行图像残影校正。本发明的方法流程简明、算法精简，具有高度的实用性和系统可集成性。
Description
技术领域
本发明属于放射影像图像处理技术领域，特别是涉及一种平板探测器以及消除平板探测器图像残影的方法。
背景技术
当前，非晶硅或非晶硒X射线平板探测器，由于其本身的物理特性，在过饱和剂量曝光时，由于电荷无法一次性释放，目标图像会存在残影的问题，且一旦出现图像残影，无法从硬件上消除，需要时间进行自然衰减，后续曝光拍摄的图像，都会存在残影的干扰。
目前的非晶硅或非晶硒X射线平板探测器，会存在图像残影这一现象，是由闪烁体(如碘化铯、琉氧化钆)余辉及TFT sensor、photo Diode、本身RC以及材料(非晶硅、非晶硒)缺陷态等原因同时造成的，会严重影响临床诊断。目前的平板探测器厂商，大多只是对图像进行本底校正、增益校正和坏点、坏线校正而已，考虑问题都是基于单一因素的影响，并没有将平板探测器看作一个小系统，对图像残影进行校正。因此，由于图像lag信息的自然衰减规律，大剂量后曝光拍摄的图像，存在残影的可能非常大，这会导致图像失真，影响临床诊断。当前，大多数企业或科研单位，主要还是通过硬件上的改进，如重新设计电路，改进FPGA的工作时序，或者改进玻璃的工艺制成，以此来消除或者降低图像残影。但是，从目前来看，这些技术手段繁琐、复杂且风险高，代价昂贵，效果也不尽人意，无法做到图像无可见残影。
也有部分企业或研究机构从图像处理技术角度出发，利用前一帧图像的具体区域特征，进行参考图像获取，并设计理论算法，如专利CN103377464A，对目标图像和参考图像进行梯度域变换，经过差分运算后获得去除残影的最终梯度图像，最后进行逆变换，恢复得到时域图像。该方法，算法比较复杂，像素计算量大，实时性差，且无法准确消除目标图像与残影信息之间相交接的边缘残影，会形成明显的带状残影。其次，该算法只是针对相邻两帧图像之间的残影消除，并未考量多次过饱和剂量曝光下的图像残影消除。对多次、隔帧残影，算法并未具有很好的鲁棒性，导致校正结果会错乱，甚至引起过校正，图像有效信息被错误去除而引起的严重失真。
鉴于上述的分析概括，从软硬件技术角度来说，消除或者降低图像残影，具备极大的改进和提高空间。
发明内容
鉴于以上所述现有技术的缺点，本发明的目的在于提供一种平板探测器以及消除平板探测器图像残影的方法，用于解决现有技术中平板探测器图像残影难以有效消除的问题。
为实现上述目的及其他相关目的，本发明提供一种消除平板探测器图像残影的方法，包括步骤：
1)使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1；
2)使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2；
3)根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数；
4)利用步骤3)中加权系数对步骤1)采集的原始暗场图像进行补偿后通过差分运算对步骤2)中采集的原始亮场图像进行图像残影校正。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤1)及步骤2)中，均还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤3)中，获得加权系数K(t1，t2)的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
K(t1,t2)＝a*f(t1,t2)+b*g(t1,t2)
其中，a和b分别为常量系数，f(t1,t2)和g(t1,t2)为一对指数函数。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤4)中，差分运算公式为：
E1′＝E1-K(t1,t2)*L1+K(t1,t2)*Avg
其中，E1′为图像残影校正后的亮场图像，E1为步骤2)中采集的原始亮场图像，K(t1,t2)为加权系数L1为步骤1)中采集的原始暗场图像，Avg为基于L1的区域特性得到的L1中部分像素点的平均值或中值。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤4)之后还包括一下步骤：
5)将最近一次曝光采集的原始暗场图像进行迭代更新；
6)使用平板探测器再次曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点；
7)依据加权系数对步骤5)中迭代更新后的原始暗场图像进行补偿后通过差分运算对步骤6)中采集的原始亮场图像进行图像残影校正；
8)重复步骤5)至步骤7)直至得到最终的图像。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤5)中，迭代更新公式为：
其中，Ls′为更新后的原始暗场图像，Ls为最近一次曝光采集的原始暗场图像，K为最近一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔系数，G1为最近一次曝光采集的原始暗场图像与上一次曝光采集的原始暗场图像之间的时间间隔系数，G2为上一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤7)中，还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。
作为本发明的消除平板探测器图像残影的方法的一种优选方案，步骤8)中，差分运算公式为：
Es′＝Es-K(t1,t2)*Ls′+K(t1,t2)*Avg
其中，Es′为图像残影校正后的亮场图像，Es为步骤7)中采集的原始亮场图像，K(t1,t2)为加权系数，Ls′Ls为步骤5)中迭代更新后的原始暗场图像，Avg为基于Ls′的区域特性得到的Ls′中部分像素点的平均值或中值。
本发明还提供一种平板探测器，所述平板探测器包括：
图像采集模块，适于采集原始亮场图像及原始暗场图像；
图像处理模块，与所述图像采集模块相连接，适于对采集的原始暗场图像进行本底校正、增益校正、坏点/坏线校正、依据加权系数进行加权系数补偿及迭代更新，并对原始亮场图像进行本底校正、增益校正、坏点/坏线校正及图像残影校正；
图像显示模块，与所述图像处理模块相连接，适于显示所述图像处理模块处理后的最终的图像。
作为本发明的平板探测器的一种优选方案，所述图像采集模块包括：X射线单元、闪烁体及TFT面板以及PCB电路，所述X射线单元、所述闪烁体及TFT面板以及所述PCB电路依次相连接。
作为本发明的平板探测器的一种优选方案，所述图像采集模块包括：
基本校正单元，与所述图像采集模块相连接，适于对采集的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正；
图像残影校正单元，与所述基本校正单元相连接，适于将基本校正后的原始暗场图像依据加权系数进行加权系数补偿及迭代更新，并依据加权系数补偿及迭代更新后的所述原始暗场图像对原始亮场图像进行图像残影校正。
作为本发明的平板探测器的一种优选方案，所述平板探测器还包括系统控制模块，所述系统控制模块与所述图像采集模块、所述图像处理模块及所述图像显示模块相连接，适于实现对所述图像采集模块、所述图像处理模块及所述图像显示模块的控制。
如上所述，本发明的消除平板探测器图像残影的方法及平板探测器，具有以下有益效果：本发明的方法不对图像进行时域变换，也不对图像进行分割、截断或特征提取，只是在曝光的时候分别提取原始亮场图像及原始暗场图像，对原始暗场图像进行加权系数补偿、迭代更新处理后，将原始亮场图像与处理后的原始暗场图像进行差分运算即可消除平板探测器图像残影，整个方法流程简明、算法精简，具有高度的实用性和系统可集成性，具有客观的市场应用价值和临床实用性、适用性；同时，该方法不但适用于动态本底校正后的图像残影消除，也适用于静态本底校正后的图像残影消除，不但消除两帧之间的图像残影，还可以彻底消除多次隔帧的图像残影。适用于推广到各个类型平板探测器的使用，解决了静态本底校正中的图像残影问题，实现临床图像无可见残影。
附图说明
图1显示为本发明实施例一中提供的消除平板探测器图像残影的方法的流程图。
图2显示为本发明实施例一中提供的消除平板探测器图像残影的方法中的原始暗场图像迭代更新的示意图。
图3显示为采用本发明实施例一中提供的消除平板探测器图像残影的方法处理前后的对比图。
图4显示为本发明实施例二中提供的平板探测器的结构框图。
元件标号说明
1 图像采集模块
11 X射线单元
12 闪烁体及TFT面板
13 PCB电路
2 图像处理模块
21 基本校正单元
22 图像残影校正单元
3 图像显示模块
4 系统控制模块
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。
请参阅图1～图4。需要说明的是，本实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
实施例一
如图1所示，本实施例提供一种消除平板探测器图像残影的方法，从成像角度来说，所谓的图像残影是指余辉增殖效应，是采集两张或两张以上原始亮场图时，当前原始亮场图像中会残留有前一次或多次的图像信息，会出现的一种伪影现象。所述消除平板探测器图像残影的方法包括以下步骤：
1)使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1；
2)使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2；
3)根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数；
4)利用步骤3)中加权系数对步骤1)采集的原始暗场图像进行补偿后通过差分运算对步骤2)中采集的原始亮场图像进行图像残影校正。
在步骤1)中，请参阅图1中的S1步骤，使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1。
作为示例，使用平板探测器进行曝光之前，初始化所述平板探测器。为了便于区分，可以将该步骤中采集的原始亮场图像记为E0，将该步骤中采集的原始暗场图像记为L1，L1中含有E0的图像残影(lag)信息。
作为示例，采集一帧所述原始亮场图像及一帧所述原始暗场图像之后，还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。上对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正均为图像的基本校正，为本领域人员所熟知的技术，这里不再累述。
需要说明的是，由于步骤1)为初始平板探测器，即曝光之前并无图像残影信息存在，因此，该步骤中采集的原始亮场图像E0无需进行图像残影校正，即可在图像显示模块中予以显示。
在步骤2)中，请参阅图1中的S2步骤，使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2。
作为示例，为了便于区分，可以将该步骤中采集的原始亮场图像记为E1，将该步骤中采集的原始暗场图像记为L2。
作为示例，采集一帧所述原始亮场图像及一帧所述原始暗场图像之后，还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。上对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正均为图像的基本校正，为本领域人员所熟知的技术，这里不再累述。
在步骤3)中，请参阅图1中的S3步骤，根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数。
作为示例，根据t1与t2的时间间隔，基于图像残影自然衰减物理特性曲线，获得加权系数K(t1，t2)的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
K(t1,t2)＝a*f(t1,t2)+b*g(t1,t2)
其中，a和b分别为常量系数，f(t1,t2)和g(t1,t2)为一对指数函数，即K(t1，t2)为符合双指数衰减的物理曲线。
在步骤4)中，请参阅图1中的S4步骤，利用步骤3)中加权系数对步骤1)采集的原始暗场图像进行补偿后通过差分运算对步骤2)中采集的原始亮场图像进行图像残影校正。
作为示例，该步骤中，差分运算公式为：
E1′＝E1-K(t1,t2)*L1+K(t1,t2)*Avg
其中，E1′为图像残影校正后的亮场图像，E1为步骤2)中采集的原始亮场图像，K(t1,t2)为加权系数，L1为步骤1)中采集的原始暗场图像，Avg为基于L1的区域特性得到的L1中部分像素点的平均值或中值，用于步骤2)中采集的原始亮场图像E1的整体灰度值恢复。
当对初始化后的所述平板探测器进行两次曝光后的图像进行图像残影校正时，到步骤4)校正过程结束，即可得到最终需要的图像。但在实际应用中，需要使用平板探测器一直进行曝光拍摄，即需要进行多次曝光，此时，步骤4)之后，还应包括如下步骤：
5)将最近一次曝光采集的原始暗场图像进行迭代更新；
6)使用平板探测器再次曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点；
7)依据加权系数对步骤5)中迭代更新后的原始暗场图像进行补偿后通过差分运算对步骤6)中采集的原始亮场图像进行图像残影校正；
8)重复步骤5)至步骤7)直至得到最终的图像。
作为示例，在步骤5)中，由于图像残影的自然衰减特性会导致第二曝光及之后的每次曝光采集的原始暗场图像中含有的图像残影信息不全，所以执行步骤5)是对后续曝光的图像进行图像残影消除的关键步骤，可以抵消后续需要图像残影校正的原始亮场图像与原始暗场图像之间的图像残影信息不对等的情况。
作为示例，步骤5)中，迭代更新公式为：
其中，Ls′为更新后的原始暗场图像，Ls为最近一次曝光采集的原始暗场图像，K为最近一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔系数，即K为图2中第二次曝光采集的原始暗场图像L2与第三次曝光采集的原始亮场图像E2之间的时间间隔系数；G1为最近一次曝光采集的原始暗场图像与上一次曝光采集的原始暗场图像之间的时间间隔系数，即G1为图2中第一次曝光采集的原始暗场图像L1与第二次曝光采集的暗场图像L2之间的时间间隔系数；G2为上一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔系数，即G2为第一次曝光采集的原始暗场图像L1与第三次曝光采集的原始亮场图像E2之间的时间间隔系数。
需要说明的是，当步骤5)承接步骤4)时，所述“最近一次曝光”即为步骤2)中对应的曝光。
需要进一步说明的是，迭代更新公式中，K为对应时间间隔的加权系数，G1及G2为基于平板探测器在饱和剂量下的空曝光，进行等间隔的多帧原始暗场图像的采集获得的基于对应时间间隔的离散G数值拟合构造出的函数。
作为示例，步骤6)中，采集一帧所述原始亮场图像及一帧所述原始暗场图像之后，还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。上对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正均为图像的基本校正，为本领域人员所熟知的技术，这里不再累述。
作为示例，步骤7)中，差分运算公式为：
Es′＝Es-K(t1,t2)*Ls′+K(t1,t2)*Avg其中，Es′为图像残影校正后的亮场图像，Es为步骤7)中采集的原始亮场图像，K(t1,t2)为加权系数，Ls′为步骤5)中迭代更新后的原始暗场图像，Avg为基于Ls′的区域特性得到的Ls′中部分像素点的平均值或中值。
作为示例，步骤8)中，可以根据实际需要设定重复步骤5)至步骤7)的次数。
需要说明的是，步骤8)中，第一次重复步骤5)至步骤7)中，重复步骤5)中对应所述的“最近一次曝光”为对应步骤6)中的曝光。以后的重复循环过程中，依次类推，此次不在穷举。
图3为采用本实施例中所述的消除平板探测器图像残影的方法对图像进行处理前后的对比图，其中图(a)为处理前的图像，(b)为处理后的图像，由图3中图(a)及图(b)可以明显看出：在处理前，图(a)中的图像中存在明显的模糊的图像残影，处理后，图(b)中的图像非常清晰，对应于图(a)中的图像残影已被完全去除。
本发明的方法不对图像进行时域变换，也不对图像进行分割、截断或特征提取，只是在曝光的时候分别提取原始亮场图像及原始暗场图像，对原始暗场图像进行加权系数补偿、迭代更新处理后，将原始亮场图像与处理后的原始暗场图像进行差分运算即可消除平板探测器图像残影，整个方法流程简明、算法精简，具有高度的实用性和系统可集成性，具有客观的市场应用价值和临床实用性、适用性；同时，该方法不但适用于动态本底校正后的图像残影消除，也适用于静态本底校正后的图像残影消除，不但消除两帧之间的图像残影，还可以彻底消除多次隔帧的图像残影。适用于推广到各个类型平板探测器的使用，解决了静态本底校正中的图像残影问题，实现临床图像无可见残影。
实施例二
请参阅图4，本发明还提供一种平板探测器，所述平板探测器包括：图像采集模块1，所述图像采集模块1适于采集原始亮场图像及原始暗场图像；图像处理模块2，所述图像处理模块2与所述图像采集模块1相连接，适于采用实施例一中所述的方法对采集的原始暗场图像进行本底校正、增益校正、坏点/坏线校正、依据加权系数进行加权系数补偿及迭代更新，并对原始亮场图像进行本底校正、增益校正、坏点/坏线校正及图像残影校正；图像显示模块3，所述图像显示模块3与所述图像处理模块2相连接，适于显示所述图像处理模块2处理后的最终的图像。
作为示例，所述图像采集模块1包括：X射线单元11、闪烁体及TFT面板12以及PCB电路13，所述X射线单元11、所述闪烁体及TFT面板12以及所述PCB电路13依次相连接。
作为示例，所述图像采集模块2包括：基本校正单元21，所述基本校正单元21与所述图像采集模块1相连接，适于对采集的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正；图像残影校正单元22，所述图像残影校正单元22与所述基本校正单元21相连接，适于采用实施例一所述的方法将基本校正后的原始暗场图像依据加权系数进行加权系数补偿及迭代更新，并依据加权系数补偿及迭代更新后的所述原始暗场图像对原始亮场图像进行图像残影校正。
作为示例，所述平板探测器还包括系统控制模块3，所述系统控制模块3与所述图像采集模块1、所述图像处理模块2及所述图像显示模块3相连接，适于实现对所述图像采集模块1、所述图像处理模块2及所述图像显示模块3的控制。
综上所述，本发明的消除平板探测器图像残影的方法及平板探测器，所述消除平板探测器图像残影的方法包括以下步骤：1)使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1；2)使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2；3)根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数；4)利用步骤3)中加权系数对步骤1)采集的原始暗场图像进行补偿后通过差分运算对步骤2)中采集的原始亮场图像进行图像残影校正。本发明的方法不对图像进行时域变换，也不对图像进行分割、截断或特征提取，只是在曝光的时候分别提取原始亮场图像及原始暗场图像，对原始暗场图像进行加权系数补偿、迭代更新处理后，将原始亮场图像与处理后的原始暗场图像进行差分运算即可消除平板探测器图像残影，整个方法流程简明、算法精简，具有高度的实用性和系统可集成性，具有客观的市场应用价值和临床实用性、适用性；同时，该方法不但适用于动态本底校正后的图像残影消除，也适用于静态本底校正后的图像残影消除，不但消除两帧之间的图像残影，还可以彻底消除多次隔帧的图像残影。适用于推广到各个类型平板探测器的使用，解决了静态本底校正中的图像残影问题，实现临床图像无可见残影。
所以，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。根据本发明的技术方案或者构思进行相应的改动或者替换，比如对目标图像和对比图像的预处理方式(校正方式)的改动，对加权系数的产生采用类似的数学公式进行计算(如指数函数、高斯函数、高次非线性函数等)以及对连续三帧或者以上曝光图像基于本发明方式进行重复迭代计算，都应该属于本发明所附有的权利要求和保护范围。
任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
Claims (11)
1.一种消除平板探测器图像残影的方法，其特征在于，包括步骤：
1)使用平板探测器进行曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t1；
2)使用平板探测器继续曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点t2；
3)根据图像残影的自然衰减规律，结合t1与t2之间的时间间隔获取加权系数；获得加权系数K(t1，t2)的方法为拟合各个(T1,T2)点的图像残影信息情况，确定其自然衰减规律符合的公式：
K(t1,t2)＝a*f(t1,t2)+b*g(t1,t2)
其中，a和b分别为常量系数，f(t1,t2)和g(t1,t2)为一对指数函数；
4)利用步骤3)中加权系数对步骤1)采集的原始暗场图像进行补偿后通过差分运算对步骤2)中采集的原始亮场图像进行图像残影校正。
2.根据权利要求1所述的消除平板探测器图像残影的方法，其特征在于：步骤1)及步骤2)中，均还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。
3.根据权利要求1所述的消除平板探测器图像残影的方法，其特征在于：步骤4)中，差分运算公式为：
E1′＝E1-K(t1,t2)*L1+K(t1,t2)*Avg
其中，E1′为图像残影校正后的亮场图像，E1为步骤2)中采集的原始亮场图像，K(t1,t2)为加权系数，L1为步骤1)中采集的原始暗场图像，Avg为基于L1的区域特性得到的L1中部分像素点的平均值或中值。
4.根据权利要求1至3中任一项所述的消除平板探测器图像残影的方法，其特征在于：步骤4)之后还包括一下步骤：
5)将最近一次曝光采集的原始暗场图像进行迭代更新；
6)使用平板探测器再次曝光，采集一帧原始亮场图像及一帧原始暗场图像，并记录所述原始亮场图像的采集时间点；
7)依据加权系数对步骤5)中迭代更新后的原始暗场图像进行补偿后通过差分运算对步骤6)中采集的原始亮场图像进行图像残影校正；
8)重复步骤5)至步骤7)直至得到最终的图像。
5.根据权利要求4所述的消除平板探测器图像残影的方法，其特征在于：步骤5)中，迭代更新公式为：
其中，Ls'为更新后的原始暗场图像，Ls为最近一次曝光采集的原始暗场图像，K为最近一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔系数，G1为最近一次曝光采集的原始暗场图像与上一次曝光采集的原始暗场图像之间的时间间隔系数，G2为上一次曝光采集的原始暗场图像与下一次曝光采集的原始亮场图像之间的时间间隔系数。
6.根据权利要求4所述的消除平板探测器图像残影的方法，其特征在于：步骤7)中，还包括对采集到的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正的步骤。
7.根据权利要求4所述的消除平板探测器图像残影的方法，其特征在于：步骤8)中，差分运算公式为：
Es′＝Es-K(t1,t2)*Ls′+K(t1,t2)*Avg
其中，Es'为图像残影校正后的亮场图像，Es为步骤7)中采集的原始亮场图像，K(t1,t2)为加权系数，Ls'为步骤5)中迭代更新后的原始暗场图像，Avg为基于Ls'的区域特性得到的Ls'中部分像素点的平均值或中值。
8.一种平板探测器，其特征在于，所述平板探测器包括：
图像采集模块，适于采集原始亮场图像及原始暗场图像；
图像处理模块，与所述图像采集模块相连接，适于对采集的原始暗场图像进行本底校正、增益校正、坏点/坏线校正、依据加权系数进行加权系数补偿及迭代更新，并对原始亮场图像进行本底校正、增益校正、坏点/坏线校正及图像残影校正；
图像显示模块，与所述图像处理模块相连接，适于显示所述图像处理模块处理后的最终的图像。
9.根据权利要求8所述的平板探测器，其特征在于：所述图像采集模块包括：X射线单元、闪烁体及TFT面板以及PCB电路，所述X射线单元、所述闪烁体及TFT面板以及所述PCB电路依次相连接。
10.根据权利要求8所述的平板探测器，其特征在于：所述图像采集模块包括：
基本校正单元，与所述图像采集模块相连接，适于对采集的原始亮场图像及原始暗场图像进行本底校正、增益校正及坏点/坏线校正；
图像残影校正单元，与所述基本校正单元相连接，适于将基本校正后的原始暗场图像依据加权系数进行加权系数补偿及迭代更新，并依据加权系数补偿及迭代更新后的所述原始暗场图像对原始亮场图像进行图像残影校正。
11.根据权利要求8至10中任一项所述的平板探测器，其特征在于：所述平板探测器还包括系统控制模块，所述系统控制模块与所述图像采集模块、所述图像处理模块及所述图像显示模块相连接，适于实现对所述图像采集模块、所述图像处理模块及所述图像显示模块的控制。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610513009.1A CN106097282B (zh) | 2016-07-01 | 2016-07-01 | 消除平板探测器图像残影的方法及平板探测器 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610513009.1A CN106097282B (zh) | 2016-07-01 | 2016-07-01 | 消除平板探测器图像残影的方法及平板探测器 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106097282A CN106097282A (zh) | 2016-11-09 |
| CN106097282B true CN106097282B (zh) | 2019-02-26 |
Family
ID=57212386
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610513009.1A Active CN106097282B (zh) | 2016-07-01 | 2016-07-01 | 消除平板探测器图像残影的方法及平板探测器 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106097282B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171765B (zh) * | 2017-12-08 | 2021-06-18 | 上海奕瑞光电子科技股份有限公司 | 平板探测器残影校正方法及其校正装置 |
| CN108172659B (zh) * | 2017-12-20 | 2019-08-09 | 上海奕瑞光电子科技股份有限公司 | 平板探测器及其残影数据表的生成方法、残影补偿校正方法 |
| CN108416817B (zh) * | 2018-02-01 | 2022-04-19 | 上海奕瑞光电子科技股份有限公司 | 一种残影校正系数自动获得方法 |
| CN109408953A (zh) * | 2018-10-23 | 2019-03-01 | 上海奕瑞光电子科技股份有限公司 | 一种适用于连续可变帧率的本底值补偿方法 |
| CN109712084B (zh) * | 2018-12-10 | 2021-01-19 | 上海奕瑞光电子科技股份有限公司 | 一种图像修复方法、图像修复系统及平板探测器 |
| CN109803100A (zh) * | 2019-01-25 | 2019-05-24 | 冠捷显示科技（厦门）有限公司 | 一种自适应消残影方法 |
| JP7200816B2 (ja) * | 2019-04-22 | 2023-01-10 | 株式会社島津製作所 | Ｘ線位相イメージング装置およびｘ線位相コントラスト画像生成方法 |
| CN110161555A (zh) * | 2019-05-28 | 2019-08-23 | 深圳市菲森科技有限公司 | 一种x射线探测器校正方法、装置、系统及存储介质 |
| CN110933256B (zh) * | 2019-12-03 | 2022-05-03 | 上海奕瑞光电子科技股份有限公司 | 图像暗场漏电流的校正方法、装置、电子终端、存储介质 |
| CN111445397B (zh) * | 2020-03-09 | 2023-07-04 | 东软医疗系统股份有限公司 | 平板探测器残影校正方法、装置、存储介质及医疗设备 |
| CN113447485A (zh) * | 2020-03-26 | 2021-09-28 | 捷普电子(新加坡)公司 | 光学检测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101510298A (zh) * | 2009-03-17 | 2009-08-19 | 西北工业大学 | 一种ct伪影的综合校正方法 |
| CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
| CN104349083A (zh) * | 2013-07-25 | 2015-02-11 | 南京普爱射线影像设备有限公司 | 一种平板探测器图像拖尾的递归校正方法 |
| CN104749648A (zh) * | 2013-12-27 | 2015-07-01 | 清华大学 | 多能谱静态ct设备 |
| CN104835125A (zh) * | 2015-05-08 | 2015-08-12 | 上海联影医疗科技有限公司 | 平板探测器的坏点校正方法 |
-
2016
- 2016-07-01 CN CN201610513009.1A patent/CN106097282B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101510298A (zh) * | 2009-03-17 | 2009-08-19 | 西北工业大学 | 一种ct伪影的综合校正方法 |
| CN104349083A (zh) * | 2013-07-25 | 2015-02-11 | 南京普爱射线影像设备有限公司 | 一种平板探测器图像拖尾的递归校正方法 |
| CN104749648A (zh) * | 2013-12-27 | 2015-07-01 | 清华大学 | 多能谱静态ct设备 |
| CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
| CN104835125A (zh) * | 2015-05-08 | 2015-08-12 | 上海联影医疗科技有限公司 | 平板探测器的坏点校正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 非晶硅平板探测器DR成像校正方法;闵吉磊 等;《失效分析与预防》;20150630;第10卷(第3期);第135页右栏，第136页左栏，图8 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106097282A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106097282B (zh) | 消除平板探测器图像残影的方法及平板探测器 | |
| CN101779961B (zh) | X射线图像诊断装置及x射线图像诊断装置的控制方法 | |
| JP5315157B2 (ja) | 情報処理装置、ライン状ノイズ低減処理方法、及びプログラム | |
| CN105100632B (zh) | 成像设备自动曝光的调整方法及装置、成像设备 | |
| CN106231210B (zh) | 一种x射线图像传感器及消除图像残影的方法 | |
| CN108172659B (zh) | 平板探测器及其残影数据表的生成方法、残影补偿校正方法 | |
| CN108833775B (zh) | 一种抗运动鬼影的hdr方法、装置及便携式终端 | |
| CN101571947A (zh) | 图像处理设备和图像处理方法 | |
| CN108596861A (zh) | 一种基于深度残差学习的ct图像金属伪影去除方法 | |
| CN102007759A (zh) | 光或放射线摄像装置 | |
| CN109709597B (zh) | 平板探测器的增益校正方法 | |
| CN106296595B (zh) | 一种平板探测器以及降低平板探测器图像残影的方法 | |
| CN105631819A (zh) | 一种ccd dr探测器的平场矫正方法及系统 | |
| CN108171765B (zh) | 平板探测器残影校正方法及其校正装置 | |
| CN106308834A (zh) | 一种x射线图像传感器及其消除图像残影的方法 | |
| JP4210464B2 (ja) | Ｘ線画像診断装置 | |
| CN106129078B (zh) | 一种x射线图像传感器及其校正图像干扰的方法 | |
| CN108110015A (zh) | X射线图像传感器、平板探测器及其图像采集校正方法 | |
| CN106206637B (zh) | 一种x射线图像传感器及校正图像干扰的方法 | |
| CN108596993A (zh) | 校正图像不饱和伪影的系统及校正方法 | |
| CN104394334B (zh) | 图像传感器在轨噪声智能抑制方法 | |
| Nguyen | Digital Radiography with a Consumer Camera: Image Denoising and Deblurring | |
| CN105962961B (zh) | 一种抑制饱和带伪影方法及系统 | |
| CN109408953A (zh) | 一种适用于连续可变帧率的本底值补偿方法 | |
| CN106341624B (zh) | 用于校正温度和漏电流的偏置模板的生成方法 |
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 Zhangjiang hi tech Industrial District Ruiqinglu No. 590 9 2 storey 202 room 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 |