Nếu bạn có thể giữ tập dữ liệu, có các công cụ R
để thực hiện việc này. Tôi có Etopo1 dưới dạng GeoTIFF, tôi nghĩ rằng đó là băng / ô từ đây mặc dù tôi có thể đã tự chuyển đổi nó từ định dạng nhị phân.
http://www.ngdc.noaa.gov/mgg/global/global.html
Đọc dữ liệu (có thể với độ phân giải giảm), tính toán độ dốc và âm mưu.
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)
Tôi sử dụng readGDAL vì tôi quen thuộc hơn với nó, nhưng bạn có thể sử dụng raster như một trình bao bọc xung quanh công cụ rgdal để xử lý giảm độ phân giải, v.v. và không yêu cầu sử dụng bộ nhớ.
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
Xem ?hist
để biết thêm các tùy chọn âm mưu.