#include <ossimNitfProjectionFactory.h>

Inheritance diagram for ossimNitfProjectionFactory:

Public Member Functions
virtual	~ossimNitfProjectionFactory ()
	destructor More...

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

virtual ossimProjection *	createProjection (const ossimString &name) const

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

virtual ossimProjection *	createProjection (ossimImageHandler *handler) const

virtual ossimObject *	createObject (const ossimString &typeName) const
	Creates an object given a type name. More...

virtual ossimObject *	createObject (const ossimKeywordlist &kwl, const char *prefix=0) const
	Creates and object given a keyword list. More...

virtual void	getTypeNameList (std::vector< ossimString > &typeList) const
	This should return the type name of all objects in all factories. More...

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 ossimNitfProjectionFactory *	instance ()
	METHOD: instance() More...

Private Member Functions
ossimProjection *	createProjectionFromHeaders (ossimNitfFileHeader fileHeader, ossimNitfImageHeader imageHeader) const

ossimProjection *	makeGeographic (const ossimNitfImageHeader *hdr, const ossimString &coordinateSysetm) const

ossimProjection *	makeUtm (const ossimNitfImageHeader *hdr, const ossimString &coordinateSysetm) const

ossimProjection *	makeEuiDistant (const ossimNitfImageHeader *hdr, const std::vector< ossimGpt > &gpts) const

ossimProjection *	makeBilinear (const ossimNitfImageHeader *hdr, const std::vector< ossimGpt > &gpts) const

void	parseGeographicString (const ossimString &geographicLocation, std::vector< ossimGpt > &gpts) const
	Parses geographic stream. More...

void	parseDecimalDegreesString (const ossimString &geographicLocation, std::vector< ossimGpt > &gpts) const
	Parses Decimal degree stream. More...

void	parseUtmString (const ossimString &utmLocation, ossim_uint32 &zone, std::vector< ossimDpt > &utmPoints) const

bool	parseMgrsString (const ossimString &mgrsLocationString, ossim_uint32 &zone, char &hemisphere, std::vector< ossimDpt > &utmPoints) const
	Gets UTM points, zone and hemisphere from IGEOLO field when ICORDS files = U which is UTM expressed in MGRS. More...

void	computeScaleInDecimalDegrees (const ossimNitfImageHeader *hdr, const std::vector< ossimGpt > &gpts, ossimDpt &scale) const
	Computes the scale in decimal degrees with scale.x being longitude, y being latitude. More...

void	computeScaleInMeters (const ossimNitfImageHeader *hdr, const std::vector< ossimDpt > &dpts, ossimDpt &scale) const
	Computes the scale in meters with scale.x being easting, y being northing. More...

bool	isSkewed (const std::vector< ossimGpt > &gpts) const

bool	isSkewed (const std::vector< ossimDpt > &gpts) const

bool	getBlockaPoints (const ossimNitfImageHeader *hdr, std::vector< ossimGpt > &gpts) const

bool	getGeolobPoints (const ossimNitfImageHeader *hdr, std::vector< ossimGpt > &gpts) const

	ossimNitfProjectionFactory ()
	Private constructor, users must go through instance() method. More...

	ossimNitfProjectionFactory (const ossimNitfProjectionFactory &obj)
	Private copy constructor to hide from use. More...

ossimNitfProjectionFactory &	operator= (const ossimNitfProjectionFactory &rhs)
	Private operator = to hide from use. More...

bool	isNitf (const ossimFilename &filename) const

Static Private Attributes
static ossimNitfProjectionFactory *	theInstance = 0

Additional Inherited Members
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 ()

Detailed Description

Definition at line 27 of file ossimNitfProjectionFactory.h.

Constructor & Destructor Documentation

◆ ~ossimNitfProjectionFactory()

ossimNitfProjectionFactory::~ossimNitfProjectionFactory ( )

virtual

destructor

Definition at line 41 of file ossimNitfProjectionFactory.cpp.

42 {

43 }

◆ ossimNitfProjectionFactory() [1/2]

ossimNitfProjectionFactory::ossimNitfProjectionFactory ( )

private

Private constructor, users must go through instance() method.

Definition at line 37 of file ossimNitfProjectionFactory.cpp.

Referenced by instance().

38 {

39 }

◆ ossimNitfProjectionFactory() [2/2]

ossimNitfProjectionFactory::ossimNitfProjectionFactory ( const ossimNitfProjectionFactory & obj )

private

Private copy constructor to hide from use.

Definition at line 1113 of file ossimNitfProjectionFactory.cpp.

1114 {

1115 }

Member Function Documentation

◆ computeScaleInDecimalDegrees()

void ossimNitfProjectionFactory::computeScaleInDecimalDegrees	(	const ossimNitfImageHeader *	hdr,
		const std::vector< ossimGpt > &	gpts,
		ossimDpt &	scale
	)		const

private

Computes the scale in decimal degrees with scale.x being longitude, y being latitude.

Parameters

hdr	The nitf image header from the currently opened nitf file.
gpts	Four geographic corner points.
scale	Scale to initialize.

Note: scale will be set to nan's on failure.

Definition at line 828 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getImageRect(), ossimIrect::height(), isSkewed(), ossimDpt::makeNan(), ossimIrect::width(), ossimDpt::x, and ossimDpt::y.

