Personalmente, realmente disfruto codificando con GDAL en Python. Una vez que lo tienes instalado y sabes cómo funciona, es bastante fácil. También con python puede utilizar algunos módulos incorporados realmente útiles 'baterías incluidas' (sistema operativo, matemática, cadena, sys, aleatorio, etc.), así como módulos de terceros como NumPy (útil para operaciones ráster) entre otros.
Para obtener una lista completa de las bibliotecas Python de SIG, consulte aquí: ¿Cuáles son las herramientas / módulos / complementos de Python cruciales en SIG? )
Una excelente introducción a los enlaces Python GDAL se puede encontrar aquí:
http://www.gis.usu.edu/~chrisg/python/2009/
Puede comenzar en un tiempo mínimo con un poco de conocimiento y lectura, por ejemplo, puede agarrar un shapefile en 2 líneas:
driver = ogr.GetDriverByName('ESRI Shapefile')
datasource = driver.Open('shapefile.shp')
Es muy fácil hacer lo que quiera al archivo una vez que haya leído la documentación ( http://www.gdal.org/annotated.html ) y sepa lo que está haciendo.
Vea aquí algunos ejemplos de lo que es posible en código simple y mínimo:
http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/