class Model(object):
def __init__(self, param):
self.c1 = ConvLayer((1, 28, 28), (16, 5, 5), (1, 0), Relu(), param)
# 4x24x24
self.p1 = PoolingLayer(self.c1.output_shape, (
2,
2,
), 2, PoolingTypes.MAX)
# 4x12x12
self.c2 = ConvLayer(self.p1.output_shape, (8, 3, 3), (1, 0), Relu(),
param)
# 4x10x10
self.p2 = PoolingLayer(self.c2.output_shape, (
2,
2,
), 2, PoolingTypes.MAX)
# 4x5x5
self.f1 = FcLayer(self.p2.output_size, 32, Relu(), param)
self.f2 = FcLayer(self.f1.output_size, 10, Softmax(), param)
def forward(self, x):
net = self.c1.forward(x)
net = self.p1.forward(net)
net = self.c2.forward(net)
net = self.p2.forward(net)
net = self.f1.forward(net)
net = self.f2.forward(net)
self.output = net
return self.output
def backward(self, y):
delta = self.output - y
delta = self.f2.backward(delta, LayerIndexFlags.LastLayer)
delta = self.f1.backward(delta, LayerIndexFlags.MiddleLayer)
delta = self.p2.backward(delta, LayerIndexFlags.MiddleLayer)
delta = self.c2.backward(delta, LayerIndexFlags.MiddleLayer)
delta = self.p1.backward(delta, LayerIndexFlags.MiddleLayer)
delta = self.c1.backward(delta, LayerIndexFlags.FirstLayer)
def update(self, learning_rate):
self.c1.update()
self.c2.update()
self.f1.update()
self.f2.update()
def save(self):
self.c1.save_parameters("c1")
self.c2.save_parameters("c2")
self.p1.save_parameters("p1")
self.p2.save_parameters("p2")
self.f1.save_parameters("f1")
self.f2.save_parameters("f2")
def load(self):
self.c1.load_parameters("c1")
self.c2.load_parameters("c2")
self.p1.load_parameters("p1")
self.p2.load_parameters("p2")
self.f1.load_parameters("f1")
self.f2.load_parameters("f2")
def try_filters(file_name):
img = cv2.imread(file_name)
# cv2 format is:G B R, change it to R G B
img1 = img[:, :, [2, 1, 0]]
#plt.imshow(img2)
#plt.show()
img2 = cv2.cvtColor(img1, cv2.COLOR_RGB2GRAY)
batch_size = 1
input_channel = 1
(height, width) = img2.shape
FH = 3
FW = 3
print(img2.shape)
data = img2.reshape((1, 1, height, width))
hp = HyperParameters_4_2(0.1,
10,
batch_size,
net_type=NetType.MultipleClassifier,
init_method=InitialMethod.Xavier,
optimizer_name=OptimizerName.Momentum)
conv = ConvLayer((1, height, width), (1, FH, FW), (1, 1), hp)
conv.initialize("know_cnn", "name")
filters = [
np.array([0, -1, 0, -1, 5, -1, 0, -1, 0]), # sharpness filter
np.array([0, 0, 0, -1, 2, -1, 0, 0, 0]), # vertical edge
np.array([1, 1, 1, 1, -9, 1, 1, 1, 1]), # surround
np.array([-1, -2, -1, 0, 0, 0, 1, 2, 1]), # sobel y
np.array([0, 0, 0, 0, 1, 0, 0, 0, 0]), # nothing
np.array([0, -1, 0, 0, 2, 0, 0, -1, 0]), # horizontal edge
np.array([0.11, 0.11, 0.11, 0.11, 0.11, 0.11, 0.11, 0.11,
0.11]), # blur
np.array([-1, 0, 1, -2, 0, 2, -1, 0, 1]), # sobel x
np.array([2, 0, 0, 0, -1, 0, 0, 0, -1])
] # embossing
filters_name = [
"sharpness", "vertical edge", "surround", "sobel y", "nothing",
"horizontal edge", "blur", "sobel x", "embossing"
]
fig, ax = plt.subplots(nrows=3, ncols=3, figsize=(9, 9))
for i in range(len(filters)):
filter = np.repeat(filters[i],
input_channel).reshape(batch_size, input_channel,
FH, FW)
conv.set_filter(filter, None)
z = conv.forward(data)
#z = normalize(z, 255)
ax[i // 3, i % 3].imshow(z[0, 0])
ax[i // 3, i % 3].set_title(filters_name[i])
ax[i // 3, i % 3].axis("off")
plt.suptitle("filters")
plt.show()
return z
class Model(object):
def __init__(self, param):
self.c1 = ConvLayer((1, 28, 28), (4, 5, 5), (1, 0), Relu(), param)
# 4x24x24
self.p1 = PoolingLayer(self.c1.output_shape, (
2,
2,
), 2, PoolingTypes.MAX)
# 4x12x12
self.f1 = FcLayer(self.p1.output_size, 32, Sigmoid(), param)
self.f2 = FcLayer(self.f1.output_size, 10, Softmax(), param)
def forward(self, x):
a_c1 = self.c1.forward(x)
a_p1 = self.p1.forward(a_c1)
a_f1 = self.f1.forward(a_p1)
a_f2 = self.f2.forward(a_f1)
self.output = a_f2
return self.output
def backward(self, y):
delta_in = self.output - y
d_f2 = self.f2.backward(delta_in, LayerIndexFlags.LastLayer)
d_f1 = self.f1.backward(d_f2, LayerIndexFlags.MiddleLayer)
d_p1 = self.p1.backward(d_f1, LayerIndexFlags.MiddleLayer)
d_c1 = self.c1.backward(d_p1, LayerIndexFlags.FirstLayer)
def update(self, learning_rate):
self.c1.update()
self.f1.update()
self.f2.update()
def save(self):
self.c1.save_parameters("c1")
self.p1.save_parameters("p1")
self.f1.save_parameters("f1")
self.f2.save_parameters("f2")
def load(self):
self.c1.load_parameters("c1")
self.p1.load_parameters("p1")
self.f1.load_parameters("f1")
self.f2.load_parameters("f2")
def conv_relu_pool():
img = cv2.imread(circle_pic)
#img2 = cv2.cvtColor(img, cv2.COLOR_RGB2GRAY)
batch_size = 1
(height, width, input_channel) = img.shape
FH = 3
FW = 3
data = np.transpose(img, axes=(2, 1, 0)).reshape(
(batch_size, input_channel, width, height))
hp = HyperParameters_4_2(0.1,
10,
batch_size,
net_type=NetType.MultipleClassifier,
init_method=InitialMethod.Xavier,
optimizer_name=OptimizerName.Momentum)
conv = ConvLayer((input_channel, width, height), (1, FH, FW), (1, 0), hp)
conv.initialize("know_cnn", "conv")
kernal = np.array([-1, 0, 1, -2, 0, 2, -1, 0, 1])
filter = np.repeat(kernal, input_channel).reshape(batch_size,
input_channel, FH, FW)
conv.set_filter(filter, None)
z1 = conv.forward(data)
z2 = Relu().forward(z1)
pool = PoolingLayer(z2[0].shape, (2, 2), 2, PoolingTypes.MAX)
pool.initialize("know_cnn", "pool")
z3 = pool.forward(z2)
fig, ax = plt.subplots(nrows=2, ncols=2, figsize=(8, 6))
ax[0, 0].imshow(img[:, :, [2, 1, 0]])
ax[0, 0].axis("off")
ax[0, 0].set_title("source:" + str(img.shape))
ax[0, 1].imshow(z1[0, 0].T)
ax[0, 1].axis("off")
ax[0, 1].set_title("conv:" + str(z1.shape))
ax[1, 0].imshow(z2[0, 0].T)
ax[1, 0].axis("off")
ax[1, 0].set_title("relu:" + str(z2.shape))
ax[1, 1].imshow(z3[0, 0].T)
ax[1, 1].axis("off")
ax[1, 1].set_title("pooling:" + str(z3.shape))
plt.suptitle("conv-relu-pool")
plt.show()