一种适用于高山峡谷地区的遥感智能解译方法与流程

1.本发明涉及遥感数据解译技术领域，特别地涉及一种适用于高山峡谷地区的遥感智能解译方法。


背景技术：

2.水资源是地球生态系统中不可或缺的重要资源，是气候变化和生态系统演化的重要影响因素。快速精准地识别水域，对水资源调查、洪水监测及防灾减灾等有重要意义。随着遥感技术的发展，高分辨率遥感成为目前水体监测的重要手段。目前，人们广泛探索了从卫星图像中提取水域的方法。已经开发了几种方法来从空间、光谱和时间特征不同的不同卫星图像中描绘水域。用卫星图像绘制森林和水域等自然资源的地图变得越来越重要。森林和水资源都受到强烈开发，定期监测对于可持续管理至关重要。在全球碳循环和气候变化中发挥关键作用的水域在时空域中被绘制出来，以分析和评估其退化和消失的程度和速度。
3.世界各地的研究人员已经使用了不同的卫星数据在空间、光谱和时间特征的变化，以产生土地利用土地覆盖的专题地图或特别强调水域。同时，已经采用了各种技术从卫星图像中提取这些特征，每种方法都有其优点和缺点。卫星数据的目视解译提供了对不同大小水域的最佳描绘，但很耗时，尤其是在处理高分辨率数据时。当水域与其他类别之间存在光谱重叠时，使用交互式自组织数据分析技术这类简单且常见的无监督分类方法，提供的结果精度和解译效率非常低。


技术实现要素：

