title: Find indices of raster cells that intersect with a polygon
body: <p>I want to get a list of indices (row,col) for all raster cells that fall within or are intersected by a polygon feature. Looking for a solution in python, ideally with gdal/ogr modules.</p>
  
  <p>Other posts have suggested rasterizing the polygon, but I would rather have direct access to the cell indices if possible.</p>