对具有给定分辨率（以公里为单位）的矢量进行光栅化

Question

我想知道如何在 R 中创建分辨率为 10km X 10km 的栅格，以覆盖整个巴西领土。

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)

创建光栅图层

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)

对具有给定分辨率（以公里为单位）的矢量进行光栅化

Raster over a vector with a given resolution in kilometers

r

r-raster