4

ここからニューオーリンズの郵便番号マップに色を付けようとして、このチュートリアルの「簡単なマップ」の手順に従っています(そのリンクからニューオーリンズの 2011 データの .shp ファイルを使用しています)。

チュートリアルのようにファイルをロードしようとすると、以下のエラーが発生します。

nolazip.shp <- readShapePoly("/PathTo/Orleans_ZCTA_2010_SP.shp", proj4string=CRS("+proj=longlat"))
Error in validityMethod(as(object, superClass)) : 
  Geographical CRS given to non-conformant data: 3820725.379655  613426.584024

このドキュメントに基づくと、このエラーは、シェープ ファイルが有効な longlat データを持つ proj4string を使用していないことを意味しているようです。

他の種類の proj4string または CRS オブジェクトを使用していますか?

これらのコマンドを実行して、CRS の出力を検索して見つけようとしましたが、何も見つかりませんでした。

    > summary(orcounty.shp)
    > str(orcounty.shp)

readShapePoly コマンドで proj4string パラメータを省略するだけで形状ファイルをインポートできますが、「いくつかの単純なマップ」セクション (唯一のセクション私は欲しい）。

  1. シェープファイルに関連付けられている proj4 文字列は何ですか? それを readShapePoly への入力として与えるにはどうすればよいですか
  2. このマップ作成方法で動作するシェープファイルをインポートする別の方法はありますか? 繰り返しますが、単に問題のあるパラメーターを除外すると、R スタジオのプロットにマップが表示されません。
4

2 に答える 2

6

を使用してこれに取り組みますreadOGR。これは投影情報を保持するため、上記の質問のようにいじる必要はありません。これは、同じシェープファイル (この米国政府のサイトからダウンロード) を読み込んで にプロットしたものと思われるものですggplot2。化粧品はおそらく片付ける必要がありますが、それでRColorBrewerスケールやその他のggplot2ことを練習できます. [編集 - に欠落しているaes呼び出しを追加geom_polygon]

# if the packages are not installed, you will have to install and load them.
install.packages("rgdal")
install.packages("ggplot2")
install.packages("scales")
library(rgdal)
library(ggplot2)
library(scales)

require(rgdal)
require(ggplot2)
require(scales)

work.dir <- "your_dirname" # your directory
                                     # no trailing slash

orl <- readOGR(work.dir, layer = "Orleans_ZCTA_2010_SP")
orl.df <- fortify(orl) # ggplot needs data frame, not spatial object

ggplot(data = orl.df, aes(x = long, y = lat, group = group)) +
    geom_polygon(aes(fill = orl.df$group)) +
    coord_equal() +
    theme(legend.position = "none")

オルレアン

于 2013-10-26T23:51:15.813 に答える