Roberts basis in matlab
l Test region:
  CR = [
5 5
  5 5];
% Coefficients:
c1 = sum(sum(W1.*CR))
% c1 = 10
c2 = sum(sum(W2.*CR))
% c2 = 0
% Do same for W3 and
W4