干货分享 | 规划用地变化分析技术分享

1.基础数据准备

根据项目实际情况，拿到数据不尽相同。仅以常见的CAD数据为例进行说明，其它格式在此不多讨论。小鸟以XX项目2010年版和2017年版土地利用规划图为例，如下图所示：

由于各地出图要求差异，加大了机器识别的困难。不论拿到的CAD数据质量如何，都要人工对其进行梳理。这里没有捷径，肯定是要花费一些时间；但请你相信，这时间花的绝对不亏！

这里不赘述处理过程和操作手法，只谈最后结果的要求：

参照新版《城市用地分类与规划建设用地标准》要求，分如下九大类用地：非建设用地（E）、居住用地（R）、公共管理与公共服务用地（A）、商业服务业设施用地（B）、工业用地（M）、物流仓储用地（W）、道路与交通设施用地（S）、公用设施用地（U）、绿地与广场用地（G）。

●  用英文代码或数字作为图层名称，将各类要素梳理到对应图层

●  进入图层的要素必须是闭合二维多段线，不能是BH填充块

●  范围单设图层，同样是闭合二维多段线，图层名：规划范围

在此仅选取部分进行示意性操作，结果如下图所示：

目标基础文件的图层示意

目标基础文件分层示意（2017土地利用）

通过以上的人工处理，配合使用“Conversion Tools｜转出至地理数据库｜要素类至要素类”命令（仅需保留Layer图层），将CAD数据转换成shpfile文件。同时在ArcMap中适当可视化调整，成果如下图所示：

   2011年土地利用      

2017年土地利用

2.变化矩阵设计

[逻辑构思] - 假设有三个类别ABC相互做差值计算：A-B、A-C、B-A、B-C、C-A、C-B。如果差值均不相等，就可由差值推算是哪两个数相减。如果将ABC这三个数字定义成三类要素，是不是就得到了一种要素互变关系呢？

思路是不是清晰了，也很容易理解：纯粹的数字游戏。接下来就从最简单的三要素互变开始，满足要求数列：1,2,4；1,2,5；1,3,6；1,3,7；1,4,8 … …

四要素互变数列：1,2,4,8；1,2,5,9；1,3,6 … …

五要素互变数列：1,2,4,8,13；1,2,5,9,14

依此类推... ...操作示例数据涉及七要素，选用如下数列：

{ 1,2,4,8,13,21,31 }

变换矩阵如下图。

数据和矩阵有了，开始用地变化比对操作吧！

 3.用地变化比对

使用“Conversion Tools｜转为栅格｜面转栅格”命令将规划用地的Shpfile文件转为栅格，参数设置如下：值字段“Layer”，像元大小“1”，其余默认，输出并保存栅格数据集命名“new” 和“old”。

使用“Spatial Analyst Tools｜重分类｜重分类”命令对转出栅格进行重分类，参数设置如下：重分类字段“LAYER”，重分类数值“A=1，B=2，E=4，G=8，R=13，U=21，S=31”，其余默认，输出并保存栅格数据集命名“new_re”和“old_re”。

“重分类”工具参数

使用“Spatial Analyst Tools｜地图代数｜栅格计算器”命令对两版栅格文件做减法运算，地图代数表达式“new_re - old_re”。计算结果如下图所示：

地图代数计算成果栅格

依据VALUE值与选用七要素互变矩阵逐个对应的原则将COUNT值填入，得到面积统计矩阵。（“COUNT”代表“面积”，像元面积取决于栅格大小）前后两版用地性质未有变化的用地面积统计：1016421平方米，其中变与未变分布情况如下图：

用地性质变化的面积统计如下图所示：

以“公共管理与公共服务用地”（A类）变化为例，阐释此表含义：

●  A类转向“商业服务业设施用地”，8699平方米

●  A类转向“绿地与广场用地”，1936平方米

●  A类转向“居住用地”，31677平方米

●  A类转向“道路与交通设施用地”，12859平方米

转载自月月小鸟，经作者授权转载

文章仅代表作者观点，与本公众号无关，版权归原作者所有

原文标题：原创丨规划用地变化分析技术分享

图文编辑：来 莱

责任编辑：张英浩

   审编：鲁嘉颐

终审： 颜子明  黄宗财  梁龙武

猜你喜欢

2. 地理学求索 | “十五周年纪念特刊”

扫描二维码，关注我们

都看到这里了，点个“在看”再走呗~👇