[gdal-dev] Calling GDAL directly from within Python

[Siebe Bosch]

Dear list,

Below you will find a question from my intern Bram. Unfortunately I could not help him out either since I only use the gdal executables.
We would highly appreciate your help calling GDAL functions from within Python.

Kind regards,
Siebe Bosch


Hi devs,

This is my first time using the Gdal library directly in Python using the Anaconda Spyder shell. I want to write a program that automatically imports GRIB files from a public ftp server and use GDAL to perform some actions. I installed GDAL with the PIP installer. Python version 2.7.13, Anaconda 32-bit, GDAL 3.0.8.
First off it worked well; GDAL opened the file and could read its information (as in raster size and coordinates (not taking the automatically importing from the server into account)). I did noticed the difference in commands (for example, gdal_Translate turns into gdal.Translate(src_ds, DestName, **kwargs)):

Import gdal
 From osgeo import osr
src_ds = gdal.Open('C:\Users\wijna_000\H11_V72_NED_SINGLE_LEVEL_00_000_GB')

raster_wkt = src_ds.GetProjection()
spatial_ref = osr.SpatialReference()
spatial_ref.ImportFromWkt(raster_wkt)
print spatial_ref.ExportToPrettyWkt()

 From here I want to transform the GRIB file into ,preferably ASCII files, but tried .tif first since it’s the default. But when I run the gdal.Translate command and do the coordinate check again it still returns as data imported from a grib-file:

tif_ds = gdal.Translate("C:\Users\wijna_000\output.tif", src_ds)

raster_wkt2 = tif_ds.GetProjection()
spatial_ref2 = osr.SpatialReference()
spatial_ref2.ImportFromWkt(raster_wkt2)
print spatial_ref2.ExportToPrettyWkt()

Return:
GEOGCS["Coordinate System imported from GRIB file", DATUM["unknown",SPHEROID["Sphere",6367470,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]

Help on internet is rare since its within the Python environment. Can anyone point me towards a better direction?

Regards, Bram


-- 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20170613/87d7b20a/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature_60pct.gif
Type: image/gif
Size: 30165 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20170613/87d7b20a/attachment-0001.gif>