Standardní výstup můžete zachytit tímto způsobem:
from subprocess import Popen, PIPE
cmd = ['gdallocationinfo', '-geoloc', '-valonly', src, lng, lat]
p = Popen(cmd, stdout=PIPE)
p.wait()
val = p.stdout.read()
Trochu jsem hloubal a našel jsem odpověď na svou otázku:
místo použití os.system je správná syntaxe pro uložení výsledku do proměnné:
result = os.popen('gdallocationinfo -valonly -wgs84 %s %s' % (lyr, loc)).read()
Myslím, že budete muset použít syntaxi podprocesu, jak je vysvětleno v tomto příspěvku:
podproces pro volání gdal z pythonu
gpspipe vrací správné časové razítko, ale špatné datum (21. prosince 1995)
Jak nainstalovat GDAL a QGIS na Raspberry Pi?