Referenced by makeEuiDistant().

 {
    if ( !hdr || isSkewed(gpts))
    {
       scale.makeNan();
       return;
    }
    ossimIrect imageRect = hdr->getImageRect();
 
    //---
    // Calculate the scale.  This assumes that the corner points are for the
    // edge of the corner pixels, not the center of the corner pixels.
    //---
    double longitudeSize  = 0.0;
    double latitudeSize = 0.0;
    if ( (gpts[1].lond() < 0.0) && (gpts[0].lond() >= 0) )
    {
       //---
       // Upper right negative(Western), upper left positive (Eastern).
       // Crossing date line maybe???
       //---
       longitudeSize = (gpts[1].lond() + 360.0) - gpts[0].lond();
    }
    else
    {
       longitudeSize = gpts[1].lond() - gpts[0].lond();
    }
 
    latitudeSize = gpts[0].latd() - gpts[2].latd();
 
    double rows = imageRect.height();
    double cols = imageRect.width();
 //   double rows = hdr->getNumberOfRows();
 //   double cols = hdr->getNumberOfCols();
        
    if (!rows || !cols)
    {
       scale.makeNan();
       return;
    }
    scale.y = latitudeSize  / rows;
    scale.x = longitudeSize / cols;
 }

◆ computeScaleInMeters()

void ossimNitfProjectionFactory::computeScaleInMeters	(	const ossimNitfImageHeader *	hdr,
		const std::vector< ossimDpt > &	dpts,
		ossimDpt &	scale
	)		const

private

Computes the scale in meters with scale.x being easting, y being northing.

Parameters

hdr	The nitf image header from the currently opened nitf file.
dpts	Four utm corner points.
scale	Scale to initialize.

Note: scale will be set to nan's on failure.

Definition at line 875 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getImageRect(), ossimIrect::height(), isSkewed(), ossimDpt::makeNan(), ossimIrect::width(), ossimDpt::x, x, ossimDpt::y, and y.

 {
    if ( !hdr || isSkewed(dpts))
    {
       scale.makeNan();
       return;
    }
    ossimIrect imageRect = hdr->getImageRect();
 
    //---
    // Calculate the scale.  This assumes that the corner points are for the
    // edge of the corner pixels, not the center of the corner pixels.
    //---
    double eastingSize  = 0.0;
    double northingSize = 0.0;
    
    eastingSize  = fabs(dpts[1].x - dpts[0].x);
    northingSize = fabs(dpts[0].y - dpts[3].y);
 
    double rows = imageRect.height();//hdr->getNumberOfRows();
    double cols = imageRect.width();//hdr->getNumberOfCols();
 
    if (!rows || !cols)
    {
       scale.makeNan();
       return;
    }
    scale.y = northingSize / rows;
    scale.x = eastingSize  / cols;
 }

◆ createObject() [1/2]

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

virtual

Creates an object given a type name.

Returns: ossimObject on success or NULL if not handled.

Implements ossimObjectFactory.

Definition at line 130 of file ossimNitfProjectionFactory.cpp.

References createProjection().

 {
    return (ossimObject*)createProjection(typeName);
 }

◆ createObject() [2/2]

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

virtual

Creates and object given a keyword list.

Returns: ossimObject on success or NULL if not handled.

Implements ossimObjectFactory.

Definition at line 136 of file ossimNitfProjectionFactory.cpp.

References createProjection().

 {
    return createProjection(kwl, prefix);
 }

◆ createProjection() [1/4]

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

virtual

Parameters

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.

Returns: ossimProjection* on success or NULL if not handled.

Implements ossimProjectionFactoryBase.

Definition at line 56 of file ossimNitfProjectionFactory.cpp.

Referenced by createObject(), and createProjection().

 {
    static const char MODULE[]="ossimNitfProjectionFactory::createProjection";
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << MODULE << " DEBUG:"
          << "\nfilename:  " << filename
          << "\nentryIdx:  " << entryIdx
          << std::endl;
    }
 
    // See if there is an external geomtry.
    ossimProjection* result = createProjectionFromGeometryFile(filename, entryIdx);
    if (result)
    {
       return result;
    }
 
    if(!isNitf(filename))
    {
       return result; // result is NULL
    }
 
    ossimRefPtr<ossimNitfFile> nitf = new ossimNitfFile();
    if (!nitf->parseFile(filename))
    {
       return result; // result is NULL
    }
 
    long imageIndex = static_cast<long>(entryIdx);
    if ( imageIndex > (nitf->getHeader()->getNumberOfImages()-1) )
    {
       return result;
    }
 
    ossimRefPtr<ossimNitfImageHeader> imageHeader = nitf->getNewImageHeader(imageIndex);
    {
       if (!imageHeader)
       {
          return result;
       }
    }
 
    result = createProjectionFromHeaders(nitf->getHeader(),
                                         imageHeader.get());
    if (traceDebug())
    {
       ossimString coordinateSystem   = imageHeader->getCoordinateSystem();
       ossimNotify(ossimNotifyLevel_DEBUG)
          << MODULE << " DEBUG:"
          << "\ncoordinateSysetm:       " << coordinateSystem
          << std::endl;
    }
 
    return result;
 }

◆ createProjection() [2/4]

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

virtual

Parameters

name	This would be a projection name.

Returns: ossimProjection* on success or NULL if not handled.

Implements ossimProjectionFactoryBase.

Definition at line 124 of file ossimNitfProjectionFactory.cpp.

 {
    return 0;
 }

◆ createProjection() [3/4]

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

virtual

Parameters

kwl	Keyword list to instantiate projection from.

Returns: ossimProjection* on success or NULL if not handled.

Implements ossimProjectionFactoryBase.

Definition at line 117 of file ossimNitfProjectionFactory.cpp.

 {
    return 0;
 }

◆ createProjection() [4/4]

ossimProjection * ossimNitfProjectionFactory::createProjection ( ossimImageHandler * handler ) const

virtual

Reimplemented from ossimProjectionFactoryBase.

Definition at line 147 of file ossimNitfProjectionFactory.cpp.

