Question

我目前正在尝试在R中的点周围创建（100m）缓冲区。我找到了两个非常有用的答案，但我仍然有问题。这里我找到了答案： Buffer (geo)spatial points in R with gbuffer Create buffer and count points in R

我的问题是，当我创建缓冲区时，我只有一个大缓冲区。

ANFR_IDF是一个数据帧。列“ID”是唯一的ID，列“coordonnees”是坐标。

ANFR_IDF$lat=as.numeric(gsub(",.*$", "", ANFR_IDF$coordonnees))
ANFR_IDF$lon=as.numeric(gsub(".*, ", "", ANFR_IDF$coordonnees))
coordinates( ANFR_IDF ) <- c( "lon", "lat" )
CRS=CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
proj4string(ANFR_IDF)=CRS

ANFR_IDF <- spTransform(ANFR_IDF,CRS)
test2=gBuffer( ANFR_IDF, width=100,id=ANFR_IDF$ID, byid=TRUE )

知道可能出现什么问题？非常感谢！

Answer 1

看起来你有一个地理坐标系，其单位为十进制度，所以当你为100参数提供width时，它会创建一个100十进制度的缓冲区（一个巨大的缓冲区！）。尝试先投射

library(sp)
ANFR_IDF_prj <- spTransform(ANFR_IDF, crs([your preferred projection]))
test2=gBuffer( ANFR_IDF_prj, width=100,id=ANFR_IDF_prj$ID, byid=TRUE )

R缓冲区地理空间点

1 个答案: