Added LEISA camera

ifeq ($(ISISROOT), $(BLANK))

.SILENT:

error:

	echo "Please set ISISROOT";

else

	include $(ISISROOT)/make/isismake.objs

endif

 No newline at end of file

/**

 * @file

 *

 *   Unless noted otherwise, the portions of Isis written by the USGS are public

 *   domain. See individual third-party library and package descriptions for

 *   intellectual property information,user agreements, and related information.

 *

 *   Although Isis has been used by the USGS, no warranty, expressed or implied,

 *   is made by the USGS as to the accuracy and functioning of such software

 *   and related material nor shall the fact of distribution constitute any such

 *   warranty, and no responsibility is assumed by the USGS in connection

 *   therewith.

 *

 *   For additional information, launch

 *   $ISISROOT/doc//documents/Disclaimers/Disclaimers.html in a browser or see

 *   the Privacy & Disclaimers page on the Isis website,

 *   http://isis.astrogeology.usgs.gov, and the USGS privacy and disclaimers on

 *   http://www.usgs.gov/privacy.html.

 */

#include "NewHorizonsLeisaCamera.h"

#include "NaifStatus.h"

#include "iTime.h"

#include "LineScanCameraGroundMap.h"

#include "LineScanCameraSkyMap.h"

#include "CameraFocalPlaneMap.h"

#include "LineScanCameraDetectorMap.h"

using namespace std;

namespace Isis {

  /**

   * Constructs a New Horizons LEISA LineScanCamera object. LEISA is technically a frame

   * type camera, but it has the etalon filer in front of it exposing each line of a frame to have a

   * different wavelength, so we treat it like a linescan camera. Each band of the ISIS cube is made 

   * by combining all the corresponding frame line numbers into that band (i.e., All the line number 

   * 1s from each frame in an observatoin are combined into band 1, all line number 2s are put into 

   * band 2, and so on) 

   * 

   * @param something Description

   *

   * @author Kristin Berry

   * @internal

   */

  NewHorizonsLeisaCamera::NewHorizonsLeisaCamera(Cube &cube) : LineScanCamera(cube) {

    // Override the SPICE error process for SPICE calls 

    NaifStatus::CheckErrors(); 

    SetFocalLength();

    SetPixelPitch();

    Pvl &lab = *cube.label();

    PvlGroup &inst = lab.findGroup("Instrument", Pvl::Traverse);

    QString expDuration = inst["ExposureDuration"];

    QString stime = inst["SpacecraftClockStartCount"];

    double m_etStart = getClockTime(stime).Et(); 

    // The line rate is set to the time between each frame since we are treating LEASA as a linescan    

    double m_lineRate = expDuration.toDouble();

    // The detector map tells us how to convert from image coordinates to

    // detector coordinates. In our case, a (sample,line) to a (sample,time)

    new LineScanCameraDetectorMap(this, m_etStart, m_lineRate);

    // The focal plane map tells us how to go from detector position

    // to focal plane x/y (distorted).  That is, (sample,time) to (x,y) and back.

    CameraFocalPlaneMap *focalMap = new CameraFocalPlaneMap(this, naifIkCode());

    focalMap->SetDetectorOrigin(128.5, 128.5);

    // Use the defualt no correction distortion map. 

    new CameraDistortionMap(this);

    // Setup the ground and sky map

    new LineScanCameraGroundMap(this);

    new LineScanCameraSkyMap(this);

    LoadCache();  

    // Check to see if there were any SPICE errors

    NaifStatus::CheckErrors();  

  }

}

/**

 * This is the function that is called in order to instantiate a 

 * NewHorizonsLeisaCamera object. 

 *

 * @param lab Cube labels

 *

 * @return Isis::Camera* NewHorizonsLeisaCamera

 * 

 */

extern "C" Isis::Camera *NewHorizonsLeisaCameraPlugin(Isis::Cube &cube) {

  return new Isis::NewHorizonsLeisaCamera(cube);

}

#ifndef NewHorizonsLeisaCamera_h

#define NewHorizonsLeisaCamera_h

/**

 * @file

 *

 *   Unless noted otherwise, the portions of Isis written by the USGS are public

 *   domain. See individual third-party library and package descriptions for

 *   intellectual property information,user agreements, and related information.

 *

 *   Although Isis has been used by the USGS, no warranty, expressed or implied,

 *   is made by the USGS as to the accuracy and functioning of such software

 *   and related material nor shall the fact of distribution constitute any such

 *   warranty, and no responsibility is assumed by the USGS in connection

 *   therewith.

 *

 *   For additional information, launch

 *   $ISISROOT/doc//documents/Disclaimers/Disclaimers.html in a browser or see

 *   the Privacy & Disclaimers page on the Isis website,

 *   http://isis.astrogeology.usgs.gov, and the USGS privacy and disclaimers on

 *   http://www.usgs.gov/privacy.html.

 */

#include "LineScanCamera.h"

namespace Isis {

  /**

   * This is the camera model for LEISA, New Hoirzon's infrared 

   * Spectrometer.  

   *

   * @ingroup SpiceInstrumentsAndCameras

   * @ingroup New Horizons

   *

   * @author 2014-09-02 Kristin Berry

   *

   * @internal 

   *  

   */

  class NewHorizonsLeisaCamera : public LineScanCamera {

    public:

      //! Create a NewHorizonsLeisaCamera object

      NewHorizonsLeisaCamera(Cube &cube);

      //! Destroys the NewHorizonsLeisaCamera object

      ~NewHorizonsLeisaCamera() {};

