आर

2015-08-23 1 views
10

में लीफलेट के साथ नासा टाइल्स मैं लीफलेट पैकेज के संबंध में कुछ मदद पूछना चाहता हूं। जब आप एक इंटरेक्टिव मानचित्र खींचते हैं, तो आप ऐसा कुछ कर सकते हैं।आर

library(leaflet) 
library(magrittr) 

m <- leaflet() %>% 
    setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>% 
    addTiles() 

m 

यदि आप एक तृतीय-पक्ष टाइल जोड़ना चाहते हैं, तो आप भी ऐसा कर सकते हैं। निम्न लिंक तृतीय-पक्ष टाइल्स (http://leaflet-extras.github.io/leaflet-providers/preview/index.html) के लिए विकल्प प्रदान करता है निम्न छवि OpenWeatherMap.Precipitation के साथ बनाई गई है।

### They work 
m %>% addProviderTiles("MtbMap") 
m %>% addProviderTiles("HikeBike.HikeBike") 
m %>% addProviderTiles("OpenWeatherMap.Precipitation") 

enter image description here

कड़ी में टाइल के कुछ विकल्प इस नासा द्वारा टाइल्स शामिल हैं। मैं उनमें से एक का उपयोग करना चाहता था। तो मैंने निम्नलिखित कोडों का प्रयास किया। दुर्भाग्य से उनमें से कोई भी काम नहीं किया।

### The default map appears, then a black layer appears on top of the default layer. 
m %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR") 
m %>% addProviderTiles("NASAGIBS.ModisTerraBands367CR") 

एकमात्र विकल्प जो काम कर रहा है वह निम्नलिखित है।

m %>% addProviderTiles("NASAGIBS.ViirsEarthAtNight2012") 

enter image description here

मेरा अगला प्रयास addTiles() के साथ कस्टम URL टेम्प्लेट का उपयोग करने के लिए किया गया था। यूआरएल उपरोक्त लिंक से है। लेकिन, यह भी सफल नहीं था; कोई त्रुटि संदेश नहीं दिखाई दिया, लेकिन टाइल में कोई बदलाव नहीं आया।

m %>%addTiles(urlTemplate = "http://map1.vis.earthdata.nasa.gov/wmts-webmerc/MODIS_Terra_CorrectedReflectance_Bands367/default/{time}/{tilematrixset}{maxZoom}/{z}/{y}/{x}.{format}", 
       tileOptions(minZoom = 1, maxZoom = 8)) 

मेरा अंतिम प्रयास निम्न था। यह डिफ़ॉल्ट मानचित्र दिखाता है, लेकिन एक अतिरिक्त टाइल या तो दिखाई नहीं दे रहा था।

leaflet() %>% 
addTiles() %>% 
setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>% 
addTiles(urlTemplate = "http://map1.vis.earthdata.nasa.gov/wmts-webmerc/MODIS_Terra_CorrectedReflectance_Bands367/default/{time}/{tilematrixset}{maxZoom}/{z}/{y}/{x}.{format}", 
     tileOptions(minZoom = 1, maxZoom = 8)) 

मेरा सवाल यह है कि यह विशेष रूप से नासा टाइल्स के साथ एक संभावित बग है या नहीं। वैकल्पिक रूप से, मुझे इन स्क्रिप्ट में किस संशोधन की आवश्यकता है? आपकी सहायता के लिए अग्रिम रूप से धन्यवाद।

अद्यतन

मैं एक website जो एक ही नासा टाइल्स का उपयोग करता पाया। मैंने NASAGIBS.ModisTerraTRUEColorCR निर्दिष्ट किया और निम्न छवि प्राप्त की। छवि दिखा रही है कि एक मेल अमेरिका से स्वीडन तक कैसे पहुंचा। जैसा कि आप देखते हैं, अमेरिका और यूरोप दोनों के लिए कोई छवि नहीं है। मुझे लगता है कि यही कारण है कि मैंने एक काला टाइल देखा। मैं जानना चाहूंगा कि क्या कोई नासा टाइल्स के कुछ विवरण जानता है। मैंने उस क्षेत्र को चुना जिसे मैं नासा छवि देख सकता हूं। लेकिन, मुझे कोई भाग्य नहीं था।

### I expected to see Japan area this time. 
foo <- leaflet() %>% 
     setView(lng = 137.37, lat = 35.93, zoom = 5) %>% 
     addTiles() 

foo %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR") 

enter image description here

अद्यतन 2

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

m <- leaflet() %>% 
    setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>% 
    addTiles() 

m %>% addProviderTiles("NASAGIBS.ModisTerraBands367CR") 

enter image description here

उत्तर

6

आपका अंतिम निष्कर्ष सही है: आप के लिए और अनुरोध के समय पर चित्र का अनुरोध क्या स्थान के आधार पर, उपग्रह या अभी तक छवि हासिल कर ली है नहीं हो सकता है। तो आप एक वास्तविक छवि या सिर्फ एक खाली हो सकता है।(यह GIBS API documentation में भी कहा गया है।)

हालांकि, आप एडप्रोवाइडरटाइल्स() के लिए 'समय' विकल्प के माध्यम से छवि का अनुरोध करने के लिए किस दिन निर्दिष्ट कर सकते हैं। निकटतम अतीत में कोई तारीख निर्दिष्ट करके, यदि आप वही पसंद करते हैं तो आप सभी स्थानों के लिए गैर-खाली छवियां प्राप्त कर सकते हैं।

इस वाक्य रचना है:

Specified today's date (or no date). Some images are not yet acquired.

यह इंग्लैंड में बादल छाए हुए हैं, जो होगा:

> library(leaflet) 
> library(magrittr) 

> m <- leaflet() %>% 
     setView(lng = 4.5, lat = 51, zoom = 1) %>% 
     addTiles() %>% 
     addProviderTiles("NASAGIBS.ModisTerraTrueColorCR", 
         options = providerTileOptions(time = "2015-08-31", opacity = 0.5)) 

> m 

लिखने के समय (2015/08/31) मैं इस परिणाम मिलता है अनुमान लगाया?

अधिकतर डेटा पहले से ही होता है, लेकिन अभी तक अलास्का के लिए कोई इमेजरी नहीं है।

Specified a date in the past. Imagery everywhere alreadycacquired.

अंत में, कारण है कि

m %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR") 

था काम नहीं था शायद: दूसरी ओर, मैं कल की तारीख निर्धारित करते हैं तो

options = providerTileOptions(time = "2015-08-30", opacity = 0.5) 

हम पूर्ण छवि मिलता है एक टाइपो की वजह से। यह

m %>% addProviderTiles("NASAGIBS.ModisTerraTrueColorCR") 
+0

इस समय विकल्प के लिए धन्यवाद! यह जानना बहुत उपयोगी है। :) – jazzurro