ggplot2

Question

का उपयोग करके एकाधिक उप-विभाजन और छेद के साथ स्पैटलियल पॉलीगॉन ड्रा करें I ggplot2 में छेद के साथ लाइब्रेरी sp से SpatialPolygons आकर्षित करना चाहते हैं। stackoverflow पर दूसरों सवाल के लिए धन्यवाद, मैं जानता हूँ कि इस घड़ी की लिखित बहुभुज साथ काम करते हुए अनुमति दी है कि: के रूप में तैयार हो
http://stackoverflow.com/questions/12047643/geom-polygon-with-multiple-hole/12051278#12051278
वास्तव में, जब एक SpatialPolygons बदलने broom::tidy का उपयोग कर (ggplot2::fortify की जगह), छेद बहुभुज घड़ी की दिशा में सहेजे जाते हैं छेद।
ggplot2 में, जिस तरह से छेद वाले बहुभुज उन्हें fill का उपयोग करके एक बार खींचने के लिए मजबूर कर रहे हैं, और दूसरी बार colour का उपयोग करके, अन्यथा आप बहुभुज को पार करने वाली रेखाएं देख सकते हैं। एकाधिक सबपोलिगन्स से निपटने पर, कुछ छेद के साथ, यह अधिक कठिन है, broom::tidy द्वारा परिभाषित बिंदु सुविधाओं का क्रम बहुभुज भरने की अनुमति नहीं दे सकता है (नीचे छवि देखें)।
क्या आप में से किसी को भी इस समस्या को हल करने के लिए समाधान का समाधान है?

library(sp) 
library(ggplot2) 

# Create two polygons: second would be a hole inside the first 
xy = cbind(
    x = c(13.4, 13.4, 13.6, 13.6, 13.4), 
    y = c(48.9, 49, 49, 48.9, 48.9) 
    ) 
hole.xy <- cbind(
    x = c(13.5, 13.5, 13.45, 13.45, 13.5), 
    y = c(48.98, 48.92, 48.92, 48.98, 48.98) 
) 

# Transform as SpatialPolygons with holes 
xy.sp <- SpatialPolygons(list(
    Polygons(list(Polygon(xy), 
       Polygon(hole.xy, hole = TRUE)), "1"), 
    Polygons(list(Polygon(xy + 0.2), 
       Polygon(xy + 0.35), 
       Polygon(hole.xy + 0.2, hole = TRUE)), "2") 
)) 

# Transform SpatialObject to be used by ggplot2 
xy.sp.l <- broom::tidy(xy.sp) 

ggplot(xy.sp.l) + 
    geom_polygon(aes(x = long, y = lat, group = id, fill = id)) 

Answer 1

एक अच्छा समय sf पैकेज के लिए "पर जाने के" करने के लिए हो सकता है:

यहाँ एक प्रतिलिपि प्रस्तुत करने योग्य उदाहरण है। sf वस्तु के साथ काम करते हुए, ggplot में बहुत आसान वास्तव में करने के लिए geom_sf ज्यामिति धन्यवाद:

library("sf") 
library("rgeos") 
sf_poly <- as(xy.sp, "sf") 
sf::st_crs(sf_poly) <- 4326 
sf_poly$id <- c(1,2) 
ggplot(sf_poly) + 
    geom_sf(aes(fill = as.factor(id))) 

Answer 2

लाइनों को जोड़ने समस्या के स्रोत को दर्शाता है। नीला "बहुभुज" खींचा जा रहा है, निचला -> ऊपरी -> छेद।

इस कोड (जो बहुत, सुरुचिपूर्ण खेद नहीं है) पथ वापस पहले टुकड़ा का प्रारंभिक बिंदु करने के लिए 3 के लिए आगे बढ़ने से पहले जाना पड़ता है।

library(dplyr) 
    extra <- xy.sp.l %>% 
     filter(piece != 1) %>% 
     group_by(id, group) %>% 
     summarise(last_pt = max(order)) 


for (n in 1:nrow(extra)) { 
    id_ex <- as.character(extra[n,"id"]) 
    x <- subset(xy.sp.l, id == id_ex & piece == 1 & order == 1) 
    x$order <- as.numeric(extra[n,"last_pt"]) + 0.5 
    xy.sp.l <- rbind(xy.sp.l,x) 
} 

xy.sp.l <- xy.sp.l[order(xy.sp.l$id, xy.sp.l$order),] 

Answer 3

इस पोस्ट में एक अच्छा सवाल है और पहले से ही महान उत्तर प्राप्त हुए। मुझे यह भी विश्वास है कि लोगों को sf ऑब्जेक्ट्स के साथ काम करना सीखना चाहिए क्योंकि यह आर में स्थानिक डेटा प्रकार की अगली पीढ़ी है। लेकिन मैं इस मामले में geom_spatialggspatial पैकेज से स्पेटियल पॉलीगॉन प्लॉट करने का विकल्प हो सकता हूं।

library(sp) 
library(ggplot2) 
library(ggspatial) 

ggplot() + 
    geom_spatial(xy.sp, aes(fill = id)) 
# Ignoring argument 'mapping' in geom_spatial.SpatialPolygons 
# Autodetect projection: assuming lat/lon (epsg 4326)