海洋学报

波段	波长范围/nm	空间分辨率/m
1	420～500	50
2	520～600	50
3	610～690	50
4	760～890	50

波段	波长范围/nm	空间分辨率/m
1	420～500	50
2	520～600	50
3	610～690	50
4	760～890	50

名称	英文全称	简写	公式	选择依据
注：NIR代表近红外波段，R代表红波段，G代表绿波段，B代表蓝波段。
归一化差值植被指数	Normalized Difference Vegetation Index	NDVI	$(NIR-R)/(NIR+R) $	对绿色植被敏感，常用于植被状态研究
比值植被指数	Ratio Vegetation Index	RVI	$R/NIR $	对浓密覆盖植被敏感
红色植被指数	Red Vegetation Index	RI	$(R-G)/(R+G) $	对土壤颜色影响的植被指数的校正
结构不敏感色素指数	Structure Insensitive Pigment Index	SIPI	$(NIR-B)/(NIR+R) $	标识植被冠层胁迫性的增加
归一化差异水分指数	Normalized Difference Water Index	NDWI	$(G-NIR)/(G+NIR) $	对土壤湿度敏感
归一化差异绿度指数	Normalized Difference Greenness Index	NDGI	$(G-R)/(G+R) $	用于检验不同活力植被形式
修改型土壤调节植被指数	Modified Soil-adjusted Vegetation Index	MSAVI	$[2NIR + 1 - \sqrt { { {\left( {2NIR + 1} \right)}^2} - 8\left( {NIR - R} \right)} ]/2$	调整土壤背景对植被指数的影响，并减少裸土的影响
增强型植被指数	Enhanced Vegetation Index	EVI	$2.5[(NIR-R)/(NIR+6R-7.5B+1)] $	同时修订土壤背景和大气噪声的影响
大气阻抗植被指数	Atmospheric Impedance Vegetation Index	ARVI	$[NIR-(2R-B)]/[NIR+(2R-B)] $	减少大气散射对植被指数的影响
CZI蓝红波段比	CZI Blue-Red Wave Segment Ratio	CBRI	$(B-R)/(B+R) $	反映不同地物在CZI波段之间差异性变化
CZI蓝绿波段比	CZI Blue-Green Wave Segment Ratio	CBGI	$(B-G)/(B+G) $	反映不同地物在CZI波段之间差异性变化
CZI近红外蓝光波段比	CZI Near Infrared-Blue Wave Segment Ratio	CNBI	$(NIR-B)/(NIR+B) $	反映不同地物在CZI波段之间差异性变化
CZI可见光光谱斜率比	CZI Visible Spectrum Slope Ratio	CVSSR	$\left( {\dfrac{{R - G}}{{650 - 560}}} \right)\bigg/\left( {\dfrac{{G - B}}{{560 - 460}}} \right)$	反映不同植被在可见光波段的光谱形态差异

名称	英文全称	简写	公式	选择依据
注：NIR代表近红外波段，R代表红波段，G代表绿波段，B代表蓝波段。
归一化差值植被指数	Normalized Difference Vegetation Index	NDVI	$(NIR-R)/(NIR+R) $	对绿色植被敏感，常用于植被状态研究
比值植被指数	Ratio Vegetation Index	RVI	$R/NIR $	对浓密覆盖植被敏感
红色植被指数	Red Vegetation Index	RI	$(R-G)/(R+G) $	对土壤颜色影响的植被指数的校正
结构不敏感色素指数	Structure Insensitive Pigment Index	SIPI	$(NIR-B)/(NIR+R) $	标识植被冠层胁迫性的增加
归一化差异水分指数	Normalized Difference Water Index	NDWI	$(G-NIR)/(G+NIR) $	对土壤湿度敏感
归一化差异绿度指数	Normalized Difference Greenness Index	NDGI	$(G-R)/(G+R) $	用于检验不同活力植被形式
修改型土壤调节植被指数	Modified Soil-adjusted Vegetation Index	MSAVI	$[2NIR + 1 - \sqrt { { {\left( {2NIR + 1} \right)}^2} - 8\left( {NIR - R} \right)} ]/2$	调整土壤背景对植被指数的影响，并减少裸土的影响
增强型植被指数	Enhanced Vegetation Index	EVI	$2.5[(NIR-R)/(NIR+6R-7.5B+1)] $	同时修订土壤背景和大气噪声的影响
大气阻抗植被指数	Atmospheric Impedance Vegetation Index	ARVI	$[NIR-(2R-B)]/[NIR+(2R-B)] $	减少大气散射对植被指数的影响
CZI蓝红波段比	CZI Blue-Red Wave Segment Ratio	CBRI	$(B-R)/(B+R) $	反映不同地物在CZI波段之间差异性变化
CZI蓝绿波段比	CZI Blue-Green Wave Segment Ratio	CBGI	$(B-G)/(B+G) $	反映不同地物在CZI波段之间差异性变化
CZI近红外蓝光波段比	CZI Near Infrared-Blue Wave Segment Ratio	CNBI	$(NIR-B)/(NIR+B) $	反映不同地物在CZI波段之间差异性变化
CZI可见光光谱斜率比	CZI Visible Spectrum Slope Ratio	CVSSR	$\left( {\dfrac{{R - G}}{{650 - 560}}} \right)\bigg/\left( {\dfrac{{G - B}}{{560 - 460}}} \right)$	反映不同植被在可见光波段的光谱形态差异

