void ossimSpectraboticsRedEdgeModel::imagingRay(const ossimDpt& image_point,
ossimEcefRay& image_ray) const
{
if(traceDebug())
{
ossimNotify(ossimNotifyLevel_DEBUG) << "ossimSpectraboticsRedEdgeModel::imagingRay: ..... entered" << std::endl;
}
ossimDpt film (image_point.x-m_calibratedCenter.x,
m_calibratedCenter.y - image_point.y); //- theRefImgPt);
// ossimDpt film (image_point-m_calibratedCenter); //- theRefImgPt);
if(m_lensDistortion.valid())
{
ossimDpt tempFilm(film.x/m_norm, film.y/m_norm);
ossimDpt filmOut;
m_lensDistortion->undistort(tempFilm, filmOut);
film.x = filmOut.x*m_norm;
film.y = filmOut.y*m_norm;
}
film.x *= m_pixelSize.x; // pixel size on the film
film.y *= m_pixelSize.y; // pixel size on the film
ossimColumnVector3d cam_ray_dir (film.x,
film.y,
-m_focalLength);
ossimEcefVector ecf_ray_dir (m_compositeMatrix*cam_ray_dir);
ecf_ray_dir = ecf_ray_dir*(1.0/ecf_ray_dir.magnitude());
image_ray.setOrigin(m_adjEcefPlatformPosition);
image_ray.setDirection(ecf_ray_dir);
}
void ossimBuckeyeSensor::imagingRay(const ossimDpt& image_point,
ossimEcefRay& image_ray) const
{
if(traceDebug()) ossimNotify(ossimNotifyLevel_DEBUG) << "ossimBuckeyeSensor::imagingRay: ..... entered" << std::endl;
ossimDpt f1 ((image_point) - theRefImgPt);
f1.x *= thePixelSize.x;
f1.y *= -thePixelSize.y;
ossimDpt film (f1 - thePrincipalPoint);
if(traceDebug())
{
ossimNotify(ossimNotifyLevel_DEBUG) << "pixel size = " << thePixelSize << std::endl;
ossimNotify(ossimNotifyLevel_DEBUG) << "principal pt = " << thePrincipalPoint << std::endl;
ossimNotify(ossimNotifyLevel_DEBUG) << "film pt = " << film << std::endl;
}
if (theLensDistortion.valid())
{
ossimDpt filmOut;
theLensDistortion->undistort(film, filmOut);
film = filmOut;
}
ossimColumnVector3d cam_ray_dir (film.x,
film.y,
-theFocalLength);
ossimEcefVector ecf_ray_dir (theCompositeMatrix*cam_ray_dir);
ecf_ray_dir = ecf_ray_dir*(1.0/ecf_ray_dir.magnitude());
image_ray.setOrigin(theAdjEcefPlatformPosition);
image_ray.setDirection(ecf_ray_dir);
if(traceDebug()) ossimNotify(ossimNotifyLevel_DEBUG) << "ossimBuckeyeSensor::imagingRay: ..... leaving" << std::endl;
}
void rspfApplanixEcefModel::imagingRay(const rspfDpt& image_point,
rspfEcefRay& image_ray) const
{
// if(traceDebug())
// {
// rspfNotify(rspfNotifyLevel_DEBUG) << "rspfApplanixEcefModel::imagingRay: ..... entered" << std::endl;
// }
rspfDpt f1 ((image_point) - theRefImgPt);
f1.x *= thePixelSize.x;
f1.y *= -thePixelSize.y;
rspfDpt film (f1 - thePrincipalPoint);
// if(traceDebug())
// {
// rspfNotify(rspfNotifyLevel_DEBUG) << "pixel size = " << thePixelSize << std::endl;
// rspfNotify(rspfNotifyLevel_DEBUG) << "principal pt = " << thePrincipalPoint << std::endl;
// rspfNotify(rspfNotifyLevel_DEBUG) << "film pt = " << film << std::endl;
// }
if (theLensDistortion.valid())
{
rspfDpt filmOut;
theLensDistortion->undistort(film, filmOut);
film = filmOut;
}
rspfColumnVector3d cam_ray_dir (film.x,
film.y,
-theFocalLength);
rspfEcefVector ecf_ray_dir (theCompositeMatrix*cam_ray_dir);
ecf_ray_dir = ecf_ray_dir*(1.0/ecf_ray_dir.magnitude());
image_ray.setOrigin(theAdjEcefPlatformPosition);
image_ray.setDirection(ecf_ray_dir);
// if(traceDebug())
// {
// rspfNotify(rspfNotifyLevel_DEBUG) << "rspfApplanixEcefModel::imagingRay: ..... leaving" << std::endl;
// }
}
void rspfBuckeyeSensor::imagingRay(const rspfDpt& image_point,
rspfEcefRay& image_ray) const
{
rspfDpt f1 ((image_point) - theRefImgPt);
f1.x *= m_pixelSize.x;
f1.y *= -m_pixelSize.y;
rspfDpt film (f1 - m_principalPoint);
if (m_lensDistortion.valid())
{
rspfDpt filmOut;
m_lensDistortion->undistort(film, filmOut);
film = filmOut;
}
rspfColumnVector3d cam_ray_dir (film.x,
film.y,
-(m_focalLength+computeParameterOffset(6)));
rspfEcefVector ecf_ray_dir (m_compositeMatrix*cam_ray_dir);
ecf_ray_dir = ecf_ray_dir*(1.0/ecf_ray_dir.magnitude());
image_ray.setOrigin(m_ecefPlatformPosition);
image_ray.setDirection(ecf_ray_dir);
}
