ossimEpsgProjectionDatabase Class Reference

Projection Database for coded projections defined in database files and specified via some coordinate reference system group:code, such as EPSG:26715. More...

#include <ossimEpsgProjectionDatabase.h>

Inheritance diagram for ossimEpsgProjectionDatabase:

Classes
class	ProjDbRecord
	Type for database record consists of EPSG code and serialized form of corresponding OSSIM projection (as a keywordlist) More...

Public Member Functions
virtual	~ossimEpsgProjectionDatabase ()
	Destructor. More...
ossimProjection *	findProjection (const ossimString &proj_spec) const
	Returns a projection corresponding to the projection specified, or NULL if no entry found. More...
ossimProjection *	findProjection (ossim_uint32 epsg_code) const
	Returns a projection corresponding to the epsg code specified, or NULL if no entry found. More...
ossim_uint32	findProjectionCode (const ossimString &projection_name) const
	Given a projection name, assigns the group (e.g., "EPSG") and code of the projection. More...
ossim_uint32	findProjectionCode (const ossimMapProjection &lost_proj) const
	Given a projection instance, this method determines the corresponding EPSG code. More...
ossimString	findProjectionName (ossim_uint32 epsg_code) const
	NECESSARY HACK ### (originally implemented by ESH) ArcMap (9.2 and less) doesn't understand the non-meters HARN codes. More...
void	getProjectionsList (std::vector< ossimString > &typeList) const
	Populates caller's list with all projections currently represented in the database. More...
size_t	numRecords () const
	ENGINEERING CODE. Used for testing. More...
ossim_uint32	getCodeFromUtmProj (const ossimUtmProjection *proj) const
	Given UTM projection, derives the associated EPSG code. This is faster than a Db lookup. More...
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

Static Public Member Functions
static ossimEpsgProjectionDatabase *	instance ()
	Instantiates singleton instance of this class: More...

Protected Types
enum	RecordFormat {   NOT_ASSIGNED =0, FORMAT_A, FORMAT_B, FORMAT_C,   CUSTOM }

Protected Member Functions
	ossimEpsgProjectionDatabase ()
	Constructor loads all Db files specified in the ossim prefs. More...
ossimMapProjection *	createProjFromFormatARecord (ProjDbRecord *record) const
	Parses the "Spadac EPSG" Db record format and produces a projection (or NULL if invalid) More...
ossimMapProjection *	createProjFromFormatBRecord (ProjDbRecord *record) const
	Parses the State Plane Db record format and produce a projection (or NULL if invalid) More...
ossimMapProjection *	createProjFromUtmCode (ossim_uint32 code) const
	HACK UTM projections as specified in the EPSG are indistinguishable from regular TM. More...
void	initialize () const
	Populates the database with contents of DB files as specified in ossim_preferences. More...
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
std::multimap< ossim_uint32, ossimRefPtr< ProjDbRecord > >	m_projDatabase
std::mutex	m_mutex

Friends
class	ossimEpsgProjectionFactory

Detailed Description

Projection Database for coded projections defined in database files and specified via some coordinate reference system group:code, such as EPSG:26715.

Definition at line 29 of file ossimEpsgProjectionDatabase.h.

Member Enumeration Documentation

RecordFormat

enum ossimEpsgProjectionDatabase::RecordFormat

protected

Enumerator
NOT_ASSIGNED
FORMAT_A
FORMAT_B
FORMAT_C
CUSTOM

Definition at line 77 of file ossimEpsgProjectionDatabase.h.

   {  
      NOT_ASSIGNED=0,
      FORMAT_A,
      FORMAT_B,
      FORMAT_C,
      CUSTOM
   };

Constructor & Destructor Documentation

~ossimEpsgProjectionDatabase()

ossimEpsgProjectionDatabase::~ossimEpsgProjectionDatabase ( )

virtual

Destructor.

Definition at line 173 of file ossimEpsgProjectionDatabase.cpp.

{
}

ossimEpsgProjectionDatabase()

ossimEpsgProjectionDatabase::ossimEpsgProjectionDatabase ( )

protected

Constructor loads all Db files specified in the ossim prefs.

Constructor loads all DB CSV files specified in the ossim prefs.

Protected as part of singleton implementation.

Definition at line 180 of file ossimEpsgProjectionDatabase.cpp.

   :
   m_projDatabase(),
   m_mutex()
{
}

Member Function Documentation

createProjFromFormatARecord()

ossimMapProjection * ossimEpsgProjectionDatabase::createProjFromFormatARecord ( ProjDbRecord *  record ) const

protected

Parses the "Spadac EPSG" Db record format and produces a projection (or NULL if invalid)

This method is invoked to parse the Db record and produce a projection (or NULL if spec invalid)

Definition at line 624 of file ossimEpsgProjectionDatabase.cpp.

