主权项 |
一种基于场景的自适应非均匀固定噪声去除方法,其特征是由以下步骤实现:步骤一,确定场景初始化;步骤二,场景自适应迭代;步骤二所描述的场景自适应迭代步骤中的迭代过程如下:假设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>Σ</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>∈</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>Σ</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>∈</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>Σ</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>∈</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>∂</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow><mrow><mo>∂</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>Σ</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>∈</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>∂</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup></mrow><mrow><mo>∂</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>Σ</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>∈</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>步骤三,场景模板更新;步骤四,生成输出图像。 |