基于多源雷达遥感技术的河流径流监测反演方法技术

一种基于多源雷达遥感技术的河流径流监测反演方法，包括以下步骤：步骤1，多源雷达遥感数据预处理：SAR图像预处理、雷达高度计数据预处理；步骤2，基于SAR图像提取河流宽度；步骤3，基于雷达高度计提取河流相对水深；步骤4，基于曼宁公式构建径流计算模型；步骤5，时间序列河流径流变化监测。本发明专利技术提供一种适用于在极端天气下的大型河流径流多源雷达遥感监测技术，从而获得精确的河流参数信息。

River runoff monitoring and inversion method based on multi-source radar remote sensing technology

【技术实现步骤摘要】
基于多源雷达遥感技术的河流径流监测反演方法
本专利技术涉及合成孔径雷达(SyntheticApertureRadar，简称SAR)图像运用、河流径流监测等
，特别设计一种基于多源雷达遥感技术的河流径流监测反演方法。
技术介绍
河流是人类重要的用水来源，与人类的生活、生产息息相关，同时也是地球水循环的重要组成部分，在区域经济发展和气候变化中发挥着至关重要的作用。近些年来，随着科技的进步和人类社会的发展，人们越来越重视对水资源的监控。目前，受全球气温上升的影响，洪涝、干旱等极端水文现象和雷电、暴雨等灾害性天气发生的可能性增加，导致水文站监测仪器发生损坏的风险也随之增加。由于雷达高度计与合成孔径雷达均具有全天时、全天候观测，具有一定的穿透能力且不受天气影响的特点，因此雷达高度计经常被用来提取河流水位，SAR图像也经常被用来提取水域面积，虽然这两种技术可以较好的提取水域参数，但是二者并没有结合起来使用。利用遥感技术监测河流径流的核心就是提高径流结果的精度。目前，利用遥感技术监测河流径流，学者已做出一些研究，但是绝大多数方法都是基于雷达高度计与光学遥感图像，并且大都在已知河底高程、河底宽度等参数已知条件下才能应用。这种方法受天气影响较大，特别是在极端环境下无法为水文站提供技术补充。另外，河底高程、河底宽度等这些参数均不易测量，大多数水文站无法获取这些数据。综上所述，如何在极端天气下对大型河流监测是一个亟待解决的问题。
技术实现思路
本专利技术的目的在于提供一种基于多源雷达遥感技术的河流径流监测反演方法，以在各种天气条件，尤其在极端天气，精确获得河流径流信息。为了实现本专利技术的目的，其提供如下技术方案：一种基于多源雷达遥感技术的河流径流监测反演方法，包括以下步骤：步骤1，多源雷达遥感数据预处理：SAR图像预处理、雷达高度计数据预处理；步骤2，基于SAR图像提取河流宽度；步骤3，基于雷达高度计提取河流相对水深；步骤4，基于曼宁公式构建径流计算模型；步骤5，时间序列河流径流变化监测。本专利技术提供一种适用于在极端天气下的大型河流径流多源雷达遥感监测技术，从而获得精确的河流参数信息。附图说明图1为本专利技术总体流程图；图2为中国黄河下游三处河段分布图；图3为图2中三处实例河段SAR图像预处理后的结果；图4为图2中为SAR图像水域分割后的结果；图5-图7为雷达高度计反演水位与实测水位对比结果；图8为河流经验参数K的拟合过程；图9-图11为本专利技术模型反演径流与实测径流对比结果。图12为列表一；图13为列表二。具体实施方式下面结合附图和具体实施方式对本专利技术做进一步详细说明。如图1所示，步骤101：多源雷达遥感数据预处理：SAR图像预处理、雷达高度计数据预处理。SAR图像预处理采用辐射定标(RadiometricCalibrate)对SAR图像进行辐射校正，再使用RefinedLee滤波减少图像中的噪声斑点，然后采用Range-Doppler校正算子完成图像的几何校正，最后对图像进行图像拉伸，提高河流与背景的对比度。雷达高度计数据预处理则首先根据数据中的经纬度信息查询对应目标区域的星下点坐标，然后将该点数据提取出来。步骤102：基于SAR图像提取河流宽度。首先，使用混合模糊图像分割算法提取河流水面，统计水面的像素点个数N。接着，根据以下公式(1)计算图像长、宽对应的实际距离Llength、Lwidth；式中，(lat1，lng1)为起点坐标，(lat2，lng2)为终点坐标，R的取值为6378137米，即WGS84坐标系下椭球的半径长度。再根据以下公式(2)利用图像的长(Lm)、宽(Ln)和其对应的实际距离Llength、Lwidth计算单位像素代表的实际面积Sreal_pixel。然后，根据公式(3)将单位像素实际面积Sreal_pixel乘以水面像素点个数N计算水面面积Sriver。Sriver＝Sreal_pixel·N(3)最后，将目标研究区域的河段长度Lriver代入以下公式(4)，计算河流宽度Wriver。式中的Lriver由GoogleEarth的测距工具测得。步骤103：基于雷达高度计提取河流相对水深。具体地，首先根据以下公式(5)将数据中卫星与地球参考椭球面的距离Altitude减去卫星与目标水面距离Range以及校正量RangeCorrected得到水面高程Height；Height＝Altitude-Range-RangeCorrected(5)其中，校正量RangeCorrected根据以下公式(6)计算，具体地，校正量RangeCorrected包括大气校正、地球物理校正，大气校正又包括干对流层修正量Rangedt、湿对流层修正量Rangewt、电离层修正量Rangei，地球物理校正又包括固体潮修正量Rangeset、极潮修正量Rangept，以及水面椭球高GeoEGM96。然后利用对应时间的历史水位数据对Height进行水准面校正，得到以中国黄海水准面为基准的河流水位；最后将近几年的历史最低水位作为相对河底高程，令河流水位减去相对河底高程，计算相对水深。步骤104：基于曼宁公式构建径流计算模型。曼宁公式已被广泛的应用于河流径流计算方面的研究。目前常见的河流断面主要有三种，分别矩形断面、梯形断面和弧形断面，而对于大型河流而言，河流的宽度W＞＞深度D，因此上述三种断面的径流均可由公式(7)计算。其中，S是水面坡度，n为曼宁粗糙系数。第一步，获取上述计算的河宽与相对水深。第二步，由于水面坡度、曼宁粗糙系数等参数在实际测量中不并容易获得，若拥有这些参数的情况下，可以直接代入至公式(7)计算径流；若没有这些参数的情况下，如公式(8)所示，这里使用河流经验参数K表示水面坡度和曼宁粗糙系数：K由实测的历史径流数据拟合得到，将相对河底高程对应的径流作为基础径流Q0，此时径流计算模型如公式(9)所示，DR为相对河流深度：最后，根据公式(9)即可以计算河流径流Q。步骤105，时间序列河流参数变化监测。具体的，将时间序列下本方法反演的河流径流结果绘制成径流过程线，获得不同时间的河流径流及其随时间的变化量，并给出分析结果。为了验证本专利技术，专利技术人针对中国黄河的部分河段实施了本专利技术。图2为中国黄河下游花园口A、夹河滩B、高村水文站C附近的河段，黄河下游是中国华北地区重要的水资源，但是下游的河床已高出大堤背河地面3-5m，比两岸平原高出更多，上述三处水文站是治理黄河关键水利工程小浪底水库后的前三处水文站，对这三处水文站的径流监测将关系到整个下游流域乃至豫、鲁、皖、苏四省的安危；图3为上述三处水文站(分别对应子图a、b、c，下同)附近河段预处理后的SAR图像(2018/08/24)；图4为SAR图像水域分割后的结果(2018/08/24)；图5-图7为雷达高度计提取河流水位与实测河流水位的对比结果(2018/01-2018/12)；图8为河流经验参数K的拟合过程；图9-图11为时间序列下本专利技术模型反演径流与实测径流的对比结果(2018/01-2018/12)；图12为一列表，给出了多时序下本专利技术模型反演的径流结果精度评价，包含决定系数、均方根误差和相对均方根误差等，该结果已经满足了传统径流监测精度要求。图13为另一列表，给出了某四个本文档来自技高网...

