matlab编程（2）

onenuister

【今天练习的是用nc月平均地面2m温度画填色图】
程序如下：
%绘制7月多年平均地面2米温度
clear;
clc;
ncdisp(strcat('F:\001\air.2m.mon.mean.nc'));
lon=ncread('F:\001\air.2m.mon.mean.nc','lon');
lat=ncread('F:\001\air.2m.mon.mean.nc','lat');
T=ncread('F:\001\air.2m.mon.mean.nc','air');
x=lon(1:192);y=lat(1:47);
Tave=squeeze(mean(T(1:192,1:47,1,7:12:475),4));
Tave1=Tave-273.15*ones(192,47);%转换为摄氏度
[C,b]=contourf(x,y,Tave1',100);%填色
b.LineStyle='none';%去除等值线
hold on
colormap(flipud(hot));%翻转色标
colorbar
caxis([0,35]);%设置色标表示范围
load coast;
plot(long,lat,'-k','LineWidth',1.8);%画世界陆地界线
hold on;
guojie=shaperead('D:\Matlab\地图\chinamap\bou2_4l.shp');
bou2_4lx=[guojie(:).X];%提取经度信息
bou2_4ly=[guojie(:).Y];%提取纬度信息
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2)%绘国界
axis([70 140 15 55]);%设置坐标轴范围
set(gca,'XTick',70:10:140)%设置坐标轴样式
set(gca,'XTickLabel',{'70°E','80°E','90°E','100°E','110°E','120°E','130°E','140°E'})
set(gca,'YTick',15:10:55)
set(gca,'YTickLabel',{'15°N','25°N','35°N','45°N','55°N'});


（出图后可右键单击色标调整颜色范围）