VoyagerCamera.cpp

Go to the documentation of this file.

#include "VoyagerCamera.h"

#include <SpiceUsr.h>

#include <QString>

#include "CameraDetectorMap.h"
#include "CameraFocalPlaneMap.h"
#include "CameraGroundMap.h"
#include "CameraSkyMap.h"
#include "FileName.h"
#include "IString.h"
#include "iTime.h"
#include "NaifStatus.h"
#include "ReseauDistortionMap.h"
#include "Spice.h"

using namespace std;

namespace Isis {
  VoyagerCamera::VoyagerCamera (Cube &cube) : FramingCamera(cube) {
    NaifStatus::CheckErrors();

    // Set the pixel pitch
    SetPixelPitch();
    SetFocalLength();
    // Find out what camera is being used, and set the focal length, altinstcode,
    // and camera
    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup ("Instrument",Pvl::Traverse);
    QString spacecraft = (QString)inst["SpacecraftName"];
    QString instId = (QString)inst["InstrumentId"];

    QString reseauFileName = "";

    // These set up which kernel and other files to access,
    if (spacecraft == "VOYAGER_1") {
      p_ckFrameId = -31100;
      p_spkTargetId = -31;
      m_spacecraftNameLong = "Voyager 1";
      m_spacecraftNameShort = "Voyager1";

      reseauFileName += "1/reseaus/vg1";

      if (instId == "NARROW_ANGLE_CAMERA") {
        reseauFileName += "na";
        m_instrumentNameLong = "Narrow Angle Camera";
        m_instrumentNameShort = "NAC";
      }
      else if (instId == "WIDE_ANGLE_CAMERA") {
        reseauFileName += "wa";
        m_instrumentNameLong = "Wide Angle Camera";
        m_instrumentNameShort = "WAC";
      }
      else {
        QString msg = "File does not appear to be a Voyager image. InstrumentId ["
            + instId + "] is invalid Voyager value.";
        throw IException(IException::User, msg, _FILEINFO_);
      }
    }
    else if (spacecraft == "VOYAGER_2") {
      p_ckFrameId = -32100;
      p_spkTargetId = -32;
      m_spacecraftNameLong = "Voyager 2";
      m_spacecraftNameShort = "Voyager2";

      reseauFileName += "2/reseaus/vg2";

      if (instId == "NARROW_ANGLE_CAMERA") {
        reseauFileName += "na";
        m_instrumentNameLong = "Narrow Angle Camera";
        m_instrumentNameShort = "NAC";
      }
      else if (instId == "WIDE_ANGLE_CAMERA") {
        reseauFileName += "wa";
        m_instrumentNameLong = "Wide Angle Camera";
        m_instrumentNameShort = "WAC";
      }
      else {
        QString msg = "File does not appear to be a Voyager image. InstrumentId ["
            + instId + "] is invalid Voyager value.";
        throw IException(IException::User, msg, _FILEINFO_);
      }
    }
    else {
      QString msg = "File does not appear to be a Voyager image. SpacecraftName ["
          + spacecraft + "] is invalid Voyager value.";
      throw IException(IException::User, msg, _FILEINFO_);
    }

    new CameraDetectorMap(this);

    // Setup focal plane map, and detector origin
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
    focalMap->SetDetectorOrigin(500.0, 500.0);

    // Master reseau location file
    reseauFileName = "$voyager" + reseauFileName + "MasterReseaus.pvl";
    FileName masterReseaus(reseauFileName);
    try {
      new ReseauDistortionMap(this, lab, masterReseaus.expanded());
    } catch (IException &e) {
      e.print();
    }

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

    // StartTime is the most accurate time available because in voy2isis the
    // StartTime is modified to be highly accurate.
    // exposure duration keyword value is measured in seconds
    double exposureDuration = inst["ExposureDuration"];
    iTime startTime;
    startTime.setUtc((QString)inst["StartTime"]);

    // set the start (shutter open) and end (shutter close) times for the image
    /*****************************************************************************
     * AS NOTED IN ISIS2 PROGRAM lev1u_vgr_routines.c:
     * StartTime (FDS count) from the labels calculated to correspond the true spacecraft
     * clock count for the frame.  The true spacecraft clock count is readout
     * time of the frame, which occurred 2 seconds after shutter close.
     *****************************************************************************/
    pair<iTime, iTime> shuttertimes = ShutterOpenCloseTimes(startTime.Et(),
                                                            exposureDuration);

    // add half the exposure duration to the start time to get the center if the image
    iTime centerTime = shuttertimes.first.Et() + exposureDuration / 2.0;
    setTime(centerTime);

    LoadCache();
    NaifStatus::CheckErrors();
  }

  
  pair<iTime, iTime> VoyagerCamera::ShutterOpenCloseTimes(double time,
                                                          double exposureDuration) {
    pair<iTime, iTime> shuttertimes;
    // To get shutter end (close) time, subtract 2 seconds from the StartTime keyword value
    shuttertimes.second = time - 2;
    // To get shutter start (open) time, take off the exposure duration from the end time.
    shuttertimes.first = shuttertimes.second.Et() - exposureDuration;
    return shuttertimes;
  }
}

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