一种用于评价城市绿化三维布局的绿度空间配置曲线构建方法
【专利摘要】本发明公开了一种用于评价城市绿化三维布局的绿度空间配置曲线构建方法，包括如下步骤：步骤1)提取单株立木识别与冠层结构信息；步骤2)基于多光谱影像反演叶面积指数，构建冠层空间结构-生态因子回归模型；步骤3)根据植被蒸腾效益扩散过程计算模型，分别计算间隔0.5m的连续高度层面上的绿辐射受益面积之和，绘制三维绿度空间曲线，并基于曲线获得绿辐射峰、绿辐射峰高度、绿辐射范围以及绿辐射集中受益范围等多项评价参数；步骤4)计算并绘制间隔0.5m的连续高度层面上的三维建筑物配置曲线；步骤5)将三维绿度空间曲线与反映各个高度上建筑面积之和的三维建筑物曲线绘制于同一数据框内，获得绿度空间配置曲线，服务于城市绿度空间评价。
【专利说明】一种用于评价城市绿化三维布局的绿度空间配置曲线构建 方法
【技术领域】
[0001] 本发明涉及一种用于评价城市绿化三维布局的绿度空间配置曲线构建方法，本发 明能够从城市绿度空间的三维垂直布局出发，量化建筑物与其邻接绿度空间的配置关系， 准确度量居民与城市绿地的接触概率。
【背景技术】
[0002] 我国城市化我国正处于快速城镇化进程中，随着城镇化建设的快速推进，城市环 境问题日益突出。快速城市化导致了空气污染蔓延、城市热岛效应加剧、水环境恶化、生物 多样性丧失等一系列环境问题，同时还引发城市核心地区土地资源紧张、功能区结构矛盾 等社会问题，对维持区域生态平衡和城市可持续发展造成极大挑战。伴随着城市化进程的 加快，加之城市生态环境问题的日益突出，人们已越来越认识到城市绿地对城市生态环境 和城市人居环境的重要价值，城市居民对城市生活环境质量要求也越来越高。
[0003] 为准确描述城市绿地空间分布以及城市居民接触城市绿地概率上的差异，一些学 者尝试建立城市绿地与城市居民间的关系模型。Ainura Nazarkulova参考地理建模空间分 析方法，利用格网法对绿地类型加权生成不同植被指数，并选取城市主要居民点做缓冲区 分析，分别计算环绕其周围间隔200m的所有环状缓冲区覆盖的绿地面积大小，从而得到指 定位置可以享受到的绿色空间总量。李小江基于LiDAR数据和多光谱数据建立了每个单体 建筑物与其周围绿度空间的接触模型，构建一系列城市绿度空间指数用于城市绿度空间研 究。
[0004] 但传统度量方法对于绿色空间生态效益的评估往往限于定性的描述，并不能准确 的描述绿度空间的生态效益，对于城市绿地在改善城市环境中的作用也只停留在感性的认 识层面，对于城市绿地的生态效应更无宏观层面上的量化的评价标准。
[0005] 鉴于以上，本发明提出三维绿度空间曲线、绿度空间与建筑物三维配置曲线评价 模型。通过城市建筑物与其邻接绿度空间的配置关系来评价人与绿地接触的概率，以及 接触概率的空间分布，以期为环境规划、城市绿化、环境保护和城市精细化管理提供参考依 据。
【发明内容】
[0006] 本发明提出一种用于评价城市绿化三维布局的绿度空间配置曲线构建方法，用于 客观定量评价建筑物与其邻接绿度空间的配置关系，以及居民与绿地的接触概率。
[0007] 本发明的目的通过以下技术步骤实现：
[0008] 步骤1)基于机载激光雷达数据与高分辨率航空影像，进行单株立木识别与冠层 结构信息提取，其中冠层结构信息主要包括树顶坐标、树高、树冠投影边界、冠层最大横截 面积、冠径、冠高、树冠体积等，进而绘制城市范围内的植被三维分布信息；
[0009] 步骤2)基于多光谱影像反演叶面积指数，并结合步骤1)所得树冠体积计算结果， 分析二者间相关系数，构建冠层空间结构-生态因子回归模型；
[0010] 步骤3)根据植被蒸腾效益扩散过程计算模型，获得城市绿辐射受益量分布图，借 助生态效益源构建方法，分别计算间隔0. 5m的连续高度层面上的绿辐射受益面积之和，经 由正态分布曲线拟合，绘制三维绿度空间曲线，并基于曲线获得绿辐射峰、绿辐射峰高度、 绿辐射范围以及绿辐射集中受益范围等多项评价参数；
[0011] 步骤4)使用步骤3)生成的建筑物平均高度、投影面积和体积等参数，分别计算并 绘制间隔〇. 5m的连续高度层面上的三维建筑物配置曲线；
[0012] 步骤5)将三维绿度空间曲线与反映各个高度上建筑面积之和的三维建筑物曲线 绘制于同一数据框内，获得三维配置曲线，通过分析该曲线的形态走势、线下面积与交点情 况，衍生出一系列评价统计参数（绿辐射匮乏层、绿辐射饱和层、绿辐射平衡点、绿度空间 有效服务比），服务于城市绿度空间评价。
[0013] 进一步所述步骤1)的具体方法为：
[0014] a)采用面向对象分类法，生成植被二值掩膜图像；b)将二维植被分布信息与数字 高度模型做叠加，构建树冠高度模型；c)在树冠高度模型基础上，开展树顶检测、树冠投影 边界识别以及树冠体积计算等信息提取步骤。
[0015] 进一步所述步骤2)的具体方法为：
[0016] a)首先对PR0SAIL模型参数进行敏感性分析；b)通过敏感性分析结果以及植被参 数实地测量确定模型的可变参数和固定参数，进而构建研究区植被叶面积反演的查找表； c)最后使用成本函数最小化方法搜索查找表进行植被叶面积指数反演；d)对叶面积指数 反演结果与树冠体积计算结果进行回归分析，获得二者的回归方程式。
[0017] 进一步所述步骤3)的具体方法为：
[0018] a)基于步骤2)所得冠层空间结构-生态因子回归模型，构建植被生态效益源；b) 根据植被蒸腾作用的扩散与衰减规律，估算植被水平蒸散梯度；c)使用步骤1)获得的数字 高度模型与植被掩膜图像，进行建筑物三维特征参数提取；d)基于上述步骤所得参数，构 建模拟植被蒸散效应的计算模型。
【专利附图】
【附图说明】
[0019] 图1为三维绿度空间曲线；
[0020] 图2为绿度空间与建筑物三维配置曲线（住宅区）；
[0021] 图3为绿度空间与建筑物三维配置曲线（商业区）；
【具体实施方式】
[0022] 下面结合附图对本发明"一种用于评价城市绿化三维布局的绿度空间配置曲线构 建方法"作进一步说明。
[0023](一）植被三维分布信息提取
[0024] 进行植被三维结构信息提取所需的数据源均为遥感数据，且全过程均可以自动化 实现。具体计算步骤如下：
[0025] 1)在数字地面模型的基础上滤除非植被信息，获得植被二值掩膜图像，将其叠加 至数字高度模型，获得树冠高度模型；
[0026] 2)基于活动窗口的局部最大值搜索法，从树冠高度模型中提取研究区范围内单株 立木的树顶位置及其高度；
[0027] 3)为获得诸如树高、冠层高度、树冠体积等结构参数，需要在树顶检测结果基础之 上进行树冠投影边界识别，此处采用伞骨算法识别树冠投影边界；
[0028] 4)最后，借助树冠最大横截面积与冠高比信息，估算椭球形树冠体积，获得单体树 冠体积分布图。
[0029](二）冠层空间结构-生态因子回归模型
[0030] 首先采用冠层辐射传输模型SAIL的耦合模型PR0SAIL反演叶面积指数，通过分析 叶面积指数与树冠体积的相关性，构建冠层空间结构-生态因子回归模型。具体计算过程 如下：
[0031] 1)通过对模型生理生化参数的敏感性进行分析，从而确定对冠层光谱反射率影响 较为敏感的参数；
[0032] 2)影响冠层光谱反射率的变动参数包括叶面积指数、平均叶倾角、叶绿素含量、胡 萝卜素素含量、等效水厚度和叶片干物质，将各个参数的变动值和固定值代入到PR0SAIL 模型中来模拟冠层光谱反射率，建立冠层叶面积指数查找表；
[0033] 3)根据大气校正后的多光谱影像反射率，在查找表中获取对应的冠层反射率，并 以代价函数最小时对应的叶面积指数作为该像元模拟的叶面积指数；
[0034] 4)利用叶面积反演结果绘制树冠体积与叶面积指数散点图，并对其尝试多种线性 和非线性回归分析，最终确定适合的冠层空间结构-生态因子回归模型。
[0035] (三）三维绿度空间曲线
[0036] 三维绿度空间曲线是由绿辐射受益量计算结果衍生出的一条用于客观描述城市 绿化空间分布的曲线。其中，绿辐射受益量的计算过程如下：
[0037] 1)将单株立木的树冠体积值与LAI反演量相结合，利用微积分定理，构建植被生 态效益源，用于量化冠层叶面积垂直分布数量；
[0038] 2)城市绿度空间由生态效益源和绿辐射空间两部分构成，其中绿辐射空间是指生 态效益源向外辐射释放的、对城市环境具有积极影响的环境因子在三维地域空间网格内的 含量及分布规模，多采用仪器测定法获得植被蒸腾作用的扩散与衰减特征，进而估算得出 植被水平蒸散梯度；
[0039] 3)城市建筑物具有一定高程，因而很容易借助高程阈值从滤除植被信息的掩膜图 像中提取建筑物，后经相关几何计算，获得单体建筑物的位置、体积等三维信息；
[0040] 4)将上述步骤所得参数代入植被蒸腾效益扩散过程计算模型，生成一系列以 〇. 5m为间隔的不同高度层上的绿辐射受益总量。
[0041] 将连续高度层上的绿辐射受益总量绘制在二维坐标轴下，即获得三维绿度空间曲 线，如图1所示，其中X轴表示高度，Y轴表示面积。不同区域的绿化格局会生成不同的曲 线形态，鉴于该曲线的对称性形如正态分布曲线，故在此使用正态分布函数对其进行拟合， 得到拟合函数如下：
【权利要求】
1. 一种用于评价城市绿化H维布局的绿度空间配置曲线构建方法，该曲线的构建过程 包括如下步骤： 步骤1)基于机载激光雷达数据与高分辨率航空影像，进行单株立木识别与冠层结构 信息提取，其中冠层结构信息主要包括树顶坐标、树高、树冠投影边界、冠层最大横截面积、 冠径、冠高、树冠体积等，进而绘制城市范围内的植被H维分布信息； 步骤2)基于多光谱影像反演叶面积指数，并结合步骤1)所得树冠体积计算结果，分析 二者间相关系数，构建冠层空间结构-生态因子回归模型； 步骤3)根据植被蒸腾效益扩散过程计算模型，获得城市绿福射受益量分布图，借助生 态效益源构建方法，分别计算间隔0. 5m的连续高度层面上的绿福射受益面积之和，经由正 态分布曲线拟合，绘制H维绿度空间曲线，并基于曲线获得绿福射峰、绿福射峰高度、绿福 射范围W及绿福射集中受益范围等多项评价参数； 步骤4)使用步骤3)生成的建筑物平均高度、投影面积和体积等参数，分别计算并绘制 间隔0. 5m的连续高度层面上的H维建筑物配置曲线； 步骤5)将H维绿度空间曲线与反映各个高度上建筑面积之和的H维建筑物曲线绘制 于同一数据框内，获得H维配置曲线，通过分析该曲线的形态走势、线下面积与交点情况， 衍生出一系列评价统计参数（绿福射匿乏层、绿福射饱和层、绿福射平衡点、绿度空间有效 服务比），服务于城市绿度空间评价。
2. 如权利要求1所述的方法，其特征在于，所述步骤1)的具体方法为： a)采用面向对象分类法，生成植被二值掩膜图像；b)将二维植被分布信息与数字高度 模型做叠加，构建树冠高度模型；C)在树冠高度模型基础上，开展树顶检测、树冠投影边界 识别W及树冠体积计算等信息提取步骤。
3. 如权利要求1所述的方法，其特征在于，所述步骤2)的具体方法为： a)首先对PROSAIL模型参数进行敏感性分析；b)通过敏感性分析结果W及植被参数实 地测量确定模型的可变参数和固定参数，进而构建研究区植被叶面积反演的查找表；C)最 后使用成本函数最小化方法搜索查找表进行植被叶面积指数反演；d)对叶面积指数反演 结果与树冠体积计算结果进行回归分析，获得二者的回归方程式。
4. 如权利要求1所述的方法，其特征在于，所述步骤3)的具体方法为： a)基于步骤2)所得冠层空间结构-生态因子回归模型，构建植被生态效益源；b)根据 植被蒸腾作用的扩散与衰减规律，估算植被水平蒸散梯度；C)使用步骤1)获得的数字高度 模型与植被掩膜图像，进行建筑物H维特征参数提取；d)基于上述步骤所得参数，构建模 拟植被蒸散效应的计算模型。
【文档编号】G06T17/30GK104463971SQ201410491578
【公开日】2015年3月25日 申请日期:2014年9月24日 优先权日:2014年9月24日
【发明者】张佳晖, 孟庆岩, 金颖, 杨健, 刘苗, 王春梅 申请人:中国科学院遥感与数字地球研究所