MinimumDifference.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "MinimumDifference.h"
#include "Chip.h"
 
namespace Isis {
 
  double MinimumDifference::MatchAlgorithm(Chip &pattern, Chip &subsearch) {
    // calculate the sampling information
 
    double diff = 0.0;
    double count = 0;
    for(double l = 1.0; l <= pattern.Lines(); l++) {
      for(double s = 1.0; s <= pattern.Samples(); s++) {
        int line = (int)l;
        int samp = (int)s;
 
        double pdn = pattern.GetValue(samp, line);
        double sdn = subsearch.GetValue(samp, line);
        if(IsSpecial(pdn)) continue;
        if(IsSpecial(sdn)) continue;
        diff += fabs(pdn - sdn);
        count++;
      }
    }
 
    return diff / count;
  }
 
  bool MinimumDifference::CompareFits(double fit1, double fit2) {
    return (fit1 <= fit2);
  }
}
 
extern "C" Isis::AutoReg *MinimumDifferencePlugin(Isis::Pvl &pvl) {
  return new Isis::MinimumDifference(pvl);
}