Matlab软件求时间序列栅格降水数据的水文年型

最新推荐文章于 2024-07-10 19:42:42 发布

qq_41424763

最新推荐文章于 2024-07-10 19:42:42 发布

阅读量653

点赞数

文章标签： matlab 数据分析

本文链接：https://blog.csdn.net/qq_41424763/article/details/123961255

版权

该博客介绍了如何利用Matlab处理时间序列栅格降水数据，通过读取地理信息、数据重塑、计算经验频率法得到水文年型，并将结果保存为新的栅格文件。主要涉及数据预处理、循环操作和条件判断等步骤。

摘要由CSDN通过智能技术生成

[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