Gdzie mogę znaleźć histogram stopnia powierzchni Ziemi?

11

Na Wikipedii znalazłem histogram wysokości powierzchni Ziemi :

Histogram wysokości

Nie daje to jednak żadnych informacji o rozkładzie ocen. Na przykład cała powierzchnia może być wykonana w całości z małych wzgórz i mieć wszędzie wysoką jakość, lub powierzchnia może być złożona w całości z doskonałych płaskowyżów, ustawiając średnią ocenę na 0 °. Oczywiście oba te scenariusze są nieprawdziwe, ale ilustruje to, jak tej informacji nie można ustalić na podstawie samego histogramu wysokości. Czy ktoś wie, gdzie mogę znaleźć podobny histogram dla oceny?

dlras2
źródło
5
W przeciwieństwie do histogramu wysokości nachylenie (nachylenie) zależy od rozdzielczości, przy której obliczane jest nachylenie. Jakiej rozdzielczości potrzebujesz? Czy potrzebujesz również stoków dna oceanu?
whuber
To bardzo dobra uwaga. Idealnie chciałbym mieć rozdzielczość około sekundy, z histogramem nie uwzględniającym dna oceanu (ale osobny byłby również dobry dla dna oceanu.) W tej chwili wezmę jednak wszystko, co mogę poprowadź mnie we właściwym kierunku.
dlras2

Odpowiedzi:

6

Jeśli możesz zdobyć zestaw danych, istnieją narzędzia Rdo tego. Mam Etopo1 jako GeoTIFF, myślę, że jest to ten lodowy / komórkowy stąd, chociaż mogłem go przekonwertować z formatu binarnego.

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

Odczytaj dane (ewentualnie ze zmniejszoną rozdzielczością), oblicz nachylenie i wykreśl.

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)

Raster działka

Korzystam z readGDAL, ponieważ jestem bardziej zaznajomiony z tym, ale możesz trzymać się rastra jako owijki wokół rgdal, aby poradzić sobie ze zmniejszaniem rozdzielczości i tak dalej, i nie wymaga użycia pamięci.

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 

Zobacz ?histwięcej opcji drukowania.

mdsumner
źródło
Generowanie własnego z danych referencyjnych wygląda na najlepszą drogę. Nigdy wcześniej nie pracowałem z żadnym GIS - GeoTIFF, Etopo1 lub R (który Googling okazał się przyjemnie trywialny) - zajmie mi to trochę czasu, ale myślę, że ustawiłeś mnie na dobrej drodze.
dlras2
1
Niestety takie podejście powoduje nieprawidłowe nachylenie, ponieważ nie wyświetla danych. Kolejną komplikacją jest to, że uzyskanie światowych stoków przy rozsądnej rozdzielczości jest ogromnym wysiłkiem. Nachylenia obliczone na ponad kilkuset metrach będą zwykle wygładzane w dół. (Zilustrowana siatka ma rozdzielczość 7 kilometrów !) Pokrycie powierzchni ziemi siatkami 100 m wymaga tysięcy siatek (każda z nich wymaga własnej projekcji dla rozsądnej dokładności) zawierających około 36 miliardów komórek. Samo zbieranie tych DEM to dużo pracy ...
whuber