Plotting contours.txt (fr)

From GHER

Jump to: navigation, search

But : regarder la correspondance entre une topographie en format gher et des contours (ascii) issus d'une autre methode.

Modifications des fichiers liés à la topographie

  • Lecture de Calvi.Y et Calvi.X (format gher)
>> y=gread('Calvi.Y');  
>> x=gread('Calvi.X');
x=x-0.0005;
y=y-0.0005;

On selectionne la dernière couche (prof)

>> Y=y(:,:,31);
>> X=x(:,:,31);

On ajoute une ligne et une colonne à l'extérieure du domaine : (About pcolor)

>> X(123,:)=8.811, X(:,48)=X(:,n); 
#n = numeros de la dernière colonne qui existe ici égale à 47
>> Y(:,48)=42.595, Y(123,:)=Y(m,:); 
#m = numeros de la dernière colonne qui existe ici égale à 122
  • Lecture de calvi.bathymetry (format gher)
>> bathy=gread('Calvi.bathymetry');

On ajoute une ligne et une colonne à l'extérieure du domaine.

>> bathy(123,:)= NaN;
>> bathy(:,48)=NaN;
  • On plot (' car on veut la transposer (format gher))
>> pcolor(X',Y',bathy'), shading flat, colorbar

Manipulations du fichier ascii de contour

  • On ouvre le fichier et creer plusieurs fichiers pour chaque Polyline.
  • On importe l'ensemble de ces fichiers dans le Workspace de Matlab
  • On convertit l'ensemble des fichiers .txt creer précédemment en format gher
>>gwrite('tc50_a.gher', 'tc50_a'),
gwrite('tc50_b.gher', 'tc50_b'),
gwrite('tc50_c.gher', 'tc50_c'),
gwrite('tc50_d.gher', 'tc50_d'),
gwrite('tc50_e.gher', 'tc50_e'),
gwrite('tc50_f.gher', 'tc50_f')
  • On plot ensuite l'ensemble des contours sur une même figure (d'où le hold on):
>>hold on
>>plot(tc50_a(:,1),tc50_a(:,2)),
plot(tc50_b(:,1),tc50_b(:,2)),
plot(tc50_c(:,1),tc50_c(:,2)),
plot(tc50_d(:,1),tc50_d(:,2)),
plot(tc50_e(:,1),tc50_e(:,2)),
plot(tc50_f(:,1),tc50_f(:,2))
Personal tools