VariableLineScanCameraDetectorMap.cpp

Go to the documentation of this file.

#include "VariableLineScanCameraDetectorMap.h"
#include "CameraFocalPlaneMap.h"
#include "iTime.h"

namespace Isis {
  VariableLineScanCameraDetectorMap::VariableLineScanCameraDetectorMap(
      Camera *parent, 
      std::vector< LineRateChange > &lineRates) :
      LineScanCameraDetectorMap(parent, 
                                lineRates[0].GetStartEt(), 
                                lineRates[0].GetLineScanRate()), 
                                p_lineRates(lineRates) {
  }


  VariableLineScanCameraDetectorMap::~VariableLineScanCameraDetectorMap() {
  }


  bool VariableLineScanCameraDetectorMap::SetDetector(const double sample,
                                                      const double line) {
    // Use the parent SetDetector for the sample, which should work fine
    if (!CameraDetectorMap::SetDetector(sample, line)) {
      return false;
    }

    // currEt is our known et time
    double currEt = p_camera->time().Et();
    int rateIndex = p_lineRates.size() - 1;

    while (rateIndex >= 0 && currEt < p_lineRates[rateIndex].GetStartEt() - 0.5) {
      rateIndex --;
    }

    if (rateIndex < 0) {
      return false;
    }

    int rateStartLine = p_lineRates[rateIndex].GetStartLine();
    double rateStartEt = p_lineRates[rateIndex].GetStartEt();
    double rate = p_lineRates[rateIndex].GetLineScanRate();

    double etDiff = currEt - rateStartEt;
    p_parentLine = etDiff / rate + (rateStartLine - 0.5);

    //std::cout << "p_parentLine = " << p_parentLine << " = " << etDiff << "/" << rate << " + " << rateStartLine << std::endl;

    SetLineRate(rate);

    return true;
  }


  bool VariableLineScanCameraDetectorMap::SetParent(const double sample, 
                                                    const double line) {
    return SetParent(sample, line, 0.0); 
  }


  bool VariableLineScanCameraDetectorMap::SetParent(const double sample, 
                                                    const double line, 
                                                    const double deltaT) {
    if (!CameraDetectorMap::SetParent(sample, line, deltaT)) {
      return false;
    } //check to make sure we're not doubling the "DeltaT" 

    p_detectorLine = p_camera->FocalPlaneMap()->DetectorLineOffset();

    int rateIndex = p_lineRates.size() - 1;

    while (rateIndex >= 0 && line < p_lineRates[rateIndex].GetStartLine() - 0.5) {
      rateIndex --;
    }

    if (rateIndex < 0) {
      return false;
    }

    int rateStartLine = p_lineRates[rateIndex].GetStartLine();
    double rateStartEt = p_lineRates[rateIndex].GetStartEt();
    double rate = p_lineRates[rateIndex].GetLineScanRate();

    double et = rateStartEt + (line - (rateStartLine - 0.5)) * rate;
    //printf("et = %.8f = %.8f + (%.3f - %i) * %.8f\n", et, rateStartEt, line, rateStartLine, rate);

    SetLineRate(rate);

    p_camera->setTime(et + deltaT);

    return true;
  }


  double VariableLineScanCameraDetectorMap::exposureDuration(const double sample,
                                                             const double line,
                                                             const int band) const {
    return lineRate(line).GetLineScanRate();
  }


  LineRateChange &VariableLineScanCameraDetectorMap::lineRate(const double line) const{
    int rateIndex = p_lineRates.size() - 1;

    while (rateIndex >= 0 && line < p_lineRates[rateIndex].GetStartLine() - 0.5) {
      rateIndex --;
    }

    if (rateIndex < 0) {
      //This may not be a good idea
      //  SetParent used to return false here
      //  exposureDuration used to throw an exception here
      rateIndex = 0;
    }

    return p_lineRates[rateIndex];
  }

}