【技术保护点】
1.一种基于多源雷达遥感技术的河流径流监测反演方法，其特征在于，包括以下步骤：步骤1，多源雷达遥感数据预处理：SAR图像预处理、雷达高度计数据预处理；步骤2，基于SAR图像提取河流宽度；步骤3，基于雷达高度计提取河流相对水深；步骤4，基于曼宁公式构建径流计算模型；步骤5，时间序列河流径流变化监测。

【技术特征摘要】
1.一种基于多源雷达遥感技术的河流径流监测反演方法，其特征在于，包括以下步骤：步骤1，多源雷达遥感数据预处理：SAR图像预处理、雷达高度计数据预处理；步骤2，基于SAR图像提取河流宽度；步骤3，基于雷达高度计提取河流相对水深；步骤4，基于曼宁公式构建径流计算模型；步骤5，时间序列河流径流变化监测。2.根据权利要求1所述的方法，其特征在于，在步骤1中，SAR图像预处理具体包括：采用辐射定标对SAR图像进行辐射校正，使用RefinedLee滤波减少图像中的噪声斑点，采用Range-Doppler校正算子完成图像的几何校正，对图像进行图像拉伸，提高河流与背景的对比度；雷达高度计数据预处理具体包括：根据数据中的经纬度信息查询对应目标区域的星下点坐标，将该点数据提取出来。3.根据权利要求1所述的方法，其特征在于，步骤2包括：使用混合模糊图像分割算法提取河流水面，统计水面的像素点个数N；根据以下公式(1)计算图像长、宽对应的实际距离Llength、Lwidth；式中，(lat1，lng1)为起点坐标，(lat2，lng2)为终点坐标，R为6378137米；根据以下公式(2)利用图像的长(Lm)、宽(Ln)和其对应的实际距离Llength、Lwidth计算单位像素代表的实际面积Sreal_pixel；根据以下公式(3)将单位像素实际面积Sreal_pixel乘以水面像素点个数N计算水面面积Sriver：Sriver＝Sreal_pixel·N(3)将目标研究区域河段长度Lriver代入以下公式(4)，计算河流...

【专利技术属性】
技术研发人员：毋琳，郭拯危，闵林，李宁，赵建辉，王宁，
申请(专利权)人：河南大学，
类型：发明
国别省市：河南,41