2015-01-21 11:29 GMT+01:00 Jäger, Frank (KRZ) <f.jae...@krz.de>: ... > gdal_calc.py -A nir.tif -B rot.tif --outfile=ndvi.tif --calc="(A - B) / (A > + B + 0.001)" > > Kostet wirklich nur einen Bruchteil der CPU. Das Ergebnis ist aber fast > schwarz mit ein paar weißen Flecken (Fahrbahnränder).
Falls auch gdal_calc.py die Formel von links nach rechts liest (wie auch r.mapcalc in GRASS GIS), ist es vielleicht notwendig, einen Faktor 1.0 vorzustellen, um eine floating point Division zu erreichen, sofern die Eingabekarten Integerwerte haben. Also (ungetestet): gdal_calc.py -A nir.tif -B rot.tif --outfile=ndvi.tif --calc=" 1.0 * (A - B) / (A + B)" (woher die 0.001 kommt, weiss ich nicht). NDVI Werte reichen von (theoretisch) -1.0 bis 1.0, sie werden aber in Mitteleuropa typischerweise nicht erreicht. Gruesse, Markus -- .................................................................... FOSSGIS 2015, Die Konferenz für Open Source GIS mit OpenData und OpenStreetMap in Münster! 11.-13. März 2015 im Schloss der Universität Münster http://www.fossgis.de/konferenz/2015/ FOSSGIS e.V, der Verein zur Förderung von Freier Software aus dem GIS-Bereich und Freier Geodaten! http://www.fossgis.de/ https://twitter.com/fossgis_eV ____________________________________________________________________ FOSSGIS-Talk-Liste mailing list FOSSGIS-Talk-Liste@fossgis.de https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
