/usr/include/ossim/projection/ossimAlbersProjection.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 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 | //*******************************************************************
// Copyright (C) 2000 ImageLinks Inc.
//
// License: See top level LICENSE.txt file.
//
// Author: Garrett Potts
//
// Description:
//
// Calls Geotrans Albers projection code.
//*******************************************************************
// $Id: ossimAlbersProjection.h 9094 2006-06-13 19:12:40Z dburken $
#ifndef ossimAlbersProjection_HEADER
#define ossimAlbersProjection_HEADER
#include <ossim/projection/ossimMapProjection.h>
//! Albers Equal Area Conic Projection
/*!
* A basic conic projection developed by H.C. Albers. Classified as a
* conic and equal area. A graphic can be found at:
* http://everest.hunter.cuny.edu/mp/pbander/albers.gif
*/
class ossimAlbersProjection : public ossimMapProjection
{
public:
ossimAlbersProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
const ossimGpt& origin = ossimGpt());
ossimAlbersProjection(const ossimEllipsoid& ellipsoid,
const ossimGpt& origin,
double stdParallel1,
double stdParallel2,
double falseEasting,
double falseNorthing);
virtual ~ossimAlbersProjection();
ossimObject *dup()const;
virtual ossimDpt forward(const ossimGpt &worldPoint) const;
virtual ossimGpt inverse(const ossimDpt &projectedPoint)const;
virtual void update();
/*!
* The passed in Argument is in degrees and will convert it to radians
* for Geotrans. Update is then called so we can pre-compute paramters.
*/
void setStandardParallel1(double degree);
/*!
* The passed in Argument is in degrees and will convert it to radians
* for Geotrans. Update is then called so we can pre-compute paramters
*/
void setStandardParallel2(double degree);
/*!
* Will allow you to set both parallels. both are passed in degrees.
*/
void setStandardParallels(double parallel1Degree,
double prallel2Degree);
/*!
* 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);
/*!
* Allows one to set all parameters for this projections. parallels are in
* degrees and easting northings are in meters.
* Update is then called so we can pre-compute paramters.
*/
void setParameters(double parallel1,
double parallel2,
double falseEasting,
double falseNorthing);
/** @returns the first standard parallel in decimal degrees. */
virtual double getStandardParallel1() const;
/** @returns the second standard parallel in decimal degrees. */
virtual double getStandardParallel2() const;
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 Albers_a; /* Semi-major axis of ellipsoid in meters */
mutable double Albers_f; /* Flattening of ellipsoid */
mutable double es; /* Eccentricity of ellipsoid */
mutable double es2; /* Eccentricity squared */
mutable double C; /* constant c */
mutable double rho0; /* height above ellipsoid */
mutable double n; /* ratio between meridians */
mutable double Albers_a_OVER_n; /* Albers_a / n */
mutable double one_MINUS_es2; /* 1 - es2 */
mutable double two_es; /* 2 * es */
/* Albers Projection Parameters */
mutable double Albers_Origin_Lat; /* Latitude of origin in radians */
mutable double Albers_Origin_Long; /* Longitude of origin in radians */
mutable double Albers_Std_Parallel_1;
mutable double Albers_Std_Parallel_2;
mutable double Albers_False_Easting;
mutable double Albers_False_Northing;
mutable double Albers_Delta_Northing;
mutable double Albers_Delta_Easting;
/*!
* The function Set_Albers_Parameters receives the ellipsoid parameters and
* projection parameters as inputs, and sets the corresponding state
* variables. If any errors occur, the error code(s) are returned by the function,
* otherwise ALBERS_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 (central meridian)
* Std_Parallel_1 : First standard parallel (input)
* Std_Parallel_2 : Second standard parallel (input)
* 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_Albers_Parameters(double a,
double f,
double Origin_Latitude,
double Central_Meridian,
double Std_Parallel_1,
double Std_Parallel_2,
double False_Easting,
double False_Northing);
/*!
* The function Get_Albers_Parameters returns the current ellipsoid
* parameters, and Albers projection parameters.
*
* 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
* Std_Parallel_1 : First standard parallel (output)
* Std_Parallel_2 : Second standard parallel (output)
* 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_Albers_Parameters(double *a,
double *f,
double *Origin_Latitude,
double *Central_Meridian,
double *Std_Parallel_1,
double *Std_Parallel_2,
double *False_Easting,
double *False_Northing)const;
/*!
* The function Convert_Geodetic_To_Albers converts geodetic (latitude and
* longitude) coordinates to Albers projection (easting and northing)
* coordinates, according to the current ellipsoid and Albers projection
* parameters. If any errors occur, the error code(s) are returned by the
* function, otherwise ALBERS_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_Albers (double Latitude,
double Longitude,
double *Easting,
double *Northing)const;
/*!
* The function Convert_Albers_To_Geodetic converts Albers projection
* (easting and northing) coordinates to geodetic (latitude and longitude)
* coordinates, according to the current ellipsoid and Albers projection
* coordinates. If any errors occur, the error code(s) are returned by the
* function, otherwise ALBERS_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_Albers_To_Geodetic(double Easting,
double Northing,
double *Latitude,
double *Longitude)const;
/*
* RTTI information defined in data_types/rtti.h
*/
TYPE_DATA
};
#endif
|