利用regress（一元回归）计算hadisst 的sst数据与nino3画b值得等值线图

18368091971

%ncdump
%regresstion analysis for enso
close all
clear
%read sst from hadisst
nc=netcdf('global sst.nc','r');
timer=nc{'TIME'}(:);
%base=datenum(1,1,1,0,0,0);
%time=timer+base;
%datestr(time/24);
lon=nc{'LON'}(:);
lat=nc{'LAT'}(:);
miss=nc{'SST'}.missing_value(:);
fill=nc{'SST'}.FillValue_(:);
sst=nc{'SST'}(:,:,:);%ncfloat('TIME', 'LAT', 'LON');time=360
sst(sst==fill)=nan;
sst(sst==miss)=nan;
sst(sst<-10)=nan;  %min sst是 -1000，去除太小值
%%
%通常下载数据直接绘图是从-180 to 180 ;将地图经度换为 0 to 360
inn=find(lon<0);
inp=find(lon>=0);
sstn=sst(:,:,inn);
sstp=sst(:,:,inp);
sstr=cat(3,sstp,sstn);
lon1=lon(inn)+360;
lonr=[lon(inp)
    lon1];
   [long,latt]=meshgrid(lonr,lat);
%%
%扣除季节性变化
tmp=reshape(sstr,[12,30,180,360]);
season=squeeze(nanmean(tmp,2));
clim=repmat(season,[30,1,1]);
ssta=sstr-clim;
%load ninoindex
nino=load('nino index.txt');
enso=nino(1:360,5);%time 1982-2011
%plot(enso);
%regress
[tt,ii,jj]=size(ssta);
cof=nan*zeros(ii,jj);
cofs=nan*zeros(ii,jj);
reg=[ones(tt,1),enso];
for i=1:ii
    for j=1:jj
        tmp=squeeze(ssta(:,i,j));
        %if~isnan(tmp)
          [b,~,~,stats]=regress(tmp,reg);  
        %end
        cof(i,j)=b(2);
        cofs(i,j)=stats(3);
    end
end
cof(cof==0)=nan;
%%
%画图
m_proj('Robinson','long',[min(lonr),max(lonr)],...
                   'lat',[min(lat)+5,max(lat)-5]);
%[X,Y]=m_112xy(lon,lat,'clip','on');
m_grid('box','off','tickdir','on',...
        'fontsize',20);
hold on
m_contourf(long,latt,cof,'linestyle','none');
m_coast('patch',[.6 .6 .6]);
hold on
m_contour(long,latt,cofs,[0.05,0.05]);
colorbar;
myc1
caxis([-0.5,0.5]);
title('regresstion enso')
%问题，怎么把地图上的lable弄好点，或没有

Nancy_qd

谢谢楼主，学习了

常山赵子龙

楼主，我在使用您的程序的时候我自己出现了一个问题，就是在regress那一行，matlab提示我tmp和reg：
Y must be a vector and must have the same number of rows as X. 请问如何解决啊？

哆啦做饿梦

感谢楼主分享

2daidai

好像很厉害的样子

luyunlong

{:5_275:}{:5_275:}{:5_275:}

lhy

太感谢了，我照着中间一段程序做出来了，是不是regression的结果，在图形上看起来应该跟correlation的差不多？

另外请问楼主ENSO指数也标准化了吗？我看你在代码中对SST去除了平均值，一般在regression中标准化就是减去均值吗，还是zscore也可以？两组数据都要标准化吗？