用matlab分水岭算法做海岸线的提取

最新推荐文章于 2024-05-06 09:30:12 发布

原创

最新推荐文章于 2024-05-06 09:30:12 发布 · 1w 阅读

96 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #算法

本文介绍了如何使用MATLAB进行海岸线提取。通过读取图像，将其转化为灰度图，然后利用sobel算子计算梯度，接着进行形态学操作如开、闭操作，并结合分水岭算法进行图像分割。最终通过一系列处理，得到清晰的海岸线图像。

先上原图：

下面是程序：

rgb = imread('你图片文件');%读取原图像
I = rgb2gray(rgb);%转化为灰度图像
figure; subplot(121)%显示灰度图像
imshow(I)
text(732,501,'Image courtesy of Corel',...
'FontSize',7,'HorizontalAlignment','right')
hy = fspecial('sobel');%sobel算子
hx = hy';
Iy = imfilter(double(I), hy, 'replicate');%滤波求y方向边缘
Ix = imfilter(double(I), hx, 'replicate');%滤波求x方向边缘
gradmag = sqrt(Ix.^2 + Iy.^2);%求摸
subplot(122); imshow(gradmag,[]), %显示梯度
title('Gradient magnitude (gradmag)')
L = watershed(gradmag);%直接应用分水岭算法
Lrgb = label2rgb(L);%转化为彩色图像
figure; imshow(Lrgb), %显示分割后的图像
title('Watershed transform of gradient magnitude (Lrgb)')
se = strel('disk', 20);%圆形结构元素
Io = imopen(I, se);%形态学开操作
figure