References A_DATUM_CODE, A_FALSE_EASTING, A_FALSE_NORTHING, A_FALSE_ORG_EASTING, A_FALSE_ORG_LAT, A_FALSE_ORG_LON, A_FALSE_ORG_NORTHING, A_NAT_ORG_LAT, A_NAT_ORG_LON, A_NAT_ORG_SCALE, A_PROJ_TYPE, A_STD_PARL_1_LAT, A_STD_PARL_2_LAT, A_UNITS, ossimEpsgProjectionDatabase::ProjDbRecord::code, ossimString::contains(), ossimEpsgDatumFactory::create(), ossimDatumFactory::create(), ossimEpsgProjectionDatabase::ProjDbRecord::csvRecord, ossimEpsgProjectionDatabase::ProjDbRecord::datumValid, decodeSexagesimalDms(), ossimDatum::ellipsoid(), ossimRefPtr< T >::get(), ossimEpsgDatumFactory::instance(), ossimDatumFactory::instance(), ossimGpt::lat, ossimGpt::lon, MTRS_PER_FT, OSSIM_FEET, OSSIM_KILOMETERS, OSSIM_METERS, OSSIM_US_SURVEY_FEET, ossimEpsgProjectionDatabase::ProjDbRecord::proj, ossimMapProjection::setDatum(), ossimNewZealandMapGridProjection::setFalseEastingNorthing(), ossimMapProjection::setOrigin(), ossimMapProjection::setPcsCode(), ossimMapProjection::setProjectionUnits(), and US_METERS_PER_FT.

Referenced by findProjection().

