/usr/include/ossim/projection/ossimTransCylEquAreaProjection.h is in libossim-dev 1.7.21-4.
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 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 | //*******************************************************************
//
// License: See top level LICENSE.txt file.
//
// Author: Garrett Potts (gpotts@imagelinks.com)
//
// Description:
//
// Calls Geotrans Transverse Cylindrical Equal Area
// projection code.
//*******************************************************************
// $Id: ossimTransCylEquAreaProjection.h 9968 2006-11-29 14:01:53Z gpotts $
#ifndef ossimTransCylEquAreaProjection_HEADER
#define ossimTransCylEquAreaProjection_HEADER
#include <ossim/projection/ossimMapProjection.h>
class ossimTransCylEquAreaProjection : public ossimMapProjection
{
public:
ossimTransCylEquAreaProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(6378137, 6356752.3142),
const ossimGpt& origin = ossimGpt());
ossimTransCylEquAreaProjection(const ossimEllipsoid& ellipsoid,
const ossimGpt& origin,
double falseEasting,
double falseNorthing,
double scaleFactor);
~ossimTransCylEquAreaProjection(){}
virtual ossimObject *dup()const{return new ossimTransCylEquAreaProjection(*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);
/*
* Set's scale and then update.
*/
void setScaleFactor(double scaleFactor);
/*!
* 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 setParameters(double falseEasting,
double falseNorthing,
double scaleFactor);
void setDefaults();
double getFalseEasting()const{return Tcea_False_Easting;}
double getFalseNorthing()const{return Tcea_False_Northing;}
double getScaleFactor()const{return Tcea_Scale_Factor;}
/*!
* 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);
private:
mutable double Tcea_a; /* Semi-major axis of ellipsoid in meters */
mutable double Tcea_f; /* Flattening of ellipsoid */
mutable double es2; /* Eccentricity (0.08181919084262188000) squared */
mutable double es4; /* es2 * es2 */
mutable double es6; /* es4 * es2 */
mutable double es; /* sqrt(es2) */
mutable double M0;
mutable double qp;
mutable double One_MINUS_es2; /* 1.0 - es2 */
mutable double One_OVER_2es; /* 1.0 / (2.0 * es) */
mutable double a0; /* es2 / 3.0 + 31.0 * es4 / 180.0 + 517.0 * es6 / 5040.0 */
mutable double a1; /* 23.0 * es4 / 360.0 + 251.0 * es6 / 3780.0 */
mutable double a2; /* 761.0 * es6 / 45360.0 */
mutable double b0; /* 3.0 * e1 / 2.0 - 27.0 * e3 / 32.0 */
mutable double b1; /* 21.0 * e2 / 16.0 - 55.0 * e4 / 32.0 */
mutable double b2; /* 151.0 * e3 / 96.0 */
mutable double b3; /* 1097.0 * e4 / 512.0 */
mutable double c0;/* 1.0 - 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 */
/* Transverse Cylindrical Equal Area projection Parameters */
mutable double Tcea_Origin_Lat; /* Latitude of origin in radians */
mutable double Tcea_Origin_Long; /* Longitude of origin in radians */
mutable double Tcea_False_Northing; /* False northing in meters */
mutable double Tcea_False_Easting; /* False easting in meters */
mutable double Tcea_Scale_Factor; /* Scale factor */
/* Maximum variance for easting and northing values for WGS 84.
*/
mutable double Tcea_Min_Easting;
mutable double Tcea_Max_Easting;
mutable double Tcea_Min_Northing;
mutable double Tcea_Max_Northing;
/*
* The function Set_Trans_Cyl_Eq_Area_Parameters receives the ellipsoid parameters and
* Transverse Cylindrical Equal Area projcetion parameters as inputs, and sets the corresponding
* state variables. If any errors occur, the error code(s) are returned by the function,
* otherwise TCEA_NO_ERROR is returned.
*
* a : Semi-major axis of ellipsoid, in meters (input)
* f : Flattening of ellipsoid (input)
* Origin_Latitude : Latitude in radians at which the (input)
* point scale factor is 1.0
* 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)
* Scale_Factor : Multiplier which reduces distances in the
* projection to the actual distance on the
* ellipsoid (input)
*/
long Set_Trans_Cyl_Eq_Area_Parameters(double a,
double f,
double Origin_Latitude,
double Central_Meridian,
double False_Easting,
double False_Northing,
double Scale_Factor);
/*!
* The function Get_Trans_Cyl_Eq_Area_Parameters returns the current ellipsoid
* parameters, Transverse Cylindrical Equal Area projection parameters, and scale factor.
*
* a : Semi-major axis of ellipsoid, in meters (output)
* f : Flattening of ellipsoid (output)
* Origin_Latitude : Latitude in radians at which the (output)
* point scale factor is 1.0
* 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)
* Scale_Factor : Multiplier which reduces distances in the
* projection to the actual distance on the
* ellipsoid (output)
*/
void Get_Trans_Cyl_Eq_Area_Parameters(double *a,
double *f,
double *Origin_Latitude,
double *Central_Meridian,
double *False_Easting,
double *False_Northing,
double *Scale_Factor)const;
/*!
* The function Convert_Geodetic_To_Trans_Cyl_Eq_Area converts geodetic (latitude and
* longitude) coordinates to Transverse Cylindrical Equal Area projection easting, and northing
* coordinates, according to the current ellipsoid and Transverse Cylindrical Equal Area projection
* parameters. If any errors occur, the error code(s) are returned by the
* function, otherwise TCEA_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_Trans_Cyl_Eq_Area (double Latitude,
double Longitude,
double *Easting,
double *Northing)const;
/*!
* The function Convert_Trans_Cyl_Eq_Area_To_Geodetic converts Transverse Cylindrical Equal Area
* projection easting and northing coordinates to geodetic (latitude and longitude) coordinates.
* coordinates, according to the current ellipsoid and Transverse Cylindrical Equal Area projection
* If any errors occur, the error code(s) are returned by the function, otherwise TCEA_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_Trans_Cyl_Eq_Area_To_Geodetic(double Easting,
double Northing,
double *Latitude,
double *Longitude)const;
TYPE_DATA
};
#endif
|