一种基于城市绿化覆盖面积测算城市绿地率的方法与流程

本发明涉及一种绿地率的测算方法，具体是一种基于城市绿化覆盖面积测算城市绿地率的方法，属于城市园林绿化成果考评测试、城市园林绿化遥感调查、城市园林绿化现状普查等空间测绘领域以及城市总体规划、城市绿地系统规划的空间规划设计领域。

背景技术：

随着空间信息技术的发展，利用遥感影像测算城市绿化覆盖率的统计方法在中国逐渐推广。遥感测算法利用多光谱波段的光谱特性识别城市植被，并自动勾绘出被植被覆盖的区域，包括屋顶绿化和乔木树冠的垂直投影区域。该方法减少了大量的实地调研工作量，并适应绿化覆盖率动态变化的特点，实现了对城市绿化覆盖率的动态监测，成为业内普遍认同的城市绿化覆盖率测算技术。

然而，城市绿地率的测算方法却没能从遥感影像中获得如绿化覆盖率一般的提升。究其原因，是由两率的统计规则差异造成的。用于计算绿化覆盖率的绿化覆盖面积是指城市中乔木、灌木、草坪等所有植被的垂直投影面积，即在遥感影像中具有绿色光谱特征的图斑面积，其中，乔灌木是按其冠幅投影面积来计算的。而用于计算绿地率的绿地面积，则是统计落在体具城市用地上的指标，即乔灌木是按其生长的种植用地来计算的。例如，一株行道树，其绿化覆盖面积是行道树的树冠面积，其绿地面积则是树冠下种植池的面积。由于遥感影像自动识别技术只能对遥感影像表层的光谱特征进行识别，不能穿过树冠层识别绿化用地层，从而使绿地率的自动提取较为困难。因此，基于遥感影像计算城市绿地率依然需要专业人员用目视解译的方法，根据自身的专业知识，在影像图中人工勾绘城市绿地面积。例如，苏国振、吴明秋等在gis操作平台上，结合遥感影像，分别对浙江省嘉兴市和象山县的绿地率进行了测算。其基本的技术流程为：将城市大比例（不小于1：500）地形图或竣工测绘图导入gis平台作为绿地率调查的源数据，对照遥感影像逐区域编辑源数据，核实并勾绘每块城市绿地，对绿地边界不准确图斑进行外业调查，通过质量检查后生成绿地率成果。该方法需对城市中的所有绿地进行边界勾绘，需要投入大量的人力物力，且耗时较长，统计效率低。

技术实现要素：

本发明要解决的技术问题是提供一种基于城市绿化覆盖面积测算城市绿地率的方法，在确保准确性的前提下，大幅减少城市绿地率的统计工作量，节约调查统计时间，解决城市两率调查工作的时效性瓶颈问题。

为了解决所述技术问题，本发明采用的技术方案是：一种基于城市绿化覆盖面积测算城市绿地率的方法，其特征在于：包括以下步骤：s01）、基于城市遥感影像测算城市绿化覆盖面积，获取城市遥感影像，采用面向对象的分类法与人工补判相结合来提取城市绿化覆盖数据，统计城市建成区绿化覆盖面积；

s02）、在遥感影像上测算城市绿地率正相关修正值，正相关修正值包括行道树绿化覆盖面积和行道树绿带面积；

s03）、在遥感影像上测算城市绿地率负相关修正值，负相关修正值包括全部或者部分纳入绿地率统计的公园水体面积和公园中未被植被覆盖的大尺度铺装面积；

s04）、基于城市绿化覆盖面积测算城市绿地率，城市绿地率计算公式为：

rgl=[sgc-sstc+lst×wst+spw+spl/2+spf]/sub，其中rgl表示城市绿地率，

sgc表示城市绿化覆盖面积，sstc表示行道树绿化覆盖面积，lst表示行道树绿带长度，wst表示行道树绿带宽度，spw表示全部纳入绿地率统计的公园水体面积，spl表示部分含水体公园的陆地面积，spf表示公园中未被植被覆盖的大尺度铺装面积，sub表示城市建成区面积。