References createProjection(), createProjectionFromHeaders(), ossimImageHandler::getCurrentEntry(), ossimNitfTileSource::getCurrentImageHeader(), ossimNitfTileSource::getFileHeader(), ossimImageHandler::getFilename(), and isNitf().

 {
    ossimNitfTileSource* nitfTileSource = dynamic_cast<ossimNitfTileSource*> (handler);
    ossimProjection* result = 0;
    if(nitfTileSource)
    {
       if(!result)
       {
          ossimNitfImageHeader* imageHeader = nitfTileSource->getCurrentImageHeader();
          if(imageHeader)
          {
             result = createProjectionFromHeaders(nitfTileSource->getFileHeader(),imageHeader);
          }
       }
    }
    else if(isNitf(handler->getFilename()))
    {
       result =  createProjection(handler->getFilename(), handler->getCurrentEntry());
    }
    return result;
 }

◆ createProjectionFromHeaders()

ossimProjection * ossimNitfProjectionFactory::createProjectionFromHeaders	(	ossimNitfFileHeader *	fileHeader,
		ossimNitfImageHeader *	imageHeader
	)		const

private

Definition at line 187 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getCoordinateSystem(), ossimNitfFileHeader::getVersion(), makeGeographic(), and makeUtm().

Referenced by createProjection().

 {
    ossimProjection* result = 0;
    ossimString version = fileHeader->getVersion();
    ossimString coordinateSystem = imageHeader->getCoordinateSystem();
 
    // Note in version 2.0 ICORDS with 'N' == NONE.  In 2.1 it is UTM North:
    if (coordinateSystem == "G" || coordinateSystem == "D")
    {
       result = makeGeographic(imageHeader, coordinateSystem);
    }
    else if( (coordinateSystem == "S") || (coordinateSystem == "U") || 
             ( version != "02.00" && (coordinateSystem == "N") ) )
    {
       result = makeUtm(imageHeader, coordinateSystem);
    }
    return result;
 }

◆ getBlockaPoints()

bool ossimNitfProjectionFactory::getBlockaPoints	(	const ossimNitfImageHeader *	hdr,
		std::vector< ossimGpt > &	gpts
	)		const

private

Parameters

hdr	The nitf image header from the currently opened nitf file.
gpts	Ground points to initialize from BLOCKA tag. This should be an empty vector.

Returns: true if BLOCKA tag was parsed.

Definition at line 670 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getTagData().

Referenced by makeGeographic().

 {
    if (!hdr)
    {
       return false;
    }
 
    ossimRefPtr<ossimNitfRegisteredTag> tag =
       hdr->getTagData(ossimString("BLOCKA"));
 
    if (!tag)
    {
       if (traceDebug())
       {
          ossimNotify(ossimNotifyLevel_DEBUG)
             << "ossimNitfProjectionFactory::getBlockaPoints DEBUG:"
             << "\nBLOCKA tag not found."
             << std::endl;
       }
       return false;
    }
 
    if (gpts.size())
    {
       gpts.clear();
    }
 
    ossimNitfBlockaTag* blockaTag = PTR_CAST(ossimNitfBlockaTag, tag.get());
    if (!blockaTag)
    {
       return false;
    }
 
    ossimDpt dpt;
    ossimGpt gpt;
 
    // Get the upper left or first row first column.
    blockaTag->getFrfcLoc(dpt);
    gpt.latd(dpt.y);
    gpt.lond(dpt.x);
    gpts.push_back(gpt);
 
    // Get the upper right or first row last column.
    blockaTag->getFrlcLoc(dpt);
    gpt.latd(dpt.y);
    gpt.lond(dpt.x);
    gpts.push_back(gpt);
 
    // Get the lower right or last row last column.
    blockaTag->getLrlcLoc(dpt);
    gpt.latd(dpt.y);
    gpt.lond(dpt.x);
    gpts.push_back(gpt);
 
    // Get the lower left or last row first column.
    blockaTag->getLrfcLoc(dpt);
    gpt.latd(dpt.y);
    gpt.lond(dpt.x);
    gpts.push_back(gpt);
 
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << "ossimNitfProjectionFactory::getBlockaPoints DEBUG:"
          << std::endl;
       for (int i=0; i<4; ++i)
       {
          ossimNotify(ossimNotifyLevel_DEBUG)
             << "gpt[" << i << "] " << gpts[i] << std::endl;
       }
    }
    
    return true;
 }

◆ getGeolobPoints()

bool ossimNitfProjectionFactory::getGeolobPoints	(	const ossimNitfImageHeader *	hdr,
		std::vector< ossimGpt > &	gpts
	)		const

private

Parameters

hdr	The nitf image header from the currently opened nitf file.
gpts	Ground points to initialize from GEOLOB tag. This should be an empty vector.

Returns: true if GEOLOB tag was parsed.

Definition at line 747 of file ossimNitfProjectionFactory.cpp.

Referenced by makeGeographic().

 {
    static const char MODULE[] = "ossimNitfProjectionFactory::getGeolobPoints";
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << MODULE << " entered...\n";
    }
    bool result = false;
 
    gpts.clear();
    
    if ( hdr )
    {
       ossimRefPtr<ossimNitfRegisteredTag> tag = hdr->getTagData(ossimString("GEOLOB"));
       if ( tag.valid() )
       {
          ossimNitfGeolobTag* geolob = dynamic_cast<ossimNitfGeolobTag*>(tag.get());
          if ( geolob )
          {
             ossimDpt gsd( geolob->getDegreesPerPixelLon(), geolob->getDegreesPerPixelLat() );
             if ( (gsd.x > 0.0 ) && (gsd.y > 0.0) )
             {
                ossim_int32 rows = hdr->getNumberOfRows();
                ossim_int32 cols = hdr->getNumberOfCols();
                if ( (rows > 0) && (cols > 0) )
                {
                   ossimGpt gpt(0.0, 0.0, 0.0);
                   ossimDpt origin( geolob->getLso(), geolob->getPso() );
 
                   if (traceDebug())
                   {
                      ossimNotify(ossimNotifyLevel_DEBUG)
                         << "origin: " << origin
                         << "\ngsd:    " << gsd << "\n";
                   }
 
                   // Note: Edge to edge here as origin is shifted makeEuiDistant method.
                   // UL:
                   gpt.latd(origin.y);
                   gpt.lond(origin.x);
                   gpts.push_back(gpt);
 
                   // UR:
                   gpt.latd(origin.y);
                   gpt.lond(origin.x + cols*gsd.x);
                   gpts.push_back(gpt);
                                
                   // LR:
                   gpt.latd(origin.y - rows*gsd.y);
                   gpt.lond(origin.x + cols*gsd.x);
                   gpts.push_back(gpt);
 
                   // LL:
                   gpt.latd(origin.y - rows*gsd.y);
                   gpt.lond(origin.x);
                   gpts.push_back(gpt);
                  
                   result = true;
                }
             }
          }
       }
       else if (traceDebug())
       {
          ossimNotify(ossimNotifyLevel_DEBUG)
             << "No GEOLOB tag found.\n";
       }
    }
 
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << MODULE << " exit status: " << (result?"true":"false") << "\n";
    }
    
    return result;
 }

