#include <ossimSpaceObliqueMercatorProjection.h>

Inheritance diagram for ossimSpaceObliqueMercatorProjection:

Public Types
enum	ossimSatelliteType { SOM_TYPE_LANDSAT_1 = 0, SOM_TYPE_LANDSAT_2 = 1, SOM_TYPE_LANDSAT_3 = 2, SOM_TYPE_LANDSAT_4 = 3, SOM_TYPE_LANDSAT_5 = 4, SOM_TYPE_LANDSAT_7 = 5 }

Public Member Functions
	ossimSpaceObliqueMercatorProjection (ossimSatelliteType type=SOM_TYPE_LANDSAT_7, double pathNumber=34, const ossimEllipsoid &ellipsoid=ossimEllipsoid())

virtual ossimObject *	dup () const

virtual ossimDpt	forward (const ossimGpt &worldPoint) const
	All map projections will convert the world coordinate to an easting northing (Meters). More...

virtual ossimGpt	inverse (const ossimDpt &projectedPoint) const
	Will take a point in meters and convert it to ground. More...

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)

virtual void	setParameters (ossimSatelliteType type, double path)

virtual void	update ()

virtual bool	operator== (const ossimProjection &projection) const
	Returns TRUE if principal parameters are within epsilon tolerance. More...

Public Member Functions inherited from ossimMapProjection
	ossimMapProjection (const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())

	ossimMapProjection (const ossimMapProjection &src)

virtual ossimGpt	origin () const

virtual ossimDpt	worldToLineSample (const ossimGpt &worldPoint) const

virtual void	worldToLineSample (const ossimGpt &worldPoint, ossimDpt &lineSample) const

virtual ossimGpt	lineSampleToWorld (const ossimDpt &projectedPoint) const

virtual void	lineSampleToWorld (const ossimDpt &projectedPoint, ossimGpt &gpt) const

virtual void	lineSampleHeightToWorld (const ossimDpt &lineSampPt, const double &heightAboveEllipsoid, ossimGpt &worldPt) const
	This is the pure virtual that projects the image point to the given elevation above ellipsoid, thereby bypassing reference to a DEM. More...

virtual void	lineSampleToEastingNorthing (const ossimDpt &liineSample, ossimDpt &eastingNorthing) const

virtual void	eastingNorthingToLineSample (const ossimDpt &eastingNorthing, ossimDpt &lineSample) const

virtual void	eastingNorthingToWorld (const ossimDpt &eastingNorthing, ossimGpt &worldPt) const

virtual double	getFalseEasting () const

virtual double	getFalseNorthing () const

virtual double	getStandardParallel1 () const
	Derived classes should implement as needed. More...

virtual double	getStandardParallel2 () const
	Derived classes should implement as needed. More...

virtual void	setPcsCode (ossim_uint32 pcsCode)

virtual ossim_uint32	getPcsCode () const
	Returns the EPSG PCS code or 32767 if the projection is a custom (non-EPSG) projection. More...

virtual ossimString	getProjectionName () const
	Returns the projection name. More...

virtual double	getA () const
	ACCESS METHODS: More...

virtual double	getB () const

virtual double	getF () const

virtual ossimDpt	getMetersPerPixel () const

virtual const ossimDpt &	getDecimalDegreesPerPixel () const
	Returns decimal degrees per pixel as an ossimDpt with "x" representing longitude and "y" representing latitude. More...

virtual const ossimDpt &	getUlEastingNorthing () const

virtual const ossimGpt &	getUlGpt () const

virtual const ossimDatum *	getDatum () const

const ossimEllipsoid &	getEllipsoid () const

const ossimGpt &	getOrigin () const

virtual bool	isGeographic () const

virtual void	applyScale (const ossimDpt &scale, bool recenterTiePoint)
	Applies scale to theDeltaLonPerPixel, theDeltaLatPerPixel and theMetersPerPixel data members (eg: theDeltaLonPerPixel *= scale.x). More...

virtual void	setEllipsoid (const ossimEllipsoid &ellipsoid)
	SET METHODS: More...

