CIFAR10数据集来源：torchvision.datasets.CIFAR10()
CIFAR10数据集是一个用于识别普适物体的小型数据集，一共包含10个类别的RGB彩色图片，图片尺寸大小为32x32，如图：
CIFAR10.png
相较于MNIST数据集，MNIST数据集是28x28的单通道灰度图，而CIFAR10数据集是32x32的RGB三通道彩色图，CIFAR10数据集更接近于真实世界的图片。
ResNet网络模型：
ResNet.png
本文采用ResNet18来构建深度网络模型:
ResNet18.png
SeNet：Squeeze-and-Excitation的缩写，特征压缩与激发，基于通道注意力
1.Squeeze-and-Excitation（SE） block 并不是一个完整的网络结构，而是一个子结构，可以嵌到其他分类或检测模型中
2.SENet网络的创新点在于关注 channel 之间的关系，希望模型可以自动学习到不同 channel 特征的重要程度
3.本质上，SeNet基于注意力机制给与每个通道不同的权重，SE模块是在 channel 维度上做 attention
具体操作过程如下：
对于每一输出通道，先 global average pool，每个通道得到 1个标量，C个通道得到C个数，
image.png
然后经过 FC(in,in/16)-ReLU-FC(in,in/16)-Sigmoid 得到 C个0~1 之间的标量，作为通道的加权，
然后原来的输出通道每个通道用对应的权重进行加权（对应通道的每个元素与权重分别相乘），得到新的加权后的特征，作者称为 feature recalibration
SeResNet.png
1. 数据集构建
每个像素点即每条数据中的值范围为0-255，有的数字过大不利于训练且难以收敛，故将其归一化到（0-1）之间
# 数据集处理
# transforms.RandomHorizontalFlip(p=0.5)---以0.5的概率对图片做水平横向翻转
# transforms.RandomCrop(32, padding=4)---填充到40*40后，再随机裁剪成32*32
transform_train = transforms.Compose([transforms.RandomHorizontalFlip(p=0.5),
transforms.RandomCrop(32, padding=4),
transforms.ToTensor(),
transforms.Normalize((0.4914, 0.4822, 0.4465), (0.2023, 0.1994, 0.2010))])
# transforms.ToTensor()---shape从(H,W,C)->(C,H,W), 每个像素点从(0-255)映射到(0-1):直接除以255
# transforms.Normalize---先将输入归一化到(0,1),像素点通过"(x-mean)/std",将每个元素分布到(-1,1)
transform = transforms.Compose([transforms.ToTensor(),
transforms.Normalize((0.4914, 0.4822, 0.4465), (0.2023, 0.1994, 0.2010))])
train_dataset = datasets.CIFAR10(root="../DataSet/cifar10", train=True, transform=transform_train,
download=True)
test_dataset = datasets.CIFAR10(root="../DataSet/cifar10", train=False, transform=transform,
download=True)
train_loader = DataLoader(train_dataset, batch_size=128, shuffle=True, num_workers=2)
test_loader = DataLoader(test_dataset, batch_size=128, shuffle=False, num_workers=2)
2.构建SE-ResNet 网络模型，最后接Softmax来处理output
1）构建 SE-Block单元
# SE-Block单元--SEblock是一个子结构，几乎可以嵌入任何一个神经网络模型之中
class SE_Block(nn.Module):
def __init__(self, input_channel, reduction=16):
super(SE_Block, self).__init__()
self.adaptive_avg_pool = nn.AdaptiveAvgPool2d(1) # 全局自适应池化
self.fc = nn.Sequential(
nn.Linear(input_channel, input_channel // reduction),
nn.ReLU(inplace=True),
nn.Linear(input_channel // reduction, input_channel),
nn.Sigmoid()
)
def forward(self, x):
b, c, h, w = x.size()
# squeeze操作:(b,c,h,w)->(b,c)
y = self.adaptive_avg_pool(x).view(b, c)
# FC获取通道注意力权重，是具有全局信息的
y = self.fc(y).view(b, c, 1, 1)
# 注意力作用每一个通道上
y = x * y.expand_as(x)
# 残差连接
return x+y
2）ResNet18_BasicBlock-残差单元
# 构建 VGGNet18 网络模型
# 1.ResNet18_BasicBlock-残差单元
class ResNet18_BasicBlock(nn.Module):
def __init__(self, input_channel, output_channel, stride, use_conv1_1):
super(ResNet18_BasicBlock, self).__init__()
# 第一层卷积
self.conv1 = nn.Conv2d(input_channel, output_channel, kernel_size=3, stride=stride, padding=1)
# 第二层卷积
self.conv2 = nn.Conv2d(output_channel, output_channel, kernel_size=3, stride=1, padding=1)
# 1*1卷积核，在不改变图片尺寸的情况下给通道升维
self.extra = nn.Sequential(
nn.Conv2d(input_channel, output_channel, kernel_size=1, stride=stride, padding=0),
nn.BatchNorm2d(output_channel)
)
self.use_conv1_1 = use_conv1_1
self.bn = nn.BatchNorm2d(output_channel)
self.relu = nn.ReLU(inplace=True)
def forward(self, x):
out = self.bn(self.conv1(x))
out = self.relu(out)
out = self.bn(self.conv2(out))
# 残差连接-(B,C,H,W)维度一致才能进行残差连接
if self.use_conv1_1:
out = self.extra(x) + out
out = self.relu(out)
return out
3. 构建损失函数和优化器
损失函数采用CrossEntropyLoss
优化器采用 SGD 随机梯度优化算法
# 构造损失函数和优化器
criterion = nn.CrossEntropyLoss()
opt = optim.SGD(model.parameters(), lr=0.01, momentum=0.8, weight_decay=5e-4)
# 动态更新学习率------每隔step_size : lr = lr * gamma
schedule = optim.lr_scheduler.StepLR(opt, step_size=10, gamma=0.6, last_epoch=-1)
4.完整代码
# -*- codeing = utf-8 -*-
# @Time : 2022/6/14 13:09
# @Software : PyCharm
import torch
from torch import nn, optim
from torch.utils.data import DataLoader
from torchvision import datasets, transforms
from matplotlib import pyplot as plt
import time
device = torch.device("cuda:0" if torch.cuda.is_available() else "cpu")
# transforms.RandomHorizontalFlip(p=0.5)---以0.5的概率对图片做水平横向翻转
# transforms.RandomCrop(32, padding=4)---填充到40*40后，再随机裁剪成32*32
transform_train = transforms.Compose([transforms.RandomHorizontalFlip(p=0.5),
transforms.RandomCrop(32, padding=4),
transforms.ToTensor(),
transforms.Normalize((0.4914, 0.4822, 0.4465), (0.2023, 0.1994, 0.2010))])
# transforms.ToTensor()---shape从(H,W,C)->(C,H,W), 每个像素点从(0-255)映射到(0-1):直接除以255
# transforms.Normalize---先将输入归一化到(0,1),像素点通过"(x-mean)/std",将每个元素分布到(-1,1)
transform = transforms.Compose([transforms.ToTensor(),
transforms.Normalize((0.4914, 0.4822, 0.4465), (0.2023, 0.1994, 0.2010))])
train_dataset = datasets.CIFAR10(root="../DataSet/cifar10", train=True, transform=transform_train,
download=True)
test_dataset = datasets.CIFAR10(root="../DataSet/cifar10", train=False, transform=transform,
download=True)
train_loader = DataLoader(train_dataset, batch_size=128, shuffle=True, num_workers=2)
test_loader = DataLoader(test_dataset, batch_size=128, shuffle=False, num_workers=2)
# SE-Block单元--SEblock是一个子结构，几乎可以嵌入任何一个神经网络模型之中
class SE_Block(nn.Module):
def __init__(self, input_channel, reduction=16):
super(SE_Block, self).__init__()
self.adaptive_avg_pool = nn.AdaptiveAvgPool2d(1) # 全局自适应池化
self.fc = nn.Sequential(
nn.Linear(input_channel, input_channel // reduction),
nn.ReLU(inplace=True),
nn.Linear(input_channel // reduction, input_channel),
nn.Sigmoid()
)
def forward(self, x):
b, c, h, w = x.size()
# squeeze操作:(b,c,h,w)->(b,c)
y = self.adaptive_avg_pool(x).view(b, c)
# FC获取通道注意力权重，是具有全局信息的
y = self.fc(y).view(b, c, 1, 1)
# 注意力作用每一个通道上
y = x * y.expand_as(x)
# 残差连接
return x+y
# ResNet18_BasicBlock-残差单元
class ResNet18_BasicBlock(nn.Module):
def __init__(self, input_channel, output_channel, stride, use_conv1_1):
super(ResNet18_BasicBlock, self).__init__()
# 第一层卷积
self.conv1 = nn.Conv2d(input_channel, output_channel, kernel_size=3, stride=stride, padding=1)
# 第二层卷积
self.conv2 = nn.Conv2d(output_channel, output_channel, kernel_size=3, stride=1, padding=1)
# 1*1卷积核，在不改变图片尺寸的情况下给通道升维
self.extra = nn.Sequential(
nn.Conv2d(input_channel, output_channel, kernel_size=1, stride=stride, padding=0),
nn.BatchNorm2d(output_channel)
)
self.use_conv1_1 = use_conv1_1
self.bn = nn.BatchNorm2d(output_channel)
self.relu = nn.ReLU(inplace=True)
def forward(self, x):
out = self.bn(self.conv1(x))
out = self.relu(out)
out = self.bn(self.conv2(out))
# 残差连接-(B,C,H,W)维度一致才能进行残差连接
if self.use_conv1_1:
out = self.extra(x) + out
out = self.relu(out)
return out
# 构建 ResNet18 网络模型
class ResNet18(nn.Module):
def __init__(self):
super(ResNet18, self).__init__()
self.conv1 = nn.Sequential(
nn.Conv2d(3, 64, kernel_size=3, stride=1, padding=1, bias=False),
nn.BatchNorm2d(64),
# nn.MaxPool2d(kernel_size=3, stride=2, padding=1),
nn.ReLU(inplace=True)
)
self.block1_1 = ResNet18_BasicBlock(input_channel=64, output_channel=64, stride=1, use_conv1_1=False)
self.se_block64 = SE_Block(64)
self.block1_2 = ResNet18_BasicBlock(input_channel=64, output_channel=64, stride=1, use_conv1_1=False)
self.block2_1 = ResNet18_BasicBlock(input_channel=64, output_channel=128, stride=2, use_conv1_1=True)
self.se_block128 = SE_Block(128)
self.block2_2 = ResNet18_BasicBlock(input_channel=128, output_channel=128, stride=1, use_conv1_1=False)
self.block3_1 = ResNet18_BasicBlock(input_channel=128, output_channel=256, stride=2, use_conv1_1=True)
self.se_block256 = SE_Block(256)
self.block3_2 = ResNet18_BasicBlock(input_channel=256, output_channel=256, stride=1, use_conv1_1=False)
self.block4_1 = ResNet18_BasicBlock(input_channel=256, output_channel=512, stride=2, use_conv1_1=True)
self.se_block512 = SE_Block(512)
self.block4_2 = ResNet18_BasicBlock(input_channel=512, output_channel=512, stride=1, use_conv1_1=False)
# 全连接层
self.FC_layer = nn.Sequential(
nn.Linear(512, 256),
nn.ReLU(inplace=True),
# 使一半的神经元不起作用，防止参数量过大导致过拟合
nn.Dropout(0.5),
nn.Linear(256, 128),
nn.ReLU(inplace=True),
nn.Dropout(0.5),
nn.Linear(128, 10)
)
self.adaptive_avg_pool2d = nn.AdaptiveAvgPool2d((1, 1))
self.relu = nn.ReLU(inplace=True)
def forward(self, x):
x = self.conv1(x)
# ResNet18-网络模型
x = self.block1_1(x)
x = self.se_block64(x)
x = self.block1_2(x)
x = self.block2_1(x)
x = self.se_block128(x)
x = self.block2_2(x)
x = self.block3_1(x)
x = self.se_block256(x)
x = self.block3_2(x)
x = self.block4_1(x)
x = self.se_block512(x)
x = self.block4_2(x)
# 平均值池化
x = self.adaptive_avg_pool2d(x)
# 数据平坦化处理，为接下来的全连接层做准备
x = x.view(x.size(0), -1)
x = self.FC_layer(x)
return x
# 初始化模型
model = ResNet18().to(device)
# 构造损失函数和优化器
criterion = nn.CrossEntropyLoss()
opt = optim.SGD(model.parameters(), lr=0.01, momentum=0.8, weight_decay=5e-4)
# 动态更新学习率------每隔step_size : lr = lr * gamma
schedule = optim.lr_scheduler.StepLR(opt, step_size=10, gamma=0.6, last_epoch=-1)
loss_list = []
# train
def train(epoch):
model.train()
start = time.time()
for epoch in range(epoch):
running_loss = 0.0
for i, (inputs, labels) in enumerate(train_loader, 0):
inputs, labels = inputs.to(device), labels.to(device)
# 将数据送入模型训练
outputs = model(inputs)
# 计算损失
loss = criterion(outputs, labels).to(device)
# 重置梯度
opt.zero_grad()
# 计算梯度，反向传播
loss.backward()
# 根据反向传播的梯度值优化更新参数
opt.step()
# 100个batch的 loss 之和
running_loss += loss.item()
loss_list.append(loss.item())
# 每一百个 batch 查看一下 平均loss
if (i + 1) % 100 == 0:
print('epoch = %d , batch = %d , loss = %.6f' % (epoch + 1, i + 1, running_loss / 100))
running_loss = 0.0
# 每一轮结束输出一下当前的学习率 lr
lr_1 = opt.param_groups[0]['lr']
print("learn_rate:%.15f" % lr_1)
schedule.step()
end = time.time()
# 计算并打印输出你的训练时间
print("time:{}".format(end - start))
# 训练过程可视化
plt.plot(loss_list)
plt.ylabel('loss')
plt.xlabel('Epoch')
plt.savefig('./SeResNet_train_img.png')
plt.show()
# Test
def verify():
model.eval()
correct = 0.0
total = 0
# 训练模式不需要反向传播更新梯度
with torch.no_grad():
print("===========================test===========================")
for inputs, labels in test_loader:
inputs, labels = inputs.to(device), labels.to(device)
outputs = model(inputs)
pred = outputs.argmax(dim=1) # 返回每一行中最大值元素索引
total += inputs.size(0)
correct += torch.eq(pred, labels).sum().item()
print("Accuracy of the network on the 10000 test images:%.2f %%" % (100 * correct / total))
print("==========================================================")
if __name__ == '__main__':
train(100)
verify()
# SeResNet: 采用se_block单元
# 使用 SeResNet 的神经网络训练 CIFAR10 数据集
