【R】世界の人口密度

世界の人口密度を表示します。データは、NASAのEarthDATAから取得します（ログインが必要です。）。

こちらのページを参考にさせていただきました。

library(raster)
library(rayshader)

pop = raster_to_matrix(raster("./gpw_v4_population_density_adjusted_to_2015_unwpp_country_totals_rev11_2020_30_min.tif"))

pal=colorRampPalette(c("red", "white"))

pop %>%
  height_shade(texture=pal(256))%>%
  plot_3d(pop, zscale = 100, solid = FALSE,
          shadowdepth = 0.3, windowsize = c(1400,700), theta = 0,zoom=0.6)

予想通り、中国、インドなどアジアは高いですね。東京がとびぬけています（笑）。Rで実行するとインタラクティブな地図になります。