Де можна знайти гістограму класу поверхні Землі?


11

Я знайшов гістограму висоти поверхні Землі у Вікіпедії :

Гістограма висоти

Однак це не дає ніякої інформації про розподіл оцінок. Наприклад, вся поверхня могла бути виготовлена ​​повністю з невеликих пагорбів і мати високий сорт скрізь, або поверхня могла бути повністю складена з ідеальних плато, ставлячи середній ступінь 0 °. Очевидно, що обидва ці сценарії не відповідають дійсності, але це ілюструє, як ця інформація не може бути визначена лише з гістограми висоти. Хтось знає, де я можу знайти подібну гістограму для класу?


5
На відміну від гістограми висоти, ступінь (ухил) залежить від роздільної здатності, при якій обчислюється схил. Яка резолюція вам потрібна? Вам також потрібні схили океанського дна?
whuber

Це дуже вдалий момент. В ідеалі, я хотів би дозволити секунду або близько того, з гістограмою, що не включає океанське дно (але окрема для океанічного дна буде також добре.) На даний момент я візьму все, що можу приведи мене в правильному напрямку.
dlras2

Відповіді:


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)

R растровий сюжет

Я використовую 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більше варіантів побудови графіків.


Створення власних даних із посилань на дані виглядає як найкращий спосіб. Ніколи раніше не працював з жодним ГІС - GeoTIFF, Etopo1 або R (що Googling виявилося приємно тривіальним) - це знадобиться мені певний час, але ви поставили мене на правильний шлях.
dlras2

1
На жаль, такий підхід отримує неправильні нахили, оскільки він не проектує дані. Ще одне ускладнення полягає в тому, що здобути по всьому світу схили при розумній резолюції - це величезні зусилля. Схили, обчислені більше декількох сотень метрів, як правило, згладжуються вниз. (Ілюстрована сітка має 7- кілометрову роздільну здатність!) Для покриття земної поверхні 100 м сітки потребують тисячі сіток (кожна з яких потребує власної проекції для розумної точності), що містить приблизно 36 мільярдів комірок. Просто збирати ці DEM - це велика робота ...
whuber
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.