/usr/include/ossim/projection/ossimEquDistCylProjection.h is in libossim-dev 1.8.16-3+b1.
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 217 | //*******************************************************************
//
// License: See top level LICENSE.txt file.
//
// Author: Garrett Potts
//
// Description:
//
// Calls Geotrans Equidistant Cylinder projection code.
//*******************************************************************
// $Id: ossimEquDistCylProjection.h 20352 2011-12-12 17:24:52Z dburken $
#ifndef ossimEquDistCylProjection_HEADER
#define ossimEquDistCylProjection_HEADER
#include <ossim/projection/ossimLlxyProjection.h>
class ossimIpt;
class OSSIMDLLEXPORT ossimEquDistCylProjection : public ossimMapProjection//public ossimLlxyProjection
{
public:
ossimEquDistCylProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
const ossimGpt& origin = ossimGpt());
ossimEquDistCylProjection(const ossimEllipsoid& ellipsoid,
const ossimGpt& origin,
double falseEasting,
double falseNorthing);
ossimEquDistCylProjection(const ossimEquDistCylProjection& rhs)
:ossimMapProjection(rhs)
// :ossimLlxyProjection(rhs)
{
*this = rhs;
}
~ossimEquDistCylProjection(){}
virtual ossimObject *dup()const
{
return new ossimEquDistCylProjection(*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);
double getFalseEasting()const{return Eqcy_False_Easting;}
double getFalseNorthing()const{return Eqcy_False_Northing;}
virtual void lineSampleHeightToWorld(const ossimDpt& lineSampPt,
const double& heightAboveEllipsoid,
ossimGpt& worldPt) const;
virtual void worldToLineSample(const ossimGpt &worldPoint,
ossimDpt& lineSample)const;
/**
* @brief Specialized worldToLineSample.
*
* Takes into consideration the image size and crossing the date line.
* Also makes local points outside of image negative or positive based on
* there relative closeness to left or right edge.
*
* @param worldPoint Ground point to convert to line sample.
* @param imageSize Size of image in pixels.
* @param lineSample Initialized by this.
*/
void worldToLineSample( const ossimGpt& worldPoint,
const ossimIpt& imageSize,
ossimDpt& lineSample ) const;
virtual bool isGeographic()const
{
return true;
}
void setDefaults();
/*!
* 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:
//---------------------GEOTRANS-------------------------------
mutable double Eqcy_a; /* Semi-major axis of ellipsoid in meters */
mutable double Eqcy_f; /* Flattening of ellipsoid */
mutable double es2; /* Eccentricity (0.08181919084262188000) squared */
mutable double es4; /* es2 * es2 */
mutable double es6; /* es4 * es2 */
mutable double Ra; /* Spherical Radius */
/* Equidistant Cylindrical projection Parameters */
mutable double Eqcy_Std_Parallel; /* Latitude of standard parallel in radians */
mutable double Cos_Eqcy_Std_Parallel; /* cos(Eqcy_Std_Parallel) */
mutable double Eqcy_Origin_Long; /* Longitude of origin in radians */
mutable double Eqcy_False_Easting;
mutable double Eqcy_False_Northing;
mutable double Eqcy_Delta_Northing;
mutable double Eqcy_Max_Easting;
mutable double Eqcy_Min_Easting;
mutable double Ra_Cos_Eqcy_Std_Parallel; /* Ra * Cos_Eqcy_Std_Parallel */
/*!
* The function Set_Equidistant_Cyl_Parameters receives the ellipsoid parameters and
* projection parameters as inputs, and sets the corresponding state
* variables. It also calculates the spherical radius of the sphere having
* the same area as the ellipsoid. If any errors occur, the error code(s)
* are returned by the function, otherwise EQCY_NO_ERROR is returned.
*
* a : Semi-major axis of ellipsoid, in meters (input)
* f : Flattening of ellipsoid (input)
* Std_Parallel : 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
* standard parallel of the projection (input)
*/
long Set_Equidistant_Cyl_Parameters(double a,
double f,
double Std_Parallel,
double Central_Meridian,
double False_Easting,
double False_Northing);
/*!
* The function Get_Equidistant_Cyl_Parameters returns the current ellipsoid
* parameters and Equidistant Cylindrical projection parameters.
*
* a : Semi-major axis of ellipsoid, in meters (output)
* f : Flattening of ellipsoid (output)
* Std_Parallel : 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
* standard parallel of the projection (output)
*/
void Get_Equidistant_Cyl_Parameters(double *a,
double *f,
double *Std_Parallel,
double *Central_Meridian,
double *False_Easting,
double *False_Northing)const;
/*!
* The function Convert_Geodetic_To_Equidistant_Cyl converts geodetic (latitude and
* longitude) coordinates to Equidistant Cylindrical projection easting, and northing
* coordinates, according to the current ellipsoid and Equidistant Cylindrical projection
* parameters. If any errors occur, the error code(s) are returned by the
* function, otherwise EQCY_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_Equidistant_Cyl (double Latitude,
double Longitude,
double *Easting,
double *Northing)const;
/*!
* The function Convert_Equidistant_Cyl_To_Geodetic converts Equidistant Cylindrical projection
* easting and northing coordinates to geodetic (latitude and longitude)
* coordinates, according to the current ellipsoid and Equidistant Cylindrical projection
* coordinates. If any errors occur, the error code(s) are returned by the
* function, otherwise EQCY_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_Equidistant_Cyl_To_Geodetic(double Easting,
double Northing,
double *Latitude,
double *Longitude)const;
TYPE_DATA
};
#endif
|