Area weighting giving too low answers

1 view (last 30 days)
Michael
Michael on 6 Aug 2014
Answered: Michael on 7 Aug 2014
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 Comments
Michael
Michael on 6 Aug 2014
AM is the pre-weighted data with same dimensions as AMWms

Sign in to comment.

Answers (1)

Michael
Michael on 7 Aug 2014
Can anyone help me?

Categories

Find more on MATLAB in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!