【Matlab为全球某指数 如NDVI等 绘制显著性】

%% 读取底图与绘制
dataread= flip(imread(chl));
    %figure;imagesc(dataread)
hang=-89.75:0.5:89.75;
lie=-179.75:0.5:179.75;
[lon,lat]=meshgrid(lie,hang);
% imagescn(lon,lat,dataread)
C = dataread.*1000;%将0-1的数据与bar的10段式对应
imagescn(lon,lat,C);

cl_map=cmocean('delta');
% cmocean('delta','pivot',0);
%天蓝-蓝色：cmocean('delta')(1:128)
%黄色-绿色flip(cmocean('delta')(129:256))
lan=colormap(cl_map(1:128,:));
lv=colormap(cl_map(129:256,:));

housan=[lan(20,:);lan(60,:);lan(110,:)];
qianwu=[lv(1,:);lv(40,:);lv(60,:);lv(80,:);lv(100,:);lv(120,:)];
col=[housan;qianwu];
colormap(col);
axis off
caxis([-3,6]);

cbh = colorbar('location','eastoutside','Limits',[-3 6],'TickLength',0.005);%,'Ticks',-3:bin_length:7-bin_length
cbh.TickLabels = {'-3','-2','-1','-0','1','2','3','4','5', '6  (10^-^3 mg/m^3)',' '};

%% sig point添加显著性点

 i = 2;  %year
    tic
    sig = ['G:\linear trend\',filelistchl(i).name];
    sigpoint= flip(imread(sig));
    StatisticallySignificant = sigpoint<0.5;
   %figure;imagesc(StatisticallySignificant)
    hold on
    [lon,lat]=meshgrid(lie,hang);
    stipple(lon,lat,StatisticallySignificant,'density',100,'color','k','marker','.','markersize',5);