三维数组的统计分析

kongfeng0824

matlab中三维数组统计分析，不知大家有何方法？
例如：
一个nc文件中的降水precip是一个三维数组，其属性包含经度、维度、时间。即precip（lon,lat,time）。这里的时间也是个点数1-384.lon是1-144、lat是1-72。

1.    现在相对precip这个变量进行概率密度函数和分布函数的统计画图。不知道该如何？
2.    先要计算出所有经纬度格子里面的超过50mm/d的天数，并返回到每一个格子之中。当然还是返回到144*72*384的三维数组之中。我用了函数调用，但是始终返回不到相同的三维数组之中。我只能返回到一个二维数组之中。想要返回到三维数组之中，数据和原来不一样，也就不能reshape了。

我用t=fpeak( B2,A,0.0001,[min(B2),max(B2),D,inf]);计算出了在实践B2范围内的所有超过50毫米的降水。如何计算期天数，并返回到类似precip的数组之中呢？因为计算出来的t的数据个数和precip肯定不一样。所以reshape肯定回不到precip的144 *72 *384了。

不知大家在做gpcp数据的时候如何处理这类问题。

当然我的想法就是想统计出每一个经纬度格点在384个时间格子中超过50mm的天数，并进行频率统计。还请大家指点！！！！@言深深 @wlzhongouc @mofangbao @ztftom @兰溪之水 @斥鷃

mofangbao

不会matlab,帮顶了

kongfeng0824

mofangbao 发表于 2013-1-25 09:08

                               
登录/注册后可看大图

不会matlab,帮顶了

多谢清风！！！！！

斥鷃

LZ要是计算频数的话返回的应该是2D的呀，每个经纬度有多少天超过50mm/d的，为什么还是3D的呢？表示fpeak函数我这个版本里面也没有……我觉得我matlab版本是不是该换了……

kongfeng0824

斥鷃 发表于 2013-1-25 09:20

                               
登录/注册后可看大图

LZ要是计算频数的话返回的应该是2D的呀，每个经纬度有多少天超过50mm/d的，为什么还是3D的呢？表示fpeak函数 ...

我用的2012b。fpeak是自己写的。是算频数。但是precip本身是三维的。所以比较难办。不知道如何下手？

斥鷃

忽略我前面的回复，我突然领悟了，是不是直接在每个格点都做概率分布密度还是对整体来做？

wlzhongouc

我想你可以用find函数找出满足条件的格点  单用reshape返还不回去

huzengyun

现在我也在处理相同类型的数据，只是做简单的线性拟合。也希望多交流

huzengyun

现在我也在处理相同类型的数据，只是做简单的线性拟合。也希望多交流

kongfeng0824

huzengyun 发表于 2013-1-25 11:19

                               
登录/注册后可看大图

现在我也在处理相同类型的数据，只是做简单的线性拟合。也希望多交流

那可以分享一下你是如何做线性拟合的吗？是对一个格点一个格点处理吗？还是全局处理？我想可以深入交流一下。