virtual void	setAB (double a, double b)

virtual void	setDatum (const ossimDatum *datum)
	Sets theDatum to datum. More...

virtual void	setOrigin (const ossimGpt &origin)
	Sets theOrigin to origin. More...

virtual void	setMetersPerPixel (const ossimDpt &gsd)

virtual void	setDecimalDegreesPerPixel (const ossimDpt &gsd)

virtual void	setUlTiePoints (const ossimGpt &gpt)

virtual void	setUlTiePoints (const ossimDpt &eastingNorthing)

virtual void	setUlEastingNorthing (const ossimDpt &ulEastingNorthing)

virtual void	setUlGpt (const ossimGpt &ulGpt)

virtual void	assign (const ossimProjection &aProjection)

virtual std::ostream &	print (std::ostream &out) const
	Prints data members to stream. More...

virtual void	computeDegreesPerPixel ()
	Computes the approximate resolution in degrees/pixel. More...

virtual void	computeMetersPerPixel ()
	This will go from the ground point and give you an approximate meters per pixel. More...

void	setMatrix (double rotation, const ossimDpt &scale, const ossimDpt &translation)

void	setMatrixScale (const ossimDpt &scale)

void	setMatrixRotation (double rotation)

void	setMatrixTranslation (const ossimDpt &translation)

void	snapTiePointTo (ossim_float64 multiple, ossimUnitType unitType)
	Utility method to snap the tie point to some multiple. More...

void	snapTiePointToOrigin ()

void	setElevationLookupFlag (bool flag)

bool	getElevationLookupFlag () const

ossimUnitType	getModelTransformUnitType () const

void	setModelTransformUnitType (ossimUnitType unit)

bool	hasModelTransform () const

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

void	setProjectionUnits (ossimUnitType units)

ossimUnitType	getProjectionUnits () const
	OSSIM considers all map projection coordinates (including false eastings and northings) to be in meters. More...

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

Public Member Functions inherited from ossimProjection
	ossimProjection ()

virtual	~ossimProjection ()

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	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

Protected Member Functions
void	seraz0 (double lam, double mult)

Protected Member Functions inherited from ossimMapProjection
virtual	~ossimMapProjection ()

void	updateFromTransform ()

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
double	a2

double	a4

double	b

double	c1

double	c3

double	q

double	t

double	u

double	w

double	p22

double	sa

double	ca

double	xj

double	rlm

double	rlm2

double	lam0

double	es

double	e

double	one_es

double	rone_es

double	a

double	thePath

ossimSatelliteType	theSatelliteType

Protected Attributes inherited from ossimMapProjection
ossimEllipsoid	theEllipsoid
	This method verifies that the projection parameters match the current pcs code. More...

ossimGpt	theOrigin

const ossimDatum *	theDatum
	This is only set if we want to have built in datum shifting. More...

ossimDpt	theMetersPerPixel
	Holds the number of meters per pixel. More...

ossimDpt	theDegreesPerPixel
	Hold the decimal degrees per pixel. More...

ossimGpt	theUlGpt
	Hold tie point in decimal degrees. More...

ossimDpt	theUlEastingNorthing
	Hold tie point as easting northing. More...

ossimDpt	theFalseEastingNorthing
	Hold the false easting northing. More...

ossim_uint32	thePcsCode
	Projection Coordinate System(PCS) code. More...

bool	theElevationLookupFlag

ossimMatrix4x4	theModelTransform

ossimMatrix4x4	theInverseModelTransform

ossimUnitType	theModelTransformUnitType

ossimUnitType	theProjectionUnits
	Linear units of the projection as indicated in the projection's specification: More...

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Detailed Description

This class implments the formulas for the ellipsoid and circular orbit. All formulas were directly taken from proj 4. Proj 4 implementation is from the following manual reference:

Map Projections - A Working Manual by John Snyder Space Oblique Mercator p. 221.

Definition at line 27 of file ossimSpaceObliqueMercatorProjection.h.