光谱指数	与NDVI之间相关系数	光谱指数	与NDVI之间相关系数
RVI	−0.985	MSAVI	0.978
RI	0.413	EVI	0.989
SIPI	0.989	CBRI	−0.004
NDWI	−0.996	CBGI	0.608
NDGI	−0.413	ARVI	−0.238
CNBI	0.996	CVSSR	−0.063

光谱指数	与NDVI之间相关系数	光谱指数	与NDVI之间相关系数
RVI	−0.985	MSAVI	0.978
RI	0.413	EVI	0.989
SIPI	0.989	CBRI	−0.004
NDWI	−0.996	CBGI	0.608
NDGI	−0.413	ARVI	−0.238
CNBI	0.996	CVSSR	−0.063

基于HY-1C CZI影像光谱指数重构数据MNF变换的红树林提取

PDF下载

梁超 ¹^,² , 刘利 ³^,^* , 刘建强 ¹^,² , 邹斌 ¹^,² , 邹亚荣 ¹^,² , 崔松雪 ¹^,²

海洋学报 | 海洋技术 2020,42(4): 104-112

收起

海洋学报 | 海洋技术 2020, 42(4): 104-112

基于HY-1C CZI影像光谱指数重构数据MNF变换的红树林提取

全屏

梁超¹^,², 刘利³^,^*, 刘建强¹^,², 邹斌¹^,², 邹亚荣¹^,², 崔松雪¹^,²

作者信息

¹ 自然资源部国家卫星海洋应用中心，北京 100081

² 自然资源部空间海洋遥感与应用研究重点实验室，北京 100081

³ 中国科学院空天信息创新研究院，北京 100094

梁超（1985-），男，陕西省咸阳市人，主要研究方向为卫星海洋遥感应用。E-mail：liangchao@mail.nsoas.org.cn

通讯作者:

^*刘利（1984-），女，江苏省淮安市人，工程师，主要从事遥感数据分析与应用研究。E-mail：liuli03rs@163.com

Extracting mangrove information using MNF transformation based on HY-1C CZI spectral indices reconstruction data

Chao Liang¹^,², Li Liu³^,^*, Jianqiang Liu¹^,², Bin Zou¹^,², Yarong Zou¹^,², Songxue Cui¹^,²

Affiliations

¹ National Satellite Ocean Application Service, Ministry of Natural Resources, Beijing 100081, China

² Key Laboratory of Space Ocean Remote Sensing and Application, Ministry of Natural Resources , Beijing 100081, China

³ Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China

出版时间: 2020-04-25 doi: 10.3969/j.issn.0253-4193.2020.04.012

文章导航

摘要

收起

本文基于广西山口国家红树林生态自然保护区的HY-1C卫星的海岸带成像仪(Coastal Zone Imager，CZI)影像，分析了红树林与一般陆地植被的光谱特征及其光谱指数的相关性，采用归一化差值植被指数(Normalized Difference Vegetation Index，NDVI)、归一化差异水分指数(Normalized Difference Water Index，NDWI)、大气阻抗植被指数(Atmospheric Impedance Vegetation Index，ARVI)及利用CZI波段构建的光谱斜率比(CZI Visible Spectrum Slope Ratio，CVSSR)4个指数替代CZI原始波段形成重构数据，基于重构数据的最小噪声分离变换(Minimum Noise Fraction Rotation，MNF)结果分量，建立决策树并实现了红树林信息的自动提取。研究结果表明：结合本文所选光谱指数重构数据及MNF变换方法，能够有效增强CZI影像上红树林与一般陆地植被的光谱差异，基于MNF变换分量建立的决策树可有效提取红树林信息，经与专家解译结果比对，本文方法面积准确率达90%以上；经随机样本点验证，总体检测精度为88%。

关键词

HY-1C / 海岸带成像仪 / 红树林 / 光谱指数 / 最小噪声分离变换

Abstract

收起

