Commit e147ef3b authored by Cédric Traizet's avatar Cédric Traizet
Browse files

BUG: fix offset formula: -mean/stddev

No related merge requests found
Showing with 7 additions and 10 deletions
+7 -10
......@@ -85,21 +85,18 @@ void LearningApplicationBase<TInputValue, TOutputValue>::TrainSharkKMeans(
auto stddevMeasurementVector = statisticsReader->GetStatisticVectorByName("stddev");
// Convert itk Variable Length Vector to shark Real Vector
shark::RealVector meanMeasurementRV(meanMeasurementVector.Size());
for (unsigned int i = 0; i<meanMeasurementVector.Size(); ++i)
{
// Substract the mean
meanMeasurementRV[i] = - meanMeasurementVector[i];
}
shark::RealVector offsetRV(meanMeasurementVector.Size());
shark::RealVector stddevMeasurementRV(stddevMeasurementVector.Size());
for (unsigned int i = 0; i<stddevMeasurementVector.Size(); ++i)
assert(meanMeasurementVector.Size()==stddevMeasurementVector.Size());
for (unsigned int i = 0; i<meanMeasurementVector.Size(); ++i)
{
stddevMeasurementRV[i] = stddevMeasurementVector[i];
// Substract the normalized mean
offsetRV[i] = - meanMeasurementVector[i]/stddevMeasurementVector[i];
}
shark::Normalizer<> normalizer(stddevMeasurementRV, meanMeasurementRV);
shark::Normalizer<> normalizer(stddevMeasurementRV, offsetRV);
centroidData = normalizer(centroidData);
}
......
Supports Markdown
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