0

Rで再生しているGISは、ggplot2やその他のパッケージを使用して緯度経度座標を使用してプロットすることであったため、これを開始する方法がわかりません。ここで行う必要があるのは、プロットのxy座標を含むHTMLファイルの形式のシェープファイルを使用するMicrostrategyの視覚化コンポーネントを使用することです（つまり、左上は0,0です）。状態レベルのファイルの例は次のとおりです。

<HTML><HEAD><TITLE>untitled</TITLE></HEAD><BODY>
<IMG SRC="" USEMAP="#myMap" WIDTH="812" HEIGHT="713" BORDER="0" />
<MAP NAME="myMap">
<AREA SHAPE="POLY" HREF="#" ALT="Texas" COORDS="299,363,299,360,....." />
</MAP></BODY></HTML>

「座標」にリストされているポイントは、その場でプロットおよび色付けされた812x713の「イメージ」に関するXポイントとYポイントです。

http://www.vdstech.com/usa-data.aspxからZip3およびZip5のshp、shx、dbfファイルがありますが、変換を開始する場所がわかりません。HTMLファイルを手作業でフォーマットするといううんざりする作業をしてもかまいません。私が立ち往生しているのはXY変換です（さびた、しばらくGISに触れていない）：

次のコードは、シェープファイルをRにインポートします

library(rgdal)
zip3 <- readOGR(dsn = '/Users/adempsey/Downloads/zip3'), layer = 'zip3')

その後、私は立ち往生していて、zip3 + xy座標をデータフレームに抽出する方法のチュートリアルを探しています。これを使用して、最終的なファイルを作成できます。

アップデート2

以下を使用して、データフレームに変換できますが、関連付けられたdbfファイルに保存されているように見える関連付けられたzip3コードをプルできません

Row long      lat order  hole piece group id
1 -151.0604 70.41873     1 FALSE     1   0.1  0
2 -150.7620 70.49722     2 FALSE     1   0.1  0

はい、これは私の現在のさびたRを超えています

update3このコードは、郵便番号をデータフレームにダンプします。zip3.codes<-as.data.frame（zip3）zip3.df <-fortify（zip3 @ polygons [[1000]]）のようなものと組み合わせることができます。特定のzip3に関連付けられているすべての行zip3.codesに置き換えられます

4

1 に答える 1

1

fastshpパッケージを使用してデータをロードできます。

 install.packages("fastshp",,"http://rforge.net")
 library(fastshp)
 s <- read.shp("zip5.shp", format="polygon")

sすべての ZIP 形状のリストになりました。xおよびコンポーネントに興味があります。yたとえば、最初の ZIP をプロットするには、単に次のようなものを使用します

 plot(s[[1]]$x, s[[1]]$y, asp=1.25)
 polygon(s[[1]]$x, s[[1]]$y, col="#eeeeee")

名前を一致させるには、 from を使用read.dbfしforeignます。

 library(foreign)
 d <- read.dbf("zip5.dbf", as.is=TRUE)
 names(s) <- d$ZIP5

?read.shp使用可能な形式の詳細については、を参照してください。NA個々のポリゴンを分離するために使用する「ポリゴン」 、「リスト」はインデックスを使用してパーツを提供します。

ところで、データセットはやや疑わしいので、TIGER/Line census ZCTA5 データ (最新は 2010 年) を調べることをお勧めします。

于 2013-02-14T04:07:19.640 に答える