发明名称 一种基于遥感图像相位相关性的图像质量判别方法
摘要 本发明公开了一种基于遥感图像相位相关性的图像质量判别方法,该方法包括以下步骤:先对遥感图像进行预处理得到遥感图像的周期分量频谱、然后向遥感图像的周期分量频谱中加入随机相位形成N幅加入随机相位形成多幅遥感图像,然后再进行离散傅里逆叶变换得到N幅时域遥感图像、再计算出N幅时域遥感图像peru1,peru2,…,peruN的全变差值和相位相关值,最后利用图像相位相关值对多幅遥感图像进行质量评判。该方法是一种无参考的数字图像质量评价方法,不需要任何参考图像数据,更适用于卫星遥感图像,而与其它无参考质量评价方法相比,本发明对各种降质因素均具有普适性,不会受到振铃、噪声等降质因素的影响。
申请公布号 CN102013100A 申请公布日期 2011.04.13
申请号 CN201010560839.2 申请日期 2010.11.26
申请人 北京空间机电研究所 发明人 满益云;阮宁娟;赵福立;张智;王殿中;鲍云飞;赵海博;许春晓
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 臧春喜
主权项 1.一种基于遥感图像相位相关性的图像质量判别方法,其特征在于包括以下步骤:(1)将一幅空间域的遥感图像u进行周期与平滑分解,分解成包含绝大多数图像信息的周期分量和在空间域中变化缓慢的平滑分量,然后将平滑分量滤除,对滤除平滑分量后的周期分量进行离散傅里叶变换得到周期分量的频谱,周期分量频谱的获取方式为:<maths num="0001"><![CDATA[<math><mrow><mo>&ForAll;</mo><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi><mo>/</mo><mrow><mo>(</mo><mn>0,0</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mover><mrow><mi>per</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow><mo>^</mo></mover><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>u</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mover><mi>v</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mrow><mrow><mn>4</mn><mo>-</mo><mn>2</mn><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;q</mi></mrow><mi>M</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;r</mi></mrow><mi>N</mi></mfrac><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mover><mrow><mi>peru</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mrow><mo>^</mo></mover><mo>=</mo><mo>|</mo><mover><mrow><mi>peru</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mrow><mo>^</mo></mover><mo>|</mo><msup><mi>e</mi><mrow><mi>i&phi;</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mrow></msup></mrow></math>]]></maths>其中<img file="FSA00000361858600014.GIF" wi="656" he="129" />且v<sub>1</sub>,v<sub>2</sub>满足:<maths num="0004"><![CDATA[<math><mrow><mo>&ForAll;</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi></mrow></math>]]></maths><img file="FSA00000361858600016.GIF" wi="1196" he="139" /><maths num="0005"><![CDATA[<math><mrow><mo>&ForAll;</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi></mrow></math>]]></maths><img file="FSA00000361858600018.GIF" wi="1235" he="153" />其中:q、r为遥感图像频域坐标;Ω/(0,0)为去除坐标原点的实数域;<img file="FSA00000361858600019.GIF" wi="202" he="110" />为周期分量的频谱;φ(q,r)为周期分量的相位;<img file="FSA000003618586000110.GIF" wi="128" he="50" />为遥感图像的频谱;v<sub>1</sub>、v<sub>2</sub>为图像像素的差值;M为遥感图像的长度,单位为像素;N为遥感图像的宽度,单位为像素;x、y为遥感图像空间域坐标;(2)向周期分量频谱的相位中加入N个随机相位形成N幅加入随机相位的遥感图像;所述向周期分量频谱的相位中加入N个随机相位的具体方式为:将相位φ增加一个随机偏移量εS,得到一个新的相位函数ψ=φ+εS,其中ε是一个固定值,S是在(-π,π)上满足独立同分布条件的均匀分布的随机变量,其中N为自然数;(3)将步骤(2)中生成的N幅加入随机相位的图像进行离散傅里叶逆变换,得到N幅时域遥感图像peru<sub>1</sub>,peru<sub>2</sub>,…,peru<sub>N</sub>;(4)计算N幅时域遥感图像peru<sub>1</sub>,peru<sub>2</sub>,…,peru<sub>N</sub>的全变差值,其中每幅时域遥感图像全变差值的计算方法为:<img file="FSA00000361858600021.GIF" wi="574" he="101" />(5)计算步骤(4)中N幅时域遥感图像全变差值TV(u<sub>N</sub>)的均值μ和方差σ,利用N幅时域遥感图像全变差值TV(u<sub>N</sub>)的均值μ和方差σ计算N幅时域遥感图像peru<sub>1</sub>,peru<sub>2</sub>,…,peru<sub>N</sub>的相位相关值IPC(u)即得到该幅空间域遥感图像u的相位相关值;<maths num="0006"><![CDATA[<math><mrow><mi>IPC</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>log</mi><mn>10</mn></msub><mi>&Phi;</mi><mrow><mo>(</mo><mfrac><mrow><mi>&mu;</mi><mo>-</mo><mi>TV</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>N</mi></msub><mo>)</mo></mrow></mrow><mi>&sigma;</mi></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中Φ为正态分布函数,即<img file="FSA00000361858600023.GIF" wi="563" he="121" />(6)利用步骤(1)至步骤(5)的计算方法计算出同一场景不同幅遥感图像的相位相关值IPC(u),对同一场景不同幅遥感图像的相位相关值IPC(u)进行比较,利用遥感图像相位相关值IPC(u)值的大小对遥感图像质量进行判别,判别方法为遥感图像相位相关值IPC(u)值越大遥感图像质量越好。
地址 100076 北京市丰台区南大红门路1号9201信箱5分箱