GradientOperator.cpp

#include "GradientOperator.h"
#include "Chip.h"

namespace Isis {
  double GradientOperator::Interest(Chip &chip) {
    double gradient = 0.0;
    double pix1, pix2;
    int height = chip.Lines();
    int width = chip.Samples();
    int offset = 0;

    while(width > 1 && height > 1) {
      for(int i = 1; i <= width; i++) {
        pix1 = chip.GetValue(i + offset, 1 + offset);
        pix2 = chip.GetValue(width - i + 1 + offset, height + offset);
        if(ValidDnValue(pix1) && ValidDnValue(pix2)) gradient += std::abs(pix1 - pix2);
      }

      for(int i = 2; i < height; i++) {
        pix1 = chip.GetValue(1 + offset, i + offset);
        pix2 = chip.GetValue(width + offset, height - i + 1 + offset);
        if(ValidDnValue(pix1) && ValidDnValue(pix2)) gradient += std::abs(pix1 - pix2);
      }

      width -= 2;
      height -= 2;
      offset += 1;
    }

    return gradient;
  }
}

extern "C" Isis::InterestOperator *GradientOperatorPlugin(Isis::Pvl &pPvl) {
  return new Isis::GradientOperator(pPvl);
}