DCGAN生成自己的图像

参考博客:https://github.com/lpty/tensorflow_tutorial/blob/master/avatarDcgan/avatar_gen.py

import os
import scipy.misc
import numpy as np
from glob import glob


class Avatar:

    def __init__(self):
        self.data_name = 'my_data'
        self.source_shape = (128, 128, 3)
        self.resize_shape = (32, 32, 3)
        self.crop = True
        self.img_shape = self.source_shape if not self.crop else self.resize_shape
        self.img_list = self._get_img_list()
        self.batch_size = 64
        self.batch_shape = (self.batch_size, ) + self.img_shape
        self.chunk_size = len(self.img_list) // self.batch_size

    def _get_img_list(self):
        path = os.path.join(os.getcwd(), self.data_name, '*.jpg')
        return glob(path)

    def _get_img(self, name):
        assert name in self.img_list
        img = scipy.misc.imread(name).astype(np.float32)
        assert img.shape == self.source_shape
        return self._resize(img) if self.crop else img

    def _resize(self, img):
        h, w = img.shape[:2]
        resize_h, resize_w = self.resize_shape[:2]
        crop_h, crop_w = self.source_shape[:2]
        j = int(round((h - crop_h) / 2.))
        i = int(round((w - crop_w) / 2.))
        cropped_image = scipy.misc.imresize(img[j:j + crop_h, i:i + crop_w], [resize_h, resize_w])
        return np.array(cropped_image) / 127.5 - 1.

    @staticmethod
    def save_img(image, path):
        scipy.misc.imsave(path, image)
        return True

    def batches(self):
        start = 0
        end = self.batch_size
        for _ in range(self.chunk_size):
            name_list = self.img_list[start:end]
            imgs = [self._get_img(name) for name in name_list]
            batches = np.zeros(self.batch_shape)
            batches[::] = imgs
            yield batches
            start += self.batch_size
            end += self.batch_size

if __name__ == '__main__':
    avatar = Avatar()
    batch = avatar.batches()
    b = next(batch)
    for num in range(len(b)):
        avatar.save_img(b[num], 'samples'+os.sep+str(num)+'.jpg')
import os
import math
import numpy as np
import tensorflow as tf
from datetime import datetime
 
from dcgan.avatar import Avatar
 
