2015-02-06 4 views
0

그래서 모양 파일 (페이지 하단에 POWIATY)이 있습니다. 내가 가지고있는 좌표로 특정 다각형을 채우고 싶습니다. 위도 : 52.599427 경도 :R : ggplot2 코드에 따라 모양 파일의 다각형 채우기

이의 말을하자, 나는에 전체 영역을 채우려 20.7572137

내가지도에 점 좌표를 설정하지만 어떻게 거기에 전체 다각형을 채우는 방법을 알아? 당신이 폴리곤 숫자 ID로 기입해야합니다 있도록

+0

내가 질문을 이해할 수 없다, 당신은 더 명확하게 설명 할 수 있습니까? – Math

+0

색상을 지정하려는 지역의 모양 파일에 정의 된 그룹 또는 지역 (또는 다른 용어)이 있습니까? – lawyeR

답변

2

그 모양 파일에 중복 된 지역 이름이 있습니다

library(rgdal) 
library(rgeos) 
library(ggplot2) 

pow <- readOGR("POWIATY.shp", "POWIATY") 
plot(pow) 

where <- over(SpatialPoints(cbind(20.7572137, 52.599427)), pow, TRUE) 
reg <- data.frame(id=rownames(where[[1]])) 

map <- fortify(pow) 

gg <- ggplot() 
gg <- gg + geom_map(map=map, data=map, 
        aes(x=long, y=lat, map_id=id), 
        fill="white", color="black", size=0.25) 
gg <- gg + geom_map(data=reg, map=map, 
        aes(fill=id, map_id=id), color="steelblue") 
gg <- gg + coord_map() 
gg <- gg + theme_bw() 
gg 

enter image description here

+0

이것은 그 것이다! 고마워요. – kodi1911

+0

나는 잠시 동안 그것을 시험했다. 이것을 다른 shapefile에 추가하려고하면 overfunction을 사용하여 "변수 identicalCRS (x, y)가 TRUE가 아닙니다"라고 표시됩니다. 또한지도의 모든 다각형에 대해 gps를 사용하여 채울 때 일부 공백이 표시됩니다. shapefile에 대한 링크가 있습니다. [link] (ftp://91.223.135.109/prng/PRNG_MIEJSCOWOSCI_SHP.zip) – kodi1911