#include <ossimRS1SarModel.h>

Inheritance diagram for ossimRS1SarModel:

Public Types
enum	ImagingMode { UNKNOWN_MODE = 0, SCN, SCW, SGC, SGF, SGX, SLC, SPG, SSG, RAW, ERS }

enum	DirectionFlag { UNKNOWN_DIRECTION = 0, ASCENDING, DESCENDING }

Public Types inherited from ossimSensorModel
enum	CovMatStatus { COV_INVALID = 0, COV_PARTIAL = 1, COV_FULL = 2 }

enum	DeriveMode { OBS_INIT =-99, EVALUATE =-98, P_WRT_X = -1, P_WRT_Y = -2, P_WRT_Z = -3 }

Public Member Functions
	ossimRS1SarModel ()

	ossimRS1SarModel (const ossimFilename &imageDir)

virtual	~ossimRS1SarModel ()

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=NULL) const
	Fulfills ossimObject base-class pure virtuals. More...

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=NULL)
	Fulfills ossimObject base-class pure virtuals. More...

virtual void	lineSampleHeightToWorld (const ossimDpt &imagePt, const double &heightAboveEllipsoid, ossimGpt &worldPt) const
	Establishes geographic 3D point given image line, sample and ellipsoid height. More...

virtual void	imagingRay (const ossimDpt &image_point, ossimEcefRay &image_ray) const
	Given an image point, returns a ray originating at some arbitrarily high point (in this model at the sensor position) and pointing towards the target. More...

virtual bool	useForward () const

virtual ossimObject *	dup () const
	image to ground faster More...

Public Member Functions inherited from ossimSensorModel
	ossimSensorModel ()

	ossimSensorModel (const ossimSensorModel &copy_this)

	ossimSensorModel (const ossimKeywordlist &geom_kwl)

const ossimSensorModel &	operator= (const ossimSensorModel &rhs)
	assignment operator More...

virtual ossimObject *	getBaseObject ()

virtual const ossimObject *	getBaseObject () const

virtual ossimGpt	origin () const

virtual ossimDpt	getMetersPerPixel () const

virtual const double &	getNominalPosError () const
	Returns the estimated Absolute horizontal position error (CE90) of the sensor model. More...

virtual const double &	getRelativePosError () const
	Returns the estimated RELATIVE horizontal position error (CE90) of the sensor model. More...

virtual void	setNominalPosError (const double &ce90)
	Assigns the absolute image position error uncertainty (abs CE90) More...

virtual void	setRelativePosError (const double &ce90)
	Assigns the relative image position error uncertainty (rel CE90) More...

virtual void	lineSampleToWorld (const ossimDpt &image_point, ossimGpt &world_point) const

virtual void	worldToLineSample (const ossimGpt &world_point, ossimDpt &image_point) const

virtual std::ostream &	print (std::ostream &out) const

void	setRefImgPt (const ossimDpt &pt)
	Sets the center line sampe of the image. More...

void	setRefGndPt (const ossimGpt &pt)
	Sets the center latitude, longitude, height of the image. More...

void	setImageRect (const ossimDrect &imageRect)

void	setGroundRect (const ossimGpt &ul, const ossimGpt &ur, const ossimGpt &lr, const ossimGpt &ll)

ossimDpt	imageSize () const

void	setImageSize (const ossimDpt &size)

virtual void	adjustableParametersChanged ()

virtual void	updateModel ()

virtual bool	insideImage (const ossimDpt &p) const

virtual bool	operator== (const ossimProjection &proj) const

const ossimString &	getImageID () const
	Access methods: More...

const ossimDrect &	getImageClipRect () const

virtual ossim_uint32	degreesOfFreedom () const

virtual bool	needsInitialState () const
	needsInitialState() More...

virtual double	optimizeFit (const ossimTieGptSet &tieSet, double *targetVariance=0)

virtual ossimDpt	getForwardDeriv (int parmIdx, const ossimGpt &gpos, double hdelta=1e-11)

virtual ossimGpt	getInverseDeriv (int parmIdx, const ossimDpt &ipos, double hdelta=1e-11)

virtual ossimSensorModel::CovMatStatus	getObsCovMat (const ossimDpt &ipos, NEWMAT::SymmetricMatrix &Cov, const ossim_float64 defPointingSigma=0.5) const
	Gives 2X2 covariance matrix of observations. More...

virtual bool	isAffectedByElevation () const
	Implementation of pure virtual ossimProjection::isAffectedByElevation method. More...

void	computeGsd ()
	This method computes the ground sample distance(gsd) and sets class attributes theGSD and theMeanGSD by doing a lineSampleHeightToWorld on four points and calculating the distance from them. More...

virtual bool	getImageGeometry (const ossimString &, const ossimString &, ossimKeywordlist &) const
	Extracts geometry info from a non-ossim key,value pair to an ossim keyword list. More...

Public Member Functions inherited from ossimProjection
	ossimProjection ()

virtual	~ossimProjection ()

virtual ossimDpt	forward (const ossimGpt &wp) const

virtual ossimGpt	inverse (const ossimDpt &pp) const

virtual void	getRoundTripError (const ossimDpt &imagePoint, ossimDpt &errorResult) const

virtual void	getRoundTripError (const ossimGpt &groundPoint, ossimDpt &errorResult) const

virtual void	getGroundClipPoints (ossimGeoPolygon &gpts) const

virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const

virtual bool	operator!= (const ossimProjection &projection) const

Public Member Functions inherited from ossimObject
	ossimObject ()

virtual	~ossimObject ()

virtual ossimString	getShortName () const

virtual ossimString	getLongName () const

virtual ossimString	getDescription () const

virtual ossimString	getClassName () const

virtual RTTItypeid	getType () const

virtual bool	canCastTo (ossimObject *obj) const

virtual bool	canCastTo (const RTTItypeid &id) const

virtual bool	canCastTo (const ossimString &parentClassName) const

virtual void	accept (ossimVisitor &visitor)

Public Member Functions inherited from ossimReferenced
	ossimReferenced ()

	ossimReferenced (const ossimReferenced &)

ossimReferenced &	operator= (const ossimReferenced &)

void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...

void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

int	referenceCount () const

Public Member Functions inherited from ossimErrorStatusInterface
	ossimErrorStatusInterface ()

virtual	~ossimErrorStatusInterface ()

virtual ossimErrorCode	getErrorStatus () const

virtual ossimString	getErrorStatusString () const

virtual void	setErrorStatus (ossimErrorCode error_status) const

virtual void	setErrorStatus () const

virtual void	clearErrorStatus () const

bool	hasError () const

Public Member Functions inherited from ossimOptimizableProjection
	ossimOptimizableProjection ()

	ossimOptimizableProjection (const ossimOptimizableProjection &source)

virtual	~ossimOptimizableProjection ()

virtual ossimOptimizableProjection &	operator= (const ossimOptimizableProjection &source)

virtual bool	setupOptimizer (const ossimString &setup)
	setupFromString() Derived classes should implement as needed. Initialize parameters needed for optimizeFit and degreesOfFreedom More...

Public Member Functions inherited from ossimAdjustableParameterInterface
	ossimAdjustableParameterInterface ()

	ossimAdjustableParameterInterface (const ossimAdjustableParameterInterface &rhs)

virtual	~ossimAdjustableParameterInterface ()

void	newAdjustment (ossim_uint32 numberOfParameters=0)

void	setAdjustmentDescription (const ossimString &description)

ossimString	getAdjustmentDescription () const

ossimString	getAdjustmentDescription (ossim_uint32 adjustmentIdx) const
	Returns adjustmentDescription of specific adjustmentInfo. More...

void	setCurrentAdjustment (ossim_uint32 adjustmentIndex, bool notify=false)

