/usr/include/ossim/projection/ossimSinusoidalProjection.h is in libossim-dev 2.2.2-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 | //*******************************************************************
//
// License: See top level LICENSE.txt file.
//
// Author: Garrett Potts (gpotts@imagelinks.com)
//
// Description:
//
// Calls Geotrans Sinusoidal projection code.
//*******************************************************************
// $Id: ossimSinusoidalProjection.h 17815 2010-08-03 13:23:14Z dburken $
#ifndef ossimSinusoidalProjection_HEADER
#define ossimSinusoidalProjection_HEADER
#include <ossim/projection/ossimMapProjection.h>
class OSSIMDLLEXPORT ossimSinusoidalProjection : public ossimMapProjection
{
public:
ossimSinusoidalProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
const ossimGpt& origin = ossimGpt());
ossimSinusoidalProjection(const ossimEllipsoid& ellipsoid,
const ossimGpt& origin,
const double falseEasting,
const double falseNorthing);
~ossimSinusoidalProjection(){}
virtual ossimObject *dup()const{return new ossimSinusoidalProjection(*this);}
virtual ossimGpt inverse(const ossimDpt &eastingNorthing)const;
virtual ossimDpt forward(const ossimGpt &latLon)const;
virtual void update();
/*!
* SetFalseEasting. The value is in meters.
* Update is then called so we can pre-compute paramters
*/
void setFalseEasting(double falseEasting);
/*!
* SetFalseNorthing. The value is in meters.
* Update is then called so we can pre-compute paramters
*/
void setFalseNorthing(double falseNorthing);
/*!
* Sets both false easting and northing values. The values are
* expected to be in meters.
* Update is then called so we can pre-compute paramters
*/
void setFalseEastingNorthing(double falseEasting, double falseNorthing);
void setDefaults();
void setCentralMeridian(double centralMeridian);
double getFalseEasting()const{return Sinu_False_Easting;}
double getFalseNorthing()const{return Sinu_False_Northing;}
/*!
* Method to save the state of an object to a keyword list.
* Return true if ok or false on error.
*/
virtual bool saveState(ossimKeywordlist& kwl,
const char* prefix=0)const;
/*!
* Method to the load (recreate) the state of an object from a keyword
* list. Return true if ok or false on error.
*/
virtual bool loadState(const ossimKeywordlist& kwl,
const char* prefix=0);
//! Returns TRUE if principal parameters are within epsilon tolerance.
virtual bool operator==(const ossimProjection& projection) const;
private:
//_______________________GEOTRANS________________________
mutable double Sinu_a; /* Semi-major axis of ellipsoid in meters */
mutable double Sinu_f; /* Flattening of ellipsoid */
mutable double es2; /* Eccentricity (0.08181919084262188000) squared */
mutable double es4; /* es2 * es2 */
mutable double es6; /* es4 * es2 */
mutable double c0; /* 1 - es2 / 4.0 - 3.0 * es4 / 64.0 - 5.0 * es6 / 256.0 */
mutable double c1; /* 3.0 * es2 / 8.0 + 3.0 * es4 / 32.0 + 45.0 * es6 / 1024.0 */
mutable double c2; /* 15.0 * es4 / 256.0 + 45.0 * es6 / 1024.0 */
mutable double c3; /* 35.0 * es6 / 3072.0 */
mutable double a0; /* 3.0 * e1 / 2.0 - 27.0 * e3 / 32.0 */
mutable double a1; /* 21.0 * e2 / 16.0 - 55.0 * e4 / 32.0 */
mutable double a2; /* 151.0 * e3 / 96.0 */
mutable double a3; /* 1097.0 * e4 / 512.0 */
/* Sinusoid projection Parameters */
mutable double Sinu_Origin_Long; /* Longitude of origin in radians */
mutable double Sinu_False_Northing; /* False northing in meters */
mutable double Sinu_False_Easting; /* False easting in meters */
/* Maximum variance for easting and northing values for WGS 84.
*/
mutable double Sinu_Max_Easting;
mutable double Sinu_Min_Easting;
mutable double Sinu_Delta_Northing;
/*!
* The function Set_Sinusoidal_Parameters receives the ellipsoid parameters and
* Sinusoidal projcetion parameters as inputs, and sets the corresponding state
* variables. If any errors occur, the error code(s) are returned by the function,
* otherwise SINU_NO_ERROR is returned.
*
* a : Semi-major axis of ellipsoid, in meters (input)
* f : Flattening of ellipsoid (input)
* Central_Meridian : Longitude in radians at the center of (input)
* the projection
* False_Easting : A coordinate value in meters assigned to the
* central meridian of the projection. (input)
* False_Northing : A coordinate value in meters assigned to the
* origin latitude of the projection (input)
*/
long Set_Sinusoidal_Parameters(double a,
double f,
double Central_Meridian,
double False_Easting,
double False_Northing);
/*!
* The function Get_Sinusoidal_Parameters returns the current ellipsoid
* parameters and Sinusoidal projection parameters.
*
* a : Semi-major axis of ellipsoid, in meters (output)
* f : Flattening of ellipsoid (output)
* Central_Meridian : Longitude in radians at the center of (output)
* the projection
* False_Easting : A coordinate value in meters assigned to the
* central meridian of the projection. (output)
* False_Northing : A coordinate value in meters assigned to the
* origin latitude of the projection (output)
*/
void Get_Sinusoidal_Parameters(double *a,
double *f,
double *Central_Meridian,
double *False_Easting,
double *False_Northing)const;
/*!
* The function Convert_Geodetic_To_Sinusoidal converts geodetic (latitude and
* longitude) coordinates to Sinusoidal projection easting, and northing
* coordinates, according to the current ellipsoid and Sinusoidal projection
* parameters. If any errors occur, the error code(s) are returned by the
* function, otherwise SINU_NO_ERROR is returned.
*
* Latitude : Latitude (phi) in radians (input)
* Longitude : Longitude (lambda) in radians (input)
* Easting : Easting (X) in meters (output)
* Northing : Northing (Y) in meters (output)
*/
long Convert_Geodetic_To_Sinusoidal (double Latitude,
double Longitude,
double *Easting,
double *Northing)const;
/*!
* The function Convert_Sinusoidal_To_Geodetic converts Sinusoidal projection
* easting and northing coordinates to geodetic (latitude and longitude)
* coordinates, according to the current ellipsoid and Sinusoida projection
* coordinates. If any errors occur, the error code(s) are returned by the
* function, otherwise SINU_NO_ERROR is returned.
*
* Easting : Easting (X) in meters (input)
* Northing : Northing (Y) in meters (input)
* Latitude : Latitude (phi) in radians (output)
* Longitude : Longitude (lambda) in radians (output)
*/
long Convert_Sinusoidal_To_Geodetic(double Easting,
double Northing,
double *Latitude,
double *Longitude)const;
TYPE_DATA
};
#endif
|