[a,R]=geotiffread('C:\Users\WH\Desktop\EPRE\PRE\Idw_PRE_1985.tif');%先导入投影信息
info=geotiffinfo('C:\Users\WH\Desktop\EPRE\PRE\Idw_PRE_1985.tif');
[m,n]=size(a);
cd=2014-1985+1;%时间跨度，根据需要自行修改
datasum=zeros(m*n,cd)+NaN;
k=1;
for year=1985:2014 %起始年份
filename=['C:\Users\WH\Desktop\EPRE\PRE\Idw_PRE_',int2str(year),'.tif'];
data=importdata(filename);
data=reshape(data,m*n,1);
datasum(:,k)=data;
k=k+1;
end
result=zeros(m*n,cd)+NaN;
for i=1:size(datasum,1)
data=datasum(i,:);
if min(data)>0 %判断是否是有效值,我这里的有效值必须大于0
valuesum=[];
for k1=1:cd
value=length(find(data>=data(k1)))/(cd+1); %经验频率法求水文年型
if value <= 0.5 %这里我想把不同的水文年型赋值
value1 = 189*0
2025年Matlab软件求时间序列栅格降水数据的水文年型
Matlab软件求时间序列栅格降水数据的水文年型a R geotiffread C Users WH Desktop EPRE PRE Idw PRE 1985 tif 先导入投影信息 info geotiffinfo C Users WH Desktop EPRE PRE Idw PRE 1985 tif m n size a cd 2014 1985 1 时间跨度 根据需要自行修改
大家好，我是讯享网，很高兴认识大家。
数据结构与算法
上一篇
2025-02-08 08:06
2025年NIS详解
下一篇
2025-02-28 14:25
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容，请联系我们，一经查实，本站将立刻删除。
如需转载请保留出处：https://51itzy.com/kjqy/119036.html