使用 Matlab 将每日数据（栅格）平均到每月平均值

Question

我有 2192 个栅格代表每日时间步长气候数据（2008 年 1 月 1 日至 2013 年 3 月 31 日），我想将其汇总为每月平均值（即给定月份和年份内的所有数据的平均值）。

我可以将数据加载到 Matlab 中，以创建大小为 110x57x2192 的 3D 矩阵。 2192 表示从 2008 年 1 月 31 日到 2013 年 3 月 31 日的时间维度。

我已经尝试过在 excel 中使用编号，但我不确定如何正确设置语法。理想情况下，月均值应替换内存中导入的每日数据。我在 Excel 中的编号有四列：日期、日期、月份和年份。忽略任何 NA 像素，像素的平均值基于给定月份的非 NA 像素。

%create monthly averages from daily data
numbering = xlsread('Numbering.xslx', B1:D2192);
month=numbering(:,2);
years=unique(year);

这是我用于编号的excel文件的前三行和后三行

2008 年 1 月 1 日 1 1 2008
2008 年 2 月 1 日 2 1 2008
2008 年 3 月 1 日 3 1 2008
29/12/2013 29 12 2013
30/12/2013 30 12 2013
31/12/2013 31 12 2013
31/12/2013 31 12 2013

Answer 1

以下是我解决问题的方法：

% Getting month breaks by detecting changes in days
day=numbering(:,2); % Day column from Excel
final_idx = find(diff(day)~=1);
initial_idx = [1 ; final_idx(1:end-1)-1];

%Number of months
nom = length(final_idx);

% Pre-allocating space for monthly mean
monthlymeans = nan(110,57,nom);

% Calculating means
for ii=1:nom
    monthlymeans(:,:,ii) = nanmean(3dmatrix(:,:,initial_idx(ii):final_idx(ii)),3);
end

请注意，我假设数据包含在一个名为 3dmatrix 的变量中，并且 NA 像素值实际上是 NaN，因为我无权访问您的原始数据。