-
Notifications
You must be signed in to change notification settings - Fork 1
/
sharppolished.m
67 lines (60 loc) · 2.03 KB
/
sharppolished.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
function [sharp,pic]=sharppolished(pic)
% clear;clc;
% pre processing
% r = medfilt2(pic(:, :, 1), [5 5]);
% g = medfilt2(pic(:, :, 2), [5 5]);
% b = medfilt2(pic(:, :, 3), [5 5]);
% pic = cat(3, r, g, b);
pic = imsharpen(pic,'Radius',7,'Amount',0.5);
pic2=pic;
% my filter
% figure;
% pic=rgb2gray(pic);
pic3=pic;
pic = wiener2(pic,[5 5]);
pic4=pic;
polished1=[-1.2 -0.8 -0.6 ;0 0 0 ;1.2 0.8 0.6 ];
polished11=[1.2 0.8 0.6 ;0 0 0 ;-1.2 -0.8 -0.6 ];
polished111=[-0.1 -0.8 -0.6 ;0 0 0;0.1 0.8 0.6 ];
polished1111=[0.1 0.8 0.6 ;0 0 0 ;-0.1 -0.8 -0.6 ];
polished2=polished1';
polished22=polished11';
polished222=polished111';
polished2222=polished1111';
polished1filter=imfilter(pic,polished1);
polished2filter=imfilter(pic,polished2);
polished11filter=imfilter(pic,polished11);
polished22filter=imfilter(pic,polished22);
polished111filter=imfilter(pic,polished111);
polished222filter=imfilter(pic,polished222);
polished1111filter=imfilter(pic,polished1111);
polished2222filter=imfilter(pic,polished2222);
polished=abs(polished1filter)+abs(polished2filter)+abs(polished11filter)+abs(polished22filter)+abs(polished111filter)+abs(polished222filter)+abs(polished1111filter)+abs(polished2222filter);
% subplot(2,2,1),imshow(pic);
% title('Original');
% subplot(2,2,2),imshow(polished1filter);
% title('Horizontal Edges');
% subplot(2,2,3),imshow(polished2filter);
% title('Vertical Edges');
% subplot(2,2,4),imshow(polished);
% title('Horizontal And Vertical Abs Edges');
% post processing
[d1,d2] = size(pic)
for i=1:d1
for j=1:d2
if polished(i,j)<60
sharp(i,j)=0;
else
sharp(i,j)=polished(i,j);
end
end
end
% figure;
% imshow(sharp);title('Edges');
% figure('units','normalized','outerposition',[0 0 1 1])
% imshow(pic2);title('RGB Smooth with medfilt2');
% figure('units','normalized','outerposition',[0 0 1 1])
% imshow(pic3);title('GrayScale Smooth');
% figure('units','normalized','outerposition',[0 0 1 1])
% imshow(pic4);title('GrayScale Smooth with wiener2');
end