r에서 Esri Shape 파일을 읽고 플로팅합니다

Question

.Shp (esri shape file)에서 R을 읽는 데 어려움이 있습니다. R에서 몇 가지 옵션을 시도했으며 ArcMap의 모양 파일을 모양 파일에서 올바르게 읽는 것으로 변환하려고했지만 아직 작동하지 않았습니다. (ArcMap에서는 단일에서 멀티 폴리곤으로 변환 된 지오메트리를 수정했습니다.

아마도 내 모양 파일에 '다각형'대신 '영역'(다중 폴리 톤)이 포함되어 있다는 사실에 대한 것이있을 것입니다 ...

플로팅을 위해 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

1 월,

당신의 모양 파일은 아마도 올바르게 읽었을 것입니다. 음모가 올바르지 않았을 가능성이 높습니다. plot () 대신 spplot ()을 시도하십시오.

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

ShapeFile이 둘 이상의 열 (가능성이 있음)이있는 경우 ZCOL = 1 (또는 다른 열)을 추가하여 하나의 패널 만 얻습니다.

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

실제로 플롯을 사용하려면 배경에 대한 색상을 설정해야합니다. 그렇지 않으면 구멍이 투명하게 인쇄되므로. PBG = "화이트"를 설정하십시오. 노력하다:

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

그건 그렇고, Dirk가 언급했다고 생각하는 목록은 R-SIG-GEO 그리고 당신은 거기에서 많은 답변을 찾고 많은 SIG 관련 질문을 할 수 있습니다.