27

米国の 50 州すべてを示すテーマ マップを作成しようとしていますが、信頼できる方法でアラスカとハワイを再配置するのに問題があります。いくつかのアイデアがありますが、どれもうまくいきません。私は今それらを実演します。

まず、データをインポートする必要があります。mapsハワイとアラスカが含まれていないため、パッケージ内のデータを使用するだけでは不十分です。

setwd(tempdir())
download.file("https://dl.dropbox.com/s/wl0z5rpygtowqbf/states_21basic.zip?dl=1", 
              "usmapdata.zip", 
              method = "curl")
# This is a mirror of http://www.arcgis.com/home/item.html?
# id=f7f805eb65eb4ab787a0a3e1116ca7e5
unzip("usmapdata.zip")

require(rgdal)
all_states <- readOGR("states_21basic/", "states")

require(ggplot2); require(maptools); require(rgeos); require(mapproj);
all_states <- fortify(all_states, region = "STATE_NAME")

ここで、いくつかのプロットの美学を定義します。

p <- ggplot() + geom_polygon( 
  aes(x=long, y=lat, group = group, fill = as.numeric(as.factor(id))), 
  colour="white", size = 0.25
) + coord_map(projection="azequalarea") + 
scale_fill_gradient(limits = c(1,50))

次に、すべての背景などを削除して、隣接していない状態を重ねたときに衝突しないようにします。

p <-   p + theme(axis.line=element_blank(),
            axis.text.x=element_blank(),
            axis.text.y=element_blank(),
            axis.ticks=element_blank(),
            axis.title.x=element_blank(),
            axis.title.y=element_blank(),
            panel.background=element_blank(),
            panel.border=element_blank(),
            panel.grid.major=element_blank(),
            panel.grid.minor=element_blank(),
            plot.background=element_blank())

ビューポートの使用

私が最初に思いついたのは、ビューポートを使用することでした:

AK <- p %+% subset(all_states, id == "Alaska") + theme(legend.position = "none")
HI <- p %+% subset(all_states, id == "Hawaii") + theme(legend.position = "none")
contiguous <- p %+% subset(all_states, id != "Alaska" & id != "Hawaii")

grid.newpage()
vp <- viewport(width = 1, height = 1)
print(contiguous, vp = vp)
subvp1 <- viewport(width = 0.25, height = 0.25, x = 0.18, y = 0.33)
print(AK, vp = subvp1)
subvp2 <- viewport(width = 0.12, height = 0.12, x = 0.32, y = 0.27)
print(HI, vp = subvp2)

ここに画像の説明を入力

これは見栄えは良いですが、サイズ変更や凡例のサイズと形状の変更など、Figure のわずかな変更に非常に敏感であるため、満足のいくものではありません。

アラスカとハワイを手動で移動する

all_states_AKHImoved <- within(all_states, {
  lat[id == "Alaska"] <- lat[id == "Alaska"] - 45
  long[id == "Alaska"] <- long[id == "Alaska"] + 40
  lat[id == "Hawaii"] <- lat[id == "Hawaii"] + 0
  long[id == "Hawaii"] <- long[id == "Hawaii"] + 70
})
p %+% all_states_AKHImoved

ここに画像の説明を入力

これは満足できるものではありません。通常、アラスカはほとんどの米国のマップで縮尺調整されていないため、非常に大きく見えるからです。また、アラスカとハワイを再配置すると、地図投影によって導入された歪みが変化します。

質問

誰もがより良いアプローチを持っていますか?

4

3 に答える 3

14

投影して変換する方法は次のとおりです。必要になるだろう：

require(maptools)
require(rgdal)

fixup <- function(usa,alaskaFix,hawaiiFix){

  alaska=usa[usa$STATE_NAME=="Alaska",]
  alaska = fix1(alaska,alaskaFix)
  proj4string(alaska) <- proj4string(usa)

  hawaii = usa[usa$STATE_NAME=="Hawaii",]
  hawaii = fix1(hawaii,hawaiiFix)
  proj4string(hawaii) <- proj4string(usa)

  usa = usa[! usa$STATE_NAME %in% c("Alaska","Hawaii"),]
  usa = rbind(usa,alaska,hawaii)

  return(usa)

}

fix1 <- function(object,params){
  r=params[1];scale=params[2];shift=params[3:4]
  object = elide(object,rotate=r)
  size = max(apply(bbox(object),1,diff))/scale
  object = elide(object,scale=size)
  object = elide(object,shift=shift)
  object
}

次に、シェープファイルを読み込みます。使用rgdal:

us = readOGR(dsn = "states_21basic",layer="states")

ここで等面積に変換し、fixup 関数を実行します。

usAEA = spTransform(us,CRS("+init=epsg:2163"))
usfix = fixup(usAEA,c(-35,1.5,-2800000,-2600000),c(-35,1,6800000,-1600000))
plot(usfix)

パラメータは、それぞれアラスカとハワイの回転、スケーリング、x および y シフトであり、試行錯誤によって取得されました。慎重に微調整してください。ハワイの縮尺パラメータを 0.99999 に変更しても、含まれる数が多いため、地球から離れてしまいました。

これを緯度経度に戻したい場合:

usfixLL = spTransform(usfix,CRS("+init=epsg:4326"))
plot(usfixLL)

ggplotしかし、 で変換を行ったので、 で変換を使用する必要があるかどうかはわかりませんspTransform。

ggplot2要塞化ビジネスをジャンプできるようになりました。あなたにとって重要かどうかはわかりませんが、usfixバージョンでは州の順序が異なっていることに注意してください。アラスカとハワイが最後の 2 つの州になりました。

于 2012-12-07T17:23:17.567 に答える
3

このように物事をシフトし始めると、アラスカとハワイをメキシコ湾のどこかにある正方形のボックスとして表すだけでよいでしょう. ハワイが何色であるかを知ることができるという追加の利点があります.

たぶん、どこまでも行って、各州の面積が等しい歪んだシステムを使用すれば、ロードアイランドを見ることができます。

グーグル画像の米国カートグラムの例は、そのようなものを示しています。ただし、シェープファイルまたはデータが含まれているものがいくつあるかはわかりません。

状態の相対的なサイズに本当に関心がありますか、それとも状態の値が何であるかを人々に示す表現が必要ですか?

于 2012-12-07T08:32:27.303 に答える