Растр по вектору с заданным разрешением в километрах

Question

Я хотел бы знать, как создать растр с разрешением 10 км х 10 км, который мог бы охватить все расширение территории Бразилии в R.

Answer 1

Вы можете сделать что-то по этим направлениям

Получить полигон для Бразилии

library(raster)
bra <- getData("GADM", country="BRA", level=1)

Спроецировать его, чтобы получить метры в единицах. Возможно, существует система координат, которая лучше подходит для ваших целей.

bcrs <- "+proj=poly +lat_0=0 +lon_0=-54 +x_0=5000000 +y_0=10000000 +ellps=GRS80"
pbra <- sp::spTransform(bra, bcrs)

Создание слоя RasterLayer

r <- raster(pbra, res=10000)

r
#class      : RasterLayer 
#dimensions : 432, 482, 208224  (nrow, ncol, ncell)
#resolution : 10000, 10000  (x, y)
#extent     : 2794657, 7614657, 6265592, 10585592  (xmin, xmax, ymin, ymax)
#crs        : +proj=poly +lat_0=0 +lon_0=-54 +x_0=5000000 +y_0=10000000 +ellps=GRS80 

Отображение

values(r) <- 1:ncell(r)
plot(r)
lines(pbra)