◆ getTypeNameList()

void ossimNitfProjectionFactory::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.

Parameters

typeList List to append names to.

Implements ossimObjectFactory.

Definition at line 142 of file ossimNitfProjectionFactory.cpp.

 {
    
 }

◆ instance()

ossimNitfProjectionFactory * ossimNitfProjectionFactory::instance ( )

static

METHOD: instance()

Returns: Singleton instance of this class.

Definition at line 45 of file ossimNitfProjectionFactory.cpp.

References ossimNitfProjectionFactory(), and theInstance.

Referenced by ossimProjectionFactoryRegistry::initializeDefaults().

 {
    if(!theInstance)
    {
       theInstance = new ossimNitfProjectionFactory();
    }
 
    return theInstance;
 }

◆ isNitf()

bool ossimNitfProjectionFactory::isNitf ( const ossimFilename & filename ) const

private

Definition at line 169 of file ossimNitfProjectionFactory.cpp.

References ossimString::c_str().

Referenced by createProjection().

 {
    std::ifstream in(filename.c_str(), ios::in|ios::binary);
    
    if(in)
    {
       char nitfFile[4];
       in.read((char*)nitfFile, 4);
       ossimString s(nitfFile, nitfFile+4);
       if ( (s == "NITF") || (s == "NSIF") )
       {
          return true;
       }
    }
 
    return false;
 }

◆ isSkewed() [1/2]

bool ossimNitfProjectionFactory::isSkewed ( const std::vector< ossimGpt > & gpts ) const

private

Returns: true if corner points are parallel to each other, else false.

Definition at line 649 of file ossimNitfProjectionFactory.cpp.

Referenced by computeScaleInDecimalDegrees(), and computeScaleInMeters().

 {
 
    return !( (gpts[0].latd() == gpts[1].latd()) &&
              (gpts[2].latd() == gpts[3].latd()) &&
              (gpts[0].lond() == gpts[3].lond()) &&
              (gpts[1].lond() == gpts[2].lond()) );
    
 }

◆ isSkewed() [2/2]

bool ossimNitfProjectionFactory::isSkewed ( const std::vector< ossimDpt > & gpts ) const

private

Returns: true if corner points are parallel to each other, else false.

Definition at line 660 of file ossimNitfProjectionFactory.cpp.

References x, and y.

 {
    return !( (dpts[0].y == dpts[1].y) &&
              (dpts[2].y == dpts[3].y) &&
              (dpts[0].x == dpts[3].x) &&
              (dpts[1].x == dpts[2].x) );
    
 }

◆ makeBilinear()

ossimProjection * ossimNitfProjectionFactory::makeBilinear	(	const ossimNitfImageHeader *	hdr,
		const std::vector< ossimGpt > &	gpts
	)		const

private

Returns: ossimBilinearProjection* or null if a problem is found.

Definition at line 617 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getNumberOfCols(), ossimNitfImageHeader::getNumberOfRows(), and ossimRefPtr< T >::release().

 {
    double rows = hdr->getNumberOfRows();
    double cols = hdr->getNumberOfCols();
    
    ossimDpt ul(0.0, 0.0);
    ossimDpt ur(cols-1.0, 0.0);   
    ossimDpt lr(cols-1.0, rows-1.0);
    ossimDpt ll(0.0, rows-1.0);
 
    ossimRefPtr<ossimBilinearProjection> proj;
    try
    {
      proj = new ossimBilinearProjection(ul,
                                       ur,
                                       lr,
                                       ll,
                                       gpts[0],
                                       gpts[1],
                                       gpts[2],
                                       gpts[3]);
    }
    catch(...)
    {
       proj = 0;
    }
    
    return proj.release();
 }

◆ makeEuiDistant()

ossimProjection * ossimNitfProjectionFactory::makeEuiDistant	(	const ossimNitfImageHeader *	hdr,
		const std::vector< ossimGpt > &	gpts
	)		const

private

Parameters

hdr	The nitf image header.
gpts	Four geographic corner points.

Returns: ossimEquDistCylProjection* or null if a problem is found.

Definition at line 570 of file ossimNitfProjectionFactory.cpp.

