আপনি যদি ডেটা সেটটি ধরে রাখতে পারেন তবে এটি করার জন্য সরঞ্জাম রয়েছে 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
আরো চক্রান্ত বিকল্পের জন্য।