diff --git a/Modules/Filtering/Polarimetry/include/otbMuellerToReciprocalCovarianceImageFilter.h b/Modules/Filtering/Polarimetry/include/otbMuellerToReciprocalCovarianceImageFilter.h index 60ac764e19540673a78c634c2e0220a116b72b81..0aa2d8665815ca14c7005fd99fc32eb9f3532b64 100644 --- a/Modules/Filtering/Polarimetry/include/otbMuellerToReciprocalCovarianceImageFilter.h +++ b/Modules/Filtering/Polarimetry/include/otbMuellerToReciprocalCovarianceImageFilter.h @@ -79,25 +79,13 @@ public: const double M12 = static_cast<double>(Mueller[1]); const double M13 = static_cast<double>(Mueller[2]); const double M14 = static_cast<double>(Mueller[3]); - //const double M21 = static_cast<double>(Mueller[4]); const double M22 = static_cast<double>(Mueller[5]); const double M23 = static_cast<double>(Mueller[6]); const double M24 = static_cast<double>(Mueller[7]); - //const double M31 = static_cast<double>(Mueller[8]); - //const double M32 = static_cast<double>(Mueller[9]); const double M33 = static_cast<double>(Mueller[10]); const double M34 = static_cast<double>(Mueller[11]); - //const double M41 = static_cast<double>(Mueller[12]); - //const double M42 = static_cast<double>(Mueller[13]); - //const double M43 = static_cast<double>(Mueller[14]); const double M44 = static_cast<double>(Mueller[15]); - /*const ComplexType hhhh(M11+M22+M12+M21, 0.0); - const ComplexType hvhv(M11+M12-M21-M22, 0.0); - const ComplexType vvvv(M11+M22-M12-M21, 0.0); - const ComplexType hhhv(M13+M23, M14+M24); - const ComplexType hhvv(-M33-M44, M43-M34); - const ComplexType hvvv(M32-M31, M41-M42);*/ const ComplexType A(0.5*(M11+M22+2*M12)); const ComplexType B(0.5*vcl_sqrt(2.0)*(M13+M23), 0.5*vcl_sqrt(2.0)*(M14+M24)); diff --git a/Modules/Filtering/Polarimetry/test/otbMuellerToReciprocalCovarianceFunctor.cxx b/Modules/Filtering/Polarimetry/test/otbMuellerToReciprocalCovarianceFunctor.cxx index 925cfdb9dfc83ec55647c9037a11d7480bf84395..db596e64ebcca9b864e8e457e095775a9fd8fdde 100644 --- a/Modules/Filtering/Polarimetry/test/otbMuellerToReciprocalCovarianceFunctor.cxx +++ b/Modules/Filtering/Polarimetry/test/otbMuellerToReciprocalCovarianceFunctor.cxx @@ -68,7 +68,6 @@ int otbMuellerToReciprocalCovarianceFunctor(int itkNotUsed(argc), char * itkNotU vcl_abs(result[5].imag()-outputFunct[5].imag()) > 1e-10 ) { - std::cout.precision(15); std::cout<<"Test gives :"<<std::endl; std::cout<<outputFunct<<std::endl; std::cout<<"Wanted results are :"<<std::endl;