Discussion :

[KAKATTAH]

Bonsoir les amis.

J´ai une serie de courbe que je genère dans une boucle ''for'' avec ''contourf'' mais il y´en une qui n´apparait pas totalement ie elle est reproduite en une fraction; pourtant, avec ''contour'' je vois bien les isolignes des autres valeurs qui ne sont pas prises en compte par ''contourf''. Et, même encore pour cette figure inachevée, quand j´utilise le ''pcolor'', elle est totalement reproduite. Comment faire pour remedier à cela car j´aimerai utiliser du contourf car elle est mieux adapté pour ces figures.

Merci d´avance

[Jerome Briot]

Montre le code que tu utilises et également un aperçu du graphique que tu obtiens.

[KAKATTAH]

Merci de votre attention;

Voici une partie du code que j'utilise et les figures qui en découlent quand j'utilise les différentes fonctions: "pcolor", "contour" ou "contourf".

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
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
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
clear all; close all; clc
 
% chemin pour les donnes
DIR = '../donees/ascii/';
 
%chemin de lon lat
DIR1 = '../dados/xls/';
 
DIR2 = '../mat/';
 
load ([DIR2,'CFI_fisica_Atol'],'T2A','dens2A','sal2A','latA','lonA')
 
pA = 1:nA; % cest pour donner la taille
iA = 1:nA; % pour donner aussi dabord la taille
 
for i = 1:nA;
    pA(i) = max(depthA(:,i)); % profondeur max
    %pA(i) = 30; % profondeur choisi
    iA(i) = find(depthA(:,i) == pA(i));% determiner lindice de cette profondeur
    %T1fA(i) = T1A(iA(i),i);
end
 
 
 
lon1 = lonA([13 12 11 10 4 2 1]);
lat1 = latA([13 12 11 10 4 2 1]);
secT11 = T1A(:,[13 12 11 10 4 2 1]);
secT12 = T2A(:,[13 12 11 10 4 2 1]);
secS11 = sal1A(:,[13 12 11 10 4 2 1]);
secS12 = sal2A(:,[13 12 11 10 4 2 1]);
secdens11 = dens1A(:,[13 12 11 10 4 2 1]);
secdens12 = dens2A(:,[13 12 11 10 4 2 1]);
 
cmin = [11.63;
    35.28;
    23.27];
cmax = [28.09;
    36.43;
    28];
cint = [1.00;
    0.05;
    0.2];
 
 
r = 6366; %rayon  de la terre
loni = -33.8167; %long
lati = -3.86667;  %lat
%Conversion des entrées en ° vers des Radians
lonir = loni*pi/180;
latir = lati*pi/180;
lon1r = lon1*pi/180;
lat1r = lat1*pi/180;
%Formule précise
%d=2*asin(sqrt( (sin((lat1-lat2)/2))^2 + cos(lat1)*cos(lat2)*(sin((lon1-lon2)/2))^2))
dp = NaN(7,1);
dpA1 = NaN(7,1);
for i = 1:7;
    dp(i) = 2*asin(sqrt(( sin((latir-lat1r(i))/2))^2) + cos(latir)*cos(lat1r(i))* ( sin((lonir-lon1r(i))/2))^2) ;
    dpA1(i) = dp(i)* r;
end
%
%
%
%
 
%figure
 
for i = 9:11;
 
    figure (i)
    colormap(jet((cmax(i-8)-cmin(i-8))/cint(i-8)));
 
    if i==9
        var2 = secT12;
        pas  = 11:1:28.09;
        var  = (' echantillon 1');
    elseif i == 10;
        var2 = secS12;
        pas  = 35:0.1:36.5;
        var  = (' echantillon 2');
    else
        var2 = secdens12;
        pas  = 23.27:0.2:28;
        var  = (' echantillon 3');
 
    end
 
    contourf(dpA1,depthA(:,1),var2,pas);
 
    %contour(dpA1,depthA(:,1),var2,pas);
 
    %  pcolor(dpA1,depthA(:,1),var2);
    %   shading flat
 
    caxis([cmin(i-8) cmax(i-8)])
    set(gca,'ylim',[1 300],'ydir','reverse','ytick',(1:20:300),'fontsize',12,'fontweight','bold')
 
    colorbar('vert','ytick',pas)
   title(strcat('Evolution vertical de ',var,'- Atoll'));
 
 
 
end

Cordialement

[Jerome Briot]

Pourrais-tu nous fournir quelques données (le fichier CFI_fisica_Atol.mat ?) pour pouvoir faire des tests ?

[KAKATTAH]

Bonsoir DUT


Merci tout d'abord de ta disponibilté!
Biensûr pour le fichier; il est en pièce jointe.
J'ai crée une nouvelle variable pour ne prendre en compte que celle qui me créait ce problème afin de ne pas vous encombrer avec le grand fichier.
Merci

[KAKATTAH]

Voici la Pièce jointe: elle n'a pas été prise en compte dans le precedent message. Cette variable est la 3eme instruction de la boucle "if" dans le code (ie dans l'instructin de "else"): Les autres variables ne me causent pas de problème.

Merci