Erdfernerkundung - Numerische Physik: Modellierung

372 KAPITEL 9. ETWAS ZUM SPIELEN
Abbildung 9.4:
Histogramme der
Bänder aus Abb. 3.48
• selbst bei einfachen arithmetischen Operationen bei der Glättung oder anderer Filterung
ist uint8 unpraktisch, da der Zahlenraum zu klein ist.
Daher bietet sich die Konversion in den konventionellen Datentyp double an:
>> X64 = double(X8) + 1
Sollen die Daten nach der Bearbeitung wieder in ein Bild überführt werden, so muss die
Conversion von double nach uint8 (oder entsprechend uint16) erfolgen gemäß
>> X8 = uint8(round(X64 - 1))
Wenn wir mit Hilfe von imwrite ein indiziertes Bild erzeugen, konvertiert MatLab die Werte
automatisch falls erforderlich.
§ 1197 Nach erfolgter Konversion lässt sich hist direkt anwenden. Die vollständige Sequenz
zur Erstellung des Histogramms ist also
>> LP1=reshape(P1(:,:),[],1); ←↪
>> LF1=double(LP1)+1; ←↪
>> hist(LF1,256) ←↪
>> axis([0 256 0 12000])
§ 1198 Die Histogramme der sieben Spektralbänder für die Szene aus Abb. 3.48 sind in
Abb. 9.4 gegeben. Der Peak bei den niedrigen Helligkeitswerten in den Kanälen 3 und höher
ist die Wasserfläche. Die Spikes im thermischen Infrarot sind kein numerisches Artefakt sondern
scheinen reale Quellen unterschiedlicher Temperatur zu reflektieren. Die relativ schmale
Verteilung der Helligkeitswerte im blauen Kanal Band 1 ist eine Folge der Dunststreuung.
9.2.4 Filtern
§ 1199 Zum Filtern empfiehlt es sich, mit den double statt den uint8 Daten zu arbeiten.
Hier wird nur das Filtern von Hand mit dem Kernel wie im Haupttext demonstriert; die
schnellere Verwendung einer Laplace-Transformation mit Hilfe des MatLab-Befehls filter
wird nicht berücksichtigt.
2. Juli 2008 c○ M.-B. Kallenrode