我有一个RasterBrick,由每天1,2和3值的积雪数据组成(1= ,2= no ,3=云模糊)。对于空间插值,我想使用数字高程模型(相同的研究区域和相同的分辨率)分别为RasterBrick 的每个层提取上、下雪线边界。上面的雪线表示所有无云像素都被归类为雪的海拔.较低的雪线识别高度,低于这个高度的所有无云像素也都是无雪的
我一直在尝试创建一个新的栅格对象,它只包含现有栅格中的几个值。我使用的是这里找到的类栅格:。# but the file is to bigs <- r[r>30 & r<33] # but the file is to bigrc <- reclass(r, c(-Inf,30,NA, 31,32, 1, 33,Inf