In this study, we first used the spectral vegetation indices such as normalized difference vegetation index (NDVI), normalized difference water index (NDWI), atmospheric impedance vegetation index (ARPI) and visible spectrum slope ratio of coastal zone imager (CVSSR) to reconstruct the HY-1C coastal zone imager (CZI) image data of the Shankou mangrove national ecosystem nature reserve in Guangxi. And then, the minimum noise fraction rotation (MNF) was used to enhance the spectral difference between mangroves and general terrestrial vegetation on the reconstructed multi-band data set. We established a decision tree based on the MNF components to achieve automatic extracting mangrove information. The results show that the spectral indices reconstruction data and its MNF transformation can effectively enhance the difference between the mangroves and the general terrestrial vegetation on CZI images, the mangrove information can be effectively extracted by our decision tree. Compared with the experts’ interpretation results, the extracted accuracy of area of our method is over 90%. The overall detection accuracy is 88% after verification by random sample points.

Key words

HY-1C satellite / CZI / mangrove / spectral indices / MNF transformation

引用本文

梁超, 刘利, 刘建强, 邹斌, 邹亚荣, 崔松雪. 基于HY-1C CZI影像光谱指数重构数据MNF变换的红树林提取. 海洋学报, 2020 , 42 (4) : 104 -112 . DOI: 10.3969/j.issn.0253-4193.2020.04.012

Chao Liang, Li Liu, Jianqiang Liu, Bin Zou, Yarong Zou, Songxue Cui. Extracting mangrove information using MNF transformation based on HY-1C CZI spectral indices reconstruction data[J]. Haiyang Xuebao, 2020 , 42 (4) : 104 -112 . DOI: 10.3969/j.issn.0253-4193.2020.04.012

正文

收起

1 引言

收起

红树林是指生长在热带、亚热带近岸潮间带上部滩涂浅滩，以红树植物为主体的常绿灌木或乔木组成的潮滩湿地木本生物群落，是陆地向海洋过渡的特殊生态系统。红树林具有防风消浪、促淤保滩、净化海水、保持生物多样性等重要功能。近几十年来，由于围海造地、围海养殖、砍伐等人为不合理开发活动等因素影响，使得红树林面积逐年锐减。红树林生态环境保护工作对长期监测红树林的分布范围、生物量、健康水平等提出了要求，其中红树林分布范围的监测是基础。利用遥感手段可以大范围、快速、重复监测红树林时空分布信息，这对于研究保护红树林生态系统具有重要意义。

国内外红树林遥感监测方面开展了大量研究，所采用的数据源通常为NOAA-AVHRR、Landsat、SPOT、Quickbird、ALOS/PALSAR等卫星以及诸如HyMap等航空高光谱遥感影像，主要方法包括人工解译、监督分类、植被指数及纹理特征、面向对象分类等^[1-14]。HY-1C是我国第一代海洋水色卫星HY-1A/B的后继星，于2018年9月7日在太原卫星发射中心发射成功，是我国首颗业务化运行、全球开机、连续观测的海洋水色卫星，载荷配置和性能指标、探测能力等均有明显改进，可获取重访周期更短、观测范围更大、精度更高的水色遥感产品。HY-1C搭载的海岸带成像仪（Coastal Zone Imager，CZI）用于获取海陆交互区域实时数据，监测近海、海岛、海岸带等信息，其幅宽、分辨率相比上一代载荷分别提升近2倍和5倍。本文拟以广西山口红树林生态自然保护区为研究区，综合利用红树林与一般陆地植被光谱特征及空间分布特征，选用适当的光谱指数重构原始数据，经最小噪声分离变换（Minimum Noise Fraction Rotation，MNF）后，建立决策树提取红树林，据此探讨HY-1C卫星CZI数据在红树林监测工作中的适用性。

2 数据与方法

收起

2.1 研究区概况

本文以广西山口红树林生态自然保护区为研究区（图1），该保护区由位于广西合浦县东南部沙田半岛东、西两侧的海域、陆域及全部滩涂组成，地处亚热带，是我国第二个国家级的红树林自然保护区，分布着发育良好、结构典型、连片较大、保存较完整的天然红树林，有红海榄、木榄、秋茄、桐花树等15种红树林植物。山口红树林国家级自然保护区具典型的大陆红树林海岸生态系统特征，红树林中还栖息着多种海洋生物和鸟类，具有重要的科学价值。

2.2 数据源

采用2018年10月31日11时25分获取的HY-1C卫星CZI传感器L1B级数据，CZI主要技术指标见表1。

2.3 研究方法

本文采用波段重构与光谱变换相结合的方法，以增强CZI影像上红树林与一般陆地植被的特征差异，提高二者可分性。图2为研究方法流程图，首先对CZI数据预处理，然后基于红树林与一般陆地植被光谱特征构建多种光谱指数，通过相关性分析筛选最佳指数组合重构原始数据，继而采用MNF进一步增强重构数据中红树林与一般陆地植被的光谱及空间分布特征差异，最后建立决策树提取红树林。

2.3.1 数据预处理

数据预处理包括几何校正和辐射定标，CZI数据采用HDF5格式，几何校正利用数据自带的RPC参数。L1B数据像元值为辐射亮度，据公式（1）转换为表观反射率^[15]：

