R: перезаписать координаты xy растрового слоя

Изменить проекцию растровых данных не так просто, как для точек (и линий, полигонов). Это потому, что если вы вычисляете новые координаты из текущих ячеек, они не будут в обычном растре.

Вы можете использовать функцию projectRaster (растровый пакет), чтобы справиться с этим.

library(raster)
test <- raster('janW1.png')

# In this case, you need to provide the correct crs to your data
# I am guessing. (this would be necessary for spatial data sets)
crs(test) <- '+proj=merc +datum=WGS84'

# you may also need to set the extent to actual coordinate values
# extent(test) <- c( , , ,) 
x <- projectRaster(test, crs='+proj=longlat +datum=WGS84') 

Кроме того, вы можете интерполировать вычисленные вами значения в новый растр. См. Примеры ?raster::interpolate.

Не могли бы вы создать растр с нуля с желаемым разрешением и пространственной протяженностью, а затем импортировать в него свои значения. Чтобы создать растр, вы можете использовать что-то вроде:

# Create a matrix of coordinates that define the limits of your raster

ex <- matrix(c(-20, -9.5, 20.5, 31.5), nrow = 2, ncol = 2, byrow = T)

# Turn those coordinates into an extent

ex <- extent(ex)

# Create a raster with the same dimensions as your original one

r <- raster(nrows = 1617, ncols = 1596)

# Set the extent of your raster

r <- setExtent(r, ex, keepres=F)

Чтобы перенести значения из вашего предыдущего растра в только что созданный растр, вы можете использовать:

test <- raster('janW1.png')

# Create vector of values from test

n <- values(test)

# Give values from test to r

values(r) <- n

Я думаю, что получил правильное разрешение из вашего кода, но вам нужно будет ввести четыре координаты вашего экстента в себя. Пустой растр должен иметь точно такое же разрешение, что и исходный растр, иначе он не будет работать. Пустой растр, который вы создаете, автоматически сохраняется в WGS84, поэтому вы можете перепроецировать его, как только вы загрузите свои данные.