Avg over a 3d array?

1 次查看（过去 30 天）
Kate
Kate 2013-7-12
Hi there!
So I have monthly gridded datasets in lat/lon/time(monthly): Example:
>> whos KoppenD_airT
Name Size Bytes Class Attributes
KoppenD_airT 360x720x361 374284800 single
>> whos KoppenDsa_airT
Name Size Bytes Class Attributes
KoppenDsa_airT 360x720x361 374284800 single
>> whos KoppenDsb_airT
Name Size Bytes Class Attributes
KoppenDsb_airT 360x720x361 374284800 single
I need to create climatologies of these regions so that I can plot monthly changes in temperature over my 30 year climate period. What is the best way to do this?
I need to aggregate my monthly data over the entire region to a single value so that I can plot all of my months.
Here's one attempt I've made that isn't working:
idx=find(Dsa);
for ii=idx
y=ind2sub(size(KoppenDsa_airT),ii);
x=squeeze(KoppenDsa_airT(y,:));
end
I get:
Index exceeds matrix dimensions.
Thanks!

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

采纳的回答

Matt J
Matt J 2013-7-12
编辑：Matt J 2013-7-12
aggregate my monthly data over the entire region to a single value
Does "aggregate" mean average? If so, this might be what you want,
mean(reshape(KoppenDsa_airT,[],361),1);
  4 个评论
显示 2更早的评论
Kate
Kate 2013-7-12
I realized my problem, I have a mask on my data so that non Koppen Dsa's are zero, I need to turn those into NaN's, thanks again for your help!
Matt J
Matt J 2013-7-12
If you turn them into NaNs, you would probably have to use nanmean()

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Climate Science and Analysis 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by