SsiCamera.cpp

Go to the documentation of this file.

#include "SsiCamera.h"

#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) {
    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);
}