विशिष्ट लंबे / लैट समन्वय के लिए प्रहरी टाइल ढूँढना


11

मैं एक टूल बना रहा हूं जो "सेंटिनल -2 पर AWS" से चित्र डाउनलोड करता है

मुझे एक विशिष्ट प्रहरी टाइल के लिए लांग / लैट निर्देशांक का अनुवाद करने में परेशानी होती है। मुझे पता है कि सेंटिनल 2 एमजीआरएस का उपयोग करता है और मैंने बस लांग / लैट को एमजीआरएस में बदलने की कोशिश की और कम परिशुद्धता के साथ मुझे एस 2 ए टाइल आईडी मिलती है। और यह काम करता है, लेकिन हमेशा नहीं।

उदाहरण के लिए, लॉन्ग / लैट 34.665,31.625 एमजीआरएस में "36 आर एक्सए" का समाधान करता है, लेकिन सेंटिनल टाइल 36 आरएक्सए मौजूद नहीं है।

लांग / लैट निर्देशांक का उपयोग करके प्रहरी टाइल का निर्धारण करने का सही तरीका क्या होगा?


1
सेंटिनल -2 टाइलिंग ग्रिड के साथ इंटरसेक्ट करें ?
user2856

हाँ, मैंने इस KML को देखा। लुकिंग महंगे होने के बाद से यह मेरा आखिरी उपाय होगा। मुझे आशा है कि मैं इसे गणना के साथ कर सकता हूं
GetUsername

3
की जाँच करें github.com/sat-utils/sentinel2-util
bugmenot123

आप Google क्लाउड में ऐसा कर सकते हैं क्योंकि उनके पास BigQuery में सभी मेटा डेटा हैं, उदाहरण के लिए SELECT * FROM [bigquery-public-data:cloud_storage_geo_index.sentinel_2_index] where north_lat < 50 and north_lat > 49 order by sensing_time LIMIT 10, लेकिन आपने AWS का उल्लेख किया है, इसलिए मुझे लगता है कि यह मदद नहीं करेगा।
जॉन पॉवेल

जवाबों:


4

MGRS में S2 विभाजन विशेष रूप से मुश्किल है जब डेटा UTM ज़ोन और टाइल ओवरलैप के बीच होता है, तो नीचे दिया गया आंकड़ा देखें। बाएं से दाएं टाइलें 20_M_QE, 20_M_RE, 21_M_SV और 21_M_TV हैं:

क्षेत्र के लिए ओपन सर्च के माध्यम से की गई खोज दो टाइलों, 20_M_RE और 21_M_TV में टूटे हुए S2 अधिग्रहण को लौटाती है:

20_M_RE

यहाँ छवि विवरण दर्ज करें

इसलिए यदि आप प्रतिनिधित्व किए गए 21_M_SV टाइल क्षेत्र में डेटा में रुचि रखते हैं तो आपको 20_M_RE का उपयोग करना चाहिए। ध्यान दें कि उदाहरणों में प्रयुक्त प्रतिनिधित्व (NGA से प्राप्त आकृति) टाइल्स के बीच ओवरलैप का प्रतिनिधित्व नहीं है।

टाइलें प्राप्त करने का सबसे सुरक्षित तरीका लगता है कि ओपनोग्राफी का उपयोग करके अपने एओआई को खोजा जाए और फिर लौटी हुई टाइलों को पुनः प्राप्त किया जाए, या बगमिनोटजैस द्वारा उल्लिखित सेंटिनल 2-उपयोग टूल का उपयोग किया जाए।

अगर आप ईएसए द्वारा प्रदान की जाने वाली केएमएल का उपयोग टाइलिंग ग्रिड परिभाषा के रूप में करने के लिए लॉन्ग / लैट टू टाइल्स का नक्शा चुनते हैं:

https://sentinel.esa.int/web/sentinel/missions/sentinel-2/data-products


4

ऐसा उपकरण पहले से मौजूद है। इसे सेंटिनेलसैट कहा जाता है और स्रोत गीथहब पर उपलब्ध है । यह एक कमांड लाइन इंटरफेस और पायथन एपीआई प्रदान करता है। यह प्रहरी 1 और 2 के साथ काम करता है। स्थानिक प्रश्न बहुभुज पर आधारित है और एक बिंदु पर नहीं है, लेकिन अन्यथा यह वही है जो आपको चाहिए।

संपादित करें: 1) आप क्वेरी फ़ंक्शन का उपयोग करके उत्पाद आईडी (या उत्पाद आईडी सूची) वापस कर सकते हैं। और कमांड लाइन से, आप -d विकल्प को छोड़ सकते हैं ताकि डाउनलोड को लॉन्च न किया जा सके (और प्रत्येक छवि के पदचिह्न बनाने के लिए -f विकल्प का उपयोग करें) 2) भू-भाग के संबंध में, आप नीचे देख सकते हैं कि एक wkt हो सकता है बजाय सीधे इस्तेमाल किया। जैसा कि @Thingumabus द्वारा उल्लेख किया गया है, आप अपने लैट / लंबे निर्देशांक के आधार पर एक छोटे बहुभुज बना सकते हैं। लेकिन वास्तव में, मैंने इसे एक बिंदु जियोजासन के साथ परीक्षण किया है और यह ठीक काम भी करता है। इसलिए अपनी बात को बहुभुज में बदलने की कोई आवश्यकता नहीं है।

# search by polygon, time, and Hub query keywords
footprint = geojson_to_wkt(read_geojson('map.geojson'))
products = api.query(footprint,
                     date = ('20151219', date(2015, 12, 29)),
                     platformname = 'Sentinel-2',
                     cloudcoverpercentage = (0, 30))

जिसमें समायोजित किया जा सकता है

# search by point, time, and Hub query keywords

products = api.query('POINT({0} {1})'.format(longitude,latitude),
                     date = ('20151219', date(2015, 12, 29)),
                     platformname = 'Sentinel-2',
                     cloudcoverpercentage = (0, 30))

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

@ user1186 एपीआई को खिलाने के लिए lat / lon के चारों ओर एक छोटे से बॉक्स का निर्माण करता है? मूर्खतापूर्ण लगता है, लेकिन आप सभी की जरूरत हो सकती है ...
geozelot

@Thingumabus, मैं ऐसा कर सकता था लेकिन API NDVI मान लौटाएगा। मैं केवल टाइल आईडी चाहता हूं, अगर एपीआई के माध्यम से ऐसा करने का कोई तरीका है, तो मुझे जानना अच्छा
लगेगा
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.