% 形态学滤波
clc,clear,close all % 清理命令区、清理工作区、关闭显示图形
warning off % 消除警告
feature jit off % 加速代码运行
tic
[filename ,pathname]=...
uigetfile({'*.bmp';'*.tif';'*.jpg';},'选择图片'); %选择图片路径
str=[pathname filename]; % 合成路径+文件名
im = imread(str); % 读图
im = imnoise(im,'gaussian',0,1e-3); % 原图像 + 白噪声
% 形态学滤波
sca = 5; % 结构元素尺寸
im_e = morphology_filter(im,sca); % 形态学滤波
figure('color',[1,1,1]),
subplot(121),imshow(im);title('原始图像')
colormap(jet) % 颜色
shading interp % 消隐
subplot(122),imshow(im_e,[]);title('形态学滤波图像')
colormap(jet) % 颜色
shading interp % 消隐
toc
function Iobrcbr = morphology_filter(im,sca)
% 形态学滤波器
%函数输入:
% im: 输入的二维图像矩阵
% sca : 结构元素尺寸代销
%函数输出:
% Iobrcbr: 形态学滤波图像
% 形态学滤波
% sca = 5; % 结构元素尺寸
se=strel('diamond',(sca-1)/2); % 形态学结构元素
Io=imopen(im,se); % 开操作
Ioc=imclose(Io,se); % 闭操作
Iobr=imreconstruct(Io,Ioc); % 结构重建
Iobrd = imdilate(Iobr,se); % 膨胀
Iobrcbr = imreconstruct(imcomplement(Iobrd),imcomplement(Iobr)); % 结构重建
Iobrcbr = imcomplement(Iobrcbr);
MATLAB---形态学滤波
于 2022-04-14 00:26:25 首次发布