Discussion :

[IFREMER]

Voici la fonction que j'ai éditée:

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
function [lon,lat,var,var2,lon2,lat2,d,m] = extract (y,m,d)
isOk=1;
isok = 1;
root = '/home9/begmeil/cersat/ftp/products/gridded/sst-pathfinder/netcdf/v5-all';
rootIma = '/home/typhon/stagiaires/jaccart/SST';

%----TEST----
if~exist(root,'dir')
       isOk=0;
       error(sprintf('ERROR, invalid directory: %s',root))
elseif~exist(rootIma,'dir')
       isOk=0;
       error(sprintf('ERROR, invalid directory: %s',rootIma))
end
 
%--NOM-DU-FICHIER--
    while m<=10
       while d<=31
		
            %--nom-du-fichier--
            fileName = sprintf('%.4d%.2d%.2d.04d2pfv50-all.nc.bz2',y,m,d);
            fileName = fullfile(root,sprintf('%.4d',y),sprintf('%.2d',m),fileName);
             %--Existance-du-fichier--
                if~exist(fileName,'file')
                isok = 0
                 end
 
           %----Lecture----
                if isok
                disp(fileName);
                nc = cl_netcdf_writer;
                nc = set_file_name (nc,fileName);
                s = read (nc);
                var = get_var_value (nc,s,'sst',1,1,1);
                lat = get_var_value (nc,s,'lat',1,1,1);
                lon = get_var_value (nc,s,'lon',1,1,1);
                var2 = var(lat>=-15 & lat<=30, lon>=-80 & lon<=20);
                end
    
            %----Sauvegarde-des-variables----
            if isok
            Name = fullfile(rootIma,sprintf('%.4d',y),sprintf('%.4d.%.2d.%.2d.mat',y,m,d));
            save(Name, 'var2');
            d = d+1;
            end
            
        end
    	m = m+1
	d=1
    end
    
end

J'aimerai juste savoir s'il existe un petit code magique ^^ pour que MATLAB charge seulement la partie de donnée qui m'intéresse (cad : var2 = var(lat>=-15 & lat<=30, lon>=-80 & lon<=20) ) plutôt que de lire tout le fichier et ensuite stocker dans var2 que les données utiles.

Je sais... je suis pénible avec ma fonction

[Jerome Briot]

Oui, mais seulement si tu connais à l'avance les numéros des lignes à lire.

[IFREMER]

Recoucou;
Comme ma fonction mettait trop de temps pour réaliser se que je voulais, je suis passé par le profiler et je me suis rendu compte que ce n'était pas simplifier le chargement du fichier que je devais faire mais réduire mes boucles.... donc en réfléchissant un peu....( si si ça m'arrive parfois ^^) je suis passé par la fonction DATE et maintenant le traitement du fichier prend 30 secondes au lieu de 5 minutes...

Merci pr toute l'aide apportée pour cette fonction ^^