CN114677401A - 一种基于极化雷达自身影像特征的水体提取方法及系统 - Google Patents
一种基于极化雷达自身影像特征的水体提取方法及系统 Download PDFInfo
- Publication number
- CN114677401A CN114677401A CN202210313593.1A CN202210313593A CN114677401A CN 114677401 A CN114677401 A CN 114677401A CN 202210313593 A CN202210313593 A CN 202210313593A CN 114677401 A CN114677401 A CN 114677401A
- Authority
- CN
- China
- Prior art keywords
- water body
- polarization
- polarization radar
- image
- threshold
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 125
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 111
- 238000000605 extraction Methods 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000002902 bimodal effect Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 description 6
- 230000005570 vertical transmission Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003809 water extraction Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000005571 horizontal transmission Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241000224466 Giardia Species 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于极化雷达自身影像特征的水体提取方法及系统,根据不同极化条件下的后向散射系数统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;计算不同极化雷达影像水体阈值之差的均值,作为阈值偏移量;针对双峰特征较弱的极化方式,在其后向散射系数上叠加阈值偏移量,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像;最后构建水体提取指数,并计算水体阈值,提取图幅范围内的水体像元。本发明不需要外部DEM即可消山体阴影的影响,实现影像图幅范围内的水体高精度提取。
Description
技术领域
本发明涉及洪涝灾害监测领域,更具体地,涉及一种基于极化雷达自身影像特征的水体提取方法及系统。
背景技术
星载合成孔径雷达系统(Synthetic Aperture Radar,SAR)成像时不依靠太阳光,使用波长较长的电磁波探测地物,受云、雨、雾等天气的影响较小,具有全天时、全天候获取地表信息的优势,能够弥补光学遥感的局限性,非常适用于对洪涝灾害的及时监测,是研究洪涝灾害的有效工具。
雷达遥感系统常用的极化电磁波有两种:水平极化(H)和垂直极化(V),传统单极化雷达只能发射和接收一种极化波,如HH(水平发射、水平接收)或VV(垂直发射、垂直接收)。在有限的波谱信息下很容易将水体与裸地、广场,以及停车场等地物混淆。而极化雷达可以同时发射和接收多种极化波,即HH,HV(水平发射、垂直接收),VH(垂直发射、水平接收)和VV,可以提供多个图像通道。因此,与传统单极化雷达相比,极化雷达可以获取更多的目标散射信息,实现更高精度的地物信息提取。
在水体提取方面,基于极化雷达影像的算法一般分为阈值分割、面向对象、机器学习三类。其中,面向对象考虑了地物的纹理、颜色特征,水体提取结果的连续性较强,不易受影像斑点噪声的影响;机器学习则因较强的特征提取能力,可以获得精度较高的提取结果。但是,成功应用上述方法的前提是高质量的训练数据和扎实的模型参数调整能力。因此,这两种方法的应用多见于学术论文中的小区域示范。
而阈值分割法因计算速度快、算法简易的优点,在诸如洪水淹没范围提取的应急响应中的具有较广泛的应用。例如,贾诗超等(2019)通过两个极化波段构建了SDWI指数来提取水体信息,李景刚等(2010)通过改进的OTSU算法实现水体提取。但是,该方法易受山体阴影等低回波地物的干扰,提取精度不高。目前普遍的改进方法是借助DEM(数字高程模型,Digital Elevation Model)数据,通过阈值分割生成山体掩膜,滤去初始的水体指数提取结果。这种做法一方面需要额外的数据输入,计算量较大;另一方面引入了第二套数据的阈值分割结果,存在新的误差。
发明内容
本发明的目的在于克服基于极化雷达影像提取水体时,阈值分割法易受山体阴影干扰的缺点,提出一种基于极化雷达自身影像特征的水体提取方法。该方法充分利用了极化雷达影像自身回波信号的统计特征,在不需要外部DEM输入的情况下,即可消山体阴影的影响,实现影像图幅范围内的水体高精度提取。
为解决上述技术问题,本发明的技术方案如下:
本发明提供一种基于极化雷达自身影像特征的水体提取方法,包括以下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ;
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
优选地,步骤S1具体包括:
S1.1:根据需要开展水体提取的时间、地点,下载相应的VV和VH极化雷达影像;
S1.2:对VV和VH极化雷达影像进行应用轨道文件校正处理;
S1.3:对VV和VH极化雷达影像进行热噪声去除处理;
S1.4:对VV和VH极化雷达影像进行斑点滤波处理;
S1.5:对VV和VH极化雷达影像进行地理编码处理;
S1.6:对VV和VH极化雷达影像进行分贝化处理,得到VV和VH两种极化方式的分贝化雷达后向散射系数,分别记为σVV和σVH。
优选地,步骤S2具体包括:
S2.1:对VV和VH极化雷达影像进行中值滤波,得到极化雷达影像后向散射系数σVV和σVH的统计直方图;
S2.2:按后向散射系数由小到大的顺序,取每N个区间(Bin)的频数最大值,设该最大值为因变量y,对应Bin所在的后向散射系数为自变量x,构建拟合函数f(x,y),以达到平滑原始直方图数据中的频数曲线,降低局部Bin之间的锯齿扰动的效果;
S2.3:根据拟合函数f(x,y)的一阶导数判断收敛性,找到频数的局部最小值y0,对应的Bin即为直方图双峰之间的谷底,进而得Bin所在的后向散射系数作为提取原始VV和VH极化雷达影像的水体阈值。
优选地,步骤S2.1中,中值滤波采用的滤波窗口大小为5×5像元。
优选地,步骤S2.2中,取N=5。
优选地,步骤S3中,计算VV和VH极化雷达影像的阈值之差,取其绝对值作为阈值偏移量tΔ。
优选地,步骤S4中,针对双峰特征较弱的VV极化雷达影像,每个像元均减去阈值偏移量tΔ,即σ’VV=σVV-tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠。
优选地,步骤S5中,比较每个像元的σ’VV和σVH大小,取每个像元σ’VV和σVH的最小值构成一景新影像,即全域最小值影像σMin。
优选地,步骤S6中,具体包括:
S6.1:构建水体提取指数PWI=lg(σMin×σMin);
S6.2:按步骤S2的方法提取水体阈值tWI;
S6.1:按PWI>tWI提取水体范围。
本发明还提供一种基于极化雷达自身影像特征的水体提取系统,包括存储器和处理器,所述存储器中包括基于极化雷达自身影像特征的水体提取方法程序,所述基于极化雷达自身影像特征的水体提取方法程序被所述处理器执行时实现如下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ;
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
与现有技术相比,本发明技术方案的有益效果是:本发明提供一种基于极化雷达自身影像特征的水体提取方法及系统,根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ;绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;最后构建水体提取指数,并计算水体阈值,提取图幅范围内的水体像元。本发明克服基于极化雷达影像提取水体时,阈值分割法易受山体阴影干扰的缺点,充分利用了极化雷达影像自身回波信号的统计特征,在不需要外部DEM输入的情况下,即可消山体阴影的影响,实现影像图幅范围内的水体高精度提取。
附图说明
图1为实施例1基于极化雷达自身影像特征的水体提取方法的流程图。
图2为哨兵1号GRD数据VV极化雷达影像图。
图3为哨兵1号GRD数据VH极化雷达影像图。
图4为σVV的统计直方图。
图5为σVH的统计直方图。
图6为修改后的VV和原始VV、VH的频数分布对比图。
图7为全域最小值影像σMin的示意图。
图8为水体提取指数PWI的直方图。
图9为通过本发明PWI方法和现有SDWI水体提取的效果对比图。
图10为图9中A、B、C、D区域的水体提取的效果对比图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
实施例1
如图1所示,本实施例提供基于极化雷达自身影像特征的水体提取方法,包括以下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ(单位:dB);
在具体实施过程中,步骤S1包括:
S1.1:如图2-3所示,根据需要开展水体提取的时间、地点,下载相应的VV(垂直发射、垂直接收)和VH(垂直发射、水平接收)极化雷达影像;本实施例以哨兵1号卫星地距多视影像(Ground Range Detected,GRD)数据为例,水体提取的地点为为珠江河口地区,水体提取的时间(成像日期)为2021年12月4日,极化方式为VV和VH。影像通过哨兵系列卫星科研数据中心网站(Sentinels Scientific Data Hub,https://scihub.copernicus.eu/dhus/#/home)下载获取,原始影像文件名为S1A_IW_GRDH_1SDV_20211204T103405_20211204T103430_040858_04D9E1_8879.zip。
S1.2:对VV和VH极化雷达影像进行应用轨道文件校正处理,替换元数据文件(.xml)中的哨兵1号卫星轨道状态信息;
S1.3:对VV和VH极化雷达影像进行热噪声去除处理,消除SAR卫星装置(例如,发射机、功率放大器、接收机)内部的热损耗,提高SAR影像的信噪比;
S1.4:对VV和VH极化雷达影像进行斑点滤波处理,使用目前较为常用的RefinedLee相干斑滤波器,去除影像相干斑噪声;
S1.5:对VV和VH极化雷达影像进行地理编码处理,使用Range-DopplerTerrainCorrection(距离多普勒法地形校正),赋予影像实际的坐标信息;
S1.6:对VV和VH极化雷达影像进行分贝化处理,将后向散射系数σ执行对数变换,即σ(dB)=10*lg(σ),得到VV和VH两种极化方式的分贝化雷达后向散射系数,分别记为σVV和σVH。
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
在具体实施过程中,步骤S2包括:
S2.1:对VV和VH极化雷达影像进行中值滤波,达到去除椒盐噪声、平滑影像的效果,得到极化雷达影像后向散射系数σVV和σVH的统计直方图,如图4-5所示;本实施例中,中值滤波的滤波窗口大小为5×5像元。
S2.2:按后向散射系数由小到大的顺序,取每5个区间(Bin)的频数最大值,设该最大值为因变量y,对应Bin所在的后向散射系数为自变量x,构建拟合函数f(x,y),以达到平滑原始直方图数据中的频数曲线,降低局部Bin之间的锯齿扰动的效果;
S2.3:根据拟合函数f(x,y)的一阶导数判断收敛性,找到频数的局部最小值y0,对应的Bin即为直方图双峰之间的谷底,进而得Bin所在的后向散射系数作为提取原始VV和VH极化雷达影像的水体阈值;本实施例中σVV的水体阈值tVV为-18.42dB,σVH的水体阈值tVH为-23.62dB。
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
在具体实施过程中,步骤S3中,计算VV和VH极化雷达影像的阈值之差,取其绝对值作为阈值偏移量tΔ。本实施例为tΔ=5.2。
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
在具体实施过程中,步骤S4中,针对双峰特征较弱的VV极化雷达影像,每个像元均减去阈值偏移量tΔ,即σ’VV=σVV-tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠。修改后的VV和原始VV、VH的频数分布对比如图6所示。
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
在具体实施过程中,步骤S5中,比较每个像元的σ’VV和σVH大小,取每个像元σ’VV和σVH的最小值构成一景新影像,即全域最小值影像σMin。全域最小值影像σMin如图7所示。
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
在具体实施过程中,步骤S6中,包括:
S6.1:构建水体提取指数PWI=lg(σMin×σMin),其直方分布如图8所示。
S6.2:按步骤S2所述的方法提取水体阈值tWI;本实施例中得到的水体阈值tW=6.384。
S6.1:按PWI>tWI提取水体范围,即对于每个像元,其PWI>tWI则该像元属于水体范围。该图9-10所示。可以看出,本方法在保留主干河流、水库水面提取效果的同时,可有效的消除山体阴影的干扰,提高极化雷达影像的水体提取精度。
实施例2
本实施例提供一种基于极化雷达自身影像特征的水体提取系统,包括存储器和处理器,所述存储器中包括基于极化雷达自身影像特征的水体提取方法程序,所述基于极化雷达自身影像特征的水体提取方法程序被所述处理器执行时实现如下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ(单位:dB);
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
在本申请所提供的实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,如:多个单元或组件可以结合,或可以集成到另一个系统,或一些特征可以忽略,或不执行。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元的间接耦合或通信连接,可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的,作为单元显示的部件可以是、或也可以不是物理单元;既可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外,在本发明各实施例中的各功能单元可以全部集成在一个处理单元中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
或者,本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
Claims (10)
1.一种基于极化雷达自身影像特征的水体提取方法,其特征在于,包括以下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ;
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
2.根据基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S1具体包括:
S1.1:根据需要开展水体提取的时间、地点,下载相应的VV和VH极化雷达影像;
S1.2:对VV和VH极化雷达影像进行应用轨道文件校正处理;
S1.3:对VV和VH极化雷达影像进行热噪声去除处理;
S1.4:对VV和VH极化雷达影像进行斑点滤波处理;
S1.5:对VV和VH极化雷达影像进行地理编码处理;
S1.6:对VV和VH极化雷达影像进行分贝化处理,得到VV和VH两种极化方式的分贝化雷达后向散射系数,分别记为σVV和σVH。
3.根据权利要求2所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S2具体包括:
S2.1:对VV和VH极化雷达影像进行中值滤波,得到极化雷达影像后向散射系数σVV和σVH的统计直方图;
S2.2:按后向散射系数由小到大的顺序,取每N个Bin的频数最大值,设该最大值为因变量y,对应Bin所在的后向散射系数为自变量x,构建拟合函数f(x,y),以达到平滑原始直方图数据中的频数曲线,降低局部Bin之间的锯齿扰动的效果;
S2.3:根据拟合函数f(x,y)的一阶导数判断收敛性,找到频数的局部最小值y0,对应的Bin即为直方图双峰之间的谷底,进而得Bin所在的后向散射系数作为提取原始VV和VH极化雷达影像的水体阈值。
4.根据权利要求3所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S2.1中,中值滤波采用的滤波窗口大小为5×5像元。
5.根据权利要求3所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S2.2中,取N=5。
6.根据权利要求3所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S3中,计算VV和VH极化雷达影像的阈值之差,取其绝对值作为阈值偏移量tΔ。
7.根据权利要求6所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S4中,针对双峰特征较弱的VV极化雷达影像,每个像元均减去阈值偏移量tΔ,即σ’VV=σVV-tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠。
8.根据权利要求7所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S5中,比较每个像元的σ’VV和σVH大小,取每个像元σ’VV和σVH的最小值构成一景新影像,即全域最小值影像σMin。
9.根据权利要求8所述的基于极化雷达自身影像特征的水体提取方法,其特征在于,步骤S6中,具体包括:
S6.1:构建水体提取指数PWI=lg(σMin×σMin);
S6.2:按步骤S2的方法提取水体阈值tWI;
S6.1:按PWI>tWI提取水体范围。
10.一种基于极化雷达自身影像特征的水体提取系统,其特征在于,包括存储器和处理器,所述存储器中包括基于极化雷达自身影像特征的水体提取方法程序,所述基于极化雷达自身影像特征的水体提取方法程序被所述处理器执行时实现如下步骤:
S1:根据需要开展水体提取的时间、地点,选择相应的极化雷达影像并完成预处理,得到不同极化条件下的后向散射系数σ;
S2:绘制不同极化条件下的后向散射系数σ的统计直方图,根据统计直方图的双峰特征,依据预设的自动阈值提取算法,定位到双峰之间的谷底位置,提取局部最小值作为原极化雷达影像的水体阈值;
S3:计算不同极化雷达影像水体阈值之差的均值,取该均值的绝对值作为阈值偏移量tΔ;
S4:针对双峰特征较弱的极化方式,在其后向散射系数σ上叠加阈值偏移量tΔ,使不同极化雷达影像的后向散射系数直方图在局部最小值处重叠;
S5:比较每个像元在不同极化方式中的后向散射系数大小,取最小值构成一景新影像σMin;
S6:构建水体提取指数,并根据步骤S2的方法计算水体阈值,提取图幅范围内的水体像元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210313593.1A CN114677401B (zh) | 2022-03-28 | 2022-03-28 | 一种基于极化雷达自身影像特征的水体提取方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210313593.1A CN114677401B (zh) | 2022-03-28 | 2022-03-28 | 一种基于极化雷达自身影像特征的水体提取方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114677401A true CN114677401A (zh) | 2022-06-28 |
CN114677401B CN114677401B (zh) | 2022-10-21 |
Family
ID=82075269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210313593.1A Active CN114677401B (zh) | 2022-03-28 | 2022-03-28 | 一种基于极化雷达自身影像特征的水体提取方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114677401B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998664A (zh) * | 2012-10-29 | 2013-03-27 | 中国科学院对地观测与数字地球科学中心 | 一种基于合成孔径雷达的水华识别方法及装置 |
CN103440489A (zh) * | 2013-09-16 | 2013-12-11 | 南京大学 | 基于像素级sar影像时间序列相似性分析的水体提取方法 |
EP2816529A2 (en) * | 2013-12-16 | 2014-12-24 | Institute of Electronics, Chinese Academy of Sciences | Automatic water area segmentation method and device for SAR image of complex terrain |
US20160125586A1 (en) * | 2014-11-03 | 2016-05-05 | The Johns Hopkins University | Ice analysis based on active and passive radar images |
CN105608414A (zh) * | 2015-12-11 | 2016-05-25 | 国网四川省电力公司电力应急中心 | 一种地表含水量分布提取方法 |
CN106600607A (zh) * | 2016-11-22 | 2017-04-26 | 武汉大学 | 一种基于水平集分割极化sar影像的水体精确提取方法 |
JP2017166949A (ja) * | 2016-03-16 | 2017-09-21 | 国立研究開発法人防災科学技術研究所 | 地形変動判定方法及び地形変動判定システム |
CN107329139A (zh) * | 2017-07-26 | 2017-11-07 | 邓少平 | 一种双时相双极化干涉sar影像水域提取方法 |
CN111259876A (zh) * | 2020-05-06 | 2020-06-09 | 中国水利水电科学研究院 | 一种基于陆表水体产品的雷达数据水体信息提取方法及系统 |
CN112102288A (zh) * | 2020-09-15 | 2020-12-18 | 北京百度网讯科技有限公司 | 水体识别和水体变化检测方法、装置、设备及介质 |
US20210149929A1 (en) * | 2019-11-20 | 2021-05-20 | University Of Connecticut | Systems and methods to generate high resolution flood maps in near real time |
CN112967308A (zh) * | 2021-02-26 | 2021-06-15 | 湖南南方水利水电勘测设计院有限公司 | 一种双极化sar影像的水陆边界提取方法及系统 |
CN113705441A (zh) * | 2021-08-27 | 2021-11-26 | 中国科学院南京地理与湖泊研究所 | 协同多光谱和sar影像的高时空分辨率地表水体提取方法 |
CN114119642A (zh) * | 2021-10-18 | 2022-03-01 | 中国科学院空天信息创新研究院 | 洪涝灾害水体提取方法及装置、设备和存储介质 |
-
2022
- 2022-03-28 CN CN202210313593.1A patent/CN114677401B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998664A (zh) * | 2012-10-29 | 2013-03-27 | 中国科学院对地观测与数字地球科学中心 | 一种基于合成孔径雷达的水华识别方法及装置 |
CN103440489A (zh) * | 2013-09-16 | 2013-12-11 | 南京大学 | 基于像素级sar影像时间序列相似性分析的水体提取方法 |
EP2816529A2 (en) * | 2013-12-16 | 2014-12-24 | Institute of Electronics, Chinese Academy of Sciences | Automatic water area segmentation method and device for SAR image of complex terrain |
US20160125586A1 (en) * | 2014-11-03 | 2016-05-05 | The Johns Hopkins University | Ice analysis based on active and passive radar images |
CN105608414A (zh) * | 2015-12-11 | 2016-05-25 | 国网四川省电力公司电力应急中心 | 一种地表含水量分布提取方法 |
JP2017166949A (ja) * | 2016-03-16 | 2017-09-21 | 国立研究開発法人防災科学技術研究所 | 地形変動判定方法及び地形変動判定システム |
CN106600607A (zh) * | 2016-11-22 | 2017-04-26 | 武汉大学 | 一种基于水平集分割极化sar影像的水体精确提取方法 |
CN107329139A (zh) * | 2017-07-26 | 2017-11-07 | 邓少平 | 一种双时相双极化干涉sar影像水域提取方法 |
US20210149929A1 (en) * | 2019-11-20 | 2021-05-20 | University Of Connecticut | Systems and methods to generate high resolution flood maps in near real time |
CN111259876A (zh) * | 2020-05-06 | 2020-06-09 | 中国水利水电科学研究院 | 一种基于陆表水体产品的雷达数据水体信息提取方法及系统 |
CN112102288A (zh) * | 2020-09-15 | 2020-12-18 | 北京百度网讯科技有限公司 | 水体识别和水体变化检测方法、装置、设备及介质 |
CN112967308A (zh) * | 2021-02-26 | 2021-06-15 | 湖南南方水利水电勘测设计院有限公司 | 一种双极化sar影像的水陆边界提取方法及系统 |
CN113705441A (zh) * | 2021-08-27 | 2021-11-26 | 中国科学院南京地理与湖泊研究所 | 协同多光谱和sar影像的高时空分辨率地表水体提取方法 |
CN114119642A (zh) * | 2021-10-18 | 2022-03-01 | 中国科学院空天信息创新研究院 | 洪涝灾害水体提取方法及装置、设备和存储介质 |
Non-Patent Citations (9)
Title |
---|
CHEN LEI等: "The Extraction of Water Body Information from TM Imagery Based on Water Quality Types", 《REMOTE SENSING FOR LAND & RESOURCES》 * |
LINAN BAO等: "Water Extraction in SAR Images Using Features Analysis and Dual-Threshold Graph Cut Model", 《REMOTE SENS》 * |
SEUNGHWAN HONG等: "Water Area Extraction Using RADARSAT SAR Imagery Combined with Landsat Imagery and Terrain Information", 《SENSORS》 * |
何颖清等: "基于高分三号雷达遥感影像的洪涝灾害监测——以郑州"7•20"特大暴雨灾害为例", 《中国水利学会2021学术年会论文集第二分册》 * |
周彬等: "基于高分三号卫星影像的水体信息提取技术研究", 《测绘与空间地理信息》 * |
孔亚男: "基于多核模型的多源遥感影像城市不透水面提取研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
崔倩等: "基于高分三号卫星数据的水体自动提取及应用", 《中国农业信息》 * |
崔舜铫等: "基于PALSAR雷达数据在南方水体信息提取研究", 《安徽农业科学》 * |
赵天杰等: "双频多极化SAR数据与多光谱数据融合的作物识别", 《地球信息科学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114677401B (zh) | 2022-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chierchia et al. | Multitemporal SAR image despeckling based on block-matching and collaborative filtering | |
IL276033B1 (en) | A method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing | |
CN110632594B (zh) | 一种长波长星载sar成像方法 | |
Liu et al. | PCA-based sea-ice image fusion of optical data by HIS transform and SAR data by wavelet transform | |
CN108230375A (zh) | 基于结构相似性快速鲁棒的可见光图像与sar图像配准方法 | |
CN108830808B (zh) | 基于相似线窗口均值补偿的星上红外图像条纹噪声去除方法 | |
CN105139396B (zh) | 一种全自动遥感影像云雾检测方法 | |
CN107689051A (zh) | 一种基于变化因子的多时相sar影像变化检测方法 | |
CN111580107A (zh) | 一种sar回波信号的射频干扰抑制方法、装置及成像方法 | |
CN114998365A (zh) | 一种基于极化干涉sar的地物分类方法 | |
CN110688923A (zh) | 一种基于哨兵1a sar数据的城市内涝风险区提取方法 | |
CN112882026A (zh) | 一种树高反演的方法、装置、计算机存储介质及终端 | |
CN114677401B (zh) | 一种基于极化雷达自身影像特征的水体提取方法及系统 | |
Kalaiyarasi et al. | A survey on: De-speckling methods Of SAR image | |
CN116626622A (zh) | 一种极化合成孔径雷达交叉定标方法 | |
CN106908791B (zh) | 基于全极化圆迹sar数据的输电线提取方法 | |
CN111123265B (zh) | 基于雷达时序图像的菠萝种植范围提取方法及装置、设备 | |
CN115147728A (zh) | 一种基于雷达数据和光学数据协同的塘坝快速识别方法 | |
CN113702933A (zh) | 一种合成孔径雷达辐射定标方法及装置 | |
Jardosh et al. | Detection of flood prone areas by flood mapping of SAR imagery | |
CN113469104A (zh) | 基于深度学习的雷达遥感影像地表水体变化检测方法和设备 | |
Lee et al. | Analysis of forest loss by Sentinel-1 SAR time series | |
Kaushik et al. | A Comparative study of pre-processing techniques of SAR images | |
Bralet et al. | Deep learning of radiometrical and geometrical SAR distorsions for image modality translations | |
CN117388812A (zh) | 一种sar影像辐射差异校正方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |