1、引　言

遥感技术已经成为植被生态系统探测的一种重要手段，为进行准确、高效的林业调查、作物估产发挥了重要作用，对地激光雷达(航空、星载)作为近年来快速发展的遥感手段，在进行植被结构参数获取、森林碳储量评估等方面的作用也日益突出(李增元 等，2016)。然而，在生物物种丰富的地区，自然植被群落通常呈现乔—灌—草立体分布，不同高度处光照、水分和养分有较大差异。因此，探测3维立体分布特性不仅是精确的生态系统监测和管理的需要，而且对病虫害早期预警、预估生态系统演化等研究也具有重要意义(王纪华 等，2004；Asner 等，2014)，但是现有的遥感手段在解决这些问题时遇到了较大的困难。

对于垂直复杂的植被分布状况，被动光学传感器接收的遥感信息通常包含了不同高度的贡献，实际上是3维立体分布状况下的2维累加，在实际成像的过程中将3维立体目标2维化，给植被参数反演带来较大的不确定性(阎广建 等，2002)。虽然多角度观测能够获取植被各个角度的光学信息，进而获取叶绿素等生化组分的垂直分布信息(赵春江 等，2006)，但是在航空和卫星尺度上很难获取“热点”和“冷点”等关键观测角度的信息，难以直接应用到卫星尺度上进行垂直分层分布反演。激光雷达虽然已经广泛应用于植被结构参数提取、电力巡线、滑坡监测等各个方面(Yang 等，2015；Gao 等，2015；庞勇 等，2005；Mallet和Bretar，2009；刘正军 等，2014)，但是主要利用了其测距特征，受到单一波段限制对光谱特征的挖掘较少。因此，将上述两种传感器进行融合受到广泛关注(Kulawardhana 等，2014；Ghosh 等，2014；Dalponte 等，2014)，但研究发现两种传感器协同无论是数据前期获取，还是后处理都存在很多缺点和误差，为进一步结合这两种传感器的优势，更直接的方法是利用一种仪器设备，直接采集带有光谱属性信息的激光雷达点云或波形数据，进行植被结构参数和生化参数提取。

高光谱激光雷达结合了被动光学的高光谱观测能力及激光雷达的垂直探测特点(刘清旺 等，2012；Hopkinson 等，2016)，具有探测冠层内部或者底部的精细结构和光谱的能力(Wing 等，2012)，因此具有较大的优点，它不依赖于太阳光，不受观测几何信息的影响；通过调整光束的出射，可以形成精确的点观测能力，减少混合像元的效应。目前已经有多个国家的实验室开展了高光谱激光雷达的研制和应用工作，如Morsdorf等人(2009)曾经利用英国爱丁堡大学4波段仪器开展了植被生理参数信息的垂直分布诊断，对利用生态过程模型模拟的不同年龄的树木进行了诊断，研究了不同树龄对植被NDVI，PRI剖面的影响。Suomalainen等人(2011)利用芬兰大地测量研究所8波段仪器研究发现，光谱的垂直分布也可以用于目标的自动探测和分类研究，垂直几何和光谱信息可以通过一次观测同时获取。Hakala等人(2012)在室内对砍伐的云松研究表明，利用不同的波段组合不仅可以对冠层枝、叶进行区分，利用其叶可以探测出冠层不同高度的叶片含水率、叶绿素等生化组分的差异。Rall等人(2004)基于2波段仪器首先提出了主动植被指数的概念，并通过低功率的激光二极管实现了水平方向的落叶植被叶片的观测，在680 nm附近的波段可以检测出植被叶绿素含量的特征变化。Gong等人(2012)利用其设计的4波段仪器研究表明其可以捕捉到精细的叶片生化组分浓度的变化。李旺等人(2014)利用中国科学院遥感与数字地球研究所研制的高光谱激光雷达样机开展了红花羊蹄甲、紫薇等两种叶片叶绿素、胡萝卜素、氮素的提取试验，与实验室化验分析比较表明其提取决策系数分别达到0.85, 0.71, 0.51。

上述研究表明，高光谱激光雷达仪器结合了高光谱和激光雷达两种传感器的特点，可以探测植被生理结构和光合作用直接相关的生理过程，在复杂地表植被特征探测过程中具有极强的科学价值和应用潜力。但是针对植被生化组分垂直分布反演方面，目前还没有基于高光谱激光雷达特定方法流程。本研究基于高光谱激光雷达开展了植被生化组分垂直探测研究，对生化组分导致的不同波段的激光反射信号的差异进行了研究，提出基于该仪器的生化组分垂直分布反演方法，并提出适应于仪器特点的定标和处理方法。

