本发明涉及一种飞机尾迹云识别方法，尤其涉及一种基于卫星云图资料的飞机尾迹云识别方法。
背景技术：
云是大气辐射系统中重要的组成部分，云覆盖率的改变会影响整个大气辐射的收支平衡，进而对气候系统产生影响。尾迹云是飞机飞行期间排出的废气与周围环境的冷空气经过混合后，由于相对湿度的增加，水汽凝结形成的一种特殊的云系。通常来说，自然状态下卷云的形成需要相对湿度达到145%-165%甚至更高，而尾迹云在弱的冰相过饱和（相对湿度100%-110%）的空气中就可以存在。这种差异意味着，在对流层高层的很大一片区域，尾迹云能够生成并长时间持续存在，从而改变高层卷云的覆盖率。另外，在尾迹云的生成和持续过程中，通过降低空气中的水汽含量，可以影响自然卷云的生成和持续，从而改变自然卷云的覆盖率及光学厚度。
目前已有研究表明，尾迹云通过自身的直接作用以及改变自然卷云的间接作用，对地球的辐射强迫已经超过了飞机排放的二氧化碳的辐射强迫。近年来，随着我国及周边地区经济和贸易的不断交流与发展，东亚地区的民用和商用航班日益增多，这些航空活动对气候的影响也在逐渐增大。可迄今为止，有关东亚地区飞机尾迹云的分布特征及其辐射强迫的研究甚少。线性尾迹云的识别对于线性尾迹云的覆盖率和光学厚度的计算至关重要。
技术实现要素：
本发明为了解决现有技术中存在的问题，提供一种本发明是通过卫星图像进行飞机尾迹云识别的方法，通过对卫星高分辨率资料的分析，准确地识别线性尾迹云，并计算尾迹云的覆盖率从而用于研究其辐射特性，将会进一步完善尾迹云对区域和全球气候影响的研究，对研究尾迹云的辐射强迫有重要意义。
为了达到上述目的，本发明提出的技术方案为：一种基于卫星资料的飞机尾迹云识别方法，可以通过以下技术方案来实现：一种基于卫星资料的飞机尾迹云识别方法，包括如下步骤：
步骤1、融合卫星资料中的红外和可见光通道的云图；
步骤2、采用显著性分析方法，在融合后的卫星云图中，识别出具有线性特征的线性尾迹云；
步骤3、利用线性尾迹云图像的像素点统计线性尾迹云的覆盖率。
对上述技术方案的进一步改进为：所述步骤1中对卫星云图的融合方案如下：
（1）采用尺度不变特征变换算子配准红外和可见光通道卫星云图；
（2）对配准后的红外和可见光通道的卫星云图进行tetrolet变换，将红外和可见光卫星云图分解为低频和高频子图；
（3）对变换后红外和可见光通道的低频子图按照区域能量设计规则融合；
（4）对变换后红外和可见光通道的高频子图基于模糊推理设计规则融合；
（5）在tetrolet变换的变换域中利用sure无偏估计方法修改tetrolet系数对（3）和（4）融合后的图像实现去噪；
（6）再对去噪后的图像进行tetrolet逆变换，得到融合后的卫星图像。
所述步骤2的实现方案为：
（1）利用卫星资料中两个不同红外通道的亮温图像，得到两个红外通道的亮温差值图像；
（2）将亮温差值图像每个像素都重新投影到兰伯特等面积方位图上，并使用渐进概率式霍夫变换来提取线性尾迹云的边缘，再用不同二元掩模检测具有这些线性的特征的图像是否是真实的线性尾迹云，从而得到识别结果；
（3）由步骤1得到的融合卫星云图作为输入图像，再利用基于对象先验知识的显著性检测模型生成融合图像的显著性图，然后通过显著性分割提取融合图像特征来完善和优化识别结果，以保存真实的线性尾迹云像素点；
（4）对识别结果进行傅里叶频谱校正，得到线性尾迹云图像。
所述卫星的两个不同红外通道的波段分别为，10.3-11.3μm和11.5-12.5μm。
所述步骤3的实现方案为：
（1）对识别出的线性尾迹云的像素点进行分析，剔除由于环境影响及仪器误差得到的错误像素点；
（2）对于不连续的线性尾迹云像素进行亚像素插值，从而得到更完整的线性尾迹云轨迹图像；还可以进一步地利用计算机算法对上述修正过的图像进行统计和分析，得到线性尾迹云的覆盖率。
所述卫星采用风云3号卫星。
本发明具有以下优点：
（1）本发明利用融合卫星云图得可见光通道和红外通道图像来提取线性尾迹云特征因子，进行准确的线性尾迹云识别，该融合方案可实现红外和可见光通道云图融合的同时，能够有效的保留红外通道的亮温数据，可见光通道的高层语义特征以及其它与线性尾迹云识别密切相关的图像底层特征信息，并通过tetrolet变换降低了噪声的影响，为后期基于融合云图的特征提取和显著性图的构建提供了更加丰富的信息。
（2）以往的线性尾迹云识别算法仅仅采用红外通道的亮温和亮温差图像来提取和识别线性尾迹云的线性特征，辅助于人工统计分析方法进行订正，错误识别率较高；本发明采用显著性分析方法，在融合后的卫星云图中，选出具有线性特征的云，剔除错误的噪音，提取线性尾迹云的线性特征信息，这样可以消除传统线性尾迹云识别算法的局限性，通过融合图像的显著性图和显著性分割，指导和完善线性尾迹云的识别，有效的去除环境和传感器因素导致的错误匹配和不连续线性尾迹云轨迹的情况，从而可得到相对准确的线性尾迹云覆盖率和光学厚度。
（3）本发明对识别得到的线性尾迹云资料，经过进一步的筛选和增补，得到更加准确的线性尾迹云像素点信息，从而可计算得到更为准确线性尾迹云覆盖率。
附图说明
图1为本发明的结构框图。
具体实施方式
下面结合附图以及具体实施例对本发明进行详细说明。
如图1所示，本发明的基于卫星云图资料的飞机尾迹云识别方法包括如下步骤：
步骤1、基于tetrolet变换结合模糊算子的红外通道和可见光通道卫星云图融合
本步骤拟将10.3-11.3μm和11.5-12.5μm波段红外和可见光的多通道融合成一幅新的云图，研究线性尾迹云的识别方案。由于成像传感器，气候条件等因素，获得的卫星云图含有不同程度的噪声。为了提高线性尾迹云的识别准确率，需要考虑多通道卫星资料融合后的噪声处理。因此，本发明采用的卫星云图的融合方案如下：
（1）采用尺度不变特征变换算子(sift,scaleinvariantfeaturetransform)配准两个红外和可见光通道卫星云图；
（2）对配准后的两个红外和可见光通道的卫星云图进行tetrolet变换，将红外和可见光卫星云图分解为低频和高频子图；
（3）对变换后各通道的低频子图按照区域能量设计规则融合；
（4）对变换后各通道的高频子图基于模糊推理设计规则融合；
（5）在tetrolet变换的变换域中利用sure无偏估计方法修改tetrolet系数对融合后的图像实现去噪；
（6）进行tetrolet逆变换，将高频和低频子图转化为一张图像，得到融合后的卫星图像。
上述融合方案在实现红外和可见光通道云图融合的同时，能够有效的保留红外通道的亮温数据，可见光通道的高层语义特征以及其它与线性尾迹云识别密切相关的图像底层特征信息，并通过tetrolet变换降低了噪声的影响，为后期基于融合云图的特征提取和显著性图的构建提供了更加丰富的信息。
步骤2、采用显著性分析方法，在融合后的卫星云图中，识别出具有线性特征的线性尾迹云；
线性尾迹云可以单独存在，但通常会出现在自然卷云的周围。由于线性尾迹云的光学厚度较薄，在卫星云图中不易直接观测得到，需要对卫星的原始数据进行处理。本发明对传统识别方法进行了改进，具体实施方案如下：
（1）利用卫星上波段为10.3-11.3μm的红外通道和波段为11.5-12.5μm的红外通道的亮温图像，得到两个红外通道的亮温差值图像；（2）为了减少视角超过50°时的失真，将亮温差图像每个像素都重新投影到兰伯特等面积方位图，并使用阈值来提取线性尾迹云的边缘，以及用不同二元掩模检测具有这些线性的特征的图像是否是真实的线性尾迹云，从而得到初步的识别结果，但其中可能包含部分误匹配的背景像素点，以及不连续的线性尾迹云信息；
（3）由步骤1得到的融合卫星云图作为输入图像，利用基于对象先验知识的显著性检测模型生成融合图像的显著性图；再通过显著性分割提取融合图像特征来完善和优化上述识别结果，保存真实的线性尾迹云像素点，剔除背景环境的误匹配点以及建立基于种子区域的生长模型完善；
（4）对识别结果进行傅里叶频谱校正；通过傅里叶变换和反变换，尽可能地保留线性的特征，消除噪音影响，从而可有效的识别线性尾迹云；
（5）以人工主观评价的线性尾迹云识别和统计结果作为基准，与传统的仅依靠亮温差图像的线性尾迹云识别算法比较，以不断完善上述图像融合和显著性特征提取的方法，并优选一种合理的参数订正方法确定上述识别匹配所需阈值；
上述线性尾迹云识别方案，可以消除传统线性尾迹云识别算法的局限性，通过融合图像的显著性图和显著性分割，指导和完善线性尾迹云的识别，有效的去除环境和传感器因素导致的错误匹配和不连续线性尾迹云轨迹的情况，从而得到相对准确的线性尾迹云覆盖率和光学厚度。
步骤3、利用线性尾迹云图像的像素点统计线性尾迹云的覆盖率；
本步骤对融合的卫星云图和通过自动识别得到的线性尾迹云定位资料，经过进一步的筛选和增补，得到更加准确的线性尾迹云像素点信息，计算得到线性尾迹云在东亚地区的覆盖率，具有实施方案如下：
（1）结合环境场以及飞行数据，对表征线性尾迹云的像素点进行分析，剔除由于环境影响及仪器误差得到的错误像素点，比如赤道地区由于对流形成的线性薄卷云会在线性尾迹云图像中产生错误像素点；
（2）由于线性尾迹云老化和仪器精度等问题往往会造成线性尾迹云像素点的断裂和缺失，本方案对于不连续的线性尾迹云像素进行亚像素插值，从而得到更完整的线性尾迹云轨迹图像；
（3）利用纠正过的图像和计算机智能算法对修改过的图像进行统计，求得线性尾迹云的覆盖率，并统计和分析覆盖率的季节及年际变化。
为了进一步提高尾迹云识别的准确性，还可采用分辨率较高的风云3号卫星的云图资料进行分析，从而准确地识别线性尾迹云，并计算尾迹云的覆盖率从而用于其辐射特性，将会进一步完善尾迹云对区域和全球气候影响的研究，对研究尾迹云的辐射强迫有重要意义。
本发明的基于卫星云图资料的飞机尾迹云识别方法不局限于上述实施例，凡采用等同替换方式得到的技术方案均落在本发明要求保护的范围内。