gistfile1.py

def gradiente_sobel(punto, imagen):
    #print punto
    pixeles = imagen.load()
    x = punto[0]
    y = punto[1]
    #print x, y
    z1 = pixeles[x-1, y-1][0] 
    z2 = pixeles[x, y-1][0]
    z3 = pixeles[x+1, y-1][0]
    z4 = pixeles[x-1, y][0]
    z5 = pixeles[x, y][0]
    z6 = pixeles[x+1, y][0]
    z7 = pixeles[x-1, y+1][0]
    z8 = pixeles[x, y+1][0]
    z9 = pixeles[x+1, y+1][0]
    Gx = ((z3)+(2*z6)+z9)-((z1)+(2*z4)+(z7))
    Gy = ((z7)+(2*z8)+z9)-((z1)+(2*z2)+(z3))
    G = sqrt((Gx**2)+(Gy**2))
    angulo = atan2(Gy, Gx)
    angulo = angulo -(pi/2)
    angulo = float("%.2f" % angulo)
    return x,y, angulo