ossimGmlSupportData Class Reference

GML support data class. More...

#include <ossimGmlSupportData.h>

Public Member Functions
	ossimGmlSupportData ()
	default constructor More...
	~ossimGmlSupportData ()
	destructor More...
bool	initialize (const ossimImageGeometry *geom, const ossimIrect &rect)
	Initializes gml block from geometry file. More...
bool	initialize (std::istream &in)
	Initializes from stream assumed to be a gml block. More...
bool	write (std::ostream &os)
	Write xml doc to stream. More...
ossimRefPtr< ossimXmlDocument >	getXmlDoc () const
bool	getImageGeometry (ossimKeywordlist &geomKwl) const
	Extracts geometry info to keyword list. More...

Private Member Functions
ossimRefPtr< ossimXmlNode >	getGmljp2V1RootNode () const
ossimRefPtr< ossimXmlNode >	getGmljp2V2RootNode () const
void	getGeoOrigin (const ossimImageGeometry *geom, const ossimIpt &ul, ossimString &originString, ossimString &offsetVector1String, ossimString &offsetVector2String) const
void	getMapOrigin (const ossimImageGeometry *geom, const ossimIpt &ul, ossimString &originString, ossimString &offsetVector1String, ossimString &offsetVector2String) const
void	getLimits (const ossimIrect &rect, ossimString &gridHighString, ossimString &gridLowString) const
void	getGeoBounds (const ossimImageGeometry *geom, const ossimIrect &rect, ossimString &upperCornerString, ossimString &lowerCornerString) const
bool	configureGmljp2V2 (ossimRefPtr< ossimXmlNode > node, const ossimImageGeometry *geom, const ossimIrect &rect)
bool	configureBounds (ossimRefPtr< ossimXmlNode > node, const ossimImageGeometry *geom, const ossimIrect &rect)
bool	getImageGeometryFromSeonsorModel (ossimKeywordlist &geomKwl) const
	For sensor model data: More...
bool	getImageGeometryFromRectifiedGrid (ossimKeywordlist &geomKwl) const
	For map projected data: More...
bool	addLineSamps (const ossimString &xpath0, bool useGmlPrefix, ossimKeywordlist &geomKwl) const
	Adds line and sample to keyword list. More...
bool	addTie (const ossimString &xpath0, bool useGmlPrefix, ossim_uint32 pcsCode, ossimKeywordlist &geomKwl) const
	Adds tie point to keyword list. More...
bool	addScale (const ossimString &xpath0, bool useGmlPrefix, ossim_uint32 pcsCode, ossimKeywordlist &geomKwl) const
	Adds tie and scale to keyword list. More...

Private Attributes
ossimRefPtr< ossimXmlDocument >	m_xmlDocument
ossim_uint32	m_pcsCodeMap
ossim_uint32	m_pcsCodeGeo
ossimString	m_srsNameStringMap
ossimString	m_srsNameStringGeo
ossimString	m_srsDimensionString
ossimString	m_axisLabelsStringMap
ossimString	m_axisLabelsStringGeo
ossimString	m_uomLabelsStringMap
ossimString	m_uomLabelsStringGeo

Detailed Description

GML support data class.

Definition at line 29 of file ossimGmlSupportData.h.

Constructor & Destructor Documentation

ossimGmlSupportData()

ossimGmlSupportData::ossimGmlSupportData

(

)

default constructor

Definition at line 40 of file ossimGmlSupportData.cpp.

   :
   m_xmlDocument(0),
   // m_use_gmljp2_version2(true), 
   m_pcsCodeMap(UNKNOWN_PCSCODE),
   m_pcsCodeGeo(UNKNOWN_PCSCODE),
   m_srsNameStringMap("http://www.opengis.net/def/crs/EPSG/0/32767"),
   m_srsNameStringGeo("http://www.opengis.net/def/crs/EPSG/0/32767"),
   m_srsDimensionString("2"),
   m_axisLabelsStringMap("X Y"),
   m_axisLabelsStringGeo("Lat Long"),
   m_uomLabelsStringMap("m m"),
   m_uomLabelsStringGeo("deg deg")
{
}

~ossimGmlSupportData()

ossimGmlSupportData::~ossimGmlSupportData

(

)

destructor

Definition at line 56 of file ossimGmlSupportData.cpp.

{
}

Member Function Documentation

addLineSamps()

bool ossimGmlSupportData::addLineSamps ( const ossimString &  xpath0, bool  useGmlPrefix, ossimKeywordlist &  geomKwl  ) const

private

Adds line and sample to keyword list.

Parameters

xpath0
useGmlPrefix	If true tack on "gml:" to paths.
geomKwl	Initialized by this.

Definition at line 1177 of file ossimGmlSupportData.cpp.

References ossimKeywordlist::add(), ossimString::afterPos(), ossimString::beforePos(), ossimString::clear(), ossimString::find(), ossimXmlDocument::findNodes(), m_xmlDocument, ossimKeywordNames::NUMBER_LINES_KW, ossimKeywordNames::NUMBER_SAMPLES_KW, status, ossimString::toInt32(), and ossimRefPtr< T >::valid().

Referenced by getImageGeometryFromRectifiedGrid().

