ossimGdalProjectionFactory Class Reference

#include <ossimGdalProjectionFactory.h>

Inheritance diagram for ossimGdalProjectionFactory:

Public Member Functions
virtual ossimProjection *	createProjection (const ossimFilename &filename, ossim_uint32 entryIdx) const
	takes a filename. More...
virtual ossimProjection *	createProjection (const ossimString &name) const
virtual ossimProjection *	createProjection (const ossimKeywordlist &kwl, const char *prefix=0) const
	Take a keywordlist. More...
virtual ossimObject *	createObject (const ossimString &typeName) const
virtual ossimObject *	createObject (const ossimKeywordlist &kwl, const char *prefix=0) const
virtual void	getTypeNameList (std::vector< ossimString > &typeList) const
virtual std::list< ossimString >	getList () const
Public Member Functions inherited from ossimProjectionFactoryBase
virtual ossimProjection *	createProjection (ossimImageHandler *handler) const
Public Member Functions inherited from ossimObjectFactory
virtual	~ossimObjectFactory ()
Public Member Functions inherited from ossimObject
	ossimObject ()
virtual	~ossimObject ()
virtual ossimObject *	dup () const
virtual ossimString	getShortName () const
virtual ossimString	getLongName () const
virtual ossimString	getDescription () const
virtual ossimString	getClassName () const
virtual RTTItypeid	getType () const
virtual bool	canCastTo (ossimObject *obj) const
virtual bool	canCastTo (const RTTItypeid &id) const
virtual bool	canCastTo (const ossimString &parentClassName) const
virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
virtual std::ostream &	print (std::ostream &out) const
	Generic print method. More...
virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const
virtual void	accept (ossimVisitor &visitor)
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 ossimGdalProjectionFactory *	instance ()

Protected Member Functions
	ossimGdalProjectionFactory ()
Protected Member Functions inherited from ossimProjectionFactoryBase
ossimProjection *	createProjectionFromGeometryFile (const ossimFilename &imageFile, ossim_uint32 entryIdx) const
	This method takes the filename and tries to find an external ".geom" file. More...
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
ossimOgcWktTranslator	theWktTranslator

Static Protected Attributes
static ossimGdalProjectionFactory *	theInstance = 0

Detailed Description

Definition at line 21 of file ossimGdalProjectionFactory.h.

Constructor & Destructor Documentation

ossimGdalProjectionFactory()

ossimGdalProjectionFactory::ossimGdalProjectionFactory ( )

inlineprotected

Definition at line 63 of file ossimGdalProjectionFactory.h.

Referenced by instance().

{}

Member Function Documentation

createObject() [1/2]

ossimObject * ossimGdalProjectionFactory::createObject ( const ossimString &  typeName ) const

virtual

Creates an object given a type name.

Implements ossimObjectFactory.

Definition at line 544 of file ossimGdalProjectionFactory.cpp.

References createProjection().

{
   return createProjection(typeName);
}

createObject() [2/2]

ossimObject * ossimGdalProjectionFactory::createObject ( const ossimKeywordlist &  kwl, const char *  prefix = 0  ) const

virtual

Creates and object given a keyword list.

Implements ossimObjectFactory.

Definition at line 549 of file ossimGdalProjectionFactory.cpp.

References createProjection().

{
   return createProjection(kwl, prefix);
}

createProjection() [1/3]

ossimProjection * ossimGdalProjectionFactory::createProjection ( const ossimFilename &  filename, ossim_uint32  entryIdx  ) const

virtual

takes a filename.

This filename can be an image file or it can also be a ossim .geom file. It could be other forms of geom files as well. The factories job will be to determine what parser to use and return a projection if successful.

Implements ossimProjectionFactoryBase.

Definition at line 48 of file ossimGdalProjectionFactory.cpp.

References ossimString::c_str(), ossimString::empty(), ossimNotify(), and ossimNotifyLevel_WARN.

Referenced by createObject(), and createProjection().

