मानचित्र टाइल निर्माण गुणों को समझना


10

मैं एक टाइल बेस मैप बनाने की कोशिश कर रहा हूं। मेरी एक छवि है जिसे मैं टाइल बनाना चाहता हूं और मैं GDAL का उपयोग कर रहा हूं।

मैंने टाइल्स बनाने में कामयाबी हासिल की लेकिन जब मुझे यह नतीजा दिखता है तो लगता है कि टाइल वाला नक्शा है:

  1. यह फ़्लिप है
  2. अनुपात सही नहीं हैं

यह मानचित्र ओवरले छवि है: नक्शा छवि

यह टाइलिंग का परिणाम है:

टाइल परिणाम

मैंने कई बार जाँच की कि ऑर्थोपोटो पर छवि के कोनों की तुलना करने से निर्देशांक अच्छे हैं।

एक और बात यह है कि गूगल पृथ्वी में समान निर्देशांक के साथ मैप छवि ठीक गूगल मैप के लिए फिट बैठता है।

गूगल अर्थ ओवरले

संपादित

gdalinfo map_overlay.png

इस आउटपुट के साथ:

चालक: PNG / पोर्टेबल नेटवर्क ग्राफिक्स फ़ाइलें: map_overlay.png का आकार 3527 है, 2494 समन्वय प्रणाली है `छवि संरचना मेटाडेटा:
INTERLEAVE = PIXEL कॉर्नर निर्देशांक: ऊपरी बाएँ (0.0, 0.0, 0.0), निचला बाएँ (0.0, 2494.0) ऊपरी दाएँ (3527.0) , 0.0) निचला अधिकार (3527.0, 2494.0) केंद्र (1763.5, 1247.0) बैंड 1 ब्लॉक = 3527x1 प्रकार = बाइट, ColorInterp = लाल मास्क फ़्लैग्स: PER_DATASET एल्पेस बैंड 2 ब्लॉक = 3527x1 प्रकार = बाइट, ColorInterp = ग्रीन मास्क फ़्लैग: PER_DATASET बैंड 3 ब्लॉक = 3527x1 प्रकार = बाइट, कलरइंटरप = ब्लू
मास्क फ़्लैग्स: PER_DATASET अल्फा बैंड 4 ब्लॉक = 3527x1 प्रकार = बाइट, कलरइंटरप = अल्फा

फिर अनुवाद करें:

gdal_translate -of VRT -a_srs EPSG: 4326 -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 31.7493769676 35.1784535489 -gcp 3527 2492 31.7431011291 35494949543 -4949543 -4349543 -4949433364943

इस आउटपुट के साथ:

इनपुट फ़ाइल का आकार 3527, 2494 है

और अंत में:

gdal2tiles.py -z 14-21 map_overlay.vrt

ओवरले बिना किसी विशेष डेटा के एक सरल। Png छवि है।

किसी भी तरह की सहायता का स्वागत है।

शनि


2
gdal_translateआपके द्वारा यहां पोस्ट किए गए मापदंडों में एक असंगतता प्रतीत होती है : आप संबद्ध (lat, lon) = (31.7431761644, 35.1680410195) दो अलग-अलग बिंदुओं (0,0) और (02424) पर करते हैं। क्या यह वास्तव में आपके द्वारा जारी किया गया आदेश है, या यह एक कॉपी-एंड-पेस्ट त्रुटि है?
whuber

जवाबों:


5

निम्न के रूप में टाइल करने से पहले मीट्रिक समन्वित छवि का उपयोग करने का प्रयास करें:

   gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif 

अपडेट करें

नोट: भले ही हम gcp को निर्दिष्ट करते हैं, लेकिन टिफ के gdal_translateकोने निर्देशांक निर्दिष्ट नहीं करेंगे।

इसके अलावा आपको केवल एक फ़ाइल होने के कारण वर्चुअल डेटासेट (.vrt) बनाने की आवश्यकता नहीं है ...

निम्नलिखित कोड का उपयोग करें:

  gdal_translate -of  GTiff  -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 
       31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643
      -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.tif

फिर

  gdalwarp -s_srs epsg:4326 -t_srs epsg:4326  map_overlay.tif warped_map_overlay.tif

और अब gdalinfo का उपयोग करने का प्रयास करें ... अगर आपको कोई समस्या नहीं दिख रही है, तो यह कमांड चलाएँ:

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

यदि आपकी टाइलें फिट नहीं हैं, तो अपडेट से पहले कोड से ऊपर चलाएं, मैंने दिया है।

आप यहाँ अधिक जानकारी प्राप्त कर सकते हैं

मुझे उम्मीद है इससे आपको मदद मिली होगी....


धन्यवाद। मुझे ऐसा कब करना चाहिए? अभी के लिए मैं 1.gdalinfo map.png 2.gdal_translate -of VRT -a_srs EPSG: 4326 -gcp .... 3. gdal2tiles.py -z 10-16 map.vrt कर रहा हूं। क्या अनुवाद के बाद मुझे इसे करने की आवश्यकता है? क्या मुझे अनुवाद विधि में ईपीएसजी को बदलना चाहिए?
शन्नोगा

यदि आपका नक्शा समन्वित नहीं है, तो gdalinfo के साथ जानकारी लेने की कोशिश करना व्यर्थ है।
अरागॉन

इसलिए मुझे इसके बजाय गालवडप का उपयोग करने की आवश्यकता है?
शन्नोगा

आपका नक्शा किसी भी प्रक्षेपण के साथ समन्वित है या नहीं? यदि हाँ, इसे रद्द करने के लिए gdalwarp का उपयोग करने का प्रयास करें, अन्यथा इसे gdal_transform के साथ समन्वित करें ...
आरागॉन

अरे अपने मरीज के लिए धन्यवाद। यदि आप सही समझते हैं तो यह नहीं है। यह एक सरल वर्गाकार चित्र है, जो मानचित्र पर आधारित है, इस पर कोई विशेष डेटा नहीं है। मैंने अपने उपयोग किए गए आदेशों के साथ प्रश्न को अपडेट किया है। आपकी मदद के लिए फिर से धन्यवाद।
शन्नोगा
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.