马上注册,结交更多好友,享用更多功能^_^
您需要 登录 才可以下载或查看,没有账号?立即注册
x
% 读取grd文件
filename = 'C:\Users\聆風祈\Desktop\论文\数据\precn-1961-2020.grd';
fileID = fopen(filename, 'r');
data = fread(fileID, [128*72, 720], 'float32');
fclose(fileID);
% 转换数据维度
data = reshape(data, [128, 72, 720]);
% 计算每年降水总量
year_num = 2020 - 1961 + 1;
annual_precipitation = zeros(128, 72, year_num);
for i = 1:year_num
start_month = (i-1) * 12 + 1;
end_month = i * 12;
annual_precipitation(:, :, i) = sum(data(:, :, start_month:end_month), 3);
end
% 输出到Excel
xlswrite('annual_precipitation.xlsx', permute(annual_precipitation, [2, 1, 3]));
大佬怎么把这个代码里面的年总和改成月总和,只需要2001到2020年即可
你好!根据你的要求,以下是修改后的MATLAB代码,用于计算2001年到2020年每个月的降水量:
% 读取grd文件
filename = 'C:\Users\聆風祈\Desktop\论文\数据\precn-1961-2020.grd';
fileID = fopen(filename, 'r');
data = fread(fileID, [128*72, 720], 'float32');
fclose(fileID);
% 转换数据维度
data = reshape(data, [128, 72, 720]);
% 选取2001-2020年的数据
start_year = 2001;
end_year = 2020;
year_num = end_year - start_year + 1;
month_num = year_num * 12;
data_2001_2020 = data(:, :, (start_year - 1961) * 12 + 1 : (end_year - 1961 + 1) * 12);
% 月降水量
monthly_precipitation = data_2001_2020;
% 输出到Excel
for i = 1:month_num
sheetname = ['month_', num2str(i)];
xlswrite('monthly_precipitation.xlsx', monthly_precipitation(:, :, i), sheetname);
end
这段代码首先从1961-2020年的数据中提取出2001-2020年的数据。然后计算每个月的降水量,并将结果保存到一个Excel文件中,每个月的数据存储在一个单独的工作表中。
|