之前做了一个三调的土地利用现状分类统计表，是按照自己的习惯做的，发现不符合土规专业的要求，于是翻开三调规程，按规范重新制作。

一、要实现的功能

如上图所示，在【规划应用】组—【三调】面板下，点击【1、土地利用现状分类面积汇总表】工具。

在弹出的工具框中，分别输入参数：

1、输入三调图层。

2，3、参数2和3是可选参数。用于需要按地块进行分区统计的情况。参数3为参数2中用于标记地块名称的字段。

4、汇总面积。工具预设了3种，【shape_area、TBMJ、TBDLMJ】。【shape_area】即投影面积，【TBMJ】即椭球面积，【TBDLMJ】为椭球面积扣去扣除面积，即净面积。

5、面积单位。工具预设了【平方米、公顷】2种。

6、输出Excel统计表的位置。

生成结果如下：

二、实现流程

因为全部代码比较长，这里只贴出主流程代码，其中部分汇总数据以及对Excel表的操作，之前的文章都有说到，可以回头翻阅一下。

// 复制嵌入资源中的Excel文件
ToolManager.CopyResourceFile(@"CCTool.Data.Excel.三调规程.土地利用现状分类面积汇总表.xlsx", excel_path);
ToolManager.CopyResourceFile(@"CCTool.Data.Excel.三调用地自转换.xlsx", folder_path + @"\三调用地自转换.xlsx");// 放到新建的数据库中
Arcpy.CreateFileGDB(folder_path, "新数据库");
string new_gdb = $@"{folder_path}\新数据库.gdb";if (fc_area_path == "")  // 如果没有分地块统计
{Arcpy.CopyFeatures(fc_path, $@"{new_gdb}\{fc_path}");
}else   // 如果分地块统计
{// 按属性分割Arcpy.SplitByAttributes(fc_area_path, new_gdb, name_field);
}// 收集分割后的地块
List list_area = ToolManager.GetFeatureClassAndTablePath(new_gdb);
// 新建的gdb列表
List gdb_list = new List();foreach (var area in list_area)
{// 裁剪新用地string area_name = area[(area.LastIndexOf(@"\") + 1)..];string new_fc = $@"{new_gdb}\{area_name}_裁剪结果";Arcpy.Clip(fc_path, area, new_fc);// 如果是计算图斑面积，就重新计算if (area_field == "TBMJ"){Arcpy.CalculateField(new_fc, "TBMJ", "!shape.geodesicarea!");}else if (area_field == "TBDLMJ"){Arcpy.CalculateField(new_fc, "TBDLMJ", "!shape.geodesicarea! * (1-!KCXS!)");}// 添加一个类一、二级字段Arcpy.AddField(new_fc, "mc_1", "TEXT");Arcpy.AddField(new_fc, "mc_2", "TEXT");// 用地编码转换ToolManager.AttributeMapper(new_fc, "DLMC", "mc_1", folder_path + @"\三调用地自转换.xlsx\一级$");ToolManager.AttributeMapper(new_fc, "DLMC", "mc_2", folder_path + @"\三调用地自转换.xlsx\二级$");// 将新用地按【坐落地代码名称】分割放入新的数据库，并收集，方便下一步处理Arcpy.CreateFileGDB(folder_path, $"{area_name}_裁剪结果");gdb_list.Add(@$"{folder_path}\{area_name}_裁剪结果.gdb");// 再分割Arcpy.SplitByAttributes(new_fc, $@"{folder_path}\{area_name}_裁剪结果.gdb", "ZLDWDM;ZLDWMC");// 收集分割后的地块List list_area_detail = ToolManager.GetFeatureClassAndTablePath($@"{folder_path}\{area_name}_裁剪结果.gdb");// 新建sheetToolManager.ExcelCopySheet(excel_path, "sheet1", area_name);string new_sheet_path = excel_path + @$"\{area_name}$";// 处理每个行政单位int start_row = 7;List row_list = new List();foreach (var area_detail in list_area_detail){// 复制行ToolManager.CopyRow(new_sheet_path, "A6:BS6", $"A{start_row}");// 记录写入行数的列表row_list.Add(start_row);// 要处理的要素名称string area_name_detail = area_detail[(area_detail.LastIndexOf(@"\") + 1)..];// 获取行政代码和行政名称string ad_name = area_name_detail[(area_name_detail.LastIndexOf(@"_") + 1)..];string ad_code = area_name_detail[1..area_name_detail.LastIndexOf(@"_")];// 汇总大、中类List list_bm = new List() { "mc_1", "mc_2" };// 单位转算【1：公顷，4：平方米】int unit_index = 0;if (unit == "平方米") { unit_index = 4; }else if (unit == "公顷") { unit_index = 1; }ToolManager.MultiStatistics(area_detail, gdb_path + @"\statistic_sd", area_field + " SUM", list_bm, "国土调查总面积", unit_index);// 将映射属性表中获取字典DictionaryDictionary dict = ToolManager.GetDictFromTable(gdb_path + @"\statistic_sd", @"分组字段", "SUM_" + area_field);// 加入行政代码和名称dict.Add("名称", ad_name);dict.Add("代码", ad_code);// 属性映射到ExcelToolManager.ExcelAttributeMapperCol(new_sheet_path, 2, start_row, dict);// 进入下一个行政单位start_row++;}// 删除空列ToolManager.DeleteNullColByList(new_sheet_path, row_list);// 删除行ToolManager.DeleteRows(new_sheet_path, new List() { 6, 2 });
}// 删除中间数据
Arcpy.Delect(gdb_path + @"\statistic_sd");
File.Delete(folder_path + @"\三调用地自转换.xlsx");
Directory.Delete(new_gdb, true);
foreach (var item in gdb_list)
{Directory.Delete(item, true);
}

三、工具文件分享

我把工具都集合成工具箱，不再单独放单个工具，可以到这里下载完整工具箱，会不断更新：

【ArcGIS Pro二次开发】：CC工具箱https://blog.csdn.net/xcc34452366/article/details/131506345PS：可以直接点击...bin\Debug\net6.0-windows\下的.esriAddinX文件直接安装。

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！