st_join renames geom column to 'geometry' if x is st_sf with only a st_sfc col.

z <- st_sf(pid=5, geom=st_sfc(st_polygon(list(matrix(c(0,0,4,0,4,4,0,4,0,0),ncol=2, byrow=T)))))
x <- st_sf(id=1:2, geom=st_sfc(st_point(1:2), st_point(2:3)))
st_join(x[,'geom'], z)

Author: faridcher

R/join.R

@@ -174,7 +174,8 @@ st_join.sf = function(x, y, join = st_intersects, ..., suffix = c(".x", ".y"),
 	if (inherits(x, "tbl_df") && requireNamespace("dplyr", quietly = TRUE))
 		st_sf(dplyr::bind_cols(x[ix,], y[unlist(i), , drop = FALSE]))
   	else
-		st_sf(cbind(as.data.frame(x)[ix,], y[unlist(i), , drop = FALSE]))	
+	        st_sf(cbind(as.data.frame(x)[ix, ,drop=FALSE], y[unlist(i), , drop = FALSE]))
+	
 }
 
 #' @export