내 목표는 특정 다각형을 병합하여 기존 모양 파일을 수정하는 것입니다.
shapefile을 가져오고 UnionSpatialPolygons 명령을 사용한 후 원하는 다각형 외곽선을 얻습니다.
그러나 이것은 이제 SpatialPolygonsDataFrame이 아니라 SpatialPolygons 객체이므로 writeOGR을 사용하여 shapefile로 내보낼 수 없습니다.
이 문제를 어떻게 해결할 수 있습니까?
답변:
이름에서 알 수 있듯이 SpatialPolygonsDataFrame은 기본적으로 데이터가 첨부 된 SpatialPolygons 객체 (속성 테이블)입니다. 데이터에는 기능이있는 행 수 이상이 있어야합니다.
library(rgdal)
ob <- SpatialPolygons(..)# Your SpatialPolygons Object
spp <- SpatialPolygonsDataFrame(ob,data=as.data.frame("yourData"),proj4string=CRS("+proj= aea > +ellps=GRS80 +datum=WGS84"))
writeOGR(spp,"shapes","testShape",driver="ESRI Shapefile",)----편집하다----
SpatialPolygonsDataFrame등을 SpatialPolygons객체 로 변환하려면 R 내의 객체 구조를 지정하면됩니다.
ob <- SpatialPolygons(spp@polygons,proj4string=spp@proj4string)IDs실행할 때 인수를 올바르게 지정하십시오 UnionSpatialPolygons. 그렇지 않으면 data로 다시 변환 할 때 인수를 지정하는 데 문제가 발생할 수 있습니다 SpatialPolygonsDataFrame.
data=as.data.frame("yourData")무엇입니까? 사용하여 내부 다각형을 용해 후 unionSpatialPolygons(...), 나는 새로운 Shape 파일로 결과를 작성하려합니다
union.
문제 :
1 : UnionSpatialPolygons의 결과는 공간 다각형입니다
2 : 결과를 다시 공간 폴리곤 데이터 프레임으로 변환하는 것은 큰 고통입니다
-ㅏ. 공간 다각형에 부착하려면 매우 정확한 데이터 프레임이 필요합니다
-비. UnionSpatialPolygons에 사용한 데이터는 출력보다 많은 행을 가지며 필요한 방식으로 형식이 지정되지 않습니다.
내 (추악한) 해결책 :
### Coerce into spatial polygon data frame with id and row name of spatial polygon
# Make a data frame that meets the requirements above:
df<- data.frame(id = getSpPPolygonsIDSlots(your.spatialpolygon))
row.names(df) <- getSpPPolygonsIDSlots(your.spatialpolygon)
# Make spatial polygon data frame
spdf <- SpatialPolygonsDataFrame(your.spatialpolygon, data =df)
# Then don't forget to make sure the projection is correct
# XXXX is your SRID
proj4string(spdf) <- CRS("+init=epsg:XXXX");
spdf <- spTransform(spdf , CRS("+init=epsg:XXXX"));