bool	setCurrentAdjustment (const ossimString &description, bool notify=false)
	Sets the current adjustment to the adjustment with a matching description. More...

void	eraseAdjustment (bool notify)

void	eraseAdjustment (ossim_uint32 idx, bool notify)

virtual void	initAdjustableParameters ()

void	resetAdjustableParameters (bool notify=false)

void	copyAdjustment (ossim_uint32 idx, bool notify)

void	copyAdjustment (bool notify=false)

void	keepAdjustment (ossim_uint32 idx, bool createCopy)

virtual void	keepAdjustment (bool createCopy=true)

const ossimAdjustableParameterInterface &	operator= (const ossimAdjustableParameterInterface &rhs)

void	removeAllAdjustments ()

ossim_uint32	getNumberOfAdjustableParameters () const

double	getAdjustableParameter (ossim_uint32 idx) const

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, bool notify=false)

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, double sigma, bool notify=false)

double	getParameterSigma (ossim_uint32 idx) const

void	setParameterSigma (ossim_uint32 idx, double value, bool notify=false)

ossimUnitType	getParameterUnit (ossim_uint32 idx) const

void	setParameterUnit (ossim_uint32 idx, ossimUnitType unit)

void	setParameterUnit (ossim_uint32 idx, const ossimString &unit)

void	setParameterCenter (ossim_uint32 idx, double center, bool notify=false)

double	getParameterCenter (ossim_uint32 idx) const

double	computeParameterOffset (ossim_uint32 idx) const

void	setParameterOffset (ossim_uint32 idx, ossim_float64 value, bool notify=false)

ossimString	getParameterDescription (ossim_uint32 idx) const

void	setParameterDescription (ossim_uint32 idx, const ossimString &descrption)

ossim_int32	findParameterIdxGivenDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

ossim_int32	findParameterIdxContainingDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

bool	isParameterLocked (ossim_uint32 idx) const

void	setParameterLockFlag (ossim_uint32 idxParam, bool flag)

bool	getParameterLockFlag (ossim_uint32 idx) const

void	lockAllParametersCurrentAdjustment ()

void	unlockAllParametersCurrentAdjustment ()

void	lockAllParameters (ossim_uint32 idxAdjustment)

void	unlockAllParameters (ossim_uint32 idxAdjustment)

void	resizeAdjustableParameterArray (ossim_uint32 numberOfParameters)

void	setAdjustment (const ossimAdjustmentInfo &adj, bool notify=false)

void	setAdjustment (ossim_uint32 idx, const ossimAdjustmentInfo &adj, bool notify=false)

void	addAdjustment (const ossimAdjustmentInfo &adj, bool notify)

void	getAdjustment (ossimAdjustmentInfo &adj) const

void	getAdjustment (ossim_uint32 idx, ossimAdjustmentInfo &adj) const

ossim_uint32	getNumberOfAdjustments () const

ossim_uint32	getCurrentAdjustmentIdx () const

void	setDirtyFlag (bool flag=true)

void	setAllDirtyFlag (bool flag=true)

bool	hasDirtyAdjustments () const

virtual void	saveCurrentAdjustmentOnly (ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))
	Saves the current active adjustment to the KWL file. More...

bool	saveAdjustments (ossimKeywordlist &kwl, const ossimString &prefix=ossimString("")) const
	Save all adjustments to the KWL file. More...

bool	loadAdjustments (const ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))

std::ostream &	print (std::ostream &out) const
	Dumps the currently active adjustment to ostream. More...

Protected Types
enum	ADJUSTABLE_PARAM_INDEXES { INTRACK_OFFSET, CRTRACK_OFFSET, RADIAL_OFFSET, LINE_SCALE, SKEW, ORIENTATION, NUM_ADJUSTABLE_PARAMS }

Protected Member Functions
void	setImagingMode (char *modeStr)

void	initFromCeos (const ossimFilename &dataDir)

void	initAdjParms ()

void	establishEphemeris ()

void	eciToEcfXform (const double &julianDay, NEWMAT::Matrix &xform) const

void	establishOrpInterp ()

void	establishOrpGrid ()

void	establishVehicleSpace ()

void	interpolatedScanORP (const ossimDpt &orp, ossimEcefPoint &orp_ecf) const

void	deallocateMemory ()

Protected Member Functions inherited from ossimSensorModel
virtual	~ossimSensorModel ()

virtual ossimDpt	extrapolate (const ossimGpt &gp) const

virtual ossimGpt	extrapolate (const ossimDpt &ip, const double &height=ossim::nan()) const

void	buildNormalEquation (const ossimTieGptSet &tieSet, NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &residue, NEWMAT::ColumnVector &projResidue, double pstep_scale)

NEWMAT::ColumnVector	getResidue (const ossimTieGptSet &tieSet)

NEWMAT::ColumnVector	solveLeastSquares (NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &r) const

NEWMAT::Matrix	invert (const NEWMAT::Matrix &m) const
	stable invert stolen from ossimRpcSolver More...

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Member Functions inherited from ossimAdjustableParameterInterface
bool	paramChanged (ossim_uint32 param_idx) const
	Returns true if specified parameter has been modified since last setAllChangeFlag(false) call. More...

void	setAllChangeFlags (bool areChanged)
	Sets all the change flags to the boolean indicated to indicate parameters are changed (TRUE) or not (FALSE). More...

void	initChangeFlags ()
	Initializes the change flags to TRUE. More...

Protected Attributes
ossimRefPtr< ossimCeosData >	theCeosData

ossimRefPtr< ossimLagrangeInterpolator >	theArpPosInterp

ossimRefPtr< ossimLagrangeInterpolator >	theArpVelInterp

ossimRefPtr< ossimLagrangeInterpolator >	theLocalOrpInterp

ImagingMode	theImagingMode

double	theIllumAzimuth

double	theIllumElevation

ossimEcefPoint	theORP

double	theRefHeight

double	theGHA

double	theEphFirstSampTime

double	theSinOrientation

double	theCosOrientation

double	theSinSkew

ossimLsrSpace	theVehicleSpace

ossimEcefVector	thePosCorrection

int	theFirstLineDay

double	theFirstLineTime

double	theTimePerLine

DirectionFlag	theDirectionFlag

double	theSrGrCoeff [6]

ossimDpt	thePixelSpacing

ossimDblGrid	theLatGrid

ossimDblGrid	theLonGrid

double	theInTrackOffset

double	theCrTrackOffset

double	theRadialOffset

double	theLineScale

double	theSkew

double	theOrientation

Protected Attributes inherited from ossimSensorModel
ossimIpt	theImageSize

ossimDpt	theSubImageOffset

ossimString	theImageID

ossimString	theSensorID

ossimDpt	theGSD

ossim_float64	theMeanGSD

ossimGpt	theRefGndPt

ossimDpt	theRefImgPt

ossimPolygon	theBoundGndPolygon

ossimDrect	theImageClipRect

ossim_float64	theRelPosError

ossim_float64	theNominalPosError

ossimDpt	theParWRTx
	Partials for current point. More...

ossimDpt	theParWRTy

ossimDpt	theParWRTz

ossimDpt	theObs
	Observations & residuals for current point. More...

ossimDpt	theResid

ossimRefPtr< ossimProjection >	theSeedFunction
	Used as an initial guess for iterative solutions and a guess for points outside the support bounds. More...

bool	theExtrapolateImageFlag

bool	theExtrapolateGroundFlag

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Additional Inherited Members
Static Public Member Functions inherited from ossimSensorModel
static void	writeGeomTemplate (ostream &os)

Detailed Description

Definition at line 37 of file ossimRS1SarModel.h.

Member Enumeration Documentation

◆ ADJUSTABLE_PARAM_INDEXES

