NewHorizonsLorriCamera.cpp

Go to the documentation of this file.

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

using namespace std;

namespace Isis {
  NewHorizonsLorriCamera::NewHorizonsLorriCamera(Cube &cube) : FramingCamera(cube) {
    m_instrumentNameLong = "Long Range Reconnaissance Imager";
    m_instrumentNameShort = "LORRI";
    m_spacecraftNameLong = "New Horizons";
    m_spacecraftNameShort = "NewHorizons";
    
    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 SPICE kernels. Version 100 of the Lorri IK uses meters for focal length, 
    // units, version 200 uses mm, and has a units keyword.
    QString unitsKey = "INS" + toString(naifIkCode()) + "_FOCAL_LENGTH_UNITS";
    QString focalLengthUnits = Spice::getString(unitsKey);
    if (focalLengthUnits != "mm") {
      QString msg = QObject::tr("SPICE keyword [%1] is expected to be mm. [%2] was found instead").
                    arg(unitsKey).arg(focalLengthUnits);
       throw IException(IException::User, msg, _FILEINFO_);      
    }
    SetFocalLength(Spice::getDouble("INS" + toString(naifIkCode()) + "_FOCAL_LENGTH"));

    // 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.
    // Uncomment the 3 lines below to use the default distortion map (i.e., no distortion)
    //    new CameraDistortionMap(this, -1);
    //    CameraDistortionMap *distortionMap = new CameraDistortionMap(this, -1);
    //    distortionMap->SetDistortion(naifIkCode());
    // Changed the SPICE keyword names to work with LORRI IK version 200
    QString e2("INS" + toString(naifIkCode()) + "_OOC_EM");
    QString e5("INS" + toString(naifIkCode()) + "_OOC_EM");
    QString e6("INS" + toString(naifIkCode()) + "_OOC_EM");
    new NewHorizonsLorriDistortionMap(this, getDouble(e2, 0), getDouble(e5, 1), 
                                      getDouble(e6, 2), -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);
}