-
Notifications
You must be signed in to change notification settings - Fork 0
/
FIltrosConvoluciones.py
169 lines (159 loc) · 9.46 KB
/
FIltrosConvoluciones.py
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
# -*- coding: utf-8 -*-
"""
Created on Sun Sep 2 14:12:49 2018
@author: FamiliaHogar
"""
from skimage import io
import matplotlib.pyplot as plt
import numpy as np
#plt.rcParams['image.cmap'] = 'gray'
matrizOriginal = (io.imread('blackWhite.png')/255.0)/255.0
i, j = matrizOriginal.shape
matrizSuavizada = np.zeros((i, j))
matrizBordeada = np.zeros((i, j))
matrizPerfilada = np.zeros((i, j))
###### Suavizado ######
for x in range(i):
for y in range(j):
if x == 0 and y == 0:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x,(y+1)] + matrizOriginal[x,y]
+ matrizOriginal[x,y] + matrizOriginal[x,y] + matrizOriginal[(x+1),y]
+ matrizOriginal[x,y])/10)
elif x == i-1 and y == 0:
matrizSuavizada[x,y] = ((matrizOriginal[(x),y] + matrizOriginal[(x),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x-1,(y+1)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y] + matrizOriginal[x,y] + matrizOriginal[(x),y]
+ matrizOriginal[x,y])/10)
elif x == i-1 and y == j-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x),y] + matrizOriginal[(x),(y)]
+ matrizOriginal[x,(y)] + matrizOriginal[x-1,(y)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x),y-1]
+ matrizOriginal[x,y])/10)
elif x == 0 and y == j-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y)]
+ matrizOriginal[x,(y)] + matrizOriginal[x,(y)] + matrizOriginal[x,y]
+ matrizOriginal[x,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1]
+ matrizOriginal[x,y])/10)
elif y == 0 and x != 0 and x != i-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x-1,(y+1)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y] + matrizOriginal[x,y] + matrizOriginal[(x+1),y]
+ matrizOriginal[x,y])/10)
elif x == i-1 and y != 0 and y != j-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x),y] + matrizOriginal[(x),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x-1,(y+1)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x),y-1]
+ matrizOriginal[x,y])/10)
elif y == j-1 and x != 0 and x != i-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y)]
+ matrizOriginal[x,(y)] + matrizOriginal[x-1,(y)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1]
+ matrizOriginal[x,y])/10)
elif x == 0 and y != 0 and y != j-1:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x,(y+1)] + matrizOriginal[x,y]
+ matrizOriginal[x,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1]
+ matrizOriginal[x,y])/10)
else:
matrizSuavizada[x,y] = ((matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
+ matrizOriginal[x,(y+1)] + matrizOriginal[x-1,(y+1)] + matrizOriginal[x-1,y]
+ matrizOriginal[x-1,y-1] + matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1]
+ matrizOriginal[x,y])/10)
###### Bordes ######
for x in range(i):
for y in range(j):
if x == 0 and y == 0:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - 2*matrizOriginal[x,y]
- matrizOriginal[x,y] + matrizOriginal[(x+1),y])/4)
elif x == i-1 and y == 0:
matrizBordeada[x,y] = ((2*matrizOriginal[(x),y] + matrizOriginal[(x),(y+1)]
- matrizOriginal[x-1,(y+1)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y] + matrizOriginal[(x),y])/4)
elif x == i-1 and y == j-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x),y] + matrizOriginal[(x),(y)]
- matrizOriginal[x-1,(y)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] + matrizOriginal[(x),y-1])/4)
elif x == 0 and y == j-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y)]
- matrizOriginal[x,(y)] - 2*matrizOriginal[x,y]
- matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1])/4)
elif y == 0 and x != 0 and x != i-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x-1,(y+1)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y] + matrizOriginal[(x+1),y])/4)
elif x == i-1 and y != 0 and y != j-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x),y] + matrizOriginal[(x),(y+1)]
- matrizOriginal[x-1,(y+1)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] + matrizOriginal[(x),y-1])/4)
elif y == j-1 and x != 0 and x != i-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y)]
- matrizOriginal[x-1,(y)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] + matrizOriginal[(x+1),y-1])/4)
elif x == 0 and y != 0 and y != j-1:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - 2*matrizOriginal[x,y]
- matrizOriginal[x,y-1] + matrizOriginal[(x+1),y-1])/4)
else:
matrizBordeada[x,y] = ((2*matrizOriginal[(x+1),y] + matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x-1,(y+1)] - 2*matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] + matrizOriginal[(x+1),y-1])/4)
###### Perfilado ######
for x in range(i):
for y in range(j):
if x == 0 and y == 0:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x,(y+1)] - matrizOriginal[x,y]
- matrizOriginal[x,y] - matrizOriginal[x,y] - matrizOriginal[(x+1),y]
+ 9*matrizOriginal[x,y])/6)
elif x == i-1 and y == 0:
matrizPerfilada[x,y] = ((-matrizOriginal[(x),y] - matrizOriginal[(x),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x-1,(y+1)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y] - matrizOriginal[x,y] - matrizOriginal[(x),y]
+ 9*matrizOriginal[x,y])/6)
elif x == i-1 and y == j-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x),y] - matrizOriginal[(x),(y)]
- matrizOriginal[x,(y)] - matrizOriginal[x-1,(y)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x),y-1]
+ 9*matrizOriginal[x,y])/6)
elif x == 0 and y == j-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y)]
- matrizOriginal[x,(y)] - matrizOriginal[x,(y)] - matrizOriginal[x,y]
- matrizOriginal[x,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x+1),y-1]
+ 9*matrizOriginal[x,y])/6)
elif y == 0 and x != 0 and x != i-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x-1,(y+1)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y] - matrizOriginal[x,y] - matrizOriginal[(x+1),y]
+ 9*matrizOriginal[x,y])/6)
elif x == i-1 and y != 0 and y != j-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x),y] - matrizOriginal[(x),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x-1,(y+1)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x),y-1]
+ 9*matrizOriginal[x,y])/6)
elif y == j-1 and x != 0 and x != i-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y)]
- matrizOriginal[x,(y)] - matrizOriginal[x-1,(y)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x+1),y-1]
+ 9*matrizOriginal[x,y])/6)
elif x == 0 and y != 0 and y != j-1:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x,(y+1)] - matrizOriginal[x,y]
- matrizOriginal[x,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x+1),y-1]
+ 9*matrizOriginal[x,y])/6)
else:
matrizPerfilada[x,y] = ((-matrizOriginal[(x+1),y] - matrizOriginal[(x+1),(y+1)]
- matrizOriginal[x,(y+1)] - matrizOriginal[x-1,(y+1)] - matrizOriginal[x-1,y]
- matrizOriginal[x-1,y-1] - matrizOriginal[x,y-1] - matrizOriginal[(x+1),y-1]
+ 9*matrizOriginal[x,y])/6)
plt.imshow(matrizOriginal,vmin=0,vmax=1)
plt.figure()
plt.imshow(matrizSuavizada,vmin=0,vmax=1)
plt.figure()
plt.imshow(matrizBordeada,vmin=0,vmax=1)
plt.figure()
plt.imshow(matrizPerfilada,vmin=0,vmax=1)
io.imsave("blackWhiteSuave.png", matrizSuavizada)
io.imsave("blackWhiteBordes.png", matrizBordeada)
io.imsave("blackWhitePerfil.png", matrizPerfilada)