enum ossimRS1SarModel::ADJUSTABLE_PARAM_INDEXES

protected

Enumerator
INTRACK_OFFSET
CRTRACK_OFFSET
RADIAL_OFFSET
LINE_SCALE
SKEW
ORIENTATION
NUM_ADJUSTABLE_PARAMS

Definition at line 142 of file ossimRS1SarModel.h.

    {
       INTRACK_OFFSET,
       CRTRACK_OFFSET,
       RADIAL_OFFSET,
       LINE_SCALE,
       SKEW,
       ORIENTATION,
       NUM_ADJUSTABLE_PARAMS  // not a parameter
    };

◆ DirectionFlag

enum ossimRS1SarModel::DirectionFlag

Enumerator
UNKNOWN_DIRECTION
ASCENDING
DESCENDING

Definition at line 60 of file ossimRS1SarModel.h.

    {
       UNKNOWN_DIRECTION = 0,
       ASCENDING,
       DESCENDING
    };

◆ ImagingMode

enum ossimRS1SarModel::ImagingMode

Enumerator
UNKNOWN_MODE
SCN
SCW
SGC
SGF
SGX
SLC
SPG
SSG
RAW
ERS

Definition at line 45 of file ossimRS1SarModel.h.

    {
       UNKNOWN_MODE = 0,
       SCN,     // ScanSAR Narrow Beam
       SCW,     // ScanSAR Wide Beam
       SGC,     // SAR Georeferenced Coarse Resolution
       SGF,     // SAR Georeferenced Fine Resolution
       SGX,     // SAR Georeferenced Extra Fine Resolution
       SLC,     // Single Look Complex
       SPG,     // SAR Precision Geocoded
       SSG,     // SAR Systematically Geocoded
       RAW,
       ERS
    };

Constructor & Destructor Documentation

◆ ossimRS1SarModel() [1/2]

ossimRS1SarModel::ossimRS1SarModel ( )

Definition at line 48 of file ossimRS1SarModel.cpp.

References ossimErrorStatusInterface::setErrorStatus().

 :  theCeosData        (0),
    theImagingMode     (UNKNOWN_MODE)
 {
    setErrorStatus();
 }

◆ ossimRS1SarModel() [2/2]

ossimRS1SarModel::ossimRS1SarModel ( const ossimFilename & imageDir )

Definition at line 59 of file ossimRS1SarModel.cpp.

    :
       theCeosData        (0),
       theImagingMode     (UNKNOWN_MODE)
 {
    static const char MODULE[] = "Constructor ossimRS1SarModel(ossimFilename)";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    clearErrorStatus();
 
    // Parse the CEOS data files:
    ossimFilename dataDirName (imageFile.path());
    initFromCeos(dataDirName);
    if(getErrorStatus() != ossimErrorCodes::OSSIM_OK)
    {
       return;
    }
    if (theImagingMode == UNKNOWN_MODE) 
       return;
 
    // Parse the image data file for the local ORPs interpolator:
    if ((theImagingMode == SCN) || (theImagingMode == SCW))
       establishOrpGrid();
    else
       establishOrpInterp();
 
    // Establish the ephemeris interpolator:
    establishEphemeris();
 
    initAdjParms();
    establishVehicleSpace();
 
    if (traceDebug())  CLOG << "returning..." << endl;
    return;
 }

◆ ~ossimRS1SarModel()

ossimRS1SarModel::~ossimRS1SarModel ( )

virtual

Definition at line 99 of file ossimRS1SarModel.cpp.

References deallocateMemory().

 {
    deallocateMemory();
 }

Member Function Documentation

◆ deallocateMemory()

void ossimRS1SarModel::deallocateMemory ( )

protected

Definition at line 1072 of file ossimRS1SarModel.cpp.

