我是初学者,我在我的问题上做了很多工作,但我仍然被卡住了......
我有一个包含世界温度值的矩阵。尺寸为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)
你知道我该怎么做吗?
非常感谢
答案 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)。单元格大小是隐式的,并根据范围和行/列计算。