r - adehabitat 包的位置数据格式

Question

我有一个这种格式的文件:

ASCII format

第一行看起来像这样:

ncols 1440
nrows 720
xllcorner -180.0
yllcorner -90
cellsize 0.25
NODATA_value -9999

基本上，我的世界在 x 方向(经度)上有 1440 个“瓷砖”，在 y 方向(纬度)上有 720 个“瓷砖”。每个“瓷砖”都是一个长度为 0.25 度的正方形。我想我的 xllcorner 和 yllcorner 是正确的。我可以在 R 中像这样绘制这张 map :

library("adehabitat")
bio1 <- import.asc("D:/ENFA/data.asc")
maps <- as.kasc(list(data = bio1))
image(maps, col = cm.colors(256), clfac = list(Aspect = cl))

map 看起来不错。

我想使用 adehabitat 包执行一些生态位因素分析 (ENFA)，但我不太确定位置数据。基本上我现在将它们作为经度和纬度，但我也可以将它们生成为“tile index”(例如，左下角的纬度为 -90 和经度 -180，因此“tile index”将为 0、0 - 对吗？ ).哪种是正确的位置数据格式？我会像这样使用 ENFA 代码:

locs <- read.table("D:/ENFA/Locs.txt", header = TRUE, sep="\t")
dataenfa1 <- data2enfa(maps, locs)
pc <- dudi.pca(dataenfa1$tab, scannf = FALSE)
enfa1 <- enfa(pc, dataenfa1$pr,scannf = FALSE)
hist(enfa1)

如有任何意见，我将不胜感激。提前致谢。

Answer 1

以经纬度形式保留坐标的问题在于，在地球上的大多数地方，经度的长度与纬度的长度不同。这可能会夸大某些方向相对于其他方向的距离，从而扭曲您的 ENFA。

特别是如果您的数据来自相对较小的区域，我建议重新表示沿 W/E x 轴和 S/N y 轴的坐标(以米为单位)。如果你所有的点都落在一个 UTM zone 内，那么您可以在 R 中使用 rgdal 包中的 project() 进行转换:

这是一个例子，taken from here :

library(rgdal)

# Make a two-column matrix, col1 = long, col2 = lat
xy <- cbind(c(118, 119), c(10, 50))
# Convert it to UTM coordinates (in units of meters)
project(xy, "+proj=utm +zone=51 ellps=WGS84")
          [,1]    [,2]
[1,] -48636.65 1109577
[2,] 213372.05 5546301 

有关如何操作空间数据的更多信息，请参阅 Bivand、Pebesma 和 Gomez-Rubio 的“使用 R 进行应用空间数据分析”。如果您需要更具体的帮助，请尝试 R-sig-Geo mailing list .

希望这对您有所帮助。