{
   bool status = false;
   
   if ( m_xmlDocument.valid() )
   {
      vector< ossimRefPtr<ossimXmlNode> > xml_nodes;

      // Number of lines & samples, for either sensor or rectified imagery:
      ossimString xpath_limits_low;
      if (useGmlPrefix)
      {
         xpath_limits_low = "/gml:limits/gml:GridEnvelope/gml:low";
      }
      else
      {
         xpath_limits_low = "/limits/GridEnvelope/low";
      }
      ossimString xpath = xpath0 + xpath_limits_low;
       m_xmlDocument->findNodes( xpath, xml_nodes );
      if ( xml_nodes.size() == 1 )
      {
         const ossimString& lowerCorner = xml_nodes[0]->getText();
         size_t spacePos = lowerCorner.find( ' ' );
         ossimString lowerXString = lowerCorner.beforePos( spacePos );
         ossimString lowerYString = lowerCorner.afterPos ( spacePos );
         ossim_uint32 lowX = lowerXString.toInt32();
         ossim_uint32 lowY = lowerYString.toInt32();

         ossimString xpath_limits_high;
         if ( useGmlPrefix )
         {
            xpath_limits_high = "/gml:limits/gml:GridEnvelope/gml:high";
         }
         else
         {
            xpath_limits_high = "/limits/GridEnvelope/high";
         }

         xpath = xpath0 + xpath_limits_high;
         xml_nodes.clear();
         m_xmlDocument->findNodes( xpath, xml_nodes );
         if ( xml_nodes.size() == 1 )
         {
            const ossimString& higherCorner = xml_nodes[0]->getText();
            size_t spacePos = higherCorner.find( ' ' );
            ossimString higherXString = higherCorner.beforePos( spacePos );
            ossimString higherYString = higherCorner.afterPos ( spacePos );
            ossim_uint32 highX = higherXString.toInt32();
            ossim_uint32 highY = higherYString.toInt32();
            
            geomKwl.add( ossimKeywordNames::NUMBER_LINES_KW,   highY - lowY + 1 );
            geomKwl.add( ossimKeywordNames::NUMBER_SAMPLES_KW, highX - lowX + 1 );
            
            status = true;
         }
      }
   }

   return status;
   
} // ossimGmlSupportData::addLineSamps

addScale()

bool ossimGmlSupportData::addScale ( const ossimString &  xpath0, bool  useGmlPrefix, ossim_uint32  pcsCode, ossimKeywordlist &  geomKwl  ) const

private

Adds tie and scale to keyword list.

Parameters

xpath0
useGmlPrefix	If true tack on "gml:" to paths.
geomKwl	Initialized by this.

Definition at line 1338 of file ossimGmlSupportData.cpp.