class AvatarModel:
 
    def __init__(self):
        self.avatar = Avatar()
        # 真实图片shape (height, width, depth)
        self.img_shape = self.avatar.img_shape
        # 一个batch的图片向量shape (batch, height, width, depth)
        self.batch_shape = self.avatar.batch_shape
        # 一个batch包含图片数量
        self.batch_size = self.avatar.batch_size
        # batch数量
        self.chunk_size = self.avatar.chunk_size
        self.mode = 'wgan-gp' #dcgan,wgan-gp
        # 噪音图片size
        self.noise_img_size = 100
        # 卷积转置输出通道数量
        self.gf_size = 64
        # 卷积输出通道数量
        self.df_size = 64
        # 训练循环次数
        self.epoch_size = 5
        # 学习率
        self.learning_rate = 0.0002
        # 优化指数衰减率
        self.beta1 = 0.5
        # 生成图片数量
        self.sample_size = 64
 
    @staticmethod
    def conv_out_size_same(size, stride):
        return int(math.ceil(float(size) / float(stride)))
 
    @staticmethod
    def linear(images, output_size, stddev=0.02, bias_start=0.0, name='Linear'):
        shape = images.get_shape().as_list()
 
        with tf.variable_scope(name):
            w = tf.get_variable("w", [shape[1], output_size], tf.float32,
                                tf.random_normal_initializer(stddev=stddev))
            b = tf.get_variable("b", [output_size],
                                initializer=tf.constant_initializer(bias_start))
            return tf.matmul(images, w) + b, w, b
 
    @staticmethod
    def batch_normalizer(x, epsilon=1e-5, momentum=0.9, train=True, name='batch_norm'):
        with tf.variable_scope(name):
            return tf.contrib.layers.batch_norm(x, decay=momentum, updates_collections=None, epsilon=epsilon,
                                                scale=True, is_training=train)
 
    @staticmethod
    def conv2d(images, output_dim, stddev=0.02, name="conv2d"):
        with tf.variable_scope(name):
            # filter : [height, width, in_channels, output_channels]
            # 注意与转置卷积的不同
            filter_shape = [5, 5, images.get_shape()[-1], output_dim]
            # strides
            # 对应每一维的filter移动步长
            strides_shape = [1, 2, 2, 1]
 
            w = tf.get_variable('w', filter_shape, initializer=tf.truncated_normal_initializer(stddev=stddev))
            b = tf.get_variable('b', [output_dim], initializer=tf.constant_initializer(0.0))
 
            conv = tf.nn.conv2d(images, w, strides=strides_shape, padding='SAME')
            conv = tf.reshape(tf.nn.bias_add(conv, b), conv.get_shape())
 
            return conv
 
    @staticmethod
    def deconv2d(images, output_shape, stddev=0.02, name='deconv2d'):
        with tf.variable_scope(name):
            # filter : [height, width, output_channels, in_channels]
            # 注意与卷积的不同
            filter_shape = [5, 5, output_shape[-1], images.get_shape()[-1]]
            # strides
            # 对应每一维的filter移动步长
            strides_shape = [1, 2, 2, 1]
 
            w = tf.get_variable('w', filter_shape, initializer=tf.random_normal_initializer(stddev=stddev))
            b = tf.get_variable('biases', [output_shape[-1]], initializer=tf.constant_initializer(0.0))
 
            deconv = tf.nn.conv2d_transpose(images, w, output_shape=output_shape, strides=strides_shape)
            deconv = tf.nn.bias_add(deconv, b)
 
            return deconv, w, b
 
    @staticmethod
    def lrelu(x, leak=0.2):
        return tf.maximum(x, leak * x)
 
    def generator(self, noise_imgs, train=True):
        with tf.variable_scope('generator'):
            # 分别对应每个layer的height, width
            s_h, s_w, _ = self.img_shape
            s_h2, s_w2 = self.conv_out_size_same(s_h, 2), self.conv_out_size_same(s_w, 2)
            s_h4, s_w4 = self.conv_out_size_same(s_h2, 2), self.conv_out_size_same(s_w2, 2)
            s_h8, s_w8 = self.conv_out_size_same(s_h4, 2), self.conv_out_size_same(s_w4, 2)
            s_h16, s_w16 = self.conv_out_size_same(s_h8, 2), self.conv_out_size_same(s_w8, 2)
 
            # layer 0
            # 对输入噪音图片进行线性变换
            z, h0_w, h0_b = self.linear(noise_imgs, self.gf_size*8*s_h16*s_w16)
            # reshape为合适的输入层格式
            h0 = tf.reshape(z, [-1, s_h16, s_w16, self.gf_size * 8])
            # 对数据进行归一化处理 加快收敛速度
            h0 = self.batch_normalizer(h0, train=train, name='g_bn0')
            # 激活函数
            h0 = tf.nn.relu(h0)
 
            # layer 1
            # 卷积转置进行上采样
            h1, h1_w, h1_b = self.deconv2d(h0, [self.batch_size, s_h8, s_w8, self.gf_size*4], name='g_h1')
            h1 = self.batch_normalizer(h1, train=train, name='g_bn1')
            h1 = tf.nn.relu(h1)
 
            # layer 2
            h2, h2_w, h2_b = self.deconv2d(h1, [self.batch_size, s_h4, s_w4, self.gf_size*2], name='g_h2')
            h2 = self.batch_normalizer(h2, train=train, name='g_bn2')
            h2 = tf.nn.relu(h2)
 
            # layer 3
            h3, h3_w, h3_b = self.deconv2d(h2, [self.batch_size, s_h2, s_w2, self.gf_size*1], name='g_h3')
            h3 = self.batch_normalizer(h3, train=train, name='g_bn3')
            h3 = tf.nn.relu(h3)
 
            # layer 4
            h4, h4_w, h4_b = self.deconv2d(h3, self.batch_shape, name='g_h4')
            return tf.nn.tanh(h4)
 
    def discriminator(self, real_imgs, reuse=False):
        with tf.variable_scope("discriminator", reuse=reuse):
            # layer 0
            # 卷积操作
            h0 = self.conv2d(real_imgs, self.df_size, name='d_h0_conv')
            # 激活函数
            h0 = self.lrelu(h0)
 
            # layer 1
            h1 = self.conv2d(h0, self.df_size*2, name='d_h1_conv')
            h1 = self.batch_normalizer(h1, name='d_bn1')
            h1 = self.lrelu(h1)
 
            # layer 2
            h2 = self.conv2d(h1, self.df_size*4, name='d_h2_conv')
            h2 = self.batch_normalizer(h2, name='d_bn2')
            h2 = self.lrelu(h2)
 
            # layer 3
            h3 = self.conv2d(h2, self.df_size*8, name='d_h3_conv')
            h3 = self.batch_normalizer(h3, name='d_bn3')
            h3 = self.lrelu(h3)
 
            # layer 4
            h4, _, _ = self.linear(tf.reshape(h3, [self.batch_size, -1]), 1, name='d_h4_lin')
 
            return tf.nn.sigmoid(h4), h4
 
