library(tidyverse)
library(tigris)
library(sf)
santacruz <- tracts("CA", "Santa Cruz")
coords_sf <- locations %>% st_as_sf(coords = c("Longitude", "Latitude"), crs=4269)
这应该具有相同的 CRS,但是当我尝试时
st_intersects(coords_sf, santacruz)
我明白了
错误:st_crs(x) == st_crs(y) 不是 TRUE
然后我尝试了
st_set_crs(santacruz, 4269)
st_set_crs(coords_sf, 4269)
st_transform(santacruz, 4269)
st_transform(coords_sf, 4269)
它不起作用。我也试过
st_transform(santacruz, crs = "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
st_transform(coords_sf, crs = "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
无论我尝试设置 CRS 并在尝试时对其进行转换
st_intersects(coords_sf, santacruz)
我明白了
错误:st_crs(x) == st_crs(y) 不是 TRUE
在这一点上,我无法判断设置 CRS 或转换或st_intersects
函数是否有问题。谢谢,