邻近 (分析)

描述

可计算输入要素与其他图层或要素类中的最近要素之间的距离和其他邻近性信息。

了解有关如何使用地理处理工具计算邻近值的详细信息

插图

按几何类型查找邻近要素

使用方法

  • 将以下字段添加到输入中。如果字段已存在,则将更新字段值。这些字段名称可以使用字段名称 (field_names) 参数进行配置。

    • NEAR_FID - 最邻近要素的 ObjectID。如果未发现邻近要素,则该值为 -1。
    • NEAR_DIST - 输入要素与邻近要素之间的距离。该值采用输入要素坐标系的线性单位,如果将 Method 参数设置为 Geodesic 且输入位于地理坐标系中,则该值将采用“米”作为单位。如果未发现邻近要素,则该值为 -1。
    • NEAR_FC - 包含邻近要素的要素类的目录路径。仅在指定了多个邻近要素时,才会将此字段添加到输出表中。如果未发现邻近要素,则该值为空字符串或为空。

    如果选中 Location 参数(在 Python 中 location 参数设置为 LOCATION),则以下字段将被添加到输入要素中。如果字段已存在,将更新字段值。字段值单位取决于为 Method 参数选择的方法。如果设置为 Planar,则该值为输入要素坐标系的线性单位。如果设置为 Geodesic,则该值在与输入要素坐标系相关联的地理坐标系中。

    字段名称描述
    NEAR_X

    邻近要素中距离输入要素最近位置的 x 坐标。如果未发现邻近要素,则该值为 -1。

    NEAR_Y

    邻近要素中距离输入要素最近位置的 y 坐标。如果未发现邻近要素,则该值为 -1。

    如果选中 Angle 参数(在 Python 中 angle 参数设置为 ANGLE),则以下字段将被添加到输入要素中。如果字段已存在,将更新字段值。

    字段描述
    NEAR_ANGLE

    连接输入要素和邻近要素的位于 FROM_XFROM_Y 位置的线的角度。如果未发现邻近要素或邻近要素与输入要素相交,则该值为 0。

  • 如果在搜索半径内未找到任何要素,则 NEAR_FIDNEAR_DIST 的值都将为 -1。

  • 输入要素和邻近要素均可为点、多点、线或面。

  • 邻近要素可包括不同形状类型(点、多点、线或面)的一个或多个要素类。

  • 同一要素类或图层可同时用作输入要素和邻近要素。这种情况下,所评估的输入要素将被排除在邻近要素候选项之外,以避免得出所有要素都与其自身最接近的结果。

  • 输入要素可以是您已执行要素选择的图层。使用工具执行操作时将使用并更新所选要素。其余要素会将新建字段(例如 NEAR_FIDNEAR_DIST)的值设置为 -1。

  • 如果有多个邻近要素与输入要素之间有相同的最短距离,则随机选择其中一个邻近要素作为最近要素。

  • 使用方法参数的平面选项时,输入要素必须采用适合距离测量值的投影,如等距投影。

    了解有关坐标系和投影的详细信息。

  • 为了显示 FROM_XFROM_YNEAR_XNEAR_Y 位置,输出表可用作创建 XY 事件图层XY 转线工具的输入。

语法

arcpy.analysis.Near(in_features, near_features, {search_radius}, {location}, {angle}, {method}, {field_names})
参数说明数据类型
in_features

输入要素可以是点、折线、面或多点类型。

Feature Layer
near_features
[near_features,...]

一个或多个包含邻近要素候选项的要素图层或要素类。邻近要素可以是点、折线、面或多点。如果指定了多个图层或要素类,则 NEAR_FC 字段将被添加到输入表中,并将存储含有找到的最近要素的源要素类的路径。同一要素类或图层可同时用作输入要素和邻近要素。

Feature Layer
search_radius
(可选)

用于搜索邻近要素的半径。如果未指定任何值,则会考虑所有邻近要素。如果指定了距离,但没有指定任何单位或将单位设置为未知,则将使用输入要素的坐标系单位。如果 Method 参数(Python 中的 method='GEODESIC')使用了 Geodesic 选项,则使用线性单位(如公里或英里)。

Linear Unit
location
(可选)

指定是否将邻近要素上最近位置的 x 和 y 坐标写入 NEAR_XNEAR_Y 字段。

  • NO_LOCATION 将不会写入位置信息。这是默认设置。
  • LOCATION 将写入位置信息。
Boolean
angle
(可选)

