更改 R 上的 matrix/raster 坐标
Change matrix/raster coordinates on R
我是初学者,我已经为我的问题做了很多工作,但我仍然卡住了...
我有一个包含世界温度值的矩阵。尺寸为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
...
这个矩阵 "looks like" 一个光栅。我的问题是我不知道如何告诉 R 矩阵分辨率。
如果我然后用 plot(raster(matrix))
栅格化和绘图,我得到这个:http://postimg.org/image/gdnblob07/
然而,我需要轴为 -180;180 (x) 和 -90;90 (y)
你知道我该怎么做吗?
非常感谢
使用 raster
包是正确的直觉,但是当您仅使用矩阵实例化它时,它无法准确获得范围。这样做就可以了:
library(raster)
r <- raster(nrow=360,ncol=720,vals=matrix)
plot(r)
raster
将使用(默认)世界边界初始化光栅对象。您可以选择指定坐标参考系统(例如,crs="+proj=longlat +datum=WGS84")和范围(ymx、ymn、xmx、xmn)。像元大小是隐式的,是根据范围和 rows/cols 计算得出的。
我是初学者,我已经为我的问题做了很多工作,但我仍然卡住了...
我有一个包含世界温度值的矩阵。尺寸为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
...
这个矩阵 "looks like" 一个光栅。我的问题是我不知道如何告诉 R 矩阵分辨率。
如果我然后用 plot(raster(matrix))
栅格化和绘图,我得到这个:http://postimg.org/image/gdnblob07/
然而,我需要轴为 -180;180 (x) 和 -90;90 (y)
你知道我该怎么做吗?
非常感谢
使用 raster
包是正确的直觉,但是当您仅使用矩阵实例化它时,它无法准确获得范围。这样做就可以了:
library(raster)
r <- raster(nrow=360,ncol=720,vals=matrix)
plot(r)
raster
将使用(默认)世界边界初始化光栅对象。您可以选择指定坐标参考系统(例如,crs="+proj=longlat +datum=WGS84")和范围(ymx、ymn、xmx、xmn)。像元大小是隐式的,是根据范围和 rows/cols 计算得出的。