Как реализовать градиентное преобразование Хафа

Я пытаюсь использовать преобразование Хафа для обнаружения краев и хотел бы использовать градиентные изображения в качестве основы. То , что я сделал до сих пор, учитывая изображение Iразмера [M,N]и ее частных производных gx, gy, чтобы вычислить угол градиента в каждом пикселе , как thetas =...