def index(self, x, y):
        col, row = [math.floor(a) for a in (~self.affine * (x, y))]
        return row, col
Given x,y in crs, return the (column, row) on the raster
dataset/ETHPy150Open perrygeo/python-rasterstats/src/rasterstats/io.py/Raster.index