Matlab 用MODIS数据画全球海表温度分布图

MODIS数据可以从下面的网址下载
https://podaac.jpl.nasa.gov/datasetlist?ids=Collections&values=MODIS_L3_SST&view=list

clear,clc
file='文件名.nc';
sst=ncread(file,'sst');
lon=ncread(file,'lon');
lat=ncread(file,'lat');
[lat,lon]=meshgrid(lat,lon);
figure(1)
选择投影方式--汉莫尔-埃托夫投影，从西经150度开始
m_proj('hammer-aitoff','clongitude',-150);
先绘制东半球海表温度图
m_pcolor(lon,lat,sst);
去掉图上网格线
shading flat;colormap(parula);
hold on
绘制西半球海表温度图
m_pcolor(lon-360,lat,sst);
shading flat;colormap(parula);
添加海岸线
m_coast('patch',[.6 1 .6]);
m_grid('xaxis','middle');
h=colorbar('h');
set(get(h,'xlabel'),'string','MODISA NSST JAN 2019');
print(gcf,'-djpeg','-r1080','global_sst.jpg');