References ossim::acosd(), ossimKeywordlist::add(), ossimString::afterPos(), ossimString::beforePos(), ossimString::c_str(), ossimString::clear(), ossimString::find(), ossimXmlDocument::findNodes(), m_xmlDocument, ossimKeywordNames::ORIGIN_LATITUDE_KW, ossimKeywordNames::PIXEL_SCALE_UNITS_KW, ossimKeywordNames::PIXEL_SCALE_XY_KW, status, ossimString::toFloat64(), ossimDpt::toString(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

Referenced by getImageGeometryFromRectifiedGrid().

{
   bool status = false;

   if ( m_xmlDocument.valid() && (pcsCode != 32767 ) )
   {
      vector< ossimRefPtr<ossimXmlNode> > xml_nodes;

      //---
      // origin:
      // Note: In GML the origin is the tie point, NOT the projection origin.
      //---

      // axis labels for rectified imagery:
      ossimString xpath_axisLabels;
      if ( useGmlPrefix )
      {
          xpath_axisLabels = "/gml:axisLabels";
      }
      else
      {
         xpath_axisLabels = "/axisLabels";
      }

      ossimString xpath = xpath0 + xpath_axisLabels;
      m_xmlDocument->findNodes( xpath, xml_nodes );
      if ( xml_nodes.size() == 1 )
      {
         ossimString axisLabelsString = xml_nodes[0]->getText();

         ossimString xpath_offsetVector;
         if ( useGmlPrefix )
         {
            xpath_offsetVector = "/gml:offsetVector";
         }
         else
         {
            xpath_offsetVector =  "/offsetVector";
         }
         
         xpath = xpath0 + xpath_offsetVector;
         xml_nodes.clear();
         m_xmlDocument->findNodes( xpath, xml_nodes );
         if ( xml_nodes.size() == 2 )
         {
            const ossimString& offsetVectorString0 = xml_nodes[0]->getText();
            const ossimString& offsetVectorString1 = xml_nodes[1]->getText();
            
            size_t spacePos0 = offsetVectorString0.find( ' ' );
            size_t spacePos1 = offsetVectorString1.find( ' ' );            
            
            ossimString firstOffsetVectorString0  = offsetVectorString0.beforePos( spacePos0 );
            ossimString secondOffsetVectorString0 = offsetVectorString0.afterPos ( spacePos0 );
            ossimString firstOffsetVectorString1  = offsetVectorString1.beforePos( spacePos1 );
            ossimString secondOffsetVectorString1 = offsetVectorString1.afterPos ( spacePos1 );
            
            // TODO: Add for rotational matrix:
            
            std::string scale_xy;
            std::string scale_units;
            if ( pcsCode == 4326 )
            {
               scale_units = "degrees";
            }
            else
            {
               scale_units = "meters";
            }

            ossimDpt scale;
            if ( ( scale_units == "degrees" ) &&
                 ( axisLabelsString == "Lat Long" ) )
            {
               scale.x = secondOffsetVectorString1.toFloat64();
               scale.y = firstOffsetVectorString0.toFloat64();
            }
            else
            {
               scale.x = firstOffsetVectorString0.toFloat64();
               scale.y = secondOffsetVectorString1.toFloat64();
            }
            if ( scale.y < 0.0 )
            {
               scale.y *= -1.0; // make positive
            }

            if ( ( scale_units == "degrees" ) && ( scale.x != scale.y ) )
            {
               // Compute/add the origin latitude of true scale.
               ossim_float64 origin_lat = ossim::acosd(scale.y/scale.x);

               geomKwl.add(  ossimKeywordNames::ORIGIN_LATITUDE_KW, origin_lat );
            }
            
            geomKwl.add( ossimKeywordNames::PIXEL_SCALE_XY_KW, scale.toString().c_str() );
            geomKwl.add( ossimKeywordNames::PIXEL_SCALE_UNITS_KW, scale_units.c_str() );                
            status = true;
         }
      }

   } // Matches: if ( m_xmlDocument.valid() && (pscCode != 32767 ) )

   return status;
   
} // End: ossimGmlSupportData::addScale( ... )

addTie()

bool ossimGmlSupportData::addTie ( const ossimString &  xpath0, bool  useGmlPrefix, ossim_uint32  pcsCode, ossimKeywordlist &  geomKwl  ) const

private

Adds tie point to keyword list.

Parameters

xpath0
useGmlPrefix	If true tack on "gml:" to paths.
geomKwl	Initialized by this.

Definition at line 1242 of file ossimGmlSupportData.cpp.

References ossimKeywordlist::add(), ossimString::afterPos(), ossimString::beforePos(), ossimString::clear(), ossimString::find(), ossimXmlDocument::findNodes(), m_xmlDocument, status, ossimString::string(), ossimKeywordNames::TIE_POINT_UNITS_KW, ossimKeywordNames::TIE_POINT_XY_KW, and ossimRefPtr< T >::valid().

Referenced by getImageGeometryFromRectifiedGrid().

{
   bool status = false;

   if ( m_xmlDocument.valid() && (pcsCode != 32767 ) )
   {
      vector< ossimRefPtr<ossimXmlNode> > xml_nodes;

      //---
      // origin:
      // Note: In GML the origin is the tie point, NOT the projection origin.
      //---

      // axis labels for rectified imagery:
      ossimString xpath_axisLabels;
      if ( useGmlPrefix )
      {
          xpath_axisLabels = "/gml:axisLabels";
      }
      else
      {
         xpath_axisLabels = "/axisLabels";
      }

      ossimString xpath = xpath0 + xpath_axisLabels;
      m_xmlDocument->findNodes( xpath, xml_nodes );
      if ( xml_nodes.size() == 1 )
      {
         ossimString axisLabelsString = xml_nodes[0]->getText();

         ossimString xpath_originPoint;

         if ( useGmlPrefix )
         {
            xpath_originPoint = "/gml:origin/gml:Point";
         }
         else
         {
            xpath_originPoint = "/origin/Point";
         }

         xpath = xpath0 + xpath_originPoint;
         xml_nodes.clear();
         m_xmlDocument->findNodes( xpath, xml_nodes );
         if ( xml_nodes.size() == 1 )
         {
            const ossimString& originString = xml_nodes[0]->getChildTextValue( ossimString( "pos" ) );
            size_t spacePos = originString.find( ' ' );
            ossimString firstOriginString  = originString.beforePos( spacePos );
            ossimString secondOriginString = originString.afterPos ( spacePos );
            
            std::string tie_point_xy;
            std::string tie_point_units;
            
            if ( pcsCode == 4326 ) // map projection
            {
               // Longitude first, e.g. (lon,lat)
               tie_point_units = "degrees";
            }
            else
            {
               tie_point_units = "meters";
            }
            
            if ( ( tie_point_units == "degrees" ) &&
                 ( axisLabelsString == "Lat Long" ) )
            {
               tie_point_xy = "(";
               tie_point_xy += secondOriginString.string();
               tie_point_xy += ",";
               tie_point_xy += firstOriginString.string();
               tie_point_xy += ")";
            }
            else
            {
               tie_point_xy = "(";
               tie_point_xy += firstOriginString.string();
               tie_point_xy += ",";
               tie_point_xy += secondOriginString.string();
               tie_point_xy += ")";
            }
            geomKwl.add( ossimKeywordNames::TIE_POINT_XY_KW, tie_point_xy.c_str() );
            geomKwl.add( ossimKeywordNames::TIE_POINT_UNITS_KW, tie_point_units.c_str() );
            status = true;
         }
      }
      
   } // Matches: if ( m_xmlDocument.valid() && (pscCode != 32767 ) )

   return status;
   
} // End: ossimGmlSupportData::addTie( ... )

configureBounds()

bool ossimGmlSupportData::configureBounds ( ossimRefPtr< ossimXmlNode >  node, const ossimImageGeometry *  geom, const ossimIrect &  rect  )

private

Definition at line 377 of file ossimGmlSupportData.cpp.

References ossimXmlNode::addAttribute(), ossimXmlNode::addChildNode(), getGeoBounds(), m_axisLabelsStringGeo, m_srsDimensionString, m_srsNameStringGeo, m_uomLabelsStringGeo, and ossimXmlAttribute::setNameValue().

Referenced by configureGmljp2V2().

{
   bool success = true;
   const ossimString BLANK = "";

   ossimString upperCornerString;         
   ossimString lowerCornerString;
   getGeoBounds( geom, rect, upperCornerString, lowerCornerString );

   ossimString path = "boundedBy";
   ossimRefPtr<ossimXmlNode> node1 = node0->addChildNode( path, BLANK );

   path = "Envelope";
   ossimRefPtr<ossimXmlNode> node1a = node1->addChildNode( path, BLANK );

   ossimRefPtr<ossimXmlAttribute> attr( 0 );
   ossimString name;

   attr = new ossimXmlAttribute();
   name = "srsName";
   attr->setNameValue( name, m_srsNameStringGeo );
   node1a->addAttribute( attr );

   attr = new ossimXmlAttribute();
   name = "axisLabels";
   attr->setNameValue( name, m_axisLabelsStringGeo );
   node1a->addAttribute( attr );

   attr = new ossimXmlAttribute();
   name = "uomLabels";
   attr->setNameValue( name, m_uomLabelsStringGeo );
   node1a->addAttribute( attr );

   attr = new ossimXmlAttribute();
   name = "srsDimension";
   attr->setNameValue( name, m_srsDimensionString );
   node1a->addAttribute( attr );

   path = "lowerCorner";
   ossimRefPtr<ossimXmlNode> node1a1 =
      node1a->addChildNode( path, lowerCornerString );

   path = "upperCorner";
   ossimRefPtr<ossimXmlNode> node1a2 =
      node1a->addChildNode( path, upperCornerString );

   return success;
}

configureGmljp2V2()

bool ossimGmlSupportData::configureGmljp2V2 ( ossimRefPtr< ossimXmlNode >  node, const ossimImageGeometry *  geom, const ossimIrect &  rect  )

private

Definition at line 183 of file ossimGmlSupportData.cpp.

References ossimXmlNode::addAttribute(), ossimXmlNode::addChildNode(), configureBounds(), ossimImageGeometry::getAsMapProjection(), getGeoOrigin(), getLimits(), getMapOrigin(), ossimMapProjection::isGeographic(), m_axisLabelsStringGeo, m_axisLabelsStringMap, m_srsNameStringGeo, m_srsNameStringMap, ossimXmlAttribute::setNameValue(), ossimIrect::ul(), and ossimRefPtr< T >::valid().

Referenced by initialize().

{
   bool success = false;

   ossimRefPtr<const ossimMapProjection> mapProj = geom->getAsMapProjection();
   if ( mapProj.valid() )
   {
      bool isGeographic = mapProj->isGeographic();
      
      const ossimString BLANK = "";
      
      ossimRefPtr<ossimXmlAttribute> attr(0);
      ossimString name;
      ossimString value;
      
      ossimString gridHighString;
      ossimString gridLowString;
      getLimits( rect, gridHighString, gridLowString );
      
      configureBounds( node0, geom, rect );
      
      ossimString path = "domainSet";
      ossimRefPtr<ossimXmlNode> node2 =
         node0->addChildNode( path, BLANK );
      
      path = "rangeSet";
      ossimRefPtr<ossimXmlNode> node3 =
         node0->addChildNode( path, BLANK );
      
      path = "File";
      ossimRefPtr<ossimXmlNode> node3a =
         node3->addChildNode( path, BLANK );
      
      path = "rangeParameters";
      ossimRefPtr<ossimXmlNode> node3a1 =
         node3a->addChildNode( path, BLANK );
      
      path = "fileName";
      ossimRefPtr<ossimXmlNode> node3a2 =
         node3a->addChildNode( path, "gmljp2://codestream" );
      
      path = "fileStructure";
      ossimRefPtr<ossimXmlNode> node3a3 =
         node3a->addChildNode( path, "inapplicable" );
      
      path = "gmlcov:rangeType";
      ossimRefPtr<ossimXmlNode> node4 =
         node0->addChildNode( path, BLANK );
      
      path = "gmljp2:featureMember";
      ossimRefPtr<ossimXmlNode> node5 =
         node0->addChildNode( path, BLANK );
      
      path = "gmljp2:GMLJP2RectifiedGridCoverage";
      ossimRefPtr<ossimXmlNode> node5a =
         node5->addChildNode( path, BLANK );

      attr  = new ossimXmlAttribute();
      name  = "gml:id";
      value = "CodeStream_0";
      attr->setNameValue( name, value );
      node5a->addAttribute( attr );

      path = "domainSet";
      ossimRefPtr<ossimXmlNode> node5a1 =
         node5a->addChildNode( path, BLANK );

      path = "RectifiedGrid";
      ossimRefPtr<ossimXmlNode> node5a1a =
         node5a1->addChildNode( path, BLANK );

      attr  = new ossimXmlAttribute();
      name  = "gml:id";
      value = "RG0001";
      attr->setNameValue( name, value );
      node5a1a->addAttribute( attr );

      attr  = new ossimXmlAttribute();
      name  = "dimension";
      value = "2";
      attr->setNameValue( name, value );
      node5a1a->addAttribute( attr );

      attr  = new ossimXmlAttribute();
      name  = "srsName";
      attr->setNameValue( name, m_srsNameStringMap );
      node5a1a->addAttribute( attr );

      path = "limits";
      ossimRefPtr<ossimXmlNode> node5a1a1 =
         node5a1a->addChildNode( path, BLANK );

      path = "GridEnvelope";
      ossimRefPtr<ossimXmlNode> node5a1a1a =
         node5a1a1->addChildNode( path, BLANK );

      path = "low";
      ossimRefPtr<ossimXmlNode> node5a1a1a1 =
         node5a1a1a->addChildNode( path, gridLowString );

      path = "high";
      ossimRefPtr<ossimXmlNode> node5a1a1a2 =
         node5a1a1a->addChildNode( path, gridHighString );

      path = "axisLabels";
      ossimRefPtr<ossimXmlNode> node5a1a2 =
         node5a1a->addChildNode( path, 
                                 mapProj->isGeographic() ? 
                                 m_axisLabelsStringGeo : m_axisLabelsStringMap );

      path = "origin";
      ossimRefPtr<ossimXmlNode> node5a1a3 =
         node5a1a->addChildNode( path, BLANK );

      path = "Point";
      ossimRefPtr<ossimXmlNode> node5a1a3a =
         node5a1a3->addChildNode( path, BLANK );

      attr  = new ossimXmlAttribute();
      name  = "gml:id";
      value = "P0001";
      attr->setNameValue( name, value );
      node5a1a3a->addAttribute( attr );

      attr  = new ossimXmlAttribute();
      name  = "srsName";
      attr->setNameValue( name, m_srsNameStringMap );
      node5a1a3a->addAttribute( attr );

      ossimString originString;
      ossimString offsetVector1String;
      ossimString offsetVector2String;
      if ( isGeographic )
      {
         getGeoOrigin( geom, rect.ul(), originString,
                       offsetVector1String, offsetVector2String );
      }
      else
      {
         getMapOrigin( geom, rect.ul(), originString,
                       offsetVector1String, offsetVector2String );
      }

      path = "pos";
      ossimRefPtr<ossimXmlNode> node5a1a3a1 =
         node5a1a3a->addChildNode( path, originString );

      path = "offsetVector";
      ossimRefPtr<ossimXmlNode> node5a1a4 =
         node5a1a->addChildNode( path, offsetVector1String );

      attr  = new ossimXmlAttribute();
      name  = "srsName";
      attr->setNameValue( name, 
                          (isGeographic ? m_srsNameStringGeo : m_srsNameStringMap ) );
      node5a1a4->addAttribute( attr );

      path = "offsetVector";
      ossimRefPtr<ossimXmlNode> node5a1a5 =
         node5a1a->addChildNode( path, offsetVector2String );

      attr  = new ossimXmlAttribute();
      name  = "srsName";
      attr->setNameValue( name, (isGeographic ? m_srsNameStringGeo : m_srsNameStringMap) );
      node5a1a5->addAttribute( attr );

      path = "rangeSet";
      ossimRefPtr<ossimXmlNode> node5a2 =
         node5a->addChildNode( path, BLANK );

      path = "File";
      ossimRefPtr<ossimXmlNode> node5a2a =
         node5a2->addChildNode( path, BLANK );

      path = "rangeParameters";
      ossimRefPtr<ossimXmlNode> node5a2a1 =
         node5a2a->addChildNode( path, BLANK );

      path = "fileName";
      ossimRefPtr<ossimXmlNode> node5a2a2 =
         node5a2a->addChildNode( path, "gmljp2://codestream" );

      path = "fileStructure";
      ossimRefPtr<ossimXmlNode> node5a2a3 =
         node5a2a->addChildNode( path, "inapplicable" );

      success = true;
   }
   
   return success;
}

getGeoBounds()

void ossimGmlSupportData::getGeoBounds ( const ossimImageGeometry *  geom, const ossimIrect &  rect, ossimString &  upperCornerString, ossimString &  lowerCornerString  ) const

private

Definition at line 1048 of file ossimGmlSupportData.cpp.

References ossimImageGeometry::getAsMapProjection(), ossimGpt::lat, ossimImageGeometry::localToWorld(), ossimGpt::lon, ossimIrect::lr(), os2, ossimIrect::ul(), and ossimRefPtr< T >::valid().

Referenced by configureBounds().

{
   if ( geom )
   {
      ossimRefPtr<const ossimMapProjection> mapProj = geom->getAsMapProjection();
      if ( mapProj.valid() )
      {
         // Get the bounding rect.  This assumes North up.
         ossimGpt ulGpt;
         ossimGpt lrGpt;
         geom->localToWorld( rect.ul(), ulGpt );
         geom->localToWorld( rect.lr(), lrGpt );
         
         std::ostringstream os;
         os.precision(15);
         os << ulGpt.lat << " " << ulGpt.lon;
         upperCornerString = os.str();
         
         std::ostringstream os2;
         os2.precision(15);
         os2 << lrGpt.lat << " " << lrGpt.lon;
         lowerCornerString = os2.str();
      }
   }
   
} // End: ossimGmlSupportData::getGeoBounds

getGeoOrigin()

void ossimGmlSupportData::getGeoOrigin ( const ossimImageGeometry *  geom, const ossimIpt &  ul, ossimString &  originString, ossimString &  offsetVector1String, ossimString &  offsetVector2String  ) const

private

Definition at line 971 of file ossimGmlSupportData.cpp.

References ossimImageGeometry::getDegreesPerPixel(), ossimGpt::lat, ossimImageGeometry::localToWorld(), ossimGpt::lon, os2, ossimDpt::x, and ossimDpt::y.

Referenced by configureGmljp2V2().

{
   if ( geom )
   {
      // Get the gsd:
      ossimDpt gsd;
      geom->getDegreesPerPixel( gsd );

      // Get tie point:
      ossimGpt tie;
      geom->localToWorld( ul, tie );

      std::ostringstream os;
      os.precision(15);
      os << tie.lat << " " << tie.lon;
      originString = os.str();

      std::ostringstream os2;
      os2.precision(15);
      os2 << -gsd.y << " " << "0.0";
      offsetVector1String = os2.str();

      std::ostringstream os3;
      os3.precision(15);
      os3 << "0.0" << " " << gsd.x;
      offsetVector2String = os3.str();
   }

} // End: ossimGmlSupportData::getGeoOrigin

getGmljp2V1RootNode()

ossimRefPtr< ossimXmlNode > ossimGmlSupportData::getGmljp2V1RootNode ( ) const

private

Definition at line 486 of file ossimGmlSupportData.cpp.

References ossimXmlNode::addAttribute(), ossimXmlAttribute::setNameValue(), and ossimXmlNode::setTag().

Referenced by initialize().

{
   ossimRefPtr<ossimXmlNode> node = new ossimXmlNode();

   ossimString os = "gml:FeatureCollection";
   node->setTag( os );

   ossimRefPtr<ossimXmlAttribute> attr(0);
   ossimString name;
   ossimString value;
   
   attr = new ossimXmlAttribute();
   name = "xmlns";
   value = "http://www.opengis.net/gml";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();   
   name = "xmlns:gml";
   value = "http://www.opengis.net/gml";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();     
   name = "xmlns:xsi";
   value = "http://www.w3.org/2001/XMLSchema-instance";
   attr->setNameValue( name, value );
   node->addAttribute( attr );
      
   attr = new ossimXmlAttribute();   
   name = "xsi:schemaLocation";
   value = "http://www.opengis.net/gml gmlJP2Profile.xsd";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   return node;
   
} // ossimGmlSupportData::getGmljp2V1Node()

getGmljp2V2RootNode()

ossimRefPtr< ossimXmlNode > ossimGmlSupportData::getGmljp2V2RootNode ( ) const

private

Definition at line 429 of file ossimGmlSupportData.cpp.

References ossimXmlNode::addAttribute(), ossimXmlAttribute::setNameValue(), and ossimXmlNode::setTag().

Referenced by initialize().

{
   ossimRefPtr<ossimXmlNode> node = new ossimXmlNode();

   ossimString os = "gmljp2:GMLJP2CoverageCollection";
   node->setTag( os );

   ossimRefPtr<ossimXmlAttribute> attr( 0 );
   ossimString name;
   ossimString value;

   attr = new ossimXmlAttribute();
   name = "gml:id";
   value = "JPEG2000_0";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();
   name = "xmlns";
   value = "http://www.opengis.net/gml/3.2";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();   
   name = "xmlns:gml";
   value = "http://www.opengis.net/gml/3.2";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();   
   name = "xmlns:gmlcov";
   value = "http://www.opengis.net/gmlcov/1.0";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();   
   name = "xmlns:gmljp2";
   value = "http://www.opengis.net/gmljp2/2.0";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   attr = new ossimXmlAttribute();     
   name = "xmlns:xsi";
   value = "http://www.w3.org/2001/XMLSchema-instance";
   attr->setNameValue( name, value );
   node->addAttribute( attr );
      
   attr = new ossimXmlAttribute();   
   name = "xsi:schemaLocation";
   value = "http://www.opengis.net/gmljp2/2.0 http://schemas.opengis.net/gmljp2/2.0/gmljp2.xsd";
   attr->setNameValue( name, value );
   node->addAttribute( attr );

   return node;
   
} // ossimGmlSupportData::getGmljp2V2Node()

getImageGeometry()

bool ossimGmlSupportData::getImageGeometry

(

ossimKeywordlist &

geomKwl

)

const

Extracts geometry info to keyword list.

Populates the keyword list with image geometry information. This method is used to relay projection/model information to users.

Parameters

geomKwl

Keyword list that will be initialized with geometry info.

Returns

true if geometry info is present, false if not.

Definition at line 544 of file ossimGmlSupportData.cpp.

References getImageGeometryFromRectifiedGrid(), and getImageGeometryFromSeonsorModel().

{
   // Try for map projected geometry first...
   bool success = getImageGeometryFromRectifiedGrid( geomKwl );
   if ( !success )
   {
      // Look for sensor model block.
      getImageGeometryFromSeonsorModel( geomKwl );
   }

   return success;
}

getImageGeometryFromRectifiedGrid()

bool ossimGmlSupportData::getImageGeometryFromRectifiedGrid ( ossimKeywordlist &  geomKwl ) const

private

For map projected data:

Parameters

geomKwl

Initialized by this.

Definition at line 1117 of file ossimGmlSupportData.cpp.

References ossimKeywordlist::add(), addLineSamps(), addScale(), addTie(), ossimString::after(), ossimString::c_str(), ossimXmlDocument::findNodes(), ossimXmlAttribute::getValue(), m_xmlDocument, ossimKeywordNames::PCS_CODE_KW, ossimKeywordNames::SRS_NAME_KW, status, ossimString::string(), ossimString::toUInt32(), ossimKeywordNames::TYPE_KW, and ossimRefPtr< T >::valid().

Referenced by getImageGeometry().

{
   bool status = false;

   if ( m_xmlDocument.valid() )
   {
      bool useGmlPrefix = false;
      
      vector< ossimRefPtr<ossimXmlNode> > xml_nodes;
      
      ossimString xpath0 = "/gmljp2:GMLJP2CoverageCollection/gmljp2:featureMember/gmljp2:GMLJP2RectifiedGridCoverage/domainSet/RectifiedGrid";
      m_xmlDocument->findNodes( xpath0, xml_nodes );
      if ( !xml_nodes.size() )
      {
         xpath0 = "/gmljp2:GMLJP2CoverageCollection/gmljp2:featureMember/gmljp2:GMLJP2RectifiedGridCoverage/gml:domainSet/gml:RectifiedGrid";
        m_xmlDocument->findNodes( xpath0, xml_nodes );
        if (  xml_nodes.size() )
        {
           useGmlPrefix = true;
        }
      }
         
      if ( xml_nodes.size() )
      {
         const ossimString SRS_NAME( "srsName" );
         const ossimString PCS_CODE_DEFINITION_STR( "http://www.opengis.net/def/crs/EPSG/0/" );
         
         const ossimRefPtr<ossimXmlAttribute> hrefAttribute =
            xml_nodes[0]->findAttribute( SRS_NAME );
         const ossimString& originSrsName = hrefAttribute->getValue();

         ossimString pcsCodeStr = originSrsName.after( PCS_CODE_DEFINITION_STR.string() );
         ossim_uint32 pcsCode = pcsCodeStr.toUInt32();

         if ( pcsCode != 32767 )
         {
            // Add the pcs code.
            geomKwl.add( ossimKeywordNames::PCS_CODE_KW, pcsCodeStr.c_str() );
            if ( pcsCode == 4326 )
            {
               ossimString srsName = "EPSG:";
               srsName += pcsCodeStr;
               geomKwl.add( ossimKeywordNames::SRS_NAME_KW, srsName.c_str() );
               geomKwl.add( ossimKeywordNames::TYPE_KW, "ossimEquDistCylProjection");
            }

            if ( addLineSamps( xpath0, useGmlPrefix, geomKwl ) )
            {
               if ( addTie( xpath0, useGmlPrefix, pcsCode, geomKwl ) )
               {
                  status = addScale( xpath0, useGmlPrefix, pcsCode, geomKwl );
               }
            }
         }
      }
   }

   return status;  
}

getImageGeometryFromSeonsorModel()

bool ossimGmlSupportData::getImageGeometryFromSeonsorModel ( ossimKeywordlist &  geomKwl ) const

private

For sensor model data:

Parameters

geomKwl

Initialized by this.

Definition at line 1096 of file ossimGmlSupportData.cpp.

References ossimXmlDocument::findNodes(), m_xmlDocument, status, and ossimRefPtr< T >::valid().

Referenced by getImageGeometry().

{
   bool status = false;

   if ( m_xmlDocument.valid() )
   {
      vector< ossimRefPtr<ossimXmlNode> > xml_nodes;
      ossimString xpath0 =  "/gmljp2:GMLJP2CoverageCollection/gmljp2:featureMember/gmljp2:GMLJP2ReferenceableGridCoverage/gml:domainSet/gmlcov:ReferenceableGridBySensorModel";

      m_xmlDocument->findNodes( xpath0, xml_nodes );
      if ( xml_nodes.size() >= 1 )
      {
         // Put sensor model code here...

         status = true;
      }
   }

   return status;
}

getLimits()

void ossimGmlSupportData::getLimits ( const ossimIrect &  rect, ossimString &  gridHighString, ossimString &  gridLowString  ) const

private

Definition at line 1078 of file ossimGmlSupportData.cpp.

References ossimIrect::hasNans(), ossimIrect::height(), and ossimIrect::width().

Referenced by configureGmljp2V2().

{
   if ( rect.hasNans() == false )
   {
      // Zero based image rect.
      gridLowString = "0 0";
         
      ossim_uint32 w = rect.width();
      ossim_uint32 h = rect.height();
      
      std::ostringstream os;
      os << (w-1) << " " << (h-1);
      gridHighString = os.str();
   }
}

getMapOrigin()

void ossimGmlSupportData::getMapOrigin ( const ossimImageGeometry *  geom, const ossimIpt &  ul, ossimString &  originString, ossimString &  offsetVector1String, ossimString &  offsetVector2String  ) const

private

Definition at line 1007 of file ossimGmlSupportData.cpp.

References ossimMapProjection::forward(), ossimImageGeometry::getAsMapProjection(), ossimImageGeometry::getMetersPerPixel(), ossimImageGeometry::localToWorld(), os2, ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

Referenced by configureGmljp2V2().

{
   if ( geom )
   {
      ossimRefPtr<const ossimMapProjection> mapProj = geom->getAsMapProjection();
      if ( mapProj.valid() != 0 )
      {
         // Get the tie point:
         ossimGpt gpt;
         geom->localToWorld( ul, gpt );

         ossimDpt tie = mapProj->forward( gpt );

         std::ostringstream os;
         os.precision(15);
         os << tie.x << " " << tie.y;
         originString = os.str();   
      }

      // Get the projected CS gsd:
      ossimDpt gsdMap;
      geom->getMetersPerPixel( gsdMap );

      std::ostringstream os2;
      os2.precision(15);
      os2 << gsdMap.x << " " << "0.0";
      offsetVector1String = os2.str();

      std::ostringstream os3;
      os3.precision(15);
      os3 << "0.0" << " " << -gsdMap.y;
      offsetVector2String = os3.str();
   }

} // End: ossimGmlSupportData::getMapBounds

getXmlDoc()

ossimRefPtr< ossimXmlDocument > ossimGmlSupportData::getXmlDoc

(

)

const

Returns

Reference to xml doc.

Definition at line 539 of file ossimGmlSupportData.cpp.

References m_xmlDocument.

{
   return m_xmlDocument;
}

initialize() [1/2]

bool ossimGmlSupportData::initialize	(	const ossimImageGeometry *	geom,
		const ossimIrect &	rect
	)

Initializes gml block from geometry file.

Parameters

geom
rect	Output rectangle (view rect).

Definition at line 78 of file ossimGmlSupportData.cpp.

References configureGmljp2V2(), ossimImageGeometry::getAsMapProjection(), getGmljp2V1RootNode(), getGmljp2V2RootNode(), ossimMapProjection::getPcsCode(), ossimXmlDocument::initRoot(), m_pcsCodeGeo, m_pcsCodeMap, m_srsNameStringGeo, m_srsNameStringMap, m_xmlDocument, ossimFilename::NIL, os2, ossimNotify(), ossimNotifyLevel_DEBUG, status, and ossimRefPtr< T >::valid().

Referenced by ossimOpjCompressor::writeGmlBox(), and ossimKakaduCompressor::writeGmlBox().

{
   static const char MODULE[] = "ossimGmlSupportData::initialize(geom)";
   
   bool status = false;
   
   if ( geom )
   {
      // ossim_uint32 code = UNKNOWN_PCSCODE; // unknown code
       
      ossimRefPtr<const ossimMapProjection> mapProj = geom->getAsMapProjection();
      if ( mapProj.valid() )
      {
         // Get the PCS code:
         m_pcsCodeMap = mapProj->getPcsCode();
         m_pcsCodeGeo = mapProj->getPcsCode();
         
         // Create an SRS Name for the map projection
         std::ostringstream os;
         os << "http://www.opengis.net/def/crs/EPSG/0/" << m_pcsCodeMap;
         m_srsNameStringMap = os.str();
         
         // Create an SRS Name for the projection datum
         std::ostringstream os2;
         os2 << "http://www.opengis.net/def/crs/EPSG/0/" << m_pcsCodeGeo;
         m_srsNameStringGeo = os2.str();
         
         m_xmlDocument = new ossimXmlDocument(ossimFilename::NIL);
         
         // if ( m_use_gmljp2_version2 == true )
         // {

         ossimRefPtr<ossimXmlNode> rootNode = getGmljp2V2RootNode();
         m_xmlDocument->initRoot( rootNode );
         status = configureGmljp2V2( rootNode, geom, rect );
         // }
         
#if 0
         else
         {
            ossimRefPtr<ossimXmlNode> rootNode = getGmljp2V1RootNode();
            m_xmlDocument->initRoot( rootNode );
            status = configureGmljp2V1( rootNode, geom, rect );
         }
#endif
         
         // cout << "gmljp2Node: " << *(gmljp2Node.get()) << endl;
         
         // cout << "xmlDoc: " << *(m_xmlDocument.get()) << endl;
      }
   }

   if ( status == false )
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << MODULE << " DEBUG Entered...\n";
      
   }

   return status;
   
} // End: ossimGmlSupportData::initialize( geom, mapProj )

