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