> -----Ursprüngliche Nachricht----- > Von: [email protected] [mailto:[email protected]] Im Auftrag > von Markus Neteler > Gesendet: Mittwoch, 21. Januar 2015 16:46 > An: Jäger, Frank (KRZ) > Cc: [email protected] > Betreff: Re: [FOSSGIS-Talk] DOP 2014 NRW +IR ...
> > Mein Fehler war .., die Formel einfach .. von "convert" auf "gdal_calc" zu > > übertragen. > > Wenn Werte zwischen -1 und +1 heraus kommen, aber z.B. Pixel-Farbwerte > > 0-255 gebraucht werden, dann erklärt das, warum fast alles schwarz bleibt. > Man koennte die Werte auch mit 100 multiplizieren bzw. auf 0-255 skalieren in > gdal_calc.py, um wieder auf den 8bit range zu kommen. > Markus Moin! Die Evolution der Formel schreitet voran: Um den theoretisch möglichen Wertebereich -1 bis +1 auf die Helligkeitswerte 0 bis 255 zu projizieren habe ich versucht: --calc="255.0/2.0 * (1.0 + (A - B) / (A + B))" Der Zuschlag +1.0 hebt den Bereich auf 0 bis 2. Der Faktor 255/2 spreizt das auf den Zielbereich. Das Ergebnis-Bild hat nur 3 Farben: schwarz, grau, weiß. Dies deutet weiter darauf hin, dass Zwischenergebnisse als Integer gerundet werden. Ein erstes ansehnliches - aber noch nicht brauchbares -Ergebnis bekomme ich wenn ich durch 1.0-Multiplikatoren überall Gleitkomma erzwinge: --calc="255.0/2.0 * (1.0 + (A * 1.0 - B * 1.0) / (A * 1.0 + B * 1.0))" Das Ergebnis entspricht noch nicht dem aus "convert -fx". Die Helligkeitswerte konzentrieren sich auf den Bereich 120 - 233. Die negativen Ergebnisse aus der Kernformel "(A - B) / (A + B)" spielen also (fast) keine Rolle. Vielleicht sollte man nur die positiven Werte mappen? Mit der folgenden Formal teste ich jetzt mal in der Fläche: --calc="255.0*((A*1.0-B*1.0)/(A*1.0+B*1.0))" Hat - für den nächsten Schritt - jemand einen zum Thema Vegetation passenden RGB-Farbbereich um die Helligkeitswerte in Farben abzubilden, möglichst in Mapserver-Syntax? MfG F. Jäger -------------- nächster Teil -------------- Ein Dateianhang mit Binärdaten wurde abgetrennt... Dateiname : smime.p7s Dateityp : application/pkcs7-signature Dateigröße : 7599 bytes Beschreibung: nicht verfügbar URL : <https://lists.fossgis.de/pipermail/fossgis-talk-liste/attachments/20150122/a1c0d363/smime.p7s> -- .................................................................... 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 [email protected] https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
