Ajout du gain, ajout de psiRB khiRB, parse PolBegin PolEnd...

Code/SARPolarimetry/otbMultiChannelsPolarimetricSynthesisFilter.h

@@ -134,7 +134,9 @@ public:
   itkGetMacro(EmissionV,bool);
   /** Set/Get Mode */    
   itkSetMacro(Mode,int);
-  itkGetMacro(Mode,int);  
+  itkGetMacro(Mode,int);
+  /** Set the gain */
+  itkSetMacro(Gain,double);  
   
   void ForceCoPolar();  
   
@@ -183,8 +185,6 @@ protected:
   void PrintSelf(std::ostream& os, itk::Indent indent) const;
   
   void Print();    
-  
-
 
 private :
   MultiChannelsPolarimetricSynthesisFilter(const Self&); //purposely not implemented
@@ -198,6 +198,9 @@ private :
   /** Khi Relechi */
   double m_KhiR;
   
+  /** Gain */
+  double m_Gain;
+  
   /** None = 0 , copolar = 1 , crosspolar = 2 */
   int m_Mode;
   

Code/SARPolarimetry/otbMultiChannelsPolarimetricSynthesisFilter.txx

@@ -39,6 +39,7 @@ MultiChannelsPolarimetricSynthesisFilter<TInputImage,TOutputImage,TFunction>
   SetEmissionH(false);
   SetEmissionV(false);  
   SetArchitectureType(0);
+  SetGain(1);  
 }
 
  /** PolarimetricSynthesisFilter
@@ -178,7 +179,7 @@ MultiChannelsPolarimetricSynthesisFilter<TInputImage,TOutputImage,TFunction>
           case 0 :
                 while( !inputIt.IsAtEnd() ) 
                 {
-                outputIt.Set( m_Functor( inputIt.Get()[0], inputIt.Get()[1], inputIt.Get()[2], inputIt.Get()[3] ) );
+                outputIt.Set( m_Gain * m_Functor( inputIt.Get()[0], inputIt.Get()[1], inputIt.Get()[2], inputIt.Get()[3] ) );
                 ++inputIt;
                 ++outputIt;
                 progress.CompletedPixel();  // potential exception thrown here
@@ -189,7 +190,7 @@ MultiChannelsPolarimetricSynthesisFilter<TInputImage,TOutputImage,TFunction>
           case 1 :
                 while( !inputIt.IsAtEnd() ) 
                 {
-                outputIt.Set( m_Functor( inputIt.Get()[0], inputIt.Get()[1], inputIt.Get()[1], inputIt.Get()[2] ) );
+                outputIt.Set( m_Gain * m_Functor( inputIt.Get()[0], inputIt.Get()[1], inputIt.Get()[1], inputIt.Get()[2] ) );
                 ++inputIt;
                 ++outputIt;
                 progress.CompletedPixel();  // potential exception thrown here
@@ -200,7 +201,7 @@ MultiChannelsPolarimetricSynthesisFilter<TInputImage,TOutputImage,TFunction>
           case 2 :
                 while( !inputIt.IsAtEnd() ) 
                 {
-                outputIt.Set( m_Functor( inputIt.Get()[0], inputIt.Get()[1], 0, 0 ) );
+                outputIt.Set( m_Gain * m_Functor( inputIt.Get()[0], inputIt.Get()[1], 0, 0 ) );
                 ++inputIt;
                 ++outputIt;
                 progress.CompletedPixel();  // potential exception thrown here
@@ -211,7 +212,7 @@ MultiChannelsPolarimetricSynthesisFilter<TInputImage,TOutputImage,TFunction>
           case 3 :
                 while( !inputIt.IsAtEnd() ) 
                 {
-                outputIt.Set( m_Functor( 0, 0, inputIt.Get()[2], inputIt.Get()[3] ) );
+                outputIt.Set( m_Gain * m_Functor( 0, 0, inputIt.Get()[2], inputIt.Get()[3] ) );
                 ++inputIt;
                 ++outputIt;
                 progress.CompletedPixel();  // potential exception thrown here

Code/SARPolarimetry/otbPolarimetricSynthesisFilter.h

@@ -97,7 +97,9 @@ public:
   /** Set/Get ArchitectureType */
   itkSetMacro(ArchitectureType,int);
   itkGetMacro(ArchitectureType,int);
-    
+  /** Set the gain */
+  itkSetMacro(Gain,double);
+  
   void SetEi(ComplexArrayType ei)
   {
        m_Ei = ei;
@@ -157,7 +159,10 @@ private:
   double m_PsiR;
   /** Khi Reflected */
   double m_KhiR;
-   
+
+  /** Gain */
+  double m_Gain;
+     
   /** None = 0 , copolar = 1 , crosspolar = 2 */
   int m_Mode;
   

Code/SARPolarimetry/otbPolarimetricSynthesisFilter.txx

@@ -34,6 +34,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
   this->SetNumberOfRequiredInputs(0);
   this->SetNumberOfInputs(4);
   SetMode(0);
+  SetGain(1);
   SetArchitectureType(0);
 }
 
@@ -327,7 +328,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
          
                 while( !inputItHH.IsAtEnd() ) 
                 {
-                    outputIt.Set( m_Functor(inputItHH.Get(), inputItHV.Get(), inputItVH.Get(), inputItVV.Get()) );
+                    outputIt.Set( m_Gain * m_Functor(inputItHH.Get(), inputItHV.Get(), inputItVH.Get(), inputItVV.Get()) );
                     ++inputItHH;
                     ++inputItHV;
                     ++inputItVH;
@@ -361,7 +362,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
 
                 while( !inputItHH.IsAtEnd() ) 
                 {
-                    outputIt.Set( m_Functor(inputItHH.Get(), inputItHV.Get(), inputItHV.Get(), inputItVV.Get()) );
+                    outputIt.Set( m_Gain * m_Functor(inputItHH.Get(), inputItHV.Get(), inputItHV.Get(), inputItVV.Get()) );
                     ++inputItHH;
                     ++inputItHV;
                     ++inputItVV;
@@ -394,7 +395,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
                 
                 while( !inputItHH.IsAtEnd() ) 
                 {
-                    outputIt.Set( m_Functor(inputItHH.Get(), inputItVH.Get(), inputItVH.Get(), inputItVV.Get()) );
+                    outputIt.Set( m_Gain * m_Functor(inputItHH.Get(), inputItVH.Get(), inputItVH.Get(), inputItVV.Get()) );
                     ++inputItHH;
                     ++inputItVH;
                     ++inputItVV;
@@ -423,7 +424,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
           
                 while( !inputItHH.IsAtEnd() ) 
                 {
-                    outputIt.Set( m_Functor(inputItHH.Get(), inputItHV.Get(), 0,0 ));
+                    outputIt.Set( m_Gain * m_Functor(inputItHH.Get(), inputItHV.Get(), 0,0 ));
                     ++inputItHH;
                     ++inputItHV;
                     ++outputIt;
@@ -450,7 +451,7 @@ PolarimetricSynthesisFilter<TInputImageHH,TInputImageHV,TInputImageVH,TInputImag
                 
                 while( !inputItVH.IsAtEnd() ) 
                 {
-                    outputIt.Set( m_Functor(0,0, inputItVH.Get(), inputItVV.Get()) );
+                    outputIt.Set( m_Gain * m_Functor(0,0, inputItVH.Get(), inputItVV.Get()) );
                     ++inputItVH;
                     ++inputItVV;
                     ++outputIt;