CRSuggest 使用教程

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 提供了一种方便的方式来处理和优化空间数据，以便更好地进行空间分析和决策制定。