I am trying to carry out the task of sending each row/polygon from a single shapefile to its own unique raster. I so far use GDAL commands through the OSGeo4W Shell like so after navigating to my directory and entering
FOR /L %i IN (1,1,250) DO gdal_rasterize -a Value -l Polygons -where "Polygon_ID=%i" -tr 100 100 Polygons.shp Polygon_raster%i.tif
What this does is loop through all the polygons (rows) in my
Polygons.shp shapefile and send them to their own unique raster, burning in the “Value” values. While this method is nice and straightforward, it does not let me do more complex tasks like multiplying each output raster by a specific value. What I specifically want to do is this:
Each polygon (row) in the shapefile contains both a “Value” column and a “Value2” column. I want to create all of my individual rasters like above, but I want to multiply each raster by the associated “Value2” value for that specific polygon. I do not know how to do this through OSGeo4W command line statements, and would much prefer doing this in python with Jupyter Notebook. However, I am having a lot of trouble directly translating my above gdal command statement into code that would work in Jupyter Notebook with python, since my command statement above uses a Windows for loop.
I would try this:
from osgeo import gdal shapefile = polygons.shp Value2 = shapefile[‘Value2’] for polygon in shapefile: raster = gdal_rasterize(‘Value’) raster * Value2 write_shapefile(raster, 'Polygon_raster%i.tif')
I want to do this in Jupyter Notebook using Python, but I am unsure how to directly translate my above command line into the necessary Python code that would work here, and particularly so that each raster would receive its own unique ID number (%i), which is straightforward with the command statement. How could I translate my command line statement into Python code that would let me loop through the rows in my shapefile, rasterize each, and multiply by its corresponding “Value2” value?
gdal.Rasterize()as opposed to running
gdal_rasterizein a command line. Yes I have tried looking through this "GDAL/OGR Python API" resource here: gdal.org/python