Calculate mean from daily to monthly values

Question

minh lan 2019-5-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/462430-calculate-mean-from-daily-to-monthly-values

回答： ABHILASH SINGH 2020-3-13

I have I have a data containing daily pricipitation data from 2001 to 2010, each year I have 1 file

lon=ncread(ncfile,'lon');
lat=ncread(ncfile,'lat');
precip=ncread(ncfile,'precip');
precip 360 x 280 x 365
lat 280x1
lon 360x1

I want to calculate the monthly mean for each year and then calculate average for 10 years.

the slow way I can do like that:

JAN{K}=nansum(precip{K}(:,:,1:31),3); %K is number of year
FEB{K}=nansum(precip{K}(:,:,32:59),3);
JAN=cat(3,JAN{:});
JAN=mean(JAN,3);

but for 2004, 2008 we have 366 day and this way is not very good to calculate

How can I do that with loop?

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

convert_to_metric 2019-5-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/462430-calculate-mean-from-daily-to-monthly-values#answer_375306

Hi minh lan,

You can consider making use of the month funciton. Take a look at the following code, perhaps it will lead you to a solution.

data=[1:365]*10; % just some example data
length_of_precip_data=length(data); % you might need to use the size function in place of length depending on how your data is organized
year=2010; % for example
start=datetime(year,1,1); 
m=month(start+caldays(1:length_of_precip_data)-1); % now you have an array the same size as your data that indicates which month an element is in
% so if you want to find just values from february:
data(m==2) % february is the second month, ie m==2