ALOS PALSAR散射总功率的面向对象林火区提取
1.
2.
Extraction of the forest fire region based on the span of ALOS PALSAR by object-oriented analysis
1.
2.
责任编辑: 陈理
收稿日期: 2018-11-30 修回日期: 2019-03-31 网络出版日期: 2019-12-15
基金资助: |
|
Received: 2018-11-30 Revised: 2019-03-31 Online: 2019-12-15
作者简介 About authors
姜德才(1990-),男,工程师,主要从事机器学习在遥感影像分类、信息提取方面的应用研究。Email:jiangdc_rs@126.com。 。
利用合成孔径雷达(synthetic aperture Radar,SAR)开展林地过火区提取主要基于单个像元开展,而以像元集合为分析单元的面向对象技术在该领域的应用研究还较少。本研究基于ALOS PALSAR影像的散射总功率数据,采用分型网络演化的多尺度分割算法提取2009年发生在美国Alaska中东部的林火灾后的过火区,并将林火区提取结果与获取的燃烧强度辅助数据进行对比,验证算法应用的有效性。研究表明,相比之前基于像元的阈值分类方法,基于本文方法火后单一时相和林火前后2时相SAR散射总功率的林火区提取精度分别提高了12.7%和15.8%。有效证明了面向对象技术能够有效应用于SAR影像的信息提取,SAR卫星遥感在森林火灾监测方面具有较好的应用潜力。
关键词:
At present, the extraction accuracy of the forest fire area by synthetic aperture Radar (SAR) is mainly limited to the analysis of single pixel. However, the application study of object-oriented technology based on pixel sets as the analysis unit is less in dealing with SAR images. In this paper, a multi-scale segmentation algorithm based on fractal net evolution approach (FNEA) was applied to the span of ALOS PALSAR images. Through the application research, the forest fire region, which happened in 2009 and was located in the Middle East of Alaska, USA, was extracted. The application validation of the algorithm was verified by comparing the experiment results with the auxiliary data of monitoring trends in burn severity (MTBS) data. The experiment results show that the classification accuracies of one-static span and two-static spans based on object-oriented analysis are improved by 12.7% and 15.8% respectively, compared with precious research. The researches show that object-oriented technology can be effectively applied to the information extraction form SAR image, and SAR technology has potential application in forest fire monitoring.
Keywords:
本文引用格式
姜德才, 李文吉, 李敬敏, 白罩峰.
JIANG Decai, LI Wenji, LI Jingmin, BAI Zhaofeng.
0 引言
林火对森林生态系统具有极大危害。林火不仅破坏植被及环境,如干扰植物群落的更替和野生动物的迁徙、改变土壤结构、降低微生物含量、造成水土流失、引起空气污染等[1],而且还威胁人类生命及财产安全。森林过火迹地的位置、面积和烈度的高精度、快速提取可为相关部门针对森林火中救治、火后灾情评估及生态恢复的决策、部署提供科学依据。
基于晴天无云的中高空间分辨率光学卫星遥感的林火区提取研究已广泛开展。该类方法主要是利用综合火烧指数(composite burn index,CBI)、归一化差值植被指数(normalized difference vegetation index,NDVI)、归一化火烧指数(normalized burn ration,NBR)和差分归一化火烧指数(differential normalized burn ration,DNBR)等光学影像遥感指数[2,3],通过人工目视解译或应用分类算法获取水体、裸地、植被以及燃烧区等地类,从而提取林火迹地; 对提取的林火迹地进一步细分还可获得燃烧烈度信息[4]。孙桂芬等[5]基于晴天无云的Landsat8影像短波红外波段提取林火迹地的精度达到了87%,而GF-1影像可见光和近红外波段提取林火迹地的精度则达到了91%。
1 研究方法
1.1 研究区及数据源
研究区位于美国Alaska中东部的Birch Creek机场西南角35 km处,火区中心点坐标为66°02'42.0″N,146°35'16.8″W,林火发生于2009年7月18日—8月5日。
图1
表1 ALOS PALSAR影像参数
Tab.1
数据类型 | 波长/cm | 像元大小(距离向× 方位向)/(m×m) | 升降轨 | 成像时间 | 像元个数(距离向× 方位向)/个 | 中心入射角/(°) |
---|---|---|---|---|---|---|
ALOS PALSAR | 23.6 L波段 | 9.37×3.58 | 降轨 | 2009/08/03 | 1 248×18 432 | 23.935 |
2009/06/18 | 1 248×18 432 | 23.941 |
1.2 极化散射总功率
全极化SAR影像数据以复数矩阵形式测量、记录4种极化状态HH,HV,VH和VV下地球分辨单元内所有分布式目标散射回波的相干叠加,构成了极化散射矩阵
通过极化协方差矩阵
式中
极化散射总功率,又称Span,由4种极化散射强度构成,包含地物丰富的纹理信息且能够有效地抑制相干噪声[14],计算公式为
1.3 FNEA方法
分类过程中,影像分割对象的合并需解决2个关键问题: ①定量表征相邻分割对象(像元)的异质性(或同质性); ②定量表征影像分割对象在合并前后的异质性(或同质性)变化。相邻分割对象间的异质性,即给定一个特定的特征空间,如果2个相邻分割对象在这个特征空间彼此接近,那么这2个对象相似。对于一个
式中f为对象特征,具体表现为光谱平均值、方差或纹理特征。特征空间距离可以通过计算所有分割对象在所有特征维的标准差
合并前后影像分割对象的异质性变化,即从单个像素对象开始,逐步合并到较大的对象上,合并对象的光谱异质性会增加。优化过程的基本目标是在每个合并过程中,最小化合并对象的异质性增加。通过影像分割对象合并到邻域对象,从而产生最小的异质性增加。因此,通过描述合并前后异质性的变化
式中:
1.4 精度评价
式中:
参照MTBS燃烧强度数据,为各样本像元赋属性值,即水体、燃烧区和未燃烧区。
1.5 技术流程
为了抑制SAR影像的相干斑噪声,在SAR影像预处理环节,对主辅影像分别进行多视处理: 距离向×方位向=1×6。多视后,主辅SAR影像的地面空间分辨率为: 距离向×方位向=23.1 m×20.7 m。高精度配准[7]的控制点误差均小于0.1个像元。对林火极化散射总功率取以10为底的对数,从而构建火后单一时相特征图
图2
2 结果与分析
提取结果如图3所示,图3(a)和图3(b)中,绿色区域为未燃烧区; 蓝色区域为水体; 红色区域为燃烧区。图3(a)为研究区火后单一时相
图3
表2 火后单一时相的混淆矩阵
Tab.2
样本像元分类 | 未燃烧区 | 水体 | 燃烧区 | 总和 |
---|---|---|---|---|
未燃烧区 | 378 | 17 | 40 | 435 |
水体 | 4 | 23 | 1 | 28 |
燃烧区 | 31 | 0 | 106 | 137 |
总和 | 413 | 40 | 147 | |
总体精度: 84.5% Kappa系数: 0.649 |
表3 林火前后2时相的混淆矩阵
Tab.3
样本像元分类 | 未燃烧区 | 水体 | 燃烧区 | 总和 |
---|---|---|---|---|
未燃烧区 | 385 | 17 | 23 | 425 |
水体 | 1 | 20 | 0 | 21 |
燃烧区 | 27 | 3 | 124 | 154 |
总和 | 413 | 40 | 147 | |
总体精度: 88.2% Kappa系数: 0.735 |
从表2和表3可知,①研究区火后单一时相
3 结论
合成孔径雷达(synthetic aperture Radar,SAR)技术可作为光学卫星遥感监测林火的有效补充手段,能够穿透烟、云、雾、雨、雹等恶劣天气,可用于林火的火中监测。本文开展了以像元集合为分析单元的面向对象技术在SAR散射总功率的林火区提取研究。以L波段全极化ALOS PALSAR为数据源,经散射总功率的计算、分型网络演化的多尺度分割、阈值分类、精度评价等过程,开展了美国Alaska中东部一处林火灾后过火区的提取研究,主要得到以下结论:
1)基于火后单一时相、林火前后2时相的SAR散射总功率开展林火区的提取结果整体与MTBS燃烧强度都较为吻合,且后者提取结果更优。
2)对比研究区SAR散射总功率基于像元阈值分类的林火分类结果,基于多尺度分割的面向对象阈值分类结果精度有较大提高,研究区火后单一时相的分类精度提高了12.7%,林火前后2时相的分类精度提高了15.8%。
3)林火研究区分类精度基本达到了基于晴天无云的常规高分可见光和红外波段影像的燃烧区和未燃烧区的分类精度。因此,极化SAR技术可作为光学卫星遥感监测林火的有效补充手段,且能适用于林火灾中的应急监测。
但是,仍有一些问题有待进一步研究,如开展多参数特征图的阈值分类,以减少未燃烧区的误分类; 深度挖掘SAR数据的极化信息,并应用于森林火灾的过火区及燃烧烈度的提取研究等。
志谢: 感谢美国国家航空航天局网站(https: //search.earthdata.nasa.gov)提供的ALOS PALSAR数据,感谢美国地质调查局的MTBS(monitoring trends in burn severity)项目网站(https: //www.mtbs.gov/direct-download)提供的燃烧强度数据。
参考文献
呼中林区火烧迹地遥感提取及林火烈度的空间分析
[J].
Burned area extraction in Huzhong forests based on remote sensing and the spatial analysis of the burned severity
[J].
L-band SAR backscatter prospects for burn severity estimation in Boreal Forests
[C]//
Polarimetric properties of burned forest areas at C- and L-band
[J].
林火烈度遥感评估指数适应性分析
[J].<p>科学合理地定量评估林火烈度,对揭示林火干扰下森林生态系统的变化,以及植被的恢复与管理具有重要意义。以美国科罗拉多大峡谷国家公园北缘的Poplar Fire为实验区,利用Landsat5 TM影像,结合实地调查的综合火烧指数(composite burn index,CBI),分析评价了归一化差值植被指数(normalized difference vegetation index,NDVI),归一化火烧指数(normalized burn ration,NBR),差分归一化植被指数(differenced normalized difference vegetation index,ΔNDVI)和差分归一化火烧指数(differenced normalized burn ration,ΔNBR)4种遥感指数对林火烈度评估的适应性。结果表明,4种遥感指数对识别不同等级林火烈度存在一定的差异。在未过火区和轻度火灾区,单一遥感指数的精度略高于差分遥感指数,其中NBR的提取精度最高,分别达到了66.7%和80%; 在中度火灾区和重度火灾区,差分遥感指数的精度高于单一遥感指数,ΔNBR的提取精度最高,分别达到了100%和90%。总体上,基于差分遥感指数的林火烈度制图精度总体高于单一遥感指数,其中ΔNBR的总体制图精度最高,达到了 86.2%。因此,ΔNBR是林火烈度分析与评估的适宜遥感指数。</p>
An adaptability analysis of remote sensing indices in evaluating fire severity
[J].
典型植被指数识别火烧迹地潜力分析
[J].
Potential analysis of typical vegetation index for identifying burned area
[J].
极化SAR林火燃烧面积提取与燃烧强度估计
[D].
Polarimetric SAR Forest Fire Scars Extraction and Burn Severity Estimation
[D].
重轨卫星InSAR技术在多云雨山区DEM生成及更新中的应用研究
[D].
Repeat-Pass Spaceborne SAR Interferometry for DEM Generation in Cloudy-Rainy Mountianous Area
[D].
林火干扰区全极化SAR影像的散射特性分析
[J].<p>目前,利用合成孔径雷达(synthetic aperture Radar,SAR)影像进行林火监测主要限于从林火前后时相的单极化通道振幅数据或者火后全极化影像振幅数据开展,而从林火对森林散射机制改变角度开展多时相全极化SAR林火监测的研究还较少。以2009年阿拉斯加地区发生的林火为例,以林火发生前后获取的Radarsat-2全极化影像为实验数据,从林火所改变的森林后向散射强度和散射机制角度,对林火发生前后各个极化通道后向散射强度、主导散射机制和去极化作用参数进行了定量分析,并对变化原因给出了解释。研究表明,对于北方森林,林火发生后林火干扰区的后向散射强度在同极化通道相比林火前增加了20%,交叉极化通道也有小幅增加,主导散射由林火前占59%的体散射变为林火后占53%的表面散射,去极化作用相比林火前减小45%。这些结论对于利用多时相全极化SAR数据开展林火燃烧面积或者燃烧强度监测具有参考价值。</p>
Analysis of fire disturbed forests scattering characteristics using polarimetric SAR image
[J].
基于面向对象方法的Sentinel-1A SAR在洪水监测中的应用
[J].
DOI:10.12082/dqxxkx.2018.170599
Magsci
[本文引用: 1]
<p>开展Sentinel-1A SAR数据在洪水淹没范围提取和水体变化监测方面的应用研究,对科学有效地管理洪涝灾害有重要意义。合成孔径雷达以其不受天气影响、能穿透云层、覆盖面积广等特点成为灾害监测的重要数据来源。面向对象的方法能有效解决影像的椒盐现象被广泛运用于信息提取研究。本文基于Sentinel-1A SAR数据,利用面向对象的方法构建洪水淹没范围提取流程,绘制灾前、灾中、灾后水体变化监测图,对比分析基于传统像元的提取方法,实现对广西临桂会仙岩溶湿地区域不同时期洪水动态监测。研究表明,Sentinel-1A SAR数据在洪水监测领域有巨大的应用潜力,相较于传统基于像元的方法,面向对象的方法能有效抑制杂斑生成,提高空间信息的利用效率,具有更好的提取精度。</p>
Flooding monitoring application based on the object-oriented method and Sentinel-1A SAR data
[J].
面向对象的高分辨率SAR图像处理及应用
[J].
The processing and application of high resolution SAR images with object-based image analysis
[J].
光谱、形状特征结合的多精度图像分割算法与应用
[J].由于高分辨率遥感图像的数据量和计算复杂性骤增,影像噪声、光谱混淆现象更为突出,这对传统主要依据于像元光谱特征的图像处理与分析方法提出了极大挑战。为此,特征更为丰富、在噪声处理和知识融合上,更具优势的面向对象的图像分析方法,逐步成为高分辨率遥感图像应用的研究热点。面向对象图像分析的第一步,以及关键一步是图像分割。本文设计了多精度图像分割算法:(1)采用降水分水岭变换进行初步分割获取图像次一级斑块,即分割亚基元;(2)设计一种可重复合并的快速图斑合并方法,进行亚基元的层次归并获得最后分割斑块,完成图像分割。在合并过程中,斑块之间的差异指标是其光谱合并代价、形状合并代价的加权和,合并结束的标志是斑块间两两合并代价超过尺度参数的平方。设置不同的尺度参数,则可实现多精度图像分割过程。实验证明,方法分割效果较好,并在算法效率上满足实际应用需求,可以开展后续图像分类、专题信息提取等工作。
Multi-resolution remotely sensed image segmentation in combination with spectral and shape features
[J].
Object-oriented and multi-scale image analysis in semantic networks
[C]//
Multiresolution segmentation:An optimization approach for high quality multiscale image segmentation
[J].
基于极化特征分析的卫星SAR林火燃烧强度估计研究
[D].
Forest Fire Severity Estimation with Satellite Polarimetric SAR Based on Polarization Feature
[D].
洪泽湖湿地纹理特征参数分析
[J].<p>应用纹理特征进行影像分类,关键在于纹理特征参数的确定。以洪泽湖湿地典型地区为研究对象,选择灰度共生矩阵进行纹理特征计算,探讨灰度共生矩阵窗口尺寸、移动步长、方向和纹理特征统计量对淡水湖泊湿地的区分能力; 然后,利用纹理特征和地物光谱特征,结合决策树方法对研究区湿地及其他主要地类进行分类,并通过混淆矩阵进行精度评价。结果表明: 研究区湿地分类中纹理特征的最佳窗口大小为3像元×3像元,方向为90°,步长为1个像元,纹理特征统计量组合为均值、熵和相关度; 分类精度为83.24%,<em>Kappa</em>为0.788,其结果验证了纹理特征参数选择的科学性和合理性。</p>
Parameter analysis of image texture of wetland in the Hongze Lake
[J].
Spatial analysis of remote sensing image classification accuracy
[J].
A note on sample size estimation for multinomial populations
[J].
Analysis of the impact of sample size,attribute variance and within-sample choicedistribution on the estimation accuracy of multinomial logit models using simulated data
[J].
/
〈 |
|
〉 |
