fichier netCDF et map avec m

fichier netCDF et map avec m_map

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:

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:

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:

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 :

Citation:

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

Montre nous le message d'erreur complet.

Que renvoi ceci ?
Code:

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:

SST = squeeze(SST);

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 :)