AlphaCube.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "AlphaCube.h"
 
#include "Cube.h"
 
using namespace std;
 
namespace Isis {
  AlphaCube::AlphaCube(Cube &cube) {
    Isis::PvlObject &isiscube = cube.label()->findObject("IsisCube");
    if(isiscube.hasGroup("AlphaCube")) {
      Isis::PvlGroup &alpha = isiscube.findGroup("AlphaCube");
      p_alphaSamples        = alpha["AlphaSamples"];
      p_alphaLines          = alpha["AlphaLines"];
      p_alphaStartingSample = alpha["AlphaStartingSample"];
      p_alphaStartingLine   = alpha["AlphaStartingLine"];
      p_alphaEndingSample   = alpha["AlphaEndingSample"];
      p_alphaEndingLine     = alpha["AlphaEndingLine"];
      p_betaSamples         = alpha["BetaSamples"];
      p_betaLines           = alpha["BetaLines"];
    }
    else {
      p_alphaSamples        = cube.sampleCount();
      p_alphaLines          = cube.lineCount();
      p_alphaStartingSample = 0.5;
      p_alphaStartingLine   = 0.5;
      p_alphaEndingSample   = (double) p_alphaSamples + 0.5;
      p_alphaEndingLine     = (double) p_alphaLines + 0.5;
      p_betaSamples         = p_alphaSamples;
      p_betaLines           = p_alphaLines;
    }
 
    ComputeSlope();
  }

  AlphaCube::AlphaCube(int alphaSamples, int alphaLines,
                       int betaSamples, int betaLines,
                       double alphaSs, double alphaSl,
                       double alphaEs, double alphaEl) {
    p_alphaSamples        = alphaSamples;
    p_alphaLines          = alphaLines;
    p_alphaStartingSample = alphaSs;
    p_alphaStartingLine   = alphaSl;
    p_alphaEndingSample   = alphaEs;
    p_alphaEndingLine     = alphaEl;
 
    p_betaSamples = betaSamples;
    p_betaLines = betaLines;
 
    ComputeSlope();
  }

  AlphaCube::AlphaCube(int alphaSamples, int alphaLines,
                       int betaSamples, int betaLines) {
    p_alphaSamples        = alphaSamples;
    p_alphaLines          = alphaLines;
    p_alphaStartingSample = 0.5;
    p_alphaStartingLine   = 0.5;
    p_alphaEndingSample   = (double) alphaSamples + 0.5;
    p_alphaEndingLine     = (double) alphaLines + 0.5;
 
    p_betaSamples = betaSamples;
    p_betaLines = betaLines;
 
    ComputeSlope();
  }

  void AlphaCube::Rehash(AlphaCube &add) {
    double sl = AlphaLine(add.AlphaLine(0.5));
    double ss = AlphaSample(add.AlphaSample(0.5));
    double el = AlphaLine(add.AlphaLine(add.BetaLines() + 0.5));
    double es = AlphaSample(add.AlphaSample(add.BetaSamples() + 0.5));
 
    p_alphaStartingLine = sl;
    p_alphaStartingSample = ss;
    p_alphaEndingLine = el;
    p_alphaEndingSample = es;
    p_betaLines = add.BetaLines();
    p_betaSamples = add.BetaSamples();
 
    ComputeSlope();
  }

  void AlphaCube::UpdateGroup(Isis::Cube &cube) {
    // If we have a mapping group we do not want to update the alpha cube
    // group as it represents the dimensions and sub-area of the raw instrument
    // cube
    Isis::PvlObject &cubeObject = cube.label()->findObject("IsisCube");
    if(cubeObject.hasGroup("Mapping")) return;
 
    // Add the labels to the pvl
    if(cubeObject.hasGroup("AlphaCube")) {
      AlphaCube temp(cube);
      temp.Rehash(*this);
      *this = temp;
 
      Isis::PvlGroup &alpha = cubeObject.findGroup("AlphaCube");
      alpha["AlphaSamples"] = toString(p_alphaSamples);
      alpha["AlphaLines"] = toString(p_alphaLines);
      alpha["AlphaStartingSample"] = toString(p_alphaStartingSample);
      alpha["AlphaStartingLine"] = toString(p_alphaStartingLine);
      alpha["AlphaEndingSample"] = toString(p_alphaEndingSample);
      alpha["AlphaEndingLine"] = toString(p_alphaEndingLine);
      alpha["BetaSamples"] = toString(p_betaSamples);
      alpha["BetaLines"] = toString(p_betaLines);
    }
    else {
      Isis::PvlGroup alpha("AlphaCube");
      alpha += Isis::PvlKeyword("AlphaSamples", toString(p_alphaSamples));
      alpha += Isis::PvlKeyword("AlphaLines", toString(p_alphaLines));
      alpha += Isis::PvlKeyword("AlphaStartingSample", toString(p_alphaStartingSample));
      alpha += Isis::PvlKeyword("AlphaStartingLine", toString(p_alphaStartingLine));
      alpha += Isis::PvlKeyword("AlphaEndingSample", toString(p_alphaEndingSample));
      alpha += Isis::PvlKeyword("AlphaEndingLine", toString(p_alphaEndingLine));
      alpha += Isis::PvlKeyword("BetaSamples", toString(p_betaSamples));
      alpha += Isis::PvlKeyword("BetaLines", toString(p_betaLines));
      cubeObject.addGroup(alpha);
    }
  }

  void AlphaCube::ComputeSlope() {
    p_lineSlope = double(p_alphaEndingLine - p_alphaStartingLine) /
                  double((p_betaLines + 0.5)       - 0.5);
    p_sampSlope = double(p_alphaEndingSample - p_alphaStartingSample) /
                  double((p_betaSamples + 0.5)       - 0.5);
  }
}