initialize() [2/2]

bool ossimGmlSupportData::initialize

(

std::istream &

in

)

Initializes from stream assumed to be a gml block.

Parameters

in

Definition at line 60 of file ossimGmlSupportData.cpp.

References m_xmlDocument, ossimXmlDocument::read(), and status.

{
   bool status = false;
   
   if ( in.good() )
   {
      m_xmlDocument = new ossimXmlDocument();
      m_xmlDocument->read( in );
      status = true;
   }
   else
   {
      m_xmlDocument = 0;
   }
   
   return status;
}

write()

bool ossimGmlSupportData::write

(

std::ostream &

os

)

Write xml doc to stream.

Definition at line 525 of file ossimGmlSupportData.cpp.

References ossimRefPtr< T >::get(), m_xmlDocument, status, and ossimRefPtr< T >::valid().

Referenced by ossimOpjCompressor::writeGmlBox(), and ossimKakaduCompressor::writeGmlBox().

{
   bool status = false;
   
   if ( m_xmlDocument.valid() )
   {
      os << *(m_xmlDocument.get());
      
      status = true;
   }

   return status;
}

Member Data Documentation

m_axisLabelsStringGeo

ossimString ossimGmlSupportData::m_axisLabelsStringGeo

private

Definition at line 164 of file ossimGmlSupportData.h.