(1)$\rho = \frac{{{\text π} \cdot {L_{\text λ} } \cdot {D^2}}}{{ESU{N_{\text λ} } \cdot {\rm{cos}}\theta }},$

式中，ρ为大气层顶表观反射率（无量纲），L_λ为大气层顶进入卫星传感器的光谱辐射亮度，单位为W/(m²·sr·μm)，D为日地距离(天文单位)，ESUN_λ为大气层顶的平均太阳光谱辐照度，单位为W/(m²·μm)，θ为太阳天顶角，相关参数可由HDF5数据集属性参量中查找。

2.3.2 光谱指数

红树林与一般陆地植被光谱特征相似，在CZI原始影像难以区分。为突出植被、水体等信息并增强红树林与一般陆地植被的差异，本文拟选用适当的光谱指数重构CZI原始数据。综合考虑大气、红树林分布环境等因素的影响，本文选用了常见的8种植被指数^[16-19]和1种水体指数，为进一步突出CZI影像上植被光谱特征，构造了3个波段比指数和1个光谱斜率比指数，其中波段比指数用来表达目标物在不同波段之间的归一化差值比，可见光光谱斜率比指数用来表达不同地物在可见光波段的光谱曲线形态差异。表2列出了本文所采用的13种光谱指数。

2.3.3 最小噪声分离变换

MNF最早由Green等^[20]提出，是在主成分变换基础上提出的一种线性变换算法。主成分变换各分量按照方差降序排列，影像主要信息集中在前几个分量，但不能保证影像质量（信噪比）也按照降序排列，MNF用以改进主成分变换在噪声消除及影像增强上的不足。Boardman和Kruse^[21]证明了MNF等效于连续两次主成分变换，设多光谱数据由加性噪声和信号两部分组成，首先对影像噪声协方差矩阵进行变换，使变换后噪声协方差矩阵为单位矩阵且波段间不相关，然后对变换后的数据集作标准主成分变换^[22]。

第一步，对多光谱影像X的噪声协方差矩阵进行对角化：

(2)$\boldsymbol{\varLambda } = {\boldsymbol{E}^{\rm{T}}}\boldsymbol{C}_{{x}}^{\rm{n}}\boldsymbol{E},$

式中，C_xⁿ为原始影像噪声协方差矩阵，Λ为其特征值组成的对角阵，E为由其特征向量组成的正交矩阵。

上式前乘以(Λ^−0.5)^T，后乘以Λ^−0.5，则

(3)$\boldsymbol{I} = {({\boldsymbol{\varLambda }^{ - 0.5}})^{\rm{T}}}{\boldsymbol{E}^{\rm{T}}}\boldsymbol{C}_{{x}}^{\rm{n}}\boldsymbol{E}{\boldsymbol{\varLambda }^{ - 0.5}},$

式中，I为单位矩阵，令F=EΛ^−0.5，则上式转换为

(4)$\boldsymbol{I} = {\boldsymbol{F}^{\rm{T}}}\boldsymbol{C}_{{x}}^{\rm{n}}\boldsymbol{F}.$

原始影像经Y=F^TX变换后的协方差矩阵为

(5)${\boldsymbol{C}_Y} = {\boldsymbol{F}^{\rm{T}}}{\boldsymbol{C}_{{x}}}\boldsymbol{F},$

式中，C_x为原始影像协方差矩阵。

第二步，对Y应用标准主成分变换：

(6)$\begin{aligned}& \boldsymbol{Z} = {\boldsymbol{B}^{\rm{T}}}\boldsymbol{Y},\\& {\boldsymbol{B}^{\rm{T}}}{\boldsymbol{C}_Y}\boldsymbol{B} = {\boldsymbol{\varLambda }_Y},\\& {\boldsymbol{B}^{\rm{T}}}\boldsymbol{B} = \boldsymbol{I},\end{aligned}$

式中，Λ_Y是C_Y特征值降序排列的对角阵，B是其特征向量组成的正交矩阵。

综上，MNF矩阵为

(7)${\boldsymbol{T}_{{\rm{MNF}}}} = \boldsymbol{FB}.$

多光谱数据经过MNF后各分量按照信噪比降序排列，影像有效信息集中在前几个分量，噪声主要存在于后面的分量中，实现了光谱增强的同时分离数据中的噪声，改善多光谱处理结果。

MNF最重要的一步是估计原始影像噪声协方差，其方法一般借助影像的空间特征，据Canty^[23]：

(8)$\boldsymbol{C}_{{x}}^{\rm{n}} \approx \frac{1}{2}\left\langle {\left( {\boldsymbol{X}\left( x \right) - \boldsymbol{X}\left( {x + {{h}}} \right)} \right){{\left( {\boldsymbol{X}\left( x \right) - \boldsymbol{X}\left( {x + {{h}}} \right)} \right)}^{\rm{T}}}} \right\rangle .$

