Re: [OSSIM] GDAL geometry

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 454-5900

I have a problem with an Erdas .img file. 
The output of gdalinfo for this image is at the end of this message. 
when I load the image using the gdal plug in and ask for the geometry, 
the function 

bool ossimGdalTileSource::getImageGeometry(ossimKeywordlist& kwl,
                                           const char* prefix)

exits at this point and I don't understand why. 
   if(ossimImageHandler::getImageGeometry(kwl, prefix))
   {
      return true;
   }

I have commented this part and I have another problem. the function read the wkt from gdal successfully
but It fails on 
   bool transOk = GDALGetGeoTransform( theDataset, geoTransform ) == CE_None;
so the resultant geometry has the projection but not the transform. 

if I comment this line
const char* wkt = GDALGetProjectionRef( theDataset );

the algorithm successfully estimates a ossimBilinearProjection

what I think the algorithm should do is to estimate geoTransform from the GCP if it fails on getting from GDALGetGeoTransform. what do you think?
also do you know why is exiting on the second line ?

Thanks in advance
Jorge Artieda

Driver: HFA/Erdas Imagine Images (.img). 
Files: c:\datos\cse05\ms\l71201032_03220040630_hrf.img
       c:\datos\cse05\ms\l71201032_03220040630_hrf.rrd
Size is 7257, 6221
Coordinate System is:
PROJCS["UTM Zone 30, Northern Hemisphere",
    GEOGCS["Undefined",
        DATUM["Undefined",
            SPHEROID["WGS 84",6378137,298.2572235629972]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-3],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GCP Projection = PROJCS["UTM Zone 30, Northern Hemisphere",GEOGCS["Undefined",DATUM["Undefined",SPHEROID["WGS 84",6378137,298.2572235629972]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]
GCP[  0]: Id=, Info=
          (0.5,0.5) -> (340787.495142927,4577896.72550787,0)
GCP[  1]: Id=, Info=
          (1451.7,0.5) -> (383478.748242927,4569362.07550787,0)
GCP[  2]: Id=, Info=
          (2902.9,0.5) -> (426170.001342927,4560827.42550787,0)
GCP[  3]: Id=, Info=
          (4354.1,0.5) -> (468861.254442927,4552292.77550787,0)
GCP[  4]: Id=, Info=
          (5805.3,0.5) -> (511552.507542927,4543758.12550787,0)
.................. more gcp here
GCP[ 35]: Id=, Info=
          (7256.5,6220.5) -> (517663.245142927,4352244.22550787,0)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 6221.0)
Upper Right ( 7257.0,    0.0)
Lower Right ( 7257.0, 6221.0)
Center      ( 3628.5, 3110.5)
Band 1 Block=64x64 Type=Byte, ColorInterp=Undefined
  Description = LANDSAT7_band1
  Min=0.000 Max=255.000 
  Minimum=0.000, Maximum=255.000, Mean=68.774, StdDev=47.123
  Overviews: 1815x1556, 908x778, 454x389, 227x195, 114x98, 57x49
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=68.773525144508
    STATISTICS_MEDIAN=86
    STATISTICS_MODE=0
    STATISTICS_STDDEV=47.123081608966
    STATISTICS_HISTONUMBINS=256
    STATISTICS_HISTOMIN=0
    STATISTICS_HISTOMAX=255
    LAYER_TYPE=athematic
    STATISTICS_HISTOBINVALUES=13351680|0|0|0|0|0|72|828|3924|5976|4464|1728|756|324|144|144|144|108|180|108|144|180|144|108|144|36|36|108|144|144|180|72|216|108|180|72|0|72|108|180|144|108|108|144|108|72|144|180|144|72|36|216|1

----- Original Message ----
From: Garrett Potts <po...@cf...>
To: Ossim users <oss...@li...>
Sent: Monday, June 23, 2008 2:50:31 PM
Subject: Re: [OSSIM] GDAL geometry

Hello:

We do support GDAL geometries.  I have a translation table in the gdal  
plugin as examples.  Just note,  I might need your help for more  
translations as you see them come up.  The code is located in:

ossim_plugins/gdal/ossimOgcWktTranslator.cpp

Note:  if you put the .type as the WKT string it should create the  
projection but,  for image models that is not enough for you still  
need tie point and pixel scale meters or decimal degrees depending on  
the input projection type.

Take care

Garrett

On Jun 23, 2008, at 6:40 AM, j artieda wrote:

> Are GDAL geometries and projections supported by OSSIM?
> Now I use gdalinfo to read metadata and write my .geom file is there  
> an easier way.
> I am developing a program which should read from several formats and  
> projections I would like to avoid writing all the particular  
> situations.
> If gdal geometries are not supported which is in your opinion the  
> best way to introduce it in OSSIM?
> thanks in advance
> Jorge Artieda
>
>
>
>
> -------------------------------------------------------------------------
> Check out the new SourceForge.net Marketplace.
> It's the best place to buy or sell services for
> just about anything Open Source.
> http://sourceforge.net/services/buy/index.php
> _______________________________________________
> www.ossim.org
> Ossim-developer mailing list
> Oss...@li...
> https://lists.sourceforge.net/lists/listinfo/ossim-developer

-------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services for
just about anything Open Source.
http://sourceforge.net/services/buy/index.php
_______________________________________________
www.ossim.org
Ossim-developer mailing list
Oss...@li...
https://lists.sourceforge.net/lists/listinfo/ossim-developer