       /**

       * CK frame ID -  - Instrument Code from spacit run on CK

       *  

       * @return @b int The appropriate instrument code for the "Camera-matrix" 

       *         Kernel Frame ID

       */

      virtual int CkFrameId() const { return (-98000); }//from NAIF_INSTRUMENT_ID in /usgs/cpkgs/isis3/data/newhorizons/kernels/ck/*.lbl or spacit R or spacit S on ck

      /** 

       *  CK Reference ID - 

       * 

       * @return @b int The appropriate code for the "Camera-matrix"  

       *         Kernel Reference ID

       */

      virtual int CkReferenceId() const { return (1); } //can get from spacit S on ck

      /** 

       *  SPK Reference ID - J2000

       *  

       * @return @b int The appropriate instrument code for the Spacecraft 

       *         Kernel Reference ID

       */

      virtual int SpkReferenceId() const { return (1); } //can get from spacit S on spk

  };

};

#endif

/**

 * @file

 *

 *   Unless noted otherwise, the portions of Isis written by the USGS are public

 *   domain. See individual third-party library and package descriptions for 

 *   intellectual property information,user agreements, and related information.

 *

 *   Although Isis has been used by the USGS, no warranty, expressed or implied,

 *   is made by the USGS as to the accuracy and functioning of such software 

 *   and related material nor shall the fact of distribution constitute any such 

 *   warranty, and no responsibility is assumed by the USGS in connection 

 *   therewith.

 *

 *   For additional information, launch

 *   $ISISROOT/doc//documents/Disclaimers/Disclaimers.html in a browser or see 

 *   the Privacy & Disclaimers page on the Isis website,

 *   http://isis.astrogeology.usgs.gov, and the USGS privacy and disclaimers on

 *   http://www.usgs.gov/privacy.html.

 */

#include <iomanip>

#include <iostream>

#include "Camera.h"

#include "CameraFactory.h"

#include "NewHorizonsLeisaCamera.h"

#include "IException.h"

#include "iTime.h"

#include "Preference.h"

#include "Pvl.h"

#include "PvlGroup.h"

using namespace std;

using namespace Isis;

void TestLineSamp(Camera *cam, double samp, double line);

int main(void) {

  Preference::Preferences(true);

  cout << "Unit Test for NewHorizonsLeisaCamera..." << endl;

  try {

    // These should be lat/lon at center of image. To obtain these numbers for a new cube/camera,

    // set both the known lat and known lon to zero and copy the unit test output "Latitude off by: "

    // and "Longitude off by: " values directly into these variables. -- in progress

    double knownLat = 0;

    double knownLon = 0; 

    Cube c("jupiterTest.cub", "r");//get a test cube!

    NewHorizonsLeisaCamera *cam = (NewHorizonsLeisaCamera *) CameraFactory::Create(c);

    cout << "FileName: " << FileName(c.fileName()).name() << endl;

    cout << "CK Frame: " << cam->instrumentRotation()->Frame() << endl << endl;

    cout.setf(std::ios::fixed);

    cout << setprecision(9);

    // Test kernel IDs

    cout << "Kernel IDs: " << endl;

    cout << "CK Frame ID = " << cam->CkFrameId() << endl;

    cout << "CK Reference ID = " << cam->CkReferenceId() << endl;

    cout << "SPK Target ID = " << cam->SpkTargetId() << endl;

    cout << "SPK Reference ID = " << cam->SpkReferenceId() << endl << endl;

    // Test all four corners to make sure the conversions are right

    //Where do these numbers come from?

    cout << "For upper left corner ..." << endl;

    TestLineSamp(cam, 1.0, 1.0);

    cout << "For upper right corner ..." << endl;

    TestLineSamp(cam, 257.0, 1.0);

    cout << "For lower left corner ..." << endl;

    TestLineSamp(cam, 1.0, 257.0);

    cout << "For lower right corner ..." << endl;

    TestLineSamp(cam, 257.0, 257.0);

    double samp = 257.0/2.0;

    double line = 257.0/2.0;

    cout << "For center pixel position ..." << endl;

    if (!cam->SetImage(samp, line)) {

      cout << "ERROR" << endl;

//      return 0;

    }

    if (abs(cam->UniversalLatitude() - knownLat) < 6E-14) {

      cout << "Latitude OK" << endl;

    }

    else {

      cout << setprecision(16) << "Latitude off by: " << cam->UniversalLatitude() - knownLat << endl;

    }

    if (abs(cam->UniversalLongitude() - knownLon) < 6E-14) {

      cout << "Longitude OK" << endl;

    }

    else {

      cout << setprecision(16) << "Longitude off by: " << cam->UniversalLongitude() - knownLon << endl;

    }

  }

  catch (IException &e) {

    e.print();

  }

}

void TestLineSamp(Camera *cam, double samp, double line) {

  bool success = cam->SetImage(samp, line);

  if (success) {

    success = cam->SetUniversalGround(cam->UniversalLatitude(), cam->UniversalLongitude());

  }

  if (success) {

    double deltaSamp = samp - cam->Sample();

    double deltaLine = line - cam->Line();

    if (fabs(deltaSamp) < 0.001) deltaSamp = 0.0;

    if (fabs(deltaLine) < 0.001) deltaLine = 0.0;

    cout << "DeltaSample = " << deltaSamp << endl;

    cout << "DeltaLine = " << deltaLine << endl << endl;

  }

  else {

    cout << "DeltaSample = ERROR" << endl;

    cout << "DeltaLine = ERROR" << endl << endl;

  }

}