Discussion :

[Djeu Kuru]

Bonjour à tous,

Je cherche à faire une carte avec une donnée satellite récupéré via un fichier netCDF.

J'ai un peu du mal à comprendre comment gérer les fichier netCDF. Ici j'ai des données pour 12 mois.

Code :

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

1
2
3
4
5
6
7
 
file_name = 'SST_1993.nc';
% open
nc = netcdf.open(file_name,'NC_NOWRITE')
SST = ncread(file_name,'temperature');
lat = ncread(file_name, 'latitude');
lon = ncread(file_name,'longitude');

size(lon) = 1440 1
size(lat) = 689 1
size(SST) = 1440 689 1 12

Je ne comprends pas trop comment utiliser ce fichier SST étant donné ces dimensions ...
Je voudrais :
1) réaliser une moyenne sur les 12 mois de mon paramètre.

pour ce faire je fais :

Code :

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

sst_mean = nanmean(SST,4);

2) faire une carte (contourf?) de ces données en rajoutant la carte du monde sur mon raster.
En utilisant le package m_map :

Code :

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

1
2
3
[Plg,Plt]=meshgrid(lon,lat);
m_proj('hammer-aitoff','clongitude',-150);
m_pcolor(Plg,Plt,set_mean);shading flat;colormap(map);

Problème :

Error using pcolor (line 59)
Matrix dimensions must agree

Pourtant : size(sst_mean) = 1440 689 ... je ne comprends pas.

Est-ce que vous pouvez m'aiguiller ?

Merci beaucoup,
DJ

[Jerome Briot]

Montre nous le message d'erreur complet.

Que renvoi ceci ?

Code :

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

1
2
3
size(Plg)
size(Plt)
size(sst_mean)

Au passage, tu peux supprimer la dimension inutile de SST (la troisième égale à 1) en faisant :

Code :

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

SST = squeeze(SST);

[Djeu Kuru]

En faite ça marche !
Il suffisait uniquement d'inverser la longitude avec la latitude.

Merci pour la fonction squeeze !

Par contre je n'arrive pas à mettre par dessus une carte des continents. Est-ce que vous avez un conseil?
Merci beaucoup