ईनार द्वारा पोस्ट की गई विधि से एक से अधिक पंक्तियां दिखाने की कोशिश करते समय, मैं दोनों को एक साथ नहीं दिखा सका। खुदाई के बाद, मैं this और this पदों पर आया। यहां दो अलग-अलग लाइनों के लिए एक छोटा कोड है।
library(geosphere)
library(leaflet)
library(dplyr)
lat_ny <- 40.73
lng_ny <- -73.9
lat_del <- 28.63
lng_del <- 77.21
lng_ca <- -121.6406
lat_ca <- 39.16414
inter1 <- gcIntermediate(c(lng_ny, lat_ny), c(lng_del, lat_del), n=10, addStartEnd=TRUE, sp = TRUE, breakAtDateLine = TRUE)
lines(inter1)
inter2 <- gcIntermediate(c(lng_ca, lat_ca), c(lng_del, lat_del), n=10, addStartEnd=TRUE, sp = TRUE, breakAtDateLine = TRUE)
lines(inter2)
inters <- c(inter1,inter2)
ll0 <- lapply(inters , function(x) `@`(x , "lines"))
ll1 <- lapply(unlist(ll0) , function(y) `@`(y,"Lines"))
Sl <- SpatialLines(list(Lines(unlist(ll1) , ID = 1)))
leaflet(Sl) %>% addTiles() %>% addPolylines()
हार्डकोड अक्षांश और देशांतर एक अच्छा विचार नहीं है, लेकिन जैसा कि मैंने केवल जगह शीर्ष 5 कनेक्ट लेने के लिए किया था, मैं एक सूची का अनुक्रमण के लिए ज्यादा समय समर्पित नहीं किया। इसके अलावा मैं अभी भी जांच कर रहा हूं कि यह चमकदार के साथ एकीकृत है या नहीं।
स्रोत
2017-05-10 06:00:39
शानदार! आप gcIntermediate को पुस्तिका में पास करने के बारे में कैसे जानते थे? –
जीसीइंटरडिएट द्वारा लौटाए गए ऑब्जेक्ट का वर्ग "स्पैटियललाइन" है। और पुस्तिका उस वर्ग को पसंद करती है (या कोई अन्य "स्थानिक ..." बहुत अधिक :-) – einar