NewHorizonsMvicFrameCamera.cpp

Go to the documentation of this file.

#include "NewHorizonsMvicFrameCamera.h"

#include <QDebug>

#include "Camera.h"
#include "CameraDetectorMap.h"
#include "CameraDistortionMap.h"
#include "CameraFocalPlaneMap.h"
#include "CameraGroundMap.h"
#include "CameraSkyMap.h"
#include "IString.h"
#include "iTime.h"
#include "NewHorizonsMvicFrameCameraDistortionMap.h"
#include "NaifStatus.h"

using namespace std;

namespace Isis {

  NewHorizonsMvicFrameCamera::NewHorizonsMvicFrameCamera(Cube &cube) : FramingCamera(cube) {
    NaifStatus::CheckErrors();

    SetFocalLength();
    SetPixelPitch();

    // Get the start time from labels
    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
    m_exposure = inst["ExposureDuration"];
    QString stime = inst["SpacecraftClockStartCount"];
    // **  TODO  **  Need an offset time added to labels at ingestion??  The 0.125 value is 
    //     the value in DELTAT00.
    double offset = 0.125;
    m_etStart = getClockTime(stime).Et() + offset;
    SpiceChar utc[30];
    et2utc_c(m_etStart, "ISOC", 3, 30, utc);
//  qDebug()<<"\n\nspacecraftClockStartCount + "<<offset<<" (offset) = "<<utc;

    // 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.  Also save times of each framelet which are stored
    // in the BandBin group.
    PvlGroup &bandBin = lab.findGroup("BandBin", Pvl::Traverse);
    PvlKeyword &origBand = bandBin["OriginalBand"];
    PvlKeyword &utcTime = bandBin["UtcTime"];
    for(int i = 0; i < origBand.size(); i++) {
      m_originalBand.push_back(toInt(origBand[i]));
      m_utcTime.push_back(utcTime[i]);
    }

    CameraDetectorMap *detectorMap = new CameraDetectorMap(this);
    detectorMap->SetDetectorSampleSumming(1);
    detectorMap->SetDetectorLineSumming(1);

    // Setup focal plane map. The class will read data from the instrument addendum kernel to pull
    // out the affine transforms from detector samp,line to focal plane x,y.
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
    focalMap->SetDetectorOrigin(2500.5, 64.5);

    // Read distortion coefficients and boresight offsets from the instrument kernels. Then
    // construct the distortion map.
    //read the distortion coefs from the NAIF Kernels
    QString naifXKey = "INS-98900_DISTORTION_COEF_X";
    QString naifYKey = "INS-98900_DISTORTION_COEF_Y";
    QString naifppKey = "INS-98900_PP_OFFSET";
    vector<double> distCoefX;
    vector<double> distCoefY;

    for (int i=0; i < 20; i++) {
      distCoefX.push_back(getDouble(naifXKey,i));
      distCoefY.push_back(getDouble(naifYKey,i));
    }

    new NewHorizonsMvicFrameCameraDistortionMap(this, distCoefX, distCoefY);

    // Setup the ground and sky map
    new CameraGroundMap(this);
    new CameraSkyMap(this);

    // Internalize all the NAIF SPICE information into memory.
    LoadCache();
    NaifStatus::CheckErrors();
  }


  void NewHorizonsMvicFrameCamera::SetBand(const int vband) {

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

    iTime time(m_utcTime[vband-1]);
    double et = time.Et();
   
    SpiceChar utc[30];
    et2utc_c(et, "ISOC", 3, 30, utc);
    Camera::setTime(et);
    pair<iTime, iTime> shuttertimes = ShutterOpenCloseTimes(et, m_exposure);

    //  Set up valid band access
    Camera::SetBand(vband);
    
  }


  pair<iTime, iTime> NewHorizonsMvicFrameCamera::ShutterOpenCloseTimes(double time, double exposureDuration) {

    return FramingCamera::ShutterOpenCloseTimes(time, exposureDuration);
  }
}


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