Необходимо наложить контуры на исходное изображение, после выделения границ различными методами

im=imread('mdb028.bmp');%загрузка изображения
imshow(im);%открытие изображения
isind(im);%проверка является ли палитровым
ans%если 1 - является палитровым
subplot(2,2,1);imshow(im);%вывыод исходного изображения и гистограммы в одно окно
title('палитровое изображение mdb028.bmp');
subplot(2,2,2);imhist(im);
title('гистограмма палитрового изображение');
[im,map]=imread('mdb028.bmp');
i=ind2gray(im,map);%преобразование палитрового в полутоновое изображение
subplot(2,2,3);imshow(i);%вывыод палитрового изображения и гистограммы в одно окно
title('полутоновое изображение mdb028.bmp');
subplot(2,2,4);imhist(i);
title('гистограмма полутонового изображение');
isgray(i);%1 - является полутоновым
ans
%эквализация гистограммы
i1=histeq(i);
figure, imshow(i1);
title('эквализированное изображение');
figure, imhist(i1);
title('гистограмма эквализированного');
%Выделение границ методом Канни
BW = edge(i1,'canny');
figure, imshow(BW);
title('Выделение границ методом Канни');
%Выделение границ методом Собеля
BW = edge(i1,'sobel');
figure, imshow(BW);
title('Выделение границ методом Собеля');
%Выделение границ методом Превита
BW = edge(i1,'prewitt');
figure, imshow(BW);
title('Выделение границ методом Превита');
%Выделение границ методом Робертса
BW=edge(i1, 'roberts'); 
figure, imshow(BW);
title('Выделение границ методом Робертса');
%Выделение границ методом лапласиан-гауссиана 
BW = edge (i1,'log');
figure, imshow(BW);
title('Выделение границ методом лапласиан-гауссиана ');
%Выделение границ методом zerocross
BW=edge(i1, 'zerocross'); 
figure, imshow(BW);
title('Выделение границ методом zerocross');

im=imread('image.jpg');
g=edge(rgb2gray(im),'sobel');
im1=im;
[n,m,k]=size(im);
for i=1:n
for j=1:m
if (g(i,j)==1)
im1(i,j,1)=255;
im1(i,j,2)=0;
im1(i,j,3)=0;
end
end
end 
imshow(im1);