function I=HaarInverse(RR,DH,DV,DO)

I1=(RR+DH+DV+DO)/2;
I2=(RR+DH-DV-DO)/2; 
I3=(RR-DH+DV-DO)/2;
I4=(RR-DH-DV+DO)/2;
I=HaarInverseAux(I1,I2,I3,I4); 