近邻分析 (Analysis)

许可等级:BasicStandardAdvanced

摘要

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

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

插图

Finding near feature by geometry type

用法

语法

Near_analysis (in_features, near_features, {search_radius}, {location}, {angle}, {method})
参数说明数据类型
in_features

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

Feature Layer
near_features
[near_features,...]

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

Feature Layer
search_radius
(可选)

用于搜索邻近要素的半径。如果未指定任何值,则会考虑所有邻近要素。如果指定了距离,但没有指定任何单位或将单位设置为未知,则将使用输入要素的坐标系单位。如果使用了 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
(可选)

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

  • PLANAR在要素间使用平面距离。这是默认设置。
  • GEODESIC在要素间使用测地线距离。这种方法考虑到了椭球体的曲率,并能够正确处理日期变更线和两极附近的数据。
String

代码实例

近邻分析 (Near) 示例 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')
近邻分析 (Near) 示例 2(独立 Python 脚本)

以下 Python 脚本演示了如何在独立脚本中使用近邻分析 (Near) 函数。

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

import arcpy
from arcpy import env

# Set workspace environment
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 ex:
    print ex.args[0]
近邻分析 (Near) 示例 3(独立 Python 脚本)

以下 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.
近邻分析 (Near) 示例 4(独立 Python 脚本)

演示使用“邻近分析”工具的派生输出进行后处理。脚本为每个邻近要素查找最接近的输入要素。

# 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]}

环境

相关主题

许可信息

ArcGIS for Desktop Basic: 否
ArcGIS for Desktop Standard: 否
ArcGIS for Desktop Advanced: 是
5/10/2014