पाइथन का उपयोग करते हुए एक यूमेटैट छवि के पिक्सेल के लिए अक्षांश और लम्बाई कैसे प्राप्त करें?

Question

मेरे पास भूगर्भीय eumetsat उपग्रह की 3712x3712 पिक्सेल आकार की छवि है। पृथ्वी के चारों ओर कुछ काला है, इस तरह कि छवि इस तरह दिखता है:

पृथ्वी के प्रत्येक पिक्सेल के लिए, मैं अपने latitute और देशांतर प्राप्त करना चाहते हैं। मुझे पता है कि pyproj है कि और मैं बहुत की तरह एक प्रक्षेपण का दृष्टांत करने में सक्षम था:

sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0') 

लेकिन पिक्सेल के latlon geting (sat(x,y,inverse=True) जहां x और y छवि में पिक्सेल के निर्देशांक हैं का प्रयोग करके) स्पष्ट रूप से प्रक्षेपण के बाद से संभव नहीं है मेरी छवि के आयाम (3712x3712) को नहीं पता है।

मुझे क्या याद आ रही है?

Answer 1

मुझे लगता है कि आप सही प्रोजेक्शन लाइब्रेरी और सेटिंग्स का उपयोग कर रहे हैं।

ठेठ पिक्सेल रिज़ॉल्यूशन (पिक्सेल प्रति किमी की) eumetsat वेबसाइट here यह 3 किमी के बारे में प्रति पिक्सल की रिपोर्ट पर सूचना दी है।

आप इसे मैरीडियन पर एक्स, वाई में लोन/लैट रूपांतरण करके और पिक्सेल गिनती से इसे विभाजित करके देख सकते हैं (-81 डिग्री, 81 डिग्री अधिकतम सीमा है, संदर्भों के लिए eumetsat साइट देखें, http://www.eumetsat.int/) :

import pyproj 
sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0') 
x,y = sat(81.299, 0, radians = False, errcheck = True) 
print (x * 2.0/3712.0)/1000.0 

आपको 2.927 का मूल्य देगा जो यूमेटैट द्वारा दी गई जानकारी को फिट करेगा।

अगला आप अपने मानचित्र के प्रसिद्ध बिंदुओं (उदाहरण के लिए तटीय विशेषताओं) के सेट को परिभाषित करके आगे बढ़ सकते हैं, अपनी एक्स/वाई स्थिति निर्धारित कर सकते हैं और अपने अक्षांश/लोन निर्देशांक ऑनलाइन देख सकते हैं। आप पिक्सेल संकल्पों की एक श्रृंखला के साथ प्रयास कर सकते हैं और जांच सकते हैं कि कौन सा सबसे अच्छा फिट बैठता है या अधिक विस्तृत दिनचर्या का उपयोग करता है।

संकल्प इस बात पर निर्भर हो सकता है कि आप भूमध्य रेखा के करीब कितने करीब हैं, here देखें)। तो आपको कई अक्षांशों पर उपरोक्त दिनचर्या का उपयोग करने की आवश्यकता हो सकती है।