交流园地

  土地资源是重要的生产资料，是人们赖以生存和发展的基础，土地利用现状和变化，不仅影响社会经济的持续发展，而且还间接影响环境的变化[1-2]。为此，人们迫切期望能快速、客观、准确地监测土地利用现状，并及时掌握变化信息，以实现土地资源的科学管理[3]。由于地形、气候和交通的影响和限制，传统生态调查是非常艰辛的工作，不仅采集数据困难，获取的结构化信息十分有限，而且测量结果受人为主观因素的影响较大。
  对于林地面积调查，过去采用的是人工测量与遥感卫星影像的方式。人工测量的方式需要耗费大量的人力和时间成本，存在效率低、精度低和更新周期冗长等问题；采用遥感卫星影像计算林地面积的方式虽然节省了人力和时间成本，但是普通卫星影像的分辨率较低，实时性相对落后，一般单位获取高分卫星影像较为困难。
  三维激光扫描和无人机技术的出现，给林业测量工作提供了新的技术方法。机载激光雷达技术可以在短时间内获取大范围的地形与地表物数据，大大降低了野外测量工作的工作量，极大的提高了数据获取的效率和精度，提高了林业资源调查的效率和准确度[4-5]。特别是在人员难以到达的地区，可以大大提高调查效率。利用获取的高精度点云数据，可以用来分析林木覆盖面积变化情况，分析林木类型、分布情况和生长阶段概况等[6]。本文利用机载激光雷达（Light Detection and Ranging，LiDAR）系统获取某区域的点云数据来调查该区域的林地面积与分布情况。
  一、机载LiDAR系统
  机载LiDAR系统以无人机为平台，传感器包括激光扫描仪（Laser Scanner,LS）、全球导航卫星系统（Global Navigation Satellite System,GNSS）和惯性测量单元（Inertial Measurement Unit,IMU）等。激光扫描仪安装在飞机的下端，用于向地面发射激光脉冲，并采集从发出到返回所用的时间、反射强度和扫描角度等信息[7-8]。利用GNSS测得的平台坐标和IMU获得的飞机姿态信息，计算得到激光所达到的地面与激光发出点的距离。GNSS和IMU用于确定飞机飞行时的位置和姿态，从而定义了距离测量值的原点，经过坐标变换，可以精确地获得激光所到达物体表面的三维坐标信息。
  机载LiDAR系统搭载的激光扫描仪具有多次回波特性，激光脉冲在穿越植被空隙时，可返回树冠、树枝、地面等多个回波数据，具有更强穿透性[9]；同时，机载LiDAR系统获取的点云数据可以达到10cm左右的精度，数据精度较高；无人机载体相对于车辆载体，其作业更加灵活，通达性更好，能够在车载不易采集和无法采集区域获取点云数据，在林业精准调查、电力线巡检和海岸带检测方面具有更强的适用性。
  二、作业方案
  在正式开展项目的全流程中，按照相关规范进行工程过程控制与管理。
  1.成立项目组。为了保证工程项目生产的顺利进行，成立项目组，由组长全面协调内外业人员的工作。
  2.工作质量控制。项目组成员严格按照质量相关要求开展数据处理和资料整理工作，严格执行相关数据处理规范；做好资料的互检与审核工作；成果资料必须经过质量审核方可交付。为达到质量管理要求需要做到：(1)切实重视成果质量，对不合格的测量数据坚决予以否定，推倒重来;(2)严格按照相关数据处理规范进行作业，并严格遵守操作规程，确保按时完成;(3)进行工序管理控制,严格抓好各道工序间的过程控制，以保证每一阶段数据的真实可靠。
  3.安全管理。严格按照安全管理体系要求进行信息安全管理，所有测量资料数据及存储均按照要求严格控制。
  4.生产计划及接收准则。按照计划对数据进行生产，生产过程中严把质量关，从数据预处理到后处理层层把关，数据提交质检小组进行检查。
  整体工作方案分为外业数据采集、内业数据解算和数据处理与分析三部分，流程图如图1所示。

  （一）外业数据采集
  外业数据采集严格按照相关的技术标准规范，测量前充分准备工作，测量中实时关注测量数据的完整性，按照计划进行数据获取，确保获取的数据完整有效。
  1.检核点测量
  使用GNSS接收机连接千寻CORS，在测区内采集一定数量的检核点，用于后期对成果进行精度检查。
  2.基准站架设
  为保证POS系统的定位精度，在航摄飞行前30分钟，需要在测区附近（距离小于30km）建立GNSS基准站，基准站上架设高精度GNSS信号接收机，该接收机与无人机POS同步记录。基准站架设时严格对中整平，量取天线高，基准站开机采集采样间隔为1HZ的静态数据。后处理中利用PPK后处理差分技术解算出无人机的精确位置。
  3.仪器安装
  仪器在不用时应存放于常温，干燥的环境中，保持机身清洁；激光头在长时间使用后，镜头上往往会沾有灰尘，需要用专用镜头纸擦拭干净。
  仪器安装过程中确保每一根连接线都连接紧，每一个连接件都连接紧。检查通电后的指示灯的情况，确保仪器正常工作。当无人机长期不用，或者无人机飞行明显的方向不正确则需要进行磁罗盘校准。检查通电后三个指示灯的情况（表1），三个灯常亮表示仪器正常工作。
  
  4.飞行航线规划
  航线设计可以在测量实施前或者现场进行规划。利用地面站软件生成相应测区的测线。
  5.机载激光扫描
  由于IMU容易受到累计误差的影响，因此在正式开始测量之前，需要对IMU进行初始化，以提高IMU的后处理精度。具体实施方法为：无人机绕“8”字飞行，观察仪器的PDOP值，当PDOP值稳定在1.5以下时，说明初始化已经完成，可以开始测量。
  无人机按照行高50米，飞行速度5m/s，航线间隔60米的规划进行作业。在飞行过程中，无人机飞行范围超出规定范围线50m，以保证飞行数据的精度以及数据的完整性。
  6.过程实时监控
  测量中可实时观察INS数据和激光数据的存储量，INS存储数据量稳定在0.01MB/s，激光存储数据量控制在2.50MB/s左右。
  （二）内业数据解算
  将无人机POS系统的数据与基站获得的GNSS数据利用Inertial Explorer软件进行联合解算，得到作业过程中无人机的准确位置信息与姿态信息。将导出的POS数据导入ScanLookPC软件，软件可以识别导入的POS数据，并可以使用POS数据解算出具有准确姿态位置的点云，满足后续作业处理要求。
  三、数据处理与分析
  （一）基于统计规律的点云去噪
  机载激光雷达点云是通过机载激光雷达扫描获得的，以离散、不规则方式分布在三维空间中的点的集合[10]。由于周围环境或者扫描仪本身特性，激光束在数据获取过程中有时会产生异常点，即噪声点。常见的噪声点主要包括高位粗差和低位粗差。高位粗差通常是因为在采集数据的过程中受到鸟类等低飞飞行物的影响，误将这些物体反射回来的信号当作被测目标的反射信号记录下来所产生的。低位粗差则是由于测量过程中的多路径误差或者扫描仪的误差导致产生的极低点，例如由于水渍导致产生低于路面的点。
  这两类噪声点以及路面上方的电力线等会影响后续处理，针对以上问题，根据点云数据的高程直方图统计规律，将高低噪声和路面上方的散乱点去除。首先将分段后的点云数据垂直投影到XOY平面，建立C列R行的规则格网,确定每个点的格网行列号；然后对每个格网内点进行高程直方图统计，绘制高程频数分布直方图，横轴为高程值，纵轴为单组高程范围的点数量。根据地物点云的高程连续性和频数差值特征，判断高程阈值Z1、Z2，保留该格网内高程值大于Z1且小于Z2的点。将高程较小或较大且数量较少的异常点删除，例如高低位噪点，路面正上方的电力线点、交通标志牌等。
  将原始点云数据进行去噪，得到的数据通过高程渲染。该数据为乡村区域的点云，地形复杂，起伏较大；林地分布散乱，且高度不均，既有多年的乔木，也有栽种不久的小树苗；同时存在低洼区域和房屋，以及零星的灯杆等。通过高程渲染，大部分树木和灯杆等被渲染为红色，道路和大棚区域被渲染为黄色区域，地面被渲染为绿色区域，低洼区域被渲染为蓝色。因此可通过高程区分得到林地范围。
  （二）分离地面点
  为减少不必要的数据量，首先采用CSF[11]方法将地面点分离，采用陡坡模式，分辨率设置为0.5m，距离阈值设置为0.3m，迭代次数设置为500，去除地面点之后的结果如图5所示。地面点被分类，树木点云仍然保留，减少了数据量，便于后期处理。
  （三）计算对比林地面积
  获取非地面点云之后，通过高程判断，利用点云处理软件在三维视图中勾画林地范围，利用软件自动计算单个多边形的面积，并将林地区域的多边形面积累加，即可得到准确的林地面积。
  为验证该技术的可行性和准确性，通过人工多次测量的手段取均值，统计不同区块的林地面积记为SP，通过软件在点云中多次计算取均值的林地面积记为S，并计算偏差P。机载LiDAR系统偏差符合工作要求，获取机载点云数据计算林地面积的测量计算方法减少了人工作业难度，提高了效率。
  四、结束语
  机载LiDAR系统是一种较为高效的新型测绘技术方法，该技术实现了测区三维点云数据的快速采集，数据精度高，测区覆盖全，对地貌地形具有三维实景还原能力。本文利用机载LiDAR系统获取某区域的机载点云数据，通过内业数据解算和后处理，通过高程判别林地范围，并利用点云处理软件在三维视图中确定林地范围并计算计算林地面积。通过实际工程试验验证了机载LiDAR技术在林地面积调查等领域具有较大优势，是一种较为完善的解决方案。