2、仪器试验与数据处理     (2.1) 高光谱激光雷达仪器

本研究利用了遥感科学国家重点实验室自行研制的32波段高光谱激光雷达样机(Sun 等，2014)，样机覆盖范围从可见光到近红外区域(图1)。该设备主要由2维扫描平台、超连续谱脉冲激光光源、固定于扫描平台的同轴发射接收系统、多通道全波形测量装置、控制中心等组成(Niu 等，2015)，其中扫描系统采用了2维转台方案PTU-D48 E(FLIR Systems Inc)，设计的水平和垂直的最小角度步长分别为0.026°和0.013°，对应10 m距离处的水平扫描位移约为4—5 mm。扫描系统与主控系统采用串口连接，遵循转台的通信协议。

超连续谱脉冲激光光源(NKT SuperK Compact)的脉冲宽度为1—2 ns，峰值功率约20 KW，脉冲频率为20—40 KHz。超连续谱激光器所发出的“白光”激光经过准直器后变成汇聚的光束(发散角小于5 mrad)，该光束经过两个反射镜(M1，M2)产生两次折射，进入望远镜光轴并发射出去(图2)。其中反射镜M1为微透反射镜，激光束中绝大部分光反射到镜筒内，少部分光透过反射镜再经过光纤投射到发射探测传感器APD1中，实现发射光束的采集，采集的信号用来触发一次测量。被探测目标的散射光经过消色差折射式望远镜(焦距400 mm，口径80 mm)进行收集，在焦点处由光纤进行收集，然后通过光栅光谱仪进行分光，产生32个独立的波段，分光后的出射光投射到相应的线性探测器(PA)上，同时实现4个波段的回波探测。探测器采用线阵PMT光电倍增阵列(Hamamatsu Photon Techniques Inc)，敏感波长为400—1000 nm，上升时间为0.6 ns。光电探测器输出信号通过高速示波器(DPO5204B，Tektronix Inc)进行采样和存储，生产高分辨率的波形回波。扫描过程由控制主机控制扫描云台产生X，Y两轴的运动。

原型系统可以通过对各个通道记录波形的后处理实现飞行时间和回波的光谱的测量。两轴转台带动光学部分进行扫描，可以获取的一系列的不同波长的回波信号。

    (2.2) 试验数据处理

试验以火炬花(Kniphofia uvaria)为扫描对象，火炬花属多年生草本植物，宜生长于疏松肥沃的沙壤土中，多丛植于草坪之中或植于假山石旁，用作配景。试验中高光谱激光雷达在距离其顶部7 m距离处进行扫描。

试验中，高光谱激光雷达扫描的水平步长为80，垂直步长为160，因此共扫描12800个点。在扫描过程中，首先对准扫描目标的中心，转台运动到扫描范围的左上角；其次向右转动，转动的角度为一个水平步长角度，每转完一个步长，要测量32次；再次，依次重复到指定位置后，向下转动一个倾斜步长的角度；然后再向左转动到指定位置后停止，依次重复最终完成面扫描。

高光谱激光雷达对地物扫描，在每个点记录32个波段回波波形信息，典型的回波波形如图3所示。

对于记录的32个波段的回波波形文件，首先进行数据清洗，包括将文件进行重命名，删除重复数据等，然后根据信号出现频率进行噪声的去除，继而进行波段分割获取每个波段的信息，并单独进行保存。对于每个波段的回波波形信号，通过发射波形和接收波形能量分布，确定发射和接收信号的阈值范围，如统计所有发射波形，对发射波形强度最大值出现的位置进行统计，从而得到近似高斯分布，以3倍的Sigma为阈值确定发射波形位置的阈值范围。根据发射波形和接收波形的阈值范围，提取出阈值距离范围的信号，通过高斯分解的最大值作为该段波形的强度并记录其位置，并利用接收强度除以发射强度进行发射强度的校正，最后依次获得32个波段的在某个位置处的校正强度值，并以点云的方式进行存储，最终通过算法提取各种参数，典型的处理流程如图4所示。

    (2.3) 高光谱激光雷达点云

根据激光雷达仪器扫描结果，经过基本数据处理，获取了高光谱激光雷达点云数据，图5为点云数据高程分布情况，顶部为距离仪器较近的位置，底部为距离较远的区域，从图5中可以看出高光谱激光雷达在40 cm范围内能够准确区分叶片的不同高度。

