HrscCamera.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "HrscCamera.h"
 
#include <string>
 
#include <QString>
 
#include "CameraDistortionMap.h"
#include "CameraFocalPlaneMap.h"
#include "iTime.h"
#include "LineScanCameraGroundMap.h"
#include "LineScanCameraSkyMap.h"
#include "NaifStatus.h"
#include "Statistics.h"
#include "VariableLineScanCameraDetectorMap.h"
 
using namespace std;
namespace Isis {
  HrscCamera::HrscCamera(Cube &cube) : LineScanCamera(cube) {
    m_instrumentNameLong = "High Resolution Stereo Camera";
    m_instrumentNameShort = "HRSC";
    m_spacecraftNameLong = "Mars Express";
    m_spacecraftNameShort = "MEX";
 
    NaifStatus::CheckErrors();
    // Setup camera characteristics from instrument and frame kernel
    SetFocalLength();
    SetPixelPitch(0.007);
    instrumentRotation()->SetFrame(-41210);
 
    // Get required keywords from instrument group
    Pvl &lab = *cube.label();
    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);
 
    ReadLineRates(lab.fileName());
 
    // Setup detector map for transform of image pixels to detector position
    new VariableLineScanCameraDetectorMap(this, p_lineRates);
    DetectorMap()->SetDetectorSampleSumming(inst["Summing"]);
 
    // Setup focal plane map for transform of detector position to
    // focal plane x/y.  This will read the appropriate CCD
    // transformation coefficients from the instrument kernel
 
    new CameraFocalPlaneMap(this, naifIkCode());
 
    QString ikernKey = "INS" + toString(naifIkCode())  + "_BORESIGHT_SAMPLE";
    double sampleBoresight = getDouble(ikernKey);
 
    ikernKey = "INS" + toString(naifIkCode())  + "_BORESIGHT_LINE";
    double lineBoresight = getDouble(ikernKey);
 
    FocalPlaneMap()->SetDetectorOrigin(sampleBoresight, lineBoresight);
 
    // Setup distortion map.  This will read the optical distortion
    // coefficients from the instrument kernel
    new CameraDistortionMap(this);
 
    // Setup the ground and sky map to transform undistorted focal
    // plane x/y to lat/lon or ra/dec respectively.
    new LineScanCameraGroundMap(this);
    new LineScanCameraSkyMap(this);
 
    LoadCache();
    NaifStatus::CheckErrors();
  }
 
 
  HrscCamera::~HrscCamera() {
  }
 
 
  int HrscCamera::CkFrameId() const {
    return (-41001);
  }
 
 
  int HrscCamera::CkReferenceId() const {
    return (1);
  }
 
 
  int HrscCamera::SpkReferenceId() const {
    return (1);
  }
 
 
  void HrscCamera::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 *HrscCameraPlugin(Isis::Cube &cube) {
  return new Isis::HrscCamera(cube);
}