首页
/ CRSuggest 使用教程

CRSuggest 使用教程

2025-04-17 12:25:45作者:薛曦旖Francesca

1. 项目介绍

crsuggest 是一个为R语言编写的开源包,旨在帮助空间分析人员为他们的空间数据确定一个合适的投影坐标参考系统(CRS)。该包通过尝试匹配输入的空间数据集与相应的坐标参考系统,以用于地图绘制或空间分析。crsuggest 使用来自 EPSG 注册表的数据,该注册表是国际油汽生产者协会的产品,并遵守其数据使用条款。

2. 项目快速启动

在R环境中安装 crsuggest 包:

install.packages("crsuggest")

或者,从GitHub获取开发版本:

remotes::install_github("walkerke/crsuggest")

3. 应用案例和最佳实践

分析未投影数据

如果你有一个地理坐标系中的点数据集,并希望绘制这些点的缓冲区,你可以使用 suggest_crs 函数来推荐合适的坐标系统。

library(tigris)
library(sf)

# 示例:获取德克萨斯州的地理标记点
tx_landmarks <- landmarks("TX")

# 由于未投影,直接使用st_buffer会报错
# landmark_buf <- st_buffer(tx_landmarks, 1)

# 使用crsuggest来推荐合适的CRS
library(crsuggest)
possible_crs <- suggest_crs(tx_landmarks)

# 查看推荐的CRS
glimpse(possible_crs)

# 选择一个CRS并重新计算缓冲区
landmarks_projected <- st_transform(tx_landmarks, 3081)
buffer_1km <- st_buffer(landmarks_projected, 1000)

定制参数

如果你已经知道你想要一个使用特定地理坐标系(例如NAD 1983,EPSG代码4269)和特定测量单位(例如美国英尺)的投影坐标系统,你可以使用 suggest_crs 函数的 gcsunits 参数来筛选。

library(tidycensus)
library(ggplot2)
library(sf)

# 示例:获取Cook县伊利诺伊州的人口统计数据
cook_age <- get_acs(geography = "tract", variables = "B01002_001", state = "IL", county = "Cook County", geometry = TRUE)

# 使用crsuggest来推荐合适的CRS
cook_crs <- suggest_crs(cook_age, gcs = 4269, units = "us-ft")

# 查看推荐的CRS
glimpse(cook_crs)

# 使用推荐的CRS进行绘图
ggplot(cook_age, aes(fill = estimate)) +
  geom_sf(color = NA) +
  coord_sf(crs = 3435) +
  scale_fill_viridis_c()

猜测数据集的CRS

如果分析师使用的数据集没有提供坐标系统信息(例如缺少 .prj 文件的Shapefile或CSV文件),可以使用 guess_crs 函数来帮助确定数据的坐标系统。

library(sf)
library(mapview)

# 示例:创建一个坐标数据框
locations <- data.frame(
  X = c(1200822.97857801, 1205015.51644983, 1202297.44383987, 1205877.68696743, 1194763.21511923, 1195463.42403192, 1199836.01037452, 1207081.96500368, 1201924.15986897),
  Y = c(1246476.31475063, 1248612.72571423, 1241479.45996392, 1243898.58428024, 1246033.7550009, 1241827.7730307, 1234691.50899912, 1251125.67808482, 1252188.4333016),
  id = 1:9
)

# 将数据框转换为sf对象
locations_sf <- st_as_sf(locations, coords = c("X", "Y"))

# 使用mapview查看数据
mapview(locations_sf)

# 使用guess_crs猜测CRS
guessed_crs <- guess_crs(locations_sf)

4. 典型生态项目

crsuggest 可以被广泛用于需要确定合适投影系统的空间数据分析项目,特别是在以下场景中:

  • 地理信息系统(GIS)开发
  • 空间数据转换和集成
  • 地图制作和可视化
  • 空间统计分析

这些场景中,crsuggest 提供了一种方便的方式来处理和优化空间数据,以便更好地进行空间分析和决策制定。

项目优选

收起
RuoYi-Vue3RuoYi-Vue3
🎉 (RuoYi)官方仓库 基于SpringBoot,Spring Security,JWT,Vue3 & Vite、Element Plus 的前后端分离权限管理系统
Vue
416
317
ohos_react_nativeohos_react_native
React Native鸿蒙化仓库
C++
90
157
openGauss-serveropenGauss-server
openGauss kernel ~ openGauss is an open source relational database management system
C++
46
115
leetcodeleetcode
🔥LeetCode solutions in any programming language | 多种编程语言实现 LeetCode、《剑指 Offer(第 2 版)》、《程序员面试金典(第 6 版)》题解
Java
50
13
Cangjie-ExamplesCangjie-Examples
本仓将收集和展示高质量的仓颉示例代码,欢迎大家投稿,让全世界看到您的妙趣设计,也让更多人通过您的编码理解和喜爱仓颉语言。
Cangjie
268
402
cherry-studiocherry-studio
🍒 Cherry Studio 是一款支持多个 LLM 提供商的桌面客户端
TSX
309
28
openHiTLSopenHiTLS
旨在打造算法先进、性能卓越、高效敏捷、安全可靠的密码套件,通过轻量级、可剪裁的软件技术架构满足各行业不同场景的多样化要求,让密码技术应用更简单,同时探索后量子等先进算法创新实践,构建密码前沿技术底座!
C
87
238
HarmonyOS-ExamplesHarmonyOS-Examples
本仓将收集和展示仓颉鸿蒙应用示例代码,欢迎大家投稿,在仓颉鸿蒙社区展现你的妙趣设计!
Cangjie
342
213
MateChatMateChat
前端智能化场景解决方案UI库,轻松构建你的AI应用,我们将持续完善更新,欢迎你的使用与建议。 官网地址:https://matechat.gitcode.com
625
73
RuoYi-Cloud-Vue3RuoYi-Cloud-Vue3
🎉 基于Spring Boot、Spring Cloud & Alibaba、Vue3 & Vite、Element Plus的分布式前后端分离微服务架构权限管理系统
Vue
85
61