maximum d'une matrice

Bjr à tous, je ss nouveau dans le forum, j'ai un petit problème dans l'analyse de mes données avec matlab, en effet, je ss entrain de rechercher le maximum d'une
matrice, cependant il y'a des matrices qui ont des valeurs égales et me donne un maximum peu satisfaisant, j'utilise la fonction max(max(matrice)), mais je trouve
un maximum qui n'est pas maximum.
le probleme est le suivant je dois trouver des peaks d'une surface rugueuse, donc le maximum doit etre le point cent des 9 points de la matrice, un maximum strictement superieur à tous les 8 points autour de lui.
merci de votre aide, car je ss vraiment amateur dans l'analyse des données

Bonjour,

Personnellement je n'ai pas bien saisi ta démarche. Pourrais tu nous montrer ce que tu as fait et nous indiquer là où le bât blesse? Dans tout les cas si tu cherches des pics, la fonction max(max) ne t'indiquera jamais plus que le maximum de la matrice?

Gooby

Merci de l'attention porté à mon problème, voici mon programme entier matlab:
je cherche des peaks sur surface rugueuse pour ce faire je crée une matrice 3X3 = 9 points, le point milieu ou central doit etre le maximum absolu des points avoisinants, dans le cas où la valeur de valeur maximum est égale à la valeur de l'un des 8 points avoisinants, alors le programme continue la recherche, la fonction max(max(matrice)) ne m'aide pas car ne me donne pas le point max absolu. ainsi je cherche comment retrouver le maximum central absolu d'une matrice 3X3, car je continuerais avec les matrices 5X5, 7X7, 9X9, 11X11.

programme entier:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72

clear all; load('Cil3_50x_onepeak'); s1=length(Zm111(1,:)); s2=length(Zm111(:,1)); k=0; w=1.4; nu1=0.41; nu2=0.41; E1=103e+9; E2=107e+9; Estar=((1-nu1^2)/E1+(1-nu2^2)/E2)^-1; Z1 = 0.3e-9; Zms = Zm111*0; Xms = Xm111; Yms = Ym111; dimframe=1; for j=dimframe:s1-2*dimframe i=dimframe; while i<s2-2*dimframe matrice = Zm111(i:i+2*dimframe,j:j+2*dimframe); massimo = max(max(matrice)); if massimo == Zm111(i+dimframe,j+dimframe) k=k+1; peaks(k,3)= Zm111(i+1,j+1); peaks(k,2)= Ym111(i+1,j+1); peaks(k,1)= Xm111(i+1,j+1); f = fittype('a*((x-x0)*cos(th)+(y-y0)*sin(th)).^2+b*(-(x-x0)*sin(th)+(y-y0)*cos(th)).^2+z0','problem',{'x0','y0'},'independent', {'x', 'y'},'dependent', 'z','coefficients', {'a','b','th','z0'}); [cfun,gof] = fit([reshape(Xm111(i:i+2*dimframe,j:j+2*dimframe),[],1),... reshape(Ym111(i:i+2*dimframe,j:j+2*dimframe),[],1)],... reshape(Zm111(i:i+2*dimframe,j:j+2*dimframe),[],1),f,... 'problem',{Xm111(i+dimframe,j+dimframe),Ym111(i+dimframe,j+dimframe)},... 'Start', [1,1,0,0],... 'Lower',[-Inf -Inf -pi/2 -Inf],... 'Upper',[0 0 0 Inf]); peaks(k,3)= Zm111(i+dimframe,j+dimframe); peaks(k,2)= Ym111(i+dimframe,j+dimframe); peaks(k,1)= Xm111(i+dimframe,j+dimframe); peaks(k,4)=cfun.a; peaks(k,5)=cfun.b; peaks(k,6)=cfun.th; peaks(k,7)=cfun.z0; Rmax = -real(1/((cfun.a+cfun.b)-sqrt((cfun.a)^2-(2*cfun.a*cfun.b)+(cfun.b)^2))); Rmin = -real(1/((cfun.a+cfun.b)+sqrt((cfun.a)^2-(2*cfun.a*cfun.b)+(cfun.b)^2))); peaks(k,8)= Rmax; peaks(k,9)= Rmin; muTmax =((16*Rmax*1e-6*w^2)/(9*Estar^2*Z1^3))^1/3; muTmin =((16*Rmin*1e-6*w^2)/(9*Estar^2*Z1^3))^1/3; peaks(k,10)= muTmax; peaks(k,11)= muTmin; rmse = gof.rmse; peaks(k,12)=rmse;   i= i+dimframe;   Ymatr = Ym111(i:i+2*dimframe, j:j+2*dimframe); Xmatr = Xm111(i:i+2*dimframe, j:j+2*dimframe); x = reshape(Xmatr, 1,[]); y = reshape(Ymatr, 1,[]); z = reshape(Zmatr, 1,[]); figure grid on hold on Zsphere=cfun(Xm111,Ym111); surf(Xm111,Ym111,Zsphere); surf(Xm111,Ym111,Zm111); plot3(x,y,z,'.','Markersize',15) zlim([-0.4 0.4]) xlim([min(x)-2.5 max(x)+2.5]) ylim([min(y)-2.5 max(y)+2.5]) caxis([-0.6 0.6]) shading interp end i=i+1; end j; end

