IsisShape.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
#include "IsisShape.h"
 
#include <QVector>
 
#include "ShapeModel.h"
#include "SurfacePoint.h"
 
using namespace std;
 
namespace Isis {
  IsisShape::IsisShape (ShapeModel* shape) {
    m_shape = shape;
  }
 
 
  SensorUtilities::Intersection IsisShape::intersect(const SensorUtilities::Vec &sensorPos, const SensorUtilities::Vec &lookVec, bool computeLocalNormal) {
    vector<double> pos = {sensorPos.x, sensorPos.y, sensorPos.z};
    vector<double> look = {lookVec.x, lookVec.y, lookVec.z};
    m_shape->intersectSurface(pos, look);
    SurfacePoint *intersection = m_shape->surfaceIntersection();
    if (computeLocalNormal && !m_shape->isDEM()) {
      QVector<double *> unusedNeighborPoints(4);
      double origin[3] = {0, 0, 0};
      unusedNeighborPoints.fill(origin);
      m_shape->calculateLocalNormal(unusedNeighborPoints);
    }
    else {
      m_shape->calculateSurfaceNormal();
    }
    vector<double> normal = m_shape->normal();
    return {
          {intersection->GetX().meters(), intersection->GetY().meters(), intersection->GetZ().meters()},
          {normal[0], normal[1], normal[2]}};
  }
}