进一步的，采用面向对象的分类法与人工补判相结合来提取城市绿化覆盖数据的具体操作为：在ecognitiondevelopertrial9.0软件平台下对遥感影像进行多尺度分割，利用图像纹理信息和波段统计信息，构造绿化覆盖分类函数ndvi，对城市绿化覆盖信息进行提取，并输出矢量图形；在arcmap软件平台下，结合绿化图斑特征，进行人工分类补判，确保绿化覆盖信息分类精度达到90%以上，构建城市建成区绿化覆盖面源数据图层；统计城市建成区绿化覆盖面积sgc。

进一步的，测算行道树绿化覆盖面积和行道树绿带面积的具体操作为：目视解译城市遥感影像中0.8-3m宽的线形行道树图斑，直接在城市绿化覆盖面源数据图层上勾绘并裁剪，构建行道树绿化覆盖面积面源数据图层，测算行道树绿化覆盖面积修正值sstc；根据勾绘的行道树绿化覆盖图层，勾绘相对应的行道树绿带线源数据图层，测算行道树绿带长度修正值lst，根据实地取样测量数据，测算城市行道树绿带平均宽度修正值wst。

进一步的，测算城市绿地率负相关修正值的具体操作为：目视解译城市遥感影像中所有含水体公园，分别勾绘公园内的水体边界和陆地边界，构建含水体公园水体面源数据图层、含水体公园陆地面源数据图层，对比每个公园中水体面积spwi与陆地面积spli的关系，将不纳入绿地率统计的范围删除，测算spw与spl修正值；目视解译城市遥感影像中所有公园未被植被覆盖的大尺度铺装，构建公园大尺度铺装面源数据图层，测算公园中未被植被覆盖的大尺度铺装面积spf。

进一步的，若spwi≤spli/2，i代表含水公园的编号，i=1,2,3……，说明该公园的水体面积小于公园陆地面积的50%，此种情况下，水体全部计入绿地率，即删除该公园的spli数据，将spwi计入绿地率统计；若spwi＞spli/2，说明该公园的水体面积大于公园陆地面积的50%，此种情况下，水体按公园陆地面积的50%计入绿地率，删除该公园的spwi数据，取spli/2计入绿地率统计；从属性表中面积求和得到城市需全部纳入绿地率统计的公园水体面积和部分纳入绿地率统计的公园水体面积即部分含水体公园的陆地面积，测算修正值spw和spl。

进一步的，对遥感影像进行多尺度分割前，将遥感影像的默认波段修改为layer1—b、layer2—g、layer3—r、layer4—g，并进行标准差拉伸，突出遥感影像中的绿色植被信息的显示效果；进行多尺度分割时，将显示植被波段layer4的权重设置为2，其余波段的权重值为1，将分割尺度参数设置为200，将shape形状因子参数设置为0.1，平滑度参数设置为0.5，进行自下而上的多尺度分割；在featureview中创建植被指数函数ndvi提取植被，ndvi的特征定义公式为（[meannir]-[meanr]）/（[meannir]+[meanr]），在ndvi更新的窗口中对植被和非植被进行区分，经过对比调试，最终确定ndvi的阈值为0.29—0.72，在此阈值内进行植被提取。

进一步的，人工补判时，去除与绿色植被色彩相近的建筑屋顶或铺地面积，补充在阴影中未被识别的植被或被云层遮挡的植被。

本发明的有益效果：本发明在深入研究了城市绿地率与绿化覆盖率的统计意义与量化关系后，提出了一个利用遥感影像识别技术快速计算城市绿地率的统计方法，在城市绿化覆盖面积的遥感影像自动识别技术的基础上，只需在gis平台测算行道树绿化覆盖面积（sstc）、行道树绿带长度（lst）与宽度（wst）、全部纳入绿地率统计的公园水体面积（spw）、部分含水体公园的陆地面积（spl）、公园中未被植被覆盖的大尺度铺装面积（spf）等6个修正值，即可得出城市绿地率的数值。该方法极大的缩小了勾绘城市绿地的工作量，提高了城市绿地率的统计效率。为园林城市测评、城市绿化横向纵向比较等工作提供了高效、便捷的技术支持。

