【ArcGIS Pro二次开发】(48)：三调土地利用现状分类面积汇总统计

之前做了一个三调三大类面积统计，有小伙伴反映太粗糙，想要一个完整的地类面积汇总表。

【ArcGIS Pro二次开发】(35)：三调三大类面积统计

本质上并没有多少难度，之前也做过类似的用地用海汇总表，于是拿出来改一改就好了。

一、要实现的功能

如上图所示，点击【三调】下的【统计三调地类】工具。

在弹出的工具框中，分别输入三调图斑，要统计的面积字段，以及要输出Excel表格的位置。

三调图层注意要确保图层中有【DLMC】字段，这是三调自带的字段。

面积字段，以投影面积为例，选择【shape_area】，如果要统计椭球面积，就选择相应的面积字段。

点击执行即可，生成结果如下：

二、实现流程

这里表格中的地类我是按【TDT 1055-2019 第三次全国国土调查技术规程】中的用地表录入的。但是实际使用过程中，发现有些用地在规范之外，比如【1104A养殖坑塘】等，目前我只能按自己碰到的情况先做一个修正，比如【1104A养殖坑塘】就归类到【1104坑塘水面】。

我对三调业务不熟，如果有熟悉业务的请多指导。

                    // 打开GDB数据库using Geodatabase gdb = new Geodatabase(new FileGeodatabaseConnectionPath(new Uri(gdb_path)));// 获取要素类using FeatureClass featureClasse = gdb.OpenDataset<FeatureClass>("fc_sd");// 逐行游标using (RowCursor rowCursor = featureClasse.Search(null, false)){while (rowCursor.MoveNext()){using (Row row = rowCursor.Current){var mc = row["DLMC"].ToString();// 赋值if (mc == "养殖坑塘"){row["DLMC"] = "坑塘水面";}else if (mc.Contains("可调整")){row["DLMC"] = mc.Replace("可调整", "");}// 保存row.Store();}}}

然后还要归纳一下大类，这里不像用地用海，可以从小类或中类直接取编码前两位归纳大类，只能做一个对照表，然后属性映射。

        // 属性映射public static void AttributeMapper(string in_data, string in_field, string map_field, string map_tabel){// 获取连接表的2个字段名string exl_field01 = GetCellFromExcel(map_tabel, "A1");string exl_field02 = GetCellFromExcel(map_tabel, "B1");List<string> fields = new List<string>() { exl_field02 };// 连接字段Arcpy.JoinField(in_data, in_field, map_tabel, exl_field01, fields);// 计算字段Arcpy.CalculateField(in_data, map_field, "!" + exl_field02 + "!");// 删除多余字段Arcpy.DeleteField(in_data, fields);}

表格像这样：

然后汇总一下【DLMC】和映射后的大类面积：

        // 汇总统计加强版public static void MultiStatistics(string in_table, string out_table, string statistics_fields, List<string> case_fields, string total_field = "合计", int unit = 0, bool is_output = false){try{List<string> list_table = new List<string>();for (int i = 0; i < case_fields.Count; i++){Arcpy.Statistics(in_table, out_table + i.ToString(), statistics_fields, case_fields[i]);    // 调用GP工具【汇总】Arcpy.AlterField(out_table + i.ToString(), case_fields[i], @"分组字段", @"分组字段");  // 调用GP工具【更改字段】list_table.Add(out_table + i.ToString());}Arcpy.Statistics(in_table, out_table + "_total", statistics_fields, "");    // // 调用GP工具【汇总】Arcpy.AddField(out_table + "_total", @"分组字段", "TEXT");    // 调用GP工具【更改字段】Arcpy.CalculateField(out_table + "_total", @"分组字段", "\"" + total_field + "\"");    // 调用GP工具【计算字段】list_table.Add(out_table + "_total");     // 加入列表// 合并汇总表Arcpy.Merge(list_table, out_table, is_output);       // 调用GP工具【合并】// 转换为公顷if (unit > 0){ChangeUnit(out_table, "SUM_Shape_Area", unit);        // 单位转换}}catch (Exception ee){MessageBox.Show(ee.Message + ee.StackTrace);return;}}// 调用
List<string> list_bm = new List<string>() { "DLMC", "mc_dl" };
ToolManager.MultiStatistics(gdb_path + @"\fc_sd", gdb_path + @"\statistic_sd", bm_field + " SUM", list_bm, "国土调查总面积", 1, true);

最后属性映射到Excel表格中，再删除0值行，这部分详细的可以看之前的用地用海汇总表，这里就不展开放代码了：

                    // 将映射属性表中获取字典DictionaryDictionary<string, string> dict = ToolManager.GetDictFromTable("statistic_sd", @"分组字段", "SUM_" + bm_field);// 属性映射大类ToolManager.ExcelAttributeMapper(excel_path + @"\sheet1$", 4, 5, dict, 5);pw.AddProcessMessage(20, time_base, "删除0值行");// 删除0值行ToolManager.Delete0Row(excel_path + @"\sheet1$", 5, 5, true);

以上便是工具的核心代码。