Discussion :

[KrokroAus]

Bonjour a tous,

Je cherche a créer une carte de valeur de réflectance mesurée par le satellite MODIS au Nord de l'Australie. Pour cela je veux utiliser un pcolor en lisant tout simplement 3 fichier texte contenant:
- des latitude
- des longitudes
- et des valeurs de réflectance correspondant aux couples de lat, long

La zone à cartographier est immense donc les fichiers txt sont lourds: 50 871 ko chacun.. et mon ordi est vieux et peu puissant.

Quand je lance MATLAB(7.0), j'ai l'erreur suivante:

??? Error using ==> figure
Attempt to write a locked object.

Est ce que quelqu'un a déjà eu ce genre de problème et est ce que c'est une erreur liée aux capacités de mon ordi?
Est-ce que les fichiers à lire sont particulièrement lourds pour un pcolor?

Merci d'avance pour toute aide!

[Invité]

Bonjour,

As-tu essayé de suivre Why does MATLAB crash when trying to produce a plot, run BENCH, or PCOLOR?

[KrokroAus]

Merci Winjerome pour ta réponse.
Oui, j'ai essayé de tapper feature('usegenericopengl',1), mais cela n'a rien changé à mon probleme: j'obtient une erruer de type : "Attempt to write a locked object" à chaque fois...

[Jerome Briot]

Essaie en diminuant la quantité de données à afficher en prenant un point sur cinq, voire sur dix ou sur cent.

[KrokroAus]

Merci Dut. Il y a t'il une façon simple de sélectionner les données tous les 1/10 ou 1/1000 ou est ce que je dois créer une boucle for -- end?

Voila mon code actuel pour créer les cartes:

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
    Lon=load([rep 'Longitude.txt']);
    Lat=load([rep 'Latitude.txt']);  
    CN=load([rep 'sur_refl_b01_1.txt']);  
    ReflB1=CN/10000/pi; % reflectance
 
    %Grid data
    nlon=1;
    for t=2:size(Lat)
        if(Lat(t)~=Lat(t-1)), break; end;
        nlon=nlon+1;
    end
    nlat=length(Lat)/nlon;
    Xp=zeros(nlat,nlon); 
    Yp=zeros(nlat,nlon);
    RB1=zeros(nlat,nlon);
 
    for j=1:nlat
        Xp(j,1:nlon)=Lon((j-1)*nlon+1:j*nlon);
        Yp(j,1:nlon)=Lat((j-1)*nlon+1:j*nlon);
        for i=1:nlon
            ReflB1t=ReflB1((j-1)*nlon+i);
            RB1(j,i) =ReflB1t; % reflectance à cartographier
        end       
    end
 
     fig=figure(1);clf %% Polynomiale - Turb
     set(fig, 'visible', 'off')
     pcolor(Xp,Yp,RB1); shading interp; 
     caxis([0 0.05]) ; axis equal tight ; colormap(jet2);colorbar('FontSize',20,'Fontweight','bold');
     saveas(gcf,[rep 'RB1_Insitu_Prof.png'],'png');

Merci d'avance,
Caroline

[Jerome Briot]

Quelles sont les dimensions des variables Xp, Yp et RB1 ?

Si ce sont des matrices MxN, tu peux faire ceci juste avant le code correspondant à l'affichage :

Code :

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

Xp = Xp(1:100:end,1:100:end);

Pareil pour les deux autres variables.