一种大范围淤地坝位置确定方法与流程

本发明涉及一种大范围淤地坝位置确定方法，是一种使用计算机技术的水利设施低成本快速调查方法。

背景技术：

淤地坝是一种重要的水土保持工程措施，在世界范围内均有广泛的应用。据统计，仅在中国的黄河流域2003-2010年就建设了8180座淤地坝。这些淤地坝对河流流域的灌溉、防洪、水土保持、减少河流泥沙等方面做出了重要的贡献。然而，随着淤地坝使用时间的延伸，淤地坝及周围流域情况的不断变化，以及早期淤地坝建设的年代还没有成熟的数字技术，众多淤地坝的数字化信息或没有或残缺不全，难以实现数字化管理。如何补齐这些缺失的淤地坝数字化信息，需要进行缜密的调查，详细的了解各个淤地坝的状态，构建各个淤地坝的数字化模型，实现淤地坝的数字化管理。然而要实现这一详细的调查，现有的方式是采用人工调查的方式，对淤地坝进行现场测量，这需要耗费大量的人力物力。尽管有一些利用数字化地形图dem进行淤地坝调查的方案，如对dem进行三角化等方式，但这些方案都存在这样或那样的缺点，如分析过程复杂，没有充分的利用dem数据，使用的计算机资源过多，使计算成本过高等。由于这些问题的存在，使利用计算机进行淤地坝调查停留在方案上，在实际中应用很少。

技术实现要素：

为了克服现有技术的问题，本发明提出了一种大范围淤地坝位置确定方法。所述的方法结合卫星图片和dem数据，快速辨认淤地坝的位置，实现大范围的淤地坝的快速确定。

本发明的目的是这样实现的：一种大范围淤地坝位置确定方法，所述方法的步骤如下：

步骤一，收集dem数据：收集需要进行淤地坝调查区域和时段的分辨率大于等于30米的dem数据；

步骤二，收集遥感影像数据：收集对应区域和时段的分辨率大于等于10米的遥感影像数据；

步骤三，提取调查区域河网，确定集水面积小于50km²的河道位置和流向：基于gis平台，利用所收集的dem数据，以0.3km²为河道提取阈值，提取与调查区域河网，并制作0.3km²~50km²汇水面积河道位置栅格文件和流向栅格文件；

步骤四，制作区域栅格坡向文件：以30米dem为基础，制作30米栅格坡向文件；

步骤五，提取可能的淤地坝位置点区域：依据河道位置栅格文件进行坡向分析，在流向上，间距n个栅格以内的坡向如果差值的绝对值在180°附近，则标定为淤地坝预选点；

步骤六，选取淤地坝预选区域：以预选点为中心点，以半径rn个栅格区域为淤地坝预选位置区域；

步骤七，影像分析：针对淤地坝预选位置区域，提取相应的遥感影像数据，并进行遥感影像纹理特征分析；

步骤八，确定淤地坝位置：将淤地坝预选位置区域内河道上，与河道方向近似垂直的直线确定为淤地坝的位置；

步骤九：输出淤地坝位置与数量统计：输出淤地坝经纬度坐标，汇总淤地坝数量。

进一步的，所述的步骤五中n的确定方法为：淤地坝的坝底宽l，栅格分辨率d，则栅格数n=int（l/d）+1。

进一步的，所述的步骤六中半径rn的确定方法为：淤地坝长度ll，则rn=int（ll/d）+1。

本发明产生的有益效果是：本发明根据dem数据快速构建各个淤地坝可能出现位置的数字模型，再根据淤地坝在卫星图片的图像纹理特征与河道垂直的特点，利用图像识别技术辨认淤地坝可能出现的位置，最终建立各个淤地坝的位置及走向。结合dem数据能够快速进行大范围的辨认，并且图像识别技术能够快速的辨别大范围的地形图像，并根据淤地坝位置的特点缩小图像识别的范围，最终实现快速的大范围淤地坝位置调查。由于所使用的资料均为免费公开的资料，并且使用的人工费用极少，计算过程在普通计算机上即可实现，因此总体成本极低。

