探讨Excel2007与ArcGis10.0结合提取小班四至界限的坐标值
2018-12-18蒋家权秦沂林
蒋家权 秦沂林
摘 要：林木采伐管理实现网上申请、审核和发证后，必须输入小班四至界限的坐标值，提取小班四至界限坐标值的工作至关重要。巧用Excel函数与ArcGis 10.0结合提取小班四至界限的坐标值，快捷且精度高，达到事倍功半的效果。
关键词：Excel2007;ArcGis 10.0;小班界;坐标值
中图分类号：S757.4 文献标识码：A DOI：10.11974/nyyjs.20180932159
启动林木采伐管理系统前，不要求填写小班四至界限的坐标值。但2015年1月1日全面启动运行全国林木采伐管理系统以来，办证时必须输入小班四至界限的坐标值，才能通过林木采伐管理系统的审核。
1 直接求得小班四至界限的点极值
启动ArcMap，将需要提取坐标值的.shp格式的小班面文件添加到工作区，打开属性表，添加Xmin、Xmax、Ymin、Ymax、X、Y等字段，在.shp格式的小班面图层处于编辑状态下，选中Xmin字段列，打开字段计算器选项，点中Python，在编辑区域内编辑Python语句round（！shape.extent.Xmin！，2）”（在英文输入状态编辑下），单击“确定”，直接求得并为2位小位的Xmin值。同方法直接求得Xmax、Ymin、Ymax值。
2 运行要素包络矩形转面、面转线、折点转点工具
因小班四至界限是指小班东西南北方向的顶点，利用平面几何中矩形对边平行原理，故采用要素包络矩形转面工具。其操作步骤是，打开Arctoolbox工具箱→数据管理工具→要素①→要素包络矩形转面，添加需要提取坐标值的shp格式的小班面文件，存储路径可不改，但文件名便于后面操作的区分，建议文件名改为矩A。运行完成后，自动添加到ArcMap工作界面中。在要素工具集中，单击“面转线”工具，选中要面转线的图层是矩A，文件名同理更名为线A。同法按①法操作，选择要素折点转點工具，添加需要提取坐标值的.shp格式的小班面文件，只更改存储文件名为点B后运行。
3 利用相交工具提取相交点
其方法是，打开Arctoolbox工具箱→分析工具→叠加分析→相交，选择线A与点B图层相交，本步骤操作关键是“输出类型”必须选择“point”，否则不是所需要的结果。存储文件更名为点C。在图层点C的属性表中，分别计算输出的相交点“X”“Y”的值。再在图层点C属性表中，选中X字段，输入Python语句“round（！X！，2）”，计算X字段列的数据并保留为2位小数。同样方法求Y字段数据。
4 巧用Excel2007逻辑函数提取小班四至界限坐标值
在Arcmap工作界面，导出属性表，其导出的是DBF文件，用Excel2007打开导出的DBF文件，另存为以.xls扩展名的Excel表后关闭。打开刚存储的Excel表，添加方位字段。输入逻辑函数公式，提取小班四至界限坐标值。
4.1 删除非四至界限的点
在工作表内“小班号”“X”“Y”的数据按高级升序排序。在I2单元格输入逻辑函数公式“=（C2=C3）+（D2=D3）”，在J2单元格输入公式“=IF（E2=C2，"西"，IF（F2=C2，"东"，IF（G2=D2，"南”，IF（H2=D2，”北"，0））））”，在K2单元格输入逻辑函数“=（J2=J3）”后，选中I2、J2、K2单元格，将鼠标指针移到单元格K2的右下角，当指针变成“+”形状时，按住鼠标左键向下拖动填充完后释放，系统自动地对各个单元格的条件进行判断，并显示出判断结果，再是筛选删除。删除在I列显示内容为“2”和“1”的所有行区域，只保留显示内容“0”的所有行区域，删除J列显示内容为“0”的所有行区域，删除K列显示内容为“TRUE”的所有行区域。再次恢复I列、K列公式，确保完全删除非四至界限点。在J1单元格内填入“方位”字段，删除I列、K列后保存工作表。
4.2 用数据透视表将复杂数据生成恰当且适合使用的汇总表
其操作方法是：选中4-1整理保存后的工作表中任意一个单元格，插入→数据透视表→创建数据透视表并设置→确定，创建一个不包含数据的数据透视表，随后把“小班号”放在行字段，“方位”放在列字段，“X”“Y”放在数据区域后形成数据透视表。再复制选择粘贴为数值，按所需整理表格。小班图斑呈三角形或斜楔形的，在4个方位中只有2～3个数据值，是数据透视表内有无数据原因。
5 运用效果
若忠县国有林场282个小班全部森林抚育，直接从ArcMap工作界面逐个提取小班四至界限的坐标值，每个小班至少1min，全部提取完需4h多。按笔者所述方法，不超过20min就完成全场282个小班的四至界限坐标值，不仅提高了工作效率，减少内业工作量，还确保数据准确性。
