本发明涉及一种基于遥感影像融合的冰川面积变化监测方法，属于对地观测与冰冻圈领域。
背景技术：
：20世纪以来全球气候变化的速度在不断加快，并在一定程度上不断加强。冰川是气候变化的指示器，对气候变化的响应极为敏感，气候变暖，冰川加速退缩，冰川泥石流将进入高潮期，可能形成冰川泥石流灾害点。冰川变化监测是全球变化和区域环境研究中不可或缺的内容之一。大部分冰川地处偏远山区，实地测量冰川面积费时费力，遥感覆盖范围广、信息量大、可连续观测的特点使它为冰川变化监测提供了快速、可靠的方法。20世纪80年代以来，卫星遥感技术在冰川变化监测中得到了广泛的应用。本世纪以来，利用遥感进行冰川监测能力随着卫星及传感器的发展不断提高。冰川面积变化的遥感监测方法主要可分为两类，分别是目视解译和计算机辅助分类法。在对研究区域先验知识比较丰富的情况下，传统目视解译方法因准确性较高，但效率很低。计算机辅助分类法为了提高解译效率，有相关学者对冰川自动提取方法进行了研究。hall等[1]最先尝试采用波段比值法区分冰雪；ye等[2]采用非监督分类方法监测了各拉丹冬冰川1969-2002年的变化；别强等[3]c采用雪盖指数方法对1960-2010年黑河流域的冰川变化进行了检测；此外，面向对象的图像信息自动分析方法也被运用到了冰川信息提取研究中[4]。不同的冰川自动提取方法均存在一定的局限性，如监督非监督分类方法可能会错分部分的积雪与岩石；雪盖指数法难以识别冰川的表碛分布区与积雪覆盖区[5]。传统冰川变化提取方法中，波段比值法的总体精度高。综上，现有的冰川变化监测方法多为间接方法，即需要先对各时期影像的整体冰川信息进行提取，然后通过计算得到冰川整体面积变化情况。这就使得设计一种直接对冰川变化区域提取进行直接提取的方法具有科学意义与应用价值。相关文献如下：[1]halldk,ormsbyjp,bindschadlerra,etal.characterizationofsnowandicereflectancezonesonglaciersusinglandsatthematicmapperdata[j].annalsofglaciology,1987,9:104-108[2]yeq,kangs,chenf,etal.monitoringglaciervariationsongeladandongmountain,centraltibetanplateau,from1969to2002usingremote-sensingandgistechnologies[j].journalofglaciology,2006,52(179):537-545[3]bieqiang,qiangwenli,wangchao,etal.monitoringtheglaciervariationintheupperreachesoftheheiheriverbasedonremotesensingin1960-2010[j].journalofglaciologyandgeocryology,2013,35(3):574—582.[4]nieyong,zhangyili,liulinshan,etal.monitoringglacierchangebasedonremotesensinginthemt.qomolangmanationalnaturepreserve,1976-2006[j].actageographicasinica,2010,65(1):13—28[5]yanlili,wangjian.studyofextractingglacierinformationfromremotesensing[j].journalofglaciologyandgeocryology,2013,35(1):110—118技术实现要素：本发明主要解决的技术问题是提供一种面向冰川变化监测的遥感影像融合方法，能够在定性分析和总体精度方面优于常用的冰川变化监测方法。针对传统冰川变换检测方法需要对不同时期研究区中冰川面积分别统计的问题，本方法针对陆地资源卫星数据提出了一种基于像元分析法的冰川变化监测分析方法，对不同时期影像预处理后进行同源数据的色调替换ihs融合，使实际冰川变化区域在融合后图像中发生颜色突变，从而实现对冰川变化的监测。为解决上述技术问题，本发明采用的一个技术方案如下：整个技术流程主要分为三个步骤，首先对两期同源landsat影像进行几何纠正和辐射校正，并将不同时期影像进行ihs分解，得到图像的亮度、色调、饱和度分量；接着通过替换色调分量获得融合图像，使冰川变化区域突显出来；最后利用决策树分类方法对冰川面积变化区域进行提取。步骤1，对两期同源landsat影像进行几何纠正和辐射校正后，通过rgb空间与ihs空间正反变换，得到图像的亮度、色调、饱和度分量；具体步骤如下：1-1)通过图像通道的dn值得到其对应的辐射亮度值，减少不同时期获得的影像的辐射亮度差异；1-2)基于外部地形数据如srtm等通过高精度几何纠正，提高不同时相的数据反映同一地物目标信息的准确度，以提高最终的监测精度。1-3)通过envi软件对遥感影像进行裁剪，获得裁剪后的覆盖冰川区域的影像，通过rgb空间与ihs空间正反变换，将两期多光谱图像变换到ihs空间，使冰川变化区域突显出来；步骤2，计算第一、二时相影像的颜色相似度，再结合两像素颜色的欧拉距离，进行新的色调分量的替换，当相似度小于一定阈值时，可认为该像元在两个时期对应的地物信息不一致，则将该像元的色调分量用第二时相该像元的色调分量值进行替换，反之，不进行替换，将i1h1's1逆变换到rgb空间，得融合图像；步骤3，颜色突变的区域为冰川减少的区域，利用决策树分类方法将变化区域提取出来，实现冰川面积变化监测。本发明的有益效果是：本发明方法较常用的分类后比较法(波段比值法、雪盖指数法、非监督分类法等)具有以下优势：(1)它能使冰川变化区域更直接直观地显示出来。传统的冰川监测方法只能得到研究区整体冰川面积变化量，而无法确定变化区域，本文方法可以将变化区域直观地显示出来，方便缩小研究区域，在分析冰川变化具体历程中有重要意义；(2)便于快速进行信息提取.传统的分类后比较法需要对研究区域所有冰川进行提取，而该方法只关注变化区域，在信息提取速度上，该方法优于传统的分类后比较法。基于改进ihs融合方法的冰川变化监测方法在以上两个方面明显优于常用的冰川变化监测法。附图说明图1是本发明的基于ihs融合的冰川变化监测技术流程图；具体实施方式下面结合附图对本发明的较佳实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。请参阅图1，本发明公开了一种面向冰川变化监测的遥感影像融合方法，包括以下步骤：1)待融合两期多光谱图像进行预处理(配准)，分别进行几何纠正和辐射校正，将两期多光谱图像变换到ihs空间；rgb空间与ihs空间正反变换可通过式(1)、(2)实现.其中，ν1、v2为中间变量，h＝tan-1(v1/v2).2)计算第一、二时相影像的颜色相似度，并根据颜色相似度得到新的色调分量。两个像素的颜色相似度公式为式(3)所示：式中，s(m1(i,j),m2(i,j))表示影像m1的第i行j列像元与影像m2的相同位置像元的颜色相似度。d为两像素的颜色的欧拉距离，可通过式(4)计算。k、l均为常数。得到颜色相似度之后，可由式(5)计算融合后图像的色调分量。3)将i1h1′s1逆变换到rgb空间，得融合图像。融合后图像中颜色突变的区域对应m1影像中是冰川而在m2中不是冰川的区域，即第二时相较第一时相冰川减少的区域。交换两个时相影像的顺序，按照上述方法进行融合，得到的融合图像中颜色突变的区域则为第二时相较第一时相冰川增加的区域，从而实现冰川变化信息的自动发现。采用改进ihs融合方法实现冰川变化信息发现之后，再利用决策树分类方法将变化区域提取出来，从而实现冰川变化监测。为了更好的说明该方法的有效性，本文也采用了传统冰川变化监测方法中总体精度最高的波段比值法对实验区冰川变化进行监测。在遥感图像处理软件envi支持下，计算可见光波段与红外波段的比值，并结合适当的阈值对两期影像的冰川进行提取，并以高分辨率遥感数据得到的结果为参考进行精度估算(表2)。表1不同方法结果对比数据源高分辨率影像低分辨率影像低分辨率影像冰川变化监测方法人工提取波段比值法基于改进ihs融合法2007年冰川面积(km2)34.3631.422013年冰川面积(km2)37.7534.40冰川面积变化量(km2)3.382.983.10精度(％)88.1191.51从表2中可以看出，针对该实验区，文中方法变化检测精度优于常用的波段比值法。以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的
技术领域：
，均同理包括在本发明的专利保护范围内。当前第1页12