[I,map]=imread('img1.bmp'); %读入图像
imshow(I,map); %显示原图像
I=double(I); %转换为double类型
[Gx,Gy]=gradient(I); % 计算梯度
G=sqrt(Gx.*Gx+Gy.*Gy); % 水平垂直差分
J=I;
K=find(G>=7); %指定灰度级
J(K)=255;
figure,imshow(J,map); % 显示处理后的图像
转自:http://zhidao.baidu.com/question/160216024.html
01 | function ruihua_Callback(hObject, eventdata, handles) |
02 | % hObject handle to Untitled_1 (see GCBO) |
03 | % eventdata reserved - to be defined in a future version of MATLAB |
04 | % handles structure with handles and user data (see GUIDATA) |
07 | [fname,fpath]=uigetfile('*.bmp;*.jpg','open a file'); |
08 | filename=[fpath,fname]; |
09 | [I,map]=imread(filename); |
12 | [Gx,Gy]=gradient(I); % 计算梯度 |
13 | G=sqrt(Gx.*Gx+Gy.*Gy); % 注意是矩阵点乘 |
38 | subplot(2,3,1),imshow(I,map); |
40 | subplot(2,3,2),imshow(J1,map); |
42 | subplot(2,3,3),imshow(J2,map); |
44 | subplot(2,3,4),imshow(J3,map); |
46 | subplot(2,3,5),imshow(J4,map); |
48 | subplot(2,3,6),imshow(J5,map); |
转自:http://www.oschina.net/code/snippet_104512_4395