Как сделать фрейм данных в растровом объекте R

Question

У меня есть фрейм данных с одним столбцом значений.Как превратить это в растровое изображение со следующими параметрами?Я также хотел бы сохранить его как растровое изображение .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 + unit = 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))

Answer 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)