Member Enumeration Documentation

◆ ossimSatelliteType

enum ossimSpaceObliqueMercatorProjection::ossimSatelliteType

Enumerator
SOM_TYPE_LANDSAT_1
SOM_TYPE_LANDSAT_2
SOM_TYPE_LANDSAT_3
SOM_TYPE_LANDSAT_4
SOM_TYPE_LANDSAT_5
SOM_TYPE_LANDSAT_7

Definition at line 30 of file ossimSpaceObliqueMercatorProjection.h.

    {
       SOM_TYPE_LANDSAT_1 = 0,
       SOM_TYPE_LANDSAT_2 = 1,
       SOM_TYPE_LANDSAT_3 = 2,
       SOM_TYPE_LANDSAT_4 = 3, 
       SOM_TYPE_LANDSAT_5 = 4,
       SOM_TYPE_LANDSAT_7 = 5
    };

Constructor & Destructor Documentation

◆ ossimSpaceObliqueMercatorProjection()

ossimSpaceObliqueMercatorProjection::ossimSpaceObliqueMercatorProjection	(	ossimSatelliteType	type = `SOM_TYPE_LANDSAT_7`,
		double	pathNumber = `34`,
		const ossimEllipsoid &	ellipsoid = `ossimEllipsoid()`
	)

Definition at line 81 of file ossimSpaceObliqueMercatorProjection.cpp.

       :
       ossimMapProjection(ellipsoid)
 {
    setParameters(type, pathNumber);
 }

Member Function Documentation

◆ dup()

virtual ossimObject* ossimSpaceObliqueMercatorProjection::dup ( ) const

inlinevirtual

Implements ossimProjection.

Definition at line 44 of file ossimSpaceObliqueMercatorProjection.h.

44 {return new ossimSpaceObliqueMercatorProjection(*this);}

ossimSpaceObliqueMercatorProjection::ossimSpaceObliqueMercatorProjection

ossimSpaceObliqueMercatorProjection(ossimSatelliteType type=SOM_TYPE_LANDSAT_7, double pathNumber=34, const ossimEllipsoid &ellipsoid=ossimEllipsoid())

Definition: ossimSpaceObliqueMercatorProjection.cpp:81

◆ forward()

ossimDpt ossimSpaceObliqueMercatorProjection::forward ( const ossimGpt & worldPoint ) const

virtual

All map projections will convert the world coordinate to an easting northing (Meters).

Implements ossimMapProjection.

Definition at line 192 of file ossimSpaceObliqueMercatorProjection.cpp.