References ossim::acosd(), computeScaleInDecimalDegrees(), ossimMapProjection::getOrigin(), ossimDpt::hasNans(), ossimGpt::lat, ossimGpt::latd(), ossimGpt::lond(), ossimMapProjection::setDecimalDegreesPerPixel(), ossimEquDistCylProjection::setOrigin(), ossimMapProjection::setUlTiePoints(), ossimDpt::x, x, and ossimDpt::y.

 {
    ossimEquDistCylProjection* proj = 0;
 
    // Get the scale.
    ossimDpt scale;
    computeScaleInDecimalDegrees(hdr, gpts, scale);
 
    if (scale.hasNans())
    {
       return proj;
    }
    // Make the projection.
    proj = new ossimEquDistCylProjection();
 
    if ( scale.x )
    {
       ossimGpt origin = proj->getOrigin();
 
       //---
       // ossimEquDistCylProjection uses the origin_latitude for meters per pixel (gsd)
       // computation.  So is not set in tiff tags, compute to achieve the proper
       // horizontal scaling.
       //---
       ossim_float64 x = (scale.y/scale.x);
       if ( x > 1.0 ) x = 1.0; // Limit to domain of arc cosine function.
       origin.lat = ossim::acosd(x);
 
       proj->setOrigin(origin);
    }
 
    // Shift the tie to the center of the pixel.
    ossimGpt tiePoint;
    tiePoint.latd(gpts[0].latd() - (scale.y/2.0));
    tiePoint.lond(gpts[0].lond() + (scale.x/2.0));
 
    // Set the tie points.
    proj->setUlTiePoints(tiePoint);
 
    // Set the scale.
    proj->setDecimalDegreesPerPixel(scale);
 
    return proj;
 }

◆ makeGeographic()

ossimProjection * ossimNitfProjectionFactory::makeGeographic	(	const ossimNitfImageHeader *	hdr,
		const ossimString &	coordinateSysetm
	)		const

private

Parameters

hdr	The nitf image header.
coordianteSystem	The coordinate system as a string.
geographicLocation	This should contain the four corner strings.

Returns: Either an ossimEquDistCylProjection or ossimBilinearProjection depending upon if the corner points line up evenly.

Definition at line 207 of file ossimNitfProjectionFactory.cpp.

References getBlockaPoints(), ossimNitfImageHeader::getGeographicLocation(), getGeolobPoints(), and ossimString::size().

Referenced by createProjectionFromHeaders().

 {
    ossimProjection* proj = 0;
 
    if (hdr)
    {
       // To hold corner points.
       std::vector<ossimGpt> gpts(0);
       std::vector<ossimGpt> geolobPts(0);
       std::vector<ossimGpt> blockaPts(0);
       std::vector<ossimGpt> igeoloPts(0);      
 
       //---
       // Look for points from the GEOLOB tag if geographic coordinate
       // system. This has an origin and scale with good precision.
       //
       // NOTES:
       // 1) More precise than the BLOCKA
       // 2) Making an assumption this is not present if rotated.
       //    Might need to test for rotation in BLOCKA and IGEOLO
       // 3) GEOLOB points are always rectangular as they are derived
       //    from a single origin.
       //---
       getGeolobPoints(hdr, geolobPts);
 
       //---
       // Look for points from the BLOCKA tag.  This may or may not be present.
       // If present since it has six digit precision use it for the points.
       //---
       getBlockaPoints(hdr, blockaPts);
             
       if ( blockaPts.empty() )
       {
          // Least precise IGEOLO field:
          ossimString geographicLocation = hdr->getGeographicLocation();
          if ( geographicLocation.size() )
          {
             if (traceDebug())
             {
                ossimNotify(ossimNotifyLevel_DEBUG)
                   << "ossimNitfProjectionFactory::makeGeographic DEBUG:"
                   << "\ngeographicLocation: " << geographicLocation
                   << std::endl;
             }
             
             if (coordinateSysetm == "G")
             {
                //---
                // If coord system is G then format is:
                // Lat = ddmmssX
                //       where d is degrees and m is minutes
                //       and s is seconds and X is either N (North) or S (South).
                // Lon = dddmmssX
                //       where d is degrees and m is minutes
                //       and s is seconds and X is either N (North) or S (South).
                //---
                parseGeographicString(geographicLocation, igeoloPts);
             }
             else if (coordinateSysetm == "D")
             {
                //---
                // If coor system is D then format is:
                // +-dd.ddd +-dd.ddd four times where + is northern hemispher and
                // - is souther hemisphere for lat and longitude
                // + is easting and - is westing.
                //---
                parseDecimalDegreesString(geographicLocation, igeoloPts);
             }
             
          } // matches: if ( geographicLocation.size() )
          
       } // matches: if ( getBlockaPoints(hdr, gpts) == false )
 
       bool isSkewedFlag = false;      
       if ( blockaPts.size() )
       {
          isSkewedFlag = isSkewed(blockaPts);
          if ( (isSkewedFlag == false) && geolobPts.size() )
          {
             gpts = geolobPts; // Not skewed and have more accurate geolob points.
          }
          else
          {
             gpts = blockaPts;
          }
       }
       else if ( igeoloPts.size() )
       {
          isSkewedFlag = isSkewed(igeoloPts);
          if ( (isSkewedFlag == false) && geolobPts.size() )
          {
             gpts = geolobPts; // Not skewed and have more accurate geolob points.
          }
          else
          {
             gpts = igeoloPts;
          }
       }
          
       if (gpts.size() == 4)
       {   
          ossimDpt scaleTest;
          computeScaleInDecimalDegrees(hdr, gpts, scaleTest);
 
          //---
          // Getting very small value above 1.0 when dividing; hence the .000001.
          // Code calling arc cosine function in makeEquDistant(...) clamps to 1.0.
          //---
          if ( !isSkewedFlag && (ossim::abs(scaleTest.y/scaleTest.x) <= 1.000001))
          {
             proj = makeEuiDistant(hdr, gpts);
          }
          else
          {
             // Image is rotated.  Make a Bilinear.
             proj = makeBilinear(hdr, gpts);
          }
       }
 
       if (traceDebug() && proj)
       {
          ossimNotify(ossimNotifyLevel_DEBUG)
             << "ossimNitfProjectionFactory::makeGeographic DEBUG:"
             << "\nUpper left corner:   " << gpts[0]
             << "\nUpper right corner:  " << gpts[1]
             << "\nLower right corner:  " << gpts[2]
             << "\nLower left corner:   " << gpts[3] << "\n";
          proj->print( ossimNotify(ossimNotifyLevel_DEBUG));
       }
       
    } // matches: if (hdr)
 
    return proj;
 }

