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

目的

  • 使用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);

版权声明:
作者:MWHLS
链接:http://panwj.top/2136.html
来源:无镣之涯
文章版权归作者所有,未经允许请勿转载。

THE END
分享
二维码
打赏
< <上一篇
下一篇>>