Discussion :

[Newenda]

Bonjour,

Voilà, j'essaye de faire un spectrogramme. Pour cela j'ai un jeu de données

Code :

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

1
2
3
4
5
6
7
temps  bande1     bande2     bande3     bande4
1        45          25          10        3 
2        58          25          5         1     
3        98          21          8         0     
4        25          18          7         0     
5        87          24          7         0    
6        47          25          4         2

bande1 correpond à la puissance spectrale entre 10 et 20Hz, la bande2 entre 20 et 35Hz, la bande 3 entre 35 et 85Hz, la bande4 entre 85 et 200Hz.

Je réalise donc une grille :

Code :

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

1
2
Xgrid = [1:6]  %%>> correspondant au temps
Ygrid = [10, 20, 35, 85, 200] %%>> correspondant au bandes spectrales.

J'aimerais maintenant faire correspondre mes données aux différentes cases de la grille, comme ça:



et remplacer mes données par des couleurs (comme colorbar, du bleu pour les valeurs les plus faibles, au rouge pour les plus grandes) pour finalement avoir une sorte de spectrogramme.

Donc j'aimerais bien un peu d'aide pour :

-faire correspondre mes données à la grille
-transformés ces données en couleur avec une échelle de couleur colorbar.

Merci bien !

[Jean Dumoncel]

Salut,

Une solution avec meshgrid et surf :
Il faut tout d'abord créer la grille, c'est à dire les valeurs de x, y et z pour chaque point de ta grille. Pour cela il faut utiliser meshgrid pour créer x et y puis tu remplis z (qui correspond à la puissance spectrale) dans le même format que x et y.

Ensuite tu utilises surf pour afficher le résultat et tu modifies view pour avoir une vue de dessus. Les couleurs par défaut seront du bleu (+faible valeur) au rouge.

Sinon je crois me souvenir d'une toolbox gratuite sur le net qui faisait beaucoup de choses en traitement du signal, mais son nom m'échappe pour le moment... Si je la retrouve, je reposte!

[Jerome Briot]

Une solution avec meshgrid et surf :

PCOLOR ou même IMAGESC peuvent également être utilisées ici

[Newenda]

Salut,

Une solution avec meshgrid et surf :

Effectivement :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
[X,Y]=meshgrid(1:6,[25 35 85 200]) 
 
valeurT1 = [45 25 10 3] ;
valeurT2 = [58 25 5 1] ;
valeurT3 = [98 21 8 0] ;
valeurT4 = [25 18 7 0] ;
valeurT5 = [87 24 7 0] ;
 
valeurT6 = [47 25 4 2] ;
Z=[valeurT1 ; valeurT2 ; valeurT3 ; valeurT4  ;valeurT5 ;valeurT6]' ;
 
surf(X,Y,Z)
view(2)

retourne bien,



Seulement et comme toujours, je n'ai pas une trentaine de données mais quelques millions...le fait de faire du SURF n'engendre t il pas des calculs inutiles pour une visualisation 2D?

Le 2ème problème est que dans mon cas, j'ai.4 bandes de fréquence, il faut donc 5 valeurs et non 4 pour les définir. Ici j'ai du zapper le 10 Hz pour la grille Y.
Y'a t il moyen de corriger cela?


En tout cas merci beaucoup.

[Caro-Line]

Comme dit par Dut tu as le même résultat avec pcolor :

Code :

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

pcolor(Xgrid,Ygrid(2:end),Z)

[Jean Dumoncel]

Quelques millions !!!

Le mieux serait d'essayer...
Sinon tu peux essayer pcolor conseillé par Dut, mais c'est juste un dérivé de surf, sauf que la vue sera au préalable ajuster donc ça permettra surement de gagner en efficacité...

Après je ne pense pas que l'on gagne quelque chose avec imagesc car la seule manière que je vois de l'utiliser ici est de remplir une matrice en fonction des indices de tes données, et donc d'augmenter considérablement le nombre de point, mais peut-être que c''est pas la meilleure manière de faire...

[Jean Dumoncel]

Le 2ème problème est que dans mon cas, j'ai.4 bandes de fréquence, il faut donc 5 valeurs et non 4 pour les définir. Ici j'ai du zapper le 10 Hz pour la grille Y.
Y'a t il moyen de corriger cela?

J'avais remarqué cela aussi, il manque une ligne et une colonne...
Tu peux les rajouter en x et y :

Code :

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

[X,Y]=meshgrid(1:7,[10 25 35 85 200]);

Et rajouter une ligne de 0 et une colonne de zeros dans z :

Code :

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

1
2
3
4
5
6
7
8
9
valeurT1 = [45 25 10 3 0] ;
valeurT2 = [58 25 5 1 0] ;
valeurT3 = [98 21 8 0 0] ;
valeurT4 = [25 18 7 0 0] ;
valeurT5 = [87 24 7 0 0] ;
valeurT6 = [47 25 4 2 0] ;
valeurT7 = [0 0 0 0 0] ;
Z=[valeurT1 ; valeurT2 ; valeurT3 ; valeurT4  ;valeurT5 ;valeurT6;valeurT7]' ;
pcolor(X,Y,Z)

[Newenda]

Effectivement PCOLOR fait la même chose et marche bien.

Merci à vous pour ces réponses variées, ça m'a permis d'en apprendre plus sur toutes ces fonctions.