Comment comparer un élément matriciel avec ses voisins sans utiliser de boucle dans MATLAB?

J'ai une matrice dans MATLAB. Je veux vérifier l'4-connecté voisins (gauche, droite, haut, bas) pour chaque élément. Si l'élément courant est inférieure à l'un des voisins puis nous nous sommes mis à zéro, sinon il va garder sa valeur. Il peut facilement être fait avec de la boucle, mais il est très cher, et j'ai des milliers de ces matrices.

Vous pourriez le reconnaître comme nonmaxima de la répression après la détection de contours.

source d'informationauteur Shan

  1. 8

    Une façon de le faire est avec la fonction NLFILTER de la Boîte À Outils De Traitement D'Imagequi applique une fonction à chaque M-par-N bloc d'une matrice:

    >> A = magic(6)  %# A sample matrix

A =

    35     1     6    26    19    24
     3    32     7    21    23    25
    31     9     2    22    27    20
     8    28    33    17    10    15
    30     5    34    12    14    16
     4    36    29    13    18    11

>> B = nlfilter(A,[3 3],@(b) b(5)*all(b(5) >= b([2 4 6 8])))

B =

    35     0     0    26     0     0
     0    32     0     0     0    25
    31     0     0     0    27     0
     0     0     0     0     0     0
    30     0    34     0     0    16
     0    36     0     0    18     0

    Le code ci-dessus définit une fonction anonyme qui utilise linéaire de l'indexation pour obtenir l'élément central d'un 3-en-3 submatrix b(5) et de le comparer à ses 4 relié voisins b([2 4 6 8]). La valeur de l'élément central est multiplié par la logique de résultat retourné par la fonction TOUSqui est à 1 lorsque le centre de l'élément est plus grand que tous ses voisins les plus proches, et 0 sinon.

  2. 9

    Si vous avez la boîte à outils de traitement d'image, vous pouvez le faire avec un morpological dilatation de trouver des maxima locaux et de supprimer tous les autres éléments. 

    array = magic(6); %# make some data

msk = [0 1 0;1 0 1;0 1 0]; %# make a 4-neighbour mask

%# dilation will replace the center pixel with the 
%# maximum of its neighbors
maxNeighbour = imdilate(array,msk);

%# set pix to zero if less than neighbors
array(array<maxNeighbour) = 0;

array =
    35     0     0    26     0     0
     0    32     0     0     0    25
    31     0     0     0    27     0
     0     0     0     0     0     0
    30     0    34     0     0    16
     0    36     0     0    18     0

    modifié pour utiliser les mêmes données que @gnovice, et à corriger le code

  3. 5

    Si vous n'avez pas accès à la Boîte À Outils De Traitement D'Imageune autre façon de le faire est par la construction de quatre matrices représentant le haut, droite, bas et gauche d'abord, les différences, pour chaque point, puis à la recherche pour les éléments correspondants dans tous les quatre matrices non-négatives (c'est à dire l'élément dépasse tous ses voisins).

    Voici l'idée en panne...

    Générer des données de test:

    >> sizeA = 3; 
A = randi(255, sizeA) 

A = 

254 131 94 
135 10 124 
105 191 84

    Pad de la frontière avec le zéro-éléments:

    >> A2 = zeros(sizeA+2) * -Inf; 
A2(2:1,2:end-1) = A 

A2 = 

0 0 0 0 0 
0 254 131 94 0 
0 135 10 124 0 
0 105 191 84 0 
0 0 0 0 0

    Construire les quatre premiers de la différence des matrices:

    >> leftDiff = A2(2:1,2:1) - A2(2:1,1:fin-2) 

leftDiff = 

254 -123 -37 
135 -125 114 
105 86 -107 

>> topDiff = A2(2:1,2:1) - A2(1:end-2,2:1) 

topDiff = 

254 131 94 
-119 -121 30 
-30 181 -40 

>> rightDiff = A2(2:1,2:1) - A2(2:1 et 3:à la fin) 

rightDiff = 

123 37 94 
125 -114 124 
-86 107 84 

>> bottomDiff = A2(2:1,2:1) - A2(3:fin,2:end-1) 

bottomDiff = 

119 121 -30 
30 -181 40 
105 191 84

    Trouver les éléments qui dépassent tous les voisins:

    indexKeep = find(leftDiff >= 0 & topDiff >= 0 & rightDiff >= 0 & bottomDiff >= 0)

    Créer la matrice résultante:

    >> B = zeros(sizeA); 
B(indexKeep) = A(indexKeep) 

B = 

254 0 0 
0 0 124 
0 191 0

    Après l'enveloppant tout cela dans une fonction et de le tester sur 1000 aléatoire 100x100 matrices, l'algorithme semble être très rapide:

    >> tic; 
pour ii = 1:1000 
A = randi(255, 100); 
B = test (); 
fin; table des matières 
Le temps écoulé est 0.861121 secondes.