Erdfernerkundung - Numerische Physik: Modellierung

9.2. BILDER IN MATLAB 371
>> P4rgb = cat(3,P4,P4,P4);
Abbildung 9.3: (Links:) Original
der Helligkeitswerte von
band 4; (rechts:) mit cat und
image erzeugtes RGB-Bild in
S/W
Das Resultat ist im rechten Teil von Abb. 9.3 gegeben; der linke Teil zeigt zum Vergleich das
‘Original’.
9.2.2 Zerfleddern kombinierter Kanäle
§ 1192 Beim Spiel mit normalen Photos, wie z.B. in Abb. 3.23 und 3.24, ist der Ausgang
bereits ein Farbphoto – das Helligkeitsmuster liegt in drei Kanälen vor. Daher müssen wir
vor einer weiteren Bearbeitung die Kanäle trennen. Als erstes schreiben wir inpig.jpg wieder
in die Variable PIG und holen uns aus dem Workspace die Information über das Format.
Nehmen wir kombi321 als Beispiel und lesen mit imread ein, so ist PIG ein dreidimensionales
Array mit dem Format 850 × 1100 × 3.
§ 1193 Aus diesem Array lassen sich die einzelnen Kanäle extrahieren; in einer vielleicht
etwas aufwendigen expliziten Form (alternativ mit For-Schleife und Index) als
>> KRED(:,:)=kombi321(:,:,1); ←↪
>> KGREEN(:,:)=kombi321(:,:,2); ←↪
>> KBLUE(:,:)=kombi321(:,:,3); ←↪
Damit sollten wir wieder die Ausgangskanäle erhalten haben: die Differenz KRED-P3 sollte
also verschwinden. Davon können wir uns z.B. mit der Sequenz
>>maximum=max(max(diff(KRED-P3)))
überzeugen.
9.2.3 Histogramme
§ 1194 Histogramme erstellt MatLab mit Hilfe des Befehls hist automatisch. Allerdings
sollten wir hist nicht auf ein Bild in voller Schönheit anwenden: MatLab verlangt als Eingabe
für hist einen Vektor oder eine Matrix; ein dreidimensionales Array ist nicht vorgesehen.
Allerdings interessieren die Histogramme auch nur für einzelne Kanäle, d.h. die Eingangsgröße
in hist ist eine Matrix.
§ 1195 Auf eine Matrix angewendet liefert hist für jede Spalte der Matrix ein Histogramm.
Um das Histogramm aller Helligkeitswerte der Matrix zu erhalten, formen wir diese mit
reshape um:
>> LP1=reshape(P1(:,:),[],1);
Auf diese Weise wird aus der 850 × 1100-Matrix ein 935000 × 1 Vektor erzeugt. Der Datentyp
wird dabei nicht verändert, d.h. die Elemente des Vektors sind weiterhin vom Typ uint8.
§ 1196 Zum weiteren Rechnen ist uint8 (ebenso wie auch uint16) aus zwei Gründen unpraktisch:
• MatLab kann intern nicht jedes Format mit jedem anderen kombinieren. So lässt sich z.B.
hist nicht auf uint8-Größen anwenden.
c○ M.-B. Kallenrode 2. Juli 2008