发明名称 一种基于欧氏距离正则化的SAR图像超分辨率重建方法
摘要 本发明公开了一种基于欧氏距离正则化的SAR图像超分辨率重建方法。首先输入时相相近的多帧SAR图像并配准,将其置于一个网格中并处理为像素均匀分布的SAR观测图像;然后在建立SAR图像退化模型的基础上,构建基于欧氏距离双边全差分正则化的高分辨率SAR图像重建公式;接着初始化观测图像,而后根据重建公式,采用一种可计算的优化方式对高分辨率SAR图像进行迭代估计;最后将估计图像作为SAR图像超分辨率重建结果输出。本发明将SAR图像重建中曼哈顿距离的双边全差分改为欧氏距离形式,使先验知识刻画更合理;为解决优化迭代过程中微分求解的困境,构建了重建公式的替代式,从而在优化求解可计算的基础上,提高了SAR图像超分辨率重建结果的准确性。
申请公布号 CN106169183A 申请公布日期 2016.11.30
申请号 CN201610551986.0 申请日期 2016.07.13
申请人 河海大学 发明人 徐枫;高红民;蒋德富;石爱业;张振;高建强
分类号 G06T5/50(2006.01)I;G06T7/00(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱桢荣
主权项 一种基于欧氏距离正则化的SAR图像超分辨率重建方法,其特征在于,包括以下步骤:1)、输入同一场景的多帧SAR图像;2)、对输入的多帧SAR图像进行空间配准;3)、根据配准结果,将所有SAR图像置于一个图像网格中,形成图像<img file="FDA0001047944450000011.GIF" wi="56" he="46" />其像素为非均匀分布;利用核回归对图像<img file="FDA0001047944450000012.GIF" wi="35" he="47" />进行处理,得到像素均匀分布的SAR观测图像y;4)、建立SAR图像的退化模型:y=Hx+n,其中,x为待估计的高分辨率SAR图像,H为已知的退化线性算子,n为独立同分布的高斯白噪声;5)、构建高分辨率SAR图像x的重建公式:<maths num="0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mo>=</mo><munder><mi>argmin</mi><mi>x</mi></munder><mrow><mo>(</mo><mrow><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><mi>H</mi><mi>x</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mi>&alpha;</mi><mrow><mi>e</mi><mrow><mo>(</mo><mi>l</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow></msup><munder><mo>&Sigma;</mo><mi>i</mi></munder><mo>|</mo><msub><mrow><mo>(</mo><msub><mo>&dtri;</mo><mrow><mi>l</mi><mi>m</mi></mrow></msub><mi>x</mi><mo>)</mo></mrow><mi>i</mi></msub><mo>|</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001047944450000013.GIF" wi="862" he="126" /></maths>其中,λ是设定的正则参数,α是设定的距离度量参数,p为设定的图像最大差分阶数,l和m分别表示像素间水平和垂直距离,e(l,m)为像素间距,<img file="FDA0001047944450000014.GIF" wi="88" he="47" />为x的lm阶的梯度,<img file="FDA0001047944450000015.GIF" wi="149" he="71" />为<img file="FDA0001047944450000016.GIF" wi="90" he="53" />在第i个像素处的模;正则项中的e(l,m)和<img file="FDA0001047944450000017.GIF" wi="146" he="78" />均表示为欧氏距离形式:<maths num="0002"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>l</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msup><mi>l</mi><mn>2</mn></msup><mo>+</mo><msup><mi>m</mi><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001047944450000018.GIF" wi="334" he="63" /></maths><maths num="0003"><math><![CDATA[<mrow><mo>|</mo><msub><mrow><mo>(</mo><msub><mo>&dtri;</mo><mrow><mi>l</mi><mi>m</mi></mrow></msub><mi>x</mi><mo>)</mo></mrow><mi>i</mi></msub><mo>|</mo><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mrow><mo>(</mo><mrow><msup><mi>&Delta;</mi><mrow><mi>l</mi><mi>h</mi></mrow></msup><mi>x</mi></mrow><mo>)</mo></mrow><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mrow><mo>(</mo><mrow><msup><mi>&Delta;</mi><mrow><mi>m</mi><mi>v</mi></mrow></msup><mi>x</mi></mrow><mo>)</mo></mrow><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001047944450000019.GIF" wi="638" he="103" /></maths>其中,Δ<sup>lh</sup>x和Δ<sup>mv</sup>x分别表示x在水平和垂直方向上的l阶和m阶差分,(Δ<sup>lh</sup>x)<sub>i</sub>和(Δ<sup>mv</sup>x)<sub>i</sub>分别表示Δ<sup>lh</sup>x和Δ<sup>mv</sup>x在第i个像素处的值;6)、通过维纳滤波,将SAR观测图像y初始化为SAR图像x<sup>(0)</sup>;7)、在已知初始SAR图像x<sup>(0)</sup>的条件下,根据步骤5)中的重建公式,对高分辨率SAR图像x进行迭代估计;8)、将步骤7)迭代估计的结果作为SAR图像x的超分辨率重建结果,并输出。
地址 211100 江苏省南京市江宁开发区佛城西路8号