摘要
生成要素类中几何问题的报告。
有效的输入格式包括 shapefile 以及存储在个人地理数据库或文件地理数据库中的要素类。而企业级地理数据库则会在上传几何时自动检查每个几何的有效性;因此,检查几何和修复几何工具无需用于企业级地理数据库。
有关几何问题、这些问题对软件造成的影响以及可能的起因的详细信息,请参阅检查和修复几何。
用法
输出表可记录发现的每个几何问题。如果未发现任何问题,则该表为空。
输出表中包含以下字段:
- CLASS - 发现问题的要素类的完整路径和名称。
- FEATURE_ID - 存在几何问题的要素的要素 ID (FID) 或对象 ID (OID)。
- PROBLEM - 有关问题的简要描述。
PROBLEM 字段将包含下列内容之一:
- Short segment - 有些线段短于与几何关联的空间参考的系统单位所允许的长度。
- Null geometry - 该要素不具有几何或者 SHAPE 字段为空。
- Incorrect ring ordering - 面的拓扑结构比较简单,但是面的环可能未正确定向(外环 - 顺时针,内环 - 逆时针)。
- Incorrect segment orientation - 各条线段的定向不一致。线段 i 的终点应该与线段 i+1 的起点相接。
- Self intersections - 面不能与自身相交。
- Unclosed rings - 环中最后一条线段的终点必须与第一条线段的起点相接。
- Empty parts - 几何具有多个部分,其中一个部分为空(没有几何)。
- Duplicate vertex - 几何的两个或多个折点坐标相同。
- Mismatched attributes - 某线段端点的 z 坐标或 m 坐标与下一条线段中与之重合的端点的 z 坐标或 m 坐标不匹配。
- Discontinuous parts - 几何的某部分由断开的或不连续的部分组成。
- Empty Z values - 几何的一个或多个折点 z 值为空(例如,NaN)。
- Bad envelope - 包络矩形与几何的坐标范围不匹配。
- Bad dataset extent - 数据集的范围不包括所有要素。
该工具识别出的问题可通过以下方式得到解决:
- 手动编辑和修复存在几何问题的要素。某些问题无法通过编辑进行修复。
- 对识别出存在几何问题的要素类运行修复几何工具。
对于点要素,只会存在空几何问题。
为方便查看报告存在几何问题的要素,可使用添加连接工具以及输入的 OBJECTID 或 FID 字段和输出表的 FEATURE_ID 字段将输入要素连接到输出表。
语法
CheckGeometry_management (in_features, out_table)
参数 | 说明 | 数据类型 |
in_features [in_feature,...] | 要检查是否存在几何问题的一个或多个要素类或要素图层。有效的输入包括 shapefile 以及存储在个人地理数据库或文件地理数据库中的要素类。 | Feature Layer |
out_table | 输出表中包含在输入要素中发现的几何问题。 | Table |
代码示例
CheckGeometry 示例 1(Python 窗口)
以下 Python 窗口脚本演示了如何在即时模式下使用 CheckGeometry 函数。
import arcpy
arcpy.env.workspace = "c:/data/data.gdb"
arcpy.CheckGeometry_management(["contours", "roads", "vegetation"], "CheckGeom_Result")
CheckGeometry 示例 2(独立脚本)
以下独立脚本通过在地理数据库的所有要素类中循环来使用 CheckGeometry 函数。
# BatchCheckGeometry.py
# Description:
# Loops through all the feature classes in a geodatabase, and generates
# a report of the problems encountered with feature geometry.
# Import modules
import arcpy
import os
# The geodatabase in which the feature classes will be checked
arcpy.env.workspace = "C:\\data\\St_Lucia.gdb"
out_table = "checkGeometryResult"
# A variable that will hold the list of all the feature classes
# inside the geodatabase
fc_list = []
# Identify all feature classes in the geodatabase
for path, dirnames, fcs in arcpy.da.Walk(arcpy.env.workspace,
datatype='FeatureClass'):
for fc in fcs:
fc_list.append(os.path.join(path, fc))
print("Running the check geometry tool on {} feature classes".format(
len(fc_list)))
arcpy.CheckGeometry_management(fc_list, out_table)
print("{} geometry problems found, see {} for details.".format(
arcpy.GetCount_management(out_table)[0], out_table))
环境
许可信息
- ArcGIS Desktop Basic: 是
- ArcGIS Desktop Standard: 是
- ArcGIS Desktop Advanced: 是