CN104680541B - 基于相位一致性的遥感图像质量评价方法 - Google Patents
基于相位一致性的遥感图像质量评价方法 Download PDFInfo
- Publication number
- CN104680541B CN104680541B CN201510112767.8A CN201510112767A CN104680541B CN 104680541 B CN104680541 B CN 104680541B CN 201510112767 A CN201510112767 A CN 201510112767A CN 104680541 B CN104680541 B CN 104680541B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- value
- mrow
- sensing images
- represent
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明公开了一种基于相位一致性的遥感图像质量评价方法，克服了现有技术仅以单个参数表征遥感图像质量优劣的不足，实现步骤为：(1)输入一幅待评价遥感图像；(2)获得遥感图像在条带噪声方向的相位一致性值；(3)对相位一致性值进行二值化处理；(4)获取条带噪声强度矩阵；(5)计算条带噪声值；(6)获得一般噪声信息；(7)计算遥感图像的一般噪声值；(8)计算遥感图像的模糊度值；(9)利用choquet模糊积分，获得遥感图像的综合评价值；(10)输出遥感图像质量的综合评价值。本发明具有无参考性、全面评价遥感图像质量的优点，能够使评价结果能更好地符合主观评价。
Description
技术领域
本发明属于图像处理技术领域，更进一步涉及遥感图像处理技术领域中的一种基于相位一致性的遥感图像质量评价方法。本发明可应用于对遥感图像的质量进行客观评价，用于保证获取遥感图像数据的有效性和准确性。
背景技术
遥感图像在成像、传输和存储的过程中，受到光学系统、大气环境以及其它因素的影响，最终会使得地面获得的图像存在条带噪声、一般噪声、模糊等各种图像质量下降的问题，遥感图像质量是比较各种遥感图像处理算法性能优劣以及优化遥感光学系统参数的重要指标，因此在遥感图像处理领域建立有效的图像质量评价机制具有重要的意义。
目前，图像质量评价一般包括主观评价和客观评价。主观评价是以人眼视觉效果作为质量评价准则，需要多名专业观测者参加评分，没有统一标准，需要多人参与，繁琐、费时费力、不能实时处理，不利于海量数据的自动判读。图像质量的客观评价方法根据需要参考图像信息量的多少可以分为全参考评价、部分参考评价和无参考评价三类。常用客观评价方法主要有峰值信噪比(Peak Signal-to-Noise Ratio,PSNR)、结构相似度(StructureSimilarity,SSIM)、信息熵、平均梯度等方法。
张飞艳、谢伟等人在论文“No-reference remote sensing image qualityassessment based on natural scene statistical in wavelet domain”(Journal ofElectronics and Information Technology,2011,33(11):2742-2747)中用遥感图像小波域统计特性，通过噪声、模糊对影像线性关系的影响方式及破坏程度的定量分析，来判定影像的噪声和模糊强度，最后，利用二者的加权和作为遥感影像的最终质量评价指标。该方法是一种无参考的图像质量评价方法，该方法仍然存在的不足是，综合评价模型比较简单，没有考虑人眼视觉的特点，评价结果不能更好的符合主观评价的结果。
北京空间机电研究所申请的专利“一种基于遥感图像相位相关性的图像质量判别方法”(公开号：CN102013100A，公开日：2011年4月13日，申请日：2010年11月26日)中公开了一种基于遥感图像相位相关性的图像质量判别方法。该方法先对遥感图像进行预处理得到遥感图像的周期分量频谱，然后向遥感图像的周期分量频谱中加入随机相位形成的N幅加入随机相位形成多幅遥感图像,再进行离散傅里叶逆变换得到N幅时域遥感图像，再计算出N幅时域遥感图像peru1,peru2,…,peruN的全变差值和相位相关值，最后利用图像相位相关值对多幅遥感图像进行质量评判。该方法是一种无参考的数字图像质量评价方法，优点是不需要任何参考图像数据，更适合于评价遥感图像。但是，该方法仍然存在的不足是，但是，该方法仍然存在的不足是，只考虑了图像的锐度，仅以单个参数表征图像质量优劣，没有综合考虑遥感图像中各个降质参数之间的相互关系，评价结果不能更好的符合主观评价的结果。
发明内容
本发明的目的是针对上述现有技术的不足，提出一种基于相位一致性的遥感图像质量评价方法。本发明综合考虑遥感图像中条带噪声、一般噪声和模糊等因素的相互关系，利用相位一致性能够检测图像中边缘、阶跃和马赫带等特征点的特性，通过相位一致性的计算方式检测出遥感图像中的条带噪声、一般噪声和模糊，利用相位一致性表征对应的强度指标并得到评价参数，最后利用Choquet非线性函数将三个评价参数结合，得到一个能够对遥感图像进行全面评价的综合评价值，评价结果能更好地符合人的主观评价。
本发明实现的具体步骤包括如下：
(1)输入一幅待评价的遥感图像；
(2)获得相位一致性值：
采用相位一致性公式，计算遥感图像在条带噪声方向的相位一致性，得到相
位一致性值；
(3)二值化处理：
采用二值化方法，对相位一致性值进行二值化处理，得到0-1模板；
(4)获得条带噪声强度矩阵：
用遥感图像矩阵乘以0-1模板，得到条带噪声，将条带噪声的像素值进行归一化处理，得到条带噪声强度矩阵；
(5)按照下式，计算遥感图像的条带噪声值：
其中，N1表示遥感图像的条带噪声值；∑表示求和操作；j表示遥感图像像素的列坐标；i表示遥感图像像素的行坐标；S表示条带噪声强度矩阵；P表示遥感图像的相位一致性值；M表示遥感图像的行数；N表示遥感图像的列数；
(6)获得一般噪声信息：
(6a)采用边缘提取方法，提取遥感图像的边缘，得到边缘图；
(6b)采用形态学中的膨胀运算，对边缘图进行膨胀，得到膨胀后的边缘图；
(6c)采用局部相位算法，计算遥感图像的局部相位值，得到局部相位图；
(6d)从膨胀后的边缘图中找出像素值为1的坐标，将找出的坐标所对应的局部相位图中像素值赋值为0，得到去掉边缘信息的局部相位图；
(7)按照下式，计算遥感图像的一般噪声值：
其中，N2表示遥感图像一般噪声的值；∑表示求和操作；p表示去掉边缘信息的局部相位图的行标号；M表示遥感图像的行数；l表示去掉边缘信息的局部相位图的列标号；N表示遥感图像的列数；P(p,l)表示去掉边缘信息的局部相位图像素点(p,l)的像素值；(p,l)表示去掉边缘信息的局部相位图像素坐标；
(8)获得遥感图像的模糊度值：
(8a)采用局部相位算法，计算遥感图像的局部相位值，得到局部相位图；
(8b)利用matlab软件中的reshape函数，将局部相位图矩阵拉伸成为一维数组，得到局部相位图的一维数组；
(8c)将局部相位图的一维数组按照从小到大进行排序；
(8d)采用权重函数，计算权重值；
(8e)按照下式，计算遥感图像的模糊度值:
其中，B表示遥感图像的模糊度值；∑表示求和操作；n表示局部相位图的一维数组从1到M×N共M×N个数的标号；M表示遥感图像的行数；N表示遥感图像的列数；Wn表示局部相位图的一维数组中第n个元素的权重值；Pn表示局部相位图的一维数组第n个元素的像素值；n表示数组从1到M×N共M×N个数的标号；
(9)计算遥感图像质量的综合评价值：
(9a)采用极差法，分别对遥感图像的条带噪声值、一般噪声值和模糊度值进行归一化处理；
(9b)将归一化处理后的遥感图像的条带噪声值、一般噪声值和模糊度值按从大到小进行排序；
(9c)根据遥感图像中条带噪声、一般噪声和模糊度在遥感图像质量中的重要程度，确定综合评价因子中各个值得的权重系数；
(9d)按照模糊测度公式，计算按降序排列后条带噪声值、一般噪声值和模糊度值的模糊测度；
(9e)按照下式，计算遥感图像质量的综合评价值：
Q＝I1g1+I2g2+I3g3
其中，Q表示遥感图像质量的综合评价值，Q越大，表示遥感图像的质量越差；I1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值；g1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值的模糊测度；I2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值；g2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值的模糊测度；I3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值；g3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值的模糊测度；
(10)输出遥感图像质量的综合评价值。
与现有技术相比，本发明具有以下优点：
第一，由于本发明在获得遥感图像的条带噪声值时，应用遥感图像的相位信息，从含有条带噪声的遥感图像中准确定位条带噪声的位置，克服了现有遥感图像评价技术中没有考虑条带噪声的问题，使得本发明具有在评价遥感图像条带噪声时结果更加精确的优点。
第二，由于本发明在遥感图像综合指标中，将导致遥感图像质量下降的条带噪声、一般噪声和模糊度综合起来，克服了现有技术单纯由人眼判决产生的主观性和仅以单个参数表征遥感图像质量优劣的片面性的问题，使得本发明在评价遥感图像质量时的结果更加全面优点。
第三，由于本发明在遥感图像综合评价指标的建立中采用Choquet模糊积分函数，充分考虑了人的视觉感知，克服了现有技术中遥感图像质量评价综合模型简单，权重系数没有考虑人眼视觉特点的问题，使得本发明在综合评价遥感图像质量时符合人的主观评价的优点。
附图说明
图1是本发明评价遥感图像质量的流程图。
具体实施方式
参照图1，对本发明具体实施方式做进一步的描述。.
步骤1，输入一幅待评价的遥感图像。
在计算机中应用matlab软件分别读取存储在计算机硬盘空间中待评价遥感图像。
步骤2，获得相位一致性值。
采用下述的相位一致性公式，计算遥感图像在条带噪声方向的相位一致性，得到相位一致性值。
其中，P(x,y)表示在遥感图像像素点(x,y)的相位一致性；P(x,y)越接近于1，表示遥感图像在像素点(x,y)的相位一致性越好；∑表示求和操作；o表示Log-Gabor小波滤波器的方向；Eo(x,y)表示对遥感图像进行在Log-Gabor小波滤波方向o上的滤波，得到的偶对称频率响应和奇对称频率响应的平方和，再对其平方和开方；Ado(x,y)表示在Log-Gabor小波滤波器的尺度d和Log-Gabor小波滤波器的方向o上遥感图像像素点(x,y)与Log-Gabor小波滤波器卷积的幅值；ε表示0.01的常数。
步骤3，二值化处理。
采用下述的二值化方法，对相位一致性值进行二值化处理，得到0-1模板。
第一步，用相位一致性值的最大值，减去相位一致性值的最小值，用其差值乘以0.146，将乘积结果作为二值化操作的阈值。
第二步，将相位一致性值中大于二值化阈值的值赋值为1，将相位一致性值中小于二值化阈值的值赋值为0，完成相位一致性值的二值化处理。
步骤4，获得条带噪声强度矩阵。
用遥感图像矩阵乘以0-1模板，得到条带噪声，将条带噪声中的每个像素值除以条带噪声的最大像素值进行归一化处理，得到条带噪声强度矩阵。
步骤5，采用下述的条带噪声公式，计算遥感图像的条带噪声值。
其中，N1表示遥感图像的条带噪声值，N1越大，遥感图像受条带噪声污染越严重，∑表示求和操作，j表示遥感图像像素的列坐标，i表示遥感图像像素的行坐标，S表示条带噪声强度矩阵，P表示遥感图像的相位一致性值，M表示遥感图像的行数，N表示遥感图像的列数。
步骤6，获得一般噪声信息。
采用下述的边缘提取方法，提取遥感图像的边缘，得到边缘图。
第一步，按照下式，计算噪声能量：
To＝2μ1+2σ1
其中，To表示在Log-Gabor小波滤波器o方向上遥感图像的噪声能量，o表示Log-Gabor小波滤波器方向，μ1表示通过Log-Gabor小波滤波器滤波后遥感图像像素的平均值，σ1表示通过Log-Gabor小波滤波器后遥感图像像素的标准差。
第二步，按照下式，计算遥感图像的边缘值：
其中，E表示遥感图像的边缘值，∑表示求和操作，o表示Log-Gabor小波滤波器方向，Eo表示在Log-Gabor小波滤波器o方向上遥感图像和Log-Gabor小波滤波器卷积的幅值，To表示在Log-Gabor小波滤波器o方向上的噪声能量，表示向下取整操作，cos(·)表示余弦操作，Φo表示在Log-Gabor小波滤波器o方向上遥感图像和Log-Gabor小波滤波器卷积后的相位值。
采用形态学中的膨胀运算，对边缘图进行膨胀，得到膨胀后的边缘图。
采用下述的局部相位算法，计算遥感图像的局部相位值，得到局部相位图。
局部相位算法的步骤如下:
第一步，按照下式，计算噪声能量：
Tk＝μ2+2σ2
其中，Tk表示在Log-Gabor小波滤波器k方向上遥感图像的噪声能量，k表示Log-Gabor小波滤波器的方向，μ2表示通过Log-Gabor小波滤波器后遥感图像像素点的平均值，σ2表示通过Log-Gabor小波滤波器后遥感图像像素点的标准差。
第二步，按照下式，计算遥感图像的局部相位值：
其中，M2表示遥感图像的局部相位值，k表示Log-Gabor小波滤波器的方向，∑表示求和操作，Ek表示在Log-Gabor小波滤波器k方向上遥感图像和Log-Gabor小波滤波器卷积后的幅值，Tk表示在Log-Gabor小波滤波器k方向上，遥感图像的噪声能量，表示向下取整操作，cos(·)表示求余弦值操作，Φk表示在Log-Gabor小波滤波器k方向上遥感图像和Log-Gabor小波滤波器卷积后的相位值。
从膨胀后的边缘图中找出像素值为1的坐标，将找出的坐标所对应的局部相位图中像素值赋值为0，得到去掉边缘信息的局部相位图。
步骤7，采用下述的一般噪声公式，计算遥感图像的一般噪声值。
其中，N2表示遥感图像一般噪声的值，N2越大，表示遥感图像受一般噪声污染越严重，∑表示求和操作，p表示去掉边缘信息的局部相位图的行标号，M表示遥感图像的行数，l表示去掉边缘信息的局部相位图的列标号，N表示遥感图像的列数，P(p,l)表示去掉边缘信息的局部相位图像素点(p,l)的像素值；(i,j)表示去掉边缘信息的局部相位图像素坐标。
步骤8，获得遥感图像的模糊度值。
采用下述的局部相位算法，计算遥感图像的局部相位值，得到局部相位图。
局部相位算法具体步骤如下：
第一步，按照下式，计算噪声能量：
Tk＝μ2+2σ2
其中，Tk表示在Log-Gabor小波滤波器k方向上遥感图像的噪声能量，k表示Log-Gabor小波滤波器的方向，μ2表示通过Log-Gabor小波滤波器后遥感图像像素点的平均值，σ2表示通过Log-Gabor小波滤波器后遥感图像像素点的标准差。
利用matlab软件中的reshape函数将局部相位图矩阵拉伸成为一维数组，得到局部相位图的一维数组。
将局部相位图的一维数组按照从小到大进行排序。
采用下述的权重函数，计算权重值。
其中，Wq表示按照升序排列后一维数组中第q个元素的权重值，exp(·)表示以自然常数e为底的指数操作，q表示数组从1到M×N共M×N个数的标号，M表示遥感图像的行数，N表示遥感图像的列数。
按照下述的模糊度公式，计算遥感图像的模糊度值：
其中，B表示遥感图像的模糊度值，∑表示求和操作，n表示局部相位图的一维数组从1到M×N共M×N个数的标号，M表示遥感图像的行数，N表示遥感图像的列数，Wn表示局部相位图的一维数组中第n个元素的权重值，Pn表示局部相位图的一维数组第n个元素的像素值；n表示数组从1到M×N共M×N个数的标号；
步骤9，计算遥感图像质量的综合评价值。
采用极差法，分别对遥感图像的条带噪声值、一般噪声值和模糊度值进行归一化处理。
将归一化处理后的遥感图像的条带噪声值、一般噪声值和模糊度值按从大到小进行排序。
根据遥感图像中条带噪声、一般噪声和模糊度在遥感图像质量中的重要程度，确定综合评价因子中各个值得的权重系数。
根据choquet模糊积分中，权重系数的确定方法，专家按照遥感图像中条带噪声、一般噪声和模糊度在遥感图像质量中的重要程度，得出语气算子，最后获得权重系数。
按照模糊测度公式，计算按降序排列后条带噪声值、一般噪声值和模糊度值的模糊测度。
模糊测度公式如下：
其中，m表示遥感图像中的条带噪声值、一般噪声值和模糊度的标号；g1(Im)表示遥感图像中的条带噪声值、一般噪声值和模糊度值按降序排列后第m位的模糊测度值，∑表示求和操作，g(Im)表示降序排列后第m位的权重系数。
利用choquet模糊积分模型，计算遥感图像质量的综合评价值。
Q＝I1g1+I2g2+I3g3
其中，Q表示遥感图像质量的综合评价值，Q越大，表示遥感图像的质量越差；I1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值，g1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值的模糊测度，I2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值，g2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值的模糊测度，I3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值，g3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值的模糊测度，
步骤10，输出遥感图像质量的综合评价值。
Claims (7)
1.一种基于相位一致性的遥感图像质量评价方法，包括如下步骤：
(1)输入一幅待评价的遥感图像；
(2)获得相位一致性值：
采用相位一致性公式，计算遥感图像在条带噪声方向的相位一致性，得到相位一致性值；
(3)二值化处理：
采用二值化方法，对相位一致性值进行二值化处理，得到0-1模板；
(4)获得条带噪声强度矩阵：
用遥感图像矩阵乘以0-1模板，得到条带噪声，将条带噪声的像素值进行归一化处理，得到条带噪声强度矩阵；
(5)按照下式，计算遥感图像的条带噪声值：
<mrow>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<mi>S</mi>
<mo>&times;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</mfrac>
</mrow>
其中，N1表示遥感图像的条带噪声值；∑表示求和操作；j表示遥感图像像素的列坐标；i表示遥感图像像素的行坐标；S表示条带噪声强度矩阵；P表示遥感图像的相位一致性值；M表示遥感图像的行数；N表示遥感图像的列数；
(6)获得一般噪声信息：
(6a)采用边缘提取方法，提取遥感图像的边缘，得到边缘图；
(6b)采用形态学中的膨胀运算，对边缘图进行膨胀，得到膨胀后的边缘图；
(6c)采用局部相位算法，计算遥感图像的局部相位值，得到局部相位图；
(6d)从膨胀后的边缘图中找出像素值为1的坐标，将找出的坐标所对应的局部相位图中像素值赋值为0，得到去掉边缘信息的局部相位图；
(7)按照下式，计算遥感图像的一般噪声值：
<mrow>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</mfrac>
</mrow>
其中，N2表示遥感图像一般噪声的值；∑表示求和操作；p表示去掉边缘信息的局部相位图的行标号；M表示遥感图像的行数；l表示去掉边缘信息的局部相位图的列标号；N表示遥感图像的列数；P(p,l)表示去掉边缘信息的局部相位图像素点(p,l)的像素值；(p,l)表示去掉边缘信息的局部相位图像素坐标；
(8)获得遥感图像的模糊度值：
(8a)采用局部相位算法，计算遥感图像的局部相位值，得到局部相位图；
(8b)利用matlab软件中的reshape函数，将局部相位图矩阵拉伸成为一维数组，得到局部相位图的一维数组；
(8c)将局部相位图的一维数组按照从小到大进行排序；
(8d)采用权重函数，计算权重值；
(8e)按照下式，计算遥感图像的模糊度值:
<mrow>
<mi>B</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>W</mi>
<mi>n</mi>
</msub>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>W</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中，B表示遥感图像的模糊度值；∑表示求和操作；n表示局部相位图的一维数组从1到M×N共M×N个数的标号；M表示遥感图像的行数；N表示遥感图像的列数；Wn表示局部相位图的一维数组中第n个元素的权重值；Pn表示局部相位图的一维数组第n个元素的像素值；n表示数组从1到M×N共M×N个数的标号；
(9)计算遥感图像质量的综合评价值：
(9a)采用极差法，分别对遥感图像的条带噪声值、一般噪声值和模糊度值进行归一化处理；
(9b)将归一化处理后的遥感图像的条带噪声值、一般噪声值和模糊度值按从大到小进行排序；
(9c)根据遥感图像中条带噪声、一般噪声和模糊度在遥感图像质量中的重要程度，确定综合评价因子中各个值得的权重系数；
(9d)按照模糊测度公式，计算按降序排列后条带噪声值、一般噪声值和模糊度值的模糊测度；
(9e)按照下式，计算遥感图像质量的综合评价值：
Q＝I1g1+I2g2+I3g3
其中，Q表示遥感图像质量的综合评价值，Q越大，表示遥感图像的质量越差；I1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值；g1表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最大值的模糊测度；I2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值；g2表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的中间值的模糊测度；I3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值；g3表示归一化处理后的条带噪声值、一般噪声值和模糊度值中的最小值的模糊测度；
(10)输出遥感图像质量的综合评价值。
2.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(2)中所述的相位一致性公式如下：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>o</mi>
</munder>
<msub>
<mi>E</mi>
<mi>o</mi>
</msub>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>o</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>d</mi>
</munder>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>&epsiv;</mi>
</mrow>
</mfrac>
</mrow>
其中，P表示遥感图像的相位一致性值；P越接近于1，表示遥感图像的相位一致性越好；∑表示求和操作；o表示Log-Gabor小波滤波器的方向；Eo表示在Log-Gabor小波滤波器方向o上遥感图像的能量值；d表示Log-Gabor小波滤波器的尺度；Ado表示在Log-Gabor小波滤波器的尺度d和Log-Gabor小波滤波器的方向o上遥感图像与Log-Gabor小波滤波器卷积的幅值；ε表示0.01的常数。
3.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(3)中所述二值化方法的具体步骤如下：
第一步，用相位一致性值的最大值减去相位一致性值的最小值，用其差值乘以0.146，将乘积结果作为二值化的阈值；
第二步，将相位一致性值中大于二值化阈值的值赋值为1，将相位一致性值中小于二值化阈值的值赋值为0，完成相位一致性值的二值化处理。
4.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(6a)中所述提取边缘方法的步骤如下：
第一步，按照下式，计算噪声能量：
To＝2μ1+2σ1
其中，To表示在Log-Gabor小波滤波器o方向上遥感图像的噪声能量；o表示Log-Gabor小波滤波器方向；μ1表示通过Log-Gabor小波滤波器滤波后遥感图像像素的平均值；σ1表示通过Log-Gabor小波滤波器后遥感图像像素的标准差；
第二步，按照下式，计算遥感图像的边缘值：
其中，E表示遥感图像的边缘值；∑表示求和操作；o表示Log-Gabor小波滤波器方向；Eo表示在Log-Gabor小波滤波器o方向上遥感图像和Log-Gabor小波滤波器卷积的幅值；To表示在Log-Gabor小波滤波器o方向上的噪声能量；表示向下取整操作，cos(·)表示余弦操作，Φo表示在Log-Gabor小波滤波器o方向上遥感图像和Log-Gabor小波滤波器卷积后的相位值。
5.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(6c)、步骤(8a)中所述局部相位算法的步骤如下：
第一步，按照下式，计算噪声能量：
Tk＝μ2+2σ2
其中，Tk表示在Log-Gabor小波滤波器k方向上遥感图像的噪声能量；k表示Log-Gabor小波滤波器的方向；μ2表示通过Log-Gabor小波滤波器后遥感图像像素点的平均值；σ2表示通过Log-Gabor小波滤波器后遥感图像像素点的标准差；
第二步，按照下式，计算遥感图像的局部相位值：
其中，M2表示遥感图像的局部相位值；k表示Log-Gabor小波滤波器的方向；∑表示求和操作；Ek表示在Log-Gabor小波滤波器k方向上遥感图像和Log-Gabor小波滤波器卷积后的幅值；Tk表示在Log-Gabor小波滤波器k方向上，遥感图像的噪声能量；表示向下取整操作，cos(·)表示求余弦值操作，Φk表示在Log-Gabor小波滤波器k方向上遥感图像和Log-Gabor小波滤波器卷积后的相位值。
6.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(8d)中所述权重函数如下：
<mrow>
<msub>
<mi>W</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>q</mi>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>0.0001</mn>
<mo>&rsqb;</mo>
</mrow>
其中，Wq表示按照升序排列后一维数组中第q个元素的权重值；exp(·)表示以自然常数e为底的指数操作；q表示数组从1到M×N共M×N个数的标号；M表示遥感图像的行数，N表示遥感图像的列数。
7.根据权利要求1所述的基于相位一致性的遥感图像质量评价方法，其特征在于，步骤(9d)中所述模糊测度公式如下：
<mrow>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>&Element;</mo>
<mo>{</mo>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>I</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>I</mi>
<mn>3</mn>
</msub>
<mo>}</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中，m表示遥感图像中的条带噪声值、一般噪声值和模糊度的标号；g1(Im)表示遥感图像中的条带噪声值、一般噪声值和模糊度值按降序排列后第m位的模糊测度值；∈表示属于符号；∑表示求和操作；g(Im)表示降序排列后第m位的权重系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510112767.8A CN104680541B (zh) | 2015-03-15 | 2015-03-15 | 基于相位一致性的遥感图像质量评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510112767.8A CN104680541B (zh) | 2015-03-15 | 2015-03-15 | 基于相位一致性的遥感图像质量评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104680541A CN104680541A (zh) | 2015-06-03 |
| CN104680541B true CN104680541B (zh) | 2018-03-13 |
Family
ID=53315534
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510112767.8A Active CN104680541B (zh) | 2015-03-15 | 2015-03-15 | 基于相位一致性的遥感图像质量评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104680541B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105389820A (zh) * | 2015-11-18 | 2016-03-09 | 成都中昊英孚科技有限公司 | 一种基于倒谱的红外图像清晰度评价方法 |
| CN105427349B (zh) * | 2015-11-18 | 2018-05-25 | 国家海洋局第一海洋研究所 | 遥感图像条带噪声定量评价方法 |
| CN105427323B (zh) * | 2015-12-07 | 2019-05-21 | 湖南大学 | 一种基于相位一致性的激光熔覆熔池边缘提取方法 |
| CN106023190B (zh) * | 2016-05-16 | 2019-05-07 | 浙江大学 | 一种自动对焦图像准焦程度的评价方法 |
| CN106408561A (zh) * | 2016-09-10 | 2017-02-15 | 天津大学 | 一种基于图像纹理特征的无参考图像质量评价方法 |
| CN106530301B (zh) * | 2016-11-30 | 2019-06-07 | 上海卫星工程研究所 | 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 |
| CN110070519B (zh) * | 2019-03-13 | 2023-07-14 | 西安电子科技大学 | 基于相位一致性的拼接图像质量测定方法、图像拼接系统 |
| CN111798421A (zh) * | 2020-06-29 | 2020-10-20 | 浙江同善人工智能技术有限公司 | 一种图像质量判断方法、设备及存储介质 |
| CN111899261A (zh) * | 2020-08-31 | 2020-11-06 | 西北工业大学 | 一种水下图像质量实时评估方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102075786A (zh) * | 2011-01-19 | 2011-05-25 | 宁波大学 | 一种图像客观质量评价方法 |
| CN102142145A (zh) * | 2011-03-22 | 2011-08-03 | 宁波大学 | 一种基于人眼视觉特性的图像质量客观评价方法 |
| CN102843572A (zh) * | 2012-06-29 | 2012-12-26 | 宁波大学 | 一种基于相位的立体图像质量客观评价方法 |
| EP2396768B1 (en) * | 2009-02-12 | 2013-04-17 | Dolby Laboratories Licensing Corporation | Quality evaluation of sequences of images |
| CN103200421A (zh) * | 2013-04-07 | 2013-07-10 | 北京理工大学 | 一种基于Curvelet变换和相位一致的无参考图像质量评价方法 |
-
2015
- 2015-03-15 CN CN201510112767.8A patent/CN104680541B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2396768B1 (en) * | 2009-02-12 | 2013-04-17 | Dolby Laboratories Licensing Corporation | Quality evaluation of sequences of images |
| CN102075786A (zh) * | 2011-01-19 | 2011-05-25 | 宁波大学 | 一种图像客观质量评价方法 |
| CN102142145A (zh) * | 2011-03-22 | 2011-08-03 | 宁波大学 | 一种基于人眼视觉特性的图像质量客观评价方法 |
| CN102843572A (zh) * | 2012-06-29 | 2012-12-26 | 宁波大学 | 一种基于相位的立体图像质量客观评价方法 |
| CN103200421A (zh) * | 2013-04-07 | 2013-07-10 | 北京理工大学 | 一种基于Curvelet变换和相位一致的无参考图像质量评价方法 |
Non-Patent Citations (3)
| Title |
|---|
| 《No-reference Remote Sensing Image Quality Assessment Using a Comprehensive Evaluation Factor》;Lin Wang等;《Proc. of SPIE》;20141231;第9124卷;全文 * |
| 《基于灰度共生矩阵的无参考模糊图像质量评价方法》;桑庆兵等;《模式识别与人工智能》;20130531;第26卷(第5期);全文 * |
| 《学习相位一致特征的无参考图像质量评价》;李朝锋等;《电子与信息学报》;20130228;第35卷(第2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104680541A (zh) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104680541B (zh) | 基于相位一致性的遥感图像质量评价方法 | |
| CN108052980B (zh) | 基于图像的空气质量等级检测方法 | |
| Zahran et al. | Automatic weld defect identification from radiographic images | |
| CN102170581B (zh) | 基于hvs的ssim与特征匹配立体图像质量评价方法 | |
| CN104156736B (zh) | 基于sae和idl的极化sar图像分类方法 | |
| WO2017024692A1 (zh) | 一种单测量节点模拟电路故障诊断方法 | |
| CN105208374A (zh) | 一种基于深度学习的无参考图像质量客观评价方法 | |
| CN106600595A (zh) | 一种基于人工智能算法的人体特征尺寸自动测量方法 | |
| CN105260998A (zh) | 基于mcmc采样和阈值低秩逼近的图像去噪方法 | |
| CN107909109A (zh) | 基于显著性和多尺度深度网络模型的sar图像分类方法 | |
| CN106780446A (zh) | 一种无参考混合失真图像质量评价方法 | |
| CN105374047B (zh) | 基于改进的双边滤波与聚类的sar图像变化检测方法 | |
| CN103310235B (zh) | 一种基于参数识别与估计的隐写分析方法 | |
| CN106570862A (zh) | 一种超分辨率重建质量的评价方法及装置 | |
| CN107343196A (zh) | 一种混合失真无参考图像质量评价方法 | |
| CN106651834A (zh) | 无参考图像的变电设备红外热像图质量评价方法及装置 | |
| CN109003265A (zh) | 一种基于贝叶斯压缩感知的无参考图像质量客观评价方法 | |
| CN107516082A (zh) | 基于自步学习的sar图像变化区域检测方法 | |
| Wang et al. | Distortion recognition for image quality assessment with convolutional neural network | |
| CN105120294A (zh) | 一种jpeg格式图像来源鉴别方法 | |
| CN107590785A (zh) | 一种基于sobel算子的布里渊散射谱图像识别方法 | |
| CN101853401B (zh) | 一种基于二维经验模态分解的多分组图像分类方法 | |
| CN104951800A (zh) | 一种面向资源开采型地区的遥感影像融合方法 | |
| CN108830829B (zh) | 联合多种边缘检测算子的无参考质量评价算法 | |
| Jobin et al. | Plant identification based on fractal refinement technique (FRT) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |