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
函数的 gcs
和 units
参数来筛选。
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
提供了一种方便的方式来处理和优化空间数据,以便更好地进行空间分析和决策制定。
热门项目推荐
相关项目推荐
- DDeepSeek-R1-0528DeepSeek-R1-0528 是 DeepSeek R1 系列的小版本升级,通过增加计算资源和后训练算法优化,显著提升推理深度与推理能力,整体性能接近行业领先模型(如 O3、Gemini 2.5 Pro)Python00
cherry-studio
🍒 Cherry Studio 是一款支持多个 LLM 提供商的桌面客户端TSX028unibest
unibest - 最好用的 uniapp 开发框架。unibest 是由 uniapp + Vue3 + Ts + Vite5 + UnoCss + WotUI 驱动的跨端快速启动模板,使用 VS Code 开发,具有代码提示、自动格式化、统一配置、代码片段等功能,同时内置了大量平时开发常用的基本组件,开箱即用,让你编写 uniapp 拥有 best 体验。TypeScript01
热门内容推荐
1 freeCodeCamp 优化测验提交确认弹窗的用户体验2 freeCodeCamp课程页面空白问题的技术分析与解决方案3 freeCodeCamp课程中Todo应用测试用例的优化建议4 freeCodeCamp课程中CSS背景与边框测验的拼写错误修复5 freeCodeCamp论坛排行榜项目中的错误日志规范要求6 freeCodeCamp全栈开发课程中"午餐选择器"项目的教学方法优化7 freeCodeCamp课程中CSS模态框描述优化分析8 freeCodeCamp JavaScript课程中十进制转二进制转换器的潜在问题分析9 freeCodeCamp 实验室项目:Event Hub 图片元素顺序优化指南10 freeCodeCamp课程中sr-only类与position: absolute的正确使用
最新内容推荐
项目优选
收起

🎉 (RuoYi)官方仓库 基于SpringBoot,Spring Security,JWT,Vue3 & Vite、Element Plus 的前后端分离权限管理系统
Vue
416
317

React Native鸿蒙化仓库
C++
90
157

openGauss kernel ~ openGauss is an open source relational database management system
C++
46
115

🔥LeetCode solutions in any programming language | 多种编程语言实现 LeetCode、《剑指 Offer(第 2 版)》、《程序员面试金典(第 6 版)》题解
Java
50
13

本仓将收集和展示高质量的仓颉示例代码,欢迎大家投稿,让全世界看到您的妙趣设计,也让更多人通过您的编码理解和喜爱仓颉语言。
Cangjie
268
402

🍒 Cherry Studio 是一款支持多个 LLM 提供商的桌面客户端
TSX
309
28

旨在打造算法先进、性能卓越、高效敏捷、安全可靠的密码套件,通过轻量级、可剪裁的软件技术架构满足各行业不同场景的多样化要求,让密码技术应用更简单,同时探索后量子等先进算法创新实践,构建密码前沿技术底座!
C
87
238

本仓将收集和展示仓颉鸿蒙应用示例代码,欢迎大家投稿,在仓颉鸿蒙社区展现你的妙趣设计!
Cangjie
342
213

前端智能化场景解决方案UI库,轻松构建你的AI应用,我们将持续完善更新,欢迎你的使用与建议。
官网地址:https://matechat.gitcode.com
625
73

🎉 基于Spring Boot、Spring Cloud & Alibaba、Vue3 & Vite、Element Plus的分布式前后端分离微服务架构权限管理系统
Vue
85
61