附图说明

图1为城市各类用地两率差值（城市绿化覆盖率-城市绿地率）分析图；

图2为道路行道树绿带示意图；

图3为遥感影像中行道树绿带两率示意图；

图4为遥感影像中城市绿化覆盖面源数据图层与行道树绿化覆盖面源数据图层示意图；

图5为行道树绿带线源数据图层示意图；

图6为spwi≤spli/2公园水体面源数据图层示意图；

图7为spwi＞spli/2公园水体面源数据图层示意图；

图8为公园中未被植被覆盖的大尺度铺装示意图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步的说明。

实施例1

本实施例公开一种基于城市绿化覆盖面积测算城市绿地率的方法，包括以下步骤：

s01）、基于城市高分辨率遥感影像测算城市绿化覆盖面积

（1）购买遥感数据

向卫星或航空遥感企业购买城市建成区遥感影像数据。首先应向企业提供城市建成区范围边界(sub)图，边界图应为西安80坐标系或2000国家大地坐标系的dwg格式文件。购买的遥感影像应满足以下要求：影像应覆盖城市的建成区范围；影像应包含全色和多光谱数据；融合后的影像空间分辨率应不低于1m；数据云量应小于20%，云、雾尽量不落在建成区范围内。

（2）测算城市绿化覆盖面积

采用面向对象的分类法与人工补判相结合来提取城市绿化覆盖数据。面向对象分类方法，是基于基元或者影像对象的分类方法。相对于传统分类方法它具有以下优势：模拟人脑解译方式，充分利用高分辨率遥感影像丰富的纹理和结构多种特征，提高信息分类精度；更适合处理空间尺度、空间分析等问题；面向对象方法可以在不同尺度空间提取特定主题的信息，提高了面向对象分类方法的分类精度。

具体可分为以下五步：

1、在ecognitiondevelopertrial9.0软件平台下导入img格式的遥感影像。

2、进入editimagelayermixing将系统默认的波段layer1—r、layer2—g、layer3—b改为layer1—b、layer2—g、layer3—r、layer4—g，并在equalizing菜单中选取standarddeviation（3.00）进行标准差拉伸，突出遥感影像中的绿色植被信息的显示效果。

3、在processtree窗口进入appendnew进行多尺度分割进程，在insertchild的子进程中插入多尺度分割算法“multiresolutionsegmentation”。在imagelayerweights中显示植被波段layer4的权重设置为2，其余波段的权重值为1，scaleparameter中将分割尺度参数设置为200。将shape形状因子参数设置为0.1，compactness平滑度参数设置为0.5，点击execute进行自下而上的多尺度分割。

4、在featureview中createnewarithmetic中创建植被指数函数ndvi提取植被，ndvi的特征定义公式为（[meannir]-[meanr]）/（[meannir]+[meanr]）。在ndvi更新的窗口中对植被和非植被进行区分。经过对比调试，最终确定植被的ndvi阈值为0.29—0.72，在此阈值内进行植被提取。

5、将ecognitiondevelopertrial9.0提取的植被矢量数据导入arcmap10.3中，结合绿化图斑特征，进行人工分类补判。去除与绿色植被色彩相近的建筑屋顶或铺地面积，补充在阴影中未被识别的植被或被云层遮挡的植被，确保绿化覆盖信息分类精度达到90%以上，构建城市建成区绿化覆盖面源数据图层（sgc），统计城市建成区绿化覆盖面积。

s02）、测算城市绿地率正相关修正值

（1）测算修正值sstc

通过对城市建成区内各类建设用地中绿化覆盖面积和绿地面积对城市两率及两率差值（绿化覆盖率-绿地率）影响的分析研究发现，城市两率差值的正相关主导因子是交通设施用地中的行道树绿带，如图1所示。因此，城市绿地率的修正应去除城市绿化覆盖面积中的行道树绿化覆盖面积，补增城市行道树绿带的绿地面积。

