dataFolder=fullfile（tempdir，“森林数据”，filesep）；数据文件=数据文件夹+“forestData.laz”;%检查文件夹和数据文件是否已存在folderExists=存在（dataFolder，“dir”)；fileExists=exist（数据文件，“文件”）;%创建一个不存在的新文件夹如果~folderExists mkdir（dataFolder）；结束%提取不存在的航空数据文件如果~fileExists解压('forestData.zip', dataFolder);结束%从文件中读取LAZ数据lasReader=lasFileReader（数据文件）；%读取点云以及相应的扫描角度信息[ptCloud，pointAttributes]=readPointCloud（lasReader，“属性”,“ScanAngle”）;%可视化输入点云图pcshow（ptCloud.Location）标题(“输入点云”)

分割地面，提取非地面和地面点groundPtsIdx = segmentGroundSMRF (ptCloud);nonGroundPtCloud =选择(ptCloud ~ groundPtsIdx);groundPtCloud =选择(ptCloud groundPtsIdx);%分别用洋红色和绿色可视化非接地点和接地点图pcshowpair（nonGroundPtCloud、groundPtCloud）标题(“分割的非地面点和地面点”)

地面点=地面云层位置；%消除沿x轴和y轴的重复点[uniqueZ，uniqueXY]=groupsummary（groundPoints（：，3），groundPoints（：，1:2），@mean）；uniqueXY=[uniqueXY{：}]；%创建插值并使用它估计每个点的地面高程F=散射三极管（双（uniqueXY），双（uniqueZ），“自然的”）;estElevation = F(双(ptCloud.Location(: 1)),双(ptCloud.Location (:, 2)));按地面归一化标高normalizedPoints = ptCloud.Location;normalizedPoints(:，3) = normalizedPoints(:，3) - estElevation;%可视化标准化点图pcshow (normalizedPoints)标题(“具有标准化高程的点云”)

%将栅格大小设置为每像素10米，将截止高度设置为2米gridSize = 10;cutOffHeight = 2;leafAngDistribution = 0.5;提取森林指标[canopyCover，gapFraction，Leaf Area Index]=helperExtractForestMetrics（标准化点，...pointAttributes.ScanAngle、gridSize cutOffHeight leafAngDistribution);可视化森林指标hForestMetrics=图；axCC=子批次（2,2,1，父项=hForestMetrics）；axCC.位置=[0.05 0.51 0.4 0.4]；图像SC（canopyCover，父项=axCC）标题（axCC，“天篷盖”)轴关颜色映射（灰色）axGF=子批次（2,2,2，父批次=hForestMetrics）；axGF.位置=[0.55 0.51 0.4 0.4]；图像分割，“家长”axGF)标题(axGF“差距分数”)轴关colormap（灰色）axLAI=子批次（2,2，[3,4]，父批次=hForestMetrics）；axLAI.Position=[0.30.02 0.40.4]；imagesc（叶面积指数，父项=axLAI）标题（axLAI，“叶面积指数”)轴关彩色地图（灰色）

%将栅格大小设置为每像素0.5米gridRes=0.5；%生成CHMcanopyModel = pc2dem (pointCloud (normalizedPoints) gridRes, CornerFillMethod =“最大值”）;%将无效和负CHM值剪裁为零canopyModel（isnan（canopyModel）| canopyModel<0）=0；%执行高斯平滑去除噪声效果H=F特殊(“高斯”，[5]，1）；canopyModel=imfilter（canopyModel，H，“复制”,“一样”）;%可视化CHM图imagesc（canopyModel）标题(“树冠高度模型”)轴关彩色地图（灰色）

%设置minTreeHeight为5 mMintreehight=5；%探测树梢[treeTopRowId，treeTopColId]=帮助检测树顶（canopyModel，gridRes，minTreeHeight）；%想象树梢图显示亮度图像(canopyModel)在…上地块（树冠、树冠、树冠），“处方”，MarkerSize=3）标题(“具有检测树顶的CHM”)轴关彩色地图(“灰色”)

%分割单株树label2D = helperSegmentTrees (canopyModel treeTopRowId、treeTopColId minTreeHeight);%标识label2D和传输标签中每个点的行和列id每点百分比rowId=ceil（（ptCloud.Location（：，2）-ptCloud.YLimits（1））/gridRes）+1；colId=ceil（（ptCloud.Location（：，1）-ptCloud.XLimits（1））/gridRes）+1；ind=sub2ind（size（label2D），rowId，colId）；label3D=label2D（ind）；%提取有效标签和对应点validSegIds=label3D~=0；ptVeg=select（ptCloud，validSegIds）；veglabel3D=label3D（有效注册表）；%为每个标签指定颜色numColors=max（veglabel3D）；colorMap=randi（[0 255]，numColors，3）/255；labelColors=label2rgb（veglabel3D，colorMap，OutputFormat=“三胞胎”）;%可视化树段图pcshow (ptVeg.Location labelColors)标题(“个人树片段”)视图（2）

%提取树属性treeMetrics = helperExtractTreeMetrics (normalizedPoints label3D);%显示前5个树段度量高级惩教主任(树木科,5);；