MATLAB: Auszug submatrix mit logischen Indizierung

Ich bin auf der Suche nach einer eleganten Lösung für dieses sehr einfache problem in MATLAB. Angenommen ich habe eine matrix 

>> M = magic(5)

M =

    17    24     1     8    15
    23     5     7    14    16
     4     6    13    20    22
    10    12    19    21     3
    11    18    25     2     9

und eine logische variable der form

I =

     0     0     0     0     0
     0     1     1     0     0
     0     1     1     0     0
     0     0     0     0     0
     0     0     0     0     0

Wenn ich versuche, die Elemente abzurufen, der M zugeordnet 1 Werte in I bekomme ich einen Spalten-Vektor 

>> M(I)

ans =

     5
     6
     7
    13

Was wäre der einfachste Weg, um die matrix [5 7 ; 6 13] von diesem logischen Indizierung?

Wenn ich weiß, dass die Form der nicht-null-Elemente von I, ich kann mit einer Umgestaltung nach der Indizierung, aber das ist keine Allgemeine Fall.

Auch, ich bin mir bewusst, dass das Standardverhalten für diese Art der Indizierung in MATLAB erzwingt Konsistenz in Bezug auf den Fall, in dem nicht-null-Werte in I nicht in form einer matrix, aber ich Frage mich, ob es gibt eine einfache Lösung für diesen speziellen Fall.

InformationsquelleAutor foglerit | 2012-07-10
  1. 11

    Dies ist ein Weg, dies zu tun. Es wird davon ausgegangen, dass alle Zeilen von I haben die gleiche Anzahl von Einsen. Es ist auch davon auszugehen, dass alle Spalten der I haben die gleiche Anzahl haben, weil Submatrix müssen rechteckig sein.

    %# Define the example data.

M = magic(5);
I = zeros(5);
I(2:3, 2:3) = 1;

%# Create the Submatrix.

Submatrix = reshape(M(find(I)), max(sum(I)), max(sum(I')));
    InformationsquelleAutor nrz
  2. 3

    Hier ist eine sehr einfache Lösung:

    T = I(any(I'),any(I));
T(:) = M(I);
    • Genau das, was ich im Begriff war, zu beantworten 🙂
    • Das funktioniert nicht immer, z.B. wenn Sie nicht symmetrisch ist.
    • Wenn du meinst die matrix I können nicht gemacht werden, nach dem Zufallsprinzip, das ist wahr, aber mit Blick auf die Frage, die ich denke, das ist eine sinnvolle Einschränkung. Wenn Sie etwas anderes meinen: ich sehe nicht, wie symetry ist hier relevant, hättest du vielleicht ein Beispiel, wo die Lösung nicht? (Matrix-Input, output, erwartete Ausgabe).
    • OK, ich war nicht präzise genug in meinem Kommentar: Wenn die Auswahl-matrix I nicht symmetrisch ist, erhalten Sie einige Probleme mit der form T. In der Tat, ich glaube, Sie können korrigieren Sie die Antwort, indem einfach die erste Zeile um: T = I(any(I'), any(I));. Versuchen Sie Ihre Lösung mit der in diesem Beispiel matrix-input: I = zeros(5);I(1:4, 2:3) = 1;
    • Es ist schön, die Begeisterung zu sehen, aber beachten Sie die Allgemeinen Nutzungsbedingungen: Wenn Sie denken, dass ein Beitrag gut ist, für ihn abstimmen. Gesprächig Kommentare machen es schwieriger für die zukünftigen Leser zum extrahieren der relevanten Informationen aus einem thread.
    InformationsquelleAutor Dennis Jaheruddin
  3. 2
    M = magic(5);
I = [ ... ];

ind = find(I); %# find indices of ones in I
[y1, x1] = ind2sub(size(M), ind(1));   %# get top-left position
[y2, x2] = ind2sub(size(M), ind(end)); %# get bottom-right position
O = M(y1:y2, x1:x2); %# copy submatrix
    InformationsquelleAutor kay
Schreibe einen Kommentar