在哪里可以找到地球表面坡度的直方图?


11

我在Wikipedia上找到了地球表面高度的直方图:

高程直方图

但是,这没有提供有关成绩分布的任何信息。例如,整个表面可以完全由小山丘制成,到处都具有较高的坡度,或者表面可以完全由完美的高原组成,平均坡度为0°。显然,这两种情况都是不正确的,但它说明了如何无法仅从海拔直方图确定此信息。有谁知道我在哪里可以找到相似的直方图?


5
与海拔直方图不同,坡度(坡度)取决于计算坡度的分辨率。您需要什么分辨率?您还需要海底斜坡吗?
ub

这是非常好的一点。理想情况下,我想要一秒钟左右的分辨率,并且直方图不包括海床(但也可以为海床单独设置一个)。但是,目前,我将采取一切可能的措施引导我朝正确的方向前进。
dlras2 2011年

Answers:


6

如果可以掌握数据集,则可以使用其中的工具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更多绘图选项,请参见。


从参考数据中生成自己的数据似乎是最好的方法。以前从未使用过任何GIS(GeoTIFF,Etopo1或R(Googling证明是微不足道的)),这需要花费我一段时间才能弄清楚,但我认为您已经把我放在了正确的轨道上。
dlras2 2011年

1
不幸的是,这种方法的斜率不正确,因为它不会投影数据。另一个复杂因素是,要以合理的分辨率获得世界范围的斜坡是一项巨大的努力。计算超过数百米的坡度将趋于向下平滑。(所示的栅格具有7 公里的分辨率!)用100m栅格覆盖地球的陆地表面需要成千上万个栅格(每个栅格都需要使用自己的投影以获得合理的精度),其中包含约360亿个像元。仅仅收集这些DEM就是很多工作……
麻烦
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.