NewHorizonsLorriCamera.cpp

Go to the documentation of this file.

#include "NewHorizonsLorriCamera.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 "NewHorizonsLorriDistortionMap.h"

using namespace std;

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

    // The LORRI focal length is fixed and is designed not to change throught the operational 
    // temperature. The NAIF code, set in the ISIS labels, will be used to read a single focal 
    // length from the SP:ICE kernels, but the units need to be changed from m to mm.
    QString fl = "INS" + toString(naifIkCode()) + "_FOCAL_LENGTH";
    double focalLength = Spice::getDouble(fl);
    focalLength *= 1000.0;
    SetFocalLength(focalLength);

    // For setting the pixel pitch, the Naif keyword PIXEL_SIZE is used instead of the ISIS
    // default of PIXEL_PITCH, so set the value directly.
    QString pp = "INS" + toString(naifIkCode()) + "_PIXEL_SIZE";
    double pixelPitch = Spice::getDouble(pp);
    pixelPitch /= 1000.0;
    SetPixelPitch(pixelPitch);

    // Since the two summing modes are handeled via different Naif codes, force the summing modes
    // to 1 for both the 1x1 and 4x4 modes
    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());

    // The boresight position recorded in the IK is zero-based and therefore needs to be adjusted 
    // for ISIS
    double boresightSample = Spice::getDouble("INS" + toString(naifIkCode()) + "_CCD_CENTER",0) + 1.0;
    double boresightLine = Spice::getDouble("INS" + toString(naifIkCode()) + "_CCD_CENTER",1) + 1.0;
    focalMap->SetDetectorOrigin(boresightSample,boresightLine);

    // Setup distortion map. Start by reading the distortion coefficient from the instrument kernel.
    // Then construct the distortion model.
//    new CameraDistortionMap(this, -1);
//    CameraDistortionMap *distortionMap = new CameraDistortionMap(this, -1);
//    distortionMap->SetDistortion(naifIkCode());
    QString e2("INS" + toString(naifIkCode()) + "_EPSILON2");
    QString e5("INS" + toString(naifIkCode()) + "_EPSILON5");
    QString e6("INS" + toString(naifIkCode()) + "_EPSILON6");
    new NewHorizonsLorriDistortionMap(this, getDouble(e2, 0), getDouble(e5, 0), 
                                      getDouble(e6, 0), -1);

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

    // The observation start time and clock count for LORRI are based on the center of the exposure.
    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
    QString clockCount = inst["SpacecraftClockStartCount"];
    double et = getClockTime(clockCount).Et();
    double exposureDuration = (double)inst["ExposureDuration"] / 1000.0;

    pair<iTime, iTime> startStop = ShutterOpenCloseTimes(et, exposureDuration);
    setTime(et);

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

  pair<iTime, iTime> NewHorizonsLorriCamera::ShutterOpenCloseTimes(double time,
                                                         double exposureDuration) {
    return FramingCamera::ShutterOpenCloseTimes(time - exposureDuration / 2.0, exposureDuration);
  }



}

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