{
   // Establish EPSG code and test for UTM (full projection is implied in the code itself -- no 
   // accessing the database). Until the database is solidified, it is probably better to do 
   // it this way:
   record->datumValid = true;
   record->proj = 0;

   // Establish the units in which the easting/northing is provided:
   ossimUnitType mapProjUnits = OSSIM_METERS;
   double mtrs_per_unit = 1.0;
   if (record->csvRecord[A_UNITS] == "US survey foot")
   {
      mapProjUnits = OSSIM_US_SURVEY_FEET;
      mtrs_per_unit = US_METERS_PER_FT;
   }
   else if (record->csvRecord[A_UNITS].contains("foot")) // THIS IS INTERNATIONAL FOOT, NOT EXACT FOR MANY INTERNATIONAL VARIETIES
   {
      mapProjUnits = OSSIM_FEET;
      mtrs_per_unit = MTRS_PER_FT;
   }
   else if (record->csvRecord[A_UNITS].contains("kilometre"))
   {
      mapProjUnits = OSSIM_KILOMETERS;
      mtrs_per_unit = 1000.0;
   }
   else if (!record->csvRecord[A_UNITS].contains("metre"))
   {
      // ### SKIP THIS MESSAGE BUT BE AWARE THAT THIS PROJECTION WON'T BE REPRESENTED IN DB ###
      //ossimNotify(ossimNotifyLevel_WARN)<<MODULE<<"EPSG:"<<pcs_code<<" units of <"
      //   <<record->csvRecord[A_UNITS]<<"> not presently supported."<<endl;
      return 0;
   }

   // First create a datum given the datum code in the record:
   ossim_uint32 gcs_code = record->csvRecord[A_DATUM_CODE].toUInt32();
   const ossimDatum* datum = ossimEpsgDatumFactory::instance()->create(gcs_code);
   if (!datum)
   {
      // Default to WGS 84 -- this may throw an exception:
      datum = ossimDatumFactory::instance()->create(ossimString("WGE"));
      record->datumValid = false;
   }
   const ossimEllipsoid* ellipsoid = datum->ellipsoid();

   ossimGpt origin(0,0,0,datum);
   ossimString proj_type = record->csvRecord[A_PROJ_TYPE];
   if (proj_type.contains("Transverse Mercator"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      double sf = record->csvRecord[A_NAT_ORG_SCALE].toDouble();
      record->proj = new ossimTransMercatorProjection(*ellipsoid, origin, fe, fn, sf);
   }
   else if (proj_type.contains("Lambert Conic Conformal (1SP)"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      record->proj = new ossimLambertConformalConicProjection(*ellipsoid, origin, origin.lat, 
         origin.lat, fe, fn);
   }
   else if (proj_type.contains("Lambert Conic Conformal (2SP)"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_FALSE_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_FALSE_ORG_LON]);
      double p1 = decodeSexagesimalDms(record->csvRecord[A_STD_PARL_1_LAT]);
      double p2 = decodeSexagesimalDms(record->csvRecord[A_STD_PARL_2_LAT]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_ORG_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_ORG_NORTHING].toDouble();
      record->proj = new ossimLambertConformalConicProjection(*ellipsoid, origin, p1, p2, fe, fn);
   }
   else if (proj_type.contains("Cassini"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      record->proj =  new ossimCassiniProjection(*ellipsoid, origin, fe, fn);
   }
   else if (proj_type.contains("Mercator (1SP)"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      double sf = record->csvRecord[A_NAT_ORG_SCALE].toDouble();
      record->proj = new ossimMercatorProjection(*ellipsoid, origin, fe, fn, sf);
   }
   else if(proj_type.contains("Popular Visualisation Sphere"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      double sf = record->csvRecord[A_NAT_ORG_SCALE].toDouble();
      record->proj = new ossimMercatorProjection(*ellipsoid, origin, fe, fn, sf);

      // Set this for saveState:
      record->proj->setPcsCode( 3785 );

   }
   else if(proj_type.contains("Pseudo Mercator"))
   {
      record->proj = new ossimGoogleProjection();

      // Set this for saveState:
      record->proj->setPcsCode( 3857 );

   }
   else if (proj_type.contains("Albers"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_FALSE_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_FALSE_ORG_LON]);
      double p1 = decodeSexagesimalDms(record->csvRecord[A_STD_PARL_1_LAT]);
      double p2 = decodeSexagesimalDms(record->csvRecord[A_STD_PARL_2_LAT]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_ORG_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_ORG_NORTHING].toDouble();
      record->proj = new ossimAlbersProjection(*ellipsoid, origin, p1, p2, fe, fn);
   }
   else if (proj_type.contains("Equidistant Cylindrical"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_STD_PARL_1_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      record->proj = new ossimEquDistCylProjection(*ellipsoid, origin, fe, fn);
   }
   else if (proj_type.contains("New Zealand Map Grid"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      ossimNewZealandMapGridProjection* nz_proj = new ossimNewZealandMapGridProjection;
      nz_proj->setOrigin(origin);
      nz_proj->setFalseEastingNorthing(fe, fn);
      record->proj = nz_proj;
   }
   else if (proj_type.contains("Polar Stereographic"))
   {
      origin.lat = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LAT]);
      origin.lon = decodeSexagesimalDms(record->csvRecord[A_NAT_ORG_LON]);
      // double fe = mtrs_per_unit*record->csvRecord[A_FALSE_EASTING].toDouble();
      // double fn = mtrs_per_unit*record->csvRecord[A_FALSE_NORTHING].toDouble();
      ossimUpsProjection* ups_proj = new ossimUpsProjection(ossimEllipsoid(), origin);
      record->proj = ups_proj;
   }
   else
   {
      // Can't handle it now. 
      //ossimNotify(ossimNotifyLevel_FATAL)<<MODULE<<"EPSG:"<<record->csvRecord[A_CODE]<<" \""<<proj_type<<"\" "
      //   "not presently supported. Returning NULL projection..."<<endl;
      return 0;
   }
  
   record->proj->setDatum(datum);
   record->proj->setPcsCode(record->code);
   record->proj->setProjectionUnits(mapProjUnits);

   return record->proj.get();
}

createProjFromFormatBRecord()

ossimMapProjection * ossimEpsgProjectionDatabase::createProjFromFormatBRecord ( ProjDbRecord *  db_record ) const

protected

Parses the State Plane Db record format and produce a projection (or NULL if invalid)

State Plane projections as specified in the EPSG are indistinguishable from regular TM.

Unfortunately OSSIM does make a distinction. For the time being, parse the code and programmatically arrive at the UTM projection.

Definition at line 797 of file ossimEpsgProjectionDatabase.cpp.

References B_CODE, B_DATUM_CODE, B_FALSE_EASTING, B_FALSE_NORTHING, B_PARAM1, B_PARAM2, B_PARAM3, B_PARAM4, B_PROJ_TYPE, B_UNITS, ossimEpsgProjectionDatabase::ProjDbRecord::code, ossimEpsgProjectionDatabase::ProjDbRecord::csvRecord, ossimEpsgProjectionDatabase::ProjDbRecord::datumValid, ossimRefPtr< T >::get(), ossimMapProjectionFactory::instance(), ossimEpsgProjectionDatabase::ProjDbRecord::name, ossimKeywordNames::PCS_CODE_KW, ossimEpsgProjectionDatabase::ProjDbRecord::proj, PTR_CAST, ossimKeywordlist::remove(), ossimMapProjection::setPcsCode(), and ossimRefPtr< T >::valid().

Referenced by findProjection().

{
   // Unfortunately, as of this writing, there is a disconnect between the Geotrans projection
   // parameters and those provided by EPSG database. In the meantime, to maintain functionality,
   // we intercept the EPSG code and programmatically arrive at the projection versus pulling
   // all the parameters out of the Db record->
   ossimStatePlaneProjectionInfo  info(db_record->name,
                                       db_record->code,
                                       db_record->csvRecord[B_PROJ_TYPE],
                                       db_record->csvRecord[B_PARAM1],
                                       db_record->csvRecord[B_PARAM2],
                                       db_record->csvRecord[B_PARAM3],
                                       db_record->csvRecord[B_PARAM4],
                                       db_record->csvRecord[B_FALSE_EASTING].toDouble(),
                                       db_record->csvRecord[B_FALSE_NORTHING].toDouble(),
                                       db_record->csvRecord[B_UNITS],
                                       db_record->csvRecord[B_DATUM_CODE]);
   
   // NOTE: In order to avoid infinite recursion with this object, we initialized the PCS code in 
   // info to NULL to insure that the projection is instantiated directly (not via this class):
   ossimKeywordlist kwl;
   info.populateProjectionKeywords(kwl);
   kwl.remove(ossimKeywordNames::PCS_CODE_KW);
   db_record->proj = 
      PTR_CAST(ossimMapProjection, ossimMapProjectionFactory::instance()->createProjection(kwl));
   if (db_record->proj.valid())
   {
      db_record->proj->setPcsCode(db_record->csvRecord[B_CODE].toUInt32());
      db_record->datumValid = true;
   }

   return db_record->proj.get();
}

createProjFromUtmCode()

ossimMapProjection * ossimEpsgProjectionDatabase::createProjFromUtmCode ( ossim_uint32  code ) const

protected

HACK

UTM projections as specified in the EPSG are indistinguishable from regular TM.

UTM projections as specified in the EPSG are indistinguishable from regular TM.

Unfortunately OSSIM does make a distinction. For the time being, parse the code and programmatically arrives at the UTM projection.

Unfortunately OSSIM does make a distinction. For the time being, parse the code and programmatically arrive at the UTM projection.

Definition at line 837 of file ossimEpsgProjectionDatabase.cpp.

References ossimDatumFactory::instance(), ossimMapProjection::setDatum(), ossimUtmProjection::setHemisphere(), ossimMapProjection::setPcsCode(), and ossimUtmProjection::setZone().

Referenced by findProjection().

{
   // This code originally found in deprecated ossimPcsCodeProjectionFactory authored by
   // David Burken.
   int type = code/100;
   int zone = code%100;
   ossimUtmProjection* proj = NULL;
   switch (type)
   {
   case 322:
      // utm, WGS72 (WGD), northern hemisphere
      if ( (zone > 0 ) && (zone < 61) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->wgs72()->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->wgs72());
         proj->setZone(zone);
         proj->setHemisphere('N');
      }
      break;

   case 323:
      // utm, WGS72 (WGD), southern hemisphere
      if ( (zone > 0) && (zone < 61) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->wgs72()->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->wgs72());
         proj->setZone(zone);
         proj->setHemisphere('S');
      }
      break;

   case 326:
      // utm, WGS84 (WGE), northern hemisphere
      if ( (zone > 0) && (zone < 61) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->wgs84()->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->wgs84());
         proj->setZone(zone);
         proj->setHemisphere('N');
      }
      break;

   case 327:
      // utm, WGS84 (WGE), southern hemisphere
      if ( (zone > 0) && (zone < 61) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->wgs84()->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->wgs84());
         proj->setZone(zone);
         proj->setHemisphere('S');
      }
      break;

   case 267:
      // utm, "NAS-C", northern hemisphere
      // Only UTM NAD27 North zones 3 to 22 are in the 267xx range...
      // 26729 through 26803 handled by state plane factory.
      //---
      if ( (code > 26702) && (code < 26723) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->
            create(ossimString("NAS-C"))->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->create(ossimString("NAS-C")));
         proj->setZone(zone);
         proj->setHemisphere('N');
      }
      break;

   case 269:
      // utm, "NAR-C", northern hemisphere
      // Only UTM NAD83 North zones 3 to 23 are in the 269xx range...
      // 26929 through 26998 handled by state plane factory.
      if ( (code > 26902) && (code < 26924) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->
            create(ossimString("NAR-C"))->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->create(ossimString("NAR-C")));
         proj->setZone(zone);
         proj->setHemisphere('N');
      }
      break;

   case 248:
      // Provisional S. American 1956, 24818 through 24880
      if ( (code > 24817) && (code < 24881) )
      {
         proj = new ossimUtmProjection(*(ossimDatumFactory::instance()->
            create(ossimString("PRP-M"))->ellipsoid()));
         proj->setDatum(ossimDatumFactory::instance()->create(ossimString("PRP-M")));
         if (zone > 60)
         {
            proj->setZone(zone - 60);
            proj->setHemisphere('S');
         }
         else
         {
            proj->setZone(zone);
            proj->setHemisphere('N');
         }
      }
      break;

   } // End of switch on code.

   if (proj)
      proj->setPcsCode(code);

   return proj;
}

