2

まず、タイトルについてお詫び申し上げますが、マップで何が起こっているかを説明する方法が完全にはわかりません (したがって、画像が含まれています)。私は ggplot2 でマッピングするのが初めてで、最初の「チュートリアル外」マップで、DC の地域ごとに国勢調査データを表示しようとしています。チュートリアルを検索すると、私が行った手順がサポートされているように見えますが、明らかにコースから外れています (このような例は他に見たことがありません)。

目標は、Census tract polys を含むシェープファイルで準備した .csv からのデータをマージすることでした。最初は、行の順序に関係があると思っていましたが、これを明示的に説明しても何も変わっていないようです。

コード：

###Combining Census data with a tract poly shapefile
library(maptools)
library(ggplot2)
library(gpclib)

setwd('~/ESRI/Trend Report Maps')

#Read data
tract<-readShapePoly('TractPly.shp')
cdata<-read.csv('census00_10.csv')

#Columns >> note that GEOID (tract) and geoid2 (cdata) are the merge keys
ntract<-names(tract)
ncdata<-names(cdata) 


#Prepare data for ggplot plotting of poly info (by conversion to DataFrame)
gpclibPermit()
tract_geom<-fortify(tract,region="GEOID")
#Note that this drops attribute and retains only spatial info. However, we don't really
#need the attribute info since we are joining it to the other dataframe

#Merge
tract_poly<-merge(tract_geom,cdata,by.x="id",by.y="geoid2")
tract_poly<-tract_poly[order(tract_poly$order),]
head(tract_poly)

workF<-ggplot(tract_poly,aes(long,lat,group-group,fill=dmed_age)) + 
        geom_polygon() +
        coord_equal() #fixes aspect ratio
workF

ここで何が起こっているのですか??

ええ、国勢調査データの視覚化は、より良い時代を迎えています。地図を表示することを許可されていない場合、ギザギザの線を無計画に並べたように見えます (誰かが紙を何度も折りたたんで、広げていく前に紙を無作為に切り取ったときと同じように見えます)。DCの外側の境界は無傷のままですが、結果として得られる形状はどれも国勢調査区のようには見えません。シェープファイルが ArcMap からエクスポートされたことを確認できます (本来あるべき場所)。

セッション情報:

> sessionInfo()
R version 2.15.3 (2013-03-01)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] grid      stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] gpclib_1.5-1    ggplot2_0.9.3.1 maptools_0.8-23 lattice_0.20-13 sp_1.0-5        foreign_0.8-52 

loaded via a namespace (and not attached):
 [1] colorspace_1.2-1   dichromat_2.0-0    digest_0.6.3       gtable_0.1.2       labeling_0.1      
 [6] MASS_7.3-23        munsell_0.4        plyr_1.8           proto_0.3-10       RColorBrewer_1.0-5
[11] reshape2_1.2.2     scales_0.2.3       stringr_0.6.2      tools_2.15.3
4

1 に答える 1

1

ただのタイプミスです。に変更group-groupするgroup=groupと、問題ないはずです。

于 2013-04-26T12:47:05.010 に答える