基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法技术

本发明专利技术公开了一种基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，属于水文遥感技术领域。通过中值滤波去除Sentinel

【技术实现步骤摘要】
基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法


[0001]本专利技术涉及一种基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，尤其适用于遥感雷达数据的洪泛湿地结构水文连通性的提取，属于水文遥感


技术介绍

[0002]水文连通性是指水体在河流、湖泊以及洪泛湿地之间的流通和交换的能力，是维持湿地生态系统服务和功能的关键性指标之一。水文连通性影响湿地的淹没范围、淹水深度以及洲滩出露时间，进而塑造湿地植物分布格局及其时间动态。水文连通性的概念因研究对象和目标的差异，其定义和度量方法有所不同。不同观测技术和方法在量化水文结构连通性和功能连通性中的应用，为后续研究提供了重要支撑。水文连通性度量方法大致有图论方法、物理或化学等参数的运移、基于水动力模型和实地调查的地统计分析和基于遥感产品提取的连通性指标等。基于图论方法的结构连通性往往指示了一种静态的水文连通性，而后两种水文水动力要素的空间分布及其与气象等环境要素的耦合关系，但有效的野外观测和调查数据成本高昂。
[0003]现有的水文连通度的度量数据源主要有：站点实测数据，以及相对静态的土地利用、地形数据和遥感数据。站点数据源的水文连通性难以实现水文联通的空间量化。基于相对固定的DEM、土地利用数据计算，得到了一个相对静态结果。而洪泛湿地中，水文特征快速变化，短期内水位快速涨落，年际差异大，静态的水文连通度无法反应这种快速的变化。基于遥感数据提取水体，具有实时性强、成本低、效率高等优点。但是常用的光学影像虽然具有丰富的光谱信息，容易受夜间、云、雨等天气因素的影响，导致数据源不稳定。

技术实现思路

