Discussion :

[alexov]

Bonjour,

Je voudrais convertir les valeurs de mes axes du tracé d'une de mes images.

En effet ils sont en case distances et je voudrais les convertir en distances.


Voici mon code, avec A qui est une matrice dépendant d'un nombre k qui varie de 0 à 200.

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
 
fig=figure       
title('Profil Distance')                                  
 
set(gca,'nextplot','replacechildren') 
 
for k=1:200
    y=20*log10(abs(fftshift(fft(Videos{1,1}{1,1}(:,k)))));       
 
    [pks loc]= findpeaks(y,'sort', 'descend');            
    pic = 10.^((pks(1:4))/10); 
 
    G(k)=(pic*(loc(1:4)'))/(sum(pic));  
 
    if k==1
       A(k,:)=(Videos{1,1}{1,1}(:,k))+2300;
    else
        A(end+1,:)=(exp(-2*pi*1i*2*((G(k)-G(1))+2300)*t1casedist)/c*(B/to*(1:4600)*Pech))...
            .*(Videos{1,1}{1,1}(:,k)).';
 
    end
end
 
imagesc(20*log10(abs(fftshift(fftshift(A)))))

Mon axe des x varie de 0 à 4600. Je pensais écrire une ligne du style :

Code :

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

x-axis = (0:(4600*t1casedist))

[Jerome Briot]

La question posée n'est pas très claire.

Essaie peut être de mettre ce bout de code à la fin du tien :

Code :

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

1
2
xt = get(gca,'xtick');
set(gca,'xticklabel',num2str(xt(:)*t1casedist)

Si ce n'est pas ça, dis nous en plus sur ce que tu entends par "case distance" et ce que représente la variable t1casedist

[alexov]

En fait t1casedist n'est pas une variable, c'est une constante, qui transforme la valeur d'une case distance en distance, c'est à dire en mètre. En fait mes axes sont ici définies en case distance, et je voudrais qu'ils le soient en mètres, c'est à dire en distance...

J'ai essayé le bout de code que vous m'avez donné mais il ne marche pazs . Les valeurs sur les axes des x ne changent pas... Et j'ai vraiment essayé plusieurs solutions.

[Invité]

Les valeurs sur les axes des x ne changent pas...

Si j'ai bien compris, il te faut les rentrer dans IMAGESC:

Code :

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

imagesc(X,Y,20*log10(abs(fftshift(fftshift(A)))))

[Jerome Briot]

Il manque une parenthèse fermante dans le code :

Code :

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

set(gca,'xticklabel',num2str(xt(:)*t1casedist))

Sinon dans ton code, il y a des parenthèses inutiles :

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
fig=figure       
title('Profil Distance')                                  

set(gca,'nextplot','replacechildren') 

for k=1:200
    y=20*log10(abs(fftshift(fft(Videos{1,1}{1,1}(:,k)))));       

    [pks loc]= findpeaks(y,'sort', 'descend');            
    pic = 10.^((pks(1:4))/10); 

    G(k)=(pic*(loc(1:4)'))/(sum(pic));  

    if k==1
       A(k,:)=(Videos{1,1}{1,1}(:,k))+2300;
    else
        A(end+1,:)=(exp(-2*pi*1i*2*((G(k)-G(1))+2300)*t1casedist)/c*(B/to*(1:4600)*Pech))...
            .*(Videos{1,1}{1,1}(:,k)).';

    end
end

imagesc(20*log10(abs(fftshift(fftshift(A)))))

[alexov]

Voici donc la fin de mon code :

Code :

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

1
2
3
4
5
6
7
fig=figure 
title('Profil Distance') 
x = get(gca,'xtick');
x=0:4600
set(gca,'xticklabel',num2str(x(:)*t1casedist))
 
imagesc(x,y,(20*log10(abs(fftshift(fftshift(A))))))

J'ai intégré vos deux bouts, de code. Ceci me donne bien des changements dans les axes, mais pas ceux que je souhaite. Ainsi l'axe des x va jusqu'à 1, alors qu'il devrait aller jusqu'à 4600*t1casedist = 1.49 à peu près... Donc je ne comprends pas trop ce que le code en fait ^^...

[Invité]

Tu as juste à spécifier les min/max:

Code :

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

imagesc([xmin xmax],[ymin ymax],(20*log10(abs(fftshift(fftshift(A))))))

[alexov]

Merci beaucoup !! Pour l'axe des x c'est parfait ... J'ai, juste un petit problème pour l'axe des y, parce que maintenant il me les mets de haut en bas lol...

Merci pour votre aide

[Invité]

C'est normal, pour les images, l'axe y est inversé

Code :

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

set(gca,'YDir','normal')