CN109212505A - 一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 - Google Patents
一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 Download PDFInfo
- Publication number
- CN109212505A CN109212505A CN201811060614.3A CN201811060614A CN109212505A CN 109212505 A CN109212505 A CN 109212505A CN 201811060614 A CN201811060614 A CN 201811060614A CN 109212505 A CN109212505 A CN 109212505A
- Authority
- CN
- China
- Prior art keywords
- point cloud
- variable
- high degree
- forest
- measured
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Electromagnetism (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于无人机多光谱高重叠度影像的森林林分特征反演方法,借助无人机机载多光谱传感器获取高重叠度影像,利用波段组合和植被指数变换获取光谱数据;并利用空中三角测量算法与点云密集匹配算法生成数字点云数据;基于归一化点云数据和光谱数据来提取待测林冠层结构特征变量;最后结合地面实测待测林林分数据和提取的待测林冠层结构特征变量构建多元回归模型,反演待测林林分特征。本发明有助于提高森林林分特征的反演精度,并可掌握待测林空间结构及动态变化规律,对于森林可持续生产经营、生态环境建模及维持区域生态平衡具有重要意义。
Description
技术领域
本发明属于森林资源管理与保护技术领域,具体涉及一种基于无人机多光谱高重叠度影像的森林林分特征反演方法。
背景技术
精准的森林林分特征提取,对于森林资源调查,生产力估测及生物多样性研究有重要意义,同时这些信息又可用于掌握森林空间结构及动态变化规律,对于森林可持续生产经营、生态环境建模及维持区域生态平衡提供数据支持。常规的林分特征提取主要依赖于地面调查方法,费时耗力,仅能够获取“点”上的数据,精度不高,且不易于在大面积上进行推广。
近年来基于多光谱高重叠影像技术进行林分特征反演研究为:Puliti等2015年在《Remote Sensing》第8卷上发表的“Inventory of Small Forest Areas Using anUnmanned Aerial System”,该研究利用装配了可见和近红外相机的固定翼无人机采集了挪威东南部北方森林的图像,并用地面参考数据拟合林分的Lorey′s高(按林分面积对林分高度的贡献加权)、优势木高度、树木株数、胸高断面积和蓄积量的线性模型;Vega等2008年在《International Journal of Remote Sensing》第11卷上发表的“Mapping canopyheight using a combination of digital stereo-photogrammetry and lidar”,该研究利用数字航空摄影测量影响生成数字表面模型,结合根据LiDAR数据获取的数字地面模型,从而建立冠层高度模型,用于森林冠层高度的测量。以上方法都是基于拼接后的二维影像提取特征,是基于单一的数据源来反演林分特征。
发明内容
发明目的:针对现有技术中的方法存在反演能力弱,反演精度不高的问题,本发明的目的在于提供一种基于无人机多光谱高重叠度影像的森林林分特征反演方法。
技术方案:为了解决上述问题,本发明所采用的技术方案如下:
一种基于无人机多光谱高重叠度影像的森林林分特征反演方法,借助无人机机载多光谱传感器获取高重叠度影像,利用波段组合和植被指数变换获取光谱特征,并利用空中三角测量算法与点云密集匹配算法生成数字点云数据,基于归一化点云数据和光谱数据来提取待测林冠层结构特征变量;最后结合地面实测待测林林分数据和提取的待测林冠层结构特征变量构建多元回归模型,反演待测林林分特征。
优选的,基于无人机多光谱高重叠度影像的森林林分特征反演方法的步骤如下:
(1)借助固定翼无人机搭载的多光谱相机采集高重叠度影像数据,并借助多旋翼无人机搭载的LiDAR传感器采集LiDAR原始点云数据;
(2)对点云数据的预处理生成数字高程模型;进一步处理获得归一化点云数据;
(3)对获取的高重叠度影像数据利用波段组合和植被指数变换提取光谱数据;
(4)基于归一化点云数据和光谱数据提取待测林林分冠层结构特征变量百分位高度、各层覆盖度和冠层体积与剖面特征变量;
(5)通过相关性分析筛选特征变量,首先筛选特征变量之间相关性低于0.6的特征变量,然后进一步筛选特征变量与各林分特征相关性高于0.6的特征变量;
(6)结合地面实测数据,以地面实测林分特征作为因变量,提取的待测林冠层结构特征变量作为自变量,分别构建多元回归模型。
优选的,所述步骤(2)中,点云数据的预处理为:首先去除LiDAR原始点云数据的噪音点,基于IDW滤波算法去除非地面点,然后通过计算每个像元内激光点高度的平均值,生成数字高程模型;
点云数据进一步处理由特征点提取和匹配、空中三角测量、三维点云加密组成,具体为:
(1)通过遥感平台拍摄研究区的数字影像,通过IMU实时记录每幅影像的姿态参数;
(2)对获取影像像对的特征点进行特征点匹配,通过光束法进行影像内外方位元素的解析并结合空中三角测量算法生成点云;
(3)在点云加密的基础上,通过加入地面控制点来对点云进行空间位置校正,并通过生成的数字高程模型对点云进行归一化处理,得到归一化点云数据。
优选的,所述数字高程模型的空间分辨率为0.5m。
优选的,所述步骤(4)林分冠层结构特征变量百分位高度、各层覆盖度和冠层体积与剖面特征变量分别为:
百分位高度变量包括4个冠层高度分布百分位数变量:H25、H50、H75、H95,冠层点云分布平均高度以上的覆盖度,冠层点云分布的变异系数;各层覆盖度变量包括点云数量在各百分数高度D3、D5、D7、D9以上的点占所有点云的百分比;冠层体积与剖面特征变量包括:Weibull函数对冠层高度分布剖面进行拟合得到2个剖面特征量Weibul α和Weibull β;冠层各结构类别的体积所占百分比,包括开放层、透光层、低光层和封闭层四个冠层结构类别。
优选的,所述步骤(6)地面实测数据:记录胸径不小于5cm的单木的树种,并测量胸径、树高和枝下高;胸径采用胸径尺测量,树高和枝下高分别用Vertex V超声和激光测高器测量;借助实测的单木信息,汇总得到样地级平均胸径、胸高断面积、Lorey′s平均树高、林分密度、蓄积量和地上生物量。
优选的,所述步骤(6)运用逐步回归法选择进入模型的自变量,即在预先给定的F水平下进行显著性检验:
t检验达不到显著水平,即p>0.1,则予以剔除;
t检验达到显著水平,即p<0.05,则予以进入。
优选的,所述步骤(6)为了降低自变量之间的相关性,通过主成分分析计算相关关系矩阵得到控制因子k,即最大特征根的平方根与最小特征根的比值,k小于30则模型进一步被选中。
优选的,采用决定系数R2、均方根误差RMSE和相对均方根误差rRMSE评价回归模型拟合的效果及估测精度:
式中,xi为某林分特征实测值;为某林分特征实测平均值;为某林分特征的模型估测值;n为样地的数量;i为某一个样地。
有益效果:与现有的技术相比,本发明具有的优点包括:
(1)本发明基于无人机平台的激光雷达技术更具有飞行成本低,数据获取便捷高效,可获取高密度点云等优势。
(2)本发明是在获取光谱信息的基础上通过立体摄影测量,从多光谱影像像对中提取三维点云,从而获取林分冠层的光谱信息和结构特征;由于林分特征与森林冠层光谱信息和结构特征有很好地相关性和机理联系,故本发明增强了林分特征反演的能力和精度,与使用其他相近遥感方法进行林分特征反演相比,其决定系数提升了10%以上。
(3)本发明同时运用光谱数据和点云数据,提取多组光谱特征和点云特征,并进行了特征变量优选,从而可以高质量反演林分特征。
(4)本发明不仅有利于特征变量的机理解释,也易于进行方法移植,即可以应用于不同树种的森林类型。
附图说明
图1从上到下分别为航空摄影测量点云条带效果图、激光雷达点云图及航空摄影测量和激光雷达点云叠加效果图;
图2为综合模型预测值与样地实测数据交叉验证散点图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合具体实施例对本发明的具体实施方式做详细的说明。
实施例1
研究区位于江苏省北部地区邳州市铁富镇,其地理位置北纬34°33′49″~34°34′23″,东经118°05′1″~~118°06′06″,属于半湿润温带季风气候;年降雨量约为903mm,最大降雨量集中在7、8月份梅雨季节;年平均温度约为13.9℃,无霜期211天;主要土壤类型为黑粘土壤,呈酸性;本地区银杏大规模种植始于1993年,总面积约5400hm2。
基于无人机多光谱高重叠度影像的森林林分特征反演方法,包括以下步骤:
(1)借助固定翼无人机采集获取的高重叠度航片的三维重建算法生成相对坐标准确的多光谱图像(蓝波段、绿波段、红波段、近红外、红边)和点云,并通过地面控制点(GCPs)进行几何校正,并借助多旋翼无人机搭载的激光雷达传感器进行激光雷达数据采集。
(2)点云数据预处理,首先去除LiDAR原始点云数据的噪音点,基于IDW滤波算法去除非地面点,然后通过计算每个像元内激光点高度的平均值,生成数字高程模型(DEM),数字高程模型(DEM)的空间分辨率为0.5m。
点云数据的进一步处理由特征点提取和匹配、空中三角测量、三维点云加密组成。首先,通过遥感平台拍摄研究区数字影像,通过IMU实时记录每幅影像的姿态参数。然后,获取影像像对的特征点并进行特征点匹配,通过光束法进行影像内外方位元素的解析并结合空中三角测量生成点云。在点云加密的基础上,通过加入地面控制点来对点云进行空间位置校正,并通过生成的数字高程模型对点云进行归一化处理,得到归一化点云数据。航空摄影测量点云条带效果图、激光雷达点云图及航空摄影测量和激光雷达点云叠加效果图如图1所示。由图1可知,通过以上方法获得的林区摄影测量点云结构完整,能够详细记录森林的三维结构信息,因此该摄影测量点云可以用来准确反演森林林分特征。
(3)将数字航空摄影测量影像与处理后的各个波段数据作为原始数据,进行波段组合、多种植被指数变换,提取多个光谱特征变量,包括5个原始单波段特征和18个植被指数特征,5个原始单波段特征是指航空摄影测量记录的红色、绿色、蓝色、红边和近红外这五个原始单波段。其中光谱特征变量已由像元尺度经均值计算转换为样地尺度,各多光谱植被指数如表1所示。
表1各多光谱植被指数
表1中的参考文献为:
[1]Y.J.Kaufman,D.Tanre.Atmospherically resistant vegetation index(ARVI)for EOS-MODIS[J].IEEE Trans.Geosci.Remote Sens.1992,30(2):261-270.
[2]Tucker C.J.Red and photographic infrared linear combinations formonitoring vegetation[J].Remote Sens.Environ.1979,8(2):127-150.
[3]A.R.Huete,H.Q.Liu,K.Batchily,W.van Leeuwen.A comparison ofvegetation indices over a global set of TM images for EOS-MODIS[J].RemoteSens.Environ.1997,59(3):440-451.
[4]Pinty B,Verstraete M.M.GEMI:A non-linear index to monitor globalvegetation from Satellites[J].Plant Ecol.1992,101(1):15-20.
[5]Gitelson A.A,Kaufman Y.J,Merzlyak M.N.Use of a green channel inremote sensing of global vegetation from EOS-MODIS[J].RemoteSens.Environ.1996,58(3):289-298.
[6]Robert E.Crippen.Calculating the vegetation index faster[J].RemoteSensing of Environment,1990.34(1):71-73.
[7]J.A.GAMON,J.S.SURFUS.Assessing leaf pigment content and activitywith a reflectometer[J].New Phytol.1999,143(1):105-117.
[8]Jiaguo Qi,A.Chehbouni,A.R.Huete,Y.H.Kerr,S.Sorooshian.A modifiedsoil adjusted vegetation index[J].Remote Sens.Environ.1994,48(2):119-126.
[9]JingMing Chen.Evaluation of vegetation indices and a modifiedsimple ratio for boreal application[J].Can.J.Remote Sens.1996,22(3):229-242.
[10]D.Haboudane,J.R.Miller,E.Pattey,P.J Zarco-Tejada,I.Strachan.Hyperspectral vegetation indices and novel algorithms forpredicting green LAI of crop canopies:modeling and validation in the contextof precision agriculture.Remote Sens.Environ.2004,90(3):337-352.
[11]J.W.Rouse,R.H.Haas,J.A.Schell,D.W.Deering,J.C Harlan.Monitoringthe Vemal Advancements and Retrogradation(Greenwave Efffect)of NatureVegetation(R).NASA/GSFC Final Report,NASA:Greenbelt,MD,USA,1974.
[12]Rondeaux G,Steven M,Baret F.Optimization of soil-adjustedvegetation indices[J].Remote Sens.Environ.1996,85(17),8369-8375.
[13]Roujean J.L,Breon F.M.Estimating PAR absorbed by vegetation frombidirectional reflectance measurements[J].Remote Sens.Environ.1995,51(3),375-384.
[14]Pearson R.L,Miller L.D.Remote Mapping of Standing Crop Biomassand Estimation of the Productivity of the Shortgrass Prairie[J].RemoteSensing of Environment.1972,45(2):7-12.
[15]HUETE A.R.A soil adjusted vegetation index(SAVI)[J].RemoteSensing of Environment.1988,25(3):295-309.
[16]C.F Jordan.Derivation of leaf area index from quality of light onthhe forest floor[J].Ecology 1969,50:663-666.
[17]Gong P,Pu R,Biging G.S,Larrieu M.R.Estimation of forest leaf areaindex using vegetation indices derived from Hyperion hyperspectral data[J].IEEE Trans.Geosci.Remote Sens.2003,41:1355-1362.
(4)基于归一化点云数据和光谱数据提取冠层结构特征变量,提取三组特征变量,即百分位高度、各层覆盖度和冠层体积与剖面特征变量;1)百分位高度变量包括:冠层高度分布百分位数(H25、H50、H75、H95),冠层点云分布平均高度以上的覆盖度(CCmean),冠层点云分布的变异系数(Hcv);2)各层覆盖度:点云数量在各百分数高度(30th、50th、70th、90th,即D3、D5、D7、D9)以上的点占所有点云的百分比;3)冠层体积与剖面特征变量:Weibull函数对冠层高度分布剖面进行拟合得到2个剖面特征量α,β(即Weibull α和Weibull β);冠层各结构类别体积占比,包括开放层,透光层,低光层和封闭层四个冠层结构类别,每个冠层结构类别的体积所占百分比(即OpenGap、Oligophotic、Euphotic、ClosedGap)。
(5)通过相关性分析筛选特征变量,即首先筛选特征变量之间相关性低于0.6的特征变量,然后进一步筛选特征变量与各林分特征相关性高于0.6的特征变量。
(6)根据历史森林资源调查数据和前期获取的卫星遥感影像数据在银杏人工林核心分布区内选取了5块1×1km的正方形大样地,然后在5个样地内按照典型选样的方法设置9块半径为15m的圆形小样地,样地的中心位置通过Trimble GeoXH6000 Handhelds手持GPS定位(结合JSCROS广域差分系统)进行定位,精度优于0.5m。样地调查过程中,对于胸径大于5cm的单木,逐一测定其胸径(胸径尺量测)、树高和枝下高(Vertex V超声波和激光测高器量测)。借助实测的单木信息,汇总得到样地级林分特征,即平均胸径、胸高断面积、Lorey′s平均树高(以每棵树的胸高断面积为权重,加权求和得到的平均树高)、林分密度、蓄积量和地上生物量。
将地面实测林分特征作为因变量,提取的冠层结构特征变量作为自变量,建立多元回归模型。运用逐步回归法选择进入模型的变量,即在预先给定的F水平下进行显著性检验,如果t检验达不到显著水平(p>0.1),则予以剔除;t检验达到显著水平(p<0.05)则予以进入。为了降低自变量之间的相关性,本方法通过主成分分析计算相关关系矩阵得到控制因子k(即最大特征根的平方根与最小特征根的比值),k小于30则模型进一步被选中。本方法采用决定系数(R2)、均方根误差(RMSE)和相对均方根误差(rRMSE)评价回归模型拟合的效果及估测精度:
式中xi为某林分特征实测值;为某林分特征实测平均值;为某林分特征模型估测值;n为样地的数量;i为某一个样地。
各林分结构参数估测模型及模型预测精度见表2,由表2可知各林分结构参数估测模型总体精度较高(R2=0.27-0.94,rRMSE=34.60-5.86%),DAP模型精度普遍高于多光谱模型。同时,集成点云和光谱特征的综合模型能够提升胸高断面积、蓄积量等林分特征的反演精度。因此,本方法对森林林分特征反演具有较高的精度,同时通过集成摄影测量影像及点云信息能够实现森林林分特征反演精度的进一步提升。
表2基于参数模型各林分特征反演交叉验证结果
综合模型预测值与样地实测数据交叉验证见图2。由图2可知,综合模型的模型预测结果与地面实测各林分特征较为接近。因此,本方法可以准确地对胸径、蓄积量等森林林分特征进行估算。
Claims (9)
1.一种基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于,借助无人机机载多光谱传感器获取高重叠度影像,利用波段组合和植被指数变换获取光谱特征,并利用空中三角测量算法与点云密集匹配算法生成数字点云数据,基于归一化点云数据和光谱数据来提取待测林冠层结构特征变量;最后结合地面实测待测林林分数据和提取的待测林冠层结构特征变量构建多元回归模型,反演待测林林分特征。
2.根据权利要求1所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:步骤如下:
(1)借助固定翼无人机搭载的多光谱相机采集高重叠度影像数据,并借助多旋翼无人机搭载的LiDAR传感器采集LiDAR原始点云数据;
(2)对点云数据的预处理生成数字高程模型,进一步处理获得归一化点云数据;
(3)对获取的高重叠度影像数据利用波段组合和植被指数变换提取光谱数据;
(4)基于归一化点云数据和光谱数据提取待测林林分冠层结构特征变量百分位高度、各层覆盖度和冠层体积与剖面特征变量;
(5)通过相关性分析筛选特征变量,首先筛选特征变量之间相关性低于0.6的特征变量,然后进一步筛选特征变量与各林分特征相关性高于0.6的特征变量;
(6)结合地面实测数据,以地面实测林分特征作为因变量,提取的待测林冠层结构特征变量作为自变量,分别构建多元回归模型。
3.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述步骤(2)中,点云数据的预处理为:首先去除LiDAR原始点云数据的噪音点,基于IDW滤波算法去除非地面点,然后通过计算每个像元内激光点高度的平均值,生成数字高程模型;
点云数据进一步处理由特征点提取和匹配、空中三角测量、三维点云加密组成,具体为:
(1)通过遥感平台拍摄研究区的数字影像,通过IMU实时记录每幅影像的姿态参数;
(2)对获取影像像对的特征点进行特征点匹配,通过光束法进行影像内外方位元素的解析并结合空中三角测量算法生成点云;
(3)在点云加密的基础上,通过加入地面控制点来对点云进行空间位置校正,并通过生成的数字高程模型对点云进行归一化处理,得到归一化点云数据。
4.根据权利要求3所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述数字高程模型的空间分辨率为0.5m。
5.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述步骤(4)林分冠层结构特征变量百分位高度、各层覆盖度和冠层体积与剖面特征变量分别为:
百分位高度变量包括4个冠层高度分布百分位数变量:H25、H50、H75、H95,冠层点云分布平均高度以上的覆盖度,冠层点云分布的变异系数;各层覆盖度变量包括点云数量在各百分数高度D3、D5、D7、D9以上的点占所有点云的百分比;冠层体积与剖面特征变量包括Weibull函数对冠层高度分布剖面进行拟合得到2个剖面特征量Weibulα和Weibullβ;冠层各结构类别的体积所占百分比,包括开放层、透光层、低光层和封闭层四个冠层结构类别。
6.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述步骤(6)地面实测数据:记录胸径不小于5cm的单木的树种,并测量胸径、树高和枝下高;胸径采用胸径尺测量,树高和枝下高分别用VertexV超声和激光测高器测量;借助实测的单木信息,汇总得到样地级平均胸径、胸高断面积、Lorey′s平均树高、林分密度、蓄积量和地上生物量。
7.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述步骤(6)运用逐步回归法选择进入模型的自变量,即在预先给定的F水平下进行显著性检验:
t检验达不到显著水平,即p>0.1,则予以剔除;
t检验达到显著水平,即p<0.05,则予以进入。
8.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:所述步骤(6)为降低自变量之间的相关性,通过主成分分析计算相关关系矩阵得到控制因子k,即最大特征根的平方根与最小特征根的比值,k小于30则模型进一步被选中。
9.根据权利要求2所述的基于无人机多光谱高重叠度影像的森林林分特征反演方法,其特征在于:采用决定系数R2、均方根误差RMSE和相对均方根误差rRMSE评价回归模型拟合的效果及估测精度:
式中,xi为某林分特征实测值;为某林分特征实测平均值;为某林分特征的模型估测值;n为样地的数量;i为某一个样地。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811060614.3A CN109212505B (zh) | 2018-09-11 | 2018-09-11 | 一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811060614.3A CN109212505B (zh) | 2018-09-11 | 2018-09-11 | 一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109212505A true CN109212505A (zh) | 2019-01-15 |
CN109212505B CN109212505B (zh) | 2019-11-05 |
Family
ID=64983868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811060614.3A Active CN109212505B (zh) | 2018-09-11 | 2018-09-11 | 一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109212505B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109754127A (zh) * | 2019-01-31 | 2019-05-14 | 浙江大学 | 基于无人机成像高光谱的水稻籽粒直链淀粉含量估测方法 |
CN110398465A (zh) * | 2019-07-06 | 2019-11-01 | 中国海洋大学 | 一种基于光谱遥感影像的养殖紫菜生物量测定方法 |
CN110530901A (zh) * | 2019-09-11 | 2019-12-03 | 中国科学院南京土壤研究所 | 融合宇宙射线中子法和无人机遥感的中小尺度土壤水监测系统及方法 |
CN110569805A (zh) * | 2019-09-10 | 2019-12-13 | 南京林业大学 | 一种基于无人机影像点云的人工林单木提取及立地质量评价的方法 |
CN111767865A (zh) * | 2020-06-30 | 2020-10-13 | 北部湾大学 | 一种利用航拍影像和激光数据反演红树林生物量的方法 |
CN112434617A (zh) * | 2020-11-26 | 2021-03-02 | 南京观微空间科技有限公司 | 一种基于多源遥感数据的森林生物量变化监测方法及系统 |
CN112698347A (zh) * | 2020-12-02 | 2021-04-23 | 北京华益瑞科技有限公司 | 地表植被参数监测装置、系统及方法 |
CN113030903A (zh) * | 2021-03-05 | 2021-06-25 | 深圳大学 | 一种基于无人机高光谱和激光雷达的苔草营养水平反演方法 |
CN113205475A (zh) * | 2020-01-16 | 2021-08-03 | 吉林大学 | 基于多源卫星遥感数据的森林高度反演方法 |
CN114965847A (zh) * | 2022-05-23 | 2022-08-30 | 西南大学 | 一种基于机载激光雷达的花椒氮累积量的无损监测方法 |
CN115629393A (zh) * | 2022-11-30 | 2023-01-20 | 北京市农林科学院智能装备技术研究中心 | 生物量估算方法、装置、系统、电子设备及存储介质 |
CN117520733A (zh) * | 2024-01-05 | 2024-02-06 | 云南师范大学 | 森林冠层高度与地理环境协变量关系确定方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867180A (zh) * | 2015-05-28 | 2015-08-26 | 南京林业大学 | 一种集成uav与lidar的森林林分特征反演方法 |
-
2018
- 2018-09-11 CN CN201811060614.3A patent/CN109212505B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867180A (zh) * | 2015-05-28 | 2015-08-26 | 南京林业大学 | 一种集成uav与lidar的森林林分特征反演方法 |
Non-Patent Citations (2)
Title |
---|
申鑫 等: ""基于高分辨率与高光谱遥感影像的北亚热带马尾松及次生落叶树种的分类"", 《植物生态学报》 * |
许子乾 等: ""集成高分辨率UAV影像与激光雷达点云的亚热带森林林分特征反演"", 《植物生态学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109754127A (zh) * | 2019-01-31 | 2019-05-14 | 浙江大学 | 基于无人机成像高光谱的水稻籽粒直链淀粉含量估测方法 |
CN110398465A (zh) * | 2019-07-06 | 2019-11-01 | 中国海洋大学 | 一种基于光谱遥感影像的养殖紫菜生物量测定方法 |
CN110569805B (zh) * | 2019-09-10 | 2023-05-05 | 南京林业大学 | 一种基于无人机影像点云的人工林单木提取及立地质量评价的方法 |
CN110569805A (zh) * | 2019-09-10 | 2019-12-13 | 南京林业大学 | 一种基于无人机影像点云的人工林单木提取及立地质量评价的方法 |
CN110530901A (zh) * | 2019-09-11 | 2019-12-03 | 中国科学院南京土壤研究所 | 融合宇宙射线中子法和无人机遥感的中小尺度土壤水监测系统及方法 |
CN110530901B (zh) * | 2019-09-11 | 2020-09-01 | 中国科学院南京土壤研究所 | 融合宇宙射线中子法和无人机遥感的中小尺度土壤水监测系统及方法 |
CN113205475A (zh) * | 2020-01-16 | 2021-08-03 | 吉林大学 | 基于多源卫星遥感数据的森林高度反演方法 |
CN113205475B (zh) * | 2020-01-16 | 2022-07-12 | 吉林大学 | 基于多源卫星遥感数据的森林高度反演方法 |
CN111767865A (zh) * | 2020-06-30 | 2020-10-13 | 北部湾大学 | 一种利用航拍影像和激光数据反演红树林生物量的方法 |
CN112434617A (zh) * | 2020-11-26 | 2021-03-02 | 南京观微空间科技有限公司 | 一种基于多源遥感数据的森林生物量变化监测方法及系统 |
CN112434617B (zh) * | 2020-11-26 | 2021-08-13 | 南京观微空间科技有限公司 | 一种基于多源遥感数据的森林生物量变化监测方法及系统 |
CN112698347A (zh) * | 2020-12-02 | 2021-04-23 | 北京华益瑞科技有限公司 | 地表植被参数监测装置、系统及方法 |
CN113030903A (zh) * | 2021-03-05 | 2021-06-25 | 深圳大学 | 一种基于无人机高光谱和激光雷达的苔草营养水平反演方法 |
CN113030903B (zh) * | 2021-03-05 | 2023-08-29 | 深圳大学 | 一种基于无人机高光谱和激光雷达的苔草营养水平反演方法 |
CN114965847A (zh) * | 2022-05-23 | 2022-08-30 | 西南大学 | 一种基于机载激光雷达的花椒氮累积量的无损监测方法 |
CN115629393A (zh) * | 2022-11-30 | 2023-01-20 | 北京市农林科学院智能装备技术研究中心 | 生物量估算方法、装置、系统、电子设备及存储介质 |
CN117520733A (zh) * | 2024-01-05 | 2024-02-06 | 云南师范大学 | 森林冠层高度与地理环境协变量关系确定方法及系统 |
CN117520733B (zh) * | 2024-01-05 | 2024-03-19 | 云南师范大学 | 森林冠层高度与地理环境协变量关系确定方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109212505B (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212505B (zh) | 一种基于无人机多光谱高重叠度影像的森林林分特征反演方法 | |
Zheng et al. | Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery | |
CN104867180B (zh) | 一种集成uav与lidar的森林林分特征反演方法 | |
Lu et al. | Mapping vegetation biophysical and biochemical properties using unmanned aerial vehicles-acquired imagery | |
CN109031344B (zh) | 一种全波形激光雷达和高光谱数据联合反演森林结构参数的方法 | |
CN110174359B (zh) | 一种基于高斯过程回归的航空高光谱影像土壤重金属浓度评估方法 | |
CN112287892B (zh) | 基于无人机高光谱和机器学习算法的乔木生物量测算方法 | |
CN110222903B (zh) | 一种基于无人机遥感的水稻产量预测方法及系统 | |
Xiao et al. | A temporally integrated inversion method for estimating leaf area index from MODIS data | |
CN111028096A (zh) | 一种天、空、地一体化数据融合的系统和方法 | |
CN110398466A (zh) | 基于遥感反演的农作物生长状态监测方法 | |
CN109061601A (zh) | 一种基于无人机激光雷达反演人工林森林结构参数的方法 | |
Motlagh et al. | Estimating and mapping forest biomass using regression models and Spot-6 images (case study: Hyrcanian forests of north of Iran) | |
CN108872964A (zh) | 基于无人机LiDAR数据的银杏人工林郁闭度提取方法 | |
CN108896021B (zh) | 基于航空摄影测量点云提取人工林林分结构参数的方法 | |
CN108981616B (zh) | 一种由无人机激光雷达反演人工林有效叶面积指数的方法 | |
CN111368817A (zh) | 一种基于地表类型进行热效应定量评价方法及系统 | |
He et al. | Using LiDAR data to measure the 3D green biomass of Beijing urban forest in China | |
CN108898070A (zh) | 一种基于无人机平台的高光谱遥感提取薇甘菊装置及方法 | |
CN109427091A (zh) | 一种基于高光谱遥感和摄影测量技术的生物量模型研建方法 | |
Fu et al. | Assessment of approaches for monitoring forest structure dynamics using bi-temporal digital aerial photogrammetry point clouds | |
Yang et al. | Forest canopy height mapping over China using GLAS and MODIS data | |
Broussard III et al. | Quantifying vegetation and landscape metrics with hyperspatial unmanned aircraft system imagery in a coastal oligohaline marsh | |
Palmer et al. | Modelling annual evapotranspiration in a semi-arid, African savanna: functional convergence theory, MODIS LAI and the Penman–Monteith equation | |
CN109212553A (zh) | 无人机LiDAR和随机森林提取银杏生物物理特性的方法 |
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 |