Detection d'un max sur un groupe de données dans un vecteur.

Version imprimable

Bonjours,

j'ai un vecteur V de taille 50 de se type:

[0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0]

j'aimerais trouver le maximum de chaques groupes de donnée.

Par groupe de donnée j'entend par la (dans cet exemple) :
(3 6 4), (3 2 6 7 6) et (7 8 2 0 6).

Je n'ai trouvé pour l'instant comme méthode que celle de "fenêtrage" càd:
Code:

1 2 3 4 5 6 taille_fenetre=10; fenetrage = reshape(V,taille_fenetre,[]); [~,indice] = max(fenetrage); indice = indice + (0:length(indice)-1)*taille_fenetre; resultat=zeros(1,length(V)); resultat(indice)=V(indice);
fenetrage.':
Code:

1 2 3 4 5 [0 0 0 3 6 4 0 0 0 0; 0 0 0 0 0 0 0 3 2 6 ; 7 6 0 0 0 0 0 0 0 0 ; 0 0 7 8 2 0 6 0 0 0 ; 0 0 0 0 0 0 0 0 0 0 ]
le problème est qu'il me renvoi les indices de 5 maximum (5 20 21 34 41) alors que je n'est que 3 groupes de donnée, du a la coupure du 2e groupe de donnée et aussi du a l'absence de groupe de donnée dans la 5ème fenêtre.

J'aimerais donc comparer les maximums de ce genre de tableau (transposé biensurs):
Code:

1 2 3 4 [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 ; 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 ; 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 ; 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 ]
Mais là je bloque sans boucle for, déjà rien que pour crée se tableau.
Mon algorithme s'appliquerait sur des grands jeu de donnée et for serait trop long...

Donc voila si vous avez des idées a me suggérer pour se problème ou même une autre méthode de détection de maximum je suis preneur :D

20/07/2012, 16h52
baalk42
ok je suis bidon...:oops:

pour le tableau(transposé):
Code:

1 2 3 4 [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 ; 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 ; 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 ; 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 ]
c'est aussi simple que ça:

Code:

fenetrage2=[fenetrage(:,1:end-1);fenetrage(:,2:end)];

maintenant il faut que je trouve comment avoir resultat.

ok le fait de parlé de mon problème me permet de mieu le posé dans mon esprit.

voila le code que je pense être juste:
Code:

1 2 3 4 5 6 7 8 9 10 taille_fenetre=10; fenetrage = reshape(V,taille_fenetre,[]); fenetrage2=[fenetrage(:,1:end-1);fenetrage(:,2:end)]; [~,indice] = max(fenetrage2); indice = indice + (0:length(indice)-1)*taille_fenetre; indice = union(indice,indice); %supprime les doublets indice(cellarray(@isnan, indice)) = []; %supprime les nan resultat=zeros(1,length(V)); resultat(indice)=V(indice);
il me reste un petit problème...

pour un groupe de donnée 4 8 6 5 7 coupé de cet façon 4 8 6 | 5 7 par reshape, j'ai peur qu'il me trouve par la suite 2 max pour un groupe de donnée...

Bonjour,

Je ne sais pas si j'ai bien compris ton problème mais une autre idée serait de faire une labelisation en composante connexe:
Code:

1 2 3 4 5 6 7 8 l = [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0]; bw = bwlabel(l); resultat = zeros(1,max(bw)); for k=1:max(bw) resultat(k) = max(l.*(bw==k)); end
Par contre je ne comprend pas pourquoi 7 8 2 0 6 n'ai q'un groupe?
Peux-tu mieux expliquer ta définition du groupe.

Rque:
Si il faut mettre un critère de proximité entre les groupe.

23/07/2012, 07h55
baalk42

Bonjour,

un groupe de données est un groupe de valeur dont les 2 extrémités sont non nul et les groupes sont séparer par plusieurs 0 (séparé par n 0 pour n>=2).

J'espère avoir mieux défini ce qu'été mes groupes de donnée, maintenant je vais étudier ton code, merci bien :) .

Bonjour,

Une solution utilisant strfind :

