<div dir="ltr"><div class="gmail_quote"><div dir="ltr">Hello,<div><br></div><div>I am trying to convert a gdal_rasterize command to a gdal.Rasterize call with the Python bindings </div><div>(Python 2.7.9 64bits, GDAL 2.1.0 compiled with libproj0 libgeos-dev on Debian 8.3).</div><div><br></div><div>The  gdal_rasterize command works perfectly  :</div><div><div><div data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>gdal_rasterize -i -b 4 -burn 0 my.geojson my.tif</div><div><br></div><div>With the information available on the web :</div><div><div><a href="http://erouault.blogspot.ie/2015/10/gdal-and-ogr-utilities-as-library.html" target="_blank">http://erouault.blogspot.ie/2015/10/gdal-and-ogr-utilities-as-library.html</a></div><div><a href="https://trac.osgeo.org/gdal/wiki/rfc59.1_utilities_as_a_library" target="_blank">https://trac.osgeo.org/gdal/wiki/rfc59.1_utilities_as_a_library</a></div><div><a href="http://gis.stackexchange.com/questions/77139/why-do-shapefiles-and-geojson-behave-differently-in-gdal-python" target="_blank">http://gis.stackexchange.com/questions/77139/why-do-shapefiles-and-geojson-behave-differently-in-gdal-python</a></div></div><div><br></div><div>I've tried this kind of script :</div><div><br></div><div><div>import os</div><div>import sys</div><div>from osgeo import gdal, ogr, osr</div></div><div><br></div><div><div>geojson_filename = "my.geojson"</div><div>driver = ogr.GetDriverByName("GeoJSON")</div><div>vector_ds = driver.Open(geojson_filename, gdal.GA_ReadOnly)</div></div><div><br></div><div><div>geotiff_filename = "my.tif"</div><div>ds = gdal.Open(geotiff_filename, gdal.GA_Update)</div><div><br></div><div>ret = gdal.Rasterize(ds, vector_ds,</div><div>            bands = [4],</div><div>            inverse = True,</div><div>            burnValues = [0])</div><div>if ret != 1:</div><div>    print 'fail'</div><div>else:</div><div>    print 'success'</div><div><br></div><div>but I kept getting this error :</div><div><div>Traceback (most recent call last):</div><div>  File "test.py", line 62, in <module></div><div>    ret = gdal.Rasterize(ds, vector_ds, bands = [4], inverse = True, burnValues = [0])</div><div>  File "/usr/local/lib/python2.7/dist-packages/osgeo/gdal.py", line 1046, in Rasterize</div><div>    return wrapper_GDALRasterizeDestName(destNameOrDestDS, srcDS, opts, callback, callback_data)</div><div>TypeError: in method 'wrapper_GDALRasterizeDestName', argument 2 of type 'GDALDatasetShadow *'</div></div><div><br></div><div>Any ideas of what could be wrong with this script ?</div><div><br></div><div>Thanks in advance</div><div><br></div><div>Cheers</div><span class="HOEnZb"><font color="#888888"><div><br></div>François ALLAIN<br><br><br></font></span></div></div></div></div></div></div>
</div></div>
</div><br></div>