Преобразование изображения в геотип с использованием R

Question

Я начал с файла формы данных точек. Затем я преобразовал этот формат в формат ppp, чтобы использовать некоторые функции spartstat.

код, чтобы сделать это приведен ниже:

sp_points_df_al <- rgdal::readOGR(dsn = shape_path, layer = "aland_points") 
# convert the spatial points data.frame to just plain ol spatial points 
sp_points_al <- as(sp_points_df_al, "SpatialPoints") 
# now convert the spatial points to ppp 
ppp_al <- as(sp_points_al, "ppp") 
# estimate the window from the points data 
wind_rr <- spatstat::ripras(ppp_al) 
ppp_al_constr <- spatstat::ppp(ppp_al$x, ppp_al$y, window = wind_rr) 

Затем я захотел сделать некоторый анализ точечных шаблонов и так использовал следующий spatstat функции:

image <- spatstat::density.ppp(ppp_al_constr, sigma = 0.004, dimyx=c(512, 512)) 

Я тогда отображение изображения используя следующее:

bias_palette <- colorRampPalette(c("blue", "magenta", "red", "yellow", "white"), bias=2, space="Lab") 
spatstat::plot.im(k, col=bias_palette(256), ribbon = FALSE) 

Что я хочу сделать дальше, это создать изображение Geotiff из выше. Я хочу сделать это, поскольку хочу наложить изображение tiff поверх некоторых векторных данных.

Мой вопрос:

Как преобразовать изображение выше в формат Geotiff ...?

Answer 1

Результат spatstat::density.ppp является объектом класса im. Если вы делаете as.matrix этого объекта, вы получаете матрицу, содержащую записи пикселей. Затем вы можете преобразовать это в любой другой тип изображения. Также as.data.frame(Frame(..)) объекта изображения будет извлекать ограничительную рамку.

Или перейти к maptools пакет, который обрабатывает spatstat объекты, и искать функцию, чтобы преобразовать spatstat класс im объект нужного класса.