Mariner10Camera.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "Mariner10Camera.h"
 
#include <iostream>
#include <iomanip>
 
#include <SpiceUsr.h>
#include <SpiceZfc.h>
#include <SpiceZmc.h>
 
#include <QString>
 
#include "CameraDetectorMap.h"
#include "CameraFocalPlaneMap.h"
#include "CameraGroundMap.h"
#include "CameraSkyMap.h"
#include "IString.h"
#include "iTime.h"
#include "FileName.h"
#include "NaifStatus.h"
#include "Pvl.h"
#include "ReseauDistortionMap.h"
 
using namespace std;
namespace Isis {
  Mariner10Camera::Mariner10Camera(Cube &cube) : FramingCamera(cube) {
    NaifStatus::CheckErrors();
 
    m_spacecraftNameLong = "Mariner 10";
    m_spacecraftNameShort = "Mariner10";
 
    //  Turn off the aberration corrections for instrument position object
    instrumentPosition()->SetAberrationCorrection("NONE");
    instrumentRotation()->SetFrame(-76000);
 
    // Set camera parameters
    SetFocalLength();
    SetPixelPitch();
 
    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
    // Get utc start time
    QString stime = inst["StartTime"];
 
    iTime startTime;
    startTime.setUtc((QString)inst["StartTime"]);
    setTime(startTime);
 
    // Setup detector map
    new CameraDetectorMap(this);
 
    // Setup focal plane map, and detector origin
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
 
    QString ikernKey = "INS" + toString((int)naifIkCode()) + "_BORESIGHT_SAMPLE";
    double sampleBoresight = getDouble(ikernKey);
    ikernKey = "INS" + toString((int)naifIkCode()) + "_BORESIGHT_LINE";
    double lineBoresight = getDouble(ikernKey);
 
    focalMap->SetDetectorOrigin(sampleBoresight, lineBoresight);
 
    // Setup distortion map which is dependent on encounter, use start time
    // MOON:  1973-11-08T03:16:26.350
    QString spacecraft = (QString)inst["SpacecraftName"];
    QString instId = (QString)inst["InstrumentId"];
    QString cam;
    if(instId == "M10_VIDICON_A") {
      cam = "a";
      m_instrumentNameLong = "Mariner 10 Vidicon A";
      m_instrumentNameShort = "VIDICON A";
    }
    else if(instId == "M10_VIDICON_B") {
      cam = "b";
      m_instrumentNameLong = "Mariner 10 Vidicon B";
      m_instrumentNameShort = "VIDICON B";
    }
    else {
      QString msg = "File does not appear to be a Mariner10 image. InstrumentId ["
          + instId + "] is invalid Mariner 10 value.";
      throw IException(IException::User, msg, _FILEINFO_);
    }
 
    QString fname = FileName("$mariner10/reseaus/mar10" + cam
                             + "MasterReseaus.pvl").expanded();
 
    try {
      new ReseauDistortionMap(this, lab, fname);
    }
    catch(IException &e) {
      string msg = "Unable to create distortion map.";
      throw IException(e, IException::Programmer, msg, _FILEINFO_);
    }
 
    // Setup the ground and sky map
    new CameraGroundMap(this);
    new CameraSkyMap(this);
 
    LoadCache();
    NaifStatus::CheckErrors();
  }
 
 
  pair<iTime, iTime> Mariner10Camera::ShutterOpenCloseTimes(double time,
                                                            double exposureDuration) {
    pair<iTime, iTime> shuttertimes;
    // To get shutter start (open) time, subtract half exposure duration
    shuttertimes.first = time - (exposureDuration / 2.0);
    // To get shutter end (close) time, add half exposure duration
    shuttertimes.second = time + (exposureDuration / 2.0);
    return shuttertimes;
  }
}
 
 
extern "C" Isis::Camera *Mariner10CameraPlugin(Isis::Cube &cube) {
  return new Isis::Mariner10Camera(cube);
}