附图说明

下面结合附图和实施例对本发明作进一步说明。

图1是本发明的实施例所述方法的流程图。

具体实施方式

实施例一：

本实施例是一种大范围淤地坝位置确定方法。为实现完全自动化快捷的淤地坝调查，尽量减少人工干预，本实施例的基本思路是：首先获取低精度dem数据，进行河网分析，根据集水面积排除大量没有可能设立淤地坝的位置，只对预选位置和周围的地域进行分析，用十分简单的计算即可以得到淤地坝的可能位置（淤地坝预选位置）。利用流域数字河网确定研究区域内河道的位置，选取淤地坝预选位置的卫星图片，进行图像识别。利用现有技术对图像进行纹理分析。淤地坝在卫星照片上的纹理特征十分有特点，在卫星图像中，淤地坝的位置影像为与河道接近垂直或垂直的纹理。利用淤地坝的这个影像特点，找到在预选区域内的与河道垂直或接近垂直的直线纹理，即可以确定为淤地坝。

分辨率大于等于30米的dem数据是一种可以免费获得的dem数据。这种免费的dem数据的分辨率和精度具有一定的水平，能够达到确定淤地坝位置的目的。例如：2015年发布的astergdemv2数据，其空间分辨率达到30乘30米，高度误差小于20m，具有95%置信度。

通过免费的dem数据，可以自动获得研究区域的河网，以及河网的形态和周围环境，因此，通过dem数据的处理，就能够在计算机系统中分辨研究区域、可能出现淤地坝的控制范围，通过这些数据结合卫星照片能够确定是否为淤地坝，藉以确定淤地坝的位置和统计淤地坝的数量。

本实施例使用的所有资料均来自免费的途径，利用淤地坝在免费卫星图片上的图像特点结合免费dem数据，使用计算机大范围快速的得到淤地坝的分布。

本实施例所述方法的具体步骤如下（流程见图1所示）：

步骤一，收集dem数据：收集需要进行淤地坝调查区域和时段的分辨率大于等于30米的dem数据。这种分辨率大于等于30米的dem数据能够在网上直接下载，完全免费。

步骤二，收集遥感影像数据：收集对应区域和时段的分辨率大于等于10米的遥感影像数据。这种分辨率大于等于10米的遥感影像数据同样可以在网上下载，完全免费。

步骤三，提取调查区域河网，确定集水面积小于50km²的河道位置和流向：基于gis平台，利用所收集的dem数据，以0.3km²为河道提取阈值，提取与调查区域河网，并制作0.3km²~50km²汇水面积河道位置栅格文件和流向栅格文件。集水面积超过50km²的河道位置应当是水库了，已经超出淤地坝的集水面积，因此不予考虑。而集水面积小于0.3km²则可能是原有的河道也不予考虑，因此阈值选取为0.3km²~50km²，并据此制作汇水面积河道位置栅格文件和流向栅格文件。

在gis软件中使用水文分析工具，先对dem数据进行填洼处理，得到无洼地的修正dem，即将栅格高程低于周边8栅格的栅格高程抬升至周边栅格高程最低值，使得能够正确计算其流向；其后，对无洼地dem进行流向计算，即根据某栅格周围8个栅格的高程差，使得当前栅格流向指向坡降最陡的栅格；第三，根据栅格流向图层计算各栅格汇流累计数，即计算某栅格上游汇入的所有栅格数量；第四，根据汇流累计数图层，提取阈值为0.3km²的模拟河网，即创建1新图层，将其中汇流累计数大于334的栅格位置设置为1，其余设置为0；第五，根据汇流累计数图层，提取阈值为50km²的模拟河网，即创建1新图层，将其中汇流累计数大于55556的栅格位置设置为1，其余设置为0；第六，提取阈值为0.3km²~50km²的模拟河网（后面称之为预选区域河网），即使用提取的阈值为0.3km²的模拟河网减去阈值为50km²的模拟河网；第七，提取阈值为0.3km²~50km²的模拟河网流向（后面称之为预选区域河网流向），即对第二步计算的流向图层进行条件筛选，只保留预选区域河网值为1的区域的流向值。预选区域河网主要用于框定可能的淤地坝位置点，即该位置点必须位于河网栅格上。预选区域河网流向配合后面制取的坡向文件用于判别可能的淤地坝位置点，即顺着流向的n个栅格的坡向差绝对值在180°附近栅格即为可能点位。

