【R】北海道の農業集落境界

G空間情報センターからダウンロードできる北海道の農業集落境界を表示してみます。

library(ggplot2)
library(sf)

shapefile_url = "http://www.machimura.maff.go.jp/shurakudata/2015/ma/MA0001_2015_2015_01.zip"

zip_file <- "zipfile.zip"
download.file(shapefile_url, destfile = zip_file) 
unzip(zip_file)          

sf_pref01_field <- read_sf("MA0001_2015_2015_01/rcom.shp") # 北海道の農業集落境界シェープファイル

ggplot() + 
  geom_sf(data = sf_pref01_field)

なかなかファイルが処理が多く、時間がかかります。。。

Add a Comment

メールアドレスが公開されることはありません。