Rでネパールのdistrict別人口密度の塗分け地図を作る

およそ1年ぶりにRの記事を書く。

今回はネパールの郡（district）別の人口密度の塗分け地図を作ってみた。

スクリプトは次の通り：

#おまじない rm(list=ls(all=TRUE)) # # library(maptools) par(family="sans") # x <- readShapePoly("c:\\rtemp\\NPL_adm\\NPL_adm3.shp") #シェープファイルの読み込み # ix <- length(x) # データのサイズ（行＝districtの数） zdata <- x$POP_2011/x$AREA # 人口密度の計算；POP_2011はdistrict別人口，AREAはdistrict面積 colnam <- rep("white", length=ix) # 塗分けのための変数 t <- c(50, 100, 500, 1000) # 塗分けのための閾値 color.set <- c("lemonchiffon", "yellow", "orange", "red", "darkred") # カラーパレット(heat color) unit <- "Persons/km^2" # zdataの単位 # # zdata（人口密度）に応じて，カラーパレットから色を選んで # 塗分けのための変数colnamに代入する for (ii in 1:ix) { if (zdata[ii] > t[4]) { colnam[ii] <- color.set[5] } else if (zdata[ii] > t[3]) { colnam[ii] <- color.set[4] } else if (zdata[ii] > t[2]) { colnam[ii] <- color.set[3] } else if (zdata[ii] > t[1]) { colnam[ii] <- color.set[2] } else { colnam[ii] <- color.set[1] } } # plot(x, col=colnam) # 人口密度で色分け # # 以下は凡例の描画用 oo <- par("usr") # プロット範囲（経度・緯度）をooに保存 ox <- oo[1] # 経度の値の小さい方を原点x座標に oy <- (oo[3] + oo[4])/2 # 緯度の値の小さい方と大きいほうの中間を原点y座標に oxl <- oo[2] - oo[1] # oyl <- oo[4] - oo[3] # par(new = T) # 現在の図に新たな図を重ね書きする for (ii in 1:5) { aa <- c(ox + oxl * 15 / 20, ox + oxl * 15 / 20, ox + oxl * 16 / 20, ox + oxl * 16 / 20) bb <- c(oy + oyl / 20 * ii, oy + oyl / 20 * (ii + 1), oy + oyl / 20 * (ii + 1), oy + oyl / 20 * ii) polygon(aa, bb, col=color.set[ii]) if (ii < 5) text(ox + oxl * 16.5 / 20, oy + oyl / 20 * (ii + 1), labels = t[ii], adj = 0) else text(ox + oxl * 15 / 20, oy + oyl / 20 * (ii + 1.5), labels = unit, adj = 0) }

詳細はここに記載。

出来上がりはこれ：

首都圏（カトマンドゥ郡，バクタプル郡，ラリトプル郡）の人口密度が非常に高いことがわかる。カトマンドゥなんか，1平方キロメートル当たり4500人以上である。

この地図では75郡だが，実は，2015年の憲法で行政区分がいろいろと変わり，現在ではナワルパラシ郡とルクム郡は東西に分割されて77郡となっている。