Ich arbeite mit Landsat-Bildern, um die Nettostrahlung zu berechnen. In meinen konvertierten Reflexionsgrad-Rasterebenen befinden sich sehr niedrige negative Werte (z. B. -0,000003). Ich möchte sicherstellen, dass meine Reflexionsgrade 0 - 1 sind, um Fehler bei meinen zukünftigen Berechnungen zu reduzieren.

Wie kann ich Rasterwerte unter 0 durch "NA" in R ersetzen? Ähnlich wie bei einer Raster-Calc-Funktion. Ich bin mir nicht sicher, wie ich ein Beispiel nennen soll, aber ich bin sicher, einer von Ihnen könnte mir helfen, ja?

Hier ist meine Gleichung für das Reflexionsvermögen von Bastiaanssen et al. (1998)

Angenommen, pb1 gibt das Reflexionsvermögen für Band 1 von Landsat an, pi = 3,14 ..., lb1 = Strahlung für Band 1, ESUN = exoatmosphärischer Wert für Band 1, dr = die relative Entfernung der Erdsonne für den Tag des Jahres.

#Calculate reflectivity per band. QC: Always 0 to 1
pb1 = (pi * lb1)/(ESUN1 * cos(solzen) * dr)

Nachdem dieses Raster erstellt wurde, möchte ich nur die pb1-Werte, die kleiner als 0 sind, auf NA setzen.

Hilfe