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
| %% Methodo : importer le fichier WOA sous forme .
clear
close all
clc
%% preparation des fichiers texte
for num=36:58 % choix du fichier
clear CTD
%% Lis le fichier station WOrl OCean atlas
[WOD,echant,C,time,heure,a,b,c,d,e ] = textread(['0',num2str(num),'.txt'],'%s %s %s %s %s %f %f %f %f %f ','delimiter','\t', 'emptyvalue', NaN);
taille=length(e);
% determination des variables
CTD=[a(1:taille) b(1:taille) c(1:taille) d(1:taille) e(1:taille)];
time=time(1:taille);
heure=heure(1:taille);
% trouve les stations
tot_l_NaN=sum(isnan(CTD),2);
x=find(tot_l_NaN==1);
clear a b c d e echant tot_l_NaN WOD C
%% MEP des coordonée temporelle trouvé
len2=length(x);
Mois=ones(len2,1)+NaN;
Annee=ones(len2,1)+NaN;
Jour=ones(len2,1)+NaN;
Min=ones(len2,1)+NaN;
Heure=ones(len2,1)+NaN;
tt = 'Progress [ 0%% ]';
h = waitbar(0,tt);
for i=1:len2
waitbar(i/(len2),h,sprintf(':D [ %2.2f%% ] \n',100*i/ (len2)));
a=time{x(i)};
b=heure{x(i)};
[X{1:3}] = strread(a,'%f%f%f','delimiter','/');
[Y{1:2}] = strread(b,'%f%f','delimiter','/:');
mois=X{1};
jour=X{2};
annee=X{3};
min=Y{2};
heure2=Y{1};
if length(mois)>0
Mois(i)=mois;
Annee(i)=annee;
Jour(i)=jour;
Min(i)=min;
Heure(i)=heure2;
end
end
close (h)
coord_time=[Mois Jour Annee Heure Min];
%% on applique les coordonée temporelles au donnée de température et
%% profondeur
[n,m]=size(coord_time);
coo=ones(length(CTD),m)+NaN;
CTD=[coo CTD];
for i=1:len2
CTD(x(i),1:5)=coord_time(i,:);
end
clear Annee Heure JOUR Mois Min X Y a annee b heure heure2 i jour len min mois time Jour
clear coo coord_time h len2 m n tt x
close all
% on efface bottom depth
CTD(:,8)=[]; |
Partager