一种基于最小熵的固定模式噪声去除方法技术

本发明专利技术公开的一种基于最小熵的固定模式噪声去除方法，属于图像处理领域。针对现有图像降噪方法在处理固定模式噪声时的不足，本发明专利技术实现方法为：首先将固定模式噪声建模为标定得到的噪声模式与一未知尺度参数的乘积，有效处理非高斯、非平稳的噪声；然后，利用最小熵原理自动估计尺度参数，保证估计的优良性；最后，从带噪声的图像中减去被尺度参数缩放后的噪声模式，得到去噪后的图像，达到有效去除图像中的固定模式噪声的目的。本发明专利技术适用图像处理等领域，能够有效校正图像中的固定模式噪声，极大改善暗光条件下的成像质量，具有较高应用价值。价值。价值。

【技术实现步骤摘要】
一种基于最小熵的固定模式噪声去除方法


[0001]本专利技术涉及一种基于最小熵的固定模式噪声去除方法，属于可见光图像处理领域。

技术介绍

[0002]受加工工艺限制，成像探测器中广泛存在像元间的不均匀，这会导致在最终的成像中存在固定模式噪声。所谓的固定模式噪声可以认为是一种特殊的“噪声”。对于给定的一个探测器来说，这种“噪声”的模式是固定不变的，强度随环境温度和曝光时间变化(关于固定模式噪声的详细介绍，见European Machine Vision Association.Standard for characterization of image sensors and cameras[J].EMVA Standard,2010,1288)。
[0003]固定模式噪声的存在极大地限制了相机在暗光下成像的能力。不同于常见的随机噪声，固定模式噪声在一次曝光内是基本恒定的，所以在一次曝光中固定模式噪声会随信号一起以相同的速度累积。增加曝光时间只会同时放大信号和噪声，不会提高最终图像的信噪比。在暗光条件下，场景的辐射强度比较低，所以最终的成像会被固定模式噪声严重污染。
[0004]现有的图像去噪方法在处理固定模式噪声时很难有好的表现。这是因为固定模式噪声往往伴随“过热像元”同时出现，这些“过热像元”导致固定模式噪声的分布函数是长尾的，不符合一般的高斯噪声假设。另外，因为加工工艺的缺陷，固定模式噪声在空间上往往存在低频的、不规则的变化，这导致固定模式噪声通常是非平稳的，不符合独立同分布的假设。非高斯、非平稳噪声的去除到目前为止依然是一个巨大的挑战。

技术实现思路