◆ makeUtm()

ossimProjection * ossimNitfProjectionFactory::makeUtm	(	const ossimNitfImageHeader *	hdr,
		const ossimString &	coordinateSysetm
	)		const

private

Parameters

hdr	The nitf image header.
coordianteSystem	The coordinate system as a string.
geographicLocation	This should contain the four corner strings.

Returns: ossimUtmProjection

Definition at line 344 of file ossimNitfProjectionFactory.cpp.

References ossimNitfImageHeader::getGeographicLocation(), parseMgrsString(), parseUtmString(), ossimString::size(), and status.

Referenced by createProjectionFromHeaders().

 {
    ossimProjection* proj = 0;
    bool isBilinear = false;
    if (hdr)
    {
       ossimString geographicLocation = hdr->getGeographicLocation();
 
       if ( geographicLocation.size() )
       {
          std::vector<ossimDpt> utmPoints;
          ossim_uint32 zone;
          ossimDpt scale;
          char hemisphere = 'N';
          
          bool status = true;
          if ( coordinateSystem == "U")
          {
             // Sets zone, hemisphere and utmPoints. Returns true on success.
             status = parseMgrsString(geographicLocation, zone, hemisphere, utmPoints);
          }
          else
          {
             // Sets zone and utmPoints.  Void return...
             parseUtmString(geographicLocation, zone, utmPoints);
             if(coordinateSystem == "S")
             {
                hemisphere = 'S';
             }
          }
          
          if ( status )
          {
             if ( traceDebug() )
             {
                std::string s;
                s.push_back(hemisphere);
                ossimNotify(ossimNotifyLevel_DEBUG)
                   << "ossimNitfProjectionFactory::makeUtm DEBUG"
                   << "\ngeo string: " << geographicLocation
                   << "\nutm zone:   " << zone
                   << "\nhemisphere: " << hemisphere
                   << std::endl;
                for (ossim_uint32 i=0; i<utmPoints.size(); ++i)
                {
                   ossimNotify(ossimNotifyLevel_DEBUG)
                      << "utmPoints[" << utmPoints[i] << std::endl;
                }
             }
             
             ossimRefPtr<ossimUtmProjection> uproj = new ossimUtmProjection;
             uproj->setHemisphere(hemisphere);
             uproj->setZone(zone);
             
             if(isSkewed(utmPoints))
             {
                std::vector<ossimGpt> gpts;
                
                // Try blocka points first as they are more accurate.
                if ( getBlockaPoints(hdr, gpts) == false )
                {
                   ossimGpt ul = uproj->inverse(utmPoints[0]);
                   ossimGpt ur = uproj->inverse(utmPoints[1]);
                   ossimGpt lr = uproj->inverse(utmPoints[2]);
                   ossimGpt ll = uproj->inverse(utmPoints[3]);
                   gpts.push_back(ul);
                   gpts.push_back(ur);
                   gpts.push_back(lr);
                   gpts.push_back(ll);
                }
 
                //---
                // Make a bilinear either from our skewed utm points or the points from the blocka
                // tag.
                //---
                proj = makeBilinear(hdr, gpts);
                if(proj) scale = proj->getMetersPerPixel();
                isBilinear = true;
                uproj = 0; // Done with utm projeciton
                
             }
             else
             {
                computeScaleInMeters(hdr, utmPoints, scale);
                
                //---
                // Assign our projection to the return "proj".
                // Use ossimRefPtr::release the so we don't delete proj when uproj
                // goes out of scope.
                //---
                proj = uproj.release(); 
             }
             
             if( (scale.hasNans() == false)&&(!isBilinear) )
             {
                //---
                // Get the tie point.
                // 
                // Look for the the BLOCKA tag which may or may not be present.
                // This has six digit precision in decimal degrees which equates to
                // about 0.11 meters (at equator) as compared to 1.0 accuaracy of the
                // IGEOLO field.
                //---
                ossimDpt tie;
                std::vector<ossimGpt> gpts;
                if ( getBlockaPoints(hdr, gpts) )
                {
                   if (traceDebug())
                   {
                      ossimNotify(ossimNotifyLevel_DEBUG)
                         << "ossimNitfProjectionFactory::makeUtm DEBUG:"
                         << "\nTie point from blocka: " << gpts[0]
                         << endl;
                   }
                   
                   tie = proj->forward(gpts[0]);
                   tie.x += scale.x/2.0;
                   tie.y -= scale.y/2.0;
                }
                else
                {
                   tie.x = utmPoints[0].x + scale.x/2.0;
                   tie.y = utmPoints[0].y - scale.y/2.0;
                }
                
                if (traceDebug())
                {
                   ossimNotify(ossimNotifyLevel_DEBUG)
                      << "ossimNitfProjectionFactory::makeUtm DEBUG:"
                      << "\nTie point: " << tie
                      << "\nScale:     " << scale
                      << endl;
                }
                
                // Set the tie and scale.
                ossimMapProjection* mproj = dynamic_cast<ossimMapProjection*>(proj);
                if ( mproj )
                {
                   mproj->setUlEastingNorthing(tie);
                   mproj->setMetersPerPixel(scale);
                }
                else // cannot cast
                {
                   if ( proj )
                   {
                      delete proj;
                      proj = 0;
                   }
                }
             }
             else if(isBilinear)
             {
                // blank
             }
             else // Scale has nans
             {
                if ( proj )
                {
                   delete proj;
                   proj = 0;
                }
             }
             
          } // matches: if (status)
 
       } // matches:  if ( geographicLocation.size() )
          
    } // matches: if (hdr)
    
    return proj;
 }

◆ operator=()

ossimNitfProjectionFactory & ossimNitfProjectionFactory::operator= ( const ossimNitfProjectionFactory & rhs )

private

Private operator = to hide from use.

Definition at line 1117 of file ossimNitfProjectionFactory.cpp.

 {
    return *this;
 }