{
   ossimKeywordlist kwl;
//    ossimRefPtr<ossimImageHandler> h = new ossimGdalTileSource;
   //std::cout << "ossimGdalProjectionFactory::createProjection: " << filename << std::endl;
   GDALDatasetH  h; 
   GDALDriverH   driverH = 0;
   ossimProjection* proj = 0;
   
   if(ossimString(filename).trim().empty()) return 0;
   
   h = GDALOpen(filename.c_str(), GA_ReadOnly);
   if(h)
   {
      driverH = GDALGetDatasetDriver( h );
      ossimString driverName( driverH ? GDALGetDriverShortName( driverH ) : "" );
      // use OSSIM's projection loader for NITF
      //
      if(driverName == "NITF")
      {
         GDALClose(h);
         return 0;
      }
      if(entryIdx != 0)
      {
         char** papszMetadata = GDALGetMetadata( h, "SUBDATASETS" );

         //---
         // ??? (drb) Should this be:
         // if ( entryIdx >= CSLCount(papszMetadata) ) close...
         //---
         if( papszMetadata&&(CSLCount(papszMetadata) < static_cast<ossim_int32>(entryIdx)) )
         {
            ossimNotify(ossimNotifyLevel_WARN) << "ossimGdalProjectionFactory::createProjection: We don't support multi entry handlers through the factory yet, only through the handler!";
            GDALClose(h);
            return 0;
         }
         else
         {
            GDALClose(h);
            return 0;
         }
      }

      ossimString wkt(GDALGetProjectionRef( h ));
      double geoTransform[6];
      bool transOk = GDALGetGeoTransform( h, geoTransform ) == CE_None;
      bool wktTranslatorOk = wkt.empty()?false:wktTranslator.toOssimKwl(wkt, kwl);
      if(!wktTranslatorOk)
      {
         ossim_uint32 gcpCount = GDALGetGCPCount(h);
         if(gcpCount > 3)
         {
            ossim_uint32 idx = 0;
            const GDAL_GCP* gcpList = GDALGetGCPs(h);
            ossimTieGptSet tieSet;
            if(gcpList)
            {
               for(idx = 0; idx < gcpCount; ++idx)
               {
                  ossimDpt dpt(gcpList[idx].dfGCPPixel,
                               gcpList[idx].dfGCPLine);
                  ossimGpt gpt(gcpList[idx].dfGCPY,
                               gcpList[idx].dfGCPX,
                               gcpList[idx].dfGCPZ);
                  tieSet.addTiePoint(new ossimTieGpt(gpt, dpt, .5));
               }

               //ossimPolynomProjection* tempProj = new ossimPolynomProjection;
               ossimBilinearProjection* tempProj = new ossimBilinearProjection;
               //tempProj->setupOptimizer("1 x y x2 xy y2 x3 y3 xy2 x2y z xz yz");
               tempProj->optimizeFit(tieSet);
               proj = tempProj;
            }
         }
      }
      if ( transOk && proj==0 )
      {
         ossimString proj_type(kwl.find(ossimKeywordNames::TYPE_KW));
         ossimString datum_type(kwl.find(ossimKeywordNames::DATUM_KW));
         ossimString units(kwl.find(ossimKeywordNames::UNITS_KW));
         if ( proj_type.trim().empty() &&
              (driverName == "MrSID" || driverName == "JP2MrSID") )
         {
            bool bClose = true;
            // ESH 04/2008, #54: if no rotation factors use geographic system
            if( geoTransform[2] == 0.0 && geoTransform[4] == 0.0 )
            {
               ossimString projTag( GDALGetMetadataItem( h, "IMG__PROJECTION_NAME", "" ) );
               if ( projTag.contains("Geographic") )
               {
                  bClose = false;

                  kwl.add(ossimKeywordNames::TYPE_KW,
                          "ossimEquDistCylProjection", true);
                  proj_type = kwl.find( ossimKeywordNames::TYPE_KW );

                  // Assign units if set in Metadata
                  ossimString unitTag( GDALGetMetadataItem( h, "IMG__HORIZONTAL_UNITS", "" ) );
                  if ( unitTag.contains("dd") ) // decimal degrees
                  {
                     units = "degrees";
                  }
                  else if ( unitTag.contains("dm") ) // decimal minutes
                  {
                     units = "minutes";
                  }
                  else if ( unitTag.contains("ds") ) // decimal seconds
                  {
                     units = "seconds";
                  }
               }
            }

            if ( bClose == true )
            {
               GDALClose(h);
               return 0;
            }
         }

         //---
         // Pixel type, "area", tie upper left corner of pixel, or "point", tie
         // center of pixel.
         //
         // Making an assumption that gdal normalize the geoTransform to area.
         // Geotiffs tagged with raster type of "pixel_is_point" are 
         // shifted.
         // (drb - 26 Oct. 2016)
         //---
         bool pixelTypeIsArea = true; // Assume area unless point is detected.

         /* removed replaced to fix world wrap issues with vrt bmng scene. 20160928 drb */         
#if 0

         //---
         // Look for AREA_OR_POINT:
         // This controls PIXEL_TYPE_KW value downstream which will do the half
         // pixel shift for us downstream.
         //---
         const char* areaOrPoint =  GDALGetMetadataItem(h, "AREA_OR_POINT", "");
         if (areaOrPoint)
         {
            if ( ossimString( areaOrPoint ).downcase().contains("point") )
            {
               pixelTypeIsArea = false;
            }
         }
         else
         {
            // Look for 
            const char* rasterTypeStr = GDALGetMetadataItem(
               h, "GEOTIFF_CHAR__GTRasterTypeGeoKey", "");
            if ( rasterTypeStr )
            {
               if ( ossimString( rasterTypeStr ).downcase().contains("point") )
               {
                  pixelTypeIsArea = false;
               }
            }
         }
         
         // Pixel-is-point of pixel-is area affects the location of the tiepoint since OSSIM is
         // always pixel-is-point so 1/2 pixel shift may be necessary:
         if( (driverName == "MrSID") || (driverName == "JP2MrSID") ||
             (driverName == "AIG") || (driverName == "VRT") )
         {
            const char* rasterTypeStr = GDALGetMetadataItem( h, "GEOTIFF_CHAR__GTRasterTypeGeoKey", "" );
            ossimString rasterTypeTag( rasterTypeStr );

            // If the raster type is pixel_is_area, shift the tie point by
            // half a pixel to locate it at the pixel center.
            if ( (driverName == "AIG") ||
                 (driverName == "VRT") || 
                 (rasterTypeTag.contains("RasterPixelIsArea")) )
            {
               geoTransform[0] += fabs(geoTransform[1]) / 2.0;
               geoTransform[3] -= fabs(geoTransform[5]) / 2.0;
            }
         }
         else
         {
            // Conventionally, HFA stores the pixel alignment type for each band. Here assume all
            // bands are the same. Consider only the first band:
            GDALRasterBandH bBand = GDALGetRasterBand( h, 1 );
            char** papszMetadata = GDALGetMetadata( bBand, NULL );
            if (CSLCount(papszMetadata) > 0)
            {
               for(int i = 0; papszMetadata[i] != NULL; i++ )
               {
                  ossimString metaStr = papszMetadata[i];
                  metaStr.upcase();
                  if (metaStr.contains("AREA_OR_POINT"))
                  {
                     ossimString pixel_is_point_or_area = metaStr.split("=")[1];
                     pixel_is_point_or_area.upcase();
                     if (pixel_is_point_or_area.contains("AREA"))
                     {
                        // Need to shift the tie point so that pixel is point:
                        geoTransform[0] += fabs(geoTransform[1]) / 2.0;
                        geoTransform[3] -= fabs(geoTransform[5]) / 2.0;
                     }
                     break;
                  }
               }
            }
         }
#endif

         kwl.remove(ossimKeywordNames::UNITS_KW);
         ossimDpt gsd(fabs(geoTransform[1]), fabs(geoTransform[5]));
         ossimDpt tie(geoTransform[0], geoTransform[3]);

         ossimUnitType savedUnitType =
            static_cast<ossimUnitType>(ossimUnitTypeLut::instance()->getEntryNumber(units));
         ossimUnitType unitType = savedUnitType;
         if(unitType == OSSIM_UNIT_UNKNOWN)
            unitType = OSSIM_METERS;
         
         if((proj_type == "ossimLlxyProjection") || (proj_type == "ossimEquDistCylProjection"))
         {
            ossimDpt halfGsd = gsd/2.0;

            //---
            // Sanity check the tie and scale for world scenes that can cause
            // a wrap issue.
            //---
            if ( (tie.x-halfGsd.x) < -180.0 )
            {
               tie.x = -180.0 + halfGsd.x;
               geoTransform[0] = tie.x;
            }
            if ( (tie.y+halfGsd.y) > 90.0 )
            {
               tie.y = 90.0 - halfGsd.y;
               geoTransform[3] = tie.y;
            }

            int samples = GDALGetRasterXSize(h);
            ossim_float64 degrees = samples * gsd.x;
            if ( degrees > 360.0 )
            {
               //---
               // If within one pixel of 360 degrees adjust it. Assume every
               // thing else is what they want.
               //---
               if ( fabs(degrees - 360.0) <= gsd.x )
               {
                  gsd.x = 360.0 / samples;
                  geoTransform[1] = gsd.x;
                  
                  // If we adjusted scale, fix the tie if it was on the edge.
                  if ( ossim::almostEqual( (tie.x-halfGsd.x), -180.0 ) == true )
                  {
                     tie.x = -180 + gsd.x / 2.0;
                     geoTransform[0] = tie.x;
                  }
               }
            }
            int lines = GDALGetRasterYSize(h);
            degrees = lines * gsd.y;
            if ( degrees > 180.0 )
            {
               //---
               // If within one pixel of 180 degrees adjust it. Assume every
               // thing else is what they want.
               //---
               if ( fabs(degrees - 180.0) <= gsd.y )
               {
                  gsd.y = 180.0 / lines;
                  geoTransform[5] = gsd.y;
                  
                  // If we adjusted scale, fix the tie if it was on the edge.
                  if ( ossim::almostEqual( (tie.y+halfGsd.y), 90.0 ) == true )
                  {
                     tie.y = 90.0 - gsd.y / 2.0;
                     geoTransform[3] = tie.y;
                  }
               }
            }
            
            // ESH 09/2008 -- Add the orig_lat and central_lon if the image
            // is using geographic coordsys.  This is used to convert the
            // gsd to linear units.

            // gsd could of changed above:
            halfGsd = gsd / 2.0;
            
            // Half the number of pixels in lon/lat directions
            int nPixelsLon = GDALGetRasterXSize(h)/2.0;
            int nPixelsLat = GDALGetRasterYSize(h)/2.0;

            // Shift from image corner to center in lon/lat
            double shiftLon =  nPixelsLon * fabs(gsd.x);
            double shiftLat = -nPixelsLat * fabs(gsd.y);

            // lon/lat of center pixel of the image
            double centerLon = (tie.x - halfGsd.x) + shiftLon;
            double centerLat = (tie.y + halfGsd.y) + shiftLat;

            kwl.add(ossimKeywordNames::ORIGIN_LATITUDE_KW,
                    centerLat,
                    true);
            kwl.add(ossimKeywordNames::CENTRAL_MERIDIAN_KW,
                    centerLon,
                    true);

            if(savedUnitType == OSSIM_UNIT_UNKNOWN)
            {
               unitType = OSSIM_DEGREES;
            }
         }
         
         kwl.add(ossimKeywordNames::PIXEL_SCALE_XY_KW,
                 gsd.toString(),
                 true);
         kwl.add(ossimKeywordNames::PIXEL_SCALE_UNITS_KW,
                 units,
                 true);
         kwl.add(ossimKeywordNames::PIXEL_TYPE_KW,
                 (pixelTypeIsArea?"area":"point"),
                 true);
         kwl.add(ossimKeywordNames::TIE_POINT_XY_KW,
                 tie.toString(),
                 true);
         kwl.add(ossimKeywordNames::TIE_POINT_UNITS_KW,
                 units,
                 true);

         //---
         // Image Model Transform:
         // 
         // Only output this if there is a rotation. There are two issues with
         // this:
         //
         // 1) When the image model tranform is set calls to
         // ossimMapProjection::setMetersPerPixel(gsd) have no effect.
         //
         // 2) Picking up the famous half pixel shift on output products.
         //
         // (drb 26 Oct. 2016)
         //---
         if ( (geoTransform[2] != 0.0) || (geoTransform[4] != 0.0) )
         {
            
            if ( pixelTypeIsArea == true )
            {
               // Not handled in ossimMapProjection::loadState so shift it here.
               geoTransform[0] += fabs(geoTransform[1]) / 2.0;
               geoTransform[3] -= fabs(geoTransform[5]) / 2.0;
            }
         
            std::stringstream matrixString;
            // store as a 4x4 matrix
            matrixString
               << ossimString::toString(geoTransform[1], 20)
               << " " << ossimString::toString(geoTransform[2], 20)
               << " " << 0 << " "
               << ossimString::toString(geoTransform[0], 20)
               << " " << ossimString::toString(geoTransform[4], 20)
               << " " << ossimString::toString(geoTransform[5], 20)
               << " " << 0 << " "
               << ossimString::toString(geoTransform[3], 20)
               << " " << 0 << " " << 0 << " " << 1 << " " << 0
               << " " << 0 << " " << 0 << " " << 0 << " " << 1;

            kwl.add(ossimKeywordNames::IMAGE_MODEL_TRANSFORM_MATRIX_KW,
                    matrixString.str().c_str(), true);
            kwl.add(ossimKeywordNames::IMAGE_MODEL_TRANSFORM_UNIT_KW,
                    units.string().c_str(), true);
         }

         //---
         // SPECIAL CASE:  ArcGrid in British National Grid
         //---
         if(driverName == "AIG" && datum_type == "OSGB_1936")
         {
            ossimFilename prj_file = filename.path() + "/prj.adf";

            if(prj_file.exists())
            {
               ossimKeywordlist prj_kwl(' ');
               prj_kwl.addFile(prj_file);

               ossimString proj = prj_kwl.find("Projection");

               // Reset projection and Datum correctly for BNG.
               if(proj.upcase().contains("GREATBRITAIN"))
               {

                  kwl.add(ossimKeywordNames::TYPE_KW,
                          "ossimBngProjection", true);

                  ossimString datum  = prj_kwl.find("Datum");

                  if(datum != "")
                  {
                     if(datum == "OGB_A")
                        datum = "OGB-A";
                     else if(datum == "OGB_B")
                        datum = "OGB-B";
                     else if(datum == "OGB_C")
                        datum = "OGB-C";
                     else if(datum == "OGB_D")
                        datum = "OGB-D";
                     else if(datum == "OGB_M")
                        datum = "OGB-M";
                     else if(datum == "OGB_7")
                        datum = "OGB-7";

                     kwl.add(ossimKeywordNames::DATUM_KW,
                             datum, true);
                  }
               }
            }
         }
      }
      
      if(traceDebug())
      {
         ossimNotify(ossimNotifyLevel_DEBUG) << "ossimGdalProjectionFactory: createProjection KWL = \n " << kwl << std::endl;
      }
      
      GDALClose(h);
      proj = ossimProjectionFactoryRegistry::instance()->createProjection(kwl);
   }

   return proj;
}

