矩阵是离散的,所以:
def gradx(Img, sx = 1.0): return sx * (Img[1: , : ] - Img[ :-1, : ])
def grady(Img, sy = 1.0): return sy * (Img[ : , 1: ] - Img[ : , : -1])