[发明专利]地震剖面图像的纹理分割方法和装置

权利要求书

1.一种地震剖面图像的纹理分割方法，所述方法包括以下步骤：

输入原始的地震剖面图像；

将所述地震剖面图像转换为灰度图像；

利用小波变换对所述灰度图像进行特征提取；

对所述灰度图像进行图像纹理边缘检测，以产生剖面边缘图；以及

基于特征提取的结果和所述剖面边缘图应用聚类算法以按照纹理类型对所述地震剖面图像进行纹理分割。

2.根据权利要求1所述的地震剖面图像的纹理分割方法，其中，所述特征提取基于定义在2n*2n大小窗口中的纹理能量宏特征E(i,j)进行：

E(i,j)=12n*2nΣk=i-n+1i+n-1Σl=j-n+1j+n-1|ω(k,l)|]]>

其中，i,j为像素坐标，n,k,l为自然数，E(i,j)为像素(i,j)的特征值，ω(k,l)表示以像素(i,j)为中心的2n*2n大小窗口的第(k,l)个小波系数。

3.根据权利要求2所述的地震剖面图像的纹理分割方法，其中，所述图像纹理边缘检测包括以下步骤：

对所述灰度图像进行滤波以消除噪声；

基于图像灰度沿X、Y两个方向的偏导数得到梯度的幅值和方向；对梯度的幅值进行非极大抑制，以定位准确的边缘；

遍历经过非极大抑制的图像，基于累计直方图计算包括高阈值和低阈值的双阈值；以及

利用双阈值对经过非极大抑制的图像进行边缘检测。

4.根据权利要求3所述的地震剖面图像的纹理分割方法，其中，对所述灰度图像进行中值滤波，具体公式为：

g(x,y)＝med{f(x-k,y-l),(k,l∈W)}

其中，g(x,y)为滤波输出的图像，f(x,y)为原始图像，med{}为中值滤波函数，W为二维模版。

5.根据权利要求3所述的地震剖面图像的纹理分割方法，其中，利用双阈值对经过非极大抑制的图像进行边缘检测包括：在经过非极大抑制的图像中，灰度值大于高阈值的像素点被确定为边缘，灰度值小于低阈值的被确定为不是边缘，对于其他的像素，如果相邻P个像素的平均值大于高阈值则被确定为边缘，否则被确定为不是边缘。