Discussion :

[Krokro00]

Bonjour à tous,

Comme le titre l'indique, je cherche à afficher les valeurs issues d'un filtre canny sur une carte .

Pour l'instant je n'ai réussi qu'avec la fonction 'contour' mais le résultat n'est pas terrible (voir fichier joint). Je souhaiterai plutot que les valeurs 1 s'affichent avec une couleur et non un contour mais je n'arrive pas à utiliser la fonction pcolor (sachant que je veux garder l'info géographique X,Y)....

Voici une partie de mon code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
TSMe=load('Poly-TSMc-06-mask.txt');
    %TSMe=load('Poly-TSMc-11-mask.txt');
    [nblon nblat]=size(TSMe)
    lat=-1.9:((-1.3)-(-1.9))/(nblat-1):-1.3;
    lon=43.9:(43.2-43.9)/(nblon-1):43.2;
    [X,Y]=meshgrid(lat,lon);

   figure(3),clf %

    hold on
    pcolor(X,Y,TSMe);shading interp
    [C,h]=contour(X,Y,TSMe,[2 4 6 8 10 12 14 16],'k'); 
    T=adaptivethreshold(TSMe,32,0.00001,0);
    C2=edge(T,'canny',0.001)
    [C,h]=contour(X,Y,C2,[1],'r'); 
    caxis([0 20]);
    colormap(jet2);
    axis equal;
     hold off

Si quelqu'un sait comment faire?
Merci d'avance!

[Krokro00]

En fait je cherche à delimiter les contours d'un panache de fleuve (ce panache étant caractérisé par de forte concentration en sédiment et donc différenciable des eaux marines pauvres en sédiment).

J'ai utilisé une première fonction pour imposer un seuil puis j 'ai fait un filtre canny pour determiner la limite du panache (normalement la zone de plus fort gradient de concentration)

Le filtre Canny me donne une réponse de type binaire en fonction du type de seuillage que j'ai choisi (1 = ok par rapport au seuil, 0 = pas ok).

Ce que je veux faire c'est afficher ces valeurs 1 sur ma carte des concentrations en sédiment.

Le problème est que quand j'utilise pcolor (pcolor(X,Y,TSMe)) j'ai une erreur de type "value must be numeric" . En utilisant str2num , j'ai une erreur sur la taille des matrices....

Du coup je ne sais pas comment faire...
En espérant que j'ai été plus claire!

[Jerome Briot]

Je ne suis pas sûr de bien te suivre mais tu peux essayer ceci : [pcolor/contourf] Superposition de graphiques

[Krokro00]

Pour mieux illustrer mon problème je cherche à superposer les gradient que j'obtient après un filtre canny (figure exemple MES canny.jpg) à ma figure exemple MES .jpg qui elle est calée géographiquement.

[Krokro00]

Oups, les images !

[Jerome Briot]

Pourrais-tu attacher le fichier Poly-TSMc-06-mask.txt (en archive zip) si les données ne sont pas confidentielles ?

[Krokro00]

le voila!!

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
    TSMe=load('Poly-TSMc-06-mask.txt');

    [nblon nblat]=size(TSMe)
    lat=-1.9:((-1.3)-(-1.9))/(nblat-1):-1.3;
    lon=43.9:(43.2-43.9)/(nblon-1):43.2;
    [X,Y]=meshgrid(lat,lon);

   figure(3),clf %

    hold on

    T=adaptivethreshold(TSMe,32,0.00001,0);
    C2=edge(double(T),'canny',0.001);
    TSMe(C2) = max(TSMe(:));
    pcolor(X,Y,TSMe);shading interp
    caxis([0 20]);

    colormap(jet2);
    axis equal;

	hold off

Et au passage :

Code :

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

1
2
3
    [nblon nblat]=size(TSMe)
    lat=-1.9:((-1.3)-(-1.9))/(nblat-1):-1.3;
    lon=43.9:(43.2-43.9)/(nblon-1):43.2;

devient :

Code :

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

1
2
3
[nblon nblat]=size(TSMe)
lat = linspace(-1.9,-1.3,nblat);
lon = linspace(43.9,43.2,nblon);

[Krokro00]

Youhou !!! J'ai testé rapido, ca à l'air de marcher!!
Je reviendrai par contre demain avec quelque questions sur tes modifs!

Et merci pour le "linspace", je ne connaissait pas!!

bonne soirée

[Jerome Briot]

Une autre possibilité :

Code :

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

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
    TSMe=load('Poly-TSMc-06-mask.txt');
 
    [nblon nblat]=size(TSMe)
    lat = linspace(-1.9,-1.3,nblat);
    lon = linspace(43.9,43.2,nblon);
    [X,Y]=meshgrid(lat,lon);
 
    figure
    pcolor(X,Y,TSMe);shading interp
 
    colormap(jet2);
    caxis([0 20]);
    axis equal;
 
    T=adaptivethreshold(TSMe,32,0.00001,0);
    C2=edge(double(T),'canny',0.001);
 
    img = cat(3,C2(end:-1:1,:),zeros(nblon,nblat,2));
 
 
    hold on
 
    image('xdata',[min(X(:)) max(X(:))],'ydata',[min(Y(:)) max(Y(:))],...
        'cdata',img,...
        'alphadata',C2(end:-1:1,:))

Les contours sont plus nets et leur couleur est plus facile à gérer (valeur RGB dans la variable img) et indépendante de la palette de couleurs jet2

[Krokro00]

Bonjour Dut et le forum,

Merci beaucoup pour ton dernier post. Ca marche très bien. Il y a par contre des formulations que je ne comprends pas vraiment ....

A quoi sert le double dans

Code :

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

 C2=edge(double(T),'canny',0.001);

a quoi sert le concatenage et que fait exactement la formule :

Code :

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

img = cat(3,C2(end:-1:1,:),zeros(nblon,nblat,2));

Merci beaucoup

[Jerome Briot]

A quoi sert le double dans

Code :

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

 C2=edge(double(T),'canny',0.001);

A priori, la fonction ADAPTIVETHRESHOLD retourne une matrice T logique (0 ou 1).
La fonction EDGE n'accepte pas ce genre de matrice, il faut donc la transformer en Double

a quoi sert le concatenage et que fait exactement la formule :

Code :

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

img = cat(3,C2(end:-1:1,:),zeros(nblon,nblat,2));

A créer une image RGB (MxNx3) dont les couleurs sont indépendantes de la pallette de couleur courante (colormap).

Ici, pour avoir les "lignes" en rouge (valeur RGB [1 0 0]), il suffit de concaténer la matrice C2 et deux matrices de même dimension avec des 0.

Il est nécessaire d'inverser les lignes de C2 car la direction des axes pour la représentation des matrices et des images est inversée... (voir la fonction AXIS et ses paramètre ij et xy)

[Krokro00]

Ok!! Merci énormement pour ton aide!