顾及时点特征的水体提取成果空间修正方法

程滔， 李广泳， 毕凯

(国家基础地理信息中心，北京 100830)

0 引言

在全国范围进行水体信息提取，需要反映某一时间节点的水体状态，以便掌握水体信息在这一时间节点的数量、分布和利用现状等情况，也利于基于多次调查，分析全国范围水体信息的动态变化特征和规律，为水体综合评价与优化配置提供决策数据支撑。因此，在全国范围开展水体信息提取，标准时点统一是必然要求。目前，遥感和地理信息系统等技术和手段是水体信息提取的主要方法，高空间分辨率遥感卫星的不断增多，为精细化的水体信息提取提供了丰富的遥感影像数据源。随着遥感、地理信息系统和互联网等技术的快速发展，水体信息提取对成果的空间精度和准确性要求也不断提升[1]。例如，第二次全国土地调查以2009年12月31日为调查标准时点，所使用的遥感影像空间分辨率为一类区优于1 m、二类区优于2.5 m、三类区和四类区优于5 m[2]。第三次全国土地调查以2019年12月31日为调查标准时点，所使用的遥感影像中： 农村土地调查全面采用优于1 m空间分辨率的航天遥感数据、城镇土地利用现状调查采用优于0.2 m空间分辨率的航空遥感数据[3]。第一次全国水利普查以2011年12月31日为普查标准时点，所使用的遥感影像根据普查任务和内容的不同有所区分，主要包括2.5 m和20 m等多种空间分辨率影像[4]。第一次全国地理国情普查(以下简称“地理国情普查”)以2015年6月30日为普查标准时点，并且在持续的年度地理国情监测中，均以6月30日为标准时点，所使用的遥感影像空间分辨率为全国优于2.5 m以及重点区域优于1 m[5]。

然而，我国地域广阔，气象条件差异性明显，在全国范围获取一遍满足标准时点要求的、统一时相的、质量可行的高空间分辨率遥感影像比较困难，这给整体时点统一要求带来了不确定性因素。对于水体信息提取而言，基于多种时相遥感影像采集水体信息，会导致部分区域影像接边处水体覆盖范围不一致，统计得到的全国水体面积数据存在一定误差。目前，较好的解决方法是先基于已获取遥感影像提取各项信息，待获取到标准时点遥感影像后，对前期调查成果进行更新，将其统一到标准时点上； 或结合外业调查的方法进行更新。但高空间分辨率的标准时点遥感影像的可获取性也存在不确定性，外业调查方法投入的人力、物力较大。相比之下，空间分辨率相对较低、时间分辨率相对较高的遥感影像比较容易获取。

因此，针对地表覆盖水体接边区域遥感影像时相不满足标准时点要求的情况，在原始成果基础上，收集获取空间分辨率相对较低、时间分辨率较高的遥感影像，开展水体信息提取； 同时，基于精细格网数字高程模型(digital elevation model，DEM)数据，利用区域生长算法，获取精细化的水体成果，对原始成果进行地理空间修正，从而保证水体提取成果既满足标准时点要求，又满足精度要求； 并以地理国情普查地表覆盖水体成果为研究实例，开展方法应用与效果验证。

1 研究方法

1.1 水体信息提取方法

水体信息提取的关键在于水体特征分析与挖掘，以及水体特征规则构建，在大部分遥感影像上，水体的纹理一般比较均匀、平滑，与周边地表覆盖物光谱差异较大。在基于遥感影像的地表覆盖信息提取中，归一化差异水体指数(normalized difference water index，NDWI)是水体信息自动提取的算法之一[6]。该算法是通过寻找水体的最强和最弱反射波段，从而增强水体对象，抑制背景地物，实现水体信息提取的目的。

据研究，水体在绿光和近红外波段上，分别表现出了强反射和强吸收特征，故NDWI公式为

NDWI=(G-N)/(G+N)，

(1)

式中G和N分别为绿光和近红外波段的亮度值。

NDWI算法在基于中、低空间分辨率遥感影像的水体信息提取中具有非常好的适用性。此外，水体在蓝光波段上表现出较强反射特性。因此，可将蓝光波段的比率值和标准方差值作为水体信息提取的其他主要判定规则[7]。

根据处理单元的粒度，水体信息提取方法可分为基于像元光谱统计的自动分类方法和面向对象语义信息的自动分类方法。前者适用于中、低空间分辨率遥感影像； 后者适用于高空间分辨率遥感影像。

1.2 区域生长算法

在基于精细格网DEM数据对水体信息进行优化处理中，区域生长算法是较适用的方法。该算法依据区域联通原理，将有相似性质的像素点合并到一起形成一个区域。具体地说，是对每一个生长区域，先指定一个种子点作为生长的起点，然后将种子点周围邻域的像素点和种子点进行对比，将具有相似性质的点合并起来继续向外生长，直到没有满足条件的像素被包括进来为止，这样一个区域的生长就完成了。区域生长算法过程有3个关键问题，即种子点选取方法[8-9]、生长过程中能将相邻像素包括进来的准则确定方法和区域生长的停止条件设置方法。

本研究在确定生长准则时，通过水体结果与DEM数据空间叠置分析，依据水体图斑内各格网DEM各项统计值确定，通过设置阈值范围，确定相邻像素是否被包括进来，基本准则为

Dh≤∂STD，

(2)

式中：Dh为相邻像素高程值与种子点高程值的差值； ∂为调整系数；STD为标准差。

