[Ubuntu] bug in YOUR version of listgeo

[Heikki Julkunen]

Btw, is there a way of getting more precision (like in the currently buggy listgeo) but make sure to include a separator ?

I'm seeing some rounding errors after doing some image processing in non gdal tools (and extracting geo data using listgeo, 
then adding it back geotifcp). With only 4 digits of precision there are slight differences from the original geo data.


On Apr 20, 2012, at 9:38 AM, Frank Warmerdam wrote:

> Folks,
> 
> I'll note that this bug is fixed in libgeotiff 1.4.0 though I notice
> that I got stuck at
> the RC level again.  I'll try to push 1.4.0 final today.
> 
> Note that libgeotiff 1.4.0 is binary incompatible with 1.3.x due to
> changes in the
> layout of some public structures.
> 
> Best regards,
> Frank
> 
> On Fri, Apr 20, 2012 at 8:44 AM, Heikki Julkunen
> <heikki.julkunen at gmail.com> wrote:
>> Guys,
>> listgeo in your unstable branch prints way to many digits so there is never a separator.
>> geotifcp then can't parse that number...
>> 
>> Please fix.
>> Also, when are you going to update to gdal1.9 ?
>> Geo-pdf is supposedly included in later versions of gdal.
>> 
>> Thanks
>> / Heikki
>> 
>> Example, check out the 2nd line in the ModelTiepointTag
>>>         -165510.298      71891.3101       0                << Ubuntu 1.2.5-1 version
>> ---
>> <          -165510.29779975571891.3100748502 0               <<< NO WHITESPACE!!
>> 
>> Entire file generated by Ubuntu 1.2.5-1 listgeo:
>> Geotiff_Information:
>>  Version: 1
>>  Key_Revision: 1.0
>>  Tagged_Information:
>>     ModelTiepointTag (2,3):
>>        0                0                0
>>        -165510.298      71891.3101       0
>>     ModelPixelScaleTag (1,3):
>>        31.735961        31.735961        0
>>     End_Of_Tags.
>>  Keyed_Information:
>>     GTModelTypeGeoKey (Short,1): ModelTypeProjected
>>     GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
>>     GTCitationGeoKey (Ascii,225): "IMAGINE GeoTIFF Support\nCopyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved\n@(#)$RCSfile: egtf.c $ $Revision: 1.9.2.4 $ $Date: 2002/01/24 16:43:13Z $\nProjection Name = San Diego TAC\nUnits = meters\nGeoTIFF Units = meters"
>>     GeographicTypeGeoKey (Short,1): GCS_NAD83
>>     ProjectedCSTypeGeoKey (Short,1): User-Defined
>>     PCSCitationGeoKey (Ascii,192): "IMAGINE GeoTIFF Support\nCopyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved\n@(#)$RCSfile: egtf.c $ $Revision: 1.9.2.4 $ $Date: 2002/01/24 16:43:13Z $\nProjection = Lambert Conformal Conic"
>>     ProjectionGeoKey (Short,1): User-Defined
>>     ProjCoordTransGeoKey (Short,1): CT_LambertConfConic_2SP
>>     ProjLinearUnitsGeoKey (Short,1): Linear_Meter
>>     ProjStdParallel1GeoKey (Double,1): 45
>>     ProjStdParallel2GeoKey (Double,1): 33
>>     ProjNatOriginLongGeoKey (Double,1): -117.25
>>     ProjFalseOriginLongGeoKey (Double,1): 0
>>     ProjFalseOriginLatGeoKey (Double,1): 33
>>     ProjFalseOriginEastingGeoKey (Double,1): 0
>>     ProjFalseOriginNorthingGeoKey (Double,1): 0
>>     End_Of_Keys.
>>  End_Of_Geotiff.
>> 
>> Projection Method: CT_LambertConfConic_2SP
>>  ProjFalseOriginLatGeoKey: 33.000000 ( 33d 0' 0.00"N)
>>  ProjFalseOriginLongGeoKey: -117.250000 (117d15' 0.00"W)
>>  ProjStdParallel1GeoKey: 45.000000 ( 45d 0' 0.00"N)
>>  ProjStdParallel2GeoKey: 33.000000 ( 33d 0' 0.00"N)
>>  ProjFalseEastingGeoKey: 0.000000 m
>>  ProjFalseNorthingGeoKey: 0.000000 m
>> GCS: 4269/NAD83
>> Datum: 6269/North American Datum 1983
>> Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
>> Prime Meridian: 8901/Greenwich (0.000000/  0d 0' 0.00"E)
>> Projection Linear Units: 9001/metre (1.000000m)
>> 
>> PROJ.4 Definition: +proj=lcc +lat_0=33.000000000 +lon_0=-117.250000000 +lat_1=45.000000000 +lat_2=33.000000000  +x_0=0.000 +y_0=0.000 +ellps=GRS80 +units=m
>> 
>> Corner Coordinates:
>> Upper Left    ( -165510.298,   71891.310)  (119d 2' 9.39"W, 33d38' 1.95"N)
>> Lower Left    ( -165510.298,  -59051.265)  (119d 0'30.96"W, 32d27'12.09"N)
>> Upper Right   (   90567.172,   71891.310)  (116d16'21.53"W, 33d38'38.97"N)
>> Lower Right   (   90567.172,  -59051.265)  (116d17'15.40"W, 32d27'48.47"N)
>> Center        (  -37471.563,    6420.023)  (117d39' 4.56"W, 33d 3'25.72"N)
>> 
>> _______________________________________________
>> UbuntuGIS mailing list
>> Ubuntu at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/ubuntu
>> http://trac.osgeo.org/ubuntugis/wiki
> 
> 
> 
> -- 
> ---------------------------------------+--------------------------------------
> I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush    | Geospatial Software Developer