如何将数据框制作成栅格对象R答案

【问题标题】：How to make data frame into raster object R如何将数据框制作成栅格对象R
【发布时间】：2019-11-12 04:52:07
【问题描述】：

我有一个包含一列值的数据框。如何使用以下参数将其制作成光栅图像？我还想将其保存为 .tif 光栅图像。

尺寸：(6,7,42) (nrow, ncol, ncell)

分辨率 (30,30) (x,y)

范围：(286695, 286905, 4620945, 4621125) (xmin, xmax, ymin, ymax)

crs: +proj=utm +zone=17 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0

值：0、65535（最小值、最大值）

感谢您的帮助！

library(raster)

#making dataframe
b1DF <- data.frame("b1" = c(18880, 19364, 20446, 21070, 23001, 22981, 23362,
 19538,18902, 19013, 18676, 20998, 23206, 24784, 21979, 20804, 20496, 19993,
 20830, 22396, 24246, 19877, 21002, 21771, 21478, 20827, 21267, 23036, 20103,
 21220, 22308, 21443, 20780, 21874, 23029, 19521, 20173, 20826, 21000, 21693,
 23388, 23935))

【问题讨论】：

标签： r r-raster

【解决方案1】：

与光栅包附带的大多数示例非常相似。

b1DF <- data.frame(b1 = c(18880, 19364, 20446, 21070, 23001, 22981, 23362,
 19538,18902, 19013, 18676, 20998, 23206, 24784, 21979, 20804, 20496, 19993,
 20830, 22396, 24246, 19877, 21002, 21771, 21478, 20827, 21267, 23036, 20103,
 21220, 22308, 21443, 20780, 21874, 23029, 19521, 20173, 20826, 21000, 21693,
 23388, 23935))

library(raster)
r <- raster(nrow=6, ncol=7, ext=extent(286695, 286905, 4620945, 4621125), crs="+proj=utm +zone=17 +datum=WGS84 +units=m")
values(r) <- b1DF$b1

r
#class      : RasterLayer 
#dimensions : 6, 7, 42  (nrow, ncol, ncell)
#resolution : 30, 30  (x, y)
#extent     : 286695, 286905, 4620945, 4621125  (xmin, xmax, ymin, ymax)
#crs        : +proj=utm +zone=17 +datum=WGS84 +units=m +ellps=WGS84 +towgs84=0,0,0 
#source     : memory
#names      : layer 
#values     : 18676, 24784  (min, max)

【讨论】：