ฉันจะหาฮิสโตแกรมของระดับพื้นผิวโลกได้จากที่ไหน

ฉันพบฮิสโตแกรมของระดับความสูงของพื้นผิวโลกบนวิกิพีเดีย :

อย่างไรก็ตามสิ่งนี้ไม่ได้ให้ข้อมูลใด ๆ เกี่ยวกับการกระจายของคะแนน ตัวอย่างเช่นพื้นผิวทั้งหมดสามารถทำจากเนินเขาเล็ก ๆ ทั้งหมดและมีระดับสูงได้ทุกที่หรือพื้นผิวอาจประกอบไปด้วยที่ราบสูงที่สมบูรณ์แบบโดยวางเกรดเฉลี่ยที่ 0 ° เห็นได้ชัดว่าทั้งสองสถานการณ์นี้ไม่จริง แต่มันแสดงให้เห็นว่าข้อมูลนี้ไม่สามารถระบุได้จากฮิสโตแกรมระดับความสูงเพียงอย่างเดียว ไม่มีใครรู้ว่าฉันสามารถหาฮิสโตแกรมที่คล้ายกันสำหรับเกรดได้หรือไม่

หากคุณได้รับชุดข้อมูลที่มีเครื่องมือใน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ตัวเลือกการพล็อตเพิ่มเติม