---

Bonne journée Goody

Si je comprends bien , tu cherches les maximum locaux, c'est bien ça ?

Si oui, et si tu possèdes l'Image Processing Toolbox, tu peux te servir de la fonction nlfilter.

Par exemple :

Code:

---

1 2 3 4 5 6 7 8 9 10 11

fun = @(x) all(x(5)>x([1:4 6:9]));   A = rand(5);   B = nlfilter(A, [3,3], fun);   C = A.*B;   idx = B(:)==1;   maxval = A(idx);

---

Ce qui donne par exemple :

Code:

---

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

A =   0.9391 0.6482 0.3479 0.3308 0.7069 0.3013 0.0252 0.4460 0.8985 0.9995 0.2955 0.8422 0.0542 0.1182 0.2878 0.3329 0.5590 0.1771 0.9884 0.4145 0.4671 0.8541 0.6628 0.5400 0.4648     B =   1 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 1 0 0 1 0 0 0     C =   0.9391 0 0 0 0 0 0 0 0 0.9995 0 0.8422 0 0 0 0 0 0 0.9884 0 0 0.8541 0 0 0     maxval =   0.9391 0.8422 0.8541 0.9884 0.9995

---

merci de la reponse apportée à mon problème, essayant de verifier ton script, j'obtiens une solution peu satisfaisante:
voici le script:

Code:

---

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

>> fun = @(x)all(x(3)> x([1:4 6:9])); a = rand(3); b = nlfilter(a,[3,3],fun); c = a.*b; idx = b(:)== 1; maxval = a(idx);   et la solution: a = 0,1455 0,5797 0,8530 0,1360 0,5498 0,6220 0,8692 0,1449 0,3509   b= 0 0 0 0 1 0 0 0 0   c= 0 0 0 0 0,5498 0 0 0 0   idx= 0 0 0 0 1 0 0 0 0   maxval = 0,5498

---

cette solution ne me satisfait pas car dans la matrice "a", maxval n'est pas le maximum absolu, nous avons autour de cette valeur centrale, autres valeurs tels que 0,5797; 0,8530; 0,6220; 0,8692 supérieurs à maxval=0,5498, je cherche maxval absolument superieur à tous les points de la matrice et , maxval doit etre au centre des points de la matrice.

Merci Dut

Il faut comparer le 5ème élément du bloc 3x3 :

Code:

---

fun = @(x) all(x(5)> x([1:4 6:9]));

---

C'est l'indice linéaire de l'élément central :

1 4 7
2 5 8
3 6 9

# Merci infiniment Dut, ton aide a été fondamentale et capitale pour la suite de l'analyse de mes données, qu'à celà ne tienne j'ai essayé d'utiliser ton script pour la détection des maximums locaux, cependant, en utilisant des données reduites c'est-à-dire 10x10, je parviens à ressortir 49 peaks, mais une fois que j'utilise des données 640x480, l'analyse est très long et à la fin me donne des messages d'erreurs.
En outre, durant les calculs, s'affiche toujours ce message en rouge "Applying Neighborhood Operation".
N'est-il pas possible de reduire le temps de calcul de peaks pour des données importantes(640x480) , car je n'arrive pas avec des données assez importantes?
ce message en rouge, peut-on le cacher durant l'analyse des données?
Merci une fois encore Dut, excellente journée.

Bonjour,

Dans le même principe, avec la fonction colfilt :

Code:

---

1 2 3 4 5

N = ... center = (ceil(N/2)-1)*N + ceil(N/2); fun = @(x) all(bsxfun(@gt,x(center,:),x([1:center-1 center+1:end],:)), 1); B = colfilt(A, [N N], 'sliding', fun);

---

Peux tu me dire que représente N, center?
Comment trouver le maximum local de matrice?
Merci Winjerome.

La taille du bloc considéré (si N=3 => 3x3 etc.)

voilà ton code fonctionne à merveilles, cependant quand j'analyse des données par exemple 480x640 élements, matlab tourne infiniment sans me donner des resultats (sans faire des calculs).
Merci Winjerome.

Il se peut que ta machine ne soit pas suffisamment puissante pour effectuer le calcul rapidement

Augmente progressivement la taille de l'image à traiter, en notant à chaque fois le temps d'exécution.
Tu pourras alors déterminer si c'est bien un problème de puissance de calcul