zwei Arrays definieren 2D-Koordinaten, als Array-Indizes, in Matlab/Oktave

Question

Ich habe ein 2D-Array, nennen Sie es "A". Ich habe zwei andere 2D-Arrays, nenne sie 'ix' und 'iy'. Ich möchte ein Ausgabearray erstellen, dessen Elemente die Elemente von A an den Indexpaaren sind, die von x_idx und y_idx bereitgestellt werden. Ich kann mit einer Schleife dazu wie folgt vor:

for i=1:nx 
    for j=1:ny 
     output(i,j) = A(ix(i,j),iy(i,j)); 
    end 
end 

Wie kann ich dies tun, ohne die Schleife? Wenn ich output = A (ix, iy), bekomme ich den Wert von A über den gesamten Bereich von (ix) X (iy).

Danke, Jason

Answer 1

Ein schneller Weg ist lineare Indizierung direkt SUB2IND ohne Aufruf zu verwenden:

output = A(size(A,1)*(iy-1) + ix) 

... denken Sie an die Matrix A als 1D-Array (spaltenweise Reihenfolge)

Answer 2

Dies ist die einzeilige Methode, die für große Matrizen

reshape(diag(A(ix(:),iy(:))),[ny nx]) 

Eine klarere und effizientere Methode wäre, sub2ind zu verwenden nicht sehr effizient ist. Ich habe für Situationen igitt Kommentar aufgenommen (wie Sie), wenn ix und iy die gleiche Anzahl von Elementen haben:

newA = A(sub2ind(size(A),ix,iy)); 

Auch nicht zu verwechseln x und y für i und j in Notation - j und x allgemein repräsentieren Spalten und i und y repräsentieren Zeilen.