遥感地学分析-地物光谱特征与遥感数字图像信息提取.ppt

第3章 遥感图像识别与信息提取遥感地学分析遥感地学分析课程小组编写课程小组编写GeographyAnalysisforRemoteSensing遥感地学分析遥感地学分析课程编写小组内容提要n3.1 遥感图像地物特征与识别q3.1.1遥感图像地物特征q3.1.2典型地物的反射光谱特征n3.2遥感图像的目视解译q3.2.1目视解译标志q3.2.2目视解译方法与步骤n3.3遥感图像的计算处理q3.3.1遥感数字图像的概念q3.3.2遥感数字图像预处理q3.3.3遥感数字图像分类处理q3.3.4遥感数字图像定量反演遥感地学分析课程编写小组3.1 遥感图像地物特征与识别3.1.1遥感图像地物特征q地物的反射光谱特性q地物的发射光谱特性q地物的透射光谱特性遥感地学分析课程编写小组3.1.1 遥感图像地物特征在可见光与近红外波段，地表物体自身的辐射几乎等于零。地物发出的波谱主要以反射太阳辐射为主。太阳辐射到达地面之后，物体除了反射作用外，还有对电磁辐射的吸收作用。电磁辐射未被吸收和反射的其余部分则是透过的部分，即：到达地面的太阳辐射能量反射能量吸收能量透射能量一般而言，绝大多数物体对可见光都不具备透射能力，而有些物体如水，对一定波长的电磁波透射能力较强，特别是对0.450.56m的蓝绿光波段，一般水体的透射深度可达1020m，清澈水体可达100m的深度。对于一般不能透过可见光的地面物体，波长5cm的电磁波却有透射能力，如超长波的透射能力就很强，可以透过地面岩石和土壤。遥感地学分析课程编写小组3.1.1 遥感图像地物特征l1、地物的反射光谱特性u反射率n地物的反射能量Pe占总入射能量Po的百分比，称为反射率u反射类型n镜面反射（Specularreflection）入射波与反射波在同一平面内，入射角与反射角相等时，所形成的反射现象n漫反射（Diffusereflection）无论入射波方向如何，其反射波分散到各个方向遥感地学分析课程编写小组3.1.1 遥感图像地物特征l方向反射：朗伯体表面实际上是一个理想化的表面。它被假定为介质是均匀的、各向同性的，并在遥感中多用以作为近似的自然表面。电磁波的反射的三种形式遥感地学分析课程编写小组3.1.1 遥感图像地物特征l2、地物的发射光谱特性u发射率地物发射电磁辐射的能力，以黑体辐射作为基准指单位面积上观测地物发射的某一波长的辐射通量密度；指与观测地物同温度下黑体的辐射通量密度；遥感地学分析课程编写小组3.1.1 遥感图像地物特征通常，根据发射率与波长的关系，将地物分为三种类型u黑体，其发射率=1，即黑体发射率对所有波长都是一个常数，并且等于1。u灰体，其发射率=常数1。即灰体的发射率始终小于1，发射率不随波长变化。u选择性辐射体，其发射率1，发射率随波长而变化。地物的发射波谱特性曲线遥感地学分析课程编写小组3.1.1 遥感图像地物特征太太太太阳阳阳阳辐辐辐辐射射射射遥感地学分析课程编写小组3.1.1 遥感图像地物特征l3、地物的透射光谱特性u透射率即地物透射的能量与入射总能量的百分率，称之为投射率=E/E100%遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征1、岩石的反射光谱特征 岩石的波谱特征是地质遥感的基础，不同的矿物成分、矿物含量、风化程度、含水状况、颗粒大小、表面的光滑程度、色泽等都会影响到其反射波谱特征。遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征2、土壤的反射光谱特征 自然状况的土壤表面的反射率没有明显的峰值和谷值，一般来说土质越细，反射率越高，有机质含量越高和含水量越高反射率越低。此外土壤的肥力也会对反射率产生影响。遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征不同含水量的土壤反射光谱曲线三种土壤反射波谱曲线比较遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征3、水体的反射光谱特征 水体的反射主要在蓝光波段，其他波段吸收都很强，特别在近红外以后水体便成为一个吸收体。光谱反射特性可能包括来自三方面的贡献：水的表面反射、水体底部物质的反射和水中悬浮物质的反射。光谱吸收和透射特性不仅与水体本身的性质有关，而且还明显地受到水中各种类型和大小的物质有机物和无机物的影响。遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征水体的反射光谱特征遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征4、植被的反射光谱特征 植被的反射波谱曲线规律性明显而独特，主要分为三个波段：可见光波段（0.40.76 m)有一个小的反射峰，位于0.55 m(绿光波段）处，两侧蓝光波段（0.45 m）和红光波段(0.67 m)则有两个吸收带；近红外波段（0.70.8 m）有一反射徒坡，至1.1 m 附近有一峰值；中红外波段（1.32.5 m）受绿色植物含水量的影响，吸收率大增，反射率下降，在1.45、1.95和2.7 m 为中心是水的吸收带，开成低谷。遥感地学分析课程编写小组3.1.2 典型地物的反射光谱特征矿区红杉林反射曲线的蓝移现象绿色植物有效光谱响应特征遥感地学分析课程编写小组3.2 遥感图像目视解译目视解译是用肉眼或借助于简单的工具如放大镜、立体镜、投影观察器等，直接由肉眼来识别图像特性，从而提取有用信息，即人把物体与图像联系起来的过程。具备的基本知识：q专业知识、q地理区域知识、q遥感系统知识。遥感地学分析课程编写小组3.2.1 目视解译标志 1、直接解译标志、直接解译标志q色调（tone）：色调是识别目标地物的基本依据。由于地物属性的差别，在遥感图像上表现出色调上的差别。一般来说由于人眼的局限性，在图像可分出16个灰度级。q颜色（colour）：由于目标地物在不同波段中反射或发射电磁辐射能量的差异性，由此而表现在彩色图像颜色的差异性。颜色也是目视解译的基本标志之一。q图型（pattern）：目标地物有规律排列而成的图形结构。遥感地学分析课程编写小组3.2.1 目视解译标志q阴影（shadow）：阴影是遥感图像上由于电磁辐能量被遮挡而产生的辐射能量减弱。由于阴影的存在，可据此地物的性质或高度，应注意的是阴影的形状与大小受到辐射能量入射角的影响。q形状（shape）：是指目标地物在遥感影像上呈现的外部轮廓。由于不同地物的顶视平面的差异，可据此判断目标地物的性质。q纹理（texture）：是指遥感图像中目标地物内部色调有规则变化造成的影影结构。q大小（size）：指遥感图像上目标地物的形状、面积与体积的度量关系。是遥感图像上测量目标地物最重要的数量特征之一。同时根据其大小可以推断地物的属性。遥感地学分析课程编写小组3.2.1 目视解译标志2、间接解译标志间接解译标志q位置（site）：指目标地物分布的地点。任何地物与其周围地理环境（geographical Environment）或生态环境(ecological Environment)总是存在着一定的空间关系，并受到周围环境的某种程度的制约。q相关布局（association）：是指多个目标地物间的空间配置关系。遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤 1、目视解译的认知过程u自下向上过程q图像信息获取q特征提取q识别证据选取u自上向下过程q特征匹配q提出假设q图像辨识遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤l2、图像解译方法q遥感资料的选择及影像处理u遥感资料的选择q资料类型选择q波段选择q时间选择q比例尺选择u遥感图像的处理q影像放大q影像数字化q图像处理遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤 目视解译的方法目视解译的方法 u直判法 是根据遥感影像目视判读直接标志，直接确定目标地物属性与范围的一种方法。u对比分析法对比分析法包括同类地物对比分析法、空间对比分析法和时相动态对比法。不仅是同一遥感影像图进行对比，而且要借助不同时相的遥感影像图之间进行对比。u综合推理法综合考虑遥感图像多种解译特征，分析、推理某种目标地物的方法。遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤u信息复合法利用专题图和地形图等信息与遥感影像图重合，根据这些辅助信息识别遥感影像图上目标地物的方法。u地理相关分析法根据地理环境中各种地理要素之间的相互依存，相互制约的关系，借助专业知识，分析和推理某种地理要素性质、类型、状况与分布的方法。遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤 3、解译步骤解译步骤u准备工作包括资料收集、分析、整理和处理u初步解译、建立解译标志包括路线路勘，制订解译对象的专业分类系统和建立解译标志u室内解译u野外验证包括解译结果校核检查，样品采集和调绘补测u成果整理包括编绘成图，资料整理和文字总结遥感地学分析课程编写小组3.2.2 遥感图像解译方法与步骤目目视视解解译译步步骤骤遥感地学分析课程编写小组3.3 遥感数字图像信息提取3.3.1遥感数字图像的概念q图象、数字图象q遥感数字图象3.3.2遥感数字图像处理q遥感数字图像的获取q遥感数字图像预处理q遥感数字图像的变换、增强和融合q遥感数字图像分析遥感地学分析课程编写小组3.3.1 遥感数字图像的概念l图象、数字图象物理世界中客观对象的一种表示u数字图像客体或可见图像的数字表述。它实际上是具有某种数值的一些点按行(横)和列(纵)排成的二维矩阵u模拟图像指空间坐标和明暗程度都连续变化的、计算机无法直接处理的图像遥感地学分析课程编写小组3.3.1 遥感数字图像的概念l遥感数字图象指以遥感方式获得的以数字形式表述的遥感影像遥感数字图像最基本的单位是像素像素的属性特征常用灰度值来表示，即该像素位置上亮暗程度的整数值遥感地学分析课程编写小组3.3.2 遥感数字图像处理l一、遥感数字图像的获取一、遥感数字图像的获取 通过接收、记录目标物的电磁波特征的仪器，即传感器获得的数字图像获取过程遥感地学分析课程编写小组3.3.2 遥感数字图像处理二、遥感数字图像预处理u1、辐射校正消除图像数据中依附在辐射亮度里的各种失真的过程称为辐射校正。完整的辐射校正包括遥感器校正、大气校正，以及太阳高度和地形校正。遥感器纠正：遥感器的设计大气辐射纠正：地形辐射纠正：需要DEM地物反射模型纠正：需要和成像时刻取得同步的地面地物光谱测量数据遥感地学分析课程编写小组3.3.2 遥感数字图像处理大气纠正方法以红外波段最低值校正可见光波段回归法相对散射模型遥感地学分析课程编写小组3.3.2 遥感数字图像处理以红外波段最低值校正可见光波段q前提假设：大气散射的影响主要在短波波段，红外波段中清洁的水体几乎不受影响，反射率值应当为0。由于散射影响，而使得水体的反射率不等于0，推定是由于受到了天空辐射项的影响。q直方图法确定q纠正方法：差值法遥感地学分析课程编写小组3.3.2 遥感数字图像处理回归法l原理l选择可见光和红外波段进行2维散点图，建立线性回归方程。遥感地学分析课程编写小组3.3.2 遥感数字图像处理相对散射模型qSTEP1:根据某个可见光波段的直方图选出黑暗地物的初始灰雾值；qSTEP2:根据此灰雾值的幅度确定大气类型（选择合适的散射模型）；qSTEP3:根据模型和初始灰雾值预测其它波段的灰雾值；qSTEP4:对每个波段进行大气纠正。遥感地学分析课程编写小组3.3.2 遥感数字图像处理地形辐射纠正q需要DEMq简单的处理方法比值法：有效消除阴影的影响。遥感地学分析课程编写小组3.3.2 遥感数字图像处理2、几何校正造成几何位置的畸变有原因：q遥感器本身引起的畸变q外部因素引起的畸变q处理过程中引起的畸变遥感地学分析课程编写小组3.3.2 遥感数字图像处理遥感器本身引起的几何畸变与遥感器的结构、特性和工作方式不同而异。这些因素主要包括：1）透镜的辐射方向畸变像差；2）透镜的切线方向畸变像差；3）透镜的焦距误差；4）透镜的光轴与投影面不正交；5）图像的投影面非平面；6）探测元件排列不整齐；7）采样速率的变化；8）采样时刻的偏差;9)扫描镜的扫描速度变化。n n 遥感器本身引起的畸变遥感器本身引起的畸变遥感地学分析课程编写小组3.3.2 遥感数字图像处理影响图像变形的外部因素包括：1）地球的曲率2）大气密度差引起的折光3）地形起伏4）地球自传5）遥感器轨道位置和姿态等n n 外部因素引起的畸变外部因素引起的畸变遥感地学分析课程编写小组3.3.2 遥感数字图像处理遥感器轨道位置和姿态引起的误差遥感器轨道位置和姿态引起的误差地球自传引起的误差遥感地学分析课程编写小组3.3.2 遥感数字图像处理地球曲率和地形起伏引起的误差地球曲率和地形起伏引起的误差遥感器轨道位置和姿态引起的误差遥感器轨道位置和姿态引起的误差遥感地学分析课程编写小组3.3.2 遥感数字图像处理 遥感图像再处理过程中产生的误差，主要是由于处理设备产生的噪声引起的。q传输、复制q光学q数字n n 处理过程中引起的畸变处理过程中引起的畸变遥感地学分析课程编写小组3.3.2 遥感数字图像处理n遥感图像的几何纠正方法遥感图像的几何纠正方法q遥感图象的几何粗处理和精处理q遥感图像的几何纠正按照处理方式分为光学纠正和数字纠正q光学纠正主要用于早期的遥感图像的处理中，现在的应用已经不多。除了对框幅式的航空照片（中心投影）可以进行比较严密的纠正以外，对于大多数动态获得的遥感影像只能进行近似的纠正遥感地学分析课程编写小组3.3.2 遥感数字图像处理l数字图象几何纠正：通过计算机对离散结构的数字图像中的每一个像元逐个进行纠正处理的方法。l几何精校正是指利用地面控制点使遥感图像的几何位置符合某种地理系统，与地图配准，并调整亮度值。也就是在遥感图像的像元与地面实际位置之间建立数学关系，将畸变图像空间中的全部像元转换到校正图像空间去。l内容一般包括两个方面：一个是图像像元空间位置的变换，另一个是像元灰度值的重采样。遥感地学分析课程编写小组3.3.2 遥感数字图像处理n几何精校正具体步骤：q选取地面控制点（GCP）q选择空间变换函数q重采样和内插最邻近内插法双线性内插三次卷积内插遥感地学分析课程编写小组3.3.2 遥感数字图像处理数字图象的纠正过程数字图象的纠正过程纠正的函数可有多种选择：多项式方法、共线方程方法、随机场内插方法等等。其中多项式方法的应用最为普遍。遥感地学分析课程编写小组3.3.2 遥感数字图像处理l三、遥感数字图像的变换、增强和融合三、遥感数字图像的变换、增强和融合 1.遥感数字图像变换u概念：图像变换指的是将图像从空间域转换到变换域例如频率域的过程。进行图像变换的目的就是为了使图像的处理过程简化。u作用：通过图像变换简单而有效的实现增强处理通过图像变换可以对图像进行特征抽取。遥感地学分析课程编写小组3.3.2 遥感数字图像处理快速傅立叶变换q目的进行数据压缩、图像的增强、特征提取q方法信号处理中的频率域分析方法q步骤选择适当的变换函数进行傅里叶变换分析变换的结果进行傅里叶逆变换遥感地学分析课程编写小组3.3.2 遥感数字图像处理KL变换q目的n减少图像波段之间的相关性，去除多余的信息，减少图像的数据量q方法n统计学中的正交变换方法q性质和特点q步骤n进行数据统计n进行主成分分析n进行旋转变换和逆变换遥感地学分析课程编写小组3.3.2 遥感数字图像处理KT变换q目的n分离和消除干扰信息n突出研究的专题信息q方法n几何中的坐标旋转方法q性质和特点q步骤n作两个波段的散点图n分析灰度值的变换特点遥感地学分析课程编写小组3.3.2 遥感数字图像处理2.遥感数字图像增强l1、反差增强图像反差增强又称对比度增强u灰度拉伸根据原图像的直方图确定需要做拉伸变换的灰度值区间，然后把这一（或一些）灰度值区间按某种直线或曲线方程关系拉伸或压缩而成为变换后的灰度值区间直方图线性拉伸直方图分段线性拉伸非线性拉伸变换遥感地学分析课程编写小组3.3.2 遥感数字图像处理u直方图均衡化q变换后的直方图接近均匀分布。即图象中每一灰度级的像元数目大致相同。q使得面积最大的地物细节得以增强，而面积小的地物与其灰度接近的地物进行合并，形成综合地物。减少灰度等级换取对比度的增大。原直方图：均衡后直方图：遥感地学分析课程编写小组原图及其直方图均衡后的图和直方图遥感地学分析课程编写小组3.3.2 遥感数字图像处理u直方图匹配q将原始图象转换为给定直方图的图象q各种拉伸以后虽然对于感兴趣的地物提高了分辨能力，但同时也造成了信息损失。遥感地学分析课程编写小组3.3.2 遥感数字图像处理l2、空间域图像增强空间域图像增强处理是应用某种数学模式直接改变图像像元灰度值的变换。这种变换与像元的坐标无关，只改变像元点的灰度深浅程度。即改变了图像的对比度，这样图像中的某些信息被突出（可能也有些信息被压抑了）达到了图像增强的目的。可分为单点处理和邻域处理遥感地学分析课程编写小组3.3.2 遥感数字图像处理u边缘增强突出目标的轮廓或边缘信息，主要是通过微分法实现的梯度法和拉普拉斯算子法u平滑滤波低频增强的空间域滤波技术。它可以滤掉由于孤立的单点噪声而引起的灰度偏差u中值滤波u定向滤波遥感地学分析课程编写小组3.3.2 遥感数字图像处理l3、频率域图像增强高通滤波（锐化）低通滤波（平滑）带通滤波（突出地物）同态滤波（改善图像质量）遥感地学分析课程编写小组3.3.2 遥感数字图像处理l4、代数运算增强q基本代数运算n差值运算两个波段的差值n比值运算两个波段的比值q应用遥感地学分析课程编写小组3.3.2 遥感数字图像处理u差值运算q利用不同地物之间光谱的特征q有利于突出目标与背景反差小的信息红外波段-红波段，突出植被实实例例B4-B3遥感地学分析课程编写小组3.3.2 遥感数字图像处理u比值运算q两个波段的比值q目的减轻地形的干扰q如果阴坡的坡度大于太阳高度角，则难于处理扩展地物，特别是植被间的光谱差异抑制大气的影响减少了数据之间的相关性遥感地学分析课程编写小组3.3.2 遥感数字图像处理l植被指数NDVITM4/TM3l其他林业上常用n4/2(3),5/1,7/3(2)n4/3(2),5/1,3/2(1)常用的比值遥感地学分析课程编写小组3.3.2 遥感数字图像处理二 应用实例NDVI遥感地学分析课程编写小组3.3.2 遥感数字图像处理TM图像的5/2突出水体遥感地学分析课程编写小组 3.3.2 遥感数字图像处理l5、彩色图像增强u假彩色密度分割将灰度按照指定的间隔分割为不同的级，对新的密度分级分别赋予不同的颜色遥感地学分析课程编写小组遥感地学分析课程编写小组3.3.2 遥感数字图像处理u彩色增强在一个输入像元的灰度级上分别进行三个独立的色变换在显示屏幕上产生一个彩色合成图像。从而达到明显的分辨和识别效果？遥感地学分析课程编写小组3.3.2 遥感数字图像处理u假彩色合成q目的综合不同波段的特征，突出研究对象的差异q原理q原则信息量最大相关性最小差异最大q方法信息量分析相关分析方差分析其他遥感地学分析课程编写小组3.3.2 遥感数字图像处理TM图像4-3-2合成实例遥感地学分析课程编写小组3.3.2 遥感数字图像处理l3.遥感数据融合u实质在统一地理坐标系中将对同一目标检测的多幅遥感图像数据采用一定的算法，生成一幅新的、更能有效表示该目标的图像信息u原理及过程预处理主要包括遥感影像的大气校正、辐射校正及空间配准数据融合根据融合目的和融合层次智能地选择合适的融合算法，将空间配准的遥感影像数据（或提取的图像特征或模式识别的属性说明）进行有机合成，得到目标的更准确表示或估计遥感地学分析课程编写小组3.3.2 遥感数字图像处理u数据融合分类q像元级融合像元级融合是一种最低水平的融合q特征级融合特征级融合是一种中等水平的融合q决策级融合决策级融合是最高水平的融合u数据融合方法IHS变换、代数法、Brovey变换、图像回归法、主成分变换（PCT）、小波变换等遥感地学分析课程编写小组3.3.2 遥感数字图像处理三级融合层次下的融合方法三级融合层次下的融合方法 像元级特征级决策级代数法熵法专家系统IHS变换表决法神经网络小波变换聚类分析Bayes估计Brovey变换Bayes估计模糊聚类法主成分变换神经网络法可靠性理论回归模型法加权平均法基于知识的融合法Kalman滤波法Dempater-shafer推理法Dempater-shafer推理法遥感地学分析课程编写小组TM741RadarSatTM和RadarSat复合遥感地学分析课程编写小组3.3.2 遥感数字图像处理l四、遥感数字图像分析四、遥感数字图像分析u图像理解从图像中提取有用信息的过程u包括图像分类图像分割特征分析u输出内容明确的数据或图形图像u特点以具体的知识和目的为基础，结果更接近于解释的成果和实际的应用遥感地学分析课程编写小组3.3.2 遥感数字图像处理图像信息获取信息处理特征提取分类 图像分类（模式识别）遥感地学分析课程编写小组3.3.2 遥感数字图像处理l1、遥感图像分类技术u原理n同类地物在相同的条件下（光照、地形等）应该具有相同或相似的光谱信息和空间信息特征。不同类的地物之间具有差异n同类地物的像元在数目较大时，其特征量的分布类型接近正态分布u方法n利用像素的灰度值之间的关系，综合考虑矢量数据，考虑属性数据。对像素的进行分类n以随机变量的统计分析为基础遥感地学分析课程编写小组3.3.2 遥感数字图像处理u特点多变量的图像分类分类在一定的比例尺上进行u原则多变量综合考虑在特征空间的位置可以用均值表示，离散程度可以用方差表示分类的实质是将特征空间分为若干子区域，每个子区域为一类遥感地学分析课程编写小组3.3.2 遥感数字图像处理u指标（特征）波段其他属性地形特征调查数据u统计量均值方差距离度量离差平方和遥感地学分析课程编写小组3.3.2 遥感数字图像处理l2、具体分类方法u非监督分类n最大似然法最大似然法n最小距离法最小距离法n遥感地学分析课程编写小组遥感地学分析课程编写小组3.3.2 遥感数字图像处理u非监督分类q特征空间识别法q系统聚类法q分裂法(isomix)q动态聚类法(k-means)遥感地学分析课程编写小组3.3.2 遥感数字图像处理n影响结果的主要因素选择的距离聚类方式n特点需要较少的人工参与结果不依赖于统计参数对于非正态分布的数据可以得到较好的结果不适合同谱异物的数据结果的稳定性较低结果可以作为监督分类的参考遥感地学分析课程编写小组3.3.2 遥感数字图像处理n非监督分类小结需要较少的人工参与不适合同谱异物的数据计算速度慢分类结果不依赖于统计参数q对于非正态分布的数据可以得到较好的结果稳定性较低可以作为监督分类的参考遥感地学分析课程编写小组3.3.2 遥感数字图像处理u数字图像分类新技术q人工神经网络分类法q模糊分类法q亚像元分类法(subpixelclassification)q其他分类法遥感地学分析课程编写小组3.4 遥感数字图像信息定量反演l遥感图像信息定量反演的原理l遥感图像信息定量反演方法l遥感地表参数反演遥感地学分析课程编写小组3.4.1 遥感图像信息定量反演的原理u遥感机理和各种前向模型q统计模型q物理模型q半经验模型u反演模型和反演策略的研究遥感地学分析课程编写小组3.4.2遥感图像信息定量反演方法u辐射传输模型（RT模型)u几何光学模型（GO模型)u几何光学辐射传翰混合模型(GORT混合模型)u计算机模拟模型遥感地学分析课程编写小组3.4.3 遥感地表参数反演u方向反射BBDF的反演u植被结构参数反演叶面积指数LAI叶倾角分布LAD