IssNACamera.cpp

Go to the documentation of this file.

// $Id: IssNACamera.cpp,v 1.6 2009/08/31 15:12:29 slambright Exp $
#include "IssNACamera.h"

#include "CameraDetectorMap.h"
#include "CameraFocalPlaneMap.h"
#include "CameraGroundMap.h"
#include "CameraSkyMap.h"
#include "IString.h"
#include "iTime.h"
#include "NaifStatus.h"
#include "RadialDistortionMap.h"
#include "iTime.h"

using namespace std;
namespace Isis {
  IssNACamera::IssNACamera(Cube &cube) : FramingCamera(cube) {
    Pvl &lab = *cube.label();

    // Get the filter wheels positions dependent focal length.
    // If we can not get the focal length for this specific filter wheel combination then
    // use the default.
    double focalLength = 0.0;
    try {
      PvlGroup bandBin = lab.findGroup("BandBin", Pvl::Traverse);
      QString key = QString("INS%1_%2_FOCAL_LENGTH").
                      arg(naifIkCode()).arg(bandBin["FilterName"][0]);
      key = key.replace("/", "_");
      focalLength = getDouble(key);   
    }
    catch (IException &firstException) {
      try {
        QString key = "INS-82360_DEFAULT_FOCAL_LENGTH";
        focalLength = getDouble(key);   
      }
      catch (IException &secondException) {
        PvlGroup bandBin = lab.findGroup("BandBin", Pvl::Traverse);
        IException finalError(IException::Unknown,
            QString("Unable to find a focal length for the requested Cassini ISS NA "
                    "filter combination [%1] or the default focal length")
                        .arg(bandBin["FilterName"][0]),
            _FILEINFO_);
        finalError.append(firstException);
        finalError.append(secondException);
        throw finalError;
      }
    }

    NaifStatus::CheckErrors();

    SetFocalLength(focalLength);
    SetPixelPitch();
    instrumentRotation()->SetFrame(Spice::getInteger("INS_" + toString(naifIkCode()) + "_FRAME_ID"));

    // Get the start time in et
    PvlGroup inst = lab.findGroup("Instrument", Pvl::Traverse);

    double et = iTime((QString)inst["StartTime"]).Et();

    // divide exposure duration keyword value by 1000 to convert to seconds
    double exposureDuration = ((double) inst["ExposureDuration"]) / 1000.0;
    pair<iTime, iTime> shuttertimes = ShutterOpenCloseTimes(et, exposureDuration);

    //correct time for center of exposure duration
    iTime centerTime = shuttertimes.first.Et() + exposureDuration / 2.0;

    // Setup detector map
    int summingMode = inst["SummingMode"];
    CameraDetectorMap *detectorMap = new CameraDetectorMap(this);
    detectorMap->SetDetectorLineSumming(summingMode);
    detectorMap->SetDetectorSampleSumming(summingMode);

    // 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
    double k1 = Spice::getDouble("INS" + toString(naifIkCode()) + "_K1");
    new RadialDistortionMap(this, k1);

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

    setTime(centerTime);
    LoadCache();
    NaifStatus::CheckErrors();
  }

  pair<iTime, iTime> IssNACamera::ShutterOpenCloseTimes(double time,
                                                        double exposureDuration) {
    return FramingCamera::ShutterOpenCloseTimes(time, exposureDuration);
  }
}

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