Якщо ви можете отримати доступ до набору даних, для цього є інструменти 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
більше варіантів побудови графіків.