findProjection() [1/2]

ossimProjection * ossimEpsgProjectionDatabase::findProjection

(

const ossimString &

spec

)

const

Returns a projection corresponding to the projection specified, or NULL if no entry found.

Returns a projection corresponding to the group (e.g., "EPSG") and PCS code provided, or NULL if no entry found.

Normally <proj_spec> takes the form of <group>: (e.g., "EPSG:26715"). Other forms accepted are (assumed EPSG), or <proj_name> to match a name in the database.

Definition at line 420 of file ossimEpsgProjectionDatabase.cpp.

References ossimString::after(), ossimString::before(), ossimEpsgProjectionDatabase::ProjDbRecord::code, ossimString::contains(), ossimDatumFactoryRegistry::create(), ossimString::downcase(), ossimProjection::dup(), ossimDatum::ellipsoid(), initialize(), ossimDatumFactoryRegistry::instance(), m_projDatabase, ossimEpsgProjectionDatabase::ProjDbRecord::name, ossimEpsgProjectionDatabase::ProjDbRecord::proj, ossimMapProjection::setDatum(), ossimMapProjection::setPcsCode(), ossimString::split(), ossimString::toUInt32(), and ossimRefPtr< T >::valid().

Referenced by ossimEpsgProjectionFactory::createProjection(), ossimWktProjectionFactory::createProjection(), and findProjectionCode().

