ArcMap

检查几何

描述

生成要素类中几何问题的报告。

有效的输入格式包括 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 字段将输入要素连接到输出表。

语法

arcpy.management.CheckGeometry(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))

环境

许可信息

  • Basic: 是
  • Standard: 是
  • Advanced: 是

相关主题