Question

关于图像处理，利用sobel算子边缘检测的Matlab程序

ps=imread('14.jpg'); %读取图像
subplot(1,3,1)
imshow(ps);
title('原图像');
ps=rgb2gray(ps);
[m,n]=size(ps); %用Sobel微分算子进行边缘检测
pa = edge(ps,'sobel');
subplot(1,3,2);
imshow(pa);
title('边缘检测得到的图像');
ps=double(ps);
for i=1:256
for j=1:256
if pa(i,j)==1
res(i,j)= g(i,j);
else
res(i,j)=ps(i,j);
end
end
end
res=uint8(res);
subplot(1,3,3);
imshow(res);
title('增强后的图像')
imwrite(res,'222.jpg');

才学的这个，不怎么会，这是别人的程序，res(i,j)= g(i,j); 这里出错了，没有定义g，边缘检测图像和原图像都可以显示，增强后的不能显示。我想问下该如何修改？
我基本一点都不会，为了做毕业设计没办法。我意思就是如何简单的改动一下这个程序？不添加额外的东西，就用sobel算子

Answer 1

这很明显啊，你没有定义g函数。看你程序的意思是先检测出边缘，因为边缘中的线都是白色的，然后就是用for 和 if语句把边缘经过处理，也就是将边缘中白色点变为g,这就要看你具体想以何种方式增强了，如果是对数变换的话，你可以先让g(:,:)=log((ps(:,:))+1)（将这条语句放在for语句的上面）;这样增强的结果就是扩展低值灰度，压缩高值灰度，当然了，还有很多其它增强方法，你可以具体再找找图像增强这方面的资料。