Area weighting giving too low answers

1 次查看（过去 30 天）
Michael
Michael 2014-8-6
回答： Michael 2014-8-7
I am trying the code below in order to area weight my data. The data looks good except it is 5 orders of magnitude lower than I want or expected. I have tried modifying a few things but can't seem to get it to work. Can anyone see where I have gone wrong?
[x,y]=meshgrid(0:1:360,-66:1:66);
x=x(1:end-1,1:end-1); y=y(1:end-1,1:end-1);
weight = cos(y(:,1));
wrepa = repmat(weight,[1 120 663]);
wrepa = permute(wrepa,[3 1 2]);
aa = areaquad(y(:,1)-(1/2),y(:,1)*0,y(:,1)+(1/2),(y(:,1)*0)+1);
repa = repmat(aa,[1 120 663]);
repa = permute(repa,[3 1 2]);
b = sum(aa)*360;
AMW = AM.*wrepa;
AMWW = AMW.*(repa/b);
AMWW = nanmean3d(AMWW,1);
AMWms = (smooth(AMWW,12));
AMWms = 663x132x120;
with values between 1 to -2 spatially.
thanks, Michael
  6 个评论
显示 4更早的评论
Michael
Michael 2014-8-6
AM is the pre-weighted data with same dimensions as AMWms

请先登录，再进行评论。

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

回答（1 个）

Michael
Michael 2014-8-7
Can anyone help me?

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by