我有两个30m x 30m的光栅文件,我想从其中采样点。在采样之前,我想从图像中删除云区。我转向R和Hijman的Raster包来完成这项任务。使用drawPoly(sp=TRUE)命令,我绘制了18个不同的多边形。该函数似乎不允许18个多边形作为一个sp对象,所以我将它们全部分开绘制。我通过lapply函数运行该列表,将它们转换为栅格(Hijman包中的rasterize函数),将多边形区域设置为NA,并将图像<e
在R中,但是我很难让它对每个栅格应用一个简单的条件语句(所有栅格层都是相同的范围,分辨率,并且是共同注册的)。我希望将所有小于每个栅格最大值的95%的像素值设置为NA。下面是我的代码: #Get max value from each raster layer
#Set all values < 95% of max