指定是否计算邻近角并将其写入输出表的 NEAR_ANGLE 字段。邻近角测量与直线(该直线连接输入要素与其最近要素的最近位置)方向之间的夹角。在 method 参数中使用 PLANAR 方法时,角度在 -180° 到 -180° 的范围内,0° 代表东,90° 代表北,180°(或 -180°)代表西,-90° 代表南。使用 GEODESIC 方法时,角度在 -180° 到 180° 的范围内,0° 代表北,90° 代表东,180°(或 -180°)代表南,-90° 代表西。

  • NO_ANGLE将不写入邻近角值。这是默认设置。
  • ANGLE将邻近角值写入 NEAR_ANGLE 字段。
Boolean
method
(可选)

确定是使用椭球体上的最短路径(测地线)还是使用地平(平面)方法。建议将测地线方法用于在不适合进行距离测量的坐标系(例如 Web 墨卡托或任何地理坐标系)中存储的数据,以及任何地理区域跨度较大的分析。

  • PLANAR要素之间将使用平面距离。这是默认设置。
  • GEODESIC要素之间将使用测地线距离。这种方法考虑到椭球体的曲率,并可以正确处理日期变更线和两极附近的数据。
String
field_names
[[property, fieldname],...]
(可选)

指定将在处理过程中添加的属性字段的名称。

如果不使用此参数或将从此参数中排除要添加的任何字段,则将使用默认字段名称。

默认情况下,将始终添加 NEAR_FIDNEAR_DIST 字段;启用位置参数(Python 中的 location)时,将添加 NEAR_XNEAR_Y 字段;启用角度参数(Python 中的 angle)时,将添加 NEAR_ANGLE 字段;使用多个输入时,将添加 NEAR_FC 字段。

Value Table

派生输出

名称说明数据类型
out_feature_class

更新后的输入要素。

要素图层

代码示例

近邻分析示例 1(Python 窗口)

以下 Python 交互式窗口脚本演示了如何在即时模式下使用 Near 函数。

import arcpy
arcpy.env.workspace = "C:/data/city.gdb" 

## find the nearest road from each house
arcpy.Near_analysis('houses', 'roads')
近邻分析示例 2(独立脚本)

以下 Python 脚本演示了如何在独立脚本中使用 Near 函数。

# Name: Near.py
# Description: Finds nearest features from input feature class to near feature class.

import arcpy

# Set workspace environment
arcpy.env.workspace = "C:/data/city.gdb"

try:
    # set local variables
    in_features = "houses"
    near_features = "parks"
    
    # find features only within search radius
    search_radius = "5000 Meters"
    
    # find location nearest features
    location = "LOCATION"
    
    # avoid getting angle of neares features
    angle = "NO_ANGLE"
    
    # execute the function
    arcpy.Near_analysis(in_features, near_features, search_radius, location, angle)
    
    # get geoprocessing messages
    print(arcpy.GetMessages())

except arcpy.ExecuteError:
    print(arcpy.GetMessages(2))
    
except Exception as err:
    print(err.args[0])
近邻分析示例 3(独立脚本)

以下 Python 脚本演示了如何将邻近角转换为方位角。

# Name: near_angle_to_azimuth.py

import arcpy

# Near tool does not calculate angle in azimuths
# This script, using the output of the tool, converts the angles
# to azimuth angles.

in_table = r"C:/data/city.gdb/near_table"
    
angles = arcpy.da.SearchCursor(in_table, ["NEAR_ANGLE"])

azimuth_angles = []

with angles as rows:
    for row in rows:
        angle = row[0]
        if angle <= 180 and angle > 90:
            azimuth_angles.append(360.0 - (angle - 90))
        else:
            azimuth_angles.append(abs(angle - 90))

# Use these azimuth angles as necessary.
近邻分析示例 4(独立脚本)

以下 Python 脚本演示了使用 Near 工具的派生输出进行后处理。脚本为每个邻近要素查找最接近的输入要素。

# Name: features_closest_to_input.py

"""
    This script finds, for each input feature, a list of near feature it is closest to.
    If near features 6, 7, 10 are closest to input feature 3 then the 
    resulting dictionary will have a list [6, 7, 10] as value for key 3
"""

import os
import arcpy

in_fc = r"C:\data\cities.gdb\cities_many"

# create a dictionary to hold the list of near ids for each input id
nearest_dict = dict()

with arcpy.da.SearchCursor(in_fc, ["OID@", "NEAR_FID"]) as rows:
    for row in rows:
        nearest_id = row[0]  # get OID value for that row
        input_id = row[1]    # get NEAR_FID value 
        
        if input_id in nearest_dict:
            # if a dict key already exists, append near id to value list for that key
            nearest_dict[input_id].append(nearest_id)
        else:
            # if the key does not exist then create a new list with near id
            # and add it to the dictionary
            nearest_dict[input_id] = [nearest_id]
            
print(nearest_dict)

# output will look like:
# {1: [13, 28], 2: [2, 9, 14, 20, 22], 3: [11, 12, 24, 25]}

许可信息

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

相关主题