干货分享 | 规划用地变化分析技术分享
The following article is from 月月小鸟 Author 月月小鸟
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｜地图代数｜栅格计算器”命令对两版栅格文件做减法运算，地图代数表达式“new_re - old_re”。计算结果如下图所示：
地图代数计算成果栅格
依据VALUE值与选用七要素互变矩阵逐个对应的原则将COUNT值填入，得到面积统计矩阵。（“COUNT”代表“面积”，像元面积取决于栅格大小）前后两版用地性质未有变化的用地面积统计：1016421平方米，其中变与未变分布情况如下图：用地性质变化的面积统计如下图所示：
以“公共管理与公共服务用地”（A类）变化为例，阐释此表含义：
● A类转向“商业服务业设施用地”，8699平方米
● A类转向“绿地与广场用地”，1936平方米
● A类转向“居住用地”，31677平方米
● A类转向“道路与交通设施用地”，12859平方米
转载自月月小鸟，经作者授权转载
文章仅代表作者观点，与本公众号无关，版权归原作者所有
原文标题：原创丨规划用地变化分析技术分享
图文编辑：来 莱
责任编辑：张英浩
审编：鲁嘉颐
终审： 颜子明 黄宗财 梁龙武
猜你喜欢
扫描二维码，关注我们
都看到这里了，点个“在看”再走呗~👇