ClipperPushBroomCamera.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "ClipperPushBroomCamera.h"
 
#include "CameraDistortionMap.h"
#include "CameraFocalPlaneMap.h"
#include "iTime.h"
#include "LineScanCameraGroundMap.h"
#include "LineScanCameraSkyMap.h"
#include "NaifStatus.h"
 
namespace Isis {
  ClipperPushBroomCamera::ClipperPushBroomCamera(Cube &cube) : LineScanCamera(cube) {
 
    m_spacecraftNameLong = "Europa Clipper";
    m_spacecraftNameShort = "Clipper";
 
    int frameCode = naifIkCode();
 
    if (frameCode == -159103) {
      m_instrumentNameLong  = "Europa Imaging System Push Broom Narrow Angle Camera";
      m_instrumentNameShort = "EIS-PBNAC";
    }
    else if (frameCode == -159104) {
      m_instrumentNameLong  = "Europa Imaging System Push Broom Wide Angle Camera";
      m_instrumentNameShort = "EIS-PBWAC";
    }
    else {
      QString msg = "Unable to construct Clipper Push Broom camera model. "
                    "Unrecognized NaifFrameCode [" + toString(frameCode) + "].";
      throw IException(IException::User, msg, _FILEINFO_);
    }
 
    NaifStatus::CheckErrors();
 
    Pvl &lab = *cube.label();
 
    PvlGroup &bandBin = lab.findGroup("BandBin", Pvl::Traverse);
    QString key = "INS" + toString(naifIkCode()) + "_" + bandBin["FilterName"][0] + "_FOCAL_LENGTH";
    SetFocalLength(Spice::getDouble(key));
 
    SetPixelPitch();
 
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
    QString startTime = inst["StartTime"];
    iTime etStart(startTime);
 
    ReadLineRates(lab.fileName());
 
    // set up detector map
    new VariableLineScanCameraDetectorMap(this, p_lineRates);
 
    // Turn off the aberration corrections for the instrument position object
    instrumentPosition()->SetAberrationCorrection("NONE");
 
    // Set up focal plane map
    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());
    // center of array (same for WAC and NAC based on XY origin in EIS_Sensor_summary.xlsx)
    focalMap->SetDetectorOrigin(2047.5, 0.5);
 
    // Can be +/- 1024 pixel offset
    // Detector offset from label is 0 - 2047 but ISIS is zero based
    // with the detector offset being +/- 1024 pixels (2048/2)
    // Double check this against the SIS
    double detectorOffset = -(toDouble(inst["DetectorOffset"]) - 1024.0);
    focalMap->SetDetectorOffset(0.0, detectorOffset);
 
    // Set up distortion map
    CameraDistortionMap *distMap = new CameraDistortionMap(this);
    distMap->SetDistortion(naifIkCode());
 
    // Set up the ground and sky map
    new LineScanCameraGroundMap(this);
    new LineScanCameraSkyMap(this);
 
    setTime(etStart.Et());
 
    LoadCache();
    NaifStatus::CheckErrors();
  }
 

  ClipperPushBroomCamera::~ClipperPushBroomCamera() {
  }
 

   int ClipperPushBroomCamera::CkFrameId() const {
     return (-159000);
   }
 

   int ClipperPushBroomCamera::CkReferenceId() const {
     return (1);
   }
 

   int ClipperPushBroomCamera::SpkReferenceId() const {
     return (1);
   }

   void ClipperPushBroomCamera::ReadLineRates(QString filename) {
     Table timesTable("LineScanTimes", filename);
 
     if(timesTable.Records() <= 0) {
       QString msg = "Table [LineScanTimes] in [";
       msg += filename + "] must not be empty";
       throw IException(IException::Unknown, msg, _FILEINFO_);
     }
 
     for(int i = 0; i < timesTable.Records(); i++) {
       p_lineRates.push_back(LineRateChange((int)timesTable[i][2],
                                            (double)timesTable[i][0],
                                            timesTable[i][1]));
     }
 
     if(p_lineRates.size() <= 0) {
       QString msg = "There is a problem with the data within the Table ";
       msg += "[LineScanTimes] in [" + filename + "]";
       throw IException(IException::Unknown, msg, _FILEINFO_);
     }
   }
}

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