- 积分
- 6
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-5-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近在进行各种干旱指数的计算,今天先在这里分享一下water storage deficit index (WSDI)的计算程序,希望大家多多交流指正。- clear;
- TWSK=load('test.txt');%2003/01-2015/12,a time series of water storage anomaly
- for i=1:13
- yK=TWSK(1+(i-1)*12:12+(i-1)*12);%y every year data
- rawK(i,:)=yK;%
- end
- for j=1:12
- damK=rawK(:,j);%13 year every day data
- meanK(j)=mean(damK);
- end
- for i=1:13
- for j=1:12
- WSDK(i,j)=rawK(i,j)-meanK(j);
- end
- end
- for i=1:13
- WSDK1(1+(i-1)*12:12+(i-1)*12)=WSDK(i,:);
- end
- WSDK1=WSDK1';
- WSDIK=zscore(WSDK1);
- save('WSDI_text.txt','WSDIK','-ascii');
复制代码
|
评分
-
查看全部评分
|