一种用于白云石信息提取的高光谱影像处理方法
【专利摘要】本发明涉及一种用于白云石信息提取的高光谱影像处理方法，保留白云石的光谱特征明显的波段，提取波段在1820nm、1985nm、2015nm、2195nm、2315nm、2380nm、2420nm的影像，进行一系列判断和计算，计算出影像范围内不同区域白云石的丰度值。这种方法可以去除其他特征不明显的波段，从而在信息提取的过程突出白云石的光谱特征，降低其他地物或噪声的影响，减少了处理的数据量，并可以用IDL程序达到用较少的人工操作实现白云石信息提取的目的，提高了白云石信息提取的精度和速度。
【专利说明】一种用于白云石信息提取的高光谱影像处理方法
【技术领域】
[0001]本发明涉及一种用于白云石信息提取的高光谱影像处理方法，特别是涉及一种减少其他地物或噪声的影响、人工操作步骤和处理数据量的用于白云石信息提取的高光谱影像处理方法。
【背景技术】
[0002]当前的高光谱遥感影像的白云石信息提取方法主要是光谱全波段匹配或是部分连续波段的光谱匹配，具体算法有光谱角、混合解调滤波等，由于地表的物质组成很少是由单一矿物组成的，这些方法在白云石信息提取的过程易受其他地物波谱或噪声的影响，提取信息精度相对较低。其次现有的光谱提取方法人工操作步骤多，增加了人为判断误差。第三是高光谱数据波段多，数据量大，现有的方法处理时间长，降低了数据处理的速度和应用规模。因此，如何在白云石信息提取的过程中减少其他地物或噪声的影响、人工操作步骤和处理数据量，成为当前高光谱遥感影像处理的前沿之一。
【发明内容】
[0003]本发明要解决的技术问题是提供一种减少其他地物或噪声的影响、人工操作步骤和处理数据量的用于白云石信息提取的高光谱影像处理方法。
[0004]为解决上述技术问题，本发明一种用于白云石信息提取的高光谱影像处理方法，包括以下步骤:
[0005]第一步，获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；
[0006]第二步，对地面反射率的影像数据进行重采样，提取波段在1820nm、1985nm、2015nm、2195nm、2315nm、2380nm、2420nm 的影像；
[0007]第三步，当1820nm波段影像的像元值大于1985nm波段影像相对应的像元值，判断该像元满足条件一；获得1820nm波段影像减去1985nm波段影像的结果影像Hl ；
[0008]第四步，当1985nm波段影像的像元值小于2015nm波段影像相对应的像元值，判断该像元满足条件二 ；获得2015nm波段影像减去1985nm波段影像的结果影像H2 ；
[0009]第五步，当2195nm波段影像的像元值大于2315nm波段影像相对应的像元值，判断该像元满足条件三；获得2195nm波段影像减去2315nm波段影像的结果影像H3 ；
[0010]第六步，当2315nm波段影像的像元值小于2380nm波段影像相对应的像元值，判断该像元满足条件四；获得2380nm波段影像减去2315nm波段影像的结果影像H4 ；
[0011]第七步，当2380nm波段影像的像元值大于2420nm波段影像相对应的像元值，判断该像元满足条件五；获得2380nm波段影像减去2420nm波段影像的结果影像H5 ；
[0012]第八步，选择同时满足条件一至条件五的像元，获得上述像元范围内的H1、H2、H3、H4、H5相加的和H。
[0013]本发明可以减少处理的数据量，SASI的波段数量为101个波段，本方法只用于7个波段，数据量减少了 94%，并且由于是计算机自动一步提取，减少了主成分变换、端元波谱的选择等操作步骤，运算速度可以提高了 13倍以上。由于去除了大部分对信息提取关系不大的波段，减少其他物质或噪声对其光谱的干扰，提高了信息提取的精度。对高光谱影像数据中白云石信息的快速提取具有较好的作用和意义。
【具体实施方式】
[0014]本发明通过对高光谱影像重采样，提取特定波段，进行一系列判断和计算，计算出影像范围内的不同区域白云石的丰度值。
[0015]具体包括以下步骤:
[0016]第一步，获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；
[0017]第二步，对地面反射率的影像数据进行重采样，提取波段在1820nm、1985nm、2015nm、2195nm、2315nm、2380nm、2420nm 的影像；
[0018]第三步，当1820nm波段影像的像元值大于1985nm波段影像相对应的像元值，判断该像元满足条件一；获得1820nm波段影像减去1985nm波段影像的结果影像Hl ；
[0019]第四步，当1985nm波段影像的像元值小于2015nm波段影像相对应的像元值，判断该像元满足条件二 ；获得2015nm波段影像减去1985nm波段影像的结果影像H2 ；
[0020]第五步，当2195nm波段影像的像元值大于2315nm波段影像相对应的像元值，判断该像元满足条件三；获得2195nm波段影像减去2315nm波段影像的结果影像H3 ；
[0021]第六步，当2315nm波段影像的像元值小于2380nm波段影像相对应的像元值，判断该像元满足条件四；获得2380nm波段影像减去2315nm波段影像的结果影像H4 ；
[0022]第七步，当2380nm波段影像的像元值大于2420nm波段影像相对应的像元值，判断该像元满足条件五；获得2380nm波段影像减去2420nm波段影像的结果影像H5 ；
[0023]第八步，选择同时满足条件一至条件五的像元，获得上述像元范围内的H1、H2、H3、H4、H5相加的和H。
[0024]H的值代表像元白云石的丰度值，H值越大，代表像元内白云石的丰度越大，SP含
量越高。
【权利要求】
1.一种用于白云石信息提取的高光谱影像处理方法，包括以下步骤: 第一步，获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据； 第二步，对地面反射率的影像数据进行重采样，提取波段在1820nm、1985nm、2015nm、2195nm、2315nm、2380nm、2420nm 的影像； 第三步，当1820nm波段影像的像元值大于1985nm波段影像相对应的像元值,判断该像元满足条件一；获得1820nm波段影像减去1985nm波段影像的结果影像Hl ； 第四步，当1985nm波段影像的像元值小于2015nm波段影像相对应的像元值，判断该像元满足条件二 ；获得2015nm波段影像减去1985nm波段影像的结果影像H2 ； 第五步，当2195nm波段影像的像元值大于2315nm波段影像相对应的像元值，判断该像元满足条件三；获得2195nm波段影像减去2315nm波段影像的结果影像H3 ； 第六步，当2315nm波段影像的像元值小于2380nm波段影像相对应的像元值,判断该像元满足条件四；获得2380nm波段影像减去2315nm波段影像的结果影像H4 ； 第七步，当2380nm波段影像的像元值大于2420nm波段影像相对应的像元值,判断该像元满足条件五；获得2380nm波段影像减去2420nm波段影像的结果影像H5 ； 第八步，选择同时满足条件一至条件五的像元，获得上述像元范围内的H1、H2、H3、H4、H5相加的和H。
【文档编号】G06T5/00GK103903225SQ201210573014
【公开日】2014年7月2日 申请日期:2012年12月25日 优先权日:2012年12月25日
【发明者】杨燕杰, 赵英俊 申请人:核工业北京地质研究院