Лично мне очень нравится программировать с помощью GDAL на Python. После того, как вы установите его и разберетесь в том, как он работает, все будет довольно просто. Также с python вы можете использовать некоторые действительно полезные встроенные модули «с батарейками» (os, math, string, sys, random и т. Д.), А также сторонние модули, такие как NumPy (полезно для растровых операций).
Полный список библиотек Python GIS см. Здесь: Какие инструменты / модули / надстройки Python имеют решающее значение в ГИС? ).
Действительно отличное введение в привязки Python GDAL можно найти здесь:
http://www.gis.usu.edu/~chrisg/python/2009/
Вы можете приступить к работе за минимальное время, немного научившись и читая, например, вы можете взять шейп-файл в 2 строки:
driver = ogr.GetDriverByName('ESRI Shapefile')
datasource = driver.Open('shapefile.shp')
Очень легко сделать с файлом все, что вы хотите, прочитав документацию ( http://www.gdal.org/annotated.html ) и узнав, что вы делаете.
Смотрите здесь некоторые примеры того, что возможно в простом и минимальном коде:
http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/