Dünya yüzeyinin derecesinin histogramını nerede bulabilirim?


11

Wikipedia'da Dünya yüzeyinin yüksekliğinin bir histogramını buldum :

Yükseklik histogramı

Ancak bu notun dağılımı hakkında herhangi bir bilgi vermez. Örneğin, tüm yüzey tamamen küçük tepelerden yapılabilir ve her yerde yüksek bir dereceye sahip olabilir veya yüzey, ortalama dereceyi 0 ° 'ye koyarak tamamen mükemmel yaylalardan oluşabilir. Açıkçası bu senaryoların her ikisi de doğru değildir, ancak bu bilginin yalnızca yükseklik histogramından nasıl belirlenemeyeceğini gösterir. Sınıf için benzer bir histogramı nerede bulabileceğimi bilen var mı?


5
Yükseklik histogramının aksine, eğim (eğim) eğimin hesaplandığı çözünürlüğe bağlıdır. Hangi çözüme ihtiyacınız var? Ayrıca okyanus tabanının eğimlerine mi ihtiyacınız var?
whuber

Bu çok iyi bir nokta. İdeal olarak, okyanus tabanını içermeyen bir histogramla bir saniye kadar çözünürlük istiyorum (ancak okyanus tabanı için ayrı bir tane de iyi olurdu.) Şu anda, yapabileceğim her şeyi alacağım beni doğru yöne götür.
dlras2

Yanıtlar:


6

Veri kümesini tutabiliyorsanız, Rbunu yapmak için araçlar vardır . Bir GeoTIFF olarak Etopo1 var, ben kendimi ikili formattan dönüştürmüş olabilir rağmen buz / hücre bir olduğunu düşünüyorum.

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

Verileri okuyun (muhtemelen düşük çözünürlükle), eğimi hesaplayın ve çizin.

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 raster grafiği

Ben daha aşina olduğum için readGDAL kullanıyorum, ancak çözünürlük azaltma ve benzeri işlemek için rgdal şeyler etrafında bir sarmalayıcı olarak raster ile yapışabilir ve bellek kullanımını gerektirmez.

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 

Daha ?histfazla çizim seçeneği için bakınız .


Başvurulan verilerden kendimi oluşturmak en iyi yol gibi görünüyor. Daha önce hiç CBS ile çalışmadıysanız - GeoTIFF, Etopo1 veya R (Googling'in hoş bir şekilde kanıtladığı kanıtlandı) - anlaması biraz zaman alacak, ama beni doğru yola koydun, sanırım.
dlras2

1
Ne yazık ki, bu yaklaşım yanlış eğimler alır, çünkü verileri yansıtmaz. Diğer bir komplikasyon, dünya çapında eğimleri makul bir çözünürlükle elde etmenin muazzam bir çaba olmasıdır. Birkaç yüz metreden uzun bir sürede hesaplanan eğimler aşağı doğru yumuşama eğiliminde olacaktır. (Gösterilen ızgara 7 kilometrelik bir çözünürlüğe sahiptir!) Dünyanın kara yüzeyini 100m ızgaralarla kaplamak, kabaca 36 milyar hücre içeren binlerce ızgara gerektirir (her biri makul doğruluk için kendi projeksiyonuna ihtiyaç duyar). Sadece bu
DEM'leri toplamak
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.