References a, a2, a4, aasin(), b, c1, c3, ca, ossimGpt::changeDatum(), ossimDatum::code(), ossimGpt::datum(), es, FORTPI, HALFPI, lam0, ossimGpt::latr(), ossimGpt::lonr(), ossimDpt::makeNan(), one_es, p22, PI_HALFPI, q, rlm, rlm2, sa, t, ossimMapProjection::theDatum, TOL, TWOPI_HALFPI, w, ossimDpt::x, xj, and ossimDpt::y.

 {
    ossimGpt gpt = worldPoint;
    
    if (theDatum)
    {
       if (theDatum->code() != worldPoint.datum()->code())
       {
          gpt.changeDatum(theDatum); // Shift to our datum.
       }
    }
 
    ossimDpt xy;
    
    int l, nn;
    double lamt, xlam, sdsq, c, d, s, lamdp, phidp, lampp, tanph,
       lamtp, cl, sd, sp, fac, sav, tanphi;
    double phi = gpt.latr();
    double lam = gpt.lonr() - lam0;
    
    if (phi > HALFPI)
       phi = HALFPI;
    else if (phi < -HALFPI)
       phi = -HALFPI;
    lampp = phi >= 0. ? HALFPI : PI_HALFPI;
    tanphi = tan(phi);
    for (nn = 0;;)
    {
       sav = lampp;
       lamtp = lam + p22 * lampp;
       cl = cos(lamtp);
       if (fabs(cl) < TOL)
          lamtp -= TOL;
       fac = lampp - sin(lampp) * (cl < 0. ? -HALFPI : HALFPI);
       for (l = 50; l; --l) {
          lamt = lam + p22 * sav;
          if (fabs(c = cos(lamt)) < TOL)
             lamt -= TOL;
          xlam = (one_es * tanphi * sa + sin(lamt) * ca) / c;
          lamdp = atan(xlam) + fac;
          if (fabs(fabs(sav) - fabs(lamdp)) < TOL)
             break;
          sav = lamdp;
       }
       if (!l || ++nn >= 3 || (lamdp > rlm && lamdp < rlm2))
          break;
       if (lamdp <= rlm)
          lampp = TWOPI_HALFPI;
       else if (lamdp >= rlm2)
          lampp = HALFPI;
    }
    if (l) {
       sp = sin(phi);
       phidp = aasin((one_es * ca * sp - sa * cos(phi) * 
                      sin(lamt)) / sqrt(1. - es * sp * sp));
       tanph = log(tan(FORTPI + .5 * phidp));
       sd = sin(lamdp);
       sdsq = sd * sd;
       s = p22 * sa * cos(lamdp) * sqrt((1. + t * sdsq)
                                        / ((1. + w * sdsq) * (1. + q * sdsq)));
       d = sqrt(xj * xj + s * s);
       xy.x = b * lamdp + a2 * sin(2. * lamdp) + a4 *
              sin(lamdp * 4.) - tanph * s / d;
       xy.y = c1 * sd + c3 * sin(lamdp * 3.) + tanph * xj / d;
    }
    else
    {
       xy.makeNan();
       return xy;
 //       xy.x = xy.y = 1.0/DBL_EPSILON; // st to large number
    }
 
    xy.x*=a;
    xy.y*=a;
    return xy;
 }

◆ inverse()

ossimGpt ossimSpaceObliqueMercatorProjection::inverse ( const ossimDpt & projectedPoint ) const

virtual

Will take a point in meters and convert it to ground.

Implements ossimMapProjection.

Definition at line 273 of file ossimSpaceObliqueMercatorProjection.cpp.

References a, a2, a4, aasin(), b, c1, c3, ca, DEG_PER_RAD, es, FORTPI, HALFPI, lam0, one_es, p22, q, rone_es, sa, t, ossimMapProjection::theDatum, TOL, u, w, ossimDpt::x, xj, and ossimDpt::y.

 {
    int nn;
    double lamt, sdsq, s, lamdp, phidp, sppsq, dd, sd, sl, fac, scl, sav, spp;
    double lam, phi;
    
    ossimDpt xy = projectedPoint;
    xy.x/=a;
    xy.y/=a;
    
    lamdp = xy.x / b;
    nn = 50;
    do {
       sav = lamdp;
       sd = sin(lamdp);
       sdsq = sd * sd;
       s = p22 * sa * cos(lamdp) * sqrt((1. + t * sdsq)
                                        / ((1. + w * sdsq) * (1. + q * sdsq)));
       lamdp = xy.x + xy.y * s / xj - a2 * sin(
          2. * lamdp) - a4 * sin(lamdp * 4.) - s / xj * (
             c1 * sin(lamdp) + c3 * sin(lamdp * 3.));
       lamdp /= b;
    } while (fabs(lamdp - sav) >= TOL && --nn);
    sl = sin(lamdp);
    fac = exp(sqrt(1. + s * s / xj / xj) * (xy.y - 
                                            c1 * sl - c3 * sin(lamdp * 3.)));
    phidp = 2. * (atan(fac) - FORTPI);
    dd = sl * sl;
    if (fabs(cos(lamdp)) < TOL)
       lamdp -= TOL;
    spp = sin(phidp);
    sppsq = spp * spp;
    lamt = atan(((1. - sppsq * rone_es) * tan(lamdp) * ca -
                 spp * sa * sqrt((1. + q * dd) * (1. - sppsq) - sppsq * u) /
                 cos(lamdp)) / (1. - sppsq * (1. + u)));
    sl = lamt >= 0. ? 1. : -1.;
    scl = cos(lamdp) >= 0. ? 1. : -1;
    lamt -= HALFPI * (1. - scl) * sl;
    
    lam = lamt - p22 * lamdp;
    if (fabs(sa) < TOL)
       phi = aasin(spp / sqrt(one_es * one_es + es * sppsq));
    else
       phi = atan((tan(lamdp) * cos(lamt) - ca * sin(lamt)) /
                     (one_es * sa));
    return ossimGpt(phi*DEG_PER_RAD,
                    (lam+lam0)*DEG_PER_RAD, 0.0, theDatum);
    
 }                        