式中，x=(i, j)^T表示影像X的一个像元坐标，h=(h1, h2)^T表示坐标的偏移量，即噪声协方差可以用原始影像与偏移影像差异的协方差来估计。

3 结果分析

收起

3.1 CZI影像上植被光谱特征分析

通过建立感兴趣区，提取红树林与一般陆地植被样本的表观反射率光谱曲线（图3），可见在CZI影像上二者均表现出典型的绿色植被光谱特征，即绿波段强反射、红波段强吸收以及近红外波段的高反射率，蓝波段由于大气散射的影响亦表现出较强的反射，而近红外波段上二者光谱曲线均较为发散，这是由于中低分辨率影像上混合像元效应的影响。红树林与一般陆地植被光谱曲线具有较高的相似性，使得二者很难在4波段CZI数据中被区分。图4为CZI影像上红树林与一般陆地植被表观反射率分布直方图，可见二者在各波段具有大致相近的直方图形态和峰值分布区间，虽然在第4波段，红树林与一般陆地植被直方图峰值位置略有差异，但如前述，红树林与陆地植被在该波段光谱曲线均较为发散，直方图上分布范围均较宽，从而存在较大的重叠区。综合考虑光谱曲线及直方图可知，单纯依靠CZI表观反射率数据很难区分红树林与一般陆地植被。

3.2 基于光谱指数重构多光谱数据

红树林与一般陆地植被在CZI原始光谱空间中较高的相似性对红树林检测带来了难度，传统的监督/非监督分类等方法提取红树林效果均不甚理想。本文通过构建对植被和水体信息敏感的光谱指数组合代替原始波段数据，以增强影像中典型植被的光谱信息及其差异性。

本文以植被遥感中最常用的归一化差异植被指数NDVI为参考基准，开展了NDVI与其他12种光谱指数的相关性分析，表3为各指数之间的相关系数。首先剔除与NDVI相关系数极低的光谱指数，如CBRI，因为极低的相关系数意味着该类指数对目标植被信息不敏感；其次，剔除与NDVI具有较高的正相关系数的光谱指数，如SIPI、CNBI、MSAVI、EVI等，该类指数在表征植被信息方面与NDVI作用相当，从避免冗余考虑予以剔除；最后，选择与NDVI具有中、低相关性及最大负相关性的光谱指数参与重构，如RI、NDWI、NDGI、CBGI、ARVI和CVSSR等，该类指数既能从不同程度反映植被信息，又避免了信息冗余，确保重构后的新光谱空间既能最大化地增强植被及水体信息，又能最大程度减少波段之间的相关性。据此原则并经多次试验，本文最终选取NDVI、NDWI、ARVI和CVSSR为最佳指数组合，并以之依次代替CZI 4个原始波段组成新的多光谱数据。

图5为重构后多光谱数据中红树林与陆地植被的直方图，对比图4可见在新光谱空间中，第1波段上红树林与陆地植被分布大致相近，但红树林分布范围更宽一些；第2波段上红树林的分布峰值位置约为−0.50，而一般陆地植被则小于−0.50；第3波段上二者分布特征最为接近；第4 波段上红树林分布峰值位置约为0.70，较之一般陆地植被偏大。这说明在重构数据光谱空间中红树林与一般陆地植被的差异性得到一定程度增强。

3.3 最小噪声分离变换

对重构后数据进行MNF，图6为变换结果的4个分量图像，根据MNF原理，变换后4个分量按照噪声水平升序排列，第1、2、3分量包含了主要的空间信息，而第4分量则包含较多的噪声信号。特别地，在MNF第3分量上，红树林与一般陆地植被表现出较明显的目视差异，前者具有更高的像元值。由图7可知，MNF后，红树林与一般植被直方图分布形态差异性得到提升，峰值分布区间的重合度进一步减小。具体表现为：第1分量上，红树林分布范围更宽，峰值位置约在9.30附近，陆地植被分布范围集中在10以上，峰值位于11.52附近；第2分量上，红树林分布范围亦较宽，峰值位置约在−1.23附近，陆地植被接近均值为−2.13的正态分布，主要范围在−3.77至1.77之间；第3分量上红树林与一般陆地植被分布差异最为明显，红树林分布中心偏右，峰值位置约为5.13，陆地植被分布峰值位置大致为3.62，且绝大部分分布在约4.36以下，而红树林一般高于此值（见图7中黑色虚线）；第4分量包含更多噪声信号，二者区分不大。这说明，基于光谱指数重构数据，MNF在分离了噪声信号的同时，由于利用了影像空间分布信息（公式（8）），进一步增强了目标地物的波谱差异，提高了红树林与一般陆地植被的可区分度。

通过上述分析，基于光谱指数重构数据MNF结果可以较容易地构建分类决策树，实现红树林提取。图8为本文构建的决策树，其中：