◆ parseDecimalDegreesString()

void ossimNitfProjectionFactory::parseDecimalDegreesString	(	const ossimString &	geographicLocation,
		std::vector< ossimGpt > &	gpts
	)		const

private

Parses Decimal degree stream.

This will initialize the vector of ossimGpt's with the four corners. If the string geographicLocation is not the correct size the vector will be cleared.

Parameters

geographicLocation	String containing corners.
gpts	Vector to initialize.

Definition at line 1066 of file ossimNitfProjectionFactory.cpp.

References ossimString::c_str(), and ossimString::toDouble().

 {
    const char* bufPtr = geographicLocation.c_str();
 
    
    ossimString ulLat(bufPtr,
                      bufPtr + 7);
    bufPtr+=7;
    ossimString ulLon(bufPtr,
                      bufPtr+8);
    bufPtr+=8;
    ossimString urLat(bufPtr,
                      bufPtr + 7);
    bufPtr+=7;
    ossimString urLon(bufPtr,
                      bufPtr+8);
    bufPtr+=8;
    ossimString lrLat(bufPtr,
                      bufPtr + 7);
    bufPtr+=7;
    ossimString lrLon(bufPtr,
                      bufPtr+8);
    bufPtr+=8;
    ossimString llLat(bufPtr,
                      bufPtr + 7);
    bufPtr+=7;
    ossimString llLon(bufPtr,
                      bufPtr+8);
 
    gpts.push_back(ossimGpt(ulLat.toDouble(), ulLon.toDouble()));
    gpts.push_back(ossimGpt(urLat.toDouble(), urLon.toDouble()));
    gpts.push_back(ossimGpt(lrLat.toDouble(), lrLon.toDouble()));
    gpts.push_back(ossimGpt(llLat.toDouble(), llLon.toDouble()));
    
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << "ossimNitfProjectionFactory::parseDecimalDegreesString DEBUG:"
          << "\nground point[" << 0 << "]:  " << gpts[0]
          << "\nground point[" << 1 << "]:  " << gpts[1]
          << "\nground point[" << 2 << "]:  " << gpts[2]
          << "\nground point[" << 3 << "]:  " << gpts[3]
          << std::endl;
    }
 }

◆ parseGeographicString()

void ossimNitfProjectionFactory::parseGeographicString	(	const ossimString &	geographicLocation,
		std::vector< ossimGpt > &	gpts
	)		const

private

Parses geographic stream.

This will initialize the vector of ossimGpt's with the four corners. If the string geographicLocation is not the correct size the vector will be cleared.

Parameters

geographicLocation	String containing corners.
gpts	Vector to initialize.

Definition at line 972 of file ossimNitfProjectionFactory.cpp.

References ossimString::c_str(), ossimDms::getDegrees(), ossimGpt::latd(), ossimGpt::lond(), ossimDms::setDegrees(), ossimDms::setLatFlag(), ossimString::size(), and SPACE.

 {
    gpts.clear();
    
    if (geographicLocation.size() != 60)
    {
       return;
    }
 
    std::string::size_type geo_index = 0;
    for (int i=0; i<4; ++i)
    {
       //---
       // We have to split up the geographicLocation string for the dms class.
       // 
       // geographicLocation = ddmmssXdddmmssX (four times).
       // "dd mm ss X" has a string length of 10
       // "ddd mm ss X" has a string length of 11
       //---
       std::string::size_type lat_index = 0;
       std::string::size_type lon_index = 0;
       const char SPACE = ' ';
       ossimString latString(10, SPACE);
       ossimString lonString(11, SPACE);
 
       // degrees latitude
       latString[lat_index++] = geographicLocation[geo_index++];
       latString[lat_index++] = geographicLocation[geo_index++];
       ++lat_index;
 
       // minutes latitude
       latString[lat_index++] = geographicLocation[geo_index++];
       latString[lat_index++] = geographicLocation[geo_index++];
       ++lat_index;
       
       // seconds latitude
       latString[lat_index++] = geographicLocation[geo_index++];
       latString[lat_index++] = geographicLocation[geo_index++];
       ++lat_index;
 
       // hemisphere
       latString[lat_index++] = geographicLocation[geo_index++];
  
       // degrees longitude
       lonString[lon_index++] = geographicLocation[geo_index++];
       lonString[lon_index++] = geographicLocation[geo_index++];
       lonString[lon_index++] = geographicLocation[geo_index++];
       ++lon_index;
 
       // minutes longitude
       lonString[lon_index++] = geographicLocation[geo_index++];
       lonString[lon_index++] = geographicLocation[geo_index++];
       ++lon_index;
 
       // seconds longitude
       lonString[lon_index++] = geographicLocation[geo_index++];
       lonString[lon_index++] = geographicLocation[geo_index++];
       ++lon_index;
 
       // hemisphere
       lonString[lon_index++] = geographicLocation[geo_index++];
 
       // Convert to decimal degrees using the dms class.
       ossimGpt gpt;
       ossimDms dms(0.0);
       dms.setLatFlag(true);
       if ( ! dms.setDegrees(latString.c_str()) )
       {
          gpts.clear();
          return;
       }
       gpt.latd(dms.getDegrees());
 
       dms.setLatFlag(false);
       if ( ! dms.setDegrees(lonString.c_str()) )
       {
          gpts.clear();
          return;
       }
       gpt.lond(dms.getDegrees());
 
       gpts.push_back(gpt);
 
       if (traceDebug())
       {
          ossimNotify(ossimNotifyLevel_DEBUG)
             << "ossimNitfProjectionFactory::parseGeographicString DEBUG:"
             << "\nground point[" << i << "]:  " << gpt
             << std::endl;
       }
    }
 }

◆ parseMgrsString()

bool ossimNitfProjectionFactory::parseMgrsString	(	const ossimString &	mgrsLocationString,
		ossim_uint32 &	zone,
		char &	hemisphere,
		std::vector< ossimDpt > &	utmPoints
	)		const

