Kombinieren Matlab Zelle Elemente zu einem einzigen Vektor

Question

Ich habe 3 480 Zeilen x 640 Spaltenbilder. Jedes Pixel enthält einen Intensitätswert.

I lese sie in Matlab mit diesem:

object_intensities = cell(1,6); 
    for pat=1:3 
     file1 = sprintf('image%d.pgm',pat-1); 
     dummy = double(imread(file1)); 
     object_intensities{pat} = dummy; 
    end 

Nun zu pat = 1 bis 3, object_intensities eine Zellenanordnung derart ist, dass jede Zelle Element object_intensities {pat} jedes Bild in einer 480x640-Matrix Array enthält .

Nun, was ich brauche, ist eine einzige Struktur, die 480 Zeilen x 640 Säule so ist, dass jede Zelle Element die drei Intensitäten jedes Bildpixel in einem einzigen Vektor wie folgt enthält:

structure{row 1, column 1} = [image1_intensity @ row1 column1 image2_intensity @ row1 column1 image3_intensity @ row1 column1] 
structure{row 1, column2} = [image1_intensity @ row1 column2 image2_intensity @ row1 column2 image3_intensity @ row1 column2] 
. 
. 
structure{row 480, column 640} = [image1_intensity @ row480 column640 image2_intensity @ row480 column640 image3_intensity @ row480 column640] 

Jede Art und Weise zu tun, dies ohne Schleifen, sondern durch Vektorisierung?

Answer 1

Die Datenstruktur, die Sie wollen, fragt nach Ärger. Verwenden Sie stattdessen eine 3D-Matrix, in der die dritte Dimension Ihre Intensitäten (in der Reihenfolge) abbildet. Es wird viel einfacher sein, mit der ganzen Linie zu arbeiten.

Ich schlage vor, Sie es wie folgt erstellen:

object_intensities = zeros(480,640,3); 
for pat=1:3 
    file1 = sprintf('image%d.pgm',pat-1); 
    object_intensities(:,:,pat) = rgb2gray(imread(file1)); 
end 

Jetzt hätte, was object_intensities{i,j} auf Ihrer Struktur gewesen ist nur object_intensities(i,j,:) und wenn Sie wirklich benötigen Sie squeeze oder permute oder (:)' nennen kann es sich um eine Reihe zu machen Vektor. Aber denken Sie nur, wie viel einfacher es jetzt so etwas wie das Finden der durchschnittlichen Intensität pro Pixel zu tun ist,

mean(object_intensities,3) 

verglichen mit dem Zellenfeld

cellfun(@mean,object_intensities,'uni',0) %// Just a loop in disguise 

Und was ist, wenn Sie die durchschnittliche pro Zeile

wollte

mean(object_intensities,2) 

Versuchen Sie dies für das Zellenfeld und Sie werden sehen, warum Sie es vermeiden wollen.

Aber wenn Sie wirklich wollen, um Ihre Struktur wirklich, dann meinen Code wie oben verwenden und mat2cell auf der 3D-Matrix nennen:

mat2cell(object_intensities, ones(size(object_intensities,1),1), ones(1,size(object_intensities,2)),3)