year=2001:2014;
lat=-89.5:1:89.5;
lon=-179.5:1:179.5;

WIND=ncread('/Users/yangyang/Aerosol_Variaiton/MODIS/ncfile/WIND_ECMWF.nc','WIND_DJF');
AOD=ncread('/Users/yangyang/Aerosol_Variaiton/MODIS/ncfile/AOD_MODIS.nc','AOD_DJF');

WIND_ECHINA=zeros(length(year),1);
AOD_ECHINA=zeros(length(year),1);

for yearID=1:length(year)
    n=0;
    for i=291:303
        for j=111:135
            if AOD(i,j,yearID)>0
                n=n+1;
                WIND_ECHINA(yearID)=WIND_ECHINA(yearID)+WIND(i,j,yearID);
                AOD_ECHINA(yearID)=AOD_ECHINA(yearID)+AOD(i,j,yearID);
            end
        end
    end
    WIND_ECHINA(yearID)=WIND_ECHINA(yearID)/n;
    AOD_ECHINA(yearID)=AOD_ECHINA(yearID)/n;
end

%WIND_ECHINA(:)=mean(mean(WIND(291:303,111:135,:),2),1);

Y=AOD_ECHINA;
X=[ones(length(year),1),WIND_ECHINA];
[b,~,~,~,stats]=regress(Y,X);
Slope=b(2);
R2=stats(1);