createProjection() [2/3]

ossimProjection * ossimGdalProjectionFactory::createProjection ( const ossimString &  name ) const

virtual

METHOD: create() Attempts to create an instance of the projection specified by name. Returns successfully constructed projection or NULL.

Implements ossimProjectionFactoryBase.

Definition at line 512 of file ossimGdalProjectionFactory.cpp.

{
   ossimString tempName = name;
   if(traceDebug())
   {
       ossimNotify(ossimNotifyLevel_WARN) << "ossimGdalProjectionFactory::createProjection: "<< name << "\n";
   }
   tempName = tempName.trim();

   if ( tempName.size() >= 6 )
   {
      ossimString testName(tempName.begin(),
                           tempName.begin()+6);
      testName = testName.upcase();
      if((testName == "PROJCS")||
         (testName == "GEOGCS"))
      {
         ossimKeywordlist kwl;
         if ( theWktTranslator.toOssimKwl(name.c_str(), kwl, "") )
         {
            if(traceDebug())
            {
               ossimNotify(ossimNotifyLevel_DEBUG)<< "ossimGdalProjectionFactory::createProjection Debug: trying to create projection \n " << kwl << std::endl;
            }
            return ossimProjectionFactoryRegistry::instance()->createProjection(kwl,"");
        }
      }
   }

   return 0;
}

