SsiCamera.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "SsiCamera.h"
 
#include <QString>
 
#include "CameraDetectorMap.h"
#include "CameraDistortionMap.h"
#include "CameraFocalPlaneMap.h"
#include "CameraGroundMap.h"
#include "CameraSkyMap.h"
#include "IString.h"
#include "iTime.h"
#include "NaifStatus.h"
#include "Pvl.h"
#include "PvlObject.h"
#include "RadialDistortionMap.h"
#include "Spice.h"
 
using namespace std;
 
namespace Isis {
  SsiCamera::SsiCamera(Cube &cube) : FramingCamera(cube) {
    m_instrumentNameLong = "Solid State Imaging System";
    m_instrumentNameShort = "SSI";
    m_spacecraftNameLong = "Galileo Orbiter";
    m_spacecraftNameShort = "Galileo";
 
    NaifStatus::CheckErrors();
    // Get the camera characteristics
    double k1;
 
    Pvl &lab = *cube.label();
    iTime removeCoverDate("1994/04/01 00:00:00");
    iTime imageDate(lab.findKeyword("StartTime", PvlObject::Traverse)[0]);
    /*
    * Change the Focal Length and K1 constant based on whether or not the protective cover is on
    * See "The Direction of the North Pole and the Control Network of Asteroid 951 Gaspra"  Icarus 107, 18-22 (1994)
    */
    if(imageDate < removeCoverDate) {
      int code = naifIkCode();
      QString key = "INS" + toString(code) + "_FOCAL_LENGTH_COVER";
      SetFocalLength(Spice::getDouble(key));
      k1 = Spice::getDouble("INS" + toString(naifIkCode()) + "_K1_COVER");
    }
    else {
      SetFocalLength();
      k1 = Spice::getDouble("INS" + toString(naifIkCode()) + "_K1");
    }
 
    SetPixelPitch();
 
    // Get the start time in et
    PvlGroup inst = lab.findGroup("Instrument", Pvl::Traverse);
 
    double et = iTime((QString)inst["StartTime"]).Et();
 
    //?????????? NEED THESE??????
    // exposure duration keyword value is measured in seconds
    double exposureDuration = ((double) inst["ExposureDuration"]);
    pair<iTime, iTime> shuttertimes = ShutterOpenCloseTimes(et, exposureDuration);
 
    // Get summation mode
    double sumMode = inst["Summing"];
 
    // Setup detector map
    CameraDetectorMap *detectorMap = new CameraDetectorMap(this);
    detectorMap->SetDetectorSampleSumming(sumMode);
    detectorMap->SetDetectorLineSumming(sumMode);
 
    // Setup focal plane map
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
 
    focalMap->SetDetectorOrigin(
      Spice::getDouble("INS" + toString(naifIkCode()) +
                       "_BORESIGHT_SAMPLE"),
      Spice::getDouble("INS" + toString(naifIkCode()) +
                       "_BORESIGHT_LINE"));
 
    // Setup distortion map
    new RadialDistortionMap(this, k1);
 
    // Setup the ground and sky map
    new CameraGroundMap(this);
    new CameraSkyMap(this);
 
    setTime(et);
    LoadCache();
    NaifStatus::CheckErrors();
  }
 
  pair<iTime, iTime> SsiCamera::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 *SsiCameraPlugin(Isis::Cube &cube) {
  return new Isis::SsiCamera(cube);
}