Referenced by ~ossimRS1SarModel().

 {
    static const char MODULE[] = "ossimRS1SarModel::deallocateMemory()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    theArpPosInterp = 0;
    theArpVelInterp = 0;
    theLocalOrpInterp = 0;
    theCeosData        = 0;
    theLatGrid.clear();
    theLonGrid.clear();
       
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ dup()

virtual ossimObject* ossimRS1SarModel::dup ( ) const

inlinevirtual

image to ground faster

Returns pointer to a new instance, copy of this.

Implements ossimProjection.

Definition at line 87 of file ossimRS1SarModel.h.

87 { return 0; } // TBR

◆ eciToEcfXform()

void ossimRS1SarModel::eciToEcfXform	(	const double &	julianDay,
		NEWMAT::Matrix &	xform
	)		const

protected

Definition at line 672 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::eciToEcfXform(acq_time)";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    // Determine time elapsed since ephemeris first sample time:
    double delta_t = acq_time - theEphFirstSampTime;
 
    // Compute GHA at time of this image line (in degrees):
    double gha = theGHA + delta_t*DEG_PER_SEC;
    if (gha >= 360.0)
       gha -=360.0;
    else if (gha < 0.0)
       gha += 360.0;
    
    // Establish rotation:
    double cos_gha = ossim::cosd(gha);
    double sin_gha = ossim::sind(gha);
 
    // Establish rotation matrix:
    xform = ossimMatrix3x3::create(cos_gha,  sin_gha,  0.0,
                                  -sin_gha,  cos_gha,  0.0,
                                       0.0,      0.0,  1.0);
    if (traceDebug())
    {
       CLOG << "DEBUG -- "
            << "\n\t acq_time = " << acq_time
            << "\n\t delta_t  = " << delta_t
            << "\n\t gha      = " << gha << endl;
    }
 
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ establishEphemeris()

void ossimRS1SarModel::establishEphemeris ( )

protected

Definition at line 579 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::establishEphemeris()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    char buf[23];
    
    // Establish pointer to platform position data record:
    const pos_data_rec* pdr = theCeosData->posDataRec();
    
    // Establish the Greenwich mean hour angle at time of first sample:
    char gha_str [23];
    strncpy(gha_str, pdr->hr_angle, 22);
    gha_str[22] = '\0';
    theGHA = atof(gha_str);
 
    // Fetch number of points in record and allocate memory:
    strncpy(buf, pdr->ndata, 4);
    buf[4] = '\0';
    int numPoints = atoi(buf);
 
    // Fetch sampling period:
    strncpy(buf, pdr->data_int, 22);
    buf[22] = '\0';
    double sampling_period = atof(buf);
 
    // Establish the first sample time's day:
    strncpy(buf, pdr->gmt_day, 4);
    buf[4] = '\0';
    int first_day =  atoi(buf);
    double time_offset = (double) (first_day - theFirstLineDay)*SEC_PER_DAY;
 
    // Establish the first sample time in seconds from beginning of day:
    strncpy(buf, pdr->gmt_sec, 22);
    buf[22] = '\0';
    theEphFirstSampTime =  atof(buf) + time_offset;
 
    int i, j;
    double sample_time = theEphFirstSampTime;
    buf[22] = '\0';
    NEWMAT::ColumnVector eciArpPos(3);
    NEWMAT::ColumnVector eciArpVel(3);
    NEWMAT::ColumnVector ecfArpPos(3);
    NEWMAT::ColumnVector ecfArpVel(3);
    NEWMAT::ColumnVector earthVel(3);
    NEWMAT::Matrix xform;
    theArpPosInterp = new ossimLagrangeInterpolator;
    theArpVelInterp = new ossimLagrangeInterpolator;
 
    // Loop over each point, converting them to vectors:
    for (i=0; i<numPoints; i++)
    {
       for (j=0; j<3; j++)
       {
          strncpy(buf, pdr->pos_vect[i].pos[j], 22);
          eciArpPos[j] = atof(buf);
          strncpy(buf, pdr->pos_vect[i].vel[j], 22);
          eciArpVel[j] = atof(buf)/1000.0;
       }
 
       // Convert pos and vel vectors from ECI to ECF:
       eciToEcfXform(sample_time, xform);
       ecfArpPos = xform*eciArpPos;
       ecfArpVel = xform*eciArpVel;
       theArpPosInterp->addData(sample_time, ecfArpPos);
 
       // Correct velocity for earth rotation:
       earthVel[0] = -ecfArpPos[1]*EARTH_ANGULAR_VELOCITY;
       earthVel[1] =  ecfArpPos[0]*EARTH_ANGULAR_VELOCITY;
       earthVel[2] =  0.0;
       ecfArpVel = ecfArpVel - earthVel;
       theArpVelInterp->addData(sample_time, ecfArpVel);
 
       sample_time += sampling_period;
    }
 
    if (traceDebug())
    {
       CLOG << "DEBUG -- "
            << "\n\t theGHA:              " << theGHA 
            << "\n\t numPoints:           " << numPoints 
            << "\n\t sampling_period:     " << sampling_period 
            << "\n\t theEphFirstSampTime: " << theEphFirstSampTime << endl;
    }
 
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ establishOrpGrid()

void ossimRS1SarModel::establishOrpGrid ( )

protected

Definition at line 928 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::establishOrpGrid()";
    if (traceDebug())  CLOG << "entering..." << endl;
    
    if (traceDebug())
    {
       CLOG << "DEBUG -- " << endl;
    }
 
    static const int      NUM_GRID_POINTS_U  = 11;
 
    desc_rec        descRec;
    pdr_prefix_rec  prefix;
    int             sizeOfDescRec = sizeof(descRec);
    int             sizeOfPrefRec = sizeof(prefix);
    int             headerSize    = sizeOfDescRec + sizeOfPrefRec;
    char            buf[] = "123456";
    ossimDpt gridSize ((double) NUM_GRID_POINTS_U, 3.0);
    ossimDpt cellSize (theImageSize.line/gridSize.u, theImageSize.samp/gridSize.v);
 
    // Because we are doing integer arithmetic, it will be necessary to
    // consider the last grid cells as being of a slightly different size:
    //theLastGridCellSize.line = theImageSize.line - 1.0 -
    //                           (theGridSize.u-2.0)*theGridCellSize.u;
    //theLastGridCellSize.samp = theImageSize.samp - 1.0 -
    //                           (theGridSize.v-2.0)*theGridCellSize.v;
    
    // Open the image data file given the image directory name and seek to first
    // SAR data record:
    FILE* fptr = fopen(theCeosData->imageFile().chars(), "r");
    if (!fptr)
    {
       CLOG << "ERROR: Could not open data file <" << theCeosData->imageFile()
            << ">" << endl;
       setErrorStatus();
       return;
    }
    fseek(fptr, theCeosData->imopDescRec()->desc.length, SEEK_SET);
 
    // Declare some constants and variables used in loop:
    int last_line    = (int) theImageSize.line - 1;
    int delta_line   = (int) cellSize.u;
    int line_number  = 0;
 
    // Allocate memory for coarse grid:
    theLatGrid.initialize(gridSize, ossimDpt(0,0), ossimDpt(1,1));
    theLonGrid.initialize(gridSize, ossimDpt(0,0), ossimDpt(1,1));
 
    // Compute the number of SAR data bytes to skip to reach prefix of next
    // record of interest:
    strncpy(buf, theCeosData->imopDescRec()->l_dataset, 6);
    int record_size = atoi(buf);
    int test_line = (int) gridSize.u - 2;
    
    // Loop over each imaging line:
    for (int u=0; u<(int)gridSize.u; u++)
    {
       if (traceDebug())
       {
          clog << "\t Processing grid line number " << u << endl;
       }
 
       // Read the image line prefix info:
       fread(&descRec, sizeOfDescRec, 1, fptr);
       fread(&prefix,  sizeOfPrefRec,  1, fptr);
 
       // Verify correct line number being read:
       if ((prefix.line_num-1) != line_number)
       {
          CLOG << "\n\tERROR: Synchronization error reading image file. "
               << "Expected line number " << line_number+1
               << " but read line number " << prefix.line_num << "." << endl;
          setErrorStatus();
      fclose(fptr); 
      fptr = 0;   
          return;
       }
       
       // Latch the imaging time for the first line collected:
       if (u == 0)
       {
          theFirstLineDay  = prefix.acq_date.acq_day;
          theFirstLineTime = (double) prefix.acq_date.acq_msec/1000.0;
       }
       
       // Read the ground points for this line:
       theLatGrid.setNode(u, 0, (double) prefix.lat_first /1.0e6);
       theLatGrid.setNode(u, 1, (double) prefix.lat_mid   /1.0e6);
       theLatGrid.setNode(u, 2, (double) prefix.lat_last  /1.0e6);
       theLonGrid.setNode(u, 0, (double) prefix.long_first/1.0e6);
       theLonGrid.setNode(u, 1, (double) prefix.long_mid  /1.0e6);
       theLonGrid.setNode(u, 2, (double) prefix.long_last /1.0e6);
       
       // Update the line number for next record, insuring we don't go passed
       // the last line:
       if (u == test_line)
       {
          delta_line  = last_line - line_number;
          line_number = last_line;
       }
       else
       {
          line_number += delta_line;
       }
 
       //***
       // Advance to read the next prefix record of interest:
       //***
       fseek(fptr, (delta_line)*record_size-headerSize, SEEK_CUR);
    }
 
    // Compute time interval between each line, handle possible day rollover:
    double lastLineTime = (double) prefix.acq_date.acq_msec/1000.0;
    if ((lastLineTime-theFirstLineTime) > (SEC_PER_DAY/2.0))
    {
       lastLineTime -= SEC_PER_DAY;
    }
    else if ((theFirstLineTime-lastLineTime) > (SEC_PER_DAY/2.0))
    {
       lastLineTime += SEC_PER_DAY;
    }
    theTimePerLine = (lastLineTime - theFirstLineTime)/line_number;
 
    if (traceDebug())
    {
       clog << "\n\t line_number (last)   = " << line_number
            << "\n\t last_line (in image) = " << last_line
            << "\n\t theFirstLineTime     = " << theFirstLineTime
            << "\n\t lastLineTime         = " << lastLineTime
            << "\n\t theTimePerLine       = " << theTimePerLine
            << endl;
    }
    
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ establishOrpInterp()

void ossimRS1SarModel::establishOrpInterp ( )

protected

Definition at line 727 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::establishOrpInterp()";
    if (traceDebug())  CLOG << "entering..." << endl;
    
    if (traceDebug())
       CLOG << "DEBUG -- " << endl;
 
    static const int      NUM_DATA_POINTS  = 11;
 
    desc_rec        descRec;
    pdr_prefix_rec  prefix;
    int             sizeOfDescRec = sizeof(descRec);
    int             sizeOfPrefRec = sizeof(prefix);
    int             headerSize    = sizeOfDescRec + sizeOfPrefRec;
    char            buf[] = "123456";
    ossimGpt localOrp (0.0, 0.0, theRefHeight);
    std::vector<double> line_numbers_list;
    std::vector<NEWMAT::ColumnVector> orpVectorList; // X, Y, Z
 
    if (traceDebug())
    {
       CLOG << "DEBUG:"
            << "\nsizeOfDescRec:  " << sizeOfDescRec
            << "\nsizeOfPrefRec:  " << sizeOfPrefRec << endl;
    }
    
    // Open the image data file given the image directory name and seek to first
    // SAR data record:
    FILE* fptr = fopen(theCeosData->imageFile().chars(), "r");
    if (!fptr)
    {
       CLOG << "ERROR: Could not open data file <" << theCeosData->imageFile()
            << ">" << endl;
       setErrorStatus();
       return;
    }
    fseek(fptr, theCeosData->imopDescRec()->desc.length, SEEK_SET);
 
    // Declare some constants and variables used in loop:
    int last_line    = (int) theImageSize.line - 1;
    int delta_line   = (int)ceil(last_line/((double)NUM_DATA_POINTS-1.0));
    int line_number  = 0;
 
    // Compute the number of SAR data bytes to skip to reach prefix of next
    // record of interest:
    strncpy(buf, theCeosData->imopDescRec()->l_dataset, 6);
    int record_size = atoi(buf);
    
    // Loop over each imaging line:
    int index;
    theLocalOrpInterp = new ossimLagrangeInterpolator;
    for (index=0; index<NUM_DATA_POINTS; index++)
    {
       if (traceDebug())
          clog << "\t Processing ORP for line number " << line_number << endl;
 
       // Read the image line prefix info:
       fread(&descRec, sizeOfDescRec, 1, fptr);
       fread(&prefix,  sizeOfPrefRec,  1, fptr);
 
       // Verify correct line number being read:
       if ((prefix.line_num-1) != line_number)
       {
          CLOG << "\n\tERROR: Synchronization error reading image file. "
               << "Expected line number " << line_number+1
               << " but read line number " << prefix.line_num << "." << endl;
          setErrorStatus();
      fclose(fptr); 
      fptr = 0; 
          return;
       }
       
       // Latch the imaging time for the first line collected:
       if (index == 0)
       {
          theFirstLineDay  = prefix.acq_date.acq_day;
          theFirstLineTime = (double) prefix.acq_date.acq_msec/1000.0;
       }
       
       // Read the local ORP for this line and convert to ECF:
       if (theDirectionFlag == ASCENDING)
       {
          localOrp.lat    = (double) prefix.lat_first  / 1.0e6;
          localOrp.lon = (double) prefix.long_first / 1.0e6;
       }
       else
       {
          localOrp.lat    = (double) prefix.lat_last  / 1.0e6;
          localOrp.lon = (double) prefix.long_last / 1.0e6;
       }
 
       theLocalOrpInterp->addData(line_number, ossimEcefPoint(localOrp).toVector());
       
       // Update the line number for next record, insuring we don't go passed
       // the last line:
       line_number += delta_line;
       if (line_number > last_line)
       {
          delta_line -= line_number - last_line;  // adjust for fseek below
          line_number = last_line;
       }
 
       // Advance to read the next prefix record of interest:
       fseek(fptr, (delta_line)*record_size-headerSize, SEEK_CUR);
    }
 
    // Compute time interval between each line, handle possible day rollover:
    double lastLineTime = (double) prefix.acq_date.acq_msec/1000.0;
    if ((lastLineTime-theFirstLineTime) > (SEC_PER_DAY/2.0))
       lastLineTime -= SEC_PER_DAY;
    else if ((theFirstLineTime-lastLineTime) > (SEC_PER_DAY/2.0))
       lastLineTime += SEC_PER_DAY;
    theTimePerLine = (lastLineTime - theFirstLineTime)/line_number;
 
    if (traceDebug())
    {
       clog << "\n\t line_number (last)   = " << line_number
            << "\n\t last_line (in image) = " << last_line
            << "\n\t theFirstLineTime     = " << theFirstLineTime
            << "\n\t lastLineTime         = " << lastLineTime
            << "\n\t theTimePerLine       = " << theTimePerLine
            << endl;
    }
    
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ establishVehicleSpace()

void ossimRS1SarModel::establishVehicleSpace ( )

protected

Definition at line 859 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::establishVehicleSpace()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    // Obtain the imaging time when the vehicle is at the ARP:
    double line = theImageSize.line/2.0;
    double arpTime = theFirstLineTime + line*theTimePerLine;
 
    // Obtain the ephemeris at this time:
    NEWMAT::ColumnVector arpPos(3);
    theArpPosInterp->interpolate(arpTime, arpPos);
    NEWMAT::ColumnVector arpVel(3);
    theArpVelInterp->interpolate(arpTime, arpVel);
 
    // Convert ephemeris to ECF:
    NEWMAT::Matrix xform;
    eciToEcfXform(arpTime, xform);
    ossimEcefPoint ecfArpPos  (xform * arpPos);
    ossimEcefVector ecfArpVel (xform * arpVel);
 
    // Need to correct the velocity vector by the earth rotational velocity:
    ossimEcefVector earthVel(-ecfArpPos.y()*EARTH_ANGULAR_VELOCITY,
                              ecfArpPos.x()*EARTH_ANGULAR_VELOCITY,
                              0.0);
    ecfArpVel = ecfArpVel - earthVel;
 
    // Establish LSR space for intrack-crosstrack-radial at vehicle:
    ossimEcefVector intrackDir (ecfArpVel);
    ossimEcefVector crtrackDir (ecfArpPos.data().cross(intrackDir.data())); 
    theVehicleSpace = ossimLsrSpace (ecfArpPos, intrackDir, crtrackDir, 0);
 
    if (traceDebug())
    {
       CLOG << "DEBUG -- "
            << "\n\t theVehicleSpace: " << theVehicleSpace << endl;
    }
 
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ imagingRay()

void ossimRS1SarModel::imagingRay	(	const ossimDpt &	image_point,
		ossimEcefRay &	image_ray
	)		const

virtual

Given an image point, returns a ray originating at some arbitrarily high point (in this model at the sensor position) and pointing towards the target.

Reimplemented from ossimSensorModel.

Definition at line 128 of file ossimRS1SarModel.cpp.

References ossimLagrangeInterpolator::interpolate(), ossimIpt::line, ossimDpt::line, ossimIpt::samp, ossimDpt::samp, SCN, SCW, theArpPosInterp, theArpVelInterp, theCosOrientation, theFirstLineTime, ossimSensorModel::theImageSize, theImagingMode, theLatGrid, theLineScale, theLocalOrpInterp, theLonGrid, theRefHeight, theSinOrientation, theSinSkew, theTimePerLine, ossimDpt::x, ossimEcefPoint::y(), and ossimDpt::y.

 {
    const char* MODULE = "ossimRS1SarModel::imagingRay()";
 
    // Apply offset, scale, skew, and rotation to image point:
    ossimDpt offset (theImageSize.line/2.0, theImageSize.samp/2.0);
    ossimDpt ip1 = image_point - offset;
    ossimDpt ip2;
    ip2.line = ip1.line*theCosOrientation*theLineScale +
       ip1.samp*(theCosOrientation*theSinSkew+theSinOrientation);
    ip2.samp = -ip1.line*theSinOrientation*theLineScale -
       ip1.samp*(theSinOrientation*theSinSkew-theCosOrientation);
 
    // Offset point back to origin of image:
    ip1 = ip2 + offset;
 
    // Given new line number, obtain the interpolated ORP:
    ossimEcefPoint localORP;
    if ((theImagingMode == SCN) || (theImagingMode == SCW))
    {
       ossimGpt gpt (theLatGrid(ip1.x, ip1.y), theLonGrid(ip1.x, ip1.y), theRefHeight);
       localORP = ossimEcefPoint(gpt);
    }
    else
    {
       NEWMAT::ColumnVector orpPos(3);
       theLocalOrpInterp->interpolate(ip1.line, orpPos);
       localORP = ossimEcefPoint(orpPos[0], orpPos[1], orpPos[2]);
    }
    
    // Establish imaging line time (zero-Doppler time):
    double arpTime = theFirstLineTime + ip1.line*theTimePerLine;
 
    NEWMAT::ColumnVector arpEph(3);
    theArpPosInterp->interpolate(arpTime, arpEph);
    ossimEcefPoint ecfArpPos (arpEph[0], arpEph[1], arpEph[2]);
 
    theArpVelInterp->interpolate(arpTime, arpEph);
    ossimEcefVector ecfArpVel (arpEph[0], arpEph[1], arpEph[2]);
 
    // Need to correct the velocity vector by the earth rotational velocity:
    ossimEcefVector earthVel 
       (-ecfArpPos.y()*EARTH_ANGULAR_VELOCITY, ecfArpPos.x()*EARTH_ANGULAR_VELOCITY, 0.0);
    ecfArpVel = ecfArpVel - earthVel;
 
    // Compute a vehicle space with Z near intrack, X normal to Z and range
    // vector to ORP (slant range normal), and Y in range direction.
    // Note the application of the position adjustable parameter to the space's
    // origin (after the range vector was established):
    ossimEcefVector rangeVector     = localORP - ecfArpPos;
    ossimEcefVector rangeNormal     = rangeVector.cross(ecfArpVel);
    // ossimEcefVector zDirection      = rangeNormal.cross(rangeVector);
    ossimLsrSpace localVehicleSpace (ossimEcefPoint(ecfArpPos + thePosCorrection), // origin
                                     rangeNormal,  // X-direction
                                     rangeVector, 0);   // Y-direction (Z not specified
    double local_orp_range = rangeVector.length();
 
    // Compute slant range distance to pixel in question:
    double slant_range = local_orp_range;
    switch (theImagingMode)
    {
    case SLC:
       if (theDirectionFlag == ASCENDING)
          slant_range += thePixelSpacing.samp*ip1.samp;
       else
          slant_range += (theImageSize.samp-1.0-ip1.samp) * thePixelSpacing.samp;
       break;
 
    case SGC:
    case SGF:
    case SGX:
    case ERS:
       {
          // Compute the slant range as a polynomial expansion given ground range
          double ground_range;
          if (theDirectionFlag == ASCENDING)
             ground_range = ip1.samp * thePixelSpacing.samp;
          else
             ground_range = (theImageSize.samp - 1.0 - ip1.samp) *
             thePixelSpacing.samp;
          //         double delta_g_i = 1.0;
          //         for (int 0=1; i<6; i++)
          double delta_g_i = ground_range;
          for (int i=1; i<6; i++)  // NOTE: not using offset (i=0) term
          {
             slant_range += theSrGrCoeff[i] * delta_g_i;
             delta_g_i *= ground_range;
          }
       }
       break;
 
    case SCN:
    case SCW:
       slant_range = local_orp_range;
       break;
 
    default:
       CLOG << "ERROR: Invalid imaging mode encountered." << endl;
       setErrorStatus();
    }
 }

◆ initAdjParms()

void ossimRS1SarModel::initAdjParms ( )

protected

Definition at line 523 of file ossimRS1SarModel.cpp.

 {
    // Adjustable model not yet implemented
 #if 0
 
    static const char* MODULE = "ossimRS1SarModel::initAdjParms()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    theInTrackOffset = (init_adj_parm[INTRACK_OFFSET]+adj_parm[INTRACK_OFFSET])
                       * adj_sigma[INTRACK_OFFSET];
    
    theCrTrackOffset = (init_adj_parm[CRTRACK_OFFSET]+adj_parm[CRTRACK_OFFSET])
                       * adj_sigma[CRTRACK_OFFSET];
    
    theRadialOffset = (init_adj_parm[RADIAL_OFFSET]+adj_parm[RADIAL_OFFSET])
                      * adj_sigma[RADIAL_OFFSET];
    
    theLineScale   = (init_adj_parm[LINE_SCALE]+adj_parm[LINE_SCALE])
                     * adj_sigma[LINE_SCALE];
    
    theSkew        = (init_adj_parm[SKEW]+adj_parm[SKEW])
                     * adj_sigma[SKEW];
    
    theOrientation = (init_adj_parm[ORIENTATION]+adj_parm[ORIENTATION]);
 
    // Initialize base-class initial adjustable parameter array:
    init_adj_parm[INTRACK_OFFSET]  = theInTrackOffset/adj_sigma[INTRACK_OFFSET];
    init_adj_parm[CRTRACK_OFFSET]  = theCrTrackOffset/adj_sigma[CRTRACK_OFFSET];
    init_adj_parm[RADIAL_OFFSET]   = theRadialOffset/adj_sigma[RADIAL_OFFSET];  
    init_adj_parm[LINE_SCALE]      = theLineScale/adj_sigma[LINE_SCALE];  
    init_adj_parm[SKEW]            = theSkew/adj_sigma[SKEW];  
    init_adj_parm[ORIENTATION]     = theOrientation/adj_sigma[ORIENTATION];  
    
    //  Initialize sensor adjustable parameter description strings:
    strncpy (adj_desc[INTRACK_OFFSET], "intrack_offset", MAX_DESC_CHARS);
    strncpy (adj_desc[CRTRACK_OFFSET], "crtrack_offset", MAX_DESC_CHARS);
    strncpy (adj_desc[RADIAL_OFFSET],  "radial_offset",  MAX_DESC_CHARS);
 //   strncpy (adj_desc[SAMP_SCALE],  "samp_scale",      MAX_DESC_CHARS);
    strncpy (adj_desc[LINE_SCALE],   "line_scale",       MAX_DESC_CHARS);
    strncpy (adj_desc[SKEW],         "image_skew",       MAX_DESC_CHARS);
    strncpy (adj_desc[ORIENTATION],  "image_orientation",MAX_DESC_CHARS);
 
    // Initialize the adj_parms (parameter corrections) to 0:
    for (int i=0; i<num_adj_parm; i++)
       adj_parm[i] = 0.0;
    
    if (traceDebug())  CLOG << "returning..." << endl;
 #endif
 }

◆ initFromCeos()

void ossimRS1SarModel::initFromCeos ( const ossimFilename & dataDir )

protected

Definition at line 340 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::initFromCeos()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    // Instantiate a CeosData object:
    theCeosData = new ossimCeosData(fname);
    
    if(theCeosData->errorStatus() != ossimErrorCodes::OSSIM_OK) return;
    const dataset_sum_rec* dsr = theCeosData->dataSetSumRec();
    const proc_parm_rec*   ppr = theCeosData->procParmRec();
    char buf[1024];
    
    // Set imaging mode:
    strncpy(buf, &(theCeosData->volDescRec()->logvol_id[11]), 3);
    buf[3] = '\0';
    setImagingMode(buf);
    
    // Determine whether ascending or descending path acquisition. This flag
    // dictates whether the first or last pixel of a line is of shorter range:
    if (theImagingMode == ERS)
    {
       theDirectionFlag = DESCENDING;
    }
    else
    {
       if (dsr->asc_des[0] == 'A')
          theDirectionFlag = ASCENDING;
       else if (dsr->asc_des[0] == 'D')
          theDirectionFlag = DESCENDING;
       else
       {
          CLOG << "ERROR: Direction Flag: " << dsr->asc_des << " not supported"
               << endl;
          theDirectionFlag = UNKNOWN_DIRECTION;
          setErrorStatus();
          if (traceDebug())  CLOG << "returning with error..." << endl;
          return;
       }
    }
 
    // Image (Product) ID:
    buf[8] = '\0';
    strncpy(buf, theCeosData->volDescRec()->product_id, 8);
    theImageID = buf;
    
    // pixel spacing in range:
    buf[16] = '\0';
    strncpy(buf, dsr->pix_spacing, 16);
    thePixelSpacing.samp = atof(buf);
    
    // pixel spacing in azimuth:
    strncpy(buf, dsr->line_spacing, 16);
    buf[16] = '\0';
    thePixelSpacing.line = atof(buf);
 
    // The ground reference point (ORP) latitude:
    ossimGpt grp;
    strncpy(buf, dsr->pro_lat, 16);
    buf[16] = '\0';
    grp.lat = atof(buf);
       
    // The ground reference point (ORP) longitude:
    strncpy(buf, dsr->pro_long, 16);
    buf[16] = '\0';
    grp.lon = atof(buf);
       
    // The ground reference point (ORP) elevation:
    strncpy(buf, dsr->terrain_h, 16);
    buf[16] = '\0';
    grp.hgt = atof(buf);
    theORP = ossimEcefPoint(grp);
    theRefHeight = grp.hgt;
    
    // Illumination elevation:
    strncpy(buf, dsr->incident_ang, 8);
    buf[8] = '\0';
    double incidence = atof(buf);
    theIllumElevation = 90.0 - incidence;
 
    // Illumination azimuth -- need to determine whether sensor is in normal
    // right-looking mode or "antarctic" (left-looking) mode:
    strncpy(buf, dsr->plat_head, 8);
    theIllumAzimuth = atof(buf);
    char sensor_orientation[] = "123456789";
    strncpy(sensor_orientation, ppr->sens_orient, 9);
    if (strcmp(sensor_orientation, "ANTARCTIC") == 0)
    {
       theIllumAzimuth -= 90.0;
       if (theIllumAzimuth < 0.0) theIllumAzimuth += 360.0;
    }
    else
    {
       theIllumAzimuth += 90.0;
       if (theIllumAzimuth >= 360.0) theIllumAzimuth -= 360.0;
    }
 
    // Determine the number of lines per image and pixels per line. NOTE: if
    // image is scan mode, the number of lines must be computed indirectly. 
    if ((theImagingMode == SCN) || (theImagingMode == SCW))
    {
       FILE* fptr = fopen(theCeosData->imageFile().chars(), "r");
       if (!fptr)
       {
          CLOG << " ERROR:\n\tCannot open CEOS image file: "
               << theCeosData->imageFile() << endl;
          setErrorStatus();
          return;
       }
       fseek(fptr, 0, SEEK_END);
       long byte_count = ftell(fptr);
       buf[6] = '\0';
       strncpy(buf, theCeosData->imopDescRec()->l_dataset, 6);
       int lengthOfRec = atoi(buf);
       theImageSize.line = (byte_count - sizeof(*theCeosData->imopDescRec()))
                            /lengthOfRec;
       fclose(fptr);
    }
    else
    {
       strncpy(buf, theCeosData->imopDescRec()->n_dataset, 6);
       buf[6] = '\0';
       theImageSize.line = atoi(buf);
    }
    strncpy(buf, theCeosData->imopDescRec()->nleft, 4);
    buf[4] = '\0';
    int left_pixels = atoi(buf);
    strncpy(buf, theCeosData->imopDescRec()->ngrp, 8);
    buf[8] = '\0';
    int data_pixels = atoi(buf);
    strncpy(buf, theCeosData->imopDescRec()->nright, 4);
    buf[4] = '\0';
    int right_pixels  = atoi(buf);
    theImageSize.samp = left_pixels + data_pixels + right_pixels;
 
    // If ground range product, need to read the ground-to-slant range
    // conversion coefficients from the proc_parm_rec. NOTE: only single-look
    // products handled properly. ScanSAR requires reading multiple coefficient
    // sets. See MDA Detailed Processing Parameter Record Description sec.3.84.
    if ((theImagingMode==SGC) ||
        (theImagingMode==SGF) ||
        (theImagingMode==SGX) ||
        (theImagingMode==ERS))
    {
       buf[16] = '\0';
       for (int i=0; i<6; i++)
       {
          strncpy(buf, ppr->srgr_coefset[0].srgr_coef[i], 16);
          theSrGrCoeff[i] = atof(buf);
       }
    }
    else
    {
       theSrGrCoeff[0] = 0.0;
       theSrGrCoeff[1] = 1.0;  // note 1.0 here (should never be accessed)
       theSrGrCoeff[2] = 0.0;
       theSrGrCoeff[3] = 0.0;
       theSrGrCoeff[4] = 0.0;
       theSrGrCoeff[5] = 0.0;
    }
    
    if (traceDebug())
    {
       CLOG << "DEBUG -- "
            << "\n\t theDirectionFlag = " << theDirectionFlag
            << "\n\t thePixelSpacing = " << thePixelSpacing
            << "\n\t theORP = " << theORP
            << "\n\t theRefHeight = " << theRefHeight
            << "\n\t theImageSize = " << theImageSize
            << "\n\t sensor_orientation = " << sensor_orientation
            << "\n\t theIllumElevation = " << theIllumElevation
            << "\n\t theIllumAzimuth = " << theIllumAzimuth
            << endl;
    }
 
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ interpolatedScanORP()

void ossimRS1SarModel::interpolatedScanORP	(	const ossimDpt &	orp,
		ossimEcefPoint &	orp_ecf
	)		const

protected

Definition at line 904 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::interpolatedScanORP(gDblPoint)";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    ossimGpt gpt;
    gpt.lat = theLatGrid(orp);
    gpt.lon = theLonGrid(orp);
 
    // Convert to ECF:
    orp_ecf = ossimEcefPoint(gpt);
    
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ lineSampleHeightToWorld()

void ossimRS1SarModel::lineSampleHeightToWorld	(	const ossimDpt &	imagePt,
		const double &	heightAboveEllipsoid,
		ossimGpt &	worldPt
	)		const

virtual

Establishes geographic 3D point given image line, sample and ellipsoid height.

Implements ossimSensorModel.

Definition at line 235 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::lineSampleHeightToWorld()";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    // static const int    MAX_NUM_ITERS = 50;
    // static const double MAX_ELEV_DIFF = 0.001;
 
    ossimEcefRay imaging_ray;
    imagingRay(image_point, imaging_ray);
    worldPt = imaging_ray.intersectAboveEarthEllipsoid(height_ellip);
          
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ loadState()

bool ossimRS1SarModel::loadState	(	const ossimKeywordlist &	kwl,
		const char *	prefix = `NULL`
	)

virtual

Fulfills ossimObject base-class pure virtuals.

Reads modeling info from KWL. Returns true if successful.

Reimplemented from ossimSensorModel.

Definition at line 107 of file ossimRS1SarModel.cpp.

References ossimErrorStatusInterface::setErrorStatus().

 {
    // NOT YET IMPLEMENTED
    setErrorStatus();
    return false;
 }

◆ saveState()

bool ossimRS1SarModel::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix = `NULL`
	)		const