规则1：NDVI > 0.50，区分植被与非植被。

规则2：MNF3 > 4.36，区分红树林与一般陆地植被。

国家卫星海洋应用中心在诸如“908”专项调查、广西海域使用本底库建设等工作中在该区域开展过多次遥感调查及现场核查，图9a是基于历史调查成果和专家经验支持下的山口红树林分布范围人工解译结果，图9b是基于本文决策树的红树林自动提取结果，可见绝大部分红树林均能被正确提取，其中专家解译的面积为9.44 km²，本文方法提取红树林的面积为9.65 km²，面积检测相对准确率可达97%。进一步在研究区内随机选取2 000个验证样本点，其中覆盖红树林样本点50个，本文方法检测出样本点44个，从而红树林总体检测精度约为88%。

4 总结

收起

本文基于HY-1C卫星CZI 4波段影像，以广西山口红树林自然保护区为研究区，采用光谱指数重构CZI数据，经MNF建立决策树，实现了研究区红树林自动提取，研究结果表明：利用NDVI、NDWI、ARVI、CVSSR 4种指数重构CZI数据，减少了地形、大气、传感器等因素对植被光谱信息的影响，且重构后数据对植被、水体等更加敏感。MNF则进一步利用了影像上红树林与陆地植被的波谱及空间分布特征，变换结果中二者差异被显著增强，最终通过建立决策树提取了红树林分布信息。本文结合光谱指数重构和MNF方法，较好地解决了在CZI原始光谱空间中红树林与陆地植被难以区分的问题。因此，HY-1C CZI数据可以有效地用于红树林空间分布监测。

HY-1C CZI数据优势在于重访周期短、幅宽大，可快速重复获取大范围中等分辨率光学影像数据，但因波段较少，缺少对研究土壤与水分情况较重要的短波红外、中波红外或热红外等波段，且波段较宽，对典型植被波谱特征细节描述能力稍弱。本文研究方法可弥补CZI数据用于红树林监测的上述不足，后续还可通过形态学方法对提取结果进一步优化，提升自动检测准确率。需要指出的是，本文决策树规则阈值对于不同地区和时相的数据可能会有所调整，目前HY-1C卫星尚处于在轨测试阶段，相信卫星正式交付后，传感器数据的定标精度和质量稳定性会得到进一步保障，有利于后续对本文方法的验证和完善。

基金

收起

国家重点研发计划（2018YFB0505001-04）。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Aschbacher J, Ofren R, Delsol J P, et al. An integrated comparative approach to mangrove vegetation mapping using advanced remote sensing and GIS technologies: preliminary results[J]. Hydrobiologia, 1995, 295: 285−294.

Green E P, Clark C D, Mumby A J, et al. Remote sensing techniques for mangrove mapping[J]. International Journal of Remote Sensing, 1998, 19(5): 935−956.

Green E P, Mumby P J, Edwards A J, et al. The assessment of mangrove areas using high resolution multispectral airborne imagery[J]. Journal of Coastal Research, 1998, 14(2): 433−443.

Wang Lanfa, Sousa W, Gong P. Integration of object-based and pixel-based classification for mapping mangroves with IKONOS imagery[J]. International Journal of Remote Sensing, 2004, 25(24): 5655−5668.

肖海燕, 曾辉, 昝启杰, 等. 基于高光谱数据和专家决策法提取红树林群落类型信息[J]. 遥感学报, 2007, 11(4): 531−537.

Xiao Haiyan, Zeng Hui, Zan Qijie, et al. Decision tree model in extraction of mangrove community information using hyperspectral image data[J]. Journal of Remote Sensing, 2007, 11(4): 531−537.

张雪红. 基于知识与规则的红树林遥感信息提取[J]. 南京信息工程大学学报, 2011, 3(4): 341−345.

Zhang Xuehong. Remote sensing information extraction of mangrove based on knowledge and rules[J]. Journal of Nanjing University of Information Science & Technology, 2011, 3(4): 341−345.

Heumann B W. Satellite remote sensing of mangrove forests: recent advances and future opportunities[J]. Progress in Physical Geography: Earth and Environment, 2011, 35(1): 87−108.

张雪红. 基于决策树方法的Landsat8 OLI影像红树林信息自动提取[J]. 国土资源遥感, 2016, 28(2): 182−187.

Zhang Xuehong. Decision tree algorithm of automatically extracting mangrove forests information from Landsat 8 OLI imagery[J]. Remote Sensing for Land & Resources, 2016, 28(2): 182−187.

Giri C. Observation and monitoring of mangrove forests using remote sensing: opportunities and Challenges[J]. Remote Sensing, 2016, 8(9): 783.

Cao Jingjing, Leng Wanchun, Liu Kai, et al. Object-based mangrove species classification using unmanned aerial vehicle hyperspectral images and digital surface models[J]. Remote Sensing, 2018, 10(1): 89.

