Pessoalmente, eu realmente gosto de codificar com GDAL em Python. Depois de instalá-lo e entender como funciona, é bem fácil. Além disso, com o python, você pode usar alguns módulos internos realmente úteis (incluindo baterias) (os, math, string, sys, random etc), bem como módulos de terceiros, como o NumPy (útil para operações de varredura), entre outros.
Para obter uma lista abrangente das bibliotecas GIS Python, consulte aqui: Quais são as ferramentas / módulos / suplementos Python cruciais no GIS? )
Uma excelente introdução às ligações Python GDAL pode ser encontrada aqui:
http://www.gis.usu.edu/~chrisg/python/2009/
Você pode avançar em um tempo mínimo com um pouco de conhecimento e leitura, por exemplo, pode pegar um arquivo shapefile em duas linhas:
driver = ogr.GetDriverByName('ESRI Shapefile')
datasource = driver.Open('shapefile.shp')
É muito fácil fazer o que quiser com o arquivo depois de ler a documentação ( http://www.gdal.org/annotated.html ) e você saber o que está fazendo.
Veja aqui alguns exemplos do que é possível no código simples e mínimo:
http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/