diff --git a/Modules/Filtering/Polarimetry/include/otbSinclairToReciprocalCircularCovarianceMatrixFunctor.h b/Modules/Filtering/Polarimetry/include/otbSinclairToReciprocalCircularCovarianceMatrixFunctor.h index 53c62b6339bd7760c3f6ec3cdedfc0ec28c548ec..6ae8b65f80b8e564dd493ff70fb61ce84aa9f06d 100644 --- a/Modules/Filtering/Polarimetry/include/otbSinclairToReciprocalCircularCovarianceMatrixFunctor.h +++ b/Modules/Filtering/Polarimetry/include/otbSinclairToReciprocalCircularCovarianceMatrixFunctor.h @@ -80,9 +80,15 @@ public: const ComplexType coef(0.5); const ComplexType j2S_hv = S_hv * ComplexType(0.0, 2.0); + const ComplexType jS_hh = S_hh * ComplexType(0.0, 1.0); + const ComplexType jS_vv = S_vv * ComplexType(0.0, 1.0); - const ComplexType Sll = coef * ( -S_hh-j2S_hv+S_vv ); + /*const ComplexType Sll = coef * ( -S_hh-j2S_hv+S_vv ); const ComplexType Slr = coef * ( -S_hh+-S_vv ); + const ComplexType Srr = coef * ( -S_hh+j2S_hv+S_vv );*/ + + const ComplexType Sll = coef * ( S_hh+j2S_hv-S_vv ); + const ComplexType Slr = coef * ( jS_hh + jS_vv ); const ComplexType Srr = coef * ( -S_hh+j2S_hv+S_vv ); //const ComplexType conjSll = vcl_conj(Sll);