Abdel-Hamid A, Dubovyk O, El-Magd I, et al. Mapping mangroves extents on the red sea coastline in Egypt using polarimetric SAR and high resolution optical remote sensing data[J]. Sustainability, 2018, 10(3): 646.

Zulfa A W, Norizah K. Remotely sensed imagery data application in mangrove forest: a review[J]. Pertanika Journal of Science & Technology, 2018, 26(3): 899−922.

Wang Le, Jia Mingming, Yin Dameng, et al. A review of remote sensing for mangrove forests: 1956-2018[J]. Remote Sensing of Environment, 2019, 231: 111223.

Pham T D, Yokoya N, Bui D T, et al. Remote sensing approaches for monitoring mangrove species, structure, and biomass: opportunities and challenges[J]. Remote Sensing, 2019, 11(3): 230.

Markham B L, Barker J L. Thematic Mapper bandpass solar exoatmospheric irradiances[J]. International Journal of Remote Sensing, 1987, 8(3): 517−523.

田庆久, 闵祥军. 植被指数研究进展[J]. 地球科学进展, 1998, 13(4): 327−333.

Tian Qingjiu, Min Xiangjun. Advances in study on vegetation indices[J]. Advance in Earth Sciences, 1998, 13(4): 327−333.

郭铌. 植被指数及其研究进展[J]. 干旱气象, 2003, 21(4): 71−75.

Guo Ni. Vegetation index and its advances[J]. Arid Meteorology, 2003, 21(4): 71−75.

Muhsoni F F, Sambah A B, Mahmudi M, et al. Comparison of different vegetation indices for assessing mangrove density using sentinel-2 imagery[J]. International Journal of Geomate, 2018, 14(45): 42−51.

Gupta K, Mukhopadhyay A, Giri S, et al. An index for discrimination of mangroves from non-mangroves using LANDSAT 8 OLI imagery[J]. MethodsX, 2018, 5: 1129−1139.

Green A A, Berman M, Switzer P, et al. A transformation for ordering multispectral data in terms of image quality with implications for noise removal[J]. IEEE Transactions on Geoscience and Remote Sensing, 1988, 26(1): 65−74.

Boardman J W, Kruse F A. Automated spectral analysis: a geological example using AVIRIS data, north Grapevine Mountains, Nevada[C]//Proceedings, ERIM Tenth Thematic Conference on Geologic Remote Sensing. Ann Arbor, MI: Environmental Research Institute of Michigan, 1994: 407−418.

Richards J A. Remote Sensing Digital Image Analysis[M]. Berlin: Springer, 2013: 186−189.

Canty M J. Image Analysis, Classification and Change Detection in Remote Sensing: with Algorithms for ENVI/IDL[M]. New York: CRC Press, 2014.

2020年第42卷第4期

PDF下载

114

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.0253-4193.2020.04.012

接收时间：2019-05-17
首发时间：2026-03-26
出版时间：2020-04-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2019-05-17
修回日期：2019-08-21

基金

国家重点研发计划（2018YFB0505001-04）。

作者信息

¹ 自然资源部国家卫星海洋应用中心，北京 100081

² 自然资源部空间海洋遥感与应用研究重点实验室，北京 100081

³ 中国科学院空天信息创新研究院，北京 100094

通讯作者:

^*刘利（1984-），女，江苏省淮安市人，工程师，主要从事遥感数据分析与应用研究。E-mail：liuli03rs@163.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.3969/j.issn.0253-4193.2020.04.012

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

波段	波长范围/nm	空间分辨率/m
1	420～500	50
2	520～600	50
3	610～690	50
4	760～890	50

波段

波长范围/nm

空间分辨率/m

420～500

520～600

610～690

760～890