virtual

Fulfills ossimObject base-class pure virtuals.

Saves modeling info to KWL. Returns true if successful.

Reimplemented from ossimSensorModel.

Definition at line 117 of file ossimRS1SarModel.cpp.

References ossimErrorStatusInterface::setErrorStatus().

 {
    // NOT YET IMPLEMENTED
    setErrorStatus();
    return false;
 }

◆ setImagingMode()

void ossimRS1SarModel::setImagingMode ( char * modeStr )

protected

Definition at line 256 of file ossimRS1SarModel.cpp.

 {
    static const char MODULE[] = "ossimRS1SarModel::setImagingMode(modeStr)";
    if (traceDebug())  CLOG << "entering..." << endl;
 
    bool illegal_mode = false;
    
    if (strcmp(modeStr, "SCN") == 0)
    {
       theImagingMode = SCN;
    }
    else if (strcmp(modeStr, "SCW") == 0)
    {
       theImagingMode = SCW;
    }
    else if (strcmp(modeStr, "SGC") == 0)
    {
       theImagingMode = SGC;
    }
    else if (strcmp(modeStr, "SGF") == 0)
    {
       theImagingMode = SGF;
    }
    else if (strcmp(modeStr, "SGX") == 0)
    {
       theImagingMode = SGX;
    }
    else if (strcmp(modeStr, "SLC") == 0)
    {
       theImagingMode = SLC;
    }
    else if (strcmp(modeStr, "SPG") == 0)
    {
       theImagingMode = SPG;
       illegal_mode = true;
    }
    else if (strcmp(modeStr, "SSG") == 0)
    {
       theImagingMode = SSG;
       illegal_mode = true;
    }
    else if (strcmp(modeStr, "RAW") == 0)
    {
       theImagingMode = RAW;
       illegal_mode = true;
    }
    else if (strcmp(modeStr, "ERS") == 0)
    {
       theImagingMode = ERS;
    }
    else 
    {
       theImagingMode = UNKNOWN_MODE;
       illegal_mode = true;
    }
 
    // Test for ERS product:
    if ((theImagingMode == UNKNOWN_MODE) && theCeosData.valid())
    {
       char buf[] = "xxx";
       strncpy(buf, &(theCeosData->textRec()->product_type[8]), 3);
       if (strcmp(buf, "ERS") == 0)
       {
          theImagingMode = ERS;
          illegal_mode = false;
       }
    }
       
    // Filter the imaging modes that are not handled:
    if (illegal_mode)
    {
       CLOG << "\n\t ERROR: The imaging mode <"
            << IMAGING_MODE_ID[(int) theImagingMode]
            << "> is currently not supported. \n" << endl;
       theImagingMode = UNKNOWN_MODE;
    }
    
    if (traceDebug())  CLOG << "returning..." << endl;
 }

