Better way to implement:
Use matrix algebra!
l Construct basis matrix
B = [
W1(:)'
W2(:)'
W3(:)'
W4(:)']
l Flatten the test area:
CRf = CR(:)
l Compute coordinate
change:
  NewCoord = B*CRf
% = (10,0,0,0)’