R将x或y坐标指定给栅格的单元格以执行计算

Question

是否有任何解决方案可以使用R将X或Y坐标分配给光栅图像的所有单元？

例如，假设我有一个包含3x3单元格的栅格。左下坐标为X = 7，Y = 15（以米为单位），X和Y方向的分辨率均为+10 m（X向右增加，Y向上增加）。然后，我想生成栅格表面，其中每个单元格都具有X和Y值，如下所示：

X raster
7  17  27
7  17  27
7  17  27

Y surface
35  35  35
25  25  25
15  15  15

更新：this是实际的栅格对象。

有没有办法做到这一点？

我尝试了套餐＆＃39; raster＆＃39;，但无法找到解决方案。

任何帮助都将不胜感激。

Answer 1

这是我发现的解决方案，受到@Carl提交的评论的启发。

假设我的栅格名为d，具有以下特征：

class       : RasterLayer
dimensions  : 59, 67, 3953  (nrow, ncol, ncell)
resolution  : 90, 90  (x, y)
extent      : 482855.6, 488885.6, 4763517, 4768827  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-3 +k=1 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : in memory
names       : subset
values      : 328, 1120  (min, max)

我首先使用x坐标值生成一个矩阵：

xm<-matrix(xFromCell(d,c(1:3953)),nrow=59,byrow=TRUE)

然后用矩阵创建了一个栅格：

x<-raster(xm,xmn=482855.6, xmx=488885.6,ymn=4763517,ymx=4768827)

最后，分配了它的预测：

projection(x)<-"+proj=tmerc +lat_0=0 +lon_0=-3 +k=1 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"

我还显示了栅格，看看结果是否正常......

plot(x)

...并阅读栅格描述

x
class       : RasterLayer
dimensions  : 59, 67, 3953  (nrow, ncol, ncell)
resolution  : 90, 90  (x, y)
extent      : 482855.6, 488885.6, 4763517, 4768827  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-3 +k=1 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : in memory
names       : layer
values      : 482900.6, 488840.6  (min, max)

我重复了这些步骤，但在第一步中使用yFromCell。