利用国土三调与林业数据寻找造林空间的方法

詹福雷

(沈阳市勘察测绘研究院有限公司，辽宁 沈阳 110000)

1 引 言

为了践行“绿水青山就是金山银山”的发展理念，推动大规模造林绿化工作，扭转城市森林覆盖率低、生态环境脆弱、土地沙化和水土流失状况，增进城市人民生态福祉，满足生态文明建设和美丽城市的迫切需求，在城市范围内寻找可造林空间势在必行[1～7]。造林空间选取必须有可靠的数据支撑，既要考虑土地利用现状，又要考虑土地规划用途等。第三次全国国土调查以国土二调成果为基础，全面细化完善全国土地利用数据，掌握翔实准确的全国国土利用现状，实现成果信息化管理与共享，满足生态文明建设、空间规划编制、供给侧改革、宏观调控、自然资源管理体制改革和统一确权登记、国土空间用途管制、国土空间生态修复、空间治理能力现代化和国土空间规划体系建设等各项工作[8～10]。三调数据完全满足造林空间的选取需求。

2 探寻造林空间

国土部门与林业部门统计的林地面积总是不一致，这种不一致导致的两部门之间数据重叠部分正是造林空间存在的区域。依据现有数据情况我们将造林空间分为：“应绿未绿”“无立木、宜林地”“非林地”“特定区域数据造林空间”这四种类型。

2.1 应绿未绿

2017年林业变更数据库中为林地，按照林业部门管理来讲这样的地块应该已经是林地，通过三调现状数据进行分析找出现状不是林地的地块，这样的造林空间我们称之为应绿未绿。

2.2 无立木、宜林地

按照林业部门对林业数据管理来说这样地块按照林地管理但现状应该没有林木，但是由于各种原因这部分地块有一些已经有林木了，故而也需要利用三调现状地类找出这部分没有林木的部分，这样的造林空间我们称之为无立木、宜林地。

2.3 非林地

这部分地块按照林业部门管理来说不按照林地管理，但是由于各种原因这部分数据县区、乡镇、村和林户有林地承包合同，这部分也可作为造林空间，这部分数据依据三调现状数据库结合林地承包合同将签订合同的地块现状没有林木的提取出来作为造林空间，这样的造林空间我们称之为非林地。

2.4 特定区域数据造林空间

特定区域指的是具有政策导向的地块如退耕还林，生态封育等，这种区域范围内有鲜明的政策导向，但是现状能否造林还要以三调现状调查为依据，这样的造林空间我们称之为特定区域数据造林空间。

2.5 总体技术流程(图1)

3 造林空间提取方法

根据本文前部分论述，我们定义了造林所需要的造林空间，接下来我们利用国土三调数据与林业数据来提取本市区域内相关的造林空间，四种造林空间的提取方法如下：

图1 总体技术流程图

3.1 应绿未绿

筛选2107年林业变更数据库中地类如下：

"C_DILEI"="111(针叶林)"OR"C_DILEI"= "112(阔叶林)"OR"C_DILEI"="113(针阔混交林)"OR"C_DILEI"="114(乔木经济林)"OR"C_DILEI"="121(疏林地)"OR"C_DILEI"="131(特规灌木林地)"OR"C_DILEI"="132(其他灌木林地)"OR"C_DILEI"="141(造林未成林地)"OR"C_DILEI"="142(封育未成林地)"OR"C_DILEI"="151(苗圃地)"，得到林业部门按照林地管理的范围，导出这部分数据。利用ArcGIS中Intersect_analysis工具将导出的图层与三调阶段性成果进行空间相交处理，之后筛选 "DLMC"="旱地"OR"DLMC"="水田"OR"DLMC"="其他草地"等土地利用现状不是林地且满足造林地理条件的相关地类，得到应绿未绿造林空间。

3.2 无立木、宜林地

筛选2107年林业变更数据库中地类如下：

"C_DILEI"="161(采伐迹地)"OR"C_DILEI"= "163(其他无立木林地)"OR"C_DILEI"="164(临时占用林地)"OR"C_DILEI"="165(毁林开垦)"OR"C_DILEI"="171(宜林荒山荒地)"OR"C_DILEI"="172(宜林沙荒地)"OR"C_DILEI"="173(其他宜林地)"，得到林业部门按照林地管理但是土地利用现状应该没有林木的范围，导出这部分数据。利用ArcGIS中Intersect_analysis工具将导出的图层与三调阶段性成果进行空间相交处理，之后筛选"DLMC"="乔木林地"OR"DLMC"="其他林地"，将这部土地利用现状已经有林木的部分删除，再筛选"DLMC"="商业服务业设施用地"OR"DLMC"="工业用地"等建设用地及地理条件不能造林的部分将其删除，得到无立木、宜林地造林空间。

3.3 非林地

收集县区、乡镇、村与林户签订的承包合同，通过图解或解析方法得到承包合同的矢量图形范围。利用ArcGIS中Intersect_analysis工具将导出的图层与三调阶段性成果进行空间相交处理，之后筛选"DLMC"="乔木林地"OR"DLMC"="其他林地"，将这部土地利用现状已经有林木的部分删除，再筛选"DLMC"="商业服务业设施用地"OR"DLMC"="工业用地"等建设用地及地理条件不能造林的部分将其删除，得到非林地造林空间。

3.4 特定区域数据造林空间

收集具有政策导向的地块如退耕还林，生态封育等数据，通过图解或解析方法得到矢量图形范围。利用ArcGIS中Intersect_analysis工具将导出的图层与三调阶段性成果进行空间相交处理，之后筛选"DLMC"="乔木林地"OR"DLMC"="其他林地"，将这部土地利用现状已经有林木的部分删除，再筛选"DLMC"= "商业服务业设施用地"OR"DLMC"="工业用地"等建设用地及地理条件不能造林的部分将其删除，得到特定区域数据造林空间。

4 造林空间数据检查

土地利用现状检查：对上述选出的地块利用最新高分辨率遥感影像核查造林空间现状是否正确，对于现状不正确的调出造林空间。

最小图斑面积检查：对上述选出的地块面积小于1亩的图形进行检查并且删除。

5 成果输出及调查底图制作

依据前文论述，选取一个乡镇的数据作为试验区。利用国土三调数据、林业数据、乡镇提供的林地承包合同数据、河流封育范围数据利用本文论述的方法快速的提取该乡镇的造林空间。以行政村为单位汇总了该乡镇造林空间统计表格如表1所示：

***乡镇造林空间统计表(单位/公顷) 表1

续表1

将上述成果分为应绿未绿，无立木、宜林地，非林地和特定区域数据四个图层，分别采用不同颜色表示，如图2所示。

图2 ***乡镇造林空间选取工作底图

利用ArcGIS输出kmz文件，提供给乡镇林业站管理人员及外业核查人员，进行实地核查形成内外业核查闭合流程。

6 结 语

本文从国土三调与林业数据出发，分析各部门对国土空间用途管制方法不同，利用ArcGIS分析数据的优势快速查找造林空间，通过一个乡镇的数据实例给出寻找造林空间的解决方案，很好地解决了造林空间寻找困难的问题，东北地区造林时节有限(4月初～5月1日前)，通过本文论述可以实现快速寻找造林目标区域满足造林时节需求[11]。