0

我正在尝试从 .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)
4

1 回答 1

1

您可以尝试flipraster包中获取功能。

library(raster)

r <- raster(nrow=18, ncol=36)
r[] <- 1:ncell(r)

par(mfrow = c(1, 2))
plot(r)
plot(flip(r, direction = "y"))

在此处输入图像描述

于 2013-11-04T13:28:38.963 回答