一种融合树冠三维结构的城市绿辐射受益量计算模型的制作方法
【专利摘要】本发明公开了一种基于融合树冠三维结构的城市绿辐射受益量计算模型，构建过程包括如下步骤：步骤1)提取城市中不同单体树种的冠层结构信息；步骤2)采用经验关系法反演所述基准图中的叶面积指数；步骤3)结合冠层结构信息与叶面积指数反演结果，构建二者之间的回归模型；步骤4)将单株立木的树冠体积值与LAI反演量相结合，获得位于冠层不同高度面上的叶面积含量；步骤5)确定所述基准图中的植被水平蒸散梯度；步骤6)提取建筑物的平均高度、投影面积和体积等参数；步骤7)结合步骤4)、步骤5)、步骤6)中得到的三个基本输入参数，构建城市绿辐射受益量计算模型，获得绿辐射受益量分布图，以及不同高度层所对应的绿辐射受益量分布特征。
【专利说明】一种融合树冠三维结构的城市绿辐射受益量计算模型
【技术领域】
[0001] 本发明涉及一种融合树冠三维结构的城市绿辐射受益量计算模型，本发明能够更 好的描述对城市环境具有积极影响的环境因子在三维地域空间网格内的含量及分布规模， 有效地揭示城市绿地生态效益服务范围。
【背景技术】
[0002] 城市绿地作为依附于城市土地资源的一种重要的城市基础设施形态，在改善城市 生态环境、促进人体健康等方面起着不可或缺的积极作用，其生态意义与环境价值的重要 性也正逐渐被人们所关注。城市绿地的生态效益功能主要体现在改善城市气候圈循环（调 节气温、增加空气湿度）；净化空气（吸收二氧化碳、释放氧气）、保持水土（涵养水源）；滞 留烟尘；降噪隔音；保护物种多样性，提升城市生态系统抗干扰能力等功能。
[0003] 目前，国外城市绿地研究多集中于改善城市环境质量、保护生物多样性、提供休闲 娱乐和生态系统管理等生态效益方面。随着人们对城市绿地功能认识的提升，以及遥感数 据源的多样化，人们不再满足于简单城市植被制图，而是考虑从多角度、多维数开展城市绿 色空间度量。
[0004] Hofmann等通过对格网中不同类型的植被赋以一定的权重系数，来度量每个格网 的绿量大小。奥地利萨尔茨堡大学利用遥感影像提取了二维植被类型信息，融合进邻接建 筑物层高、直线距离等要素，为反映公民福祉与城市绿色生态提供了一种新方法，并认为城 市绿地的邻接程度是改善居民生活条件的重要因素。Gupta等对研究区进行土地利用分类 及建筑物高度信息，利用格网法分别计算单个格网中建筑物密度、高建筑物密度和绿地的 邻接程度等参量并赋权重系数，计算研究区的城市绿色空间指数。20世纪80年代以来，我 国一直以城市人均公共绿地面积、绿地率和绿化覆盖率这3项指标来指导我国城市绿地建 设。尽管这些指标能从数量上说明一些问题，却难以在空间上描述绿地景观的空间布局，从 而真正反映绿地对人民服务的有效性。
[0005] 绿地是实际存在三维空间的几何实体，其产生的生态效应及其对环境的动态调节 理应也是三维的，考虑到这一问题，本发明试图建立一种融合树冠三维结构的城市绿辐射 受益量计算模型，用于定量描述绿辐射受益量在指定高度层上的含量及分布，量化绿地的 生态效益服务范围，以期更好地指导城市绿地规划建设。
【发明内容】
[0006] 本发明提出一种融合树冠三维结构的城市绿辐射受益量计算模型，用于描述不同 缓冲范围内蒸腾效益因子在指定高度层上的含量及分布特征，揭示城市绿地生态效益服务 范围。
[0007] 本发明的目的通过以下技术步骤实现：
[0008] 步骤1)利用机载激光雷达数据，辅以高分辨率航空影像，提取城市中不同单体树 种的冠层结构信息，具体包括树顶位置、树高、树冠投影轮廓、冠径、冠高、树冠体积等参数， 绘制冠层结构信息基准图；
[0009] 步骤2)采用经验关系法反演所述基准图中的叶面积指数，得到叶面积指数反演 结果分布图；
[0010] 步骤3)结合冠层结构信息与叶面积指数反演结果，构建二者之间的回归模型，获 得基于冠层体积信息的叶面积指数反演模型；
[0011] 步骤4)将单株立木的树冠体积值与LAI反演量相结合，经由微积分计算，获得位 于冠层不同高度面上的叶面积含量，绘制植被生态效益源分布图；
[0012] 步骤5)根据植被水平蒸散影响的最大有效距离，以及不同方位上单位距离内植 被温湿度效应的变化梯度观测值，确定所述基准图中的植被水平蒸散梯度；
[0013] 步骤6)对LiDAR数据进行预处理，生成数字高度模型与植被掩膜图像，提取建筑 物的平均高度、投影面积和体积等参数，绘制建筑物高度分布图与三维建模结果；
[0014] 步骤7)基于植物生态场理论，结合步骤4)、步骤5)、步骤6)中得到的三个基本输 入参数，构建一种融合树冠三维结构的城市绿辐射受益量计算模型，获得绿辐射受益量分 布图，以及不同高度层所对应的绿辐射受益量分布特征。
[0015] 进一步所述步骤1)的具体方法为：
[0016] a)采用基于活动窗口的局部最大值搜索法进行树顶位置与树高信息的提取；b) 在树顶检测结果基础之上，利用伞骨算法识别树冠投影边界；c)根据树冠的形态特征，使 用基于伸缩圆与冠高比的树冠体积计算方法，引入冠层孔隙率修正系数，获得树冠体积计 算结果。
[0017] 进一步所述步骤2)的具体方法为：
[0018] a)对多光谱遥感影像进行波段计算，得到NDVI分布图；b)采用表征冠层辐射传 输过程的物理模型PR0SAIL模拟获得用于拟合模型参数的叶面积指数；C)利用一定数量的 LAI参考值与NDVI数据绘制散点图，计算得到经验参数，即获得LAI遥感估算模型。
[0019] 进一步所述步骤6)的具体计算过程为：
[0020] a)使用高程阈值从经过植被掩膜后的图像提取建筑物，得到粗略的建筑物提取结 果；b)采用中值滤波（Median Filter)消除边界"锯齿现象"，得到边缘整齐的建筑物轮廓； c)经由分区统计和几何计算，得到单个独立建筑物的高度、面积、体积等信息。
【专利附图】
【附图说明】
[0021] 图1为树冠三维结构提取结果；
[0022] 图2为叶面积指数反演结果；
[0023] 图3为树冠三维简化结构示意图；
[0024] 图4为城市建筑物高度分布图与三维建模结果；
[0025] 图5为总绿辐射受益量分布图；
[0026] 图6为I. 7m绿辐射受益量分布图；
[0027] 图7为12m绿辐射受益量分布图；
[0028] 图8为28. 5m绿辐射受益量分布图；
【具体实施方式】
[0029] 下面结合附图对本发明"一种融合树冠三维结构的城市绿辐射受益量计算模型" 作进一步说明。
[0030] 本发明试图建立一种模拟植被蒸腾扩散过程的计算模型，可用于从二维平面与垂 直结构角度描述绿辐射空间的数量及分布特征，界定城市绿度空间的有效生态服务范围， 为环境规划、城市绿化指导提供参考。
[0031] (一）单株立木冠层结构信息提取
[0032] 城市绿度空间由生态效益源和绿辐射空间两部分构成，其中生态效益源是指可以 进行能量转换、发挥生理功能和生态服务效益的客观体，在本实例中特指树冠。本发明首先 基于机载激光雷达数据，辅以高分辨率航空影像，提取单株立木冠层结构信息，具体实现步 骤如下：
[0033] 1)基于单次回波点云数据构建树冠高度模型，在数字地面模型的基础上滤除非植 被信息，获得包含仅植被分布及高度信息的基准图；
[0034] 2)基于活动窗口的局部最大值搜索法，从所述基准图中提取单株立木的顶点（即 局部最商点）；
[0035] 3)在树顶检测结果基础之上，采用"伞骨算法"进行树冠投影边界识别，分离单体 树木并确定各自树冠的边界位置，得到树冠边界识别结果；
[0036] 4)结合树顶检测结果与树冠边界识别结果，基于伸缩圆与冠高比的树冠体积计算 方法生成单体树冠体积分布图；
[0037] 5)通过计算图像平均梯度反演冠层孔隙率，修正因树冠内部枝叶间存在的大量空 隙而造成的计算误差，进一步提升步骤4)中的树冠体积反演精度。
[0038] 经由上述步骤得到树冠三维结构提取结果，如图1所示。
[0039](二）叶面积指数反演
[0040] 采用经验关系方法反演植被指数。植被指数的经验关系法需要输入3个关键参 数：植被指数、经验关系形式和用于模型参数拟合的真实LAI数据。鉴于归一化差值植被指 数（NDVI)为最常用的植被指数，且比值形式的NDVI在突出植被信息的同时减小了冠层阴 影、土壤背景、大气污染和角度效益等影响。故本发明以NDVI为自变量，LAI为因变量，建 立估算模型：
【权利要求】
1. 一种融合树冠H维结构的城市绿福射受益量计算模型，该模型的构建方法包 括如下步骤： 步骤1)获取机载激光雷达数据，结合高分辨率航空影像，提取城市中不同单体树种的 冠层结构信息，具体包括树顶位置、树高、树冠投影轮廓、冠径、冠高、树冠体积参数，绘制冠 层结构信息基准图； 步骤2)采用经验关系法反演所述基准图中的叶面积指数，得到叶面积指数反演结果 分布图； 步骤3)结合冠层结构信息与叶面积指数反演结果，构建二者之间的回归模型，获得基 于冠层体积信息的叶面积指数反演模型； 步骤4)将单株立木的树冠体积值与LAI反演量相结合，经由微积分计算，获得位于冠 层不同高度面上的叶面积含量，绘制植被生态效益源分布图； 步骤5)根据植被水平蒸散影响的最大有效距离，W及不同方位上单位距离内植被温 湿度效应的变化梯度观测值，确定所述基准图中的植被水平蒸散梯度； 步骤6)对LiDAR数据进行预处理，生成数字高度模型与植被掩膜图像，提取建筑物的 平均高度、投影面积和体积等数，绘制建筑物高度分布图与H维建模结果； 步骤7)基于植物生态场理论，结合步骤4)、步骤5)、步骤6)中得到的H个基本输入参 数，构建一种融合树冠H维结构的城市绿福射受益量计算模型，获得绿福射受益量分布图， W及不同高度层所对应的绿福射受益量分布特征。
2. 如权利要求1所述的方法，其特征在于，所述步骤1)的具体方法为： a)采用基于活动窗口的局部最大值搜索法进行树顶位置与树高信息的提取；b)在树 顶检测结果基础之上，利用伞骨算法识别树冠投影边界（将伞滑算法识别树冠投影边界的 具体内容进行补充，并加入到说明书的相应部分）；c)根据树冠的形态特征，使用基于伸缩 圆与冠高比的树冠体积计算方法，引入冠层孔隙率修正系数，获得树冠体积计算结果。
3. 如权利要求1所述的方法，其特征在于，所述步骤2)的具体方法为： a)对多光谱遥感影像进行波段计算，得到NDVI分布图；b)采用表征冠层福射传输过程 的物理模型PROSAIL模拟获得用于拟合模型参数的叶面积指数；C)利用一定数量的LAI参 考值与NDVI数据绘制散点图，计算得到经验参数，即获得LAI遥感估算模型。
4. 如权利要求1所述的方法，其特征在于，所述步骤6)的具体方法为： a)使用高程阔值从经过植被掩膜后的图像提取建筑物，得到粗略的建筑物提取结果； b)采用中值滤波（Median Filter)消除边界"银齿现象"，得到边缘整齐的建筑物轮廓；C) 经由分区统计和几何计算，得到单个独立建筑物的高度、面积、体积信息。
【文档编号】G06F19/00GK104462741SQ201410444227
【公开日】2015年3月25日 申请日期:2014年9月3日 优先权日:2014年9月3日
【发明者】孟庆岩, 顾行发, 张佳晖, 吴俊 , 占玉林, 金颖 申请人:中国科学院遥感与数字地球研究所