Hello I am trying to patch process a folder full of ASCs to contours using Gdal from python. here is my code;
import os import glob PATH1 = raw_input("folder containing ASCs path please:") inputlistB =  for filename in glob.glob(os.path.join(PATH1, '*h_Max.ASC')): print filename inputlistB.append(filename) print inputlistB for filepath in inputlistB: src = filepath dst = filepath+'b.shp' os.system('gdalcontour -fl -888 -inodata -f "ESRI Shapefile "' + src + " " + dst)
The ipython console returns this;
folder containing ASCs path please:C:\D\OUT
showing that the inputlist functions as intended. However no new gis layers are generated in my target folder. If I take the statement ran by os.system and copy it with file paths into command line, it works. Though it does crash after producing the shapefile. If I just run the final line of the for loop os.system(... the Ipython console simply returns :1
Any idea why gdal would crash after producing a shapefile / successfully running?
Any ideas as to why my code does not produce any results despite no python errors?
Any alternative / better solutions to this problem?
also I use the file path input and loop in a lot of programs, feel free to steal it if this is news to you.