পাইথনের একমাত্র সরঞ্জাম দিয়ে জাতীয় উচ্চতা ডেটাসেট (আরকগ্রিড / গ্রিডফ্লোট / আইএমজি) পড়া?

10

আমি ইউএসজিএস দ্বারা সরবরাহিত জাতীয় উচ্চতা ডেটাসেট থেকে উচ্চ নির্ভুলতার উচ্চতা (1/3 - 1/9 আর্কসেকেন্ড রেজোলিউশন) ডেটা পেয়েছি। এটি আইএমজি, আর্কগ্রিড এবং গ্রিডফ্লোট ফর্ম্যাট আসে। আমি জানি আর্কগ্রিড একটি প্রদেয় সফ্টওয়্যার প্যাকেজের সাথে মিলে যায় তবে আমি অবাধে উপলব্ধ সরঞ্জামগুলি ব্যবহারের সাথে আঁকড়ে থাকার চেষ্টা করছি।

আমার জিপিএস ডেটা রয়েছে যা আমি স্থল স্তরের সাথে সম্পর্ক স্থাপন করার চেষ্টা করছি।

এমন কোনও অজগর লাইব্রেরি রয়েছে যা আমাকে আইএমজি, আর্কগ্রিড, বা গ্রিডফ্লোট ফর্ম্যাটে কোনও ডেটা ফাইল ব্যবহার করে জিপিএস ডেটা উচ্চতায় রূপান্তর করতে দেবে?

python elevation esri-grid-format

— abnowack
সূত্র

7

অজগরে সরাসরি আইএমজি ফাইলের সাথে কাজ করা জিডিএল বাইন্ডিংগুলির সাথে সোজা। উদাহরণস্বরূপ, আপনি সরাসরি একটি NumPy অ্যারে ডেটা পড়তে পারেন:

from osgeo import gdal
geo = gdal.Open('imgn36w100_11.img')
arr = geo.ReadAsArray()
print repr(arr)
array([[ 744.31896973,  743.68762207,  743.1116333 , ...,  550.42498779,
         553.77813721,  556.18640137],
       [ 744.22955322,  743.66082764,  743.05273438, ...,  552.05706787,
         554.81365967,  557.55877686],
       [ 744.0133667 ,  743.49041748,  743.00061035, ...,  553.0123291 ,
         555.78076172,  558.01312256],
       ...,
       [ 568.70880127,  567.33666992,  566.56170654, ...,  447.68035889,
         447.68804932,  447.65426636],
       [ 568.01116943,  566.95739746,  564.23382568, ...,  447.6696167 ,
         447.71224976,  447.62734985],
       [ 565.62896729,  562.65325928,  560.78759766, ...,  447.67129517,
         447.67529297,  447.65179443]], dtype=float32)

আইএমজি ফর্ম্যাট ডেটা প্লট করার আরও সম্পূর্ণ উদাহরণের জন্য, এই স্ক্রিপ্টটি দেখুন যা নীচের চিত্রটি উত্পন্ন করেছে। আপনার জিপিএস ডেটাটি উচ্চতায় রূপান্তর করার জন্য আপনাকে ফলাফলযুক্ত নুমপি অ্যারে নমুনা করতে হবে।

এখানে চিত্র বর্ণনা লিখুন

— Kelsey
সূত্র

4

জাতীয় মানচিত্র দর্শক একটি ভিন্ন সফ্টওয়্যার প্ল্যাটফর্ম ব্যবহার করতেন, যা টাইলস নির্বাচনের পরে ডেটা জিওটিআইএফএফ-তে রূপান্তর করা সম্ভব করেছিল। কিছু আঞ্চলিক ডেটা সরবরাহকারী জিওটিআইএফএফকেও পছন্দ করে না। তবে এটি আপনার কোনও উপকারে আসতে পারে না।

তবে, সুসংবাদটি হ'ল gdalএই তিনটি ফর্ম্যাটের সাথে কাজ করতে সক্ষম। আর্কগ্রিড হ'ল ইএসআরআই বাইনারি ফর্ম্যাট, গ্রিডফ্লাট ইউএসজিএসের ডেম ডেটা ফর্ম্যাট এবং আইডিজির স্বত্বাধিকারী এরদাস, যারা ইমেজিন সফ্টওয়্যার তৈরি করে। আমি সম্ভবত আরকগ্রিড বা আইএমজি ফর্ম্যাটটি ডাউনলোড করে জিওটিআইএফএফ-তে রূপান্তর করব: ইউএসজিএসের নিজেরাই gdal ডাউনলোডের ভিত্তিতে রূপান্তর স্ক্রিপ্টগুলির একটি বান্ডিল রয়েছে ( জিপ ফাইলের সরাসরি লিঙ্ক )। সেখানে আপনার প্রয়োজনীয় সংস্থানগুলি খুঁজে পাওয়া উচিত।

— chryss
সূত্র

লিঙ্কটির জন্য ধন্যবাদ এবং জিডিএল উল্লেখ করার জন্য, এটি খুব দরকারী বলে মনে হচ্ছে

— অবনাক করুন

0

কেলসির উত্তর আমার পক্ষে প্রথমে কার্যকর হয়নি কারণ .img ফাইলটি -inf শূন্য স্থানের মান হিসাবে ব্যবহার করে।

আমি লাইনটি যুক্ত করে এটি সমাধান করেছি:

topo[topo<-1000] = 0

— vatlark
সূত্র

1

এটি মূল প্রশ্নের কোনও উত্তর নয়, এটি সম্ভবত স্বীকৃত উত্তরের একটি সম্পাদনা হওয়া উচিত

— nmtoken