参考:https://ww2.mathworks.cn/help/images/correcting-nonuniform-illumination.html
总的步骤为:dilation——从原始图像 I
中减去背景逼近图像 background——提升对比度——二值化操作——去除噪点
下面介绍几个代码中用到的函数:(下面的代码适用于灰度值图像或者二值化图像,如果是彩色图像,可以更改SE为非扁平的,即三维的)
1.strel:该函数返回一个flat SE(structure element):SE和卷积核的区别是其有中心点,它是进行dilation(扩张)和erosion(腐蚀)的基本组成单元。
在matlab中,flat SE和 NonFlat SE(由offsetstrel创建
)包括下面几类:
参考:https://ww2.mathworks.cn/help/images/ref/strel.html#bu7pnv2-1
2.imopen函数:对灰度或者二值图像执行形态学开运算(dilation),返回形态学开运算后的图像,该函数需要输入两个参数,一个是图像,一个是SE(structure element)。
参考:https://ww2.mathworks.cn/help/images/ref/imopen.html
dilation的过程可以参考我之前的博文https://blog.csdn.net/weixin_44719615/article/details/105987122
3.imadjust函数:
将灰度图像 J
= imadjust(I
)I
中的强度值映射到 J
中的新值。默认情况下,imadjust
对所有像素值中最低的 1% 和最高的 1% 进行饱和处理。此运算可提高输出图像 J
的对比度。
参考:https://ww2.mathworks.cn/help/images/ref/imadjust.html
4. imbinarize函数:对图像进行二值化操作
参考:https://ww2.mathworks.cn/help/images/ref/imbinarize.html
5.bwareaopen函数:BW2
= bwareaopen(BW
,P
)从二值图像 BW
中删除少于 P
个像素的所有连通分量(对象),并生成另一个二值图像 BW2
。此运算称为面积开运算。
参考:https://ww2.mathworks.cn/help/images/ref/bwareaopen.html
%matlab代码
I = imread(fullfile('D://','fuji.jpg'));
I = rgb2gray(I);
imshow(I)
d = imdistline(gca);
delete(d)
SE = strel('square', 50 )
background = imopen(I,SE);
imshow(background)
I2 = I - background;
imshow(I2)
I3 = imadjust(I2);
imshow(I3)
bw = imbinarize(I3);
bw = bwareaopen(bw,50);
imshow(bw)