名称	英文全称	简写	公式	选择依据
注：NIR代表近红外波段，R代表红波段，G代表绿波段，B代表蓝波段。
归一化差值植被指数	Normalized Difference Vegetation Index	NDVI	$(NIR-R)/(NIR+R) $	对绿色植被敏感，常用于植被状态研究
比值植被指数	Ratio Vegetation Index	RVI	$R/NIR $	对浓密覆盖植被敏感
红色植被指数	Red Vegetation Index	RI	$(R-G)/(R+G) $	对土壤颜色影响的植被指数的校正
结构不敏感色素指数	Structure Insensitive Pigment Index	SIPI	$(NIR-B)/(NIR+R) $	标识植被冠层胁迫性的增加
归一化差异水分指数	Normalized Difference Water Index	NDWI	$(G-NIR)/(G+NIR) $	对土壤湿度敏感
归一化差异绿度指数	Normalized Difference Greenness Index	NDGI	$(G-R)/(G+R) $	用于检验不同活力植被形式
修改型土壤调节植被指数	Modified Soil-adjusted Vegetation Index	MSAVI	$[2NIR + 1 - \sqrt { { {\left( {2NIR + 1} \right)}^2} - 8\left( {NIR - R} \right)} ]/2$	调整土壤背景对植被指数的影响，并减少裸土的影响
增强型植被指数	Enhanced Vegetation Index	EVI	$2.5[(NIR-R)/(NIR+6R-7.5B+1)] $	同时修订土壤背景和大气噪声的影响
大气阻抗植被指数	Atmospheric Impedance Vegetation Index	ARVI	$[NIR-(2R-B)]/[NIR+(2R-B)] $	减少大气散射对植被指数的影响
CZI蓝红波段比	CZI Blue-Red Wave Segment Ratio	CBRI	$(B-R)/(B+R) $	反映不同地物在CZI波段之间差异性变化
CZI蓝绿波段比	CZI Blue-Green Wave Segment Ratio	CBGI	$(B-G)/(B+G) $	反映不同地物在CZI波段之间差异性变化
CZI近红外蓝光波段比	CZI Near Infrared-Blue Wave Segment Ratio	CNBI	$(NIR-B)/(NIR+B) $	反映不同地物在CZI波段之间差异性变化
CZI可见光光谱斜率比	CZI Visible Spectrum Slope Ratio	CVSSR	$\left( {\dfrac{{R - G}}{{650 - 560}}} \right)\bigg/\left( {\dfrac{{G - B}}{{560 - 460}}} \right)$	反映不同植被在可见光波段的光谱形态差异

名称

英文全称

简写

公式

选择依据

注：NIR代表近红外波段，R代表红波段，G代表绿波段，B代表蓝波段。

归一化差值植被指数

Normalized Difference Vegetation Index

NDVI

$(NIR-R)/(NIR+R) $

对绿色植被敏感，常用于植被状态研究

比值植被指数

Ratio Vegetation Index

RVI

$R/NIR $

对浓密覆盖植被敏感

红色植被指数

Red Vegetation Index

$(R-G)/(R+G) $

对土壤颜色影响的植被指数的校正

结构不敏感色素指数

Structure Insensitive Pigment Index

SIPI

$(NIR-B)/(NIR+R) $

标识植被冠层胁迫性的增加

归一化差异水分指数

Normalized Difference Water Index

NDWI

$(G-NIR)/(G+NIR) $

对土壤湿度敏感

归一化差异绿度指数

Normalized Difference Greenness Index

NDGI

$(G-R)/(G+R) $

用于检验不同活力植被形式

修改型土壤调节植被指数

Modified Soil-adjusted Vegetation Index

MSAVI

$[2NIR + 1 - \sqrt { { {\left( {2NIR + 1} \right)}^2} - 8\left( {NIR - R} \right)} ]/2$

调整土壤背景对植被指数的影响，并减少裸土的影响

增强型植被指数

Enhanced Vegetation Index

EVI

$2.5[(NIR-R)/(NIR+6R-7.5B+1)] $

同时修订土壤背景和大气噪声的影响

大气阻抗植被指数

Atmospheric Impedance Vegetation Index

ARVI

$[NIR-(2R-B)]/[NIR+(2R-B)] $

减少大气散射对植被指数的影响

CZI蓝红波段比

CZI Blue-Red Wave Segment Ratio

CBRI

$(B-R)/(B+R) $

反映不同地物在CZI波段之间差异性变化

CZI蓝绿波段比

CZI Blue-Green Wave Segment Ratio

CBGI

$(B-G)/(B+G) $

反映不同地物在CZI波段之间差异性变化

CZI近红外蓝光波段比

CZI Near Infrared-Blue Wave Segment Ratio

CNBI

$(NIR-B)/(NIR+B) $

反映不同地物在CZI波段之间差异性变化

CZI可见光光谱斜率比

CZI Visible Spectrum Slope Ratio

CVSSR

$\left( {\dfrac{{R - G}}{{650 - 560}}} \right)\bigg/\left( {\dfrac{{G - B}}{{560 - 460}}} \right)$

反映不同植被在可见光波段的光谱形态差异

光谱指数	与NDVI之间相关系数	光谱指数	与NDVI之间相关系数
RVI	−0.985	MSAVI	0.978
RI	0.413	EVI	0.989
SIPI	0.989	CBRI	−0.004
NDWI	−0.996	CBGI	0.608
NDGI	−0.413	ARVI	−0.238
CNBI	0.996	CVSSR	−0.063

光谱指数

与NDVI之间相关系数

光谱指数

与NDVI之间相关系数

RVI

−0.985

MSAVI

0.978

0.413

EVI

0.989

SIPI

0.989

CBRI

−0.004

NDWI

−0.996

CBGI

0.608

NDGI

−0.413

ARVI

−0.238

CNBI

0.996

CVSSR

−0.063