u2d_mag.py

import tensorflow as tf
from keras.layers import Input, Dense, BatchNormalization, Reshape, Conv2D, add, LeakyReLU, concatenate, Lambda, Dropout
from keras.models import Model
from keras.callbacks import TensorBoard, Callback
from keras import backend as K
from keras import optimizers
import scipy.io as sio
import numpy as np
import math
import time
import os
from keras.callbacks import ModelCheckpoint, LearningRateScheduler
from keras.callbacks import ReduceLROnPlateau
from keras.optimizers import Adam


envir = 'indoor'  # 'indoor' or 'outdoor'
# image params
img_height = 32
img_width = 32
img_channels = 1
img_total = img_height * img_width * img_channels
# network params
residual_num = 2
encoded_dim = img_total

def lr_schedule(epoch):
    """Learning Rate Schedule

    Learning rate is scheduled to be reduced after 80, 120, 160, 180 epochs.
    Called automatically every epoch as part of callbacks during training.

    # Arguments
        epoch (int): The number of epochs

    # Returns
        lr (float32): learning rate
    """
    step = 500
    lr = 1e-3
    if epoch > step+180:
        lr *= 0.5e-3
    elif epoch > step+160:
        lr *= 1e-3
    elif epoch > step+120:
        lr *= 5e-2
    elif epoch > step+80:
        lr *= 1e-1
    print('Learning rate: ', lr)
    return lr

class LossHistory(Callback):
    def on_train_begin(self, logs={}):
        self.losses_train = []
        self.losses_val = []

    def on_batch_end(self, batch, logs={}):
        self.losses_train.append(logs.get('loss'))

    def on_epoch_end(self, epoch, logs={}):
        self.losses_val.append(logs.get('val_loss'))

# Build the autoencoder model of u2d_mag
def residual_network(x, residual_num, encoded_dim):
    def add_common_layers(y):
        y = BatchNormalization()(y)
        y = LeakyReLU()(y)
        return y

    def residual_block_decoded(y):
        shortcut = y
        y = Conv2D(8, kernel_size=(3, 3), padding='same', data_format='channels_first')(y)
        y = add_common_layers(y)

        y = Conv2D(16, kernel_size=(3, 3), padding='same', data_format='channels_first')(y)
        y = add_common_layers(y)

        y = Conv2D(1, kernel_size=(3, 3), padding='same', data_format='channels_first')(y)
        y = BatchNormalization()(y)

        y = add([shortcut, y])
        y = LeakyReLU()(y)

        return y

    x = Conv2D(1, (3, 3), padding='same', data_format="channels_first")(x)
    x = add_common_layers(x)

    x = Reshape((img_total,))(x)
    x = Dropout(0.5)(x)
    encoded = Dense(encoded_dim, activation='linear')(x)
    x = Dropout(0.5)(encoded)

    x = Dense(img_total, activation='linear')(x)
    x = Reshape((img_channels, img_height, img_width,))(x)
    for i in range(residual_num):
        x = residual_block_decoded(x)

    x = Conv2D(1, (3, 3), activation='sigmoid', padding='same', data_format="channels_first")(x)

    return x


image_tensor = Input(shape=(img_channels, img_height, img_width))
auxiliary_input = Input(shape=(img_channels, img_height, img_width))
network_output = residual_network(image_tensor, residual_num, encoded_dim)
autoencoder = Model(inputs=[image_tensor], outputs=[network_output])
autoencoder.compile(optimizer='adam', loss='mse',metrics=['mse'])
print(autoencoder.summary())

# Data loading
if envir == 'indoor':
    mat = sio.loadmat('data/indoor53/Data100_Htrainin_mag.mat')
    x_train = mat['HD_train']
    x_train_up = mat['HU_train']
    mat = sio.loadmat('data/indoor53/Data100_Hvalin_mag.mat')
    x_val = mat['HD_val']
    x_val_up = mat['HU_val']

    x_test = x_val
    x_test_up = x_val_up

elif envir == 'outdoor':
    mat = sio.loadmat('data/urban3/Data100_Htrainin_mag.mat')
    x_train = mat['HD_train']
    x_train_up = mat['HU_train']
    mat = sio.loadmat('data/urban3/Data100_Hvalin_mag.mat')
    x_val = mat['HD_val']
    x_val_up = mat['HU_val']

    x_test = x_val
    x_test_up = x_val_up

x_train = x_train.astype('float32')
x_train_up = x_train_up.astype('float32')
x_val = x_val.astype('float32')
x_val_up = x_val_up.astype('float32')
x_test = x_test.astype('float32')
x_test_up = x_test_up.astype('float32')

x_train = np.reshape(x_train, (len(x_train), img_channels, img_height, img_width))  # adapt this if using `channels_first` image data format
x_train_up = np.reshape(x_train_up, (len(x_train_up), img_channels, img_height, img_width))
x_val = np.reshape(x_val, (len(x_val), img_channels, img_height, img_width))  # adapt this if using `channels_first` image data format
x_val_up = np.reshape(x_val_up, (len(x_val_up), img_channels, img_height, img_width))
x_test = np.reshape(x_test, (len(x_test), img_channels, img_height, img_width))  # adapt this if using `channels_first` image data format
x_test_up = np.reshape(x_test_up, (len(x_test_up), img_channels, img_height, img_width))





file = 'u2d_mag_' + (envir) + '_dim' + str(encoded_dim) + time.strftime('_%m_%d')
path = 'result_u2d/TensorBoard_%s/1' % file

# Prepare model model saving directory.
save_dir = os.path.join(os.getcwd(), 'result_u2d/saved_models_u2d_mag')
model_name = '%s_model.h5' % file
if not os.path.isdir(save_dir):
    os.makedirs(save_dir)
filepath = os.path.join(save_dir, model_name)

# Prepare callbacks for model saving and for learning rate adjustment.
checkpoint = ModelCheckpoint(filepath=filepath,
                             monitor='val_loss',
                             verbose=1,
                             save_best_only=True)

lr_scheduler = LearningRateScheduler(lr_schedule)

lr_reducer = ReduceLROnPlateau(factor=np.sqrt(0.1),
                               cooldown=0,
                               patience=5,
                               min_lr=0.5e-6)

history = LossHistory()

callbacks = [checkpoint, lr_reducer, lr_scheduler, history, TensorBoard(log_dir=path)]

autoencoder.fit([x_train_up], x_train,
                epochs=700,
                batch_size=200,
                shuffle=True,
                validation_data=([x_val_up], x_val),
                callbacks=callbacks)


autoencoder.load_weights(filepath)

# Testing data
tStart = time.time()
x_hat = autoencoder.predict([x_test_up])
tEnd = time.time()
print("It cost %f sec" % ((tEnd - tStart) / x_test.shape[0]))


x_test_mag = np.reshape(x_test[:, 0, :, :], (len(x_test), -1))
x_hat_mag = np.reshape(x_hat[:, 0, :, :], (len(x_hat), -1))

power = np.sum(abs(x_test_mag) ** 2, axis=1)
mse = np.sum(abs(x_test_mag - x_hat_mag) ** 2, axis=1)

# Here is the MSE for magnitude
print("In " + envir + " environment")
print("When dimension is", encoded_dim)
print("MSE is ", 10 * math.log10(np.mean(mse)))
filename = "result_u2d/decoded_%s.csv" % file
x_hat1 = np.reshape(x_hat, (len(x_hat), -1))
np.savetxt(filename, x_hat1, delimiter=",")