Avoid loops with mapping toolbox "setltln"
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
Robert Jenkins
il 4 Ago 2015
Risposto: Sean de Wolski
il 4 Ago 2015
Hi all, I'm looking for a clever way to optimize the following chunk of code. The function is to get a grid of lat, lon coordinates using the reference matrix R, and the indices of Z.
However, with such huge dimensions, this nested for loop is painfully slow! How can I avoid for loops in this case?
[Z, R] = arcgridread('biloxi_ms.asc');
tempLAT=zeros(9721,10801);
tempLON=zeros(9721,10801);
for i =1:9721
for j = 1:10801
[tempLAT(i,j), tempLON(i,j)] = setltln(Z, R, i, j);
end
end
Best, Rob
0 Commenti
Risposta accettata
Sean de Wolski
il 4 Ago 2015
You could build a grid and call setltln on all of it. This will consume a bit more memory since you'll need to store ii, jj. A hybrid approach might be to run a loop over smaller meshgrids.
[ii, jj] = meshgrid(1:9271,1:10801)
[tempLAT, tempLON]= setltln(Z, R, ii, ij);
Also, parallel computing could help here if you have the parallel computing toolbox. Simply, turn the outer loop into a parfor.
0 Commenti
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Execution Speed in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!