使用p参数对灰度图进行图像分割(Matlab)

文章首发及后续更新：https://mwhls.top/2136.html
新的更新内容请到mwhls.top查看。
无图/无目录/格式错误/更多相关请到上方的文章首发页面查看。

目录

1. 目的

2. 原理

3. 实验结果

4. 代码

目的

• 使用p参数算法分割灰度图像。

原理

1. 求灰度直方图
2. 根据灰度直方图，将灰度直方图分割，以数量处于中间的灰度值作为阈值（但因为图像不是双峰灰度直方图，因此是取了前20%）。
3. 设置大小与灰度图大小一致的二值数组，其中，灰度图像素的灰度值大于阈值，则数组对应位置为1，否则为0.
4. 将数组与灰度图相乘，可以得到一个像素灰度值只在0或(阈值, 255)间的图像，即分割后的图像。

实验结果

代码

clear all;
ima = imread('2.jpg');
grayIma = rgb2gray(ima);
[height, width] = size(grayIma);
histArray = zeros(1, 256);
array = zeros(height, width);
 
for row = 1:1:height
    for col = 1:1:width
        temp = grayIma(row, col) + 1;
        histArray(1, temp) = histArray(1, temp) + 1;
    end
end
mid = sum(histArray)/5;
sum = 0;
for pos = 1:1:256
    sum = sum + histArray(1,pos);
    if(sum > mid)
        divisionValue = pos;
        break;
    end
end
 
for row = 1:1:height
    for col = 1:1:width
        if(grayIma(row, col) < divisionValue)
            array(row, col) = 0;
        else
            array(row, col) = 1;
        end
    end
end
 
subplot(2,2,1);
imshow(grayIma);
subplot(2,2,3);
bar(histArray);
 
grayIma2Division = grayIma;
for row = 1:1:height
    for col = 1:1:width
        grayIma2Division(row, col) = grayIma(row, col) * array(row, col);
    end
end
 
 
subplot(2,2,2);
imshow(grayIma2Division);
subplot(2,2,4);
imhist(grayIma2Division);