Matlab中的多尺度形态图像简化

时间:2013-11-15 06:37:45

标签: matlab image-processing morphological-analysis mathematical-morphology multiscaleimage

大家好我正在尝试使用Dorini论文“多尺度形态图像简化”的SMMT操作符。由于在没有订阅的情况下无法访问页面,我在此处发布相关详细信息:

enter image description here enter image description here enter image description here

请注意,我将相关文章的部分内容发布为图片。我不知道如何在stackoverflow.com中编写方程式。我想使用这个SMMT操作符作为图像处理中的预处理步骤。我写下的代码如下:

clc;clear all;close all;
tic
I=imread('handwritten.jpg');
I=I(:,:,1);

dim=11 ;
HEIGHT=zeros(dim,dim);
sigma=1/10;
for i=-floor(dim/2):floor(dim/2)
    for j=-floor(dim/2):floor(dim/2)
       HEIGHT(i+ceil(dim/2),j+ceil(dim/2))=-(1/sigma).*max(abs(i),abs(j));
    end
end
NHOOD=ones(dim,dim);
se = strel('arbitrary',NHOOD,HEIGHT);
se

IM1 = imdilate(I,se,'same');
IM2 = imerode(I,se,'same');
figure;
subplot(2,2,1),imshow(I)
subplot(2,2,2),imshow(IM1)
subplot(2,2,3),imshow(IM2)

II = I;
for i=1:1
    phi1 = imdilate(II,se,'same');
    phi2 = imerode(II,se,'same');
    for j=1:size(I,1)
        for k=1:size(I,2)
            if ((phi1(j,k)-II(j,k))<(II(j,k)-phi2(j,k)))
                II(j,k) = phi1(j,k);
            elseif ((phi1(j,k)-II(j,k))==(II(j,k)-phi2(j,k)))                
                II(j,k) = II(j,k);
            else
                II(j,k) = phi1(j,k);
            end
        end
    end
end
IM3=II;
subplot(2,2,4),imagesc(IM3,[0 255]);colormap('gray');axis off;
toc

代码的结果应该是这样的:(再次来自论文): enter image description here

我的结果如下: - enter image description here

我的实施是否正确?可以进一步改进??任何建议都会有所帮助。在此先感谢您的帮助! 有关详细信息,请参阅Dorini Free Access

1 个答案:

答案 0 :(得分:0)

这些实现在全球范围内看起来不错,但有一点澄清:

if ((phi1(j,k)-II(j,k))<(II(j,k)-phi2(j,k)))
                II(j,k) = phi1(j,k);

我想在这里你想要在比较IIphi1时更改变量I中的图片值:

if ((phi1(j,k)-I(j,k))<(I(j,k)-phi2(j,k)))
                II(j,k) = phi1(j,k);