更改R上的矩阵/栅格坐标

时间:2015-06-02 21:38:43

标签: r matrix coordinates raster

我是初学者,我在我的问题上做了很多工作,但我仍然被卡住了......

我有一个包含世界温度值的矩阵。尺寸为360 * 720。值为-180,-90(经度,纬度)至180,90,分辨率为0.5。

> head(matrix)

NaN.     NaN..1 NaN..2 NaN..3 NaN..4, [...] NaN..718 NaN..719
[1,]      Na      Na      Na      Na           Na      Na
[2,]      Na      Na      Na      Na           Na      Na
[3,]      Na      Na      Na      Na           Na      Na
[4,]      Na      Na      Na      Na           Na      Na
[5,]      Na      Na      Na      Na           Na      Na
[6,]   -1.6634 -1.6634 -1.6634 -1.6634       -1.6634 -1.7469
...

此矩阵“看起来像”一个栅格。我的问题是我不知道如何告诉R矩阵分辨率。

如果我然后使用plot(raster(matrix))进行栅格化和绘图,我会得到:http://postimg.org/image/gdnblob07/ 然而,我需要轴为-180; 180(x)和-90; 90(y)

你知道我该怎么做吗?

非常感谢

1 个答案:

答案 0 :(得分:1)

使用raster包是正确的直觉,但是当你用矩阵实例化它时,它不能准确地获得范围。这样就可以了:

library(raster)
r <- raster(nrow=360,ncol=720,vals=matrix)
plot(r)

raster将使用(默认情况下)世界边界初始化栅格对象。您可以选择指定坐标参考系统(例如,crs =“+ proj = longlat + datum = WGS84”)和范围(ymx,ymn,xmx,xmn)。单元格大小是隐式的,并根据范围和行/列计算。