esri 모양의 파일 읽기 및 플로팅 R

Question

R에 .shp (esri 모양 파일)을 읽는 데 어려움을 겪고 있습니다. R에서 여러 옵션을 시도하고 ArcMap의 모양 파일을 셰이프 파일을 올바르게 읽었지만 아무 것도 작동하지 않았습니다. (ArcMap에서 형상을 수정하고 단일 또는 다중 다각형 등을 아마도 필요하지 않거나 관련성이없는 것으로 바꿨습니다)

아마도 내 모양 파일에 'regions' 'polygons'...

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

두 가지 plot() 명령은 동일한 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 월,

모양 파일이 올바르게 읽혔습니다. 작의가 올바르지 않은 경우가 많습니다. spplot() 대신 플롯()을보십시오 :

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

당신의 Shape 파일은 (가능성이) 하나 이상의 열이있는 경우, zcol = 1 (또는 다른 열) 그래서 당신이 얻을 하나 개의 패널을 추가

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

플롯을 실제로 사용하려면 배경색을 설정해야합니다. 그렇지 않으면 투명하게 구멍이 인쇄되기 때문입니다. 따라서 pbg = "white"로 설정하십시오. 나는 더크는 R-sig-Geo 언급 된 생각, 그런데 목록을

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

그리고 당신이 많은 해답을 찾을 수 많은 시그마 관련 질문을 할 수 있습니다보십시오.