Commit aef9d31b authored by Julien Malik's avatar Julien Malik
Browse files

BUG: use the right multiplication order

parent 03f565d6
......@@ -52,6 +52,7 @@ MNFImageFilter< TInputImage, TOutputImage, TNoiseImageFilter, TDirectionOfTransf
m_CovarianceEstimator = CovarianceEstimatorFilterType::New();
m_NoiseCovarianceEstimator = CovarianceEstimatorFilterType::New();
m_Transformer = TransformFilterType::New();
m_Transformer->MatrixByVectorOn();
}
template <class TInputImage, class TOutputImage,
......
......@@ -48,6 +48,7 @@ PCAImageFilter< TInputImage, TOutputImage, TDirectionOfTransformation >
m_CovarianceEstimator = CovarianceEstimatorFilterType::New();
m_Transformer = TransformFilterType::New();
m_Transformer->MatrixByVectorOn();
m_Normalizer = NormalizeFilterType::New();
}
......
......@@ -130,6 +130,7 @@ void VCAImageFilter<TImage>::GenerateData()
otbMsgDevMacroVCA( "Apply dimensionality reduction" )
// Xd = Ud.'*M;
typename MatrixImageFilterType::Pointer mulUd = MatrixImageFilterType::New();
mulUd->MatrixByVectorOn();
mulUd->SetInput(this->GetInput());
mulUd->SetMatrix(UdT);
mulUd->UpdateOutputInformation();
......@@ -179,6 +180,7 @@ void VCAImageFilter<TImage>::GenerateData()
vnl_matrix<PrecisionType> UdT = Ud.transpose();
typename MatrixImageFilterType::Pointer mulUd = MatrixImageFilterType::New();
mulUd->MatrixByVectorOn();
mulUd->SetInput(normalize->GetOutput());
mulUd->SetMatrix(UdT);
mulUd->UpdateOutputInformation();
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment