技術特征：

1.一種邊緣檢測方法，其特征在于包括以下步驟：

步驟1：為了在去除噪聲的同時保留邊緣信息，對輸入圖像I₀進行基于各向異性熱擴散方程的預處理，得到預處理圖像I₁；

步驟2：在預處理圖像I₁中初始化活動輪廓C和迭代次數T₀；

步驟3：構造能量函數F，F由膨脹因子F_b、平滑因子F_s和收縮因子F_c三個部分組成，F＝F_b+F_s+F_c；

步驟4：迭代計算步驟2中活動輪廓C的所有像素點的能量函數F，直到迭代次數為T₀；

步驟5：將第T₀次迭代后的活動輪廓C_T作為邊緣輸出。

2.根據權利要求1所述的一種邊緣檢測方法，其特征在于所述的膨脹因子F_b的迭代表達式采用如下公式：

式(1)，為第t次迭代后的膨脹因子，為第t+1次迭代后的膨脹因子，為像素點(x,y)在第t次迭代后的灰度取值，sign為像素點(x,y)處的梯度方向，DIL_d為膨脹操作符，ERO_d為腐蝕操作符，T₁為灰度閾值。

3.根據權利要求1所述的一種邊緣檢測方法，其特征在于所述的平滑因子F_s的迭代表達式采用如下公式：

式(2)，為第t次迭代后的平滑因子，為第t+1次迭代后的平滑因子，符號為數學形態學的開運算操作符，SE₁為長度為三個像素點的直線結構元素集合，SE₂為包含于8鄰域的結構元素集合，T₂為灰度閾值。

4.根據權利要求1所述的一種邊緣檢測方法，其特征在于所述的收縮因子F_c的迭代表達式采用如下公式：

式(3)，為第t次迭代后的收縮因子，為第t+1次迭代后的收縮因子，T₃為灰度閾值。

5.根據權利要求1所述的一種邊緣檢測方法，其特征在于基于各向異性熱擴散方程的預處理方法采用如下公式：

$I_{x,y}^{t+1}=I_{x,y}^t+\mathrm{\λ}\left({\mathrm{cN}}_{x,y}{\▿}_N\right(I_{x,y}^t)+{\mathrm{cS}}_{x,y}{\▿}_S(I_{x,y}^t)+{\mathrm{cE}}_{x,y}{\▿}_E(I_{x,y}^t)+{\mathrm{cW}}_{x,y}{\▿}_W(I_{x,y}^t\left)\right)---\left(4\right)$

式(4)，為像素點(x,y)在第t次迭代后的灰度取值，為像素點(x,y)在第t+1次迭代后的灰度取值，k和λ均為平滑控制系數；和分別表示在上、下、左、右四個方向上對像素點(x,y)求偏導；

和分別表示上、下、左、右四個方向上的導熱系數。