◆ loadState()

bool ossimSpaceObliqueMercatorProjection::loadState	(	const ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)

virtual

Method to the load (recreate) the state of an object from a keyword list. Return true if ok or false on error.

Reimplemented from ossimMapProjection.

Definition at line 358 of file ossimSpaceObliqueMercatorProjection.cpp.

References ossimKeywordlist::find().

 {
    bool result = true;
 
    const char* path = kwl.find(prefix, PATH_KW);
    const char* type = kwl.find(prefix,  SATELLITE_TYPE_KW);
 
    if(path)
    {
       thePath = ossimString(path).toDouble();
    }
    else
    {
       result = false;
    }
    if(type)
    {
       theSatelliteType = (ossimSatelliteType)ossimString(type).toInt();
    }
    else
    {
       result = false;
    }
 
    setParameters(theSatelliteType, thePath);
 
    if(result)
    {
       result = ossimMapProjection::loadState(kwl, prefix);
    }
    else
    {
       ossimMapProjection::loadState(kwl, prefix);
    }
 
    return result;
 }

◆ operator==()

bool ossimSpaceObliqueMercatorProjection::operator== ( const ossimProjection & projection ) const

virtual

Returns TRUE if principal parameters are within epsilon tolerance.

Reimplemented from ossimMapProjection.

Definition at line 400 of file ossimSpaceObliqueMercatorProjection.cpp.

References ossim::almostEqual(), thePath, and theSatelliteType.

 {
    if (!ossimMapProjection::operator==(proj))
       return false;
 
    const ossimSpaceObliqueMercatorProjection* p =
       dynamic_cast<const ossimSpaceObliqueMercatorProjection*>(&proj);
    if (!p) return false;
 
    if (theSatelliteType != p->theSatelliteType) return false;
    if (!ossim::almostEqual(thePath,p->thePath)) return false;
 
    return true;
 }

◆ saveState()

bool ossimSpaceObliqueMercatorProjection::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)		const

virtual

Method to save the state of an object to a keyword list. Return true if ok or false on error.

Reimplemented from ossimMapProjection.

Definition at line 340 of file ossimSpaceObliqueMercatorProjection.cpp.

References ossimKeywordlist::add().

 {
    kwl.add(prefix,
            PATH_KW,
            thePath,
            true);
    
    kwl.add(prefix,
            SATELLITE_TYPE_KW,
            (int)theSatelliteType,
            true);
    return ossimMapProjection::saveState(kwl, prefix);
 }

◆ seraz0()

void ossimSpaceObliqueMercatorProjection::seraz0	(	double	lam,
		double	mult
	)

protected

Definition at line 165 of file ossimSpaceObliqueMercatorProjection.cpp.

References a2, a4, b, c1, c3, ca, p22, q, RAD_PER_DEG, sa, sq, t, w, and xj.

Referenced by setParameters().

 {
     double sdsq, h, s, fc, sd, sq, d__1;
 
     lam *= RAD_PER_DEG;
     sd = sin(lam);
     sdsq = sd * sd;
     s = p22 * sa * cos(lam) * sqrt((1. + t * sdsq) / ((
         1. + w * sdsq) * (1. + q * sdsq)));
     d__1 = 1. + q * sdsq;
     h = sqrt((1. + q * sdsq) / (1. + w * sdsq)) * ((1. + 
         w * sdsq) / (d__1 * d__1) - p22 * ca);
     sq = sqrt(xj * xj + s * s);
     b += fc = mult * (h * xj - s * s) / sq;
     a2 += fc * cos(lam + lam);
     a4 += fc * cos(lam * 4.);
     fc = mult * s * (h + xj) / sq;
     c1 += fc * cos(lam);
     c3 += fc * cos(lam * 3.);
    
 }

