ERA5网格百分位数极端降水阈值提取

slc83328317

原始数据为ERA5的58年的日数据，需要提取90%分位数阈值。以我的数据为例（分全天、白天和夜间），41行，41列，21184天。提取阈值的时候，需要对非0降水值进行提取。


阈值提取
clc;clear   %本程序实现提取阈值90th
load('C:\Users\admin\Desktop\ERA5日总降水量\Allday_tp.mat')                    %导入数据
% load('C:\Users\admin\Desktop\ERA5日总降水量\Day_time_tp.mat')
% load('C:\Users\admin\Desktop\ERA5日总降水量\Night_time_tp.mat')
Threshold_90=zeros(41,41);   %先设置一个阈值的空矩阵，一会把结果写进去


for i=1:41
    for j=1:41
iGrid=Night_time_tp(:,i,j);     %取第i行，第j列的21184天的时间序列
iGrid(iGrid==0)=[];         %把0值去掉，不参与百分位阈值的计算
C=prctile(iGrid,90);        %提取非0的90th分位数
Threshold_90(i,j)=C;        %把i，j处的阈值写入阈值矩阵
disp(i)
disp(j)
    end
end

Threshold_allday=Threshold_90;

si18

可以上传一份示例数据吗 感谢

xuwenwen123

您好，这个数据书什么类型的啊

slc83328317

xuwenwen123 发表于 2022-4-2 19:35
您好，这个数据书什么类型的啊

类型？不明白唉，你说的是数据格式吗  NC的