CN111242224A - 一种基于无人机提取分类样本点的多源遥感数据分类方法 - Google Patents
一种基于无人机提取分类样本点的多源遥感数据分类方法 Download PDFInfo
- Publication number
- CN111242224A CN111242224A CN202010046232.6A CN202010046232A CN111242224A CN 111242224 A CN111242224 A CN 111242224A CN 202010046232 A CN202010046232 A CN 202010046232A CN 111242224 A CN111242224 A CN 111242224A
- Authority
- CN
- China
- Prior art keywords
- classification
- sample points
- remote sensing
- classified
- data set
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于无人机提取分类样本点的多源遥感数据分类方法,包括:从无人机航拍相片中均匀提取分类样本点,并对每类样本点进行准备标定;获取分类遥感数据集,并对遥感数据集进行影像处理,根据分类遥感影像数据集对分类样本点进行地理空间定位;分类遥感数据集包括:微波数据Sentinel‑1数据集、多光谱Sentinel‑2数据集、基于Sentinel‑2数据集的植被指数数据集、数字高程模型数据集;通过地理空间信息定位后的分类样本点,利用随机森林分类模型,获得分类结果。本发明基于无人机提取分类样本点的多源遥感数据随机森林分类方法,可以快速、有效、廉价地实现地表类型分类制图过程;同时剔除边缘分类样本点的影响后,分类精度明显提高,特别是kappa系数的精度更优。
Description
技术领域
本发明涉及遥感数据分类技术领域,更具体的涉及一种基于无人机提取分类样本点的多源遥感数据分类方法。
背景技术
全球喀斯特地貌面积较大,相当一部分的全球人口水源依赖于喀斯特区域的含水层。喀斯特生态系统十分脆弱,特别容易受到环境变化的侵袭,导致区域内地表植被发生破坏,进而造成其地表景观退化为裸土区域,甚至退化为岩石区域,而这种石漠化现象又是一种严重的生态系统短期不可逆过程。在我国西南喀斯特地区石漠化面积较大,其中作为岩溶中心的贵州省,其表层土壤在1974年到2001年间退化为石漠化区域的速度较快,不过这种趋势在最近20年开始转为良性,许多地区的植被开始变得比以前更绿。尽管如此,对喀斯特区域,特别是位于岩溶中心贵州省的长期监测依然不容忽视。
随着多源遥感数据的发展,遥感影像在时空分辨率和光谱分辨率方面都有极大提升,特别是针对喀斯特区域的植被动态和地物类型监测研究越来越成熟。现有的地表类型分类方法日益精准,但是作为任意一种分类模型必要输入条件的野外实测分类样本点获取较为困难。特别是在较大尺度范围内,如果通过传统野外调查法收集分类样本点,人力、物力和时间成本的花费极高,严重阻碍大范围地表分类研究的发展。
发明内容
本发明实施例提供一种基于无人机提取分类样本点的多源遥感数据分类方法,用以解决上述背景技术中提出的问题。
本发明实施例提供一种基于无人机提取分类样本点的多源遥感数据分类方法,包括:
从无人机航拍相片中均匀提取分类样本点,并对每类样本点进行准备标定;其中,准备标定的样本点类型包括:农田和草地、林地和灌木、空地和裸地、道路、建筑物;
获取分类遥感数据集,所述分类遥感数据集包括:微波数据Sentinel-1数据集、多光谱Sentinel-2数据集、基于Sentinel-2数据集的植被指数数据集、数字高程模型数据集;
对遥感数据集进行处理,获得分类遥感影像数据集;并根据分类遥感影像数据集对分类样本点进行地理空间定位;
通过地理空间信息定位后的分类样本点,利用随机森林分类模型,获得分类结果。
进一步地,所述从无人机航拍相片中提取分类样本点,包括:
通过目视解译方法,从无人机航拍相片影像中均匀提取分类样本点。
进一步地,所述从无人机航拍相片中提取分类样本点,包括:
剔除不同地表类型边缘处的样本点。
进一步地,基于10m分辨率,采用SNAP软件对Sentinel-1数据集进行轨道修正、热噪声去除、辐射校正、散斑滤波、距离-多普勒地形校正处理,获得VV极化影像数据集和VH极化影像数据集。
进一步地,所述Sentinel-2数据集包含13个波段数据,涵盖可见光、近红外和短波红外光谱波段;采用Sen2Cor软件对Sentinel-2数据集进行地形校正、大气校正和辐射校正处理,获得除第10波段外的12层影像数据集,并对12层影像数据集重采样至10m分辨率。
进一步地,所述植被指数数据集包括:NDVI、EVI、SAVI,计算公式如下所示:
NDVI=(NIR–Red)/(NIR+Red)
EVI=2.5×(NIR-Red)/(NIR+6.0Red–7.5Blue+1)
SAVI=(NIR-Red)(1+L)/(NIR+Red+L)
式中,NIR、Red和Blue分别对应近红外、红波段和蓝波段的数据;L为土壤调节系数,由实际区域条件确定;NIR、Red和Blue波段的数据分别对应Sentinel-2数据集的第8波段、第4波段和第2波段的数据。
进一步地,所述土壤调节系数L=0.5。
进一步地,所述DEM数据集采用SRTM DEM数据集,对SRTM DEM数据集重采样至10m分辨率后,获取高程DEM影像数据集、坡度slope影像数据集、坡向aspect影像数据集、剖面曲率profile curvature影像数据集。
进一步地,所述随机森林分类模型,包括:
利用readOGR()和brick()命令读取分类样本点影像和分类遥感数据集于R语言环境中;
利用如下代码搭建随机森林分类模型;
rf<-randomForest(lc~b1+b2+b3+b4+b5+b6+b7+b8+b9+b8a+b11+b12,
data=rois,
ntree=500,
importance=TRUE)
其中,b1~b12为随机森林分类模型中的参数层影像,不同数据集对应不同的参数层影像;
利用tuneRF()和randomForest()命令完成随机森林分类模型的调参训练;
利用writeRaster()命令对分类结果进行出图,生成分类结果影像。
进一步地,所述地表类型分类结果的准确率指标包括:总体准确率OA和Kappa系数,计算公式如下:
OA=(TP+TN)/(TP+FN+FP+TN)
式中,TP为真正,即被随机森林分类模型分类正确的正样本;FN为假负,即被随机森林分类模型分类错误的正样本;FP为假正,即被随机森林分类模型分类错误的负样本;TN为真负,即被随机森林分类模型分类正确的负样本;OA为总体分类精度,即分类正确的样本个数占所有样本个数的比例。
Kappa=(Po-Pe)/(1-Pe)
式中,Po为对角线单元中观测值的总和,即总体分类精度OA;Pe为对角线单元中期望值的总和;Kappa为评价一致性的测量值,表示分类与完全随机的分类产生错误减少的比例。
本发明实施例提供一种基于无人机提取分类样本点的多源遥感数据分类方法,与现有技术相比,其有益效果如下:
本发明中基于无人机提取分类样本点的多源遥感数据随机森林分类方法,可以快速、有效、廉价地实现地表类型分类制图过程,同时也可为今后上万、上百万级的海量样本点提取过程提供技术支持和方法基础。但是涉及混合像元的分类样本点存在问题,在剔除边缘分类样本点(混合像元)的影响后,分类精度明显提高,特别是kappa系数的精度更优。因此,在以后的相关研究中进行布点时,尽量避免提取边缘处的分类样本点,应选取地物类型均匀一致的区域提取样本点。该方法能够有效区分枯萎植被和裸露土地,即使仅使用可见光波段组合生成的影像,也可以参考无人机影像很方便的分辨各类地表类型;可以扩大分类样本点的采集时间,不只局限于植物生长最旺季(如7到9月);该方法也减少对植物生长季运算过程的时间消耗,不需利用长时间序列的植被研究数据反演整个植被物候过程来完成分类。
附图说明
图1为本发明实施例提供的随机森林分类结果;
图2a为本发明实施例提供的S2数据集的分类结果混淆矩阵图;
图2b为本发明实施例提供的S2&VI数据集的分类结果混淆矩阵图;
图2c为本发明实施例提供的S2&VI&DEM数据集的分类结果混淆矩阵图;
图2d为本发明实施例提供的S2&VI&S1数据集的分类结果混淆矩阵图;
图2e为本发明实施例提供的b3&b2&b4&b6数据集的分类结果混淆矩阵图;
图3a为本发明实施例提供的S2数据集的基尼指数图;
图3b为本发明实施例提供的S2&VI数据集的基尼指数图;
图3c为本发明实施例提供的S2&VI&DEM数据集的基尼指数图;
图3d为本发明实施例提供的S2&VI&S1数据集的基尼指数图;
图3e为本发明实施例提供的b3&b2&b4&b6数据集的基尼指数图;
图4a为本发明实施例提供的数据集S2数据集的随机森林分类结果和混淆矩阵;
图4b为本发明实施例提供的数据集S2&VI数据集的随机森林分类结果和混淆矩阵;
图4c为本发明实施例提供的数据集S2&VI&S1数据集的随机森林分类结果和混淆矩阵;
图5为本发明实施例提供的枯萎植被和裸地区分图;
图6为本发明实施例提供的细小、杂乱斑块图;
图7为本发明实施例提供的一种基于无人机提取分类样本点的多源遥感数据分类方法流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图7,本发明实施例提供一种基于无人机提取分类样本点的多源遥感数据分类方法,该方法包括:
步骤S1,从无人机航拍相片中均匀提取分类样本点,并对每类样本点进行准备标定;其中,准备标定的样本点类型包括:农田和草地、林地和灌木、空地和裸地、道路、建筑物。
步骤S2,获取分类遥感数据集,分类遥感数据集包括:微波数据Sentinel-1数据集、多光谱Sentinel-2数据集、基于Sentinel-2数据集的植被指数数据集、数字高程模型数据集。
步骤S3,对遥感数据集进行处理,获得分类遥感影像数据集;并根据分类遥感影像数据集对分类样本点进行地理空间定位。
步骤S4,通过地理空间信息定位后的分类样本点,利用随机森林分类模型,获得分类结果。
对于上述步骤S1~S4的具体说明如下:
详细的地表类型分类及监测是我国西南喀斯特区域关键的防治石漠化手段。目前多种多样的遥感数据广泛应用于地物分类制图研究中,但是稀少的野外实测样本点一直以来都是准确、有效感知地表类型的技术瓶颈之一。因此,本发明利用无人机航拍相片提取大量野外实测样本点的方法,完成地表类型分类过程,为后期的廉价海量分类样本点提取提供切实可行的方法。
本发明中均匀提取分布于研究区内的982个分类样本点,利用包括Sentinel-1/2(S1/2)数据集、S2数据计算获得的植被指数(Vegetation index,VI)数据集和数字高程模型(Digital Elevation Model,DEM)数据集的遥感影像资料,然后借助随机森林分类模型完成研究区地表类型分类。上述遥感数据集不仅涵盖可见光数据,还包括近红外、短波红外和微波光谱波段数据资料。分类后的结果显示,除包含DEM数据集的分类结果外(总体精度和Kappa系数分别为74.54%和61.73%),其他4种数据集组合(仅用S2数据集、S2&VI数据集、S2&VI&DEM数据集,以及4个S2波段b3&b2&b4&b6组成的数据集)的制图总体精度(Overall Accuracy,OA)和Kappa系数精度都在75%和65%以上。此外,在不考虑位于边缘的样本点(即混合像元)后,分类制图结果具有更好的鲁棒,3种分类精度最高的数据集(仅用S2数据集、S2&VI数据集和S2&VI&DEM数据集)OA和Kappa系数精度都提高至85%和79%以上。特别是kappa系数的精度更优,提高了近15%。以上研究结果能够为喀斯特地区的地表类型制图提供准确而有效的技术手段和方法支持。
另外,本发明还对分类制图的边缘分类样本点(混合像元)、枯萎植被和裸露土地的有效区分能力,并着重强调航拍影像分类采样点自动化过程实现的必要性,展望其将会成为未来研究的一个热点方向。
需要说明的是,现有的遥感数据源种类繁多,分类方法也各式各样。为充分展示利用无人机提取海量、廉价分类样本点的可能性,选取开源数据集作为分类的遥感数据源,如Sentinel-1/2和SRTM DEM。上述数据的优势是提供免费的,涵盖可见光、近红外、短波红外和微波范围的影像资料,为分类方法研究提供了广泛的光谱数据。同时,分类方法选用对遥感数据分类具有较优精度的随机森林,相比其他机器学习分类方法(除硬件要求极高的深度学习外),随机森林具有更好的鲁棒。
研究区概括
本发明实施例的研究区地处贵州省威宁县北部,位于104.100°E-104.118°E和27.179°N-27.191°N之间,面积1.7km×1.4km,采集无人机照片120张,提取分类采样点982个。
UAV航拍相片及分类样本点数据获取
无人机航拍相片采集于2018年4月21日,使用大疆精灵4拍摄,共采集120张,飞行高度近300米,相片分辨率为1200万像素,使用FragMAP软件完成操作拍摄过程。采样点通过ArcGIS软件,设置均匀分布于研究区内的矢量点完成,共计982个。如果利用传统的生态样方框完成本实验中近千个点位、甚至后续研究中上万、上百万点位的植被调查,所需的人工和时间成本花费极大,因此本发明利用无人机影像完成分类样本点数据的提取过程。
仅用开源遥感影像来目视解译,即使分辨率已达10m,仍难以清晰分辨各类地表类型。举例来说,研究区内某些裸露土地的颜色与生长不茂盛的树林颜色极为相近,十分容易解译错误为林地,但是如果有无人机航拍数据作为参考,地表类型就不会解译错误。所有分类样本点都是通过目视解译,将无人机影像定位至Sentinel-2影像上,然后完成样本点提取过程。
遥感数据源
Sentinel-1/2数据来源于欧空局(http://scihub.copernicus.ed/)。Sentinel-2(S2)数据共包含13个波段数据,涵盖可见光、近红外和短波红外光谱波段,其中有5个近红外波段的数据可应用于植被相关的研究。经过Sen2Cor软件处理后,完成地形校正、大气校正和辐射校正等基础影像处理过程,最后得到除第10波段外的12层影像数据集,全部重采样至10m分辨率与第2(蓝)、3(绿)、4(红)和8(近红外)波段分辨率一致。Sentinel-1(S1)GRD数据(C波段,VV和VH极化)分辨率也为10m,使用SNAP软件进行轨道修正、热噪声去除、辐射校正、散斑滤波和距离-多普勒地形校正操作后得到VV和VH2层数据影像。本发明所使用的S2和S1数据分别获取于2018年4月17日和2018年4月20日。DEM数据使用SRTM DEM,重采样至10m分辨率后计算获取高程(DEM)、坡度(slope)、坡向(aspect)和剖面曲率(profilecurvature)4层数据影像。经过处理,上述所有遥感影像和无人机航拍数据全部使用WGS_1984_UTM_Zone_48N投影。
植被指数的计算
为提高分类精度,引入植被指数数据(Vegetation indices,VI),主要包括NDVI(Normalized Difference Vegetation Index)EVI(Enhanced Vegetation Index)和SAVI(Soil-Adjusted Vegetation Index),其计算公式如下所示:
NDVI=(NIR–Red)/(NIR+Red) (1)
EVI=2.5×(NIR-Red)/(NIR+6.0Red–7.5Blue+1) (2)
SAVI=(NIR-Red)(1+L)/(NIR+Red+L) (3)
式中NIR、Red和Blue分别对应近红外、红和蓝波段的数值,L是土壤调节系数,由实际区域条件确定,一般情况都采用L=0.5来完成运算。其中,NIR、Red和Blue波段数据分别对应S2数据的第8、4和2波段结果。
随机森林分类模型的构建
随机森林是一种组成式的监督分类法,以决策树为基础,实现对多决策树的集成。随机森林方法在遥感数据分类研究中应用较为广泛,该分类模型从原始的训练数据集中采取有放回的抽样(bagging)方法完成子数据集的构造过程。该过程中不同子数据集的元素可以重复,同一子数据集中的元素也可以重复。同时,因其引入两个随机性属性(样本随机,特征随机),所以分类结果不容易陷入过拟合。随机森林特征重要性的大小正相关于该特征对森林中每棵树的贡献大小,当平均该特征对每个树的贡献之后,得到基尼指数(Giniindex)。此外,袋外数据(Out Of Bag,OOB)错误率可以作为评估指标来衡量特征集贡献大小,通常优先选择袋外误差率最低的特征集。
随机森林分类模型通过R语言实现。
首先,需要加载的工具包包括randomForest、raster、rgdal、lattice、ggplot2、caret和e1071。
然后利用readOGR()和brick()命令读取分类样本点影像和用于分类的遥感基础数据集于R语言环境中。
第三步,利用下述代码搭建随机森林模型。
rf<-randomForest(lc~b1+b2+b3+b4+b5+b6+b7+b8+b9+b8a+b11+b12,
data=rois,
ntree=500,
importance=TRUE)
其中,b1到b12为本次分类模型中的参数层影像,不同数据集对应于不同的参数层影像。如本示例针对S2数据集,对应权利5所述Sen2Cor软件对其地形校正、大气校正和辐射校正处理后的12个参数层影像。而S2&VI不仅包括S2数据集,还包括权利6所示的NDVI、EVI、SAVI3个植被参数层影像。S2&VI&DEM不仅包括S2和权利6中的3个植被参数层影像,还包括权利8所示的高程、坡度、坡向和剖面曲率4个地形参数影像。S2&VI&S1包括S2和权利6中的3个植被参数层影像,还包括权利4所述的VV和VH2个极化参数层影像。而b3&b2&b4&b6仅包括S2数据集的第3、2、4和6这4个波段参数层影像。
第四步,利用tuneRF()和randomForest()命令完成模型的调参训练。
最后,利用writeRaster()命令对分类结果进行出图,生成分类结果影像。
本发明实施例中,主要区分5种地表类型,如表1所示。
表1实验研究区地表类型分类样本点概况
验证分类结果的准确率指标主要包括总体准确率(Overall Accuracy,OA)和Kappa系数两种,其计算公式如下:
OA=(TP+TN)/(TP+FN+FP+TN) (4)
式中,TP为真正,即被模型分类正确的正样本;FN为假负,即被模型分类错误的正样本;FP为假正,即被模型分类错误的负样本;TN为真负,即被模型分类正确的负样本;OA为总体分类精度,即分类正确的样本个数占所有样本个数的比例。
Kappa=(Po-Pe)/(1-Pe) (5)
式中,Po为对角线单元中观测值的总和,也就是总体分类精度OA;Pe为对角线单元中期望值的总和;Kappa为评价一致性的测量值,表示分类与完全随机的分类产生错误减少的比例。
依据上述数据和随机森林分类模型完成的制图结果如图1所示,除S2&VI&DEM数据集的分类结果较差外(表2),其余分类结果都有相对稳定的结果,OA在75%以上,Kpppa系数在65%以上,同时各分类结果的空间分布位置基本一致。其中,2个数据集S2&VI的分类精度最高,但是OOB的数值以3个数据集S2&VI&S1的分类结果为最高。本发明的研究结果基本呈现数据资料越多,分类精度越高的规律。不过,本发明中DEM及其相关计算结果的引入会明显产生噪声,导致分类结果降低。因此,分类研究也应注意筛选变量,以免冗余变量的引入产生误差,导致分类精度降低。
表2分类结果精度评价指标表
根据图2a~2e可以发现,分类样本点越多的地表类型,分类精度越高;分类样本点越少的地表类型,最终分类精度越低。特别是具有最少采样本点的建筑物分类,在测试集的5种分类结果中都没能正确区分。
图3a~3e展示了各分类数据集中,不同遥感图层的基尼指数值。在前4种数据集中,S2数据的第2、3、4和6波段都具有较高的基尼指数数据值,因此仅用上述4层数据波段完成分类后发现,其分类后的空间分布形式基本同其他多数据图层分类结果相似(图1)。特别是仅用以上4层数据波段完成的分类结果精度略低于最高分类数据集的结果精度,远高于引入冗余数据集(DEM)后的分类精度(表2和图2a~2e)。以上结果也说明,数据降维的操作方式,在大数据量处理时会有效的节约时间,同时又保证较高的分类精度。
边缘分类样本点(混合像元)对分类结果准确率的影响
本发明中分类样本点是按照均匀分布规则设置的,因此有接近1/3的样本点位于不同地表类型的边缘处(即混合像元),如表1所示。在此,通过剔除表1中318个边缘样本点,使所有样本点只是位于一致性极高区域后的分类精度。将位于一致性较高地表类型的664个样本点输入分类模型中分类结果精度最高的3个数据集进行分类,结果如图4a~4c所示,其分类精度存在明显提升(表3),OA和Kappa都在85%和79%以上。因此,剔除边缘分类样本点(混合像元)对分类精度的提高十分重要,本发明中对OA和Kappa的提升,一个接近10%,一个接近15%。
表3数据集S2、S2&VI和S2&VI&S1分类结果精度评价指标表
基于无人机提取分类样本点的分类方法能够有效识别枯萎植被
常规目视解译,特别是没有无人机航拍相片作为参考资料时,可见光遥感影像目视解译对枯萎林地和某些裸地的区分相对困难,如图5A所示从S2的合成图像上判读枯萎林地和裸地的颜色十分接近。不过,本发明利用无人机提取分类样本点完成分类后,对枯萎林地和裸地的区分会准确很多,如图5B所示。特别是图5B从上到下第二个圈中对枯萎林地的分类,基本可以有效地衔接旁边的旺盛林地。其中,图5A中从上到下圈内地表类型依次为森林、森林、裸地;图5B中从上到下圈内地表类型依次为森林、森林、裸地。
细小斑块和杂乱斑块
图6展示了分类对细小、杂乱斑块的区分效果。图6中地表类型为森林的圈中的分类效果基本可以满足常规地表类型分类需求,不过该结果也存在部分椒盐效应现象,而且部分道路也存在不连贯现象。此外,需要着重关注的是对建筑物区分效果很差,如图6红色圆圈中的房屋完全没有区分出来,在图5中更多的建筑物中也仅区分出部分样本。造成这种现象的原因主要有两方面:1)分类样本点过少,如表1所示,建筑物的分类样本点仅有4个,而且全部是边缘样本点;2)建筑物面积小,多数个体难以覆盖一个2×2的像素单元。根据实际分类制图需求,如果确实需要对这类个数少、面积小的地表类型进行分类,建议人为加大采样样本点,而不单纯依靠均匀布点所设置的样本点。此外,如果条件允许,考虑使用更高分辨率的遥感影像(一般非开源)进行分类制图,加强对小面积样本的识别能力。其中,图6A为S2合成图像上的细小斑块和杂乱斑块;图6A为本发明利用无人机提取分类样本点完成分类后的细小斑块和杂乱斑块;图6A中从上到下圈内地表类型依次为森林、森林、裸地;图6B中从上到下圈内地表类型依次为森林、森林、裸地。
以上公开的仅为本发明的几个具体实施例,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,包括:
从无人机航拍相片中均匀提取分类样本点,并对每类样本点进行准备标定;其中,准备标定的样本点类型包括:农田和草地、林地和灌木、空地和裸地、道路、建筑物;
获取分类遥感数据集,所述分类遥感数据集包括:微波数据Sentinel-1数据集、多光谱Sentinel-2数据集、基于Sentinel-2数据集的植被指数数据集、数字高程模型数据集;
对遥感数据集进行处理,获得分类遥感影像数据集;并根据分类遥感影像数据集对分类样本点进行地理空间定位;
通过地理空间信息定位后的分类样本点,利用随机森林分类模型,获得分类结果。
2.如权利要求1所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述从无人机航拍相片中提取分类样本点,包括:
通过目视解译方法,从无人机航拍相片影像中均匀提取分类样本点。
3.如权利要求1或2所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述从无人机航拍相片中提取分类样本点,包括:
剔除不同地表类型边缘处的样本点。
4.如权利要求1所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,基于10m分辨率,采用SNAP软件对微波数据Sentinel-1数据集进行轨道修正、热噪声去除、辐射校正、散斑滤波、距离-多普勒地形校正处理,获得VV极化影像数据集和VH极化影像数据集。
5.如权利要求1所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述多光谱Sentinel-2数据集包含13个波段数据,涵盖可见光、近红外和短波红外光谱波段;采用Sen2Cor软件对多光谱Sentinel-2数据集进行地形校正、大气校正和辐射校正处理,获得除第10波段外的12层影像数据集,并对12层影像数据集重采样至10m分辨率。
6.如权利要求1或5所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述植被指数数据集包括:NDVI、EVI、SAVI,计算公式如下所示:
NDVI=(NIR–Red)/(NIR+Red)
EVI=2.5×(NIR-Red)/(NIR+6.0Red–7.5Blue+1)
SAVI=(NIR-Red)(1+L)/(NIR+Red+L)
式中,NIR、Red和Blue分别对应近红外、红波段和蓝波段的数据;L为土壤调节系数,由实际区域条件确定;NIR、Red和Blue波段的数据分别对应Sentinel-2数据集的第8波段、第4波段和第2波段的数据。
7.如权利要求6所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述土壤调节系数L=0.5。
8.如权利要求1所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述DEM数据集采用SRTM DEM数据集,对SRTM DEM数据集重采样至10m分辨率后,获取高程DEM影像数据集、坡度slope影像数据集、坡向aspect影像数据集、剖面曲率profilecurvature影像数据集。
9.如权利要求1所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述随机森林分类模型,包括:
利用readOGR()和brick()命令读取分类样本点影像和分类遥感数据集于R语言环境中;
利用如下代码搭建随机森林分类模型;
rf<-randomForest(lc~b1+b2+b3+b4+b5+b6+b7+b8+b9+b8a+b11+b12,
data=rois,
ntree=500,
importance=TRUE)
其中,b1~b12为随机森林分类模型中的参数层影像,不同数据集对应不同的参数层影像;
利用tuneRF()和randomForest()命令完成随机森林分类模型的调参训练;
利用writeRaster()命令对分类结果进行出图,生成分类结果影像。
10.如权利要求1或9所述的基于无人机提取分类样本点的多源遥感数据分类方法,其特征在于,所述地表类型分类结果的准确率指标包括:总体准确率OA和Kappa系数,计算公式如下:
OA=(TP+TN)/(TP+FN+FP+TN)
式中,TP为真正,即被随机森林分类模型分类正确的正样本;FN为假负,即被随机森林分类模型分类错误的正样本;FP为假正,即被随机森林分类模型分类错误的负样本;TN为真负,即被随机森林分类模型分类正确的负样本;OA为总体分类精度,即分类正确的样本个数占所有样本个数的比例;
Kappa=(Po-Pe)/(1-Pe)
式中,Po为对角线单元中观测值的总和,即总体分类精度OA;Pe为对角线单元中期望值的总和;Kappa为评价一致性的测量值,表示分类与完全随机的分类产生错误减少的比例。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010046232.6A CN111242224B (zh) | 2020-01-16 | 2020-01-16 | 一种基于无人机提取分类样本点的多源遥感数据分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010046232.6A CN111242224B (zh) | 2020-01-16 | 2020-01-16 | 一种基于无人机提取分类样本点的多源遥感数据分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111242224A true CN111242224A (zh) | 2020-06-05 |
CN111242224B CN111242224B (zh) | 2021-07-20 |
Family
ID=70869466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010046232.6A Active CN111242224B (zh) | 2020-01-16 | 2020-01-16 | 一种基于无人机提取分类样本点的多源遥感数据分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111242224B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680427A (zh) * | 2020-06-16 | 2020-09-18 | 中国气象科学研究院 | 一种极度稀疏植被区域负土壤调节因子计算方法 |
CN111898503A (zh) * | 2020-07-20 | 2020-11-06 | 中国农业科学院农业资源与农业区划研究所 | 基于云覆盖遥感影像和深度学习的作物识别方法和系统 |
CN112084991A (zh) * | 2020-09-18 | 2020-12-15 | 中国农业科学院农业资源与农业区划研究所 | 基于多源遥感时序影像和卷积神经网络的作物早期识别方法 |
CN112597855A (zh) * | 2020-12-15 | 2021-04-02 | 中国农业大学 | 农作物倒伏程度识别方法及装置 |
CN112818880A (zh) * | 2021-02-05 | 2021-05-18 | 郑州科技学院 | 基于深度学习的航拍图像植被提取与分类方法 |
CN113009485A (zh) * | 2021-03-10 | 2021-06-22 | 安徽皖南烟叶有限责任公司 | 一种基于改进植被指数的遥感烟田识别方法 |
CN113421273A (zh) * | 2021-06-30 | 2021-09-21 | 中国气象科学研究院 | 林草搭配信息遥感提取方法及装置 |
CN114494198A (zh) * | 2022-01-26 | 2022-05-13 | 自然资源部第一航测遥感院(陕西省第五测绘工程院) | 一种综合多源信息的玉米干旱损害范围提取方法 |
CN115205688A (zh) * | 2022-09-07 | 2022-10-18 | 浙江甲骨文超级码科技股份有限公司 | 茶树种植区域提取方法及系统 |
CN115965812A (zh) * | 2022-12-13 | 2023-04-14 | 桂林理工大学 | 无人机影像对湿地植被物种和地物分类的评估方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318270A (zh) * | 2014-11-21 | 2015-01-28 | 东北林业大学 | 一种基于modis时间序列数据的土地覆盖分类方法 |
CN105046673A (zh) * | 2015-07-13 | 2015-11-11 | 哈尔滨工业大学 | 基于自学习的高光谱图像和可见光图像融合分类方法 |
CN105787457A (zh) * | 2016-03-08 | 2016-07-20 | 浙江工商大学 | 一种modis卫星集成dem提高植被分类遥感精度的估算方法 |
CN107389036A (zh) * | 2017-08-02 | 2017-11-24 | 珠江水利委员会珠江水利科学研究院 | 一种结合无人机影像的大空间尺度植被覆盖度计算方法 |
CN108020211A (zh) * | 2017-12-01 | 2018-05-11 | 云南大学 | 一种无人机航空摄影估测入侵植物生物量的方法 |
US10325370B1 (en) * | 2016-05-31 | 2019-06-18 | University Of New Brunswick | Method and system of coregistration of remote sensing images |
CN110132237A (zh) * | 2019-05-05 | 2019-08-16 | 四川省地质工程勘察院 | 一种城市地面变形灾害早期识别的方法 |
CN110533052A (zh) * | 2019-09-16 | 2019-12-03 | 贵州省草业研究所 | 一种协同遥感影像的航拍相片植被信息提取方法 |
-
2020
- 2020-01-16 CN CN202010046232.6A patent/CN111242224B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318270A (zh) * | 2014-11-21 | 2015-01-28 | 东北林业大学 | 一种基于modis时间序列数据的土地覆盖分类方法 |
CN105046673A (zh) * | 2015-07-13 | 2015-11-11 | 哈尔滨工业大学 | 基于自学习的高光谱图像和可见光图像融合分类方法 |
CN105787457A (zh) * | 2016-03-08 | 2016-07-20 | 浙江工商大学 | 一种modis卫星集成dem提高植被分类遥感精度的估算方法 |
US10325370B1 (en) * | 2016-05-31 | 2019-06-18 | University Of New Brunswick | Method and system of coregistration of remote sensing images |
CN107389036A (zh) * | 2017-08-02 | 2017-11-24 | 珠江水利委员会珠江水利科学研究院 | 一种结合无人机影像的大空间尺度植被覆盖度计算方法 |
CN108020211A (zh) * | 2017-12-01 | 2018-05-11 | 云南大学 | 一种无人机航空摄影估测入侵植物生物量的方法 |
CN110132237A (zh) * | 2019-05-05 | 2019-08-16 | 四川省地质工程勘察院 | 一种城市地面变形灾害早期识别的方法 |
CN110533052A (zh) * | 2019-09-16 | 2019-12-03 | 贵州省草业研究所 | 一种协同遥感影像的航拍相片植被信息提取方法 |
Non-Patent Citations (1)
Title |
---|
张飞 等: ""基于无人机低空遥感与卫星遥感的洛宁县中药资源种植面积估算研究"", 《中国中药杂志》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680427A (zh) * | 2020-06-16 | 2020-09-18 | 中国气象科学研究院 | 一种极度稀疏植被区域负土壤调节因子计算方法 |
CN111898503A (zh) * | 2020-07-20 | 2020-11-06 | 中国农业科学院农业资源与农业区划研究所 | 基于云覆盖遥感影像和深度学习的作物识别方法和系统 |
CN111898503B (zh) * | 2020-07-20 | 2021-02-26 | 中国农业科学院农业资源与农业区划研究所 | 基于云覆盖遥感影像和深度学习的作物识别方法和系统 |
CN112084991A (zh) * | 2020-09-18 | 2020-12-15 | 中国农业科学院农业资源与农业区划研究所 | 基于多源遥感时序影像和卷积神经网络的作物早期识别方法 |
CN112597855A (zh) * | 2020-12-15 | 2021-04-02 | 中国农业大学 | 农作物倒伏程度识别方法及装置 |
CN112597855B (zh) * | 2020-12-15 | 2024-04-16 | 中国农业大学 | 农作物倒伏程度识别方法及装置 |
CN112818880B (zh) * | 2021-02-05 | 2022-09-30 | 郑州科技学院 | 基于深度学习的航拍图像植被提取与分类方法 |
CN112818880A (zh) * | 2021-02-05 | 2021-05-18 | 郑州科技学院 | 基于深度学习的航拍图像植被提取与分类方法 |
CN113009485A (zh) * | 2021-03-10 | 2021-06-22 | 安徽皖南烟叶有限责任公司 | 一种基于改进植被指数的遥感烟田识别方法 |
CN113421273A (zh) * | 2021-06-30 | 2021-09-21 | 中国气象科学研究院 | 林草搭配信息遥感提取方法及装置 |
CN114494198A (zh) * | 2022-01-26 | 2022-05-13 | 自然资源部第一航测遥感院(陕西省第五测绘工程院) | 一种综合多源信息的玉米干旱损害范围提取方法 |
CN115205688A (zh) * | 2022-09-07 | 2022-10-18 | 浙江甲骨文超级码科技股份有限公司 | 茶树种植区域提取方法及系统 |
CN115965812A (zh) * | 2022-12-13 | 2023-04-14 | 桂林理工大学 | 无人机影像对湿地植被物种和地物分类的评估方法 |
CN115965812B (zh) * | 2022-12-13 | 2024-01-19 | 桂林理工大学 | 无人机影像对湿地植被物种和地物分类的评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111242224B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111242224B (zh) | 一种基于无人机提取分类样本点的多源遥感数据分类方法 | |
AU2020101054A4 (en) | A Multi-source Remote Sensing Data Classification Method Based On the Classification Sample Points Extracted By the UAV | |
Bastin et al. | Aboveground biomass mapping of African forest mosaics using canopy texture analysis: toward a regional approach | |
Kamal et al. | Assessment of multi-resolution image data for mangrove leaf area index mapping | |
US9940514B2 (en) | Automated geospatial image mosaic generation with multiple zoom level support | |
CN105678281B (zh) | 基于光谱和纹理特征的地膜覆盖农田遥感监测方法 | |
US9042674B2 (en) | Automated geospatial image mosaic generation | |
Pasher et al. | Multivariate forest structure modelling and mapping using high resolution airborne imagery and topographic information | |
CN108932521B (zh) | 一种基于深度学习的农作物分类方法及系统 | |
CN107527014A (zh) | 县级作物种植面积遥感统计抽样调查方案设计方法 | |
CN111626269B (zh) | 一种实用的大空间范围滑坡提取方法 | |
CN113221765B (zh) | 一种基于数字相机影像有效像元的植被物候期提取方法 | |
CN113850139B (zh) | 一种基于多源遥感的森林年际物候监测方法 | |
Manakos et al. | Comparison between atmospheric correction modules on the basis of worldview-2 imagery and in situ spectroradiometric measurements | |
CN112052757B (zh) | 火烧迹地信息提取方法、装置、设备和存储介质 | |
Xu et al. | Mapping impervious surface fractions using automated Fisher transformed unmixing | |
Ma et al. | Two-step constrained nonlinear spectral mixture analysis method for mitigating the collinearity effect | |
CN113076796B (zh) | 喀斯特石漠化的遥感制图方法及装置 | |
Bektas Balcik et al. | Determination of magnitude and direction of land use/land cover changes in Terkos Water Basin, Istanbul | |
Sakieh et al. | An integrated spectral-textural approach for environmental change monitoring and assessment: analyzing the dynamics of green covers in a highly developing region | |
CN115830464A (zh) | 基于多源数据的高原山地农业大棚自动提取方法 | |
Brigante et al. | USE OF MULTISPECTRAL SENSORS WITH HIGH SPATIAL RESOLUTION FOR TERRITORIAL AND ENVIRONMENTAL ANALYSIS. | |
CN116246272A (zh) | 针对国产卫星多光谱图像质量标记的云雪区分方法 | |
CN115082812A (zh) | 一种农业景观非农生境绿地斑块提取方法及其相关设备 | |
Taha | Assessment of urbanization encroachment over Al-Monib island using fuzzy post classification comparison and urbanization metrics |
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 |