发明名称 一种基于场景的自适应非均匀固定噪声去除方法
摘要 本发明描述了一种基于场景的自适应非均匀固定噪声去除方法,其特征是由以下步骤实现:第一,确定场景初始化;第二,场景自适应迭代;第三,场景模板更新;第四,生成输出图像。本发明适用于非均匀固定噪声去除。该方法能适应场景的变化,不产生鬼影,收敛速度快,解决了现有算法收敛问题,鬼影问题。
申请公布号 CN104166960B 申请公布日期 2017.02.15
申请号 CN201410356021.7 申请日期 2014.07.24
申请人 中国科学院光电技术研究所 发明人 王辉;周进;雷涛
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 孟卜娟;卢纪
主权项 一种基于场景的自适应非均匀固定噪声去除方法,其特征是由以下步骤实现:步骤一,确定场景初始化;步骤二,场景自适应迭代;步骤二所描述的场景自适应迭代步骤中的迭代过程如下:假设x是存在非均匀性噪声的图像,经过算法处理后输出为y,对于像元(i,j),它们的关系为:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>*</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>b</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow>]]></math><img file="FDA0001127107840000011.GIF" wi="398" he="62" /></maths>其中,n表示第n帧图像,<img file="FDA0001127107840000012.GIF" wi="62" he="70" />和<img file="FDA0001127107840000013.GIF" wi="59" he="70" />分别为像元(i,j)的增益和偏移校正系数,<img file="FDA0001127107840000014.GIF" wi="61" he="63" />为是第n帧存在非均匀性噪声的图像中位置(i,j)的像素值,像元(i,j)期望输出<img file="FDA0001127107840000015.GIF" wi="67" he="63" />采用(i,j)相似邻域的空间均值,相似邻域用符号T表示,含义是距离相近,内容相似的区域,这样输出值由网络中的隐含层计算得到:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mi>t</mi></mfrac><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>T</mi></mrow></msub><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow>]]></math><img file="FDA0001127107840000016.GIF" wi="364" he="118" /></maths>其中,t表示相似邻域中像素总数;误差函数如下:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>e</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>M</mi></mrow></msub><msup><mrow><mo>(</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>*</mo><msub><mi>M</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>M</mi></mrow></msub><msup><mrow><mo>(</mo><msubsup><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><msubsup><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>b</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>*</mo><msub><mi>M</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow>]]></math><img file="FDA0001127107840000017.GIF" wi="1358" he="87" /></maths>其中,M表示图像更新模板,也是在迭代过程中计算的邻域;对其求关于<img file="FDA0001127107840000018.GIF" wi="62" he="63" />和<img file="FDA0001127107840000019.GIF" wi="66" he="62" />的偏导,结果如下:<maths num="0004"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow><mrow><mo>&part;</mo><msubsup><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><mo>=</mo><mn>2</mn><mo>*</mo><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>M</mi></mrow></msub><msubsup><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>*</mo><mrow><mo>(</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011271078400000110.GIF" wi="686" he="142" /></maths><maths num="0005"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow><mrow><mo>&part;</mo><msubsup><mi>b</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><mo>=</mo><mn>2</mn><mo>*</mo><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>M</mi></mrow></msub><mrow><mo>(</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011271078400000111.GIF" wi="598" he="143" /></maths>步骤三,场景模板更新;步骤四,生成输出图像。
地址 610209 四川省成都市双流350信箱