4.针对上述现有技术中的问题，本技术提出了一种适用于高山峡谷地区的遥感智能解译方法，使用至少5个波段的数据，基于线性预测函数分类，所述线性预测函数结合了特征向量，这些特征向量与权重相结合以构建预测函数；在找到特征向量并初始化权重后，计算加权和；所述加权和作为二元分类器输出函数的输入，利用所述二元分类器的输出函数，将所述输入“x”映射到一个值“y”，使用硬限功能执行分类过程，在执行分类过程时使用临界值；所述临界值取决于所述特征向量和所述权重，选择物体像素显示的加权和的最大可能值，将每个像素获取的加权和与临界值进行比较，从而完成分类。
5.优选地，特征向量通过找出作为物体特征的描述符来定义。
6.优选地，权重取决于对象的空间和时间属性。
7.优选地，每个向量的权重函数是关于时间和位置的函数gi(r,t)，通过统计大量不同时间、不同地点、不同物体的反射率值的变化数据，拟合出每个向量的权重函数gi(r,t)，通常情况下权重函数gi(r,t)的数值与标准时间严格相关，如果遇到多个像素点反射率实际测量值偏离历史同期数据时，允许计算机根据反射率实际测量值自动匹配另一时间段的权重函数曲线，以适应特例情况。
8.优选地，特征向量包括短波红外波段物体的反射辐射值和反射波段中不同物体的
反射率值的变化规律。
9.优选地，根据反射波段中不同物体的反射率值的变化规律，对绿地、建筑物和非绿地的智能提取。
10.上述技术特征可以各种适合的方式组合或由等效的技术特征来替代，只要能够达到本发明的目的。
11.本发明提供的一种适用于高山峡谷地区的遥感智能解译方法，与现有技术相比，至少具备有以下有益效果：本发明对基于遥感技术获取的大量数据进行智能、高效解译，综合考虑了多种与被测对象属性相关的特征，兼顾了时间、空间、特例水文气象等特征，使得高山峡谷地区的遥感数据解译更加智能、更加准确。遥感数据量庞杂繁巨，按照土地类型的变化情况，优选关注区域，才能获得更加高效、优质的解译结果。
附图说明
12.在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
13.图1显示了本发明的感知器模型示意图；
14.图2显示了本发明的水域和其他土地特征的反射率值示意图；
15.图3显示了反射波段中不同物体的反射率值示意图；
具体实施方式
16.下面将结合附图对本发明作进一步说明。
17.本发明提供了一种适用于高山峡谷地区的遥感智能解译方法，使用至少5个波段的数据，即蓝色波段、绿色波段、红色波段、近红外波段和短波红外波段，基于线性预测函数分类，该函数结合了特征向量。特征向量是表示对象的属性；特征向量越多，分类过程就越容易。这些特征向量与权重相结合以构建预测函数。利用二元分类器，将输入“x”映射到一个值“y”使用硬限制功能(如图1所示)。
18.w为权重向量，w
·
x为点积。要使用感知器执行分类，第一步是通过找出作为水域特征的描述符来定义特征向量。在当前的研究中，仅使用基于像素的信息。一旦定义了向量特征，下一步就是初始化每个向量的权重。这些权重取决于对象的空间和时间属性。水域显示的波段反射值在不同季节会有所不同，因为水位增加、沉积物沉积和其他杂质。因此，每个向量的权重函数是关于时间和位置的函数gi(r,t)，通过统计大量不同时间、不同地点、不同物体的反射率值的变化数据，拟合出每个向量的权重函数gi(r,t)，通常情况下权重函数gi(r,t)的数值与标准时间严格相关，如果遇到多个像素点反射率实际测量值偏离历史同期数据时，允许计算机根据反射率实际测量值自动匹配另一时间段的权重函数曲线，以适应特例情况(数十年一遇的干旱、洪水、严寒等)。
19.在找到特征向量并初始化权重后，计算加权和。这个加权和作为输出函数的输入。输出函数是一个硬限制二元函数，在执行分类过程时它使用临界值。临界值取决于特征向量和权重。在这种情况下，我们选择水域像素显示的加权和的最大可能值，并使用极值来计算该最大可能值。为每个像素获取的加权和与临界值进行比较，这样就完成了分类。
20.特征向量一：短波红外波段水的低反射辐射
21.在遥感图像中捕获的物体的反射辐射取决于物体吸收的电磁辐射的程度；也就是
说，它吸收的越多，反射的就越少。水域在红外区域吸收更多。因此，水在红外和短波红外波段的反射辐射低于其他物体，如绿地、建筑物和非绿地等。非绿地包括裸露的土地、道路、沙地等，水域和地物部分像素的反射率值如图2所示。
22.进一步从图3可以清楚地看出，水域在所有红外波段(包括第五波段b5，即短波红外波段)中都显示出明显的低反射辐射值。
23.特征向量二：反射波段中不同物体的反射率值的变化规律
24.第五波段b5中水的反射率小于第二波段b2中水的反射，但不适用于其他对象，如非绿地、绿地或建筑物。因此，通过应用这个条件，可以很容易地提取水域。
25.第一波段b1中非绿地的反射率小于第二波段b2中非绿地的反射，但不适用于其他对象，如水域、绿地或建筑物。因此，通过应用这个条件，可以很容易地提取非绿地。
26.第三波段b3中绿地的反射率小于第四波段b4中绿地的反射，但不适用于其他对象，如水域、非绿地或建筑物。因此，通过应用这个条件，可以很容易地提取绿地。
27.第二波段b2中建筑物的反射率小于第一波段b1中建筑物的反射，且第二波段b2中建筑物的反射率小于第三波段b3中建筑物的反射，但不适用于其他对象，如水域、非绿地或绿地。因此，通过应用这个条件，合理设置第二波段b2与第三波段b3的差值临界值，可以很容易地提取建筑物。所述差值临界值是关于时间和位置的函数gi(r,t)，通过统计大量不同时间、不同地点、不同物体的反射率值的变化数据，拟合出差值临界值函数gi(r,t)，不同的差值临界值函数gi(r,t)之间相互关联，例如，通常情况下差值临界值函数gi(r,t)的数值与标准时间严格相关，如果遇到多个反射率实际测量值偏离历史同期数据时，允许计算机根据反射率实际测量值自动匹配另一时间段的差值临界值函数曲线，以适应特例情况(数十年一遇的干旱、洪水、严寒等)。
28.特征向量三：像素点位置的相对海拔高度与邻近区域的反射率值
29.高山峡谷地区的水域通常位于山体之间，两侧山体与峡谷共同构成天然的容器，同一容器内水域通常位于相对海拔较低的区域，解译过程中如若计算机判定同一山体同一连续坡面中海拔较高处有水域存在，那么推断同一山体同一连续坡面中海拔较低处也有水域存在；同理，解译过程中如若计算机判定同一山体同一连续坡面中海拔较低处没有水域存在，那么推断同一山体同一连续坡面中海拔较高处也没有水域存在。
30.此外，在同一个地区中，根据该区域不同土地覆盖类型对其生态系统服务功能状况的影响和作用，可以科学合理地对土地覆盖类型的生态级别进行分类。但对于不同类型的区域类型来说，不同的土地覆盖变化类型对其生态系统服务功能变化的影响和作用不同。因此，本技术针对不同的区域类型，分别对各种土地覆盖类型的转换趋势进行了计算。利用转换矩阵表征各种土地类型转换趋势，表达式如表1所示。
31.表1.土地利用转换矩阵的描述
[0032][0033]
其中，t1和t2为不同的提取期；ty
11-ty
nn
是不同的土地利用类型；tr
1-trn是t1到t2的转移区；ty
n
到ty
n
表示t1和t2期间地面n的总面积。
[0034]
计算某一土地利用类型的转移面积占转移前总面积的比例，代表转移期内的单一活跃度，表达式如下：lj＝(s
j-s
jj
)/s
j t，其中，j为某一土地利用类型；sj是转移前的面积；s
jj
为未转移面积，单位为km2；t代表持续时间。lj值反映变化前后两个时期之间某类型面积的净变化。计算机自动筛选lj值较大的区域作为第一优先关注的目标区域。
[0035]
为客观反映土地利用的动态变化，引入综合活跃度，即流转后某类土地总面积与流转前的比值，计算公式如下：uj＝(sa sc)/s
j t，其中，sa为某一类型转换为其他类型的面积，sc表示其他土地利用类型转换为该类型的面积。uj值大时，表明土地利用综合活跃程度大，土地利用流转频繁。计算机自动筛选(lj uj)值较大的区域作为第二优先关注的目标区域。
[0036]
为了反映区域特征并减少信息冗余，选择了5个空间度量来描述研究的区域。
[0037]
表2.本发明的区域特征描述
[0038][0039][0040]
以斑点数、斑点平均面积、斑点密度、周长面积指数和最大斑点比例五个指标作为雷达图的五个要素，自动生成对应的雷达图，按照雷达图中五边形面积值，从大到小排列图像顺序，计算机自动筛选五边形面积值较大的区域作为第三优先关注的目标区域。
[0041]
虽然在本文中参照了特定的实施方式来描述本发明，但是应该理解的是，这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是，可以对示例性的实施例进行
许多修改，并且可以设计出其他的布置，只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是，可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是，结合单独实施例所描述的特征可以使用在其他所述实施例中。