参考：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)