createProjection() [3/3]

ossimProjection * ossimGdalProjectionFactory::createProjection ( const ossimKeywordlist &  kwl, const char *  prefix = 0  ) const

virtual

Take a keywordlist.

Implements ossimProjectionFactoryBase.

Definition at line 479 of file ossimGdalProjectionFactory.cpp.

References ossimKeywordlist::add(), createProjection(), ossimKeywordlist::find(), ossimObject::getClassName(), ossimProjection::loadState(), ossimProjection::saveState(), and ossimKeywordNames::TYPE_KW.

{
   const char *lookup = keywordList.find(prefix, ossimKeywordNames::TYPE_KW);
   if(lookup&&(!ossimString(lookup).empty()))
   {
      ossimProjection* proj = createProjection(ossimString(lookup));
      if(proj)
      {
         // make sure we restore any passed in tie points and meters per pixel information
         ossimKeywordlist tempKwl;
         ossimKeywordlist tempKwl2;
         proj->saveState(tempKwl);
         tempKwl2.add(keywordList, prefix, true);
         tempKwl.add(prefix, tempKwl2,  true);
         tempKwl.add(prefix, ossimKeywordNames::TYPE_KW, proj->getClassName(), true);
     proj->loadState(tempKwl);
         if(traceDebug())
         {
            tempKwl.clear();
            proj->saveState(tempKwl);
            if(traceDebug())
            {
               ossimNotify(ossimNotifyLevel_DEBUG)<< "ossimGdalProjectionFactory::createProjection Debug: resulting projection \n " << tempKwl << std::endl;
            }
         }

         return proj;
      }
   }
   return 0;
}

