- 积分
- 6155
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
整理一下,方便使用。
1.求CMIP5数据多格点斜率
inPath='/data/NEX-GDDP/outdata/r85-min/';
inFiles=dir([inPath,'*.xlsx']);
fs=length(inFiles);
% x=1950:1:2005;
x=2006:1:2099;
for i=1:fs
file=([inPath,inFiles(i).name]);
data=xlsread(file,'1','C1:CR3963');
k=size(data,1);
for j =1:k
y=data(j,:);
[a m b]=regression(x,y);
m=m*10;
M(:,j)=m;
end
N=M';
Nn(:,i)=N;
end
2.求max,min,mean
a=xlsread('/data/NEX-GDDP/outdata/nex.xlsx','85-max');
A=a(:,2:22);
[max_A,index]=max(A,[],2);
[min_A,index]=min(a,[],2);
B=mean(A,2);
C=[max_A min_A B];
3.批量读取arcgis处理结果dbf
addpath('/home/syrus/桌面/model/code/');
% data1='/home/syrus/桌面/model/Js/80.dbf';
% b11=dbfread(data1);
inPath='/home/syrus/桌面/model/Mw/';
a=dir([inPath,'*.dbf']);
k=length(a);
for i=1:k
data=([inPath,'/',a(i).name]);
b=dbfread(data);
b1=b(:,7);
B(:,:,i)=b1(:,:);
end
|
|