private

Gets UTM points, zone and hemisphere from IGEOLO field when ICORDS files = U which is UTM expressed in MGRS.

Parameters

mgrsLocationString	String from IGEOLO field.
zone	Initialized in method.
hemisphere	Initialized in method.
utmPoints	Initialized in method.

Returns: True on success, false on error.

Definition at line 518 of file ossimNitfProjectionFactory.cpp.

References Convert_OSSIM_MGRS_To_UTM(), n, ossimString::size(), and ossimString::substr().

Referenced by makeUtm().

 {
    bool result = false; // Start false.
    
    //---
    // From spec:
    // UTM expressed in MGRS use the format zzBJKeeeeennnnn (15 characters).
    // Assumption Zone and hemisphere same for all corners.
    //---
    if ( mgrsLocationString.size() >= 60 )
    {
       // Split the location string into four separate ones.
       std::vector<std::string> mgrsStr(4); // Corner strings.
       mgrsStr[0] = mgrsLocationString.substr(0, 15);
       mgrsStr[1] = mgrsLocationString.substr(15, 15);
       mgrsStr[2] = mgrsLocationString.substr(30, 15);
       mgrsStr[3] = mgrsLocationString.substr(45, 15);
       
       utmPoints.resize(4);
       long z = 0;
       ossim_float64 e=0.0;
       ossim_float64 n=0.0;
       
       result = true; // Set to true.
 
 
       //---
       // Convert each string to Easting Northing.  This also sets zone hemisphere.
       // Method takes long for zone.
       //---
       for (ossim_uint32 i = 0; i < 4; ++i)
       {
          if ( Convert_OSSIM_MGRS_To_UTM(mgrsStr[i].c_str(), &z, &hemisphere, &e, &n) == 0 )
          {
             utmPoints[i].x = e;
             utmPoints[i].y = n;
          }
          else
          {
             result = false; // Geotrans code errored on string.
             break;
          }
       }
       if (result) zone = static_cast<ossim_uint32>(z); // Set the zone.
    }
 
    return result;
 }

◆ parseUtmString()

void ossimNitfProjectionFactory::parseUtmString	(	const ossimString &	utmLocation,
		ossim_uint32 &	zone,
		std::vector< ossimDpt > &	utmPoints
	)		const

private

Definition at line 909 of file ossimNitfProjectionFactory.cpp.

References ossimString::begin(), ossimString::push_back(), ossimString::toDouble(), and ossimString::toUInt32().

Referenced by makeUtm().

 {
    ossim_uint32 idx = 0;
    ossimString z;
    ossimString east;
    ossimString north;
 
    
    z    = ossimString(utmLocation.begin() + idx,
                    utmLocation.begin() + idx + 2);
    idx += 2;
    east = ossimString(utmLocation.begin() + idx,
                       utmLocation.begin() + idx + 6);
    idx += 6;
    north = ossimString(utmLocation.begin() + idx,
                        utmLocation.begin() + idx + 7);
    idx += 7;
 
    utmPoints.push_back(ossimDpt(east.toDouble(),
                                 north.toDouble()));
    
    z    = ossimString(utmLocation.begin() + idx,
                    utmLocation.begin() + idx + 2);
    idx += 2;
    east = ossimString(utmLocation.begin() + idx,
                       utmLocation.begin() + idx + 6);
    idx += 6;
    north = ossimString(utmLocation.begin() + idx,
                        utmLocation.begin() + idx + 7);
    idx += 7;
    utmPoints.push_back(ossimDpt(east.toDouble(),
                                 north.toDouble()));
 
    z    = ossimString(utmLocation.begin() + idx,
                    utmLocation.begin() + idx + 2);
    idx += 2;
    east = ossimString(utmLocation.begin() + idx,
                       utmLocation.begin() + idx + 6);
    idx += 6;
    north = ossimString(utmLocation.begin() + idx,
                        utmLocation.begin() + idx + 7);
    idx += 7;
    utmPoints.push_back(ossimDpt(east.toDouble(),
                                 north.toDouble()));
 
    z    = ossimString(utmLocation.begin() + idx,
                    utmLocation.begin() + idx + 2);
    idx += 2;
    east = ossimString(utmLocation.begin() + idx,
                       utmLocation.begin() + idx + 6);
    idx += 6;
    north = ossimString(utmLocation.begin() + idx,
                        utmLocation.begin() + idx + 7);
    idx += 7;
    utmPoints.push_back(ossimDpt(east.toDouble(),
                                 north.toDouble()));
 
    zone = z.toUInt32();
 }

Member Data Documentation

◆ theInstance

ossimNitfProjectionFactory * ossimNitfProjectionFactory::theInstance = 0

staticprivate

Definition at line 266 of file ossimNitfProjectionFactory.h.

Referenced by instance().

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ~ossimNitfProjectionFactory()

◆ ossimNitfProjectionFactory() [1/2]

◆ ossimNitfProjectionFactory() [2/2]

Member Function Documentation

◆ computeScaleInDecimalDegrees()

◆ computeScaleInMeters()

◆ createObject() [1/2]

◆ createObject() [2/2]

◆ createProjection() [1/4]

◆ createProjection() [2/4]

◆ createProjection() [3/4]

◆ createProjection() [4/4]

◆ createProjectionFromHeaders()

◆ getBlockaPoints()

◆ getGeolobPoints()

◆ getTypeNameList()

◆ instance()

◆ isNitf()

◆ isSkewed() [1/2]

◆ isSkewed() [2/2]

◆ makeBilinear()

◆ makeEuiDistant()

◆ makeGeographic()

◆ makeUtm()

◆ operator=()

◆ parseDecimalDegreesString()

◆ parseGeographicString()

◆ parseMgrsString()

◆ parseUtmString()

Member Data Documentation

◆ theInstance