本发明涉及模式识别中基于谱分解打分的区域标注方法，特别适用于modis卫星遥感图像的冰海云概率预测和区域标注问题。
背景技术：
以南极科考为重点的极地科学考察目前已有51个国家参与，包括世界大多数发达国家和主要发展中国家，它关系着全球变化和人类的未来，也是一个国家综合国力、高科技水平在国际舞台上的展现和角逐，在政治、科学、经济、外交、军事等方面都有其深远意义与重大影响。随着科学技术的发展，近年来，各国重视利用最先进的技术，为极地科考提供保障，用技术为科考保驾护航。2014年1月，在中国第30次南极考察队完成对俄罗斯“绍卡利斯基院士”号被困52名乘客救援工作后，“雪龙”船准备回撤时受阻被困在密集浮冰区，五天后才成功驶出乱冰区，进入清水区中航行。如何利用遥感图像，有效判断出冰海云区域，找到冰间水道，为极地考察船只规划合理航线，是现在科考面临的重要问题之一。到目前为止，遥感图像分割和分类主要采用选取特征阈值和目视解译的方法进行处理。自动进行遥感图像处理及区域预测，为极地科考提供有效的冰海云区域标注辅助，减少科考船只航行过程中受到浮冰阻碍的影响，是亟待解决的问题。
要达到自动进行modis卫星遥感图像概率预测与区域标注的目的，需要识别图像中有效区分冰海云区域的纹理信息，实验发现，不同区域图像块谱特征分布、综合光密度有明显差异，其中海区域图像块谱能量集中，光密度低；厚冰和厚云区域图像块谱能量相对集中，光密度高；薄冰和薄云区域谱能量分散，光密度分布广；且存在冰(云)间水道的局部区域的近似秩存在各向异性，不同方向上计算出的近似秩方差较大。利用这些特性，我们采用局部谱分解的方法提取图像块特征，进而训练分类模型对冰海云分布进行概率预测，从而达到区域标注的目的。
技术实现要素：
发明目的：目前的modis卫星遥感图像分割和分类主要采用选取特征阈值和目视解译的方法进行处理，这些方法需要大量的人工步骤，如何对modis卫星遥感图像进行自动处理，得到图像的区域标注，是亟待解决的问题。针对这个问题，本发明提供一种基于局部谱分解打分的modis卫星遥感图像标注方法，具体来说，首先收集一定量的modis卫星遥感图像并人工标注出待分类预测的区域，然后对图像进行处理、分割和提取特征并训练分类器，最后在实际使用中利用训练完毕的分类器对待预测的图像块进行分类，进而得到完整的modis卫星遥感图像区域标注。
技术方案：一种基于局部谱分解打分的modis卫星遥感图像标注方法，包括人工标注训练样本步骤，样本特征提取与分类器训练步骤，modis卫星遥感图像概率预测与区域标注步骤：
所述人工标注训练样本步骤具体为：
步骤100，根据要待标注的modis卫星遥感图像的时间和经纬度的大致范围，收集历年相同时间段，相近范围的大量modis卫星遥感图像；
步骤101，对收集的modis卫星遥感图像频段1做极地方位投影转换，并做同态滤波和直方图均衡化处理；
步骤102，对步骤101处理后的modis卫星遥感图像频段1，在modis29工具辅助下，初步勾勒出图像中的云区域，冰区域和海区域
步骤103，对步骤102处理后的modis卫星遥感图像频段1，分别针对云区域和冰区域，人工鉴别厚冰、厚云、薄冰、薄云；
步骤104，对步骤101处理后的modis卫星遥感图像频段1，随机生成一组坐标点，提取以该点为中心的20*20大小的图像块，将步骤103中人工鉴别出的标记作为提取出的图像块的标记，并将厚冰区域和厚云区域归为同一类，薄冰区域和薄云区域归为同一类，海区域单独归为一类，从而得到大量带标记的20*20大小的频段1图像块；
步骤105，对收集的modis卫星遥感图像频段2，重复步骤101到步骤104，生成带标记的20*20频段2图像块，并与带标记的20*20频段1图像块一一对应，2个一组作为训练数样本。
所述样本特征提取与分类器训练的具体步骤为：
步骤200，以步骤104中提取并手工标记的20*20频段1图像块作为输入数据；
步骤201，将20*20的频段1图像旋转7个方向，分别是0°，15°，30°，45°，60°，75°，90°，得到对应的7张图像块，对每张图像块的灰度矩阵进行谱分解，每个矩阵得到其对应的一组特征值；
步骤202，针对步骤201中生成的一个矩阵的全部特征值，把这些特征值按降序排列，按顺序累加特征值，直到当前特征值之和达到所有特征值之和的90％为止，记录当前所需特征值的个数，记作ni，ni即为所处理矩阵对应的近似秩；
步骤203，针对步骤202中生成的近似秩，按序拼接7个方向近似秩的值，得到20*20频段1图像块对应的一个长度为7的向量；
步骤204，针对20*20频段2图像块，重复步骤200-203，生成20*20频段2图像块对应的一个长度为7的向量，并与步骤203中生成的向量拼接，作为近似秩特征
步骤205，提取频段1与频段2图像块的综合光密度的均值、方差，作为光密度特征，与近似秩特征拼接，获得一对图像块的完整特征表示；
步骤206，利用人工标注的训练样本的特征表示和标记，使用rbf-svm训练分类模型，选取最佳分类参数，得到分类器c。
所述modis卫星遥感图像冰海云概率预测与区域标注的具体步骤为：
步骤300，收集待预测的modis卫星遥感图像
步骤301，对modis卫星遥感图像频段1和频段2做极地方位投影转换，并做同态滤波和直方图均衡化处理；
步骤302，对步骤301处理后的modis卫星遥感图像频段1和频段2，每隔10个像素点为中心，截取20*20区域，针对该区域提取近似秩特征和光密度特征，获得一对图像块的完整特征表示；
步骤303，使用已训练完成的分类器c对每个图像块的完整特征表示进行分类，得到图像块的标记的概率分布，作为图像块的标注打分；
步骤304，训练区域分类器，再结合前向插值，得到每个像素点的概率预测结果。对步骤301处理后的modis卫星遥感图像中每个像素点，将包含该点的小图像块的概率预测值进行平均，得到该点最终的海、厚冰厚云、薄冰薄云的概率预测结果；
步骤305，根据对每个像素点的概率预测，生成热力图，即把每点是云海冰的概率用颜色表现出来，越蓝表示是该点是海的可能性越大，越安全，越红表示该点是冰的可能性越大，越危险。
附图说明
图1为本发明的人工标注训练样本步骤流程图；
图2为本发明的样本特征提取与分类器训练步骤流程图；
图3为本发明的modis卫星遥感图像概率预测与区域标注步骤流程图。
具体实施方式
下面结合具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
人工标注训练样本阶段的工作流程如图1所示。首先收集极地区域的modis卫星遥感图像(步骤100)；然后对其进行极地方位投影转换，并做同态滤波和直方图均衡化处理(步骤101)；初步勾勒出图像中冰、海、云的区域(步骤103)；对云区域和冰区域，仔细鉴别区出厚冰、薄冰和厚云、薄云(步骤104)；对频段1图像提取大量20*20图像块，同时记录图像块的标记(步骤105)；对频段2图像提取大量20*20图像块，同时记录图像块的标记，与同位置的频段1图像块配对(步骤106)。
样本特征提取与分类器训练的工作流程如图2所示。以步骤一生成的图像块对作为输入(步骤200)，对频段1图像块进行7个方向旋转，分别是0°，15°，30°，45°，60°，75°，90°，得到对应的7张图像块，对其灰度矩阵做谱分解(步骤201)；对于各个矩阵的谱分解结果，计算其90％能量谱个数，作为近似秩(步骤202)；拼接各个方向近似秩的值(步骤203)；对配对的频段2图像块做相同处理，并与频段1图像块近似秩整合得到近似秩特征(步骤204)；对频段1和频段2图像块的综合光密度的均值、方差，作为光密度特征，与近似秩特征拼接，获得一对图像块的完整特征表示(步骤205)；利用人工标注的训练样本的特征表示和标记，使用rbf-svm训练分类模型，选取最佳分类参数，得到分类器c(步骤206)
modis卫星遥感图像概率预测与区域标注的工作流程如图2所示。首先以待预测modis卫星遥感图像作为输入(步骤300)，对其进行极地方位投影转换，并做同态滤波和直方图均衡化处理(步骤301)；对处理后的modis卫星遥感图像频段1和频段2提取大量20*20图像块，对提取的图像块进行特征化(步骤302)；将图像块的特征表示应用于训练好的分类器c，得到图像块的概率预测结果，作为该图像块中每个像素的预测结果(步骤303)；对待遇测的modis图像上每个像素的各个预测结果取均值，整合出完整图像预测结果(步骤304)；利用完整的逐像素预测结果，生成航运危险程度热力图(步骤305)。