contour_line.m

关于偏微分方程进行图像处理的电子书,对于学习动态轮廓实现图像的分割非常有用

function contour_line(image_name,N)
%求一图像的等高线
%输入参数：image_name－－－图像文件名
%        ：N－－－－－－－－阀值（0～255）
if nargin<1 
     disp('请输入文件名!!!');
     return;
elseif nargin<2
     N=128;
end

image_input=imread(image_name);
subplot(1,2,1);
imshow(image_input);
title('原图像');
[size_m,size_n]=size(image_input);

% 求图像的等高线
contour=zeros(size_m,size_n);

for row=1:size_m
    for col=1:size_n
        if image_input(row,col)==N
            contour(row,col)=1;
        end
        contour(row,col)=1-contour(row,col);
    end
end

subplot(1,2,2);
imshow(contour,[]);
title(strcat('图像的等高线(',num2str(N),')'));