一种基于卫星遥感和DEM数据的水库蓄水量计量方法技术

本发明专利技术公开了一种基于卫星遥感和DEM数据的水库蓄水量计量方法，属于蓄水量计量域。本发明专利技术包括步骤：S1建立基于卫星遥感影像和水文站监测数据的水库水位‑水面面积关系曲线；S2基于数字高程模型进行DEM插补，获取无遥感影像资料的水库水位相应的水面面积；S3结合DEM插补建立基于DEM的水库水位‑水面面积关系曲线；S4利用相关关系回归计算方程对S1获得的水库水位‑水面面积关系曲线以及S3获得的水库水位‑水面面积关系曲线合并，得到完整的水库水位‑水面面积关系曲线；S5由S4得到的完整的水库水位‑水面面积关系曲线进一步推求水库水位‑水面面积‑水库库容关系曲线。本发明专利技术方法显著提高了水库蓄水状态监测精度。

【技术实现步骤摘要】
一种基于卫星遥感和DEM数据的水库蓄水量计量方法
本专利技术涉及水利工程
，特别涉及一种基于卫星遥感和DEM数据的水库蓄水量计量方法。
技术介绍
地表水体蓄水量监测是水资源评价中的一个重要环节，蓄水量变化也是水资源管理的重要参考指标。传统地表水体蓄水量的估算方法主要依赖于地面水文站点的监测数据，通过水库特征曲线查算水库蓄水量。而实际上，由于水文站点分布的有限性的不均匀性，区域蓄水量的精确计算存在一定的困难。水库在建成后多年运行造成库区的淤积，改变了原有的水位-库容关系，或者水库在建成时受技术条件等方面的限制，并不能获得真实、准确的水位-库容关系曲线，因此，水库建成或运行多年后，有必要对水库的水位-库容曲线进行复核，掌握水库准确的水位-库容关系曲线，在水库、大坝安全的前提下，最大限度地发挥水库地效益。
技术实现思路
为了弥补现有技术的不足，解决现有技术不能获得真实准确的水位-库容关系的问题，本专利技术提供了一种基于卫星遥感和DEM数据的水库蓄水量计量方法。本专利技术的技术方案为：一种基于卫星遥感和DEM数据的水库蓄水量计量方法，包括步骤：S1建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线；S2基于数字高程模型进行DEM插补，获取无遥感影像资料的水库水位相应的水面面积；S3结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；S4利用相关关系回归计算方程对S1获得的水库水位-水面面积关系曲线以及S3获得的水库水位-水面面积关系曲线合并，得到完整的水库水位-水面面积关系曲线；S5由S4得到的完整的水库水位-水面面积关系曲线进一步推求水库水位-水面面积-水库库容关系曲线。作为一种实施方式，S1中，获取水库不同时间下的卫星遥感影像水面面积信息，并获取对应的水文站监测水位数据，利用线性回归建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线(Ⅰ)；S1(h)＝α1h2+β1h+γ1(Ⅰ)其中，S1(h)为基于卫星遥感影像面积S1关于水位h的函数，h表示水位，α1、β1、γ1为遥感影像水库水位-水面面积关系曲线修正系数。作为一种实施方式，S3中，获取水库的DEM数据信息和对应的监测水位信息，并结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；S2(h)＝α2h2+β2h+γ2(Ⅱ)其中，S2(h)为基于DEM数据面积S2关于水位h的函数，h表示水位，α2、β2、γ2为DEM水库水位-水面面积关系曲线修正系数。作为一种实施方式，S4中，利用相关关系回归计算方程对曲线(Ⅰ)和曲线(Ⅱ)进行合并，得到完整的水库水位-水面面积关系曲线(Ⅲ)；S(h)＝Ah2+Bh+C(Ⅲ)B＝θβ1+(1-θ)β20<θ<1C＝εγ1+(1-ε)γ20<ε<1。作为一种实施方式，S5中，根据S4获得的水库水位-水面面积关系曲线(Ⅲ)建立水库水位-水面面积查算表，按0.01m的间隔计算水库水位相应的面积，以整个水域由若干紧密排列的四棱柱构成，得到相应水位差下的库容：Vi＝hi[kSi+(1-k)Si+1]0<k<1(Ⅳ)式中，Vi表示第i个四棱柱的体积，hi为第i个四棱柱的高程，Si表示第i个四棱柱对应的水面面积，Si+1表示第i+1个四棱柱对应的水面面积，k为修正系数。作为一种实施方式，水库蓄水量，即，累计库容按照式(Ⅴ)计算；式中，V表示水库蓄水量，Vi表示水库第i部分的库容，V0表示死库容。作为一种实施方式，S2中，通过四个临近点的中心点进行数字高程模型DEM的插补，定义插求点为Q，插求点Q的四个临近点为(i,j)、(i,j+1)、(i+1,j)、(i+1,j+1)，(i,j)、(i,j+1)的中点为Q0，则Q0的高程值Z(Q0)＝0.5[Z(i,j)+Z(i,j+1)]；(i+1,j)、(i+1,j+1)的中点为Q1，则Q1的高程值Z(Q1)＝0.5[Z(i+1,j)+Z(i+1,j+1)]，通过式(Ⅵ)求得插求点Q的高程Z(Q)；Z(Q)＝Z(Q0)(1-t)+Z(Q1)t(0≤t≤1)(Ⅵ)，t为修正系数。作为一种实施方式，S1中，在提取水库卫星遥感影像水面面积之前，先对遥感影像进行辐射校正和几何校正，然后利用绿光波段和近红外波段组合构建归一化差异水体指数NDWI，设置阈值提取水体信息，统计水体像元个数，根据栅格面积得到水面面积；归一化差异水体指数NDWI计算方法如式(Ⅶ)所示：NDWI＝(ρGREEN-ρNIR)/(ρGREEN+ρNIR)(VII)。作为一种实施方式，所述辐射校正包括辐射定标和大气校正，用于消除设备误差和大气干扰，得到真实的反射率数据；所述几何校正采用多项式函数校正方法，选取控制点纠正影像的几何畸变。本专利技术的有益效果为：将遥感数据提取的水体参数作为综合数据的一部分，建立遥感数据提取水体的水面面积与地面监测水位数据的同化模型，更新水库特征曲线，显著提高水库蓄水状态监测精度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术基于卫星遥感和DEM数据的水库蓄水量计量方法的流程图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。如图1所示，一种基于卫星遥感和DEM数据的水库蓄水量计量方法，包括步骤：S1对遥感影像进行预处理，提取遥感影像提取水体面积，同时提取对应的水文监测数据的水位信息，建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线；S2基于数字高程模型进行DEM插补，获取无遥感影像资料的水库水位相应的水面面积；S3提取DEM相关数据，并结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；S4利用相关关系回归计算方程对S1获得的水库水位-水面面积关系曲线以及S3获得的水库水位-水面面积关系曲线合并，得到完整的水库水位-水面面积关系曲线；S5由S4得到的完整的水库水位-水面面积关系曲线进一步推求水库水位-水面面积-水库库容关系曲线。实施例1本实施例所用数据源包括：多光谱遥感数据：为了满足时间、尺度和成本方面的要求，采用重访周期较短的中分辨率卫星影像数据，如LandsatTM数据、HJ数据，以及重访周期较长、中高分辨率的高分系列卫星数据，分本文档来自技高网...