步骤四，制作区域栅格坡向文件：以30米dem为基础，制作30米栅格坡向文件。

采用gis软件，以30米dem为基础，计算栅格坡向。即以某栅格为中心的3×3的窗口内的9个栅格高程为基础，根据一定的算法计算当前栅格的坡面朝向，其值为0-360°，0表示正北，90表示正东，180表示正南，270表示正西。该文件用于提取可能的淤地坝位置点区域。

步骤五，提取可能的淤地坝位置点区域：依据河道位置栅格文件进行坡向分析，在流向上，间距n个栅格以内的坡向如果差值的绝对值在180°附近，则标定为淤地坝预选点。n的选取一般可以是3个栅格，格栅选取得太少则不够精确，太多则增大计算量，不够经济。栅格数的确定方法可以是：淤地坝的坝底宽lm，栅格分辨率dm，则栅格数n=int（l/d）+1。

步骤六，选取淤地坝预选区域：以预选点为中心点，以半径rn个栅格区域为预选区域。以该位置为中心点，以半径rn个栅格区域为可能的淤地坝位置区域（淤地坝预选位置区域）。rn的具体确定方法可以是：淤地坝长度llm，栅格分辨率dm，则rn=int（ll/d）+1。

步骤七，影像分析：针对淤地坝预选位置区域，提取相应的遥感影像数据，并进行遥感影像纹理特征分析。对淤地坝预选位置的区域的卫星图像（遥感影像数据）进行分析，是否有直线纹理。这种直线纹理在卫星图像中能够十分容易地从众多自然线条中判别出来。

纹理分析可以使用纹理图像的方向场估计的方法，也可以使用其他方法进行。

步骤八，确定淤地坝位置：将淤地坝预选位置区域内河道上，与河道方向近似垂直的直线确定为淤地坝的位置。在对淤地坝预选位置区域找到直线纹理还不能完全确定是否为淤地坝，还要经过判别，即对直线纹理的方向与河道方向的关系进行判别，若与河道方向存在近似垂直的纹理（夹角在90°左右），则确认垂直纹理区域即为淤地坝。

步骤九：输出淤地坝位置与数量统计：输出淤地坝经纬度坐标，汇总淤地坝数量。对研究区域的核对进行上述研究，确定所有在河道上的淤地坝，之后总结淤地坝的数量。

实施例二：

本实施例是实施例一的改进，是实施例一关于第五步骤的细化。本实施例所述的步骤五中n的确定方法为：淤地坝的坝底宽l，栅格分辨率d，则栅格数n=int（l/d）+1。

本实施例所确定的是使用多少的栅格n进行坡面分析。这一参数与淤地坝的坝底宽和栅格的分辨率有关。淤地坝的坝底宽实际是淤地坝在水平面上的投影宽度，也就是淤地坝的最大面积，而格栅分辨率是水平面上格栅所覆盖的面积，公式的物理意义是，以最少的格栅数覆盖淤地坝宽度。

实施例三：

本实施例是实施例一的改进，是实施例一关于第五步骤的细化。本实施例所述的步骤六中半径rn的确定方法为：淤地坝长度ll，则rn=int（ll/d）+1。

本实施例中公式的物理意义是以最少的格栅数覆盖淤地坝的长度。这样，本实施例所述公式与实施例二中的公式结合，其物理意义就是用最少的格栅覆盖整个淤地坝。

最后应说明的是，以上仅用以说明本发明的技术方案而非限制，尽管参照较佳布置方案对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案（比如调查的具体方式、各种淤地坝的形式、步骤的先后顺序等）进行修改或者等同替换，而不脱离本发明技术方案的精神和范围。