{
   //std::cout << "ossimEpsgProjectionDatabase::findProjection: entered with spec " << spec << "\n";
   ossimProjection* proj = 0;

   // Use the CRS code to access the database. The spec should be <group>:<code> where <group> is 
   // "EPSG" (the only group handled here):
   ossim_uint32 spec_code;
   ossimString spec_group ("epsg"); // default if only integer code provided
   if (spec.contains(":"))
   {
      spec_group = spec.before(":");
      spec_code = spec.after(":").toUInt32();
      spec_group = spec_group.downcase();
   }
   else
   {
      spec_code = spec.toUInt32();
   }

   // Presently only EPSG database is handled:
   if ((spec_code != 0) && (spec_group == "epsg"))
      return findProjection(spec_code);

   // The spec is probably a projection name. Need to search Db by the projection name. 
   // Search database for entry. The spec may use different delimiters than
   // the DB so need to split the strings and compare the words:
   ossimString separators ("_ /()");
   vector<ossimString> split_spec = spec.split(separators, true);
   vector<ossimString> split_db_name;
   ossimRefPtr<ossimMapProjection> map_proj = 0;
   if (m_projDatabase.empty())
      initialize();
   std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter = m_projDatabase.begin();
   while ((db_iter != m_projDatabase.end()) && !proj)
   {
      ossimRefPtr<ProjDbRecord> db_record = db_iter->second;
      if ( db_record.valid() )
      {
         split_db_name.clear();
         db_record->name.split(split_db_name, separators, true);
         if (split_spec == split_db_name)
         {
            // We may already have instantiated this projection, in which case just return its copy.
            // Otherwise, create the projection from the EPSG code that corresponds to the name:
            if (db_record->proj.valid())
               proj = (ossimMapProjection*) db_record->proj->dup();
            else
               proj = findProjection(db_record->code);
            return proj;
         }
      }
      ++db_iter;
   }
    
   // No hit? Could be that just a datum was identified, in which case we need a simple 
   // Platte Carree:
   const ossimDatum* datum = ossimDatumFactoryRegistry::instance()->create(spec);
   if (datum)
   {
      ossimEquDistCylProjection* proj = new ossimEquDistCylProjection(*(datum->ellipsoid()));
      proj->setDatum(datum);
      proj->setPcsCode(spec_code);
      return proj;
   }

   return 0;
}

findProjection() [2/2]

ossimProjection * ossimEpsgProjectionDatabase::findProjection

(

ossim_uint32

epsg_code

)

const

Returns a projection corresponding to the epsg code specified, or NULL if no entry found.

Returns a projection corresponding to the group (e.g., "EPSG") and PCS code provided, or NULL if no entry found.

Definition at line 332 of file ossimEpsgProjectionDatabase.cpp.

References createProjFromFormatARecord(), createProjFromFormatBRecord(), createProjFromUtmCode(), ossimEpsgProjectionDatabase::ProjDbRecord::csvFormat, ossimEpsgProjectionDatabase::ProjDbRecord::csvRecord, ossimEpsgProjectionDatabase::ProjDbRecord::datumValid, ossimProjection::dup(), FORMAT_A, FORMAT_B, ossimRefPtr< T >::get(), initialize(), m_mutex, m_projDatabase, NOT_ASSIGNED, ossimEpsgProjectionDatabase::ProjDbRecord::proj, ossimMapProjection::setPcsCode(), and ossimRefPtr< T >::valid().

{
   ossimMapProjection* proj = 0;

   // Quick check for bogus EPSG:
   if ((epsg_code == 0) || (epsg_code == 32767))
      return 0;

   // Check for Google projection:
   else if (( epsg_code==3857) || (epsg_code == 900913))
   {
      proj = new ossimGoogleProjection();

      // Set this for saveState:
      proj->setPcsCode( 3857 );
      /*
      const ossimDatum* datum = ossimDatumFactory::instance()->create(ossimString("6055"));
      ossimMercatorProjection* merc_proj = new ossimMercatorProjection();
      ossimGpt origin(0.0,0.0,0.0, datum);
      merc_proj->setFalseEasting(0.0);
      merc_proj->setFalseNorthing(0.0);
      merc_proj->setOrigin(origin); // Also sets the projections datum to the origin's datum
      merc_proj->update();
      merc_proj->setPcsCode(900913);
      proj = merc_proj;
      */
   }

   else
   {
      // Search database for entry:
      m_mutex.lock();
      if (m_projDatabase.empty())
      {
         initialize();
      }
      m_mutex.unlock();

      std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter = 
         m_projDatabase.find(epsg_code);
      if (db_iter != m_projDatabase.end())
      {
         // See if a projection has already been created for this entry:
         m_mutex.lock();
         ossimRefPtr<ProjDbRecord> db_record = db_iter->second;
         if ( db_record.valid() )
         {
            if (db_record->proj.valid())
               proj = (ossimMapProjection*) db_record->proj->dup();
            else
            {
               // Try decoding the EPSG code before accessing DB:
               proj = createProjFromUtmCode(epsg_code);
               if (proj)
               {
                  db_record->proj = proj;
                  db_record->datumValid = true;
               }
               else if (db_iter->second->csvFormat == FORMAT_A)
               {
                  proj = createProjFromFormatARecord( db_record.get() );
               }
               else if (db_iter->second->csvFormat == FORMAT_B)
               {
                  proj = createProjFromFormatBRecord( db_record.get() );
               }
               
               if (proj)
               {
                  // To save allocated memory, get rid of the original CSV entry since a real 
                  // projection is now represented in the database:
                  db_record->csvRecord.clear();
                  db_record->csvFormat = NOT_ASSIGNED;
               }
            }

         } // Matches: if  ( db_record.valid() )
         m_mutex.unlock();
      }
   }

   return proj;
}