#     @staticmethod
    def loss_graph(self,real_logits, fake_logits, fake_imgs, real_imgs):
        if self.mode == 'dcgan':
            # 生成器图片loss
            # 生成器希望判别器判断出来的标签为1
            gen_loss = tf.reduce_mean(tf.nn.sigmoid_cross_entropy_with_logits(logits=fake_logits, labels=tf.ones_like(fake_logits)))
            # 判别器识别生成器图片loss
            # 判别器希望识别出来的标签为0
            fake_loss = tf.reduce_mean(tf.nn.sigmoid_cross_entropy_with_logits(logits=fake_logits, labels=tf.zeros_like(fake_logits)))
            # 判别器识别真实图片loss
            # 判别器希望识别出来的标签为1
            real_loss = tf.reduce_mean(tf.nn.sigmoid_cross_entropy_with_logits(logits=real_logits, labels=tf.ones_like(real_logits)))
            # 判别器总loss
            dis_loss = tf.add(fake_loss, real_loss)
            clip_disc_weights = None
            
            return gen_loss, fake_loss, real_loss, dis_loss
        
        elif self.mode == 'wgan-gp':
            lamda = 10
            gen_loss = -tf.reduce_mean(fake_logits)
            dis_loss = tf.reduce_mean(fake_logits) - tf.reduce_mean(real_logits)
            print(fake_imgs.shape)
            print(real_imgs.shape)
            #优化器
            alpha = tf.random_uniform(shape=[self.batch_size,1,1,1],minval=0.,maxval=1.)  
            interpolates = alpha*fake_imgs + (1-alpha)*real_imgs
            print(interpolates.shape)
            _prob,_output = self.discriminator(interpolates,reuse=True)
            gradients = tf.gradients(_output,[interpolates,])[0]  
            slopes = tf.sqrt(tf.reduce_sum(tf.square(gradients),axis=[1,2,3]))  
            gradient_penalty = tf.reduce_mean((slopes-1.)**2)  
            dis_loss += lamda * gradient_penalty  
            
            clip_disc_weights = None
            
            return gen_loss, dis_loss
 
    @staticmethod
    def optimizer_graph(gen_loss, dis_loss, learning_rate, beta1):
        # 所有定义变量
        train_vars = tf.trainable_variables()
        # 生成器变量
        gen_vars = [var for var in train_vars if var.name.startswith('generator')]
        # 判别器变量
        dis_vars = [var for var in train_vars if var.name.startswith('discriminator')]
        # optimizer
        # 生成器与判别器作为两个网络需要分别优化
        gen_optimizer = tf.train.AdamOptimizer(learning_rate=learning_rate, beta1=beta1).minimize(gen_loss, var_list=gen_vars)
        dis_optimizer = tf.train.AdamOptimizer(learning_rate=learning_rate, beta1=beta1).minimize(dis_loss, var_list=dis_vars)
        return gen_optimizer, dis_optimizer
 
    def train(self):
        # 真实图片
        real_imgs = tf.placeholder(tf.float32, self.batch_shape, name='real_images')
        # 噪声图片
        noise_imgs = tf.placeholder(tf.float32, [None, self.noise_img_size], name='noise_images')
 
        # 生成器图片
        fake_imgs = self.generator(noise_imgs)
 
        # 判别器
        real_outputs, real_logits = self.discriminator(real_imgs)
        fake_outputs, fake_logits = self.discriminator(fake_imgs, reuse=True)
 
        # 损失
        if self.mode == 'dcgan':
            gen_loss, fake_loss, real_loss, dis_loss = self.loss_graph(real_logits, fake_logits,fake_imgs,real_imgs)
        elif self.mode == 'wgan-gp':
            gen_loss, dis_loss = self.loss_graph(real_logits, fake_logits,fake_imgs,real_imgs)
        # 优化
        gen_optimizer, dis_optimizer = self.optimizer_graph(gen_loss, dis_loss, self.learning_rate, self.beta1)
 
        # 开始训练
        saver = tf.train.Saver()
        step = 0
        # 指定占用GPU比例
        # tensorflow默认占用全部GPU显存 防止在机器显存被其他程序占用过多时可能在启动时报错
        gpu_options = tf.GPUOptions(per_process_gpu_memory_fraction=0.8)
        with tf.Session(config=tf.ConfigProto(gpu_options=gpu_options)) as sess:
            sess.run(tf.global_variables_initializer())
            for epoch in range(self.epoch_size):
                batches = self.avatar.batches()
                for batch_imgs in batches:
                    # generator的输入噪声
                    noises = np.random.uniform(-1, 1, size=(self.batch_size, self.noise_img_size)).astype(np.float32)
                    if self.mode == 'dcgan':
                        # 优化
                        _ = sess.run(dis_optimizer, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                        _ = sess.run(gen_optimizer, feed_dict={noise_imgs: noises})
                        _ = sess.run(gen_optimizer, feed_dict={noise_imgs: noises})
                    if self.mode == 'wgan-gp':
                        for i in range(5):
                            _ = sess.run(dis_optimizer, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                        _ = sess.run(gen_optimizer, feed_dict={noise_imgs: noises})
                        
                    step += 1
                    print(datetime.now().strftime('%c'), epoch, step)
                    
                if self.mode == 'dcgan':
                    # 每一轮结束计算loss
                    # 判别器损失
                    loss_dis = sess.run(dis_loss, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                    # 判别器对真实图片
                    loss_real = sess.run(real_loss, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                    # 判别器对生成器图片
                    loss_fake = sess.run(fake_loss, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                    # 生成器损失
                    loss_gen = sess.run(gen_loss, feed_dict={noise_imgs: noises})
                    
                    print(datetime.now().strftime('%c'), ' epoch:', epoch, ' step:', step, ' loss_dis:', loss_dis,
                          ' loss_real:', loss_real, ' loss_fake:', loss_fake, ' loss_gen:', loss_gen)
                else:
                    loss_dis = sess.run(dis_loss, feed_dict={real_imgs: batch_imgs, noise_imgs: noises})
                    loss_gen = sess.run(gen_loss, feed_dict={noise_imgs: noises})
                    
                    print(datetime.now().strftime('%c'), ' epoch:', epoch, ' step:', step, ' loss_dis:', loss_dis,
                           ' loss_gen:', loss_gen)
            saver.save(sess,'Save_wgan/dcgan.ckpt' )
#             model_path = os.getcwd() + os.sep + "avatar.model"
#             saver.save(sess, model_path, global_step=step)
 
    def gen(self):
        # 生成图片
        noise_imgs = tf.placeholder(tf.float32, [None, self.noise_img_size], name='noise_imgs')
        sample_imgs = self.generator(noise_imgs, train=False)
        saver = tf.train.Saver()
        with tf.Session() as sess:
            sess.run(tf.global_variables_initializer())
            #saver.restore(sess, tf.train.latest_checkpoint("Save_wgan"))
            saver.restore(sess,'Save_wgan/wgan20.ckpt')
            sample_noise = np.random.uniform(-1, 1, size=(self.sample_size, self.noise_img_size))
            samples = sess.run(sample_imgs, feed_dict={noise_imgs: sample_noise})
        for num in range(len(samples)):
            self.avatar.save_img(samples[num], 'samples'+os.sep+str(num)+'.jpg')
if __name__ == '__main__':
    avatar = AvatarModel()
    avatar.train()
#     avatar.gen()

 

你可能感兴趣的:(Gan,深度学习)