मुझे पृथ्वी की सतह के ग्रेड का हिस्टोग्राम कहां मिल सकता है?

Question

मुझे विकिपीडिया पर पृथ्वी की सतह के उत्थान का एक हिस्टोग्राम मिला है :

हालाँकि यह ग्रेड वितरण के बारे में कोई जानकारी नहीं देता है। उदाहरण के लिए, पूरी सतह को पूरी तरह से छोटी पहाड़ियों से बनाया जा सकता है और हर जगह एक उच्च ग्रेड होता है, या सतह को पूरी तरह से पूर्ण पठारों से बनाया जा सकता है, औसत ग्रेड को 0 ° पर रखा जा सकता है। जाहिर है कि ये दोनों परिदृश्य असत्य हैं, लेकिन यह दर्शाता है कि यह जानकारी केवल ऊंचाई हिस्टोग्राम से निर्धारित नहीं की जा सकती है। क्या किसी को पता है कि मुझे ग्रेड के लिए एक समान हिस्टोग्राम कहां मिल सकता है?

Comments

ऊंचाई हिस्टोग्राम के विपरीत, ग्रेड (ढलान) उस संकल्प पर निर्भर करता है जिस पर ढलान की गणना की जाती है। आपको किस संकल्प की आवश्यकता है? क्या आपको समुद्र तल के ढलान की भी आवश्यकता है?

— whuber

यह बहुत अच्छी बात है। आदर्श रूप से, मैं एक या दूसरे का एक संकल्प चाहूंगा, जिसमें हिस्टोग्राम समुद्र तल सहित नहीं है (लेकिन समुद्र तल के लिए एक अलग भी अच्छा होगा।) फिलहाल, मैं कुछ भी ले सकता हूं। मुझे सही दिशा में ले जाने के लिए।

— dlras2

Answer 1

यदि आप डेटा सेट को पकड़ सकते हैं Rतो ऐसा करने के लिए उपकरण हैं । मेरे पास Etopo1 एक GeoTIFF के रूप में है, मुझे लगता है कि यह यहाँ से बर्फ / सेल है, हालांकि मैंने इसे बाइनरी प्रारूप से खुद को परिवर्तित किया होगा।

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

डेटा पढ़ें (संभवतः कम संकल्प के साथ), ढलान और साजिश की गणना करें।

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)

मैं readGDAL का उपयोग करता हूं क्योंकि मैं इसके साथ अधिक परिचित हूं, लेकिन आप रिजोल्यूशन को कम करने के लिए rgdal सामान के चारों ओर रैपर के रूप में छड़ी कर सकते हैं और इतने पर स्मृति उपयोग की आवश्यकता नहीं है।

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, या आर (जो कि Googling ने सुखद रूप से तुच्छ साबित कर दिया है) से पहले कभी भी जीआईएस के साथ काम नहीं किया है - मुझे यह पता लगाने में थोड़ा समय लगेगा, लेकिन आपने मुझे सही रास्ते पर सेट किया है, मुझे लगता है।

— dlras2

दुर्भाग्य से, इस दृष्टिकोण को गलत ढलान मिलता है, क्योंकि यह डेटा को प्रोजेक्ट नहीं करता है। एक और जटिलता यह है कि एक उचित संकल्प पर दुनिया भर में ढलान प्राप्त करने के लिए एक बहुत बड़ा प्रयास है। कुछ सौ मीटर से अधिक दूरी पर गणना की जाने वाली ढलानों को नीचे की ओर चिकना किया जाएगा। (सचित्र ग्रिड में 7 किलोमीटर का रिज़ॉल्यूशन है!) 100 मीटर ग्रिड के साथ पृथ्वी की भूमि की सतह को कवर करने के लिए लगभग 36 बिलियन कोशिकाओं वाले हजारों ग्रिड (जिनमें से प्रत्येक को उचित सटीकता के लिए इसके स्वयं के प्रक्षेपण की आवश्यकता होती है) की आवश्यकता होती है। बस बहुत काम इन Dems एकत्रित कर रहा है ...

— whuber