ポイント データのシェープ ファイルから始めました。次に、この形式を ppp 形式に変換して、いくつかの spatstat 関数を使用しました。
これを行うコードを以下に示します。
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 形式に変換するにはどうすればよいですか?