我正在使用Landsat图像来计算净辐射.在我转换的反射光栅图层中有非常低的负值(例如:-0.000003).我想确保我的反射率为0 - 1,以减少我未来计算中的错误.
如何在R中将小于0的栅格值替换为"NA".类似于栅格计算功能.我不确定如何提供一个例子,但我确定你们其中一个可以帮助我,是吗？
这是我从Bastiaanssen等人得到的反射率的等式.(1998)
假设,pb1表示Landsat的带1的反射率,pi = 3.14 ...,lb1 =带1的辐射率,ESUN =带1的外大气值,dr =一年中的相对地球太阳距离.
#Calculate reflectivity per band. QC: Always 0 to 1
pb1 = (pi * lb1)/(ESUN1 * cos(solzen) * dr)
创建此栅格后,我想要做的就是将小于0的pb1值设置为NA.
救命？