在arcmap10.3软件平台下，对城市行道树覆盖图斑进行目视解译，解译原则为：线形、连续、乔木，绿带长度不低于3m，宽度在0.8-3m之间，如图2、3所示。图3中被圈部分为行道树绿化覆盖面积及其统计范围，带点直线部分为行道树绿带绿地面积及其统计范围。在城市绿化覆盖图层勾选符合条件的行道树绿带后，执行相交（intersect）、剪裁(clip)命令，得到行道树绿化覆盖面源数据图层，如图4所示，测算修正值sstc。

（2）测算修正值lst和wst

城市行道树绿带的面积可通过公式sst=lst×wst求得，在arcmap10.3软件平台下，勾画与第一步sstc图层中行道树覆盖面相对应的行道树线要素，构建行道树绿带线源数据图层，如图5所示，从属性表中找到线的长度数据计算行道树绿带总长度，测算修正值lst；根据《城市道路绿化规划与设计规范》（cjj75-97）要求和实地观察测量数据，求城市行道树绿带平均宽度，测算修正值wst。

s03）、测算城市绿地率负相关修正值

（1）测算修正值spw和spl

如图1所示，城市两率差值的负相关主导因子是公园绿地中的水体与铺装。由于水体对公园绿地的景观塑造和休闲游憩发挥了重要作用，因此允许城市绿地率统计时纳入部分水体。作为绿地率统计的一种特殊情况，公园中水体与陆地的面积对比计算是不可省略的。

具体操作为：目视解译含水体公园，勾绘公园绿地内的水体，构建公园水体面源数据图层（spw）；勾绘含水体公园绿地内的陆地，构建含水体公园陆地面源数据图层（spl）；测算每个含水体公园中spwi(i代表含水公园的编号，i=1,2,3……)与spli的关系：若spwi≤spli/2，如图6所示，则该公园的绿地率面积统计需在绿化覆盖面积的基础上增加公园水体面积，此种情况的公园删除其spli数据，保留spwi数据；若spwi＞spli/2，如图7所示，则该公园的绿地率面积统计需在绿化覆盖面积的基础上增加部分公园水体面积，即公园陆地面积的1/2（spli/2），此种情况的公园删除其spwi数据，保留spli数据；从属性表中面积求和得到城市需全部纳入绿地率统计的公园水体面积和部分纳入绿地率统计的公园水体面积即部分含水体公园的陆地面积，测算修正值spw和spl。

（2）测算修正值spf

基于遥感影像解译的城市绿化覆盖面积不包含公园绿地中未被植物覆盖的大广场面积，而城市绿地率计算规则要求公园绿地中的铺装应全部纳入绿地率统计。因此，在公园绿地面积统计时需在绿化覆盖面积的基础上增加公园中未被植被覆盖的大尺度铺装面积。目视解译城市公园绿地中未被植被覆盖的大尺度铺装，单独勾绘公园铺装面源数据图层，测算修正值spf。

s04）、基于城市绿化覆盖面积测算城市绿地率

通过gis属性表统计sgc和sub基础数据，以及sstc、lst、wst、spw、spl、sgf等6个修正数据，代入城市绿地率修正公式rgl=[sgc-sstc+lst×wst+spw+spl/2+spf]/sub测算城市绿地率。

本实施例在深入研究了城市绿地率与绿化覆盖率的统计意义与量化关系后，提出了一个利用遥感影像识别技术快速计算城市绿地率的统计方法，在城市绿化覆盖面积的遥感影像自动识别技术的基础上，只需在gis平台测算行道树绿化覆盖面积（sstc）、行道树绿带长度（lst）与宽度（wst）、全部纳入绿地率统计的公园水体面积（spw）、部分含水体公园的陆地面积（spl）、公园中未被植被覆盖的大尺度铺装面积（spf）等6个修正值，即可得出城市绿地率的数值。该方法极大的缩小了勾绘城市绿地的工作量，提高了城市绿地率的统计效率。为园林城市测评、城市绿化横向纵向比较等工作提供了高效、便捷的技术支持。

以上描述的仅是本发明的基本原理和优选实施例，本领域技术人员根据本发明做出的改进和替换，属于本发明的保护范围。