读取和绘图中的R ESRI形状文件

Question

我有我试图R中的几个选项，并试图形状文件转换ArcMap中的东西，在造型上的文件，但工作没有正确地读取困难的.SHP（ESRI形状文件）读入R.然而。 （在ArcMap我纠正的几何形状，从单个转换为MULTIPOLYGON等这可能是不必要的或不相关的）

这可能与我的形状文件中包含“区”（多多边形）的事实的东西，而不是“多边形” ...

如何正确读出中的R为绘制该类型形文件？ （它看起来像在ArcMap正常形状）

在ArcMap中的形状文件看起来像这样： http://bfast.r-forge.r-project.org/arcmapshape.jpg

（示出了具有其他多边形内的多边形的形状文件）

在该R它看起来像这样：

（示出了一个形状文件，其中一些多边形被错误填充的）

我用下面的代码中R：

require(maptools)
require(rgdal)

newproj <- "+proj=utm +zone=55 +south +ellps=GRS80 +units=m"
shape<- readShapeSpatial(pdir, proj4string = CRS(newproj),repair=TRUE,force_ring=T,verbose=TRUE) # without any errors
plot(shape, col="gray",border="blue", axes=TRUE)

# via rgdal
folder <- c("spatial")
lyr <- c("clipped_forest_mga")
shp <- readOGR(dsn=folder,layer=lyr)
plot(shp, col="gray",border="blue", axes=TRUE)

这两个曲线图（）命令给了相同的R-结果。没有发生错误。只有下面的消息

OGR data source with driver: ESRI Shapefile 
Source: "P:/Victoria_DSE/BFAST_spatial/vector/PLM_excl_fire03_09_GDA94LL/mgaz94z55/clipped_EG", layer: "clipped_forest_mga"
with 1 features and 4 fields
Feature type: wkbMultiPolygon with 2 dimensions

这又如何解决？

Answer 1

一月，

您shape文件可能已被正确读取。它更可能是绘图是不正确的。尝试spplot（）而不是图（）：

spplot(shp, col.regions="gray", col="blue")

如果您的shape文件有一个以上的列（这是可能的），添加ZCOL = 1（或任何其它列）这样就可以获得只有一个面板

spplot(shp, zcol=1, col.regions="gray", col="blue")

如果你真的想使用你必须设置背景颜色的情节，因为否则将打印透明的孔。因此，设置PBG =“白色”。尝试：

plot(shp, col="gray", border="blue", axes=TRUE, pbg="white")

顺便说，该列表我认为德克指的是 R-SIG -GeO ，你可以找到很多有答案，问很多SIG相关的问题。