Question

我意识到这个问题的标题与另一个问题相似，但是那里的解决方案对我不起作用。

我有一个光栅文件。我有一个纬度矩阵和一个经度矩阵。我想将栅格文件的x和y坐标替换为两个矩阵中的坐标。

我在这里上传了示例文件。一共有三个文件，一个是栅格文件（GeoTIFF），两个是经度和纬度矩阵。

exec_program:  program.cpp
    g++ -o program program.cpp -lm -std=c++11 -O3 -I/usr/include/python3.6m/
    ./program
    rm program

https://guides.rubyonrails.org/action_cable_overview.html#in-app

执行此操作的简单方法是什么？我可以在RasterBrick上使用相同的方法吗？

我尝试了很多事情，但找不到如何替换坐标。

任何帮助将不胜感激。

Answer 1

类似的事情可能会有所帮助：

library(raster)
library(rgdal)

# starting raster:
r<- raster(matrix(1:100, nrow = 10))

# new coordinates:
x <- 51:60
y <- 21:30
xy <- expand.grid(x, y)

# create new raster via SpatialPixelsDataFrame, copying the values from the original raster but supplying the new coordinates
r2 <- raster(SpatialPixelsDataFrame(xy, data.frame(values(r))))

# plot to have a look:
par(mfrow = c(1,2))
plot(r)
plot(r2)

请注意，栅格以xy坐标为像素中心，给出的范围可能会让您感到惊讶：

extent(r2)

但是您可以轻松调整以适应任何需要

如何替换R中栅格的XY坐标？

1 个答案: