地空同步试验高光谱影像特征提取与分类
四川测绘地理信息局测绘技术服务中心,成都 610081
Feature extraction and classification of hyperspectral image with ground-sky synchronization test
Surveying and Mapping Technology Service Center, Sichuan Surveying and Mapping Geographic Information Bureau, Chengdu 610081, China
责任编辑: 陈理
收稿日期: 2018-08-20 修回日期: 2019-01-4 网络出版日期: 2019-09-15
基金资助: |
|
Received: 2018-08-20 Revised: 2019-01-4 Online: 2019-09-15
作者简介 About authors
廖小露(1986-),男,硕士研究生,工程师,主要研究方向为高光谱遥感技术应用、应急测绘保障和灾情地理空间分析。Email:liaoxlj@163.com.。 。
针对地空同步的高光谱特征提取与反演理论研究匮乏的现状,结合“类内密度最大,类间距离最大”原则,探索不同谱段对不同地物识别的可分性和重要性选择方法,并提出改进的投影寻踪分类方法,实现不同地物在特征波谱的加权投影寻踪。以山西省临汾市某苗圃为实验区,同步采集实验区不同地物光谱与PHI高光谱影像,依托实测光谱,围绕先整体最优后局部最优策略,构建不同分类规则,并将其应用于PHI高光谱影像分类,实现植被与非植被的信息提取,并细分为10余种不同地类。结果表明,420 nm,520 nm,570 nm,610 nm,660 nm,690 nm,715 nm及810 nm等8个谱段为不同植被的重要可分性谱段,分类结果具有层次感强、植被轮廓清晰、避免阴影影响等优点。
关键词:
In view of the lack of theoretical research on the extraction and inversion of hyperspectral features in ground-sky synchronization, the authors, in combination with the principle of “maximum density within class and maximum distance between classes”, studied the separability and importance selection of the spectra for different ground objects in different spectral regions, proposed an improved projection pursuit classification method, and realized the projection pursuit method based on weighted feature band. In the case study, the spectral and PHI hyperspectral images of different ground objects in the experimental area were collected synchronously and, combined with the measured spectra on the ground, constructed different classification rules around the strategy of overall optimization and local optimization. It was applied to the classification of PHI hyperspectral image to extract the information of vegetation and non-vegetation and to subdivide more than ten different kinds of ground objects. The results show that 8 spectral areas are important separability bands of different vegetation: 420 nm, 520 nm, 570 nm , 610 nm, 660 nm, 690 nm, 715 nm, and 810 nm. The classification results have the advantages of strong stratification, clear outline of vegetation and avoiding shadow influence.
Keywords:
本文引用格式
廖小露, 刘嘉, 周兴霞.
LIAO Xiaolu, LIU Jia, ZHOU Xingxia.
0 引言
当前,机器学习、计算机视觉和模式识别理论与方法已成为高光谱影像分类的重要技术手段。国内外学者为充分挖掘高光谱数据隐含的丰富信息和特征,发展了核变换技术[4,5,6]、特征挖掘技术[7,8]、半监督学习和主动学习[9,10]、光谱-空间分类[11]、稀疏表达[12]和多分类器集成[13,14]等一系列解决方案[15,16]。高光谱特征提取与反演方法处理思想多源于多光谱研究,存在以相对较粗层次的方法和思想解决较细层次问题的不足,如何在应用实践中充分发挥高光谱数据精细光谱、图谱合一等特点,开展针对高光谱影像“空间维、光谱维、空间-光谱维”综合知识的高光谱数据处理方法研究,选取高光谱影像的特征谱段以突出感兴趣区域,对有效地进行高光谱数据处理、分析及信息提取至关重要[17]。
本文从光谱分析的角度,探讨高分高光谱成像应用下 “地-空”一体的高光谱遥感成像应用实验,通过获取、分析采集的地面光谱特征,从已知光谱中选取能有效区分不同地物的特征谱段,采用改进的投影寻踪方法计算投影向量,构建分类规则,并结合航空高光谱影像开展分类研究,探索地面光谱与航空影像同步的应用方法。
1 光谱特征分析与提取
1.1 波段可分性与重要性分析
式中:
以
式中Z表示某一波段的像素值,0<i≤n,0<j≤n,n为波段数。
则类内密度
式中:
1.2 改进的投影寻踪方法
投影寻踪方法是将样本在整体上按照一定的约束条件实现样本最大差异化的投影[24],例如在局部投影点尽可能密集成若干个点团; 而在整体上投影点团之间尽可能分散。其默认各波谱范围上反射率相对于不同波段具有相同的重要性。但在实际应用中,不同地物在不同波段上具有不同的特性; 不同波段对不同地物的识别具有不同的响应特征,因此本文提出了改进的投影寻踪方法。
该方法首先根据不同波谱范围,讨论其对不同地物识别的重要性程度,再给出不同波段在不同地物识别中的权重系数,并在此基础上,对各波段反射率进行加权处理。具体步骤如下:
1)按波谱范围计算不同波段权重系数。该步骤的目的是突出波谱范围内不同地物的特征,以增强不同谱段的地物识别能力,提高可区分性,在式(1)的基础上,确定不同波谱范围上地物识别权重系数,权重系数的大小根据第i波段的
式中ω(ai)为权重系数。
2)反射率加权处理。由不同波段
式中:
3)投影寻踪分类。以投影寻踪方法为基础,构建指标使得投影后的值在局部投影点尽可能密集,且整体投影点团之间尽可能分散。结合目标函数式(1),实现非线性优化,约束条件为
改进的投影寻踪方法,本质上讲就是对数据分别进行2次不同方向投影的“综合”,第一次投影是针对各评价指标相对于评价目标的重要程度而进行的; 后一次加权是为了尽量拉开各被评价对象之间的整体差异而进行的。2次投影的背景是截然不同的,前者的投影方向由“功能驱动”决定,即波段重要性,而后者是由“差异驱动”生成。
2 应用实验验证
实验区位于山西省临汾市某苗圃,范围大小为400 m×200 m。实验以地面实测数据为基础,在进行不同谱段上地物的可分性与重要性分析基础上,确定分类特征谱段范围,并对波段进行加权处理,拉大波段间的差异性。结合改进的投影寻踪分析方法,以“地物逐类识别,逐层剔除”为目标,构建基于地面实测光谱的不同谱段不同地物识别分类规则,即完成具有监督样本信息的分类规则构建; 最后结合PHI航空遥感试验数据,应用训练的分类规则,对采样到的与航空影像具有相同成像光谱范围和光谱分辨率的地面光谱数据计算投影向量,实现不同地物信息的提取。
2.1 数据采集与预处理
野外光谱实验中,对实验区范围内有代表性的不同地物进行选择,并确定具有群体特性的植株作为观测目标。测量中以10个光谱为一个采样组,计算3组测量值的平均值作为该点的光谱反射率,并采用S-G去噪方法[26],对平均后的光谱曲线进行去噪处理,获取实测地物光谱曲线。实验区主要地物类型有灌木、冬青、柳树、龙爪槐、柿子树、伏地猪草、玉米、红叶林、桃树、道路以及土壤等。
同步开展PHI高光谱影像获取,PHI航摄仪成像波段为450 ~ 1 000 nm,共计128个波段,成像相对高度为1 500 m,影像空间分辨率为1.5 m。其数据预处理主要包括光谱重采样、大气校正以及影像去噪。通过设定波长信息与波谱半幅值全宽(full width half maximum,FWNM),应用FWHM间隔的高斯模型,对光谱进行重采样,将地面样本光谱重采样到与影像相同的波谱范围与光谱分辨率。由于航空影像飞行高度较低,受大气影响较小,采用基于统计的大气校正方法,设真实值X与辐亮度值Y之间存在的函数关系为
式中k和b为转换系数。根据实地采集光谱数据(主要为土壤、柏油路和水泥路),计算地面反射率,便可计算k和b,再反求其他点校正后的光谱曲线。最后应用S-G去噪方法对高光谱影像进行去噪处理。预处理后影像如图1所示。
图1
2.2 特征光谱选择
运用“类内密度最大,类间距离最大”指标的波段可分性与重要性分析方法,开展地面实测植被光谱的自动特征选择。局部散点密度的窗宽参数R依距离取2个波段最大距离的1.5倍,并与PHI高光谱影像谱段范围一致,分析350~1 000 nm范围上的光谱重要性,如图2所示。
图2
图2
植被谱段范围350~1 000 nm的波谱重要性分析
Fig.2
Spectral importance analysis of spectral range 350~1 000 nm
结合植被波谱重要性分析,计算其导数值并归一化处理后,确定极值点对应的特征谱段。
2)另一较大可分性波段在520 nm处,区间范围为500~540 nm处,该谱段亦是土壤、道路与植被可分性特征,该范围可分性主要由叶绿素的反射产生,可有效区分植被中的非绿色植被。
3)570 nm,610 nm,660 nm,810 nm,890 nm和930 nm范围均为较明显的分类特征谱段。570 nm附近光谱曲线上绿光波段反射峰右侧有一急速下滑趋势,红叶林的该特征不太明显,其值基本为0,与其他植被特性相差较远,其中伏地猪草在该处为负值,最为明显。
应用“类内密度最大,类间距离最大”指标,自动计算波谱指标重要性,且结合PHI影像成像光谱范围以及信噪比较高谱段,确定8个特征波段范围区间,即420 nm附近、520 nm附近、570 nm附近、610 nm附近、660 nm附近、690 nm附近、715 nm附近及810 nm附近。波段重要性分析的特征谱段与光谱一阶微分结果的特征谱段高度吻合。
2.3 基于实测光谱数据分类规则建立
2.3.1 全谱段整体最优分类分析
在确定的8个谱段范围上对11个细分地类进行投影分析,结果见表1。
表1 实测光谱8个谱段投影结果
Tab.1
细分地类 | 投影均值 | 细分地类 | 投影均值 | |
---|---|---|---|---|
冬青 | 0.603 | 玉米 | 0.744 | |
柳树 | 0.760 | 红叶林 | 0.663 | |
龙爪槐 | 0.688 | 桃树 | 0.896 | |
柿子树 | 0.714 | 土壤 | 0.087 | |
伏地猪草 | 1.232 | 道路 | 0.014 | |
灌木 | 0.844 |
土壤和道路投影均值最低,而其他地类投影均值相对于土壤和道路较明显,呈数量级差距,可优先区分植被与非植被; 但不同种类植被间可分性仍不强,需在全谱段整体最优识别基础上,继续围绕光谱重要性分析,开展不同特征谱段上的局部最优投影,实现植被的精细化分类。2.3.2 分谱段局部最优投影分析以不同特征谱段为局部波谱分析单元,分别在不同谱段上对不同植被进行投影计算,分析投影后地物可识别的最佳投影谱段。并在此基础上,通过逐次寻优,实现地物最佳识别的波谱范围提取,逐次剔除已识别植被,分谱段、分层次循环投影实现不同植被识别。2.3.3 基于实测光谱数据建立地物分类规则基于以上分析,构建结合投影寻踪方法的实验区地物分类规则,为PHI高光谱影像分类提供依据,如图3所示。
图3
图3
基于实测光谱数据构建地物分类规则
Fig.3
Classification rules based on measured spectral data
2.4 结合地面光谱的PHI航空影像分类
将地面光谱数据重采样到与航空影像具有相同光谱成像范围及光谱分辨率,再运用由实测数据构建的分类规则(图3),对重采样后的地面光谱数据按相应谱段计算投影,并应用此投影方向,计算PHI高光谱影像在此投影上的结果,以此逐次提取分类后地物。分类思路为先应用典型样本训练规则,再应用规则指导分类地物,其数学意义在于,在大样本统计下应用少量的典型样本可以估计整体分布。图4(a)为全部8个投影方向实现不同种类地物识别的结果,即植被与非植被(道路,土壤)分类提取结果,图4(b)为计算得到的归一化植被指数(normalized difference vegetation index,NDVI)图像。图4中2个分类结果整体相似,植被在影像上表现出亮色调,而道路和土壤明显得到抑制,整个影像中表现为明显的暗色调。因投影寻踪分析中约束条件的影响,投影分类的结果层次更强,提取的植被明显呈现出层次感(亮度等级不同),植被与土壤、道路等反差增大(投影分析结果值域无固定范围,且明显较NDVI图像值域范围大),可分性更强。投影寻踪分析中植被个体特征更加明显,较之于NDVI提取的结果,能有效地避免太阳阴影的影响,NDVI图像中植被多为成片分布,而投影寻踪结果能更好地区分不同物种,植被轮廓更加清晰,纹理特征更加丰富,例如高大茂密的植被特征明显突出。基于改进的投影寻踪分析方法,PHI航空影像分类结果如图5所示。
图4
图4
8个谱段范围投影图像及NDVI图像
Fig.4
Projection image of eight spectral range and NDVI image
图5
以实地采样样本类型与位置为依据,对提取结果开展定量分析,分类结果如表2所示。在11种地类中,有5种地物分类准确率高达90%以上,3种地物分类准确率达80%以上,分类总体精度达85%以上,算法分类结果准确度高。
表2 分类结果精度分析
Tab.2
序号 | 样本情况 | 分类情况 | |||
---|---|---|---|---|---|
地物类型 | 样本数 量/个 | 正确分类 数量/个 | 错误分类 数量/个 | 准确 率/% | |
1 | 冬青 | 9 | 8 | 1 | 88.89 |
2 | 玉米 | 12 | 10 | 2 | 83.33 |
3 | 柳树苗 | 9 | 7 | 2 | 77.78 |
4 | 红叶林 | 3 | 3 | 0 | 100 |
5 | 龙爪槐 | 6 | 5 | 1 | 83.33 |
6 | 桃树 | 3 | 2 | 1 | 66.67 |
7 | 柿子树 | 9 | 7 | 2 | 77.78 |
8 | 伏地猪草 | 3 | 3 | 0 | 100 |
9 | 灌木 | 12 | 11 | 1 | 91.67 |
10 | 道路 | 10 | 10 | 0 | 100 |
11 | 土壤 | 12 | 11 | 1 | 91.67 |
3 结论
围绕“地-空”一体的高光谱遥感成像应用,以特征选择和精细化分类为目的,讨论了基于实测光谱的PHI高光谱影像不同地物类型提取方法,结论如下:
1)应用“类内密度最大,类间距离最大”指标,结合高光谱PHI影像成像范围与信噪比较高波段,自动计算不同波段在地物识别中的可分性和重要性,确定研究区内作物在420 nm,520 nm,570 nm,610 nm,660 nm,690 nm,715 nm及810 nm 等8个谱段附近具备强可分性。该结论与地物光谱一阶微分结果以及同类研究高度吻合。
2)采用改进的投影寻踪分类方法进行PHI高光谱影像分类,研究区分类精度达85%以上,识别效果同NDVI相比,具有层次感强、植被轮廓更加清晰、有效避免阴影影响的特点,适用于不同植被精细分类应用。
但是,改进的投影寻踪分析方法的应用需以知道场景内基本地物分布为前提,若假设前提不能完全满足时,如何有效运用该方法仍有待进一步探讨。
参考文献
中国高光谱遥感的前沿进展
[J].
Current progress of hyperspectral remote sensing in China
[J].
我国高光谱遥感的发展
[N].
The development of hyperspectral remote sensing in China
[N].
Representative multiple kernel learning for classification in hyperspectral imagery
[J].
Combining support vector machines and Markov random fields in an intergrated framework for contextual image classification
[J].
高光谱影像地物智能化提取理论与方法研究
[D].
Research on the Theory and Method for Intelligent Extraction of Objects in Hyperspectral Images
[D].
Feature mining for hyperspectral image classification
[J].
高光谱遥感数据特征提取与特征选择方法研究
[D].
Research on Feature Extraction and Feature Selection for Hyperspectral Remote Sensing Data
[D].
联合丰度信息与空谱特征的高光谱影像分类研究
[D].
Hyperspectral Image Classification with Abundance Information and Spatial Spectrum Characteristics
[D].
基于信息熵的自训练半监督高光谱遥感影像分类研究
[D].
Research on Self-Training Semi-Supervised Classification for Hyperspectral Remote Sensing Image Based on Information Entropy
[D].
高光谱遥感影像空谱特征提取与分类方法研究
[D].
Research on Spectral-Spatial Feature Extraction and Classification Methods for Hyperspectral Remote Sensing Imagery
[D].
Simultaneous sparse graph embedding for hyperspectral image classification
[J].
Multiple classifier systems in remote sensing:From basics to recent developments
[C]//
基于机载PHI高光谱数据的森林优势树种分类研究
[J].
Classification of forest species using airborne PHI hyperspectral data
[J].
高光谱遥感影像分类研究进展
[J].
Review of hyperspectral remote sensing image classification
[J].
高光谱图像处理与信息提取前沿
[J].
Advancement of hyperspectral image processing and information extraction
[J].
高光谱遥感数据最佳波段选择方法研究
[D].
Study on the Optimum Band Selection Method for Hyperspectral Remote Sensing Data
[D].
基于类内类间距离的模糊聚类算法及图像分割应用
[D].
Fuzzy Clustering Algorithm Based on Intra-Class Distance and Its Application in Image Segmentation
[D].
基于类内类间离散度的分类器设计方法
[J].
Classifier design based on between-class and within-class scatter
[J].
基于类内类间距离的模糊C-均值聚类分割算法
[J].
Fuzzy C-means clustering segmentation algorithm based on intra-class and inter-class distance
[J].
Review of image segmentation techniques:Including pre and post processing operations
[J].
基于投影寻踪学习网络算法的植物群落高分遥感分类研究
[J].
DOI:10.3724/SP.J.1047.2016.00124
URL
Magsci
[本文引用: 1]
<p>传统的植物群落调查方法主要是野外样地调查和抽样统计,其对于地形复杂的区域难以做到对数据的全面调查;将遥感技术应用于植物群落调查,可实现数据的全面获取,以及对植物群落的快速分类。在深圳市植物群落野外样地调查的基础上,本文应用高分辨率Pléiades影像,结合光谱、地形及纹理信息,采用投影寻踪学习网络的方法,实现了深圳市东部地区植物分类。在实验中,选取人工林和次生林中典型群落样本,将投影寻踪与学习网络算法结合应用于植被分类,通过分类结果与经典监督分类方法比较表明,该算法应用于植物群落分类是可行的;并且该算法分类精度高,更新速度快,能满足深圳市重点项目基本生态控制线专项调查的要求。</p>
Research on classification of plant community using projection pursuit learning network algorithm on high resolution remote sensing images
[J].
上海城市植被光谱反射特征分析
[J].地物光谱特征是遥感机理的重要内容, 也是遥感应用研究的重要依据。 城市植被作为城市中重要的生命支持系统, 对于维持城市生态系统平衡起着核心的作用。 文章以上海城市植被为研究对象, 利用美国ASD公司生产的FieldSpec3便携式地物波谱仪, 通过对2009年10月金边黄杨, 金丝桃, 圆柏, 麦冬, 珊瑚树, 小叶黄杨等野外和室内光谱数据采集和处理, 获得了这6种城市植被的高光谱数据。 采用Savitzky Golay平滑方法进行光谱噪声去除和利用微分光谱技术进行植被环境背景影响消除, 分别从植被冠层光谱和叶片光谱以及微分光谱, 对这6种城市植被的光谱特征以及光谱曲线的差异和变化规律进行分析。 可为今后进一步研究城市植被理化性能、 参数反演、 植被分类、 植被调查以及城市环境监测等遥感应用服务。
Analysis on urban vegetations reflectance characteristics in Shanghai
[J].
基于PHI高光谱影像的植被光谱特征应用研究
[J].
Application of vegetation spectral features based on PHI imaging data
[J].
/
〈 |
|
〉 |
