1.导入设置好的经纬度数据和色标条(另一篇文章: MATLAB保存遥感影像经纬度文件和色标条)
load('jingweidu.mat','lat','lon'); %load经纬度数据
load('mycolormap.mat'); %load色标条
mycolormap(end,:) =[0.5,0.5,0.5]; %最后为灰色
mycolormap(1,:) =[0,0,0]; %最大为白色
2.设置显示的值范围
tmp=data; %数据矩阵
a=0; %值的显示空间最小值
b=4.8; %值的显示空间最大值
3.设置窗口、背景值显示问题等
h=figure;
set(h,'position',[100 100 900 700]); %设置窗口大小
tmp(mask==0) = b; %把背景显示为白色
imagesc(lon,lat,tmp,[a,b]); %显示,值只显示[a,b]范围之间的
4.设置引用色标条
colormap(mycolormap);
bar=colorbar;
5.设置横纵坐标、标题等
set(get(bar,'title'),'string','mg/m^3');
str=get(gca,'xticklabel'); %获得句柄值
strtxt=strcat(str,'°E'); %字符串上加上°N'
set(gca,'xticklabel',strtxt,'FontSize',22,'FontName','宋体'); %重新设定横坐标%
str=get(gca,'yticklabel'); %获得句柄值
strtxt=strcat(str,'°N'); %字符串上加上°N'
set(gca,'yticklabel',strtxt,'FontSize',22,'FontName','宋体'); %重新设定横坐标
axis xy; %笛卡尔坐标系,原点在左下
title( '(a) ','fontsize',22,'FontName','宋体','position',[120 36.5]);