This is my first question in the forum, so please apologize if I’m missing to follow some rules.

I’m trying to run gdal_polygonize.py to polygonize a raster file with this code,

print("# 4- CONVERT RASTER DATA TO VECTOR")
input_layer = output
print("input_layer: ", input_layer)
output = vec_3857 + "BUILT_%s_%s_3857.shp"%(fua,date)
print("output to be generated: ", output)
name = "BUILT_%s_%s_3857"%(fua,date)
subprocess.call('"C:\OSGeo4W\apps\gdal\share\bash-completion\completions\gdal_polygonize.py" %s -f "ESRI Shapefile" %s' %(input_layer, output),shell=True)

The input_layer is already where the code is looking for it. But I’m not understanding the reason for that gdal_polygonize.py is not generating the output shapefile. There’s no error message indicating something wrong when running subprocess.call(...).

I tried to address the issue with solutions provided in this entry, /a/40157820/17237490, but I was unsuccessful. I think something

Any help is much appreciated! Thank you.