◆ useForward()

virtual bool ossimRS1SarModel::useForward ( ) const

inlinevirtual

useForward() return true when it's better (more accurate / fast) to use forward (from ground to image) than inverse(from image to ground)

Implements ossimOptimizableProjection.

Definition at line 84 of file ossimRS1SarModel.h.

84 { return false; }

Member Data Documentation

◆ theArpPosInterp

ossimRefPtr<ossimLagrangeInterpolator> ossimRS1SarModel::theArpPosInterp

protected

Definition at line 102 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theArpVelInterp

ossimRefPtr<ossimLagrangeInterpolator> ossimRS1SarModel::theArpVelInterp

protected

Definition at line 103 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theCeosData

ossimRefPtr<ossimCeosData> ossimRS1SarModel::theCeosData

protected

Definition at line 101 of file ossimRS1SarModel.h.

◆ theCosOrientation

double ossimRS1SarModel::theCosOrientation

protected

Definition at line 113 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theCrTrackOffset

double ossimRS1SarModel::theCrTrackOffset

protected

Definition at line 133 of file ossimRS1SarModel.h.

◆ theDirectionFlag

DirectionFlag ossimRS1SarModel::theDirectionFlag

protected

Definition at line 120 of file ossimRS1SarModel.h.

◆ theEphFirstSampTime

