Wo finde ich ein Histogramm der Erdoberfläche?

11

Ich habe auf Wikipedia ein Histogramm der Höhe der Erdoberfläche gefunden :

Höhenhistogramm

Dies gibt jedoch keine Auskunft über die Verteilung der Note. Zum Beispiel könnte die gesamte Oberfläche vollständig aus kleinen Hügeln bestehen und überall einen hohen Grad aufweisen, oder die Oberfläche könnte vollständig aus perfekten Plateaus bestehen, wobei der durchschnittliche Grad bei 0 ° liegt. Natürlich sind diese beiden Szenarien nicht wahr, aber es zeigt, wie diese Informationen nicht allein aus dem Höhenhistogramm ermittelt werden können. Weiß jemand, wo ich ein ähnliches Histogramm für die Note finden kann?

data dlras2
quelle
5
Im Gegensatz zum Höhenhistogramm hängt der Grad (Steigung) von der Auflösung ab, mit der die Steigung berechnet wird. Welche Auflösung benötigen Sie? Benötigen Sie auch Hänge des Meeresbodens?
whuber
Das ist ein sehr guter Punkt. Im Idealfall möchte ich eine Auflösung von ungefähr einer Sekunde mit einem Histogramm ohne Meeresboden (aber ein separates für den Meeresboden wäre auch gut). Im Moment werde ich jedoch alles nehmen, was ich kann leite mich in die richtige Richtung.
dlras2

Antworten:

6

Wenn Sie auf den Datensatz zugreifen können, stehen hierfür Tools Rzur Verfügung. Ich habe Etopo1 als GeoTIFF, ich denke, es ist das Eis / die Zelle von hier, obwohl ich es möglicherweise selbst aus dem Binärformat konvertiert habe.

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

Lesen Sie die Daten (möglicherweise mit reduzierter Auflösung), berechnen Sie die Steigung und das Diagramm. 

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)

R-Rasterplot

Ich verwende readGDAL, da ich damit besser vertraut bin, aber Sie können Raster als Wrapper um das rgdal-Material verwenden, um die Reduzierung der Auflösung usw. zu handhaben, und es ist keine Verwendung im Speicher erforderlich.

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

Siehe ?histfür mehr Plotten Optionen.

mdsumner
quelle
Das Generieren meiner eigenen aus den referenzierten Daten scheint der beste Weg zu sein. Nachdem ich noch nie mit einem GIS gearbeitet habe - GeoTIFF, Etopo1 oder R (was Googeln als angenehm trivial erwiesen hat) -, werde ich eine Weile brauchen, um es herauszufinden, aber Sie haben mich auf den richtigen Weg gebracht, denke ich.
dlras2
1
Leider erhält dieser Ansatz die falschen Steigungen, da die Daten nicht projiziert werden. Eine weitere Komplikation besteht darin, dass es eine enorme Anstrengung ist, weltweite Steigungen mit einer angemessenen Auflösung zu erhalten . Steigungen, die über mehr als einige hundert Meter berechnet werden, werden tendenziell nach unten geglättet. (Das abgebildete Gitter hat eine Auflösung von 7 Kilometern !) Um die Landoberfläche der Erde mit 100-Meter-Gittern zu bedecken, sind Tausende von Gittern erforderlich (von denen jedes seine eigene Projektion für eine angemessene Genauigkeit benötigt), die ungefähr 36 Milliarden Zellen umfassen. Nur diese DEMs zu sammeln ist eine Menge Arbeit ...
whuber