NewHorizonsLeisaCamera.cpp

Go to the documentation of this file.

#include <QDebug>
#include <QVector>

#include "NewHorizonsLeisaCamera.h"
#include "NaifStatus.h"
#include "iTime.h"
#include "LineScanCameraGroundMap.h"
#include "LineScanCameraSkyMap.h"
#include "CameraFocalPlaneMap.h"
#include "LineScanCameraDetectorMap.h"

using namespace std;

namespace Isis {

  NewHorizonsLeisaCamera::NewHorizonsLeisaCamera(Cube &cube) : LineScanCamera(cube) {

    // Override the SPICE error process for SPICE calls 
    NaifStatus::CheckErrors(); 

    SetFocalLength();
    SetPixelPitch();

    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
    QString expDuration = inst["ExposureDuration"];

    QString stime = inst["SpacecraftClockStartCount"];
    double m_etStart = getClockTime(stime).Et(); 
    
    // The line rate is set to the time between each frame since we are treating LEASA as a linescan    
    double m_lineRate = expDuration.toDouble();

    // The detector map tells us how to convert from image coordinates to
    // detector coordinates. In our case, a (sample,line) to a (sample,time)
    new LineScanCameraDetectorMap(this, m_etStart, m_lineRate);

    // The focal plane map tells us how to go from detector position
    // to focal plane x/y (distorted).  That is, (sample,time) to (x,y) and back.
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
    focalMap->SetDetectorOrigin(128.5, 128.5);

    // Pull out the focal plane map affine coefficients so we can use them to adjust when the band
    // is changed. The coefficients as read from the iak are only valid for band 2. The constant
    // terms need to be multiplied by band-1 and then put back into the focal plane map
    m_origTransl.append(focalMap->TransL()[0]);
    m_origTransl.append(focalMap->TransL()[1]);
    m_origTransl.append(focalMap->TransL()[2]);

    m_origTranss.append(focalMap->TransS()[0]);
    m_origTranss.append(focalMap->TransS()[1]);
    m_origTranss.append(focalMap->TransS()[2]);

    m_origTransx.append(focalMap->TransX()[0]);
    m_origTransx.append(focalMap->TransX()[1]);
    m_origTransx.append(focalMap->TransX()[2]);

    m_origTransy.append(focalMap->TransY()[0]);
    m_origTransy.append(focalMap->TransY()[1]);
    m_origTransy.append(focalMap->TransY()[2]);

    // If bands have been extracted from the original image then we need to read the band bin group
    // so we can map from the cube band number to the instrument band number
    PvlGroup &bandBin = lab.findGroup("BandBin", Pvl::Traverse);
    PvlKeyword &orgBand = bandBin["OriginalBand"];
    for (int i = 0; i < orgBand.size(); i++) {
      m_originalBand.append(toInt(orgBand[i]));
    }

    // Use the defualt no correction distortion map. 
    new CameraDistortionMap(this);

    // Setup the ground and sky map
    new LineScanCameraGroundMap(this);
    new LineScanCameraSkyMap(this);

    LoadCache();  

    // Check to see if there were any SPICE errors
    NaifStatus::CheckErrors();  
  }


  void NewHorizonsLeisaCamera::SetBand(const int vband) {

    // Lookup the original band from the band bin group.
    if (vband > (int) m_originalBand.size()) {
     QString msg = QObject::tr("Band number out of array bounds in NewHorizonsLeisaCamera::SetBand legal "
                                "bands are [1-%1], input was [%2]").
                    arg(m_originalBand.size()).arg(vband);
      throw IException(IException::Programmer, msg, _FILEINFO_);
    }
    int band;
    band = m_originalBand[vband-1];

    // Get the affine coefficients from the focal plane map and adjust the constant terms to
    // provide the correct Y/Line offset for this band
    QVector<double> temp;

    temp.append(m_origTransl[0] * (band - 1));
    temp.append(m_origTransl[1]);
    temp.append(m_origTransl[2]);
    this->FocalPlaneMap()->SetTransL(temp);
    temp.clear();


    temp.append(m_origTranss[0] * (band - 1));
    temp.append(m_origTranss[1]);
    temp.append(m_origTranss[2]);
    this->FocalPlaneMap()->SetTransS(temp);
    temp.clear();

    temp.append(m_origTransx[0] * (band - 1));
    temp.append(m_origTransx[1]);
    temp.append(m_origTransx[2]);
    this->FocalPlaneMap()->SetTransX(temp);
    temp.clear();

    temp.append(m_origTransy[0] * (band - 1));
    temp.append(m_origTransy[1]);
    temp.append(m_origTransy[2]);
    this->FocalPlaneMap()->SetTransY(temp);

  }
}


extern "C" Isis::Camera *NewHorizonsLeisaCameraPlugin(Isis::Cube &cube) {
  return new Isis::NewHorizonsLeisaCamera(cube);
}