Discussion :

[tobalito]

Bonjour,

Je cherche à extraire les valeurs d'un volume (

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

map

) autour du voxel (

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

row,column,tube

), étant donné un voisinage de côté K (

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

neighbour

). Mon code actuel est brutal :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
function values=collectNeighbourValue3D(map,row,column,tube,neighbour)
 
limitLoop=(neighbour-1)/2;
l=1;
values=zeros(neighbour^3,1);
for i=(-limitLoop:1:limitLoop)
    for j=(-limitLoop:1:limitLoop),
        for k=(-limitLoop:1:limitLoop),
        values(l)=map(row+i,column+j,tube+k);
        l=l+1;
        end
    end
end

Je possede une liste de voxels (1 millions) dont je desire extraire les valeurs des voisins, et mon petit doigt me dit qu'il serait possible de faire cela sans boucle pour le rendre plus rapide.
Ma recherche sur divers forums n'a pas été fructueuse, bien que ce probléme paraisse assez banal.

Merci d'avance pour vos idées.

[Jean Dumoncel]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
limitLoop=(neighbour-1)/2;
values=map(row-limitLoop:row+limitLoop,column-limitLoop:column+limitLoop,tube-limitLoop:tube+limitLoop);

Il faudrait peut-être rajouter une condition sur neighbour pour vérifier que c'est un nombre pair.

[tobalito]

Simple et efficace.
Merci