GRiB 形式の JMA のデータセットを使用して、MATLAB で 2015 年 4 月の全球海面温度をプロットしようとしています。nctoolbox と m_map ツールボックスもインストールしました。
以下は私のコードです：
!wget http://ds.data.jma.go.jp/tcc/tcc/products/elnino/cobesst/gpvdata/2010-2019/sst201504.grb
nc=ncgeodataset('sst201504.grb')
nc.variables %to check the variable names in this file
lat=double('lat');
lon=double('lon');
sst=double(squeeze('Water_temperature_depth_below_sea'));
m_proj('miller','lat',[min(lat(:)) max(lat(:))],...'lon',[min(lon(:)) max(lon(:))])
m_pcolor(lon,lat,sst);
ただし、m-pcolor 関数を使用すると、次のエラー メッセージが生成されます。
Error using pcolor (line 53)
Color data input must be a matrix.
Error in m_pcolor (line 53)
[h]=pcolor(X,Y,data,varargin{:});
ただし、次のコードを使用して海岸線とグリッド線をプロットすることはできますが、色付きの温度異常はありません。
m_coast;
m_grid;
コードで何かを見逃していませんか? latとlonは1x3倍精度配列で、sstは1x33倍精度配列です。