【技术保护点】
1.一种基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于，包括步骤：/nS1建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线；/nS2基于数字高程模型进行DEM插补，获取无遥感影像资料的水库水位相应的水面面积；/nS3结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；/nS4利用相关关系回归计算方程对S1获得的水库水位-水面面积关系曲线以及S3获得的水库水位-水面面积关系曲线合并，得到完整的水库水位-水面面积关系曲线；/nS5由S4得到的完整的水库水位-水面面积关系曲线进一步推求水库水位-水面面积-水库库容关系曲线。/n

【技术特征摘要】
1.一种基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于，包括步骤：
S1建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线；
S2基于数字高程模型进行DEM插补，获取无遥感影像资料的水库水位相应的水面面积；
S3结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；
S4利用相关关系回归计算方程对S1获得的水库水位-水面面积关系曲线以及S3获得的水库水位-水面面积关系曲线合并，得到完整的水库水位-水面面积关系曲线；
S5由S4得到的完整的水库水位-水面面积关系曲线进一步推求水库水位-水面面积-水库库容关系曲线。


2.如权利要求1所述基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于：S1中，获取水库不同时间下的卫星遥感影像水面面积信息，并获取对应的水文站监测水位数据，利用线性回归建立基于卫星遥感影像和水文站监测数据的水库水位-水面面积关系曲线(Ⅰ)；
S1(h)＝α1h2+β1h+γ1(Ⅰ)
其中，S1(h)为基于卫星遥感影像面积S1关于水位h的函数，h表示水位，α1、β1、γ1为遥感影像水库水位-水面面积关系曲线修正系数。


3.如权利要求2所述基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于：S3中，获取水库的DEM数据信息和对应的监测水位信息，并结合DEM插补建立基于DEM的水库水位-水面面积关系曲线；
S2(h)＝α2h2+β2h+γ2(Ⅱ)
其中，S2(h)为基于DEM数据面积S2关于水位h的函数，h表示水位，α2、β2、γ2为DEM水库水位-水面面积关系曲线修正系数。


4.如权利要求3所述基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于：S4中，利用相关关系回归计算方程对曲线(Ⅰ)和曲线(Ⅱ)进行合并，得到完整的水库水位-水面面积关系曲线(Ⅲ)；
S(h)＝Ah2+Bh+C(Ⅲ)



B＝θβ1+(1-θ)β20<θ<1
C＝εγ1+(1-ε)γ20<ε<1。


5.如权利要求4所述基于卫星遥感和DEM数据的水库蓄水量计量方法，其特征在于：S5中，根据S4获得的水库水位-水面面积关...

【专利技术属性】
技术研发人员：孙启玉，刘玉峰，邓青，李红超，荆志铎，
申请(专利权)人：山东锋士信息技术有限公司，
类型：发明
国别省市：山东;37