Referenced by configureBounds(), and configureGmljp2V2().

m_axisLabelsStringMap

ossimString ossimGmlSupportData::m_axisLabelsStringMap

private

Definition at line 163 of file ossimGmlSupportData.h.

Referenced by configureGmljp2V2().

m_pcsCodeGeo

ossim_uint32 ossimGmlSupportData::m_pcsCodeGeo

private

Definition at line 159 of file ossimGmlSupportData.h.

Referenced by initialize().

m_pcsCodeMap

ossim_uint32 ossimGmlSupportData::m_pcsCodeMap

private

Definition at line 158 of file ossimGmlSupportData.h.

Referenced by initialize().

m_srsDimensionString

ossimString ossimGmlSupportData::m_srsDimensionString

private

Definition at line 162 of file ossimGmlSupportData.h.

Referenced by configureBounds().

m_srsNameStringGeo

ossimString ossimGmlSupportData::m_srsNameStringGeo

private

Definition at line 161 of file ossimGmlSupportData.h.

Referenced by configureBounds(), configureGmljp2V2(), and initialize().

m_srsNameStringMap

ossimString ossimGmlSupportData::m_srsNameStringMap

private

Definition at line 160 of file ossimGmlSupportData.h.

Referenced by configureGmljp2V2(), and initialize().

m_uomLabelsStringGeo

ossimString ossimGmlSupportData::m_uomLabelsStringGeo

private

Definition at line 166 of file ossimGmlSupportData.h.

Referenced by configureBounds().

m_uomLabelsStringMap

ossimString ossimGmlSupportData::m_uomLabelsStringMap

private

Definition at line 165 of file ossimGmlSupportData.h.

m_xmlDocument

ossimRefPtr<ossimXmlDocument> ossimGmlSupportData::m_xmlDocument

private

Definition at line 155 of file ossimGmlSupportData.h.

Referenced by addLineSamps(), addScale(), addTie(), getImageGeometryFromRectifiedGrid(), getImageGeometryFromSeonsorModel(), getXmlDoc(), initialize(), and write().

---

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

• ossimGmlSupportData.h
• ossimGmlSupportData.cpp