findProjectionCode() [1/2]

ossim_uint32 ossimEpsgProjectionDatabase::findProjectionCode

(

const ossimString &

proj_name

)

const

Given a projection name, assigns the group (e.g., "EPSG") and code of the projection.

Returns integer EPSG code if match was found or 0 if not found.

Returns "<group>:<code>" if match was found or empty string if not found.

Definition at line 494 of file ossimEpsgProjectionDatabase.cpp.

References ossimEpsgProjectionDatabase::ProjDbRecord::code, ossimRefPtr< T >::get(), initialize(), m_projDatabase, ossimEpsgProjectionDatabase::ProjDbRecord::name, and ossimRefPtr< T >::valid().

Referenced by ossimOgcWktTranslator::fromOssimKwl(), and ossimMapProjection::getPcsCode().

{
   if (m_projDatabase.empty())
      initialize();
   std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter = m_projDatabase.begin();
   while (db_iter != m_projDatabase.end())
   {
      ossimRefPtr<ProjDbRecord> db_record = db_iter->second.get();
      if ( db_record.valid() )
      {
         if (db_record->name == proj_name)
            return (db_record->code);
      }
      ++db_iter;
   }
      
   return 0;
}

findProjectionCode() [2/2]

ossim_uint32 ossimEpsgProjectionDatabase::findProjectionCode

(

const ossimMapProjection &

lost_proj

)

const

Given a projection instance, this method determines the corresponding EPSG code.

Obviously this is only needed if the projection does not have its PCS code assigned (it is NULL). This happens when the projection was constructed with full parameters instead of the EPSG code. Returns integer EPSG code if match was found or 0 if not found.

Definition at line 520 of file ossimEpsgProjectionDatabase.cpp.

References ossimEpsgProjectionDatabase::ProjDbRecord::code, ossimDatum::epsgCode(), findProjection(), ossimRefPtr< T >::get(), ossimObject::getClassName(), getCodeFromUtmProj(), ossimMapProjection::getDatum(), initialize(), m_projDatabase, ossimEpsgProjectionDatabase::ProjDbRecord::proj, and ossimRefPtr< T >::valid().

{
   ossimString lost_type (lost_proj.getClassName());

   // Shortcut for EPSG:4326 (WGS-85 geographic rectangular -- very common):
   if ((lost_type == "ossimEquDistCylProjection") && (lost_proj.getDatum()->epsgCode() == 6326))
      return 4326;

   ossim_uint32 found_code = 0;
   if (lost_type == "ossimUtmProjection")
   {
      found_code = getCodeFromUtmProj(dynamic_cast<const ossimUtmProjection*>(&lost_proj));
      if (found_code)
         return found_code;
   }

   if (m_projDatabase.empty())
      initialize();
   ossimString lookup;
   std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter =
      m_projDatabase.begin();
   while ((db_iter != m_projDatabase.end()) && (found_code == 0))
   {
      ossimRefPtr<ProjDbRecord> db_record = db_iter->second;
      if ( db_record.valid() )
      {
         // Has a projection already been created for this db iter?
         if (!db_record->proj.valid())
         {
            // No projection has been created yet for this DB entry. 
            // NOTE: THIS IS VERY SLOW BECAUSE WE ARE INSTANTIATING EVERY PROJECTION IN THE DB!!!
            db_record->proj = dynamic_cast<ossimMapProjection*>(findProjection(db_record->code));
         }
         if (db_record->proj.valid() && (*(db_record->proj.get()) == lost_proj))
         {
            found_code = db_record->code;
            
            // Hack to remap projection code 4087 to 4326 (which is not really a projection 
            // code but other packages like to see 4326 for geographic projections.
            // Hacked under protest (OLK, 08/2010)
            if (found_code == 4087)
               found_code = 4326;
         }
      }
      ++db_iter;
   }
   return found_code;
}

findProjectionName()

ossimString ossimEpsgProjectionDatabase::findProjectionName

(

ossim_uint32

epsg_code

)

const

NECESSARY HACK ### (originally implemented by ESH)

ArcMap (9.2 and less) doesn't understand the non-meters HARN codes.

This method provides acess to the projection name given the EPSG code.

This method provides acess to the projection name given the EPSG code. The name provides enough info for recognizing a non-meters HARN-based projection.

The name provides enough info for recognizing a non-meters HARN-based projection.

Definition at line 573 of file ossimEpsgProjectionDatabase.cpp.

References initialize(), and m_projDatabase.

Referenced by ossimOgcWktTranslator::fromOssimKwl().

{
   if (m_projDatabase.empty())
      initialize();

   ossimString name ("");
   std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter = 
       m_projDatabase.find(epsg_code);
   
   if (db_iter != m_projDatabase.end())
      name = db_iter->second->name;
   
   return name;
}

getCodeFromUtmProj()

ossim_uint32 ossimEpsgProjectionDatabase::getCodeFromUtmProj

(

const ossimUtmProjection *

proj

)

const

Given UTM projection, derives the associated EPSG code. This is faster than a Db lookup.

Definition at line 950 of file ossimEpsgProjectionDatabase.cpp.

