পৃথিবীর পৃষ্ঠের গ্রেডের হিস্টোগ্রামটি কোথায় পাব?

Question

আমি উইকিপিডিয়ায় পৃথিবীর পৃষ্ঠের উচ্চতা বৃদ্ধির একটি হিস্টোগ্রাম পেয়েছি :

তবে এটি গ্রেড বিতরণ সম্পর্কে কোনও তথ্য দেয় না। উদাহরণস্বরূপ, পুরো পৃষ্ঠটি সম্পূর্ণ ছোট ছোট পাহাড় দিয়ে তৈরি হতে পারে এবং সর্বত্র একটি উচ্চ গ্রেড থাকতে পারে, বা পৃষ্ঠটি পুরোপুরি নিখুঁত মালভূমি দ্বারা গঠিত হতে পারে, গড় গ্রেড 0 at রেখে ° স্পষ্টতই এই উভয় পরিস্থিতি অসত্য, তবে এটি বর্ণনা করে যে কীভাবে এই তথ্যটি একাই উচ্চতা হিস্টোগ্রাম থেকে নির্ধারণ করা যায় না। কেউ কি জানতে পারে আমি গ্রেডের জন্য একই রকম হিস্টোগ্রামটি কোথায় পাব?

Comments

উচ্চতা হিস্টগ্রামের বিপরীতে, গ্রেড (opeাল) যে রেজোলিউশনে slাল গণনা করা হয় তার উপর নির্ভর করে। আপনার কি রেজোলিউশন দরকার? আপনার কি সমুদ্রের তল ?ালু প্রয়োজন?

— হোবার

এটি একটি খুব ভাল পয়েন্ট। আদর্শভাবে, আমি সমুদ্রের তলকে অন্তর্ভুক্ত না করে একটি হিস্টোগ্রাম সহ একটি দ্বিতীয় বা তার মতো একটি রেজোলিউশন চাই (তবে সমুদ্রের তলগুলির জন্য পৃথক পৃথকটিও ভাল হবে)) তবে এই মুহুর্তে, আমি যা কিছু করতে পারি তা গ্রহণ করব আমাকে সঠিক পথে নিয়ে যেতে।

— dlras2

Answer 1

আপনি যদি ডেটা সেটটি ধরে রাখতে পারেন তবে এটি করার জন্য সরঞ্জাম রয়েছে R। আমার জিওটিআইএফএফ হিসাবে ইটোপো 1 রয়েছে, আমি মনে করি এটি এখান থেকে বরফ / কোষের এক হলেও যদিও আমি এটি বাইনারি ফর্ম্যাট থেকে নিজেকে রূপান্তর করেছি।

http://www.ngdc.noaa.gov/mgg/global/global.html

ডেটা পড়ুন (সম্ভবত হ্রাসযুক্ত রেজোলিউশন সহ), opeাল এবং প্লটটি গণনা করুন।

library(rgdal)
library(raster)

## orig dims, reduced 4-fold (choose divisor to suit your needs / system)
x <- readGDAL("Etopo1.tif", output.dim = c(10800, 21600)/4)

## convert to raster format for calculations
r <- raster(x)

g <- slopeAspect(r, out = "slope", unit = "degrees")

## plot histogram
hist(g)

আমি রিডজিডাল ব্যবহার করি যেহেতু আমি এর সাথে আরও বেশি পরিচিত তবে আপনি আরজিডাল স্টাফের চারপাশে একটি র‍্যাপার হিসাবে রাস্টার দিয়ে আটকে রাখতে পারেন যাতে রেজোলিউশন হ্রাস করতে পারে এবং এ জাতীয় কিছু করা যায় না, এবং স্মৃতি ব্যবহারের প্রয়োজন হয় না।

g
class       : RasterLayer 
dimensions  : 2700, 5400, 14580000  (nrow, ncol, ncell)
resolution  : 0.06666667, 0.06666667  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 
values      : in memory
min value   : 0 
max value   : 38.11677 

দেখুন ?histআরো চক্রান্ত বিকল্পের জন্য।

Comments

রেফারেন্স করা ডেটা থেকে আমার নিজের তৈরি করা দেখতে যাওয়ার সর্বোত্তম উপায় বলে মনে হচ্ছে। জিওটিআইএফএফ, ইটোপো 1, বা আর (যা গুগলিং আনন্দদায়ক তুচ্ছ প্রমাণ করেছে) - এর সাথে এর আগে কখনও কোন জিআইএস নিয়ে কাজ করেনি, তা আবিষ্কার করতে আমার কিছুটা সময় লাগবে, তবে আপনি আমাকে সঠিক পথে রেখেছেন, আমি মনে করি।

— dlras2

দুর্ভাগ্যক্রমে, এই পদ্ধতির ভুল opালু পায়, কারণ এটি ডেটা প্রজেক্ট করে না। আরেকটি জটিলতা হ'ল যুক্তিসঙ্গত রেজোলিউশনে বিশ্বব্যাপী slালু প্রাপ্তি একটি বিশাল প্রচেষ্টা। কয়েক শ মিটারেরও বেশি জুড়ে গণনা করা opালগুলি নীচের দিকে মসৃণ হবে। (সচিত্র গ্রিডটিতে 7 কিলোমিটার রেজোলিউশন রয়েছে!) 100 মিটার গ্রিড দিয়ে পৃথিবীর স্থলভাগকে ingেকে ফেলার জন্য হাজার হাজার গ্রিড (যার প্রত্যেকটির যুক্তিসঙ্গত নির্ভুলতার জন্য নিজস্ব প্রক্ষেপণ প্রয়োজন) প্রায় 36 বিলিয়ন কোষ সমন্বিত রয়েছে। শুধু কাজ অনেক এই DEMs সংগ্রহ করছে ...

— whuber