我需要一个线性回归来计算经验参数.L1是光栅图像,格式为.tif.L2也是预先计算的光栅图像.两个图像具有相同数量的列和行.
公式为:L1 = a + b*L2,其在R中翻译为:
lm(L1 ~ L2)
在第二个式I以后需要一个第二b.
我现在面临的问题是,两个栅格都包含NA值,我不知道如何为线性回归构建函数.我对R不熟悉所以我坚持这个可能相当简单的问题.我想我必须使用calc,但不知道怎么做.
编辑:到目前为止我有这个代码:
s = stack(L1,L2)
fun = function(x) {if (is.na(x[1])) { NA } else {lm(x[1] ~ x[2])$coefficients[2]}}
但是,计算需要很长时间才能得出结果
有没有办法改变R中NA值的颜色？我已经屏蔽了光栅堆栈中的一些值，并且想知道是否可以将屏蔽值的颜色从白色更改为黄色？我正在使用plotRGB 来绘制我的蒙版光栅文件。