I = imread( 'sar.tif' );
%I = rgb2gray(rgb);
subplot(5,3,1),imshow(I);
k=medfilt2(I,[5,5]);
subplot(5,3,2),imshow(k);
title('5*5中值滤波图像');
%f=imread('tuxiang1.tif');
%subplot(1,2,1),imshow(f);
%title('原图像');
g1=histeq(k,256);
subplot(5,3,3),imshow(g1);
title('直方图匹配');
%g2=histeq(k2,256);
%subplot(2,2,2),imshow(g2);
%title('5*5直方图匹配');
%k=medfilt2(f,[5,5]);
%k2=medfilt2(f,[5,5]);
%j=imnoise(f,'gaussian',0,0.005);
%subplot(1,3,3),imshow(k2);
%title('5*5中值滤波图像');
hy = fspecial( 'sobel' );
hx = hy;
Iy = imfilter(double(g1), hy, 'replicate' );
Ix = imfilter(double(g1), hx, 'replicate' );
gradmag = sqrt(Ix.^2 + Iy.^2);
subplot(5,3,4), imshow(gradmag,[ ]), title( 'gradmag' );
L = watershed(gradmag);
Lrgb = label2rgb(L);
subplot(5,3,5), imshow(Lrgb), title( 'Lrgb' );
se = strel( 'disk' , 9);
Io = imopen(g1, se);
subplot(5,3,6), imshow(Io), title( 'Io' )
Ie = imerode(g1, se);
Iobr = imreconstruct(Ie, g1);
subplot(5,3,7), imshow(Iobr), title( 'Iobr' );
Ioc = imclose(Io, se);
subplot(5,3,8), imshow(Ioc), title( 'Ioc' );
Iobrd = imdilate(Iobr, se);
Iobrcbr = imreconstruct(imcomplement(Iobrd), imcomplement(Iobr));
Iobrcbr = imcomplement(Iobrcbr);
subplot(5,3,9), imshow(Iobrcbr), title( 'Iobrcbr' );
fgm = imregionalmax(Iobrcbr);
subplot(5,3,10), imshow(fgm), title( 'fgm' );
I2 = g1; I2(fgm) = 255;
subplot(5,3,11),imshow(I2), title( 'fgm superimposed on original image' );
se2 = strel(ones(5,5)); I3 = g1; I3(fgm) = 255;
subplot(5,3,12) ,imshow(I3);
title( 'fgm4 superimposed on original image' );
bw = im2bw(Iobrcbr, graythresh(Iobrcbr));
subplot(5,3,13) , imshow(bw), title( 'bw' );
D = bwdist(bw); DL = watershed(D);
bgm = DL == 0;
subplot(5,3,14) , imshow(bgm), title( 'bgm' );
gradmag2 = imimposemin(gradmag, bgm | fgm);
L = watershed(gradmag2);
I4 = g1;
I4(imdilate(L == 0, ones(3, 3)) | bgm | fgm) = 255;
figure, imshow(I4);
title( 'Markers and object boundaries superimposed on original image' );
Lrgb = label2rgb(L, 'jet' , 'w' , 'shuffle' );
figure, imshow(Lrgb);
title( 'Lrgb' );
figure, imshow(I), hold on
himage = imshow(Lrgb);
set(himage, 'AlphaData' , 0.3);
title( 'Lrgb superimposed transparently on original image' );
