[发明专利]一种带补偿策略的火山口形SAR图像边缘检测方法

权利要求书

1.一种带补偿策略的火山口形SAR图像边缘检测方法，其特征在于，包括以下步骤：

步骤1，确定SAR图像，所述SAR图像的大小为a×b，所述SAR图像中每一个像素点分别对应通过一个火山口形滑动窗口滤波器的窗函数进行滤波，进而计算得到SAR图像中每一个像素点处的边缘响应值；其中，a表示SAR图像的长，b表示SAR图像的宽；

步骤2，根据SAR图像中每一个像素点处的边缘响应值，确定强候选边缘图和待确定边缘图；

步骤3，计算待确定边缘图中每一个待确定像素点的比例响应强度值；

步骤4，根据待确定边缘图中每一个待确定像素点的比例响应强度值，得到B'个弱边缘像素点和B”个非边缘像素点；B'和B”分别为大于或等于0的正整数；

步骤5，计算得到补偿后的B'个弱边缘像素点；

步骤6，根据强候选边缘图、补偿后的B'个弱边缘像素点和B”个非边缘像素点，得到SAR图像的边缘强度图；

步骤7，根据所述SAR图像的边缘强度图，得到SAR图像边缘图，所述SAR图像边缘图为一种带补偿策略的火山口形SAR图像边缘检测结果；其中，所述SAR图像边缘图的大小与所述SAR图像的大小相同。

2.如权利要求1所述的一种带补偿策略的火山口形SAR图像边缘检测方法，其特征在于，步骤1的子步骤为：

1.1计算SAR图像中每一个像素点通过火山口形滑动窗口滤波器的窗函数滤波后的计算值，其中像素点(x',y')通过对应火山口形滑动窗口滤波器的窗函数滤波后的计算值为W_CSW(x',y')，其表达式为：

其中，1≤x'≤a，1≤y'≤b，a表示SAR图像的长，b表示SAR图像的宽，Γ(·)表示伽马函数，α表示第一可调参数，α>1；β表示第二可调参数，β>0；然后将x'的值初始化为1，将y'的值初始化为1；

1.2将像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数平均划分为k个互不重叠的子窗口，记为w₁,…,w_k；

1.3计算得到像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数滤波到第i个子窗口w_i的局部加权均值为m_i(x',y')，其计算公式为：

其中，第i个子窗口w_i的大小为c×d，c表示第i个子窗口w_i的长，d表示第i个子窗口w_i的宽，0<c<a，0<d<b；第i个子窗口包括c×d个像素点，W_CSW(x_i,y_i)表示像素点(x_i,y_i)通过对应火山口形滑动窗口滤波器的窗函数滤波后的计算值，1≤x_i≤c，1≤y_i≤d；将I(x',y')记为SAR图像中像素点(x',y')处的灰度值，I(x'-x_i,y'-y_i)表示SAR图像中像素点(x'-x_i,y'-y_i)处的灰度值，S_CSW表示a×b个火山口形滑动窗口滤波器的窗函数横截面积，-∞<x”<+∞，-∞<y”<+∞，d·表示·的微分；

1.4令i的值分别从1取至k，重复1.3，进而分别得到像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数滤波到第1个子窗口w₁的局部加权均值m₁(x',y')至像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数滤波到第k个子窗口w_k的局部加权均值m_k(x',y')，然后计算得到像素点(x',y')处的边缘响应值ECM_CSW(x',y')，其表达式为：

其中，j＝(i+1)mod k，mod表示为取模计算，m_i(x',y')表示像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数滤波到第i个子窗口w_i的局部加权均值，m_j(x',y')像素点(x',y')对应的火山口形滑动窗口滤波器的窗函数滤波到第j个子窗口w_j的局部加权均值；然后将i的值初始化为1；

1.5令1≤x'≤a，1≤y'≤b，重复1.3和1.4遍历SAR图像中a×b个像素点，直到得到SAR图像中每一个像素点处的边缘响应值，即得到a×b个像素点处的边缘响应值。