R: ggplot2 заполняет многоугольники в шейп-файле с помощью коордов

Question

Итак, у меня есть файл формы (с именем POWIATY внизу). Я хочу заполнить конкретные полигоны с помощью коордов, которые у меня есть.

Скажем, я хочу, чтобы заполнить весь регион: Lat: 52,599427 долготу: 20,7572137

Я знаю, как установить точку COORDS на карте, но как заполнить весь полигон в нем?

Answer 1

Есть дублирующие названия региона в этом шейп, так что вы должны будете заполнить полигоном числовым идентификатором:

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