Как рассчитать значения наклона и соотношения сторон, указав в качестве входных данных данные о широте, долготе и высоте?

Question

Поскольку я не очень разбираюсь в расчете наклона, может ли кто-нибудь помочь мне с кодом, который может рассчитывать значения наклона и пропорции?

У меня есть значения широты, долготы и высоты в файле.Есть ли способ получить эти значения?

Я много искал и узнал, что в R есть функция "terrain", которая дает наклон и соотношение сторон.Но у меня везде есть значения нан.

Answer 1

используя некоторые фиктивные данные:

library(raster)

long <- rep(rep(seq(12,36,0.5),41))
lat <-rep(seq(32,52,0.5), each=49)
dat <- rnorm(2009, 26.5, 44.0)

data <- data.frame(long, lat, dat)

rast <- rasterFromXYZ(data)
crs(rast) <- "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs "

slope <- terrain(rast,opt = 'slope', unit = 'degrees')
aspect <- terrain(rast, opt = 'aspect', unit = 'degrees')
flowdir <- terrain(rast, opt = 'flowdir')

spplot(stack(rast, slope, aspect, flowdir), scales = list(draw = T)) # the legend is not correct - only for demonstration purposes