CGAN 和 DCGAN
import torch
import torch.nn as nn
import torch.optim as optim
from torchvision import datasets, transforms
import numpy as np
import matplotlib.pyplot as plt
# 基本參數
z_dim = 100
batch_size = 128
learning_rate = 0.0002
total_epochs = 30
# 使用GPU訓練,可以在菜單 "代碼執行工具" -> "更改運行時類型" 里進行設置
device = torch.device("cuda:0" if torch.cuda.is_available() else "cpu")
# 加載MNIST數據集
dataloader = torch.utils.data.DataLoader(
datasets.MNIST('./data', train=True, download=True,
transform=transforms.Compose([transforms.ToTensor(),transforms.Normalize((0.1307,), (0.3081,))])
), batch_size, shuffle=False, drop_last=True)
CGAN
-
下面分別是判別器和生成器的網絡結構可以看出網絡結構非常簡單,具體如下:
- 生成器:(784 + 10) ==> 512 ==> 256 ==> 1
- 判別器:(100 + 10) ==> 128 ==> 256 ==> 512 ==> 784
可以看出,去掉生成器和判別器那 10 維的標簽信息,和普通的GAN是完全一樣的。下面是網絡的具體實現代碼:
class Discriminator(nn.Module): '''全連接判別器,用於1x28x28的MNIST數據,輸出是數據和類別''' def __init__(self): super(Discriminator, self).__init__() self.model = nn.Sequential( nn.Linear(28*28+10, 512), nn.LeakyReLU(0.2, inplace=True), nn.Linear(512, 256), nn.LeakyReLU(0.2, inplace=True), nn.Linear(256, 1), nn.Sigmoid() ) def forward(self, x, c): x = x.view(x.size(0), -1) validity = self.model(torch.cat([x, c], -1)) return validity class Generator(nn.Module): '''全連接生成器,用於1x28x28的MNIST數據,輸入是噪聲和類別''' def __init__(self, z_dim): super(Generator, self).__init__() self.model = nn.Sequential( nn.Linear(z_dim+10, 128), nn.LeakyReLU(0.2, inplace=True), nn.Linear(128, 256), nn.BatchNorm1d(256, 0.8), nn.LeakyReLU(0.2, inplace=True), nn.Linear(256, 512), nn.BatchNorm1d(512, 0.8), nn.LeakyReLU(0.2, inplace=True), nn.Linear(in_features=512, out_features=28*28), nn.Tanh() ) def forward(self, z, c): x = self.model(torch.cat([z, c], dim=1)) x = x.view(-1, 1, 28, 28) return x
-
定義相關模型
# 初始化構建判別器和生成器 discriminator = Discriminator().to(device) generator = Generator(z_dim=z_dim).to(device) # 初始化二值交叉熵損失 bce = torch.nn.BCELoss().to(device) ones = torch.ones(batch_size).to(device) zeros = torch.zeros(batch_size).to(device) # 初始化優化器,使用Adam優化器 g_optimizer = optim.Adam(generator.parameters(), lr=learning_rate) d_optimizer = optim.Adam(discriminator.parameters(), lr=learning_rate)
-
開始訓練
# 開始訓練,一共訓練total_epochs for epoch in range(total_epochs): # torch.nn.Module.train() 指的是模型啟用 BatchNormalization 和 Dropout # torch.nn.Module.eval() 指的是模型不啟用 BatchNormalization 和 Dropout # 因此,train()一般在訓練時用到, eval() 一般在測試時用到 generator = generator.train() # 訓練一個epoch for i, data in enumerate(dataloader): # 加載真實數據 real_images, real_labels = data real_images = real_images.to(device) # 把對應的標簽轉化成 one-hot 類型 tmp = torch.FloatTensor(real_labels.size(0), 10).zero_() real_labels = tmp.scatter_(dim=1, index=torch.LongTensor(real_labels.view(-1, 1)), value=1) real_labels = real_labels.to(device) # 生成數據 # 用正態分布中采樣batch_size個隨機噪聲 z = torch.randn([batch_size, z_dim]).to(device) # 生成 batch_size 個 ont-hot 標簽 c = torch.FloatTensor(batch_size, 10).zero_() c = c.scatter_(dim=1, index=torch.LongTensor(np.random.choice(10, batch_size).reshape([batch_size, 1])), value=1) c = c.to(device) # 生成數據 fake_images = generator(z,c) # 計算判別器損失,並優化判別器 real_loss = bce(discriminator(real_images, real_labels), ones) fake_loss = bce(discriminator(fake_images.detach(), c), zeros) d_loss = real_loss + fake_loss d_optimizer.zero_grad() d_loss.backward() d_optimizer.step() # 計算生成器損失,並優化生成器 g_loss = bce(discriminator(fake_images, c), ones) g_optimizer.zero_grad() g_loss.backward() g_optimizer.step() # 輸出損失 print("[Epoch %d/%d] [D loss: %f] [G loss: %f]" % (epoch, total_epochs, d_loss.item(), g_loss.item()))
-
用隨機噪聲生成一組圖像,看看CGAN的效果:
#用於生成效果圖 # 生成100個隨機噪聲向量 fixed_z = torch.randn([100, z_dim]).to(device) # 生成100個one_hot向量,每類10個 fixed_c = torch.FloatTensor(100, 10).zero_() fixed_c = fixed_c.scatter_(dim=1, index=torch.LongTensor(np.array(np.arange(0, 10).tolist()*10).reshape([100, 1])), value=1) fixed_c = fixed_c.to(device) generator = generator.eval() fixed_fake_images = generator(fixed_z, fixed_c) plt.figure(figsize=(8, 8)) for j in range(10): for i in range(10): img = fixed_fake_images[j*10+i, 0, :, :].detach().cpu().numpy() img = img.reshape([28, 28]) plt.subplot(10, 10, j*10+i+1) plt.imshow(img, 'gray')
DCGAN
-
下面分別是判別器和生成器的網絡結構,和之前類似,只是使用了卷積結構。
class D_dcgan(nn.Module): '''滑動卷積判別器''' def __init__(self): super(D_dcgan, self).__init__() self.conv = nn.Sequential( # 第一個滑動卷積層,不使用BN,LRelu激活函數 nn.Conv2d(in_channels=1, out_channels=16, kernel_size=3, stride=2, padding=1), nn.LeakyReLU(0.2, inplace=True), # 第二個滑動卷積層,包含BN,LRelu激活函數 nn.Conv2d(in_channels=16, out_channels=32, kernel_size=3, stride=2, padding=1), nn.BatchNorm2d(32), nn.LeakyReLU(0.2, inplace=True), # 第三個滑動卷積層,包含BN,LRelu激活函數 nn.Conv2d(in_channels=32, out_channels=64, kernel_size=3, stride=2, padding=1), nn.BatchNorm2d(64), nn.LeakyReLU(0.2, inplace=True), # 第四個滑動卷積層,包含BN,LRelu激活函數 nn.Conv2d(in_channels=64, out_channels=128, kernel_size=4, stride=1), nn.BatchNorm2d(128), nn.LeakyReLU(0.2, inplace=True) ) # 全連接層+Sigmoid激活函數 self.linear = nn.Sequential(nn.Linear(in_features=128, out_features=1), nn.Sigmoid()) def forward(self, x): x = self.conv(x) x = x.view(x.size(0), -1) validity = self.linear(x) return validity class G_dcgan(nn.Module): '''反滑動卷積生成器''' def __init__(self, z_dim): super(G_dcgan, self).__init__() self.z_dim = z_dim # 第一層:把輸入線性變換成256x4x4的矩陣,並在這個基礎上做反卷機操作 self.linear = nn.Linear(self.z_dim, 4*4*256) self.model = nn.Sequential( # 第二層:bn+relu nn.ConvTranspose2d(in_channels=256, out_channels=128, kernel_size=3, stride=2, padding=0), nn.BatchNorm2d(128), nn.ReLU(inplace=True), # 第三層:bn+relu nn.ConvTranspose2d(in_channels=128, out_channels=64, kernel_size=3, stride=2, padding=1), nn.BatchNorm2d(64), nn.ReLU(inplace=True), # 第四層:不使用BN,使用tanh激活函數 nn.ConvTranspose2d(in_channels=64, out_channels=1, kernel_size=4, stride=2, padding=2), nn.Tanh() ) def forward(self, z): # 把隨機噪聲經過線性變換,resize成256x4x4的大小 x = self.linear(z) x = x.view([x.size(0), 256, 4, 4]) # 生成圖片 x = self.model(x) return x
-
定義相關模型
# 構建判別器和生成器 d_dcgan = D_dcgan().to(device) g_dcgan = G_dcgan(z_dim=z_dim).to(device) def weights_init_normal(m): classname = m.__class__.__name__ if classname.find('Conv') != -1: torch.nn.init.normal_(m.weight.data, 0.0, 0.02) elif classname.find('BatchNorm2d') != -1: torch.nn.init.normal_(m.weight.data, 1.0, 0.02) torch.nn.init.constant_(m.bias.data, 0.0) # 使用均值為0,方差為0.02的正態分布初始化神經網絡 d_dcgan.apply(weights_init_normal) g_dcgan.apply(weights_init_normal) # 初始化優化器,使用Adam優化器 g_dcgan_optim = optim.Adam(g_dcgan.parameters(), lr=learning_rate) d_dcgan_optim = optim.Adam(d_dcgan.parameters(), lr=learning_rate) # 加載MNIST數據集,和之前不同的是,DCGAN輸入的圖像被 resize 成 32*32 像素 dcgan_dataloader = torch.utils.data.DataLoader( datasets.MNIST('./data', train=True, download=True, transform=transforms.Compose([transforms.Resize(32), transforms.ToTensor(),transforms.Normalize((0.5,), (0.5,))]) ), batch_size, shuffle=True, drop_last=True)
-
開始訓練
# 開始訓練,一共訓練 total_epochs for e in range(total_epochs): # 給generator啟用 BatchNormalization g_dcgan = g_dcgan.train() # 訓練一個epoch for i, data in enumerate(dcgan_dataloader): # 加載真實數據,不加載標簽 real_images, _ = data real_images = real_images.to(device) # 用正態分布中采樣batch_size個噪聲,然后生成對應的圖片 z = torch.randn([batch_size, z_dim]).to(device) fake_images = g_dcgan(z) # 計算判別器損失,並優化判別器 real_loss = bce(d_dcgan(real_images), ones) fake_loss = bce(d_dcgan(fake_images.detach()), zeros) d_loss = real_loss + fake_loss d_dcgan_optim.zero_grad() d_loss.backward() d_dcgan_optim.step() # 計算生成器損失,並優化生成器 g_loss = bce(d_dcgan(fake_images), ones) g_dcgan_optim.zero_grad() g_loss.backward() g_dcgan_optim.step() # 輸出損失 print ("[Epoch %d/%d] [D loss: %f] [G loss: %f]" % (e, total_epochs, d_loss.item(), g_loss.item()))
-
用一組隨機噪聲輸出圖像,看看DCGAN的效果:
#用於生成效果圖 # 生成100個隨機噪聲向量 fixed_z = torch.randn([100, z_dim]).to(device) g_dcgan = g_dcgan.eval() fixed_fake_images = g_dcgan(fixed_z) plt.figure(figsize=(8, 8)) for j in range(10): for i in range(10): img = fixed_fake_images[j*10+i, 0, :, :].detach().cpu().numpy() img = img.reshape([32, 32]) plt.subplot(10, 10, j*10+i+1) plt.imshow(img, 'gray')