[0004]针对现有技术的不足之处，提供一种基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，先判定长期河道位置，并计算洪泛湿地范围内不同像元与长期河道连通时长作为量化水文连通性指标；在不依赖成本较高的水文要素及气象要素监测数据的情况下，实现洪泛湿地结构水文连通度的提取；不受天气影响，动态的洪泛湿地水文连通度；其计算量小，处理效率高，实时性强。
[0005]为达到上述目的，本专利技术基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，其特征在于步骤如下：
[0006]步骤1、收集洪泛湿地区域的Sentinel
‑
1雷达数据影像，对Sentinel
‑
1雷达数据影像进行预处理，去除Sentinel
‑
1雷达数据影像中影响分类的噪声点；
[0007]步骤2、根据每景Sentinel
‑
1雷达数据影像的灰度特性，将Sentinel
‑
1雷达数据影像分成背景和前景两部分，计算每景Sentinel
‑
1雷达数据影像的灰度频率直方图，依据使类间方差最大化原则，得到水体、非水体分类的最优阈值K，利用最优阈值K对Sentinel
‑
1雷达数据影像进行逐个像元分类，得到分类后的水体/非水体数据；将水体赋值为1，非水体赋值为0；
[0008]步骤3、将水体/非水体数据进行空间叠加分析得到预设时间段内水体淹没概率，将水体淹没概率小于95％的像元组成的区域删除，将水体淹没概率大于95％的像元组成的区域标注为长期河道水体，简称长期河道；
[0009]步骤4、为去除洪泛湿地中短期淹没的水体，将水体/非水体数据转为矢量数据，判断矢量数据影像内每个矢量水体是否与步骤3判断为长期河道水体的区域相交，如部分相交则判断为与长期河道连通水体；不相交则为阻隔水体并移除；得到与长期河道连通的水体空间分布数据，通过空间叠加计算后，最终得到预设时间范围内与长期河道连通水体时长，作为量化洪泛湿地水文结构连通性的指标。
[0010]进一步，通过待计算像素及其周围3x3邻域内像素点的像素值进行排序，然后将位于中间位置的像素值作为当前像素点的像素值，让周围的像素值接近真实值从而消除孤立的噪声点。
[0011]进一步，获取水体数据时，当背景和前景之间的类间方差越大，则构成图像的两部分的差别越大，使类间方差最大的分割意味着错分概率最小；
[0012]类间方差最大的阈值最优选择公式为：
[0013][0014]其中，k
*
是区分水体和非水体的最佳阈值，ω0(k)和ω1(k)是被阈值k分为水体和非水体的概率，μ0(k)和μ1(k)是水体和非水体在阈值处的平均值，δ2是总方差。
[0015]有益效果：本专利技术的目的是基于高时空分辨率遥感雷达数据，融合地学空间分析方法，提出一种洪泛湿地水文结构连通性时空动态监测的方法。洪泛湿地淹没范围具有高度动态性，受制于数据和方法的限制难以于精细尺度上进行时空动态变化监测。本专利技术突破了常规数据频率低、覆盖范围窄的局限性，实现洪泛湿地复杂水文特征的时间和空间量化。
[0016]本专利技术的优点：1)本专利技术基于免费开源数据计算，提高了方法可应用范围，避免了传统实地监测调查的人力和经济成本；2)基于雷达数据提取的水文连通度指数可以得到高时空分辨率的水文连通度数据，有效反应洪泛湿地中水文的高度动态特性。
附图说明
[0017]图1为本专利技术基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法的流程示意图；
[0018]图2为本专利技术的实施例中利用中值滤波处理雷达数据影像的噪声示意图；
[0019]图3为本专利技术实施例中Sentinel
‑
1雷达数据影像极化数据频率分布直方图；
[0020]图4为本专利技术实施中Sentinel
‑
1雷达数据分类示意图，图中A.Sentinel
‑
1雷达数据影像极化数据B.水体、非水体分类结果；
[0021]图5为本专利技术实施中长期河道水体示意图；
[0022]图6为本专利技术实施中判断水体与长期河道是否连通概念图；
[0023]图7为本专利技术实施中水文连通性提取示意图。
具体实施方式
[0024]以下结合附图对本专利技术的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本专利技术，并不用于限定本专利技术。
[0025]如图1所示，本专利技术的基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，其步骤如下：
[0026]步骤1：去除Sentinel
‑
1雷达数据影像的噪声
[0027]Sentinel
‑
1雷达数据影像通常具有斑点噪声，会降低影像质量，影响分类精度。水体中有细小的噪声，在分类过程中被误分为陆地。采用中值滤波处理，通过待计算像素及其周围3x3邻域内像素点的像素值进行排序，然后将位于中间位置的像素值作为当前像素点的像素值，让周围的像素值接近真实值从而消除孤立的噪声点。如图，待计算像素为噪点，值为2。将待计算像素点周围3x3邻域内像素点按从小到大排列为{2,16,25,27,28,36,44,45,80}，中位数为28，将28赋值给待计算像素，具体见图2所示。
[0028]步骤2本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于遥感雷达数据的洪泛湿地结构水文连通性的提取方法，其特征在于步骤如下：步骤1、收集洪泛湿地区域的Sentinel
‑
1雷达数据影像，对Sentinel
‑
1雷达数据影像进行预处理，去除Sentinel
‑
1雷达数据影像中影响分类的噪声点；步骤2、根据每景Sentinel
‑
1雷达数据影像的灰度特性，将Sentinel
‑
1雷达数据影像分成背景和前景两部分，计算每景Sentinel
‑
1雷达数据影像的灰度频率直方图，依据使类间方差最大化原则，得到水体、非水体分类的最优阈值K，利用最优阈值K对Sentinel
‑
1雷达数据影像进行逐个像元分类，得到分类后的水体/非水体数据；将水体赋值为1，非水体赋值为0；步骤3、将水体/非水体数据进行空间叠加分析得到预设时间段内水体淹没概率，将水体淹没概率小于95％的像元组成的区域删除，将水体淹没概率大于95％的像元组成的区域标注为长期河道水体，即长期河道；步骤4、为去除洪泛湿地中...

【专利技术属性】
技术研发人员：滕佳昆，陈继龙，王冉，杨慧，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：