double ossimRS1SarModel::theEphFirstSampTime

protected

Definition at line 111 of file ossimRS1SarModel.h.

◆ theFirstLineDay

int ossimRS1SarModel::theFirstLineDay

protected

Definition at line 117 of file ossimRS1SarModel.h.

◆ theFirstLineTime

double ossimRS1SarModel::theFirstLineTime

protected

Definition at line 118 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theGHA

double ossimRS1SarModel::theGHA

protected

Definition at line 110 of file ossimRS1SarModel.h.

◆ theIllumAzimuth

double ossimRS1SarModel::theIllumAzimuth

protected

Definition at line 106 of file ossimRS1SarModel.h.

◆ theIllumElevation

double ossimRS1SarModel::theIllumElevation

protected

Definition at line 107 of file ossimRS1SarModel.h.

◆ theImagingMode

ImagingMode ossimRS1SarModel::theImagingMode

protected

Definition at line 105 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theInTrackOffset

double ossimRS1SarModel::theInTrackOffset

protected

Definition at line 132 of file ossimRS1SarModel.h.

◆ theLatGrid

ossimDblGrid ossimRS1SarModel::theLatGrid

protected

Definition at line 126 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theLineScale

double ossimRS1SarModel::theLineScale

protected

Definition at line 135 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theLocalOrpInterp

