【ArcGIS Pro实操第10期】统计某个shp文件中不同区域内的站点数

统计某个shp文件中不同区域内的站点数

  • [方法 1:使用"空间连接 (Spatial Join)"工具](#方法 1:使用“空间连接 (Spatial Join)”工具)
  • [方法 2:使用"点计数 (Point Count)"工具](#方法 2:使用“点计数 (Point Count)”工具)
  • [方法 3:通过"选择 (Select by Location)"统计](#方法 3:通过“选择 (Select by Location)”统计)
  • [方法 4:通过"Python 脚本 (ArcPy)"实现](#方法 4:通过“Python 脚本 (ArcPy)”实现)
  • 参考

在 ArcGIS Pro 中,如果需要统计某个shp文件中不同区域内的站点数,可以通过以下步骤实现。通常,这种任务涉及两个数据集:

1、站点数据:点状数据,例如站点坐标。

2、区域数据:多边形数据,例如行政区划、分区等。

目标是计算每个区域(多边形)中包含的站点数量。

方法 1:使用"空间连接 (Spatial Join)"工具

步骤如下:

1、加载数据:

打开 ArcGIS Pro 并加载两个图层:

  • 站点图层(点数据),如站点位置的 shp 文件。
  • 区域图层(多边形数据),如行政区划的 shp 文件。

2、打开"空间连接 (Spatial Join)"工具:

在"分析"选项卡中,选择"工具 (Tools)"。

搜索并打开 "Spatial Join" 工具。

3、设置参数:

  • 目标要素(Target Features):选择 区域图层(多边形数据)。
  • 连接要素(Join Features):选择 站点图层(点数据)。
  • 连接操作(Join Operation):选择 One-To-One(默认)。
  • 匹配选项(Match Option):选择 Intersects(确保统计所有位于多边形内的点)。

4、运行工具:

运行工具后,输出的多边形图层将包含一个新字段,表示每个区域内的站点数(通常是 Join_Count 字段)。

5、结果:

打开输出图层的属性表,可以在 Join_Count 字段中查看每个区域的站点数量。

方法 2:使用"点计数 (Point Count)"工具

如果仅需要统计点的数量,可以使用专门的点计数工具。

步骤如下:

1、加载数据:

打开 ArcGIS Pro 并加载两个图层:

  • 站点图层(点数据)。
  • 区域图层(多边形数据)。

2、打开"点计数 (Summarize Within)"工具:

在"分析"选项卡中,选择"工具 (Tools)"。

搜索并打开 "Summarize Within" 工具。

3、设置参数:

  • 输入多边形(Input Polygons):选择 区域图层。
  • 输入点(Input Summary Features):选择 站点图层。
  • 输出字段:选择一个字段(如站点 ID 或其他唯一标识),用作计数依据。

示例如下:

4、运行工具:

运行工具后,输出的多边形图层将包含一个新字段,表示每个区域内的点数量。

5、结果:

打开输出图层的属性表,可以查看每个区域的点数量。

方法 3:通过"选择 (Select by Location)"统计

如果不需要生成新的图层,只是想快速统计每个区域的站点数量,可以使用"选择 (Select by Location)"功能。

步骤如下:

1、加载数据:

同样加载站点图层和区域图层。

2、打开"按位置选择 (Select by Location)"工具:

在"地图"选项卡中,选择"按位置选择 (Select by Location)"工具。

3、设置参数:

  • 输入要素(Input Features):选择 站点图层。
  • 选择要素(Selecting Features):选择 区域图层。
  • 空间关系(Spatial Relationship):选择 完全包含 (Completely Within) 或 相交 (Intersects)。

4、统计结果:

在属性表中查看被选中的点的数量。
分别对每个多边形执行上述操作,记录每个区域的站点数量。

方法 4:通过"Python 脚本 (ArcPy)"实现

使用 ArcPy 脚本可以自动化此过程,更适合批量处理。

脚本示例:

以下代码统计每个多边形内的站点数,并将结果写入区域图层的属性表。

python 复制代码
import arcpy

# 输入文件路径
point_layer = r"path_to_point_layer.shp"  # 站点图层路径
polygon_layer = r"path_to_polygon_layer.shp"  # 区域图层路径
output_layer = r"path_to_output_layer.shp"  # 输出图层路径

# 执行空间连接
arcpy.analysis.SpatialJoin(
    target_features=polygon_layer,
    join_features=point_layer,
    out_feature_class=output_layer,
    join_operation="JOIN_ONE_TO_ONE",
    match_option="INTERSECT"
)

print("空间连接完成,输出保存到:", output_layer)

运行该脚本后,output_layer 将包含一个新字段(Join_Count),表示每个区域的站点数量。

参考

相关推荐
杨超越luckly2 天前
HTML应用指南:利用GET请求获取全国奈雪的茶门店位置信息
大数据·前端·python·arcgis·信息可视化·html
zhou_x_b3 天前
解决栅格数据裁剪矢量数据问题两种方法,ArcGIS解决与PYTHON解决
arcgis
维维180-3121-14553 天前
ArcGIS水文及空间分析与SWMM融合协同在城市排水防涝领域中的应用
arcgis·水文·内涝
GIS小小研究僧7 天前
ArcGIS Pro+PS 实现地形渲染效果图
arcgis·gis·qgis·地理信息
杨超越luckly7 天前
Python应用指南:使用PyKrige包实现ArcGIS的克里金插值法
python·算法·arcgis·信息可视化·克里金法
WangYan20228 天前
Python+ArcGIS+AI蒸散发与GPP估算|Penman-Monteith模型|FLUXNET数据处理|多源产品融合|专业科研绘图与可视化等
arcgis·蒸散发·光合作用·植被生产力估算
清纯世纪9 天前
Arcgis 10.7 矢量的分区统计
arcgis
猫头虎9 天前
什么是 npm、Yarn、pnpm? 有什么区别? 分别适应什么场景?
前端·python·scrapy·arcgis·npm·beautifulsoup·pip
杨超越luckly10 天前
ArcGISPro应用指南:使用ArcGIS Pro创建与优化H3六边形网格
arcgis·信息可视化·数据可视化·shp·h3网格
柳晓黑胡椒10 天前
Cesium源码打包
arcgis