Band math с calc и / или оверлеем, возвращающим однозначные изображения

Question

Когда я пытаюсь запустить банд математику, результатом всегда будет изображение цвета, а значения min и max сильно отличаются от предсказанных. Я не нашел здесь ни одного вопроса, который бы показал эту проблему. Я работал таким образом

r.stack <- stack("path to raster file"))

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

prj <- "+proj=utm +zone=23 +south +datum=WGS84 +units=m"
r <- raster(res=11.47, ext=extent(c(301496,  323919, 9888968, 9913982)), crs=prj, vals=NA
r.stack <- resample(r.stack, r)

После этого изображения имеют такую ​​конфигурацию:

> class       : RasterBrick 
> dimensions  : 2181, 1955, 4263855, 4  (nrow, ncol, ncell, nlayers)
> resolution  : 11.47, 11.47  (x, y)
> extent      : 301496, 323919.8, 9888966, 9913982  (xmin, xmax, ymin, ymax)
>coord. ref. : +proj=utm +zone=23 +south +datum=WGS84 +units=m +ellps=WGS84 +towgs84=0,0,0 
>data source : in memory
>names       : l.1, l.2, l.3, l.4
>min values  :       -36.12217,           -45.12768,   -46.30455,                       -35.26328 
>max values  :      10.567671,         4.050200,         3.878345,                       11.613799

и затем используйте функцию ниже для расчета

f <- function(x){
      (x[[2]])/(x[[1]])
}

s <- r.stack[[c(1,2)]]
r2 <- calc(s, f)

и я тоже запускаю оверлей с забавой

f <- function(x,y){
     y/x
}
r2 <- overlay(r.stack[[1]], r.stack[[2]], fun= f)

Любой из методов приводит к изображению с одним значением

Я пропустил несколько шагов?

Answer 1

Вот ваш код с некоторыми примерами данных (без этого сложно ответить на вопросы). Я немного упростил одну функцию, но результаты те же.

library(raster)
b <- brick(system.file("external/rlogo.grd", package="raster"))
b <- b/10 + 1

f <- function(x){ x[2]/ x[1] }
s <- b[[c(1,2)]]
r1 <- calc(s, f)

f <- function(x,y){ y / x }
r2 <- overlay(b[[1]], b[[2]], fun= f)

Или просто

r3 <- b[[2]] / b[[1]]   
r3
#class       : RasterLayer 
#dimensions  : 77, 101, 7777  (nrow, ncol, ncell)
#resolution  : 1, 1  (x, y)
#extent      : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=merc +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
#data source : in memory
#names       : layer 
#values      : 0.7692308, 1.7  (min, max)

r1 и r2 одинаковы.

Причина, по которой вы получаете «один цвет», состоит в том, что большинство значений близко к 1, но есть несколько больших выбросов; возможно из-за деления на число от -1 до 1? Это может проиллюстрировать это:

q <- quantile(r3, c(0.1, 0.9))
d <- clamp(r3, q[1], q[2])
plot(d)

И посмотрите на крайности

i <- which.max(r3)
b[i][,2:1]