ossimRefPtr<ossimLagrangeInterpolator> ossimRS1SarModel::theLocalOrpInterp

protected

Definition at line 104 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theLonGrid

ossimDblGrid ossimRS1SarModel::theLonGrid

protected

Definition at line 127 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theOrientation

double ossimRS1SarModel::theOrientation

protected

Definition at line 137 of file ossimRS1SarModel.h.

◆ theORP

ossimEcefPoint ossimRS1SarModel::theORP

protected

Definition at line 108 of file ossimRS1SarModel.h.

◆ thePixelSpacing

ossimDpt ossimRS1SarModel::thePixelSpacing

protected

Definition at line 122 of file ossimRS1SarModel.h.

◆ thePosCorrection

ossimEcefVector ossimRS1SarModel::thePosCorrection

protected

Definition at line 116 of file ossimRS1SarModel.h.

◆ theRadialOffset

double ossimRS1SarModel::theRadialOffset

protected

Definition at line 134 of file ossimRS1SarModel.h.

◆ theRefHeight

double ossimRS1SarModel::theRefHeight

protected

Definition at line 109 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theSinOrientation

double ossimRS1SarModel::theSinOrientation

protected

Definition at line 112 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theSinSkew

double ossimRS1SarModel::theSinSkew

protected

Definition at line 114 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theSkew

double ossimRS1SarModel::theSkew

protected

Definition at line 136 of file ossimRS1SarModel.h.

◆ theSrGrCoeff

double ossimRS1SarModel::theSrGrCoeff[6]

protected

Definition at line 121 of file ossimRS1SarModel.h.

◆ theTimePerLine

double ossimRS1SarModel::theTimePerLine

protected

Definition at line 119 of file ossimRS1SarModel.h.

Referenced by imagingRay().

◆ theVehicleSpace

ossimLsrSpace ossimRS1SarModel::theVehicleSpace

protected

Definition at line 115 of file ossimRS1SarModel.h.

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ ADJUSTABLE_PARAM_INDEXES

◆ DirectionFlag

◆ ImagingMode

Constructor & Destructor Documentation

◆ ossimRS1SarModel() [1/2]

◆ ossimRS1SarModel() [2/2]

◆ ~ossimRS1SarModel()

Member Function Documentation

◆ deallocateMemory()

◆ dup()

◆ eciToEcfXform()

◆ establishEphemeris()

◆ establishOrpGrid()

◆ establishOrpInterp()

◆ establishVehicleSpace()

◆ imagingRay()

◆ initAdjParms()

◆ initFromCeos()

◆ interpolatedScanORP()

◆ lineSampleHeightToWorld()

◆ loadState()

◆ saveState()

◆ setImagingMode()

◆ useForward()

Member Data Documentation

◆ theArpPosInterp

◆ theArpVelInterp

◆ theCeosData

◆ theCosOrientation

◆ theCrTrackOffset

◆ theDirectionFlag

◆ theEphFirstSampTime

◆ theFirstLineDay

◆ theFirstLineTime

◆ theGHA

◆ theIllumAzimuth

◆ theIllumElevation

◆ theImagingMode

◆ theInTrackOffset

◆ theLatGrid

◆ theLineScale

◆ theLocalOrpInterp

◆ theLonGrid

◆ theOrientation

◆ theORP

◆ thePixelSpacing

◆ thePosCorrection

◆ theRadialOffset

◆ theRefHeight

◆ theSinOrientation

◆ theSinSkew

◆ theSkew

◆ theSrGrCoeff

◆ theTimePerLine

◆ theVehicleSpace