지구 표면 등급의 히스토그램은 어디에서 찾을 수 있습니까?

Wikipedia 에서 지구 표면의 고도에 대한 막대 그래프를 찾았습니다 .

고도 히스토그램

그러나 이것은 등급 분포에 대한 정보를 제공하지 않습니다. 예를 들어, 전체 표면은 전적으로 작은 언덕으로 만들어져 어디에서나 높은 등급을 받거나 표면이 완전한 고원으로 구성되어 평균 등급이 0 °가 될 수 있습니다. 분명히이 두 시나리오는 모두 사실이 아니지만, 고도 히스토그램에서만이 정보를 확인할 수없는 방법을 보여줍니다. 학년에 대한 비슷한 히스토그램을 어디서 찾을 수 있는지 아는 사람이 있습니까?

data

— dlras2
소스

표고 히스토그램과 달리 경사 (기울기)는 기울기가 계산되는 해상도에 따라 다릅니다. 어떤 해결책이 필요합니까? 해저 경사도 필요합니까?

— whuber

아주 좋은 지적입니다. 이상적으로는 해저를 포함하지 않는 막대 그래프로 1 초 정도의 해상도를 원합니다 (그러나 해저에 대한 별도의 것도 좋습니다). 그러나 지금은 할 수있는 모든 것을 취할 것입니다 나를 올바른 방향으로 인도하십시오.

— dlras2

데이터 세트를 확보 할 수있는 경우이를 수행 할 수있는 도구가 있습니다 R. 나는 GeoTIFF로 Etopo1을 가지고 있는데, 바이너리 형식에서 직접 변환했을 수도 있지만 여기에서 얼음 / 셀이라고 생각합니다.

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더 많은 플로팅 옵션을 참조하십시오 .

— mdsumner
소스

참조 된 데이터에서 내 자신을 생성하는 것이 가장 좋은 방법처럼 보입니다. GeoTIFF, Etopo1 또는 R (Googling은 매우 사소한 것으로 판명 되었음)과 같은 GIS와 함께 일한 적이 없었습니다.

— dlras2

불행하게도,이 방법은 데이터를 투영하지 않기 때문에 잘못된 기울기를 얻습니다. 또 다른 합병증은 합리적인 해상도로 전세계 경사면을 얻는 것이 엄청난 노력이라는 것입니다. 수백 미터 이상에서 계산 된 경사는 아래쪽으로 부드럽게되는 경향이 있습니다. (그림의 격자는 7 킬로미터의 해상도를가집니다!) 100m 격자로 지구의 지표면을 덮으려면 약 360 억 개의 셀로 구성된 수천 개의 격자 (각각 합리적인 정확도를 위해 자체 투영이 필요)가 필요합니다. 이러한 DEM을 수집하는 것은 많은 작업입니다 ...

— whuber