◆ setParameters()

void ossimSpaceObliqueMercatorProjection::setParameters	(	ossimSatelliteType	type,
		double	path
	)

virtual

Definition at line 93 of file ossimSpaceObliqueMercatorProjection.cpp.

References ossimEllipsoid::a(), a, a2, a4, b, c1, c3, ca, e, ossimEllipsoid::eccentricity(), es, lam0, one_es, p22, PI, q, RAD_PER_DEG, rlm, rlm2, rone_es, sa, seraz0(), SOM_TYPE_LANDSAT_1, SOM_TYPE_LANDSAT_2, SOM_TYPE_LANDSAT_3, SOM_TYPE_LANDSAT_4, SOM_TYPE_LANDSAT_5, SOM_TYPE_LANDSAT_7, t, ossimMapProjection::theEllipsoid, thePath, theSatelliteType, TWOPI, u, update(), w, and xj.

 {
    thePath = path;
    theSatelliteType = type;
    double lam = 0.0;
    double alf = 0.0;
    double esc = 0.0;
    double ess = 0.0;
    
    e  = theEllipsoid.eccentricity();
    es = e*e;
    one_es = 1- es;
    rone_es = 1.0/one_es;
    a   = theEllipsoid.a();
    
    switch(theSatelliteType)
    {
    case SOM_TYPE_LANDSAT_1:
    case SOM_TYPE_LANDSAT_2:
    case SOM_TYPE_LANDSAT_3:
    {
       lam0 = RAD_PER_DEG*(128.87 - (360.0/251.0)*path);
       p22 = 103.2669323;
       alf = RAD_PER_DEG * 99.092;
       break;
    }
    case SOM_TYPE_LANDSAT_4:
    case SOM_TYPE_LANDSAT_5:
    case SOM_TYPE_LANDSAT_7:
    {
       
       lam0 = RAD_PER_DEG*(129.305582487 - (360.0/233.0)*path);
       p22 = 98.8841202;
       alf = RAD_PER_DEG * 98.2;
       break;
    }
    }
    p22 /= 1440.;
    sa = sin(alf);
    ca = cos(alf);
    if (fabs(ca) < 1e-9)
       ca = 1e-9;
    esc = es * ca * ca;
    ess = es * sa * sa;
    w = (1. - esc) * rone_es;
    w = w * w - 1.;
    q = ess * rone_es;
    t = ess * (2. - es) * rone_es * rone_es;
    u = esc * rone_es;
    xj = one_es * one_es * one_es;
    rlm = PI * (1. / 248. + .5161290322580645);
    rlm2 = rlm + TWOPI;
    a2 = a4 = b = c1 = c3 = 0.;
    seraz0(0.0, 1.0);
    for (lam = 9.; lam <= 81.0001; lam += 18.)
       seraz0(lam, 4.0);
    for (lam = 18; lam <= 72.0001; lam += 18.)
       seraz0(lam, 2.0);
    seraz0(90.0, 1.0);
    a2 /= 30.0;
    a4 /= 60.0;
    b /= 30.0;
    c1 /= 15.0;
    c3 /= 45.0;
 
    update();
 }

◆ update()

void ossimSpaceObliqueMercatorProjection::update ( )

virtual

Reimplemented from ossimMapProjection.

Definition at line 327 of file ossimSpaceObliqueMercatorProjection.cpp.

Referenced by setParameters().

 {
 //   theMetersPerPixel.x    = 1.0;
 //   theMetersPerPixel.y    = 1.0;
    
 //   theUlEastingNorthing.x = 0.0;
 //   theUlEastingNorthing.y = 0.0;
 }