Isis 3.0 Object Programmers' Reference

Home

---

AlphaCube.cpp

Go to the documentation of this file.

#include "AlphaCube.h"

using namespace std;
namespace Isis {
  AlphaCube::AlphaCube (Isis::Pvl &pvl) {
    Isis::PvlObject &isiscube = pvl.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 {
    Isis::PvlGroup &dims = isiscube.FindGroup("Dimensions",Isis::Pvl::Traverse);
      p_alphaSamples        = dims["Samples"];
      p_alphaLines          = dims["Lines"];
      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::Pvl &pvl) {
    // 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 &cube = pvl.FindObject("IsisCube");
    if (cube.HasGroup("Mapping")) return;
    
    // Add the labels to the pvl
    if (cube.HasGroup("AlphaCube")) {
      AlphaCube temp(pvl);
      temp.Rehash(*this);
      *this = temp;
  
      Isis::PvlGroup &alpha = cube.FindGroup("AlphaCube");
      alpha["AlphaSamples"] = p_alphaSamples;
      alpha["AlphaLines"] = p_alphaLines;
      alpha["AlphaStartingSample"] = p_alphaStartingSample;
      alpha["AlphaStartingLine"] = p_alphaStartingLine;
      alpha["AlphaEndingSample"] = p_alphaEndingSample;
      alpha["AlphaEndingLine"] = p_alphaEndingLine;
      alpha["BetaSamples"] = p_betaSamples;
      alpha["BetaLines"] = p_betaLines;
    }
    else {
      Isis::PvlGroup alpha("AlphaCube");
      alpha += Isis::PvlKeyword("AlphaSamples",p_alphaSamples);
      alpha += Isis::PvlKeyword("AlphaLines",p_alphaLines);
      alpha += Isis::PvlKeyword("AlphaStartingSample",p_alphaStartingSample);
      alpha += Isis::PvlKeyword("AlphaStartingLine",p_alphaStartingLine);
      alpha += Isis::PvlKeyword("AlphaEndingSample",p_alphaEndingSample);
      alpha += Isis::PvlKeyword("AlphaEndingLine",p_alphaEndingLine);
      alpha += Isis::PvlKeyword("BetaSamples",p_betaSamples);
      alpha += Isis::PvlKeyword("BetaLines",p_betaLines);
      cube.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);
  }
}