ApolloPanoramicCamera.cpp

#include "ApolloPanoramicCamera.h"

#include "ApolloPanIO.h"
#include "ApolloPanoramicDetectorMap.h"

#include "CameraDistortionMap.h"
#include "CameraFocalPlaneMap.h"
#include "IException.h"
#include "IString.h"
#include "iTime.h"
#include "LineScanCameraDetectorMap.h"
#include "LineScanCameraGroundMap.h"
#include "LineScanCameraSkyMap.h"
#include "PvlGroup.h"
#include "PvlKeyword.h"

using namespace std;
namespace Isis {
    ApolloPanoramicCamera::ApolloPanoramicCamera(Isis::Cube &cube) : Isis::LineScanCamera(cube) {
      // Set up the camera info from ik/iak kernels
      SetFocalLength(610.0);  //nominal (uncalibrated) focal length in mm from "Apollo 15 SIM Bay
                              // Photographic Equipment and Mission Summary" August, 1971
      SetPixelPitch(0.005);   //internally all images are modeled as if they have 5 micron pixels

      double  constantTimeOffset = 0.0,
              additionalPreroll = 0.0,
              additiveLineTimeError = 0.0,
              multiplicativeLineTimeError = 0.0;


      //following keywords in InstrumentAddendum file
      QString ikernKey = "INS" + toString((int)naifIkCode()) + "_CONSTANT_TIME_OFFSET";
      constantTimeOffset = getDouble(ikernKey);

      ikernKey = "INS" + toString((int)naifIkCode()) + "_ADDITIONAL_PREROLL";
      additionalPreroll = getDouble(ikernKey);

      ikernKey = "INS" + toString((int)naifIkCode()) + "_ADDITIVE_LINE_ERROR";
      additiveLineTimeError = getDouble(ikernKey);

      ikernKey = "INS" + toString((int)naifIkCode()) + "_MULTIPLI_LINE_ERROR";
      multiplicativeLineTimeError = getDouble(ikernKey);

      Pvl &lab = *cube.label(); 
      PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
      QString stime = (QString)inst["StartTime"];  
      SpiceDouble etStart;
      str2et_c(stime.toAscii().data(), &etStart);
      stime = (QString) inst["StopTime"];
      SpiceDouble etStop;
      str2et_c(stime.toAscii().data(), &etStop);
      iTime isisTime( (QString) inst["StartTime"]);
      
      // Get other info from labels
      // line exposure duration, sec/mm
      double lineRate = ( (double) inst["LineExposureDuration"] )*0.005;    

      lineRate *= 1.0 + multiplicativeLineTimeError;
      lineRate += additiveLineTimeError;
      etStart += additionalPreroll * lineRate;
      etStart += constantTimeOffset;

      setTime(isisTime);

      // Setup detector map
      //note (etStart+etStop)/2.0 is the time in the middle of image 
      //  (line = 0 after interior orientation)
      ApolloPanoramicDetectorMap *detectorMap = 
          new ApolloPanoramicDetectorMap((Camera *)this, 
                                         (etStart+etStop)/2.0, 
                                         (double)lineRate, &lab);
      //interior orientation residual stats
      m_residualMean = detectorMap->meanResidual();
      m_residualMax = detectorMap->maxResidual();
      m_residualStdev = detectorMap->stdevResidual();

      detectorMap->SetDetectorSampleSumming(1.0);
      detectorMap->SetStartingDetectorSample(0.0);
      // Setup focal plane map
      PvlGroup &kernel = lab.findGroup("Kernels", Pvl::Traverse);
      CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, (int) kernel["NaifFrameCode"]);

      //  Retrieve boresight location from instrument kernel (IK) (addendum?)
      double sampleBoreSight = 0.0;  //Presently no NAIF keywords for this sensor
      double lineBoreSight = 0.0;  //Presently no NAIF keywords for this sensor

      focalMap->SetDetectorOrigin(sampleBoreSight, lineBoreSight);
      focalMap->SetDetectorOffset(0.0, 0.0);

      // Setup distortion map
      new CameraDistortionMap(this, -1.0);
      //distMap->SetDistortion(naifIkCode());    Presently no NAIF keywords for this sensor

      //Setup the ground and sky map
      new LineScanCameraGroundMap(this);
      new LineScanCameraSkyMap(this);

      PvlGroup &instP = lab.findGroup("Kernels", Pvl::Traverse);
      m_CkFrameId = toInt(instP["NaifFrameCode"][0]);
      m_CkFrameId = -int(-m_CkFrameId/1000)*1000;

      LoadCache();
    }
}// end Isis namespace

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