通过非SF列内部连接两个SF对象 [英] Inner-Joining two sf objects by non sf column
本文介绍了通过非SF列内部连接两个SF对象的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
CHECK_JOIN(x,y)出错: y应为data.frame;对于空间连接,请使用st_joinfalse
可重现示例如下:
df1 <- data.frame(
var = c("a", "b", "c"),
lon1 = c(20,35,45),
lat1 = c(50,10,15)
) %>% st_as_sf(coords = c("lon1", "lat1"), dim = "XY") %>%
st_set_crs(4326)
df2 <- data.frame(
var = c("a", "b", "c"),
lon2 = c(15,25,35),
lat2 = c(5,10,15)
) %>% st_as_sf(coords = c("lon2", "lat2"), dim = "XY") %>%
st_set_crs(4326)
df <- inner_join(df1, df2, by = "var")
我不喜欢丢弃几何体,因为我认为那样会在以后弄乱我的结果,但欢迎任何解决方案
推荐答案
如果您只想执行非空间联接,但将几何列向前推进,则可以首先"停用"它们(例如as.data.frame()),联接,然后"重新激活"要激活的几何列。现在,您的data.frame中有两个SFC列,但"geometry.x"是活动的列。
df <- inner_join(df1 %>% as.data.frame(), df2 %>% as.data.frame(), by = "var")
df %<>% st_sf(sf_column_name = 'geometry.x')
> str(df)
Classes ‘sf’ and 'data.frame': 3 obs. of 3 variables:
$ var : Factor w/ 3 levels "a","b","c": 1 2 3
$ geometry.x:sfc_POINT of length 3; first list element: Classes 'XY', 'POINT', 'sfg' num [1:2] 20 50
$ geometry.y:sfc_POINT of length 3; first list element: Classes 'XY', 'POINT', 'sfg' num [1:2] 15 5
- attr(*, "sf_column")= chr "geometry.x"
- attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA
..- attr(*, "names")= chr "var" "geometry.y"
> st_crs(df)
Coordinate Reference System:
EPSG: 4326
proj4string: "+proj=longlat +datum=WGS84 +no_defs"
这篇关于通过非SF列内部连接两个SF对象的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!
查看全文