CN110647786A - 一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 - Google Patents
一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 Download PDFInfo
- Publication number
- CN110647786A CN110647786A CN201810677193.2A CN201810677193A CN110647786A CN 110647786 A CN110647786 A CN 110647786A CN 201810677193 A CN201810677193 A CN 201810677193A CN 110647786 A CN110647786 A CN 110647786A
- Authority
- CN
- China
- Prior art keywords
- grassland
- livestock
- sheep
- vegetation
- pasture
- 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
- 244000144972 livestock Species 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005516 engineering process Methods 0.000 title claims abstract description 16
- 244000025254 Cannabis sativa Species 0.000 title claims description 34
- 239000002028 Biomass Substances 0.000 claims abstract description 54
- 241001494479 Pecora Species 0.000 claims abstract description 49
- 241001465754 Metazoa Species 0.000 claims abstract description 46
- 244000038280 herbivores Species 0.000 claims abstract description 37
- 238000011156 evaluation Methods 0.000 claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 238000000605 extraction Methods 0.000 claims abstract description 4
- 238000009304 pastoral farming Methods 0.000 claims description 25
- 239000004459 forage Substances 0.000 claims description 24
- 241000283690 Bos taurus Species 0.000 claims description 14
- 238000011835 investigation Methods 0.000 claims description 12
- 241000196324 Embryophyta Species 0.000 claims description 7
- 231100000331 toxic Toxicity 0.000 claims description 6
- 230000002588 toxic effect Effects 0.000 claims description 6
- 239000003086 colorant Substances 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 241000283707 Capra Species 0.000 claims description 3
- 241000283074 Equus asinus Species 0.000 claims description 3
- 241001331845 Equus asinus x caballus Species 0.000 claims description 3
- 238000012790 confirmation Methods 0.000 claims description 3
- 235000021050 feed intake Nutrition 0.000 claims description 3
- 230000037406 food intake Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 230000000153 supplemental effect Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 9
- 230000007547 defect Effects 0.000 abstract 1
- 235000015278 beef Nutrition 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 5
- 229930002875 chlorophyll Natural products 0.000 description 4
- 235000019804 chlorophyll Nutrition 0.000 description 4
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 241000209504 Poaceae Species 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000002542 deteriorative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 244000144980 herd Species 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- 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/188—Vegetation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mining & Mineral Resources (AREA)
- General Business, Economics & Management (AREA)
- Marine Sciences & Fisheries (AREA)
- Human Resources & Organizations (AREA)
- Agronomy & Crop Science (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，可有效克服传统卫星遥感不能监测干枯植被和食草动物信息的缺点，面向枯草期草原开展草畜平衡监测，包括以下步骤：基于坡度‑角度‑距离限制的地形提取法，对无人机LIDAR点云数据进行分类，根据分类得到的地面点和植被点，计算非生长季的草地植被覆盖度和冠层高度，实现非生长季草地地上生物量反演；基于面向对象和监督分类法，获得可利用牧草比例、草地合理利用率等参数，计算理论载畜量；根据食草动物种类、数量、和每种动物的羊单位换算公式，计算实际载畜量；最后，根据理论载畜量和实际载畜量，计算超载率、剩余牧草可利用天数及需补饲牧草量，完成非生长季草畜平衡状况评估。
Description
技术领域
本发明涉及无人机LIDAR点云及航拍影像数据处理技术领域，具体涉及一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法
背景技术
我国是除澳大利亚外草地面积最大的国家，草原面积约占国土面积的41.2％，但长期以来受自然因素和人为因素双重作用，尤其是近几十年来超载过牧的影响，我国退化草原面积已达90％以上，严重退化草原占50％(侯向阳,2014)，且每年还以200万公顷的速度增加。为保护、建设和合理利用草原，1985年国家颁布了《草原法》，开始对草原实行“以草定畜、草畜平衡”制度，并对严重退化、沙化、盐碱化、石漠化和生态脆弱区的草原，实行禁牧、休牧制度。近年来，又开始实施“生态草牧业”(方精云等,2016)、“国家公园”(付梦娣等,2017)等生态文明建设重大举措，对草原可持续发展起到积极作用。但现有的草畜平衡监测技术发展缓慢，不足以支撑相关工作的顺利和高效开展，存在以下问题：
(1)现有技术难以直接监测非生长季草地生物量，主要通过饲草总量与每月的家畜饲草需求量估算获得。我国草地饲草集中产生于夏秋季，而家畜最缺草的季节往往是在冬春季(徐敏云,2014)。以往基于光学卫星数据的生长季草地生物量监测方法，过分依赖植被的叶绿素含量，而非生长季草地植被干枯，直接使用基于光学数据的监测方法估算生物量的精度低，不能满足非生长季的草畜平衡评估需要(Feng et al.,2017)。为了精确计算草地合理载畜量，学者们发展了考虑生物量年际变化的草畜平衡核定方法(李青丰,2011)，基于经验模型法将生物量最大月份的卫星遥感数据与地面调查结果建立关系，估算一年中生物量最大月份值，并根据生物量和家畜数量估算当年的草畜平衡(李梦娇等,2016)，或根据饲草剩余量与家畜饲草需求量，估算每月的草畜平衡(李青丰,2011)。该方法缩短了牧草剩余量核定周期(从5年缩短到1年或1月)，可对草畜平衡进行逐年或逐月预报，有助于牧民提前准备饲草，但受家畜数量月季变化、牧草损失等影响，每月的真实饲草存量可能与估算值差别很大。如，赵晨光and李青丰(2016)认为冬季可利用饲草量仅为生物量最大月份值的1/4。因此需探索使用LIDAR等可探测植被垂直结构信息的数据进行非生长季草地生物量反演。
(2)食草动物数据更新频率低。现有的家畜数量主要通过入户抽样调查方法统计获得，效率低，且由于牧户的积极性和抽样调查方法本身的不确定性，使得该方法获取的家畜数据存在较大不确定性。考虑到食草动物数量受自然灾害和人为影响变化较大，且统计数据缺乏空间位置信息，不适于精细尺度的草畜平衡评估。无人机高分影像等具有探测范围广、获取周期短、时效性强、成本低等优点，却没有用于食草动物识别与现实载畜量估算。
传统航空遥感平台主要为有人机，可搭载可见光、多光谱、LIDAR等多种传感器，进行生物量反演。与卫星遥感相比，航空遥感数据分辨率更高，更适合中小尺度的植被参数便捷提取，亦可为卫星遥感产品提供验证数据(Tang et al.,2014)。但有人机数据获取成本高、操作复杂，应用范围有限。而近年来兴起的无人机，可搭载普通彩色相机、多光谱相机、LIDAR传感器等，飞行高度可低至10～100m，获得厘米级、甚至毫米级精度的遥感数据，具有成本低和操作灵活等优势，已在草地植被覆盖度监测(Laliberte et al.,2010)、植株密度计算(Hardin and Jackson,2005)、草地类型自动识别等方面有所应用(Laliberte andRango,2011)。无人机携带LIDAR传感器能够不依赖植被叶绿素信息，低成本获取植被三维空间结构信息，进而结合地面实测数据开展草地地上生物量反演，且可同时携带普通光学相机获取厘米级分辨率的食草动物影像，并估算草场的现实载畜量。然而，至今尚未见有学者将无人机LIDAR和普通相机有机结合，用于非生长季的草、畜监测及草畜平衡评估研究，相关的监测评估理论和方法研究体系也亟待探索。
发明内容
针对现有草畜平衡评估技术中存在的不能对非生长季生物量和食草动物数据有效监测的问题，本发明提出一种适合非生长季的草畜平衡评估方法，其使用无人机LIDAR数据直接探测草地冠层高度和盖度，进而结合地面实测数据反演非生长季的草地生物量；然后，基于面向对象和监督分类法，对预设区的无人机真彩色影像进行解译，获取食草动物数据；最后，根据非生长季的草地生物量计算得到理论载畜量和根据食草动物数据计算得到现实载畜量，开展非生长季草畜平衡状况评估，以显著提高非生长季草畜平衡评估的精度和效率。
为了解决上述问题，本发明实例提出了一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，包括以下步骤：
步骤1、基于坡度-角度-距离限制的地形提取法，对预设区的无人机LIDAR点云数据进行监督分类，将点云分类为地面点、植被点、其他点，进而计算非生长季的草地植被覆盖度和冠层高度等植被参数，并结合地面实测数据，构建生物量反演模型，以反演非生长季草地地上生物量；
步骤2、基于面向对象和监督分类法，根据颜色、纹理、尺寸、隶属关系等，结合专家知识，对预设区的无人机航拍影像进行分割与解译，提取食草动物种类和数量等信息，并对解译精度精细评估；
步骤3、根据得到的地上生物量数据，结合无人机遥感、文献或实地调查方法获得的可利用牧草比例(不含毒杂草)、草地合理利用率等参数，以及草场确权等数据，生成可利用牧草数据，计算理论载畜量；
步骤4、根据食草动物种类、数量、采食范围、月数和每种动物的羊单位换算公式，推算食草动物羊单位，计算实际载畜量；
步骤5、根据理论载畜量和实际载畜量，计算超载率、剩余牧草可利用天数及需补饲牧草量，完成非生长季草畜平衡状况评估。
进一步地，所述步骤1中的草地冠层高度、覆盖度、生物量通过以下方法获取：
草地冠层高度按下式(1)计算：
CHM(i,j)＝DVSM(i,j)-DEM(i,j) (1)
式中，DVSM(i,j)代表由植被点拟合生成的数字植被表面模型中(i,j)点的海拔高度；DEM(i,j)代表由地面点拟合生成的数字地形模型中(i,j)点的海拔高度；
植被覆盖度为求解点窗口Wd内植被点个数占全部点的比例，按下式(2)计算：
FVC(i,j)＝(N_veg(u,v))/(N_(all)(u,v)),u,v∈W_d(i,j) (2)
其中，Nveg(u,v)为窗口Wd内植被点个数，Nall(u,v)为窗口Wd内包括植被、地面、其他点在内的所有点个数，Wd为边长d的正方形窗口，d通常可取值为0.5-1m。
草地地上生物量由从LIDAR点云提取的草地冠层高度和(或)植被覆盖度与实测生物量之间的拟合获得，按下式(3)计算：
Bio(i,j)＝f(CHM(i,j),FVC(i,j)，b) (3)
其中，b为待定参数；拟合函数f为线性函数(Bio＝b1x+b0)、对数型拟合函数(Bio＝b1lnx+b0)、指数函数或幂指数函数
中表现最好的拟合函数(实测值与拟合值间的相关系数R2最大，均方根误差RMSE最小)，x代表草地冠层高度CHM(i,j)和植被覆盖度FVC(i,j)两个自变量。
进一步地，所述步骤2中的食草动物遥感识别精度按下式(4)计算：
PA(C_i)＝(N_right(C_i))/(N_all(C_i)) (4)
其中，Nright(Ci)为准确识别为Ci类食草动物的对象数，Nall(Ci)为Ci类食草动物的地面调查对象数。
食草动物遥感识别误识率按下式(5)计算：
PA(C_i)＝(N_wrong(C_i))/(N_all(C_i)) (5)
式中，Nwrong(Ci)为遥感解译过程中错误识别为Ci类食草动物的对象数。
食草动物遥感识别漏识率按下式(6)计算：
PA(C_i)＝(N_undetected(C_i))/(N_all(C_i)) (6)
式中，Nundetected(Ci)为遥感解译过程中漏识的Ci类食草动物的对象数。
进一步地，所述步骤3中的理论载畜量，即草场合理载畜量的食草兽类动物羊单位按下式(7)计算：
N_R＝F_rest/(I×D_plan) (7)
其中，I为羊单位饲草日采食量，Dplan为计划放牧或刈割牧草饲喂天数(一年按365天计算，非生长季依据当地实际情况计算，本发明实例按9个月计算)，F为放牧或刈割草地可利用(剩余)标准干草量，按下式(8)计算：
Frest＝Y×E×U×H (8)
式中，Y为草地地上生物量，E为可食牧草比率(不含毒杂草，割草地不超过90％，放牧草地不超过100％)，U为放牧或刈割草地合理利用率(通常为50％)，H为草地标准干草折算系数(取决于草地类型和草地质量)。以上E、U、H通过查询文献和实地测量获得。
进一步地，所述步骤4中预设区的实际载畜量通过下式(9)确定：
其中，m为从航拍数据中提取的畜种类数，NPi为从航拍数据中提取的第i种畜类头只数，kpi为第i种畜类的羊单位折算系数(绵羊/山羊折算为1羊单位，黄牛/马/骡子为5羊单位，牦牛为4羊单位，乳牛为6.5羊单位，驴为3羊单位)，NP为m种畜类的羊单位总体数，即实际载畜量。
进一步地，所述步骤5中预设区的草地潜/超载率通过下式(10)确定：
W＝(N_R-N_P)/N_R (10)
式中，NR为草地合理载畜羊单位，NP为实际载畜羊单位。当W＝0，表明草地区域合理承载量与现存食草动物利用量达到平衡，表明草地利用适度。
当W>0，草地还有载畜潜力，则剩余牧草可利用天数通过下式(11)确定：
D_rest＝F_rest/(N_P×I) (11)
当W<0，表明草地已超载，需补饲牧草量通过下式(12)确定：
Fadding＝NP×I×Dplan-Frest (12)
式中，Dplan为计划放牧或刈割牧草饲喂天数(主要参考历年预设区的非生长季天数确定，本发明实例按9个月计算)。
本发明的上述技术方案的有益效果如下：
1、本发明实例提出的基于无人机LIDAR数据的生物量反演方法不依赖植被叶绿素信息，能够通过测量草地冠层高度和盖度来更准确地估算非生长季期间任一时间节点的草地生物量，较传统光学遥感有不可比拟的精度优势。实验表明，生物量反演的相对精度达85.9％，R2＝0.35，RMSE＝81.9g·m-2。
2、本发明实例提出的基于无人机航拍数据的食草动物数量提取方法，无需地面调查，即可实现对现实载畜量的准确估算，总体识别精度达89.71％，而传统基于地面调查的方法虽然准确率也很高，但受调查人员的主客观因素影响较大，成本更高，且有些调查区域地面难以进入时，准确度难以保证。
3、本发明实例采用无人机搭载LIDAR和真彩色相机，对预设区域的非生长季期内的草畜平衡进行有效遥感监测评估，结果可靠、调查效率高、可操作性强。该发明特别适合对精细尺度的非生长季草畜平衡评估，还可为卫星遥感产品提供验证数据，对防止广大偏远地区的草原超载过牧，减少非生长季家畜掉膘、甚至死亡等风险有重要意义；同时，在确保春季牧草留茬高度的基础上，可防止草场退化，预防和控制沙尘灾害等。
附图说明
图1为本发明实例的基于无人机LIDAR航测技术的非生长季草畜平衡评估流程示意图；
图2为本发明实例的基于无人机LIDAR数据的草地生物量反演流程示意图。
图3为本发明实例的基于无人机真彩色影像的食草动物数量估算流程示意图。
图4为本发明实例的肉牛无人机真彩色影像实例(主要有黑黄白三色)。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实例进行详细描述。
本发明提供一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，分别基于无人机LIDAR点云和真彩色影像实现非生长季的草地生物量估算和食草动物解译，进而对非生长季草地的草畜平衡进行评估，步骤如下：
步骤1、基于坡度-角度-距离限制的地形提取法，对预设区的无人机LIDAR点云数据进行监督分类，将点云分类为地面点、植被点、其他点，进而计算非生长季的草地植被覆盖度和冠层高度等植被参数，并结合地面实测数据，构建生物量反演模型，以反演非生长季草地地上生物量；
草地冠层高度按下式(1)计算：
CHM(i,j)＝DVSM(i,j)-DEM(i,j) (1)
式中，DVSM(i,j)代表由植被点拟合生成的数字植被表面模型中(i,j)点的海拔高度；DEM(i,j)代表由地面点拟合生成的数字地形模型中(i,j)点的海拔高度；
植被覆盖度为求解点窗口Wd内植被点个数占全部点的比例，按下式(2)计算：
FVC(i,j)＝(N_veg(u,v))/(N_(all)(u,v)),u,v∈W_d(i,j) (2)
其中，Nveg(u,v)为窗口Wd内植被点个数，Nall(u,v)为窗口Wd内包括植被、地面、其他点在内的所有点个数，Wd为边长d的正方形窗口，d通常可取值为0.5-1m。
草地地上生物量由从LIDAR点云提取的草地冠层高度和/或植被覆盖度与实测生物量之间的拟合获得，按下式(3)计算：
Bio(i,j)＝f(CHM(i,j),FVC(i,j)，b) (3)
其中，b为待定参数；拟合函数f为线性函数(Bio＝b1x+b0)、对数型拟合函数(Bio＝b1lnx+b0)、指数函数
或幂指数函数
中表现最好的拟合函数(实测值与拟合值间的相关系数R2最大，均方根误差RMSE最小)，x代表草地冠层高度CHM(i,j)和植被覆盖度FVC(i,j)两个自变量。
步骤2、基于面向对象和监督分类法，根据颜色、纹理、尺寸、隶属关系等，结合专家知识，对预设区的无人机航拍影像进行分割与解译，提取食草动物种类和数量等信息，并对解译精度精细评估；
食草动物遥感识别精度按下式(4)计算：
PA(C_i)＝(N_right(C_i))/(N_all(C_i)) (4)
其中，Nright(Ci)为准确识别为Ci类食草动物的对象数，Nall(Ci)为Ci类食草动物的地面调查对象数。
食草动物遥感识别误识率按下式(5)计算：
PA(C_i)＝(N_wrong(C_i))/(N_all(C_i)) (5)
式中，Nwrong(Ci)为遥感解译过程中错误识别为Ci类食草动物的对象数。
食草动物遥感识别漏识率按下式(6)计算：
PA(C_i)＝(N_undetected(C_i))/(N_all(C_i)) (6)
式中，Nundetected(Ci)为遥感解译过程中漏识的Ci类食草动物的对象数。
步骤3、根据得到的地上生物量数据，结合无人机遥感、文献或实地调查方法获得的可利用牧草比例(不含毒杂草)、草地合理利用率等参数，以及草场确权等数据，生成可利用牧草数据，计算理论载畜量；
草场合理载畜量的食草兽类动物羊单位按下式(7)计算：
N_R＝(F_rest)/(I×D_plan) (7)
其中，I为羊单位饲草日采食量，Dplan为计划放牧或刈割牧草饲喂天数(一年按365天计算，非生长季依据当地实际情况计算，本发明实例按9个月计算)，F为放牧或刈割草地可利用(剩余)标准干草量，按下式(8)计算：
Frest＝Y×E×U×H (8)
式中，Y为草地地上生物量，E为可食牧草比率(不含毒杂草，割草地不超过90％，放牧草地不超过100％)，U为放牧或刈割草地合理利用率(通常为50％)，H为草地标准干草折算系数(取决于草地类型和草地质量)。以上E、U、H通过查询文献和实地测量获得。
步骤4、根据食草动物种类、数量、采食范围、月数和每种动物的羊单位换算公式，推算食草动物羊单位，计算实际载畜量；
实际载畜量通过下式(9)确定：
其中，m为从航拍数据中提取的畜种类数，NPi为从航拍数据中提取的第i种畜类头只数，kpi为第i种畜类的羊单位折算系数(绵羊/山羊折算为1羊单位，黄牛/马/骡子为5羊单位，牦牛为4羊单位，乳牛为6.5羊单位，驴为3羊单位)，NP为m种畜类的羊单位总体数，即实际载畜量。
步骤5、根据理论载畜量和实际载畜量，计算超载率、剩余牧草可利用天数及需补饲牧草量，完成非生长季草畜平衡状况评估。
草地潜/超载率通过下式(10)确定：
W＝(N_R-N_P)/N_R (10)
式中，NR为草地合理载畜羊单位，NP为实际载畜羊单位。当W＝0，表明草地区域合理承载量与现存食草动物利用量达到平衡，表明草地利用适度。
当W>0，草地还有载畜潜力，则剩余牧草Frest可利用天数通过下式(11)确定：
D_rest＝F_rest/(N_P×I) (11)
当W<0，表明草地已超载，需补饲牧草量通过下式(12)确定：
Fadding＝NP×I×Dplan-Frest (12)
式中，Dplan为计划放牧或刈割牧草饲喂天数(主要参考历年预设区的非生长季天数确定，本发明实例按9个月计算)。
图1展示了本发明实例的基于无人机LIDAR航测技术的非生长季草畜平衡评估流程。
下面利用一套无人机LIDAR点云和真彩色影像数据，对本发明——一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法——的效果和效率进行了测试，并和传统基于地面调查的方法进行了对比。本发明实例中，评估在内蒙古呼伦贝尔草原站的肉牛放牧样地进行。该放牧样地有6个放牧梯度，每个放牧梯度有3个重复。W/M/E0-5代表每个小区分别放养0、2、3、4、6、8头200-300kg的标准牛，理论上每个小区之间由于有栅栏，肉牛不能穿越(实际上个别会穿越栅栏)，因而每个小区内的草地冠层高度、盖度、生物量都明显不同，每个小区面积约90亩。无人机LIDAR数据由八旋翼无人机LI-AIR搭载Velodyne公司的HDL-32E传感器，于2015年8月底获取。无人机相机数据由组装固定翼无人机搭载CanonPowerShot G16。具体参数见表2。
无人机LIDAR点云和真彩色影像数据的相关参数如表1和表2所示。
表1无人机Lidar传感器及数据获取参数
表2无人机相机及数据获取参数
上述基于无人机LIDAR点云对非生长季草地生物量反演时，在完成对草地冠层高度和植被覆盖度的基础上，需找到最适合生物量反演的模型。本发明实例在使用线性、非线性模型，分析单独基于草地冠层高度和植被覆盖度，以及组合使用草地冠层高度和植被覆盖度构建生物量模型的基础上，发现组合草地冠层高度和植被覆盖度可更有效的反演草地地上生物量。实验中采集地面样方90个，其中随机抽取1/4作为验证数据，其余用作训练模型。生物量反演模型为：Bio＝35.09CHM-0.047FVC+8.71，R2＝0.341，RMSE＝82.1g·m-2。若使用仅单一指标，建议使用草地冠层高度，Bio＝34.79CHM+7.08，R2＝0.34，RMSE＝82.1g·m-2。图2展示了本发明实例的基于无人机LIDAR数据的草地生物量反演流程。表3展示了每个放牧样区的生物量反演结果及理论载畜量核定结果。结果表明：8月底每个小区的草地地上生物量为5216.22～16383.51kg，理论载畜量为4.76-14.96羊单位,总计理论载畜量153.04羊单位。
随后，基于面向对象和监督分类法，根据颜色、纹理、尺寸、隶属关系等，结合专家知识，对预设区的无人机航拍影像进行分割与解译，提取了食草动物种类和数量等信息。图4为本发明实例的肉牛无人机真彩色影像实例(主要有黑黄白三色)。表4展示了食草动物遥感解译及现实载畜量羊单位折算结果。由表4可知，本放牧样地只有肉牛一种食草动物家畜，除部分误识外，大多数小区的遥感识别较为准确，总体识别精度达89.71％，误识率为2.94％，漏识率为5.88％。考虑到地面调查与无人机影像获取时间存在一定时间差，受牧草存量和质量影响，部分肉牛会跨越栅栏到其他小区，地面调查也存在一定误差，因此遥感识别精度会更高。
最后，根据理论载畜量和实际载畜量，计算超载率、剩余牧草可利用天数及需补饲牧草量，完成非生长季草畜平衡状况评估。表5展示了呼伦贝尔放牧样地非生长季草畜平衡评估结果。结果表明，18个样区中，仅5个样区未超载，其余13个小区均已超载，超载率达23.45(M1)～681.59％(W5)。按整个放牧样地计算，平均超载133.49％，剩余牧草可利用数为5.14个月，需补饲56.65吨标准干草。
需要注意的是，本发明实例的评估方法在内蒙古呼伦贝尔草原站的肉牛放牧样地进行理论验证，不代表内蒙古呼伦贝尔草原站肉牛放牧平台在非生长季会继续放牧或刈割牧草。
表3生物量反演结果及理论载畜量核定
注：每个羊单位每日按消耗1.8kg干草，按放牧或刈割牧草饲喂9个月计算。
表4食草动物遥感解译结果及现实载畜量羊单位折算
注：本样地只有肉牛一种食草家畜。每头200-300kg的标准肉牛按4羊单位折算。
表5呼伦贝尔放牧样地非生长季草畜平衡评估结果
注：当W>0表示潜载，反之超载。
对比结果表明，本发明实例所提出的一种基于无人机LIDAR数据的生物量反演方法，不依赖植被叶绿素信息，能够通过测量草地冠层高度和盖度来更准确地估算非生长季期间任一时间节点的草地生物量，生物量反演结果相对于最大生物量的精度达85.9％，R2＝0.35，RMSE＝81.9g·m-2，较传统光学遥感有不可比拟的精度优势，可提供高精度的非生长季草地理论载畜量数据。本发明实例提出的基于无人机航拍数据的食草动物数量提取方法，无需地面调查，即可实现对现实载畜量的准确估算，总体识别精度达89.71％，可作为地面食草动物调查的一种替代手段，具有成本低、结果可靠、野外工作量小等优势。
此外，根据上述方式获得理论载畜量和实际载畜量，进而计算超载率、剩余牧草可利用天数及需补饲牧草量，可实现非生长季草畜平衡状况快速精准评估，对非生长季草地的牧草存量合理利用和及时补饲，防止草地退化、非生长季家畜掉膘、甚至死亡等有重要理论和现实意义，在精准草牧业、国家生态文明建设等领域具有可预期的应用前景。发明方法可供牧户、草原管理者、防沙治沙办等相关人员参考。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，包括：
步骤1、基于坡度-角度-距离限制的地形提取法，对预设区的无人机LIDAR点云数据进行监督分类，将点云分类为地面点、植被点、其他点，进而计算非生长季的草地植被覆盖度和冠层高度等植被参数，并结合地面实测数据，构建生物量反演模型，以反演非生长季草地地上生物量；
步骤2基于面向对象和监督分类法，根据颜色、纹理、尺寸、隶属关系等，结合专家知识，对预设区的无人机航拍影像进行分割与解译，提取食草动物种类和数量等信息，并对解译精度精细评估；
步骤3、根据得到的地上生物量数据，结合无人机遥感、文献或实地调查方法获得的可利用牧草比例(不含毒杂草)、草地合理利用率等参数，以及草场确权等数据，生成可利用牧草数据，计算理论载畜量；
步骤4、根据食草动物种类、数量、采食范围、月数和每种动物的羊单位换算公式，推算食草动物羊单位，计算实际载畜量；
步骤5、根据理论载畜量和实际载畜量，计算超载率、剩余牧草可利用天数及需补饲牧草量，完成非生长季草畜平衡状况评估。
2.根据权利要求1所述的一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，所述步骤1中的草地冠层高度、覆盖度、生物量通过以下方法获取：
根据公式(1)计算位于(i,j)点的草地冠层高度；
CHM(i,j)＝DVSM(i,j)-DEM(i,j) (1)
其中，DVSM(i,j)代表由植被点拟合生成的数字植被表面模型中(i,j)点的海拔高度；DEM(i,j)代表由地面点拟合生成的数字地形模型中(i,j)点的海拔高度；
根据公式(2)计算植被覆盖度，所述植被覆盖度为求解点窗口Wd内植被点个数占全部点的比例：
其中，Nveg(u,v)为窗口Wd内植被点个数，Nall(u,v)为窗口Wd内包括植被、地面、其他点在内的所有点个数，Wd为边长d的正方形窗口，d通常可取值为0.5-1m；
根据公式(3)计算草地地上生物量，通过从LIDAR点云提取的草地冠层高度和(或)植被覆盖度与实测生物量之间的拟合获得：
Bio(i,j)＝f(CHM(i,j),FVC(i,j)，b) (3)
其中，b为待定参数；拟合函数f为线性函数(Bio＝b1x+b0)、对数型拟合函数(Bio＝b1lnx+b0)、指数函数
或幂指数函数
中表现最好的拟合函数(实测值与拟合值间的相关系数R2最大，均方根误差RMSE最小)，x代表草地冠层高度CHM(i,j)、植被覆盖度FVC(i,j)等自变量。
3.根据权利要求1所述的一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，所述步骤2中的食草动物遥感识别精度按下式(4)计算：
其中，Nright(Ci)为准确识别为Ci类食草动物的对象数，Nall(Ci)为Ci类食草动物的地面调查对象数；
食草动物遥感识别误识率按下式(5)计算：
式中，Nwrong(Ci)为遥感解译过程中错误识别为Ci类食草动物的对象数。
食草动物遥感识别漏识率按下式(6)计算：
式中，Nundetected(Ci)为遥感解译过程中漏识的Ci类食草动物的对象数。
4.根据权利要求1所述的一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，所述步骤3中的理论载畜量羊单位按下式(7)计算：
其中，I为羊单位饲草日采食量，Dplan为计划放牧或刈割牧草饲喂天数(一年按365天计算，非生长季依据当地实际情况计算，本发明实例按9个月计算)，F为放牧或刈割草地可利用(剩余)标准干草量，按下式(8)计算：
Frest＝Y×E×U×H (8)
式中，Y为草地地上生物量，E为可食牧草比率(不含毒杂草，割草地不超过90％，放牧草地不超过100％)，U为放牧或刈割草地合理利用率(通常为50％)，H为草地标准干草折算系数(取决于草地类型和草地质量)。以上E、U、H通过查询文献和实地测量获得。
5.根据权利要求1所述的一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，所述步骤4中预设区的实际载畜量通过下式(9)确定：
其中，m为从航拍数据中提取的畜种类数，NPi为从航拍数据中提取的第i种畜类头/只数，kpi为第i种畜类的羊单位折算系数(绵羊/山羊折算为1羊单位，黄牛/马/骡子为5羊单位，牦牛为4羊单位，乳牛为6.5羊单位，驴为3羊单位)，NP为m种畜类的羊单位总体数，即实际载畜量。
6.根据权利要求1所述的一种基于无人机LIDAR航测技术的非生长季草畜平衡评估方法，其特征在于，所述步骤5中预设区的草地潜/超载率通过下式(10)确定：
式中，NR为草地合理载畜羊单位，NP为实际载畜羊单位。当W＝0，表明草地区域合理承载量与现存食草动物利用量达到平衡，表明草地利用适度；
当W>0，草地还有载畜潜力，则剩余牧草Frest可利用天数通过下式(11)确定：
当W<0，表明草地已超载，需补饲牧草量通过下式(12)确定：
Fadding＝NP×I×Dplan-Frest (12)
式中，Dplan为计划放牧或刈割牧草饲喂天数(主要参考历年预设区的非生长季天数确定，本发明实例按9个月计算)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810677193.2A CN110647786B (zh) | 2018-06-27 | 2018-06-27 | 一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810677193.2A CN110647786B (zh) | 2018-06-27 | 2018-06-27 | 一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110647786A true CN110647786A (zh) | 2020-01-03 |
| CN110647786B CN110647786B (zh) | 2022-04-26 |
Family
ID=68988759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810677193.2A Active CN110647786B (zh) | 2018-06-27 | 2018-06-27 | 一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110647786B (zh) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111582246A (zh) * | 2020-06-10 | 2020-08-25 | 河西学院 | 一种基于高寒草甸草地产草量估算放牧率的方法和系统 |
| CN111666855A (zh) * | 2020-05-29 | 2020-09-15 | 中国科学院地理科学与资源研究所 | 基于无人机的动物三维参数提取方法、系统及电子设备 |
| CN111695606A (zh) * | 2020-05-25 | 2020-09-22 | 中国科学院东北地理与农业生态研究所 | 一种多类型城市绿地分类方法 |
| CN112216109A (zh) * | 2020-10-12 | 2021-01-12 | 无锡臻永科技有限公司 | 采用数据库存储的信息分析平台及方法 |
| CN112734224A (zh) * | 2021-01-06 | 2021-04-30 | 中国科学院地理科学与资源研究所 | 一种草食性野生动物造成牧场损失的评估方法 |
| CN113222453A (zh) * | 2021-05-28 | 2021-08-06 | 中国农业科学院草原研究所 | 一种基于北斗导航的天然草地草畜平衡评估方法、系统、设备及存储介质 |
| CN113885060A (zh) * | 2021-09-26 | 2022-01-04 | 中国农业科学院草原研究所 | 基于无人机遥感技术的放牧强度监测方法 |
| CN115049248A (zh) * | 2022-06-10 | 2022-09-13 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种草地绿色承载力测算方法 |
| CN116310798A (zh) * | 2023-02-13 | 2023-06-23 | 中国林业科学研究院资源信息研究所 | 天然草地合理载畜量高精度遥感估算方法 |
| CN116414931A (zh) * | 2023-02-21 | 2023-07-11 | 中国科学院地理科学与资源研究所 | 一种国家公园家畜分布密度空间化的方法 |
| EP4245133A1 (en) * | 2022-03-16 | 2023-09-20 | Digitanimal, S.L. | Characterization of pasture for improved and sustainable grazing and feeding management of livestock |
| CN116993058A (zh) * | 2023-09-26 | 2023-11-03 | 中国农业科学院北京畜牧兽医研究所 | 放牧草地精准补饲方法和装置 |
| CN117456364A (zh) * | 2023-11-09 | 2024-01-26 | 中国农业科学院农业信息研究所 | 一种基于SfM和草地高度因子的草地地上生物量估测方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105842707A (zh) * | 2015-01-15 | 2016-08-10 | 兰州大学 | 基于无人机遥感图像的草地地上生物量测算方法及装置 |
| CN106199627A (zh) * | 2016-09-14 | 2016-12-07 | 中国农业科学院农业资源与农业区划研究所 | 草地植被参数获取方法 |
| US20180049407A1 (en) * | 2016-08-22 | 2018-02-22 | International Business Machines Corporation | Unmanned aerial vehicle for determining geolocation foraging zones |
-
2018
- 2018-06-27 CN CN201810677193.2A patent/CN110647786B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105842707A (zh) * | 2015-01-15 | 2016-08-10 | 兰州大学 | 基于无人机遥感图像的草地地上生物量测算方法及装置 |
| US20180049407A1 (en) * | 2016-08-22 | 2018-02-22 | International Business Machines Corporation | Unmanned aerial vehicle for determining geolocation foraging zones |
| CN106199627A (zh) * | 2016-09-14 | 2016-12-07 | 中国农业科学院农业资源与农业区划研究所 | 草地植被参数获取方法 |
Non-Patent Citations (2)
| Title |
|---|
| JINGLI HE: "Study on threshold expansion of the carrying capacity in desert grassland within the environmental capacity of grass and water resources", 《2011 INTERNATIONAL CONFERENCE ON REMOTE SENSING, ENVIRONMENT AND TRANSPORTATION ENGINEERING》 * |
| 钱拴琴: "青藏高原载畜能力及草畜平衡状况研究", 《生态经济》 * |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111695606B (zh) * | 2020-05-25 | 2022-10-04 | 中国科学院东北地理与农业生态研究所 | 一种多类型城市绿地分类方法 |
| CN111695606A (zh) * | 2020-05-25 | 2020-09-22 | 中国科学院东北地理与农业生态研究所 | 一种多类型城市绿地分类方法 |
| CN111666855A (zh) * | 2020-05-29 | 2020-09-15 | 中国科学院地理科学与资源研究所 | 基于无人机的动物三维参数提取方法、系统及电子设备 |
| CN111666855B (zh) * | 2020-05-29 | 2023-06-30 | 中国科学院地理科学与资源研究所 | 基于无人机的动物三维参数提取方法、系统及电子设备 |
| CN111582246A (zh) * | 2020-06-10 | 2020-08-25 | 河西学院 | 一种基于高寒草甸草地产草量估算放牧率的方法和系统 |
| CN112216109A (zh) * | 2020-10-12 | 2021-01-12 | 无锡臻永科技有限公司 | 采用数据库存储的信息分析平台及方法 |
| CN112734224A (zh) * | 2021-01-06 | 2021-04-30 | 中国科学院地理科学与资源研究所 | 一种草食性野生动物造成牧场损失的评估方法 |
| CN113222453A (zh) * | 2021-05-28 | 2021-08-06 | 中国农业科学院草原研究所 | 一种基于北斗导航的天然草地草畜平衡评估方法、系统、设备及存储介质 |
| CN113222453B (zh) * | 2021-05-28 | 2024-01-16 | 中国农业科学院草原研究所 | 一种基于北斗导航的天然草地草畜平衡评估方法、系统、设备及存储介质 |
| CN113885060A (zh) * | 2021-09-26 | 2022-01-04 | 中国农业科学院草原研究所 | 基于无人机遥感技术的放牧强度监测方法 |
| EP4245133A1 (en) * | 2022-03-16 | 2023-09-20 | Digitanimal, S.L. | Characterization of pasture for improved and sustainable grazing and feeding management of livestock |
| WO2023175095A1 (en) * | 2022-03-16 | 2023-09-21 | Digitanimal, S.L. | Characterization of pasture for improved and sustainable grazing and feeding management of livestock |
| CN115049248A (zh) * | 2022-06-10 | 2022-09-13 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种草地绿色承载力测算方法 |
| CN116310798A (zh) * | 2023-02-13 | 2023-06-23 | 中国林业科学研究院资源信息研究所 | 天然草地合理载畜量高精度遥感估算方法 |
| CN116414931A (zh) * | 2023-02-21 | 2023-07-11 | 中国科学院地理科学与资源研究所 | 一种国家公园家畜分布密度空间化的方法 |
| CN116414931B (zh) * | 2023-02-21 | 2024-05-14 | 中国科学院地理科学与资源研究所 | 一种国家公园家畜分布密度空间化的方法 |
| CN116993058A (zh) * | 2023-09-26 | 2023-11-03 | 中国农业科学院北京畜牧兽医研究所 | 放牧草地精准补饲方法和装置 |
| CN116993058B (zh) * | 2023-09-26 | 2024-05-03 | 中国农业科学院北京畜牧兽医研究所 | 放牧草地精准补饲方法和装置 |
| CN117456364A (zh) * | 2023-11-09 | 2024-01-26 | 中国农业科学院农业信息研究所 | 一种基于SfM和草地高度因子的草地地上生物量估测方法及系统 |
| CN117456364B (zh) * | 2023-11-09 | 2024-04-26 | 中国农业科学院农业信息研究所 | 一种基于SfM和草地高度因子的草地地上生物量估测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110647786B (zh) | 2022-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110647786B (zh) | 一种基于无人机lidar航测技术的非生长季草畜平衡评估方法 | |
| Higgins et al. | Vegetation sampling and measurement | |
| Reagan | Habitat selection in the three-toed box turtle, Terrapene carolina triunguis | |
| Xu et al. | Classification method of cultivated land based on UAV visible light remote sensing | |
| Radtke et al. | An evaluation of overhead laser scanning to estimate herbage removals in pasture quadrats | |
| Song et al. | Estimating effective leaf area index of winter wheat using simulated observation on unmanned aerial vehicle-based point cloud data | |
| Jurišić et al. | Current status and perspective of remote sensing application in crop management | |
| Mason et al. | Measurement of habitat predictor variables for organism-habitat models using remote sensing and image segmentation | |
| Chen et al. | 3D model construction and ecological environment investigation on a regional scale using UAV remote sensing | |
| Hoffmeister et al. | 3D terres-trial laser scanning for field crop modelling | |
| Liu et al. | Prediction of trees diameter at breast height based on unmanned aerial vehicle image analysis | |
| McCarthy et al. | Field evaluation of automated site-specific irrigation for cotton and perennial ryegrass using soil-water sensors and Model Predictive Control | |
| Bretas et al. | Precision livestock farming applied to grazingland monitoring and management—A review | |
| Li et al. | Quantifying grazing intensity in China using high temporal resolution MODIS data | |
| CN114358615A (zh) | 一种动态草畜平衡监测方法 | |
| Klaver et al. | Associating seasonal range characteristics with survival of female white‐tailed deer | |
| Fertu et al. | Monitoring the vegetation of agricultural crops using drones and remote sensing-comparative presentation. | |
| Lussem et al. | Feasibility study of directly georeferenced images from low-cost unmanned aerial vehicles for monitoring sward height in a long-term experiment on grassland | |
| Yang et al. | Flowering leafy spurge (Euphorbia esula) detection using unmanned aerial vehicle imagery in biological control sites: Impacts of flight height, flight time and detection method | |
| Hema et al. | First assessment of age and sex structures of elephants by using dung size analysis in a West African savannah | |
| Yang et al. | Field monitoring of fractional vegetation cover based on UAV low-altitude remote sensing and machine learning | |
| Meghraoui et al. | Towards 3D CNN for precise crop yield estimation using multimodal remote sensing data: case study of wheat in Morocco | |
| Swe et al. | Using digital photography to monitor changes in biocrusts and ground cover in a savanna rangeland | |
| Edwards | Limited Movement of Individuals in a Population of the Migratory Grasshopper, Melanoplus bilituratus (Walker)(Acrididae) at Kamloops, British Columbia1 | |
| Zainuddin et al. | Utilizing quadcopter as lars image platform to determine the paddy spectral and growth parameter |
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 |