Area weighting giving too low answers

2 次查看（过去 30 天）

Michael 2014-8-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/146497-area-weighting-giving-too-low-answers

回答： 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));