【ArcGIS Pro二次开发】(29)：村庄规划生成空间功能结构调整表

根据现在村规成果要求，【空间功能结构调整表】是必需的。

以福建省为例，它长这样：

下面就来实现从现状用地和规划用地导出这样的Excel表格。

一、要实现的功能

如上图所示，点击【汇总村庄空间功能结构调整表】工具，选择【现状用地】和【规划用地】图层，选择要保存输出的Excel文件路径，再点击【执行】即可。

生成结果如下：

PS:注意输入的【现状用地】和【规划用地】的规范要求

一、【现状用地】

输入的【现状用地】图斑需要注意要有2个字段：【JQDLBM】和【CZCSXM】。

【JQDLBM】是村规入库的标准代码，即基期地类编码。

【CZCSXM】是三调自带的字段，用来表示三调用地的城镇村属性。201、201A是城市建设用地，202、202A是建制镇建设用地，203、203A是村庄建设用地。这里提取【201、201A、202、202A】作为城镇用地。

二、【规划用地】

输入的【规划用地】图斑需要注意要有2个字段：【GHDLBM】和【SSBJLX】。

【GHDLBM】和【SSBJLX】都是村规入库的标准代码。

【GHDLBM】即规划地类编码。

【SSBJLX】用来表示规划用地的城镇村属性，字段值为【Z】则表示该地块为城镇用地。这里为了照顾输入时的误差，字段值为小写的【z】也可以被认得。

二、实现流程

工具的实现流程和我的上一篇文章很像：

【ArcGIS Pro二次开发】(28)：用地图斑导出用地用海汇总表

流程长但不复杂，更重要的是Excel模板和输入【现状用地、规划用地】的准备。

1、准备Excel模板

一、【空间功能结构调整表】模板

首先是要输出表的模板，需要输入的单元格包括面积、比重、和面积增减。这里的大部分面积值由用地汇总而来，部分合计值、比重以及面积增减值就在Excel里用公式计算。另外，还需要一个参照列，用来属性映射。

二、【用地用海代码_村庄功能】属性映射表

这个表格是作为用地用海代码向村庄功能映射的参照表，将用地用海代码转换成村庄功能代码。

2、主要代码

实现功能的主要流程如下：

            // 获取默认数据库string init_gdb = Project.Current.DefaultGeodatabasePath;// 获取默认文件夹位置string init_path = Project.Current.HomeFolderPath;// 获取参数string yd_xz = combox_fc_xz.Text;           // 现状用地string yd_gh = combox_fc_gh.Text;         // 规划用地string excel_GN = init_path + @"\用地用海代码_村庄功能.xlsx";    // 村庄功能映射Excelstring excel_TZB = textExcelPath.Text;      // 输出空间功能结构调整表Excelstring output_table_xz = init_gdb + @"\output_table_xz";string output_table_gh = init_gdb + @"\output_table_gh";string statistics_fields = "Shape_Area SUM";string field_xz = "JQDLBM";string field_gh = "GHDLBM";string field_GN = "村庄功能";// 复制模板Excel和输出结果ExcelToolManager.CopyResourceFile(@"CCTool.Excel.用地用海代码_村庄功能.xlsx", excel_GN);ToolManager.CopyResourceFile(@"CCTool.Excel.【模板】空间功能结构调整表.xlsx", excel_TZB);// 现状、规划用地功能映射GN_Mapper(yd_xz, field_xz, field_GN, excel_GN);GN_Mapper(yd_gh, field_gh, field_GN, excel_GN);// 统计面积ToolManager.MultiStatistics(yd_xz, output_table_xz, statistics_fields, new List<string>() { field_GN }, "总计", 1, true);ToolManager.MultiStatistics(yd_gh, output_table_gh, statistics_fields, new List<string>() { field_GN }, "总计", 1, true);// 删除中间字段Arcpy.DeleteField(yd_xz, field_GN);Arcpy.DeleteField(yd_gh, field_GN);// 属性映射至ExcelToolManager.ExcelAttributeMapper(excel_TZB + @"\sheet1$", 13, 6, "output_table_xz", @"分组字段", "SUM_Shape_Area", 5);ToolManager.ExcelAttributeMapper(excel_TZB + @"\sheet1$", 13, 8, "output_table_gh", @"分组字段", "SUM_Shape_Area", 5);// 删除0值行ToolManager.Delete0Row(excel_TZB + @"\sheet1$", new List<int>() {6, 8}, 5);// 删除参照列ToolManager.DeleteColumns(excel_TZB + @"\sheet1$", new List<int>() { 13 });// 删除中间数据Arcpy.Delect(new List<string>() { output_table_xz, output_table_gh });

其中【GN_Mapper】方法用于计算城村庄功能字段：

        // 村庄功能映射public async void GN_Mapper(string yd, string field_bm, string field_GN, string excel_GN){// 汇总统计现状用地Arcpy.AddField(yd, field_GN, "TEXT");     // 添加一个参照字段// 属性映射ToolManager.AttributeMapper(yd, field_bm, field_GN, excel_GN + @"\sheet1$");// 计算城镇用地var map = MapView.Active.Map;FeatureLayer init_layer = map.FindLayers(yd)[0] as FeatureLayer;await QueuedTask.Run(() =>{using (ArcGIS.Core.Data.Table table = init_layer.GetTable()){using (RowCursor rowCursor = table.Search(null, false)){TableDefinition tableDefinition = table.GetDefinition();while (rowCursor.MoveNext()){using (Row row = rowCursor.Current){if (field_bm == "JQDLBM"){// 获取valuevar va = row["CZCSXM"].ToString();// 赋值if (va == "202" || va == "202A" || va == "201" || va == "201A"){row[field_GN] = "城镇用地";}}else if (field_bm == "GHDLBM"){// 获取valuevar va = row["SSBJLX"].ToString();// 赋值if (va == "Z" || va == "z"){row[field_GN] = "城镇用地";}}row.Store();}}}}});}

流程里涉及到的其它功能包括：复制Excel、属性映射、统计面积、Excel属性映射、删除0值行等，可以在以往的文章里查到，具体可以看上一篇文章，或在文章末尾提供的工程文件里查看。实在也没什么好讲的。

PS:这次提供的工程文件和别的工具合在一起了。代码里的功能嵌套越来越多，分离开来反而很麻烦，就将整个工具包一块分享了。

三、工程文件分享

最后，放上工程文件的链接：

StatisticVillagehttps://pan.baidu.com/s/1nsjUHqRqfogHrjis5pXDXw?pwd=dwsv

PS：可以直接点击...bin\Debug\net6.0-windows\下的.esriAddinX文件直接安装。