当前位置: 首页 > news >正文

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

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

  • 方法 1:使用“空间连接 (Spatial Join)”工具
  • 方法 2:使用“点计数 (Point Count)”工具
  • 方法 3:通过“选择 (Select by Location)”统计
  • 方法 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 脚本可以自动化此过程,更适合批量处理。

脚本示例:

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

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),表示每个区域的站点数量。

参考

http://www.lryc.cn/news/494274.html

相关文章:

  • JavaScript中类数组对象及其与数组的关系
  • 基础入门-Web应用架构搭建域名源码站库分离MVC模型解析受限对应路径
  • C#:时间与时间戳的转换
  • QT的exec函数
  • Css—实现3D导航栏
  • 树莓集团:以人工智能为核心,打造数字化生态运营新典范
  • 2024年首届数证杯 初赛wp
  • 2017 NHOI小学(C++)
  • 【一维DP】【三种解法】力扣983. 最低票价
  • 【头歌实训:递归实现斐波那契数列】
  • IntelliJ IDEA配置(mac版本)
  • CSAPP Cache Lab(缓存模拟器)
  • 【机器学习】机器学习的基本分类-监督学习-逻辑回归-对数似然损失函数(Log-Likelihood Loss Function)
  • 51c自动驾驶~合集35
  • 网络安全体系与网络安全模型
  • antd table 自定义表头过滤表格内容
  • Elasticsearch实战:从搜索到数据分析的全面应用指南
  • BEPUphysicsint定点数3D物理引擎介绍
  • 宠物领养平台构建:SpringBoot技术路线图
  • 解决Flink读取kafka主题数据无报错无数据打印的重大发现(问题已解决)
  • python自动化测开面试题汇总(持续更新)
  • 1-1 Gerrit实用指南
  • docker如何安装redis
  • 省级新质生产力数据(蔡湘杰版本)2012-2022年
  • 【游资悟道】-作手新一悟道心法
  • Diffusion中的Unet (DIMP)
  • 编译以前项目更改在x64下面时报错:函数“PVOID GetCurrentFiber(void)”已有主体
  • 【AIGC】大模型面试高频考点-数据清洗篇
  • 当测试时间与测试资源有限时,你会如何优化测试策略?
  • 基于R语言森林生态系统结构、功能与稳定性分析与可视化