由于缺乏商业软件支撑，数据处理存在难度，因此本研究在Microsoft Visual Studio 2010集成开发环境中，利用C++语言[10]，自主研发了地表覆盖水体优化计算软件。通过设计区域生长算法，结合地表覆盖水体结果和精细格网DEM数据的特点，解决了该算法的3个关键问题。开发过程中，调用了开源栅格空间数据转换库(geospatial data abstraction library，GDAL)[11]，利用其栅格空间数据处理能力，实现栅格数据的读出、空间计算和写入。软件主界面如图1所示。

图1 地表覆盖水体优化计算软件主界面

该软件的研发提高了数据处理的自动化水平和成果质量，实现了将水体成果优化至高空间分辨率水平的目的，满足了研究中数据处理工作对软件工具的需求。

2 研究区概况及数据源

以地理国情普查地表覆盖水体成果为研究实例，选取了长江流域的局部区域作为研究区。该研究区位于省级行政区划交界处，面积为37.26 km2。影像接边两侧均为WorldView-2影像，空间分辨率为0.5 m； 东北部影像时相为2014年1月4日，其余部分影像时相为2013年1月19日，两者距离地理国情普查标准时点均较远； 不符合地理国情普查标准时点要求。使得图2(a)黄色矩形框范围，即影像接边处水体范围变化明显。经统计，研究区地理国情普查地表覆盖水体成果总面积为11.91 km2，范围如图2(a)红线所示。另外，研究区属于山地、高山地地形，高程范围为417.90～1 311.63 m。从图2(b)的三维视图可以看出，研究区沿河流方向呈四周地势高、中间地势低的特点。DEM数据来源于地理国情普查精细化生产[12]，格网单元尺寸为2 m，等高距为5 m，是基于1∶ 1万数字线划图(digital line graphic，DLG)数据线性内插生成，高程中误差为1.5 m，数据源现势性为2012年12月。

(a) 地理国情普查影像及水体成果 (b) 三维视图

本研究收集了研究区符合地理国情普查标准时点要求的Landsat8遥感影像数据，数据来源于地理空间数据云网站(http: //www.gscloud.cn/)，空间分辨率为15 m，影像时相为2015年6月29日，经波段组合和彩色增强处理后的结果如图3(a)所示。虽然此影像时相最接近地理国情普查时点，但影像云量较大，对水体信息提取具有一定的影响。因此，另收集了一景2014年1月17日Landsat 8遥感影像数据，如图3(b)所示。

经数据配准、对比分析发现，2幅影像水体分布范围吻合度极好，因此，研究中采用2014年1月17日影像代替2015年6月29日影像开展研究。

3 结果与分析

3.1 基于Landsat8遥感影像数据的水体信息提取

基于2014年1月17日Landsat8遥感影像数据计算得到NDWI指数数据(图4)，反演得到水体提取结果。

图4 NDWI指数计算结果

将水体提取结果分别于2幅Landsat8遥感影像数据叠合(图5)，发现2幅影像水体分布范围吻合度极好，也进一步证明了采用2014年1月17日影像代替2015年6月29日影像开展研究的可行性。

(a) 2015年6月29日 (b) 2014年1月17日

3.2 基于2 m格网DEM数据的水体成果精度优化

由于研究区属于山地、高山地地形，高程落差较大。因此，考虑水体落差因素，研究中将最大图斑进行分段处理，共计分为13段，提取出各段几何中心点13个，作为待定种子点。

通过与DEM空间叠置分析，提取出各段几何中心点的高程值，同时，统计出各段高程最小值、最大值、平均值、值域范围和标准差等统计值。

通过对各段几何中心点的高程值与各项统计值进行对比分析，对5个待定种子点的空间位置进行了修正，最终确定的种子点空间分布范围如图6所示。

图6 水体种子点空间分布范围

采用8邻域算子，通过图形区域生长，对水体种子点进行迭代生长计算，∂系数取1.5，并进行二值化处理，得出基于精细格网DEM的区域生长结果，如图7所示。

图7 基于精细格网DEM的水体区域生长结果

3.3 结果对比分析

地理国情普查水体成果、基于Landsat8遥感影像数据提取的水体结果和区域生长结果数据对比显示如图8所示。

图8 地理国情普查水体成果、基于Landsat8遥感影像数据提取的水体结果和区域生长结果数据对比

从图8可以看出，研究区区域生长结果与基于Landsat8遥感影像数据提取的水体结果图斑吻合度较高，两者与地理国情普查水体成果之间都存在一定差异，通过水体优化处理，原始成果得到了修正。

对研究区水体成果优化的各项指标进行统计，结果如表1所示。

表1 水体成果修正各项指标统计结果

从表1可以得出，研究区水体成果空间范围相对于标准时点监测修正了17.97%，通过2 m格网DEM数据对基于15 m空间分辨率Landsat8遥感影像数据提取的水体结果进行精度优化，研究区水体成果空间范围优化了1.51%。通过该方法，实现了水体成果空间范围修正。

4 结论

1)采用较高时间分辨率遥感影像与精细格网DEM数据协同的方式，开展水体提取成果空间修正，能够对大范围的标准时点高空间分辨率遥感影像覆盖形成较好的影像补充，该方法保证了水体提取成果的现势性满足标准时点要求； 同时，精细格网DEM的应用，能够使得基于低空间分辨率遥感影像采集的信息精度保持在较高的水平。从而在大尺度上形成最接近标准时点状态的调查成果，缩小了由于影像数据源的差异造成的调查成果误差，使统计结果更加客观、合理。

2)通过选取地理国情普查地表覆盖水体成果开展方法应用与效果验证，验证了该方法的适用性； 地表覆盖水体优化计算软件的研发，提高了数据处理的自动化水平和成果质量，满足了数据处理工作的需要，提出的方法能够在水体提取成果优化中推广应用。