Jak zaimplementować gradientową transformację Hougha

Próbuję użyć transformacji Hougha do wykrywania krawędzi i chciałbym wykorzystać obrazy gradientowe jako podstawę. Co zrobiłem do tej pory, ze względu na obraz Iwielkości [M,N]i jej pochodnych cząstkowych gx, gyjest obliczenie kąta gradientu w każdym pikselu jako thetas = atan(gy(x,y) ./ gx....

image-processing edge-detection image-processing computer-vision image-registration discrete-signals noise bpsk snr demodulation bpsk multipath synchronization timing image-processing filters algorithms edge-detection sampling demodulation bpsk synchronization timing fft fourier-transform delay audio speech-recognition soft-question discrete-signals discrete-signals autocorrelation frequency computer-vision