R - 栅格化 SpatialPolygonsDataFrame - 变量/属性被删除答案

【问题标题】：R - Rasterize a SpatialPolygonsDataFrame - variables/attributes droppedR - 栅格化 SpatialPolygonsDataFrame - 变量/属性被删除
【发布时间】：2016-11-12 14:06:10
【问题描述】：

每当我栅格化我的 SpatialPolygonsDataFrame 时，我都会丢失属性/数据信息部分。命令“rasterize”来自R中的“raster”包。

我有以下 RasterLayer（名为“raster1”）

class       : RasterLayer 
dimensions  : 6000, 4800, 28800000  (nrow, ncol, ncell)
resolution  : 0.00025, 0.00015  (x, y)
extent      : 8.699875, 9.899875, 46.69993, 47.59993  (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:3857 +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs

我的 SpatialPolygonsDataFrame（名为“bw1”）具有以下属性

class       : SpatialPolygonsDataFrame 
features    : 7663 
extent      : 980075.6, 1076577, 5908811, 6023151  (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:3857 +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs 
variables   : 2
names       :    A,  B 
min values  : 3231, 11 
max values  : 3955, 19

当我使用命令时

bw1_raster<-rasterize(bw1,raster1,fun='last',field=c("A","B")

我得到了新对象“bw1_raster”：

class       : RasterLayer 
dimensions  : 6000, 4800, 28800000  (nrow, ncol, ncell)
resolution  : 0.00025, 0.00015  (x, y)
extent      : 8.699875, 9.899875, 46.69993, 47.59993  (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:3857 +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs 
data source : in memory
names       : layer     
values      : NA, NA  (min, max)

如何获得具有“A”和“B”属性/图层的栅格对象？

【问题讨论】：

标签： r spatial r-raster

【解决方案1】：

问题是 SpatialPolygons 的坐标参考系统与 RasterLayer 的坐标参考系统不匹配。 看起来它们匹配，但这可能是因为您将 RasterLayer 的 crs 更改为不正确的值。 RasterLayer 几乎可以肯定有一个 longlat crs，而不是 merc。

您需要为其分配正确的 crs（或者至少不要将其更改为错误的 crs！）。

根据您想要实现的目标，您可以使用 spTransform 将 SpatialPolygons 转换为 RasterLayer 的 crs，然后重试。

或者，您可以按照这些方式做一些事情

library(raster)
r <- raster(bw1, res=10000)
r <- rasterize(bw1, r)

【讨论】：