Code:

1
2
3
4
5
6
7
8
9
10
11
V = [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0];
n = 2; % nombre de zéros minimum séparant les données
idx = [zeros(1,n) V zeros(1,n)]~=0;
debut = strfind(idx,[false(1,n) true]); % indices marquant le début de chaque groupe
fin = strfind(idx,[true false(1,n)])-n; % indices marquant la fin de chaque groupe
 
data = arrayfun(@(from, to) V(from:to) , debut,fin,'UniformOutput',false); % groupes
data{:}
 
maxima = arrayfun(@(from, to) max(V(from:to)) , debut,fin) % max de chaque groupe
% ou : cellfun(@(d) max(d) , data)

Citation:

Envoyé par baalk42

Bonjour,

un groupe de données est un groupe de valeur dont les 2 extrémités sont non nul et les groupes sont séparer par plusieurs 0 (séparé par n 0 pour n>=2).

J'espère avoir mieux défini ce qu'été mes groupes de donnée, maintenant je vais étudier ton code, merci bien :) .

D'accord dans ce cas un petit changement:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 n = 2; l = [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0]; tic s = imclose(l>0, strel('line', 2*n+1, 1) ); bw = bwlabel(s); resultat = zeros(1,max(bw)); for k=1:max(bw) resultat(k) = max(l.*(bw==k)); end toc resultat
Marche très bien ici. Les fonctions utilisées sont en générale pour le traitement d'image. Mais je remarque ici que le temps d'exécution reste court.

A comparer à celle de Winjerome ( qui à mon sens reste plus adapté dans votre cas):
Pour un jeu de données de 500000:

Code:

V = rand(1,50000).*(rand(1,50000)<0.3);

Winjerome: Elapsed time is 0.108319 seconds.
Gakusei: Elapsed time is 0.552182 seconds.

23/07/2012, 14h58
baalk42

Un grand merci a tous les 2 pour vos solutions.

Par contre Gakusei mon MATLAB ne reconnait les fonctions imclose, strel et bwlabel. Je n'ai pas Image Processing Toolbox sur ma version MATLAB...

@Winjerome:

T'as fonction règle mon dernier petit problème (c'est super !) :

Citation:

pour un groupe de donnée 4 8 6 5 7 coupé de cet façon 4 8 6 | 5 7 par reshape, j'ai peur qu'il me trouve par la suite 2 max pour un groupe de donnée...

Par contre n'étant sur matlab que depuis peu, je n'arrive pas à comprendre et à trouver une explication sur l'help de Matlab ou sur internet à ce terme:

Code:

@(from, to)

Aurais tu un liens, un FAQ ou ton explication qui me permettais de comprendre ?

Et enfin je tien à m'excuser se n'est pas le max de chaque groupe que je cherche mais les indices du max de chaque groupe :oops:.
Donc je me permet de modifier ton code:
Code:

1 2 3 4 5 6 7 8 9 10 11 V = [0 0 0 3 6 4 0 0 0 0 0 0 0 0 0 0 0 3 2 6 7 6 0 0 0 0 0 0 0 0 0 0 7 8 2 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0]; n = 2; % nombre de zéros minimum séparant les données idx = [zeros(1,n) V zeros(1,n)]~=0; debut = strfind(idx,[false(1,n) true]); % indices marquant le début de chaque groupe fin = strfind(idx,[true false(1,n)])-n; % indices marquant la fin de chaque groupe data = arrayfun(@(from, to) V(from:to) , debut,fin,'UniformOutput',false); % groupes data{:} [~,idxDMax] = arrayfun(@(from, to) max(V(from:to)) , debut,fin) ;% max de chaque groupe idxDMax=idxDMax+debut-1;

23/07/2012, 18h35
Invité

Cette syntaxe permet de créer ce que l'on appelle une fonction anonyme (ou anonymous Function en anglais). (from , to) représentent ses entrées (qui sont tour à tour tirées de debut et fin fournis à la fonction arrayfun), et V(from:to) ce qu'elle retourne.
Tu pourras trouver plus de détails dans ce topic de la FAQ.