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:[email protected]]<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:[email protected]]<http://www.jbagroup.co.uk/>
[cid:[email protected]] <http://www.jbarisk.com/>
[cid:[email protected]] <http://www.jbaenergy.com/>
[cid:[email protected]] <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:[email protected]]<http://www.jbaconsulting.com/nce-winner-2010>
[cid:[email protected]] <http://twitter.com/JBAConsulting>
<<inline: [email protected]>>
<<inline: [email protected]>>
<<inline: [email protected]>>
<<inline: [email protected]>>
<<inline: [email protected]>>
<<inline: [email protected]>>
<<inline: [email protected]>>
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