[0005]针对现有图像降噪方法在处理固定模式噪声时的不足，本专利技术旨在提供一种基于最小熵的固定模式噪声去除方法，达到去除图像中固定模式噪声的目的。
[0006]为达到以上目的，本专利技术采用以下技术方案：
[0007]本专利技术公开的一种基于最小熵的固定模式噪声去除方法，基于探测器噪声的物理特性，将固定模式噪声建模为噪声模式和一尺度参数的乘积，有效处理非高斯、非平稳的噪声；然后利用最小熵原理自动估计尺度参数，保证估计的优良性。从带噪声的图像中减去被尺度参数缩放后的噪声模式之后，得到去噪后的图像，达到有效去除图像中的固定模式噪声的目的。
[0008]本专利技术公开的一种基于最小熵的固定模式噪声去除方法，包含以下步骤：
[0009]步骤一：通过暗场标定得到噪声模式F。
[0010]所述F为一m
×
n矩阵，其中m为图像高度，n为图像宽度。具体的标定方式为：在暗室中，将相机镜头盖盖上，然后采集图像，采集到的图像即为噪声模式F。为了抑制随机噪声，采集多张图像，然后平均，最终得到噪声模式F。
[0011]步骤二：读取带噪声的图像Y，图像Y为一m
×
n矩阵。
[0012]步骤三：初始化尺度参数α＝0。
[0013]步骤四：定义辅助变量Z。
[0014]所述辅助变量Z的定义为
[0015]Z＝Y
‑
α
·
F,
ꢀꢀꢀꢀ
(1)
[0016]显然Z也为一m
×
n矩阵。
[0017]步骤五：更新α。
[0018]所述α按如下方式进行更新
[0019][0020]式中δ为一实数，控制梯度下降的步长，取值反比于噪声模式的标准差；g(α)为Z的熵。
[0021]步骤六：重复步骤四
‑
五，直到α收敛。
[0022]记步骤五更新前的尺度参数为α
old
，更新后的尺度参数为α
new
，α收敛的判据为|α
new
‑
α
old
|<∈，其中∈为一小量。
[0023]步骤七：减去固定模式噪声得到去噪后的图像X，达到去除图像中固定模式噪声的目的。
[0024]所述X由下式计算得到
[0025]X＝Y
‑
α
·
F,
ꢀꢀꢀꢀ
(3)
[0026]显然，X也为一个m
×
n矩阵。
[0027]有益效果：
[0028]1、本专利技术公开的一种基于最小熵的固定模式噪声去除方法，将固定模式噪声建模为标定得到的噪声模式与一尺度参数的乘积，能够有效处理非高斯、非平稳的噪声。
[0029]2、本专利技术公开的一种基于最小熵的固定模式噪声去除方法，使用熵最小化方法估计未知的尺度参数，保证估计的优良性。
[0030]3、本专利技术公开的一种基于最小熵的固定模式噪声去除方法，在实现有益效果1、2的基础上，能够有效去除图像中的固定模式噪声，极大改善暗光条件下成像的质量。
附图说明
[0031]图1是本专利技术公开的一种基于最小熵原理的固定模式噪声去除方法的流程图。
[0032]图2是使用本专利技术固定模式噪声去除方法去噪前后图像视觉质量对比。
[0033]其中(a)(c)为去除固定模式噪声之前的图像；(b)(d)为去除固定模式噪声之后的图像。
具体实施方式
[0034]下面将结合附图和实施例对本专利技术加以详细说明。同时也叙述了本专利技术技术方案解决的技术问题及有益效果，需要指出的是，所描述的实施例仅旨在便于对本专利技术的理解，而对其不起任何限定作用。
[0035]实施例1：
[0036]本实施例公开的一种基于最小熵的固定模式噪声去除方法，针对现有图像去噪方
法在处理固定模式噪声方面的不足。本实验的硬件测试条件为：Inter i7 6700，8G RAM，Matlab 2016a；所用的相机为PointGrey灰度相机，曝光时间设置为300ms。δ设置为0.1。
[0037]本实施例首先将固定模式噪声建模为噪声模式和一尺度参数的乘积，然后利用最小熵原理自动估计尺度参数。从带噪声图像中减去被尺度参数缩放后的噪声模式即可有效去除固定模式噪声。本实施例的流程图如图1所示。
[0038]本实施例公开的一种基于最小熵原理的固定模式噪声去除方法，包含以下步骤：
[0039]步骤一：通过暗场标定得到噪声模式F。
[0040]所述F为一512
×
640矩阵，其中512为图像高度，640为图像宽度。具体的标定方式为：在暗室中，将相机镜头盖盖上，然后采集图像，采集到的图像即为噪声模式F。为了抑制随机噪声，采集100张图像，然后平均，最终得到噪声模式F。
[0041]步骤二：读取带噪声的图像Y。
[0042]所述Y代表带噪声的图像，为一512
×
640矩阵。
[0043]步骤三：初始化尺度参数α＝0。
[0044]步骤四：定义辅助变量Z。
[0045]所述辅助变量Z的定义为
[0046]Z＝Y
‑
α
·
F,本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于最小熵的固定模式噪声去除方法，其特征在于：包括如下步骤，步骤一：通过暗场标定得到噪声模式F；步骤二：读取带噪声的图像Y，图像Y为一m
×
n矩阵；步骤三：初始化尺度参数α＝0；步骤四：定义辅助变量Z；步骤五：更新α；步骤六：重复步骤四
‑
五，直到α收敛；步骤七：减去固定模式噪声得到去噪后的图像X，达到去除图像中固定模式噪声的目的。2.如权利要求1所述的一种基于最小熵的固定模式噪声去除方法，其特征在于：步骤一的实现方法为，所述F为一m
×
n矩阵，其中m为图像高度，n为图像宽度；具体的标定方式为：在暗室中，将相机镜头盖盖上，然后采集图像，采集到的图像即为噪声模式F；为了抑制随机噪声，采集多张图像，然后平均，最终得到噪声模式F。3.如权利要求1所述的一种基于最小熵的固定模式噪声去除方法，其特征在于：步骤四中所述辅助变量Z的定义为Z＝Y
‑
α
·
F,Z也为一m<...

【专利技术属性】
技术研发人员：黄华，宋凌飞，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：