References ossimDatum::code(), ossimMapProjection::getDatum(), ossimUtmProjection::getHemisphere(), and ossimUtmProjection::getZone().

Referenced by findProjectionCode(), and ossimUtmProjection::getPcsCode().

{
   if (proj == NULL)
      return 0;

   char hemisphere = proj->getHemisphere();
   ossim_uint32 zone = proj->getZone();
   ossimString datum_code = proj->getDatum()->code();
   
   if ((zone < 1 ) || (zone > 60))
      return 0;

   ossim_uint32 epsg_code = zone;
   if ((hemisphere == 'N') && (datum_code == "WGD"))
      epsg_code += 32200;
   
   else if ((hemisphere == 'S') && (datum_code == "WGD"))
      epsg_code += 32300;

   else if ((hemisphere == 'N') && (datum_code == "WGE"))
      epsg_code += 32600;

   else if ((hemisphere == 'S') && (datum_code == "WGE"))
      epsg_code += 32700;

   else if ((hemisphere == 'N') && (datum_code == "NAS-C") && (zone > 2) && (zone < 23))
      epsg_code += 26700;

   else if ((hemisphere == 'N') && (datum_code == "NAR-C") && (zone > 2) && (zone < 24))
      epsg_code += 26900;

   else if ((hemisphere == 'N') && (datum_code == "PRP-M"))
      epsg_code += 24800;

   else if ((hemisphere == 'S') && (datum_code == "PRP-M"))
      epsg_code += 24800 + 60;
   
   else
   {
      //---
      // Use a projection code that does not imply a datum.
      // See section "6.3.3.2 Projection Codes" for definition.
      //---
      if  ( hemisphere == 'N' ) // Northern hemisphere.
      {
         epsg_code += 16000;
      }
      else // Southern hemisphere.
      {
         epsg_code += 16100;
      }
   }

   return epsg_code;
}

getProjectionsList()

void ossimEpsgProjectionDatabase::getProjectionsList

(

std::vector< ossimString > &

list

)

const

Populates caller's list with all projections currently represented in the database.

The format is: EPSG: "<Name as appears in Db>" There may be duplicate entries if the projection is duplicated in the various Db files.

Definition at line 594 of file ossimEpsgProjectionDatabase.cpp.

References ossimEpsgProjectionDatabase::ProjDbRecord::code, initialize(), m_projDatabase, ossimEpsgProjectionDatabase::ProjDbRecord::name, ossimString::push_back(), ossimString::toString(), and ossimRefPtr< T >::valid().

Referenced by ossimEpsgProjectionFactory::getTypeNameList().

{
   if (m_projDatabase.empty())
      initialize();

   std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> >::iterator db_iter = m_projDatabase.begin();
   while (db_iter != m_projDatabase.end())
   {
      ossimRefPtr<ProjDbRecord> db_record = db_iter->second;
      if ( db_record.valid() )
      {
         ossimString record ("EPSG:");
         record += ossimString::toString(db_record->code);
         record += "  \"";
         record += db_record->name;
         record += "\"";
         list.push_back(record);
      }
      ++db_iter;
   }
   return;
}

initialize()

void ossimEpsgProjectionDatabase::initialize ( ) const

protected

Populates the database with contents of DB files as specified in ossim_preferences.

Definition at line 190 of file ossimEpsgProjectionDatabase.cpp.

References ossimString::chars(), ossimString::clear(), ossimFilename::dirCat(), ossimKeywordlist::findKey(), getline(), ossimKeywordlist::getSubstringKeyList(), ossimPreferences::instance(), ossimFilename::isReadable(), ossimFilename::isRelative(), ossimPreferences::preferencesKWL(), ossimString::size(), ossimString::string(), and ossimString::trim().

Referenced by findProjection(), findProjectionCode(), findProjectionName(), and getProjectionsList().