根据2.2节介绍的数据处理方法，对32个波段的激光雷达回波波形进行拟合，以高斯拟合的峰值作为该波段在该点处的回波强度值，将此回波强度值除以发射强度值从而得到32个波段的光谱定标回波强度信息，图6分别是800 nm，670 nm，572 nm和523 nm处激光雷达回波强度点云，从图6中可以看出，红色和绿色火炬花叶片对800 nm激光都具有较强的反射，红色叶片对670 nm激光具有较强的反射，而绿色叶片则反射较弱，572 nm和523 nm的激光点云也在绿色叶片处存在明显的差别。

3、方法与结果     (3.1) 生化组分化验值与高光谱激光雷达关系

扫描完成后，将叶片进行破坏性采样，并在实验室中采用萃取法进行化验，图7为火炬花不同位置叶片的高光谱激光雷达回波强度曲线。对采样样本试验室内化验得到叶片叶绿素a、叶绿素b、类胡萝卜素、凯氏氮、水分等生化组分含量，将测量生化组分与激光雷达光谱建立关系，在研究中主要利用NDVI和PRI与植被生化组分建立关系，建模总样本数为18个，包括9个红色叶片和9个绿色叶片样本，相关关系如表1。

${\rm{NDVI}} = \left( {{\rho _{800}} - {\rho _{670}}} \right)/\left( {{\rho _{800}} + {\rho _{670}}} \right)$

(1)

${\rm{PRI}} = \left( {{\rho _{572}} - {\rho _{523}}} \right)/\left( {{\rho _{572}} + {\rho _{523}}} \right)$

(2)

从图7中可以看出，红色叶片与绿色叶片的高光谱激光雷达回波强度曲线有明显的不同，红色叶片670 nm附近红光波段具有较强的反射，而绿色叶片在此波段区域存在明显的波谷。研究表明波谷的存在是由于叶绿素在此波段有较强的吸收作用，因此从图7中可以看出高光谱激光雷达强度曲线能够反映植被叶绿素含量的大小。

    (3.2) 高光谱激光雷达植被指数垂直分布

对点云的光谱反射强度进行处理，得到激光雷达植被指数点云数据，图8是获取的NDVI点云和PRI点云。基于植被指数点云数据，计算了植被指数的垂直分布图，图8(b)和图8(d)为对应的植被指数沿着高度分布的曲线，NDVI明显在顶部红色叶片处(7.1 m)的值偏低，而在中部7.2 m左右值较高。PRI植被指数则在垂直方向上变化不明显。

    (3.3) 生化组分垂直分布

基于高光谱激光雷达植被指数与生化组分的关系(表1)，可以将植被指数点云根据公式进行计算从而得到各种生化参数点云。本研究计算了火炬花叶绿素a和类胡萝卜素(单位：mg/g)生化组分含量沿高度的变化，如图9所示。从图9可看出，在顶部(7.1 m)处，叶片的生化组分含量较低，叶绿素a普遍低于0.5，胡萝卜素低于0.2，而在中部叶片处(7.2 m)，生化组分含量较高，这与植被红色叶片与绿色叶片在垂直方向上分布一致，即顶部的红色叶片叶绿素、胡萝卜素含量低，而中部的绿色叶片生化参数含量高。

4、结　论

本文利用研究室自行研制的对地高光谱激光雷达仪器，以火炬花为例开展了室内扫描，针对获取的含有32个波段的激光雷达点云数据提出了数据处理的基本流程。然后根据生化组分与光谱的关系，获取生化组分(叶绿素a，胡萝卜素)在3维空间的垂直分布，与实验室化验结果具有较好的一致性，研究表明新型的对地高光谱激光雷达仪器可以为精确的植被3维建模和定量遥感反演提供支持。

高光谱激光雷达仪器可以详细测量不同位置处的光谱，可以对不同位置处的生化组分进行详细的测定，避免了结构参数对生化组分反演的影响，相比于多角度被动光学观测能够提供更有效的观测手段。相比于被动的高光谱遥感，该仪器能够获取不同位置的高光谱信息，从而避免不同高度处光谱信息的干扰，相比于激光雷达传感器，能够提供丰富的光谱信息，从而可以与地物属性建立联系。

由于实验室化验样本量较少，生化组分提取方面，本研究只利用了单一的植被指数与生化组分建立关系，没有完全利用32个波段的反射强度信息进行生化组分垂直分布反演算法的深入研究。同时，由于超连续谱激光器的光源强度较弱，因此其探测距离较短，造成其回波信号较弱，给波形分解及后续的点云生成带来较大的难度，限制了其大范围的应用。但是研究表明，随着硬件水平的提高和更多试验的开展，高光谱激光雷达仪器可以获得更多更高质量的数据，在植被生理生化特征探测方面具有极大的应用前景。

志 谢 此次试验样本化验得到了北京农业质量标准与检测技术研究中心的大力协助，在此表示衷心的感谢！