Canny edge的程式，我是依據此網頁改的 www.pages.drexel.edu/~nk752/cannyTut2.html#Step 1
程式執行時會發生 atan2:DOMAIN error 錯誤訊息!
發現是因為程式中

[code cpp]
thisAngle = (atan2(Gx,Gy)/3.14159) * 180.0; // Calculate actual direction of edge
[/code]

Gx與Gy同時為零，此時函式 atan2 就會發生錯誤。

假設一影像的灰階 p1 p2 p3
p4 p5 p6
p7 p8 p9

若 p1、p2、p3、p4、p6、p7、p8、p9灰階值皆相同，經sobel mask 後，Gx與Gy不就會為零嗎!
這樣我該如何計算 thisAngle ?

請問有人知道應該怎麼做嗎?

附加檔案：47d0eef0bbe9a_canny.rar

[code cpp]
if ( Gx == 0 && Gy == 0 )
thisAngle = 0;
else
thisAngle = (atan2(Gx,Gy)/3.14159) * 180.0;
[/code]