我正在尝试从 .csv 文件创建空间点数据框。它(通常)有效,但输出是镜像的,所以我用它创建的栅格地图是“颠倒的”。我在这里使用了这种方法(创建 RasterLayer)来创建栅格。
我无法弄清楚问题是什么。任何帮助将非常感激!
base <- read.csv(file="Ancylus_Cleaned_Coordinates.csv", header=T, sep="\t")
head(base)
#latitude longitude species
#1 40.55000 -8.16000 Afluviatilis
#2 40.22000 -8.15000 Afluviatilis
#3 54.43584 -8.12788 Afluviatilis
spdf <- SpatialPointsDataFrame( base[ c("latitude" , "longitude") ], data = data.frame( base$species ), proj4string = CRS("+proj=longlat +datum=WGS84") )
r <- raster(extent(spdf))
r
然后我设置了光栅的分辨率
res(r) <- 1
展开
r <- extend(r, extent(r)+5)
绘制我的子样本
acsel <- gridSample(base, r, n=1)
并绘制结果
p <- rasterToPolygons(r)
plot(p, border='gray')
points(afc)