VariableLineScanCameraDetectorMap.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#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];
  }
 
}