getList()

std::list< ossimString > ossimGdalProjectionFactory::getList ( ) const

virtual

METHOD: getList() Returns list of all projections represented by this factory:

Definition at line 559 of file ossimGdalProjectionFactory.cpp.

{
   std::list<ossimString> result;

   return result;
}

getTypeNameList()

void ossimGdalProjectionFactory::getTypeNameList ( std::vector< ossimString > &  typeList ) const

virtual

This should return the type name of all objects in all factories. This is the name used to construct the objects dynamially and this name must be unique.

Implements ossimObjectFactory.

Definition at line 555 of file ossimGdalProjectionFactory.cpp.

{
}

instance()

ossimGdalProjectionFactory * ossimGdalProjectionFactory::instance ( )

static

METHOD: instance() Instantiates singleton instance of this class:

Definition at line 38 of file ossimGdalProjectionFactory.cpp.

References ossimGdalProjectionFactory(), and theInstance.

Referenced by ossimSharedLibraryFinalize(), and ossimSharedLibraryInitialize().

{
   if(!theInstance)
   {
      theInstance = new ossimGdalProjectionFactory;
   }

   return (ossimGdalProjectionFactory*) theInstance;
}

Member Data Documentation

theInstance

ossimGdalProjectionFactory * ossimGdalProjectionFactory::theInstance = 0

staticprotected

Definition at line 65 of file ossimGdalProjectionFactory.h.

Referenced by instance().

theWktTranslator

ossimOgcWktTranslator ossimGdalProjectionFactory::theWktTranslator

protected

Definition at line 64 of file ossimGdalProjectionFactory.h.

---

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

• ossimGdalProjectionFactory.h
• ossimGdalProjectionFactory.cpp