{
   // Fetch filenames of all projection DB files from the share directory specified in
   // ossim_preferences:

   // Optional ossim share dir:
   ossimFilename share_dir = ossimPreferences::instance()->
      preferencesKWL().findKey( std::string( "ossim_share_directory" ) );
   
   ossimString regEx =  ossimString("^epsg_database_file[0-9]+");
   vector<ossimString> keys = 
      ossimPreferences::instance()->preferencesKWL().getSubstringKeyList(regEx);
   vector<ossimString>::const_iterator i = keys.begin();

   // Create only once outside the loop:
   ossimFilename epsg_path;
   ossimFilename db_name;
   ossimString group_id;
   ossimString format_id;
   ossimString line;

   // Loop over each file and read contents into memory:
   while ( i != keys.end() )
   {
      db_name.clear();
      epsg_path = ossimPreferences::instance()->preferencesKWL().findKey( (*i).string() );
      if ( epsg_path.size() )
      {
         if ( !epsg_path.isRelative() )
         {
            //---
            // example:
            // epsg_database_file0:/usr/share/ossim/projection/ossim_epsg_projections-v7_4.csv
            //---
            db_name = epsg_path;
         }
         else if ( share_dir.size() )
         {
            //---
            // example:
            // ossim_share_dir: /usr/share/ossim
            // epsg_database_file0: projection/ossim_epsg_projections-v7_4.csv
            //---
            db_name = share_dir.dirCat( epsg_path );

            //---
            // This block is for backwards compatibility.
            // Try tacking "projection" onto share dir.
            //---
            if ( !db_name.isReadable() )
            {
               db_name = share_dir.dirCat( ossimFilename("projection") );
               db_name = db_name.dirCat( epsg_path);

               // Lastly: Try tacking "ossim/projection" onto share dir.
               if ( !db_name.isReadable() )
               {
                  db_name = share_dir.dirCat( ossimFilename("ossim/projection") );
                  db_name = db_name.dirCat( epsg_path);
               }
            }
         }
      }
      ++i;
      
      if (!db_name.isReadable())
      {
         continue;
      }

      // Open the DB file:
      std::ifstream db_stream (db_name.chars());
      bool good_file = false;
      if (db_stream.good())
      {
         // Format specification implied in file's magic number:
         std::getline(db_stream, format_id.string());
         format_id.trim();
         if ((format_id == EPSG_DB_FORMAT_A) || 
             (format_id == STATE_PLANE_FORMAT_B) ||
             (format_id == SPCS_EPSG_MAP_FORMAT_C))
            good_file = true;
      }
      if (!good_file)
      {
         ossimNotify(ossimNotifyLevel_WARN)<<"ossimEpsgProjectionDatabase::initialize() -- "
            "Encountered bad database file <"<<db_name<<">. Skipping this file."<<endl;
         db_stream.close();
         continue;
      }

      // The file is good. Skip over the column descriptor line:
      std::getline(db_stream, line.string());

      // Loop to read all data records:
      while (!db_stream.eof())
      {
         ossimRefPtr<ProjDbRecord> db_record = new ProjDbRecord;
         std::getline(db_stream, line.string());
         db_record->csvRecord = line.explode(","); // ONLY CSV FILES CONSIDERED HERE
         if (db_record->csvRecord.size())
         {
            // Check if primary EPSG database format A:
            if (format_id == EPSG_DB_FORMAT_A)
            {
               db_record->code = db_record->csvRecord[A_CODE].toUInt32();
               db_record->name = db_record->csvRecord[A_NAME];
               db_record->csvFormat = FORMAT_A;
            }

            // Check if State Plane (subset of EPSG but handled differently until projection 
            // geotrans-EPSG disconnect is resolved. 
            else if (format_id == STATE_PLANE_FORMAT_B)
            {
               db_record->code = db_record->csvRecord[B_CODE].toUInt32();
               db_record->name = db_record->csvRecord[B_NAME];
               db_record->csvFormat = FORMAT_B;
            }

            // This format is for Ming-special State Plane Coordinate System coded format.
            // This format is simply a mapping from SPCS spec name (OSSIM-specific) to EPSG code.
            // Note that no proj is instantiated and no KWL is populated. Only name and EPSG mapped
            // code is saved.
            else if (format_id == SPCS_EPSG_MAP_FORMAT_C)
            {
               db_record->code = db_record->csvRecord[C_CODE].toUInt32();
               db_record->name = db_record->csvRecord[C_NAME];
               db_record->csvFormat = FORMAT_C;
            }

            m_projDatabase.insert(make_pair(db_record->code, db_record));
         }
      }

      db_stream.close();
   } // end of while loop over all DB files
}

instance()

ossimEpsgProjectionDatabase * ossimEpsgProjectionDatabase::instance ( )

static

Instantiates singleton instance of this class:

Implements singleton pattern.

Definition at line 160 of file ossimEpsgProjectionDatabase.cpp.

Referenced by ossimEpsgProjectionFactory::createProjection(), ossimWktProjectionFactory::createProjection(), ossimOgcWktTranslator::fromOssimKwl(), ossimUtmProjection::getPcsCode(), ossimMapProjection::getPcsCode(), ossimGeoTiff::getPcsUnitType(), ossimTiffWriter::getPcsUnitType(), and ossimEpsgProjectionFactory::getTypeNameList().

{
   static ossimEpsgProjectionDatabase inst;

   return &inst;
   //if(!m_instance)
   //   m_instance = new ossimEpsgProjectionDatabase;
   //return m_instance;
}

numRecords()

size_t ossimEpsgProjectionDatabase::numRecords ( ) const

inline

ENGINEERING CODE. Used for testing.

Definition at line 71 of file ossimEpsgProjectionDatabase.h.

{ return m_projDatabase.size(); }

Friends And Related Function Documentation

ossimEpsgProjectionFactory

friend class ossimEpsgProjectionFactory

friend

Definition at line 31 of file ossimEpsgProjectionDatabase.h.

Member Data Documentation

m_mutex

std::mutex ossimEpsgProjectionDatabase::m_mutex

mutableprotected

Definition at line 126 of file ossimEpsgProjectionDatabase.h.

Referenced by findProjection().

m_projDatabase

std::multimap<ossim_uint32, ossimRefPtr<ProjDbRecord> > ossimEpsgProjectionDatabase::m_projDatabase

mutableprotected

Definition at line 125 of file ossimEpsgProjectionDatabase.h.

Referenced by findProjection(), findProjectionCode(), findProjectionName(), and getProjectionsList().

---

The documentation for this class was generated from the following files:

• ossimEpsgProjectionDatabase.h
• ossimEpsgProjectionDatabase.cpp