使用pbsmapping和shapefiles帮助在r中绘制地理数据

使用 O'Reilly's Data Mashups in R 犹他州盐湖县的一个形状文件中发现了一些地址,作为灵感来源,我正试图将其绘制出来。 here .

我有数据帧地质表:

> geoTable
         address        Y         X EID
1    130 E 300 S 40.76271 -111.8872   1
2    875 E 900 S 40.74992 -111.8660   2
3   2200 S 700 E 40.72298 -111.8714   3
4    702 E 100 S 40.76705 -111.8707   4
5 177 East 200 S 40.76518 -111.8859   5
6    702 3rd ave 40.77264 -111.8683   6
7   2175 S 900 E 40.72372 -111.8652   7
8   803 E 2100 S 40.72556 -111.8680   8

我已经将它强制为一个EventData对象:

> addressEvents<-as.EventData(geoTable,projection=NA)
> addressEvents
         address        Y         X EID
1    130 E 300 S 40.76271 -111.8872   1
2    875 E 900 S 40.74992 -111.8660   2
3   2200 S 700 E 40.72298 -111.8714   3
4    702 E 100 S 40.76705 -111.8707   4
5 177 East 200 S 40.76518 -111.8859   5
6    702 3rd ave 40.77264 -111.8683   6
7   2175 S 900 E 40.72372 -111.8652   7
8   803 E 2100 S 40.72556 -111.8680   8

所以看起来我已经准备好了所有需要设计的东西,但它不起作用。当我加载shapefile并使用

addPoints(addressEvents,col="red",cex=.5)

我在看一个空的形状文件。此外,当我尝试对我的EventData对象运行findpolys时,它返回空值。

> findPolys(addressEvents,myShapeFile)
NULL

我怎样才能做到这一点?我能够毫无问题地完成O'Reilly教程,并且很难确定我在这里哪里出错。我不知道它是形状文件,我的数据框架,还是其他什么。

以下是用于导入数据和形状文件的命令

slc<-read.table('~/utah.txt',sep=',',header=TRUE,strip.white=TRUE,stringsAsFactors=FALSE)

myShapeFile<-importShapefile("/Users/neil/Downloads/SGID93_DEMOGRAPHIC_CensusTracts2000/SGID93_DEMOGRAPHIC_CensusTracts2000",readDBF=TRUE)