[gdal-dev] Difference in output using gdal_translate from FWTools shell and cmd window

Anne Brookes Anne.Brookes at jbaconsulting.com
Thu Dec 1 06:21:09 EST 2011


Hi,
Thx for your reply.  I'm not entirely sure what you mean about access to the proj and epsg libraries.  proj does seem to be a command that can be run from any cmd window.
I have checked for the existence of pcs.csv and there seem to be numerous versions (also pcs.override.csv) in various locations, including my FWTools2.0.4\data; under my GRASS6.4 installation folders and under ArcGIS installation folders.  I have tried running the gdalinfo commands both with a GDAL_DATA environment variable (pointing to FWTools2.0.4\data) and without this env variable but it didn't seem to make any difference.

What proj and epsg libraries are used by gdal and where should they be?
Thanks,
Anne

Anne Brookes
Senior Analyst





http://www.jbaconsulting.com/
South Barn
Broughton Hall
Skipton
North Yorkshire
BD23 3AE
T: +44 (0)1756 799919 F: 0845 8627772
The JBA Family
http://www.jbagroup.co.uk/ http://www.jbarisk.com/ http://www.jbaenergy.com/ http://www.jbatrust.org/ http://www.jbaconsulting.com/

This email is covered by JBA Consulting's http://www.jbaconsulting.com/disclaimer JBA Consulting is the trading name of Jeremy Benn Associates Limited, registered in England, company number 03246693, South Barn, Broughton Hall, Skipton, North Yorkshire, BD23 3AE.
http://www.jbaconsulting.com/nce-winner-2010  http://twitter.com/JBAConsulting





From: gdal-dev-bounces at lists.osgeo.org [mailto:gdal-dev-bounces at lists.osgeo.org] On Behalf Of Frank Broniewski
Sent: 01 December 2011 11:00
To: gdal-dev at lists.osgeo.org
Subject: Re: [gdal-dev] Difference in output using gdal_translate from FWTools shell and cmd window

Hi,

are you sure that the same libraries (proj, epsg et. al.) are available
under both cirumstances? Do you have access to proj and epsg codes in
plain cmd?
That's where I would start to investigate ...

btw. I think the osgeo4w package is more recent than the fwtools build,
so might also have a look at that one


Frank

Am 01.12.2011 11:36, schrieb Anne Brookes:
> Hello,
> When I run gdalinfo and gdal_translate commands from the FWTools Shell (on Windows) the output is different to when I run the same commands from an ordinary cmd window (or from a python script). The difference is in the PROJCS, DATUM and SPHEROID metadata.
>
> Here is the (correct) output from the FWTools shell command window:
> C:\Program Files\FWTools2.0.4>gdalinfo : C:\JFLOW8\testgdal\D108787.tif
>
> Driver: GTiff/GeoTIFF
> Files: C:\JFLOW8\testgdal\D108787.tif
> Size is 597, 597
> Coordinate System is:
> PROJCS["ETRS89 / ETRS-LAEA",
>      GEOGCS["ETRS89",
>          DATUM["European_Terrestrial_Reference_System_1989",
>              SPHEROID["GRS 1980",6378137,298.2572221010002,
>                  AUTHORITY["EPSG","7019"]],
>              AUTHORITY["EPSG","6258"]],
>          PRIMEM["Greenwich",0],
>          UNIT["degree",0.0174532925199433],
>          AUTHORITY["EPSG","4258"]],
>      UNIT["metre",1,
>          AUTHORITY["EPSG","9001"]],
>      AUTHORITY["EPSG","3035"]]
> Origin = (4039990.000000000000000,2950770.000000000000000)
> Pixel Size = (10.000000000000000,-10.000000000000000)
> Metadata:
>    AREA_OR_POINT=Area
> Image Structure Metadata:
>    INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 4039990.000, 2950770.000)
> Lower Left  ( 4039990.000, 2944800.000)
> Upper Right ( 4045960.000, 2950770.000)
> Lower Right ( 4045960.000, 2944800.000)
> Center      ( 4042975.000, 2947785.000)
> Band 1 Block=597x3 Type=Float32, ColorInterp=Gray
>    NoData Value=-9999
>
>
> And here is the output from the cmd window :
> C:\Document and Settings\jflownw>  " C:\Program Files\FWTools2.0.4\bin\gdalinfo" C:\JFLOW8\testgdal\D108787.tif
>
> Driver: GTiff/GeoTIFF
> Files: C:\JFLOW8\testgdal\D108787.tif
> Size is 597, 597
> Coordinate System is:
> PROJCS["unnamed",
>      GEOGCS["ETRS89",
>          DATUM["unknown",
>              SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
>          PRIMEM["Greenwich",0],
>          UNIT["degree",0.0174532925199433]],
>      UNIT["metre",1,
>          AUTHORITY["EPSG","9001"]],
>      AUTHORITY["EPSG","3035"]]
> Origin = (4039990.000000000000000,2950770.000000000000000)
> Pixel Size = (10.000000000000000,-10.000000000000000)
> Metadata:
>    AREA_OR_POINT=Area
> Image Structure Metadata:
>    INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 4039990.000, 2950770.000)
> Lower Left  ( 4039990.000, 2944800.000)
> Upper Right ( 4045960.000, 2950770.000)
> Lower Right ( 4045960.000, 2944800.000)
> Center      ( 4042975.000, 2947785.000)
> Band 1 Block=597x3 Type=Float32, ColorInterp=Gray
>    NoData Value=-9999
>
> Please can anyone tell me why this is happening and what I need to do to be able to get the full DATUM and SPHEROID metadata when running the GDAL commands from outside of the FWTools shell, eg. from a python script or another cmd window?
> I have been puzzling over this for a while and it is now becoming quite urgent.
> Thanks
> Anne Brookes
> Senior Analyst
>
>
>
>
>
>
>
> [cid:imagea52157.png at 0b641363.58fb4337]<http://www.jbaconsulting.com/>
>
>
> South Barn
> Broughton Hall
> Skipton
> North Yorkshire
> BD23 3AE
> T: +44 (0)1756 799919 F: 0845 8627772
>
>
> The JBA Family
>
> [cid:imageb7288b.png at 04236f77.d2b04ae4]<http://www.jbagroup.co.uk/>   [cid:image840d25.png at a2776b7f.c5f34cb4]<http://www.jbarisk.com/>    [cid:image98f3a8.png at fecc7dc0.dc434408]<http://www.jbaenergy.com/>    [cid:image8ba6f0.png at b44205ef.33e84574]<http://www.jbatrust.org/>
> www.jbaconsulting.com<http://www.jbaconsulting.com/>
>
>
>
>
>
> This email is covered by JBA Consulting's email disclaimer<http://www.jbaconsulting.com/disclaimer>   JBA Consulting is the trading name of Jeremy Benn Associates Limited, registered in England, company number 03246693, South Barn, Broughton Hall, Skipton, North Yorkshire, BD23 3AE.
>
>
> [cid:image836b39.png at 13d0faf1.51b34be6]<http://www.jbaconsulting.com/nce-winner-2010>   [cid:imagea20f2e.png at e453a025.69b74420]<http://twitter.com/JBAConsulting>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev


--
Frank BRONIEWSKI

METRICO s.à r.l.
géomètres
technologies d'information géographique
rue des Romains 36
L-5433 NIEDERDONVEN

tél.: +352 26 74 94 - 28
fax.: +352 26 74 94 99
http://www.metrico.lu
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev



More information about the gdal-dev mailing list