def classifier_one(dataset, model_name, x_target, y_target, len_x_target):
NUM_CLASSES = 10
fp = open("classifier_result.txt", 'a')
print("dataset: {}; model name: {} Evaluation start.".format(
dataset, model_name))
_, X_dim, width, height, channels, len_x_train, x_train, y_train, len_x_test, x_test, y_test = data_loader(
dataset)
graph = tf.Graph()
with graph.as_default():
session_conf = tf.ConfigProto(allow_soft_placement=True,
log_device_placement=False)
sess = tf.Session(config=session_conf)
with sess.as_default():
#input placeholder
input_shape = [None, width, height, channels]
filter_sizes = [5, 5, 5, 5, 5]
hidden = 128
n_filters = [channels, hidden, hidden * 2, hidden * 4]
X = tf.placeholder(tf.float32,
shape=[None, width, height, channels])
Y = tf.placeholder(tf.float32, shape=[None, NUM_CLASSES])
#discriminator variables
W1 = tf.Variable(he_normal_init([5, 5, channels, hidden // 2]))
W2 = tf.Variable(he_normal_init([5, 5, hidden // 2, hidden]))
W3 = tf.Variable(he_normal_init([5, 5, hidden, hidden * 2]))
W4 = tf.Variable(xavier_init([4 * 4 * hidden * 2, NUM_CLASSES]))
b4 = tf.Variable(tf.zeros(shape=[NUM_CLASSES]))
var_C = [W1, W2, W3, W4, b4]
global_step = tf.Variable(0, name="global_step", trainable=False)
C_real_logits = classifier(X, var_C)
C_loss = tf.reduce_mean(
tf.nn.softmax_cross_entropy_with_logits_v2(
labels=Y, logits=C_real_logits))
num_batches_per_epoch = int((len_x_target - 1) / mb_size) + 1
C_solver = tf.contrib.opt.AdamWOptimizer(
weight_decay=1e-4, learning_rate=1e-3, beta1=0.5,
beta2=0.9).minimize(C_loss,
var_list=var_C,
global_step=global_step)
timestamp = str(int(time.time()))
sess.run(tf.global_variables_initializer())
x_temp = np.append(x_test, x_test[:mb_size], axis=0)
y_temp = np.append(y_test, y_test[:mb_size], axis=0)
best_accuracy = 0.0
for it in range(num_batches_per_epoch * 500):
X_mb, Y_mb = next_batch(mb_size, x_target, y_target)
_, C_curr = sess.run([C_solver, C_loss],
feed_dict={
X: X_mb,
Y: Y_mb
})
if it % 100 == 0:
print('Iter: {}; C_loss: {:.4};'.format(it, C_curr))
if it % 1000 == 0:
predictions = []
for jt in range(len_x_test // mb_size + 1):
Xt_mb, Yt_mb = next_test_batch(jt, mb_size, x_temp,
y_temp)
_, C_pred = sess.run([C_solver, C_real_logits],
feed_dict={
X: Xt_mb,
Y: Yt_mb
})
if len(predictions) == 0:
predictions = C_pred
else:
predictions = np.append(predictions,
C_pred,
axis=0)
predictions = predictions[:len_x_test]
predictions = np.argmax(predictions, axis=1)
correct_y = np.argmax(y_test, axis=1)
correct_predictions = sum(predictions == correct_y)
accuracy = correct_predictions / float(len_x_test)
print('Iter: {}; accuracy: {:.4}; best accuracy: {:.4}'.
format(it, accuracy, best_accuracy))
if accuracy > best_accuracy:
best_accuracy = accuracy
print("dataset: {} model name: {} with best accuracy: {:.4} fin.".
format(dataset, model_name, best_accuracy))
print("dataset: {} model name: {} with best accuracy: {:.4} fin.".
format(dataset, model_name, best_accuracy),
file=fp)
fp.close()
sess.close()
return
if dataset == 'celebA' or dataset == 'lsun':
n_filters = [channels, hidden, hidden * 2, hidden * 4, hidden * 8]
else:
n_filters = [channels, hidden, hidden * 2, hidden * 4]
X = tf.placeholder(tf.float32, shape=[None, width, height, channels])
Z_noise = tf.placeholder(tf.float32, shape=[None, z_dim])
A_true_flat = X
#autoencoder variables
var_G = []
var_H = []
var_A = []
#discriminator variables
W1 = tf.Variable(he_normal_init([5, 5, channels, hidden // 2]))
W2 = tf.Variable(he_normal_init([5, 5, hidden // 2, hidden]))
W3 = tf.Variable(he_normal_init([5, 5, hidden, hidden * 2]))
if dataset == 'celebA' or dataset == 'lsun':
W4 = tf.Variable(he_normal_init([5, 5, hidden * 2, hidden * 4]))
W5 = tf.Variable(xavier_init([4 * 4 * hidden * 4, 1]))
b5 = tf.Variable(tf.zeros(shape=[1]))
var_D = [W1, W2, W3, W4, W5, b5]
else:
W4 = tf.Variable(xavier_init([4 * 4 * hidden * 2, 1]))
b4 = tf.Variable(tf.zeros(shape=[1]))
var_D = [W1, W2, W3, W4, b4]
global_step = tf.Variable(0, name="global_step", trainable=False)
G_sample, A_sample, z_original, z_noised, epsilon_layer, delta_layer, z_noise = eddp_autoencoder(
def generate_one(dataset,model_name, z_dim,USE_DELTA):
NUM_CLASSES = 10
mb_size, X_dim, width, height, channels,len_x_train, x_train, y_train, len_x_test, x_test, y_test = data_loader(dataset)
graph = tf.Graph()
with graph.as_default():
session_conf = tf.ConfigProto(allow_soft_placement=True, log_device_placement=False)
sess = tf.Session(config=session_conf)
with sess.as_default():
#input placeholder
input_shape=[None, width, height, channels]
filter_sizes=[5, 5, 5, 5, 5]
hidden = 128
n_filters=[channels, hidden, hidden*2, hidden*4]
X = tf.placeholder(tf.float32, shape=[None, width, height,channels])
Z_S = tf.placeholder(tf.float32, shape=[None, z_dim])
Z_noise = tf.placeholder(tf.float32, shape=[None, z_dim])
Y = tf.placeholder(tf.float32, shape=[None, NUM_CLASSES])
A_true_flat = X
#generator variables
var_G = []
#discriminator variables
W1 = tf.Variable(he_normal_init([5,5,channels, hidden//2]))
W2 = tf.Variable(he_normal_init([5,5, hidden//2,hidden]))
W3 = tf.Variable(he_normal_init([5,5,hidden,hidden*2]))
W4 = tf.Variable(xavier_init([4*4*hidden*2, 1]))
b4 = tf.Variable(tf.zeros(shape=[1]))
var_D = [W1,W2,W3,W4,b4]
#classifier variables
W4_c = tf.Variable(xavier_init([4*4*hidden*2, NUM_CLASSES]))
b4_c = tf.Variable(tf.zeros(shape=[NUM_CLASSES]))
var_C = [W1,W2,W3,W4_c,b4_c]
var_D_C = [W1,W2,W3,W4,b4,W4_c,b4_c]
global_step = tf.Variable(0, name="global_step", trainable=False)
G_sample, G_zero, z_original, z_noise, z_noised = generator(input_shape,
n_filters,
filter_sizes,
X,
Z_noise,
var_G,
z_dim,
Z_S,
USE_DELTA)
D_real, D_real_logits = discriminator(X, var_D)
D_fake, D_fake_logits = discriminator(G_sample, var_D)
C_real_logits = classifier(X, var_C)
C_fake_logits = classifier(G_sample, var_C)
D_real_loss = tf.nn.sigmoid_cross_entropy_with_logits(
logits=D_real_logits, labels=tf.ones_like(D_real))
D_fake_loss = tf.nn.sigmoid_cross_entropy_with_logits(
logits=D_fake_logits, labels=tf.zeros_like(D_fake))
D_S_loss = tf.reduce_mean(D_real_loss) + tf.reduce_mean(D_fake_loss)
D_C_loss = tf.reduce_mean(
tf.nn.softmax_cross_entropy_with_logits_v2(labels = Y,
logits = C_real_logits))
G_zero_loss = tf.reduce_mean(tf.pow(X - G_zero,2))
G_S_loss = tf.reduce_mean(
tf.nn.sigmoid_cross_entropy_with_logits(logits=D_fake_logits,
labels=tf.ones_like(D_fake)))
G_C_loss = tf.reduce_mean(
tf.nn.softmax_cross_entropy_with_logits_v2(labels = Y,
logits = C_fake_logits))
D_loss = D_S_loss + D_C_loss
G_loss = G_S_loss + G_C_loss + G_zero_loss
#sensitivity estimation
latent_max = tf.reduce_max(z_original, axis = 0)
latent_min = tf.reduce_min(z_original, axis = 0)
tf.summary.image('Original',X)
tf.summary.image('fake',G_sample)
tf.summary.image('fake_zero', G_zero)
tf.summary.scalar('D_loss', D_loss)
tf.summary.scalar('D_S_loss',D_S_loss)
tf.summary.scalar('D_C_loss',D_C_loss)
tf.summary.scalar('G_zero_loss',G_zero_loss)
tf.summary.scalar('G_S_loss',G_S_loss)
tf.summary.scalar('G_C_loss',G_C_loss)
tf.summary.scalar('G_loss',G_loss)
tf.summary.histogram('z_original',z_original)
tf.summary.histogram('z_noise',z_noise)
tf.summary.histogram('z_noised',z_noised)
merged = tf.summary.merge_all()
num_batches_per_epoch = int((len_x_train-1)/mb_size) + 1
A_solver = tf.contrib.opt.AdamWOptimizer(weight_decay=1e-4,learning_rate=1e-4,beta1=0.5, beta2=0.9).minimize(G_zero_loss,var_list=var_G, global_step=global_step)
D_solver = tf.contrib.opt.AdamWOptimizer(weight_decay=1e-4,learning_rate=1e-4,beta1=0.5, beta2=0.9).minimize(D_loss,var_list=var_D_C, global_step=global_step)
G_solver = tf.contrib.opt.AdamWOptimizer(weight_decay=1e-4,learning_rate=1e-4,beta1=0.5, beta2=0.9).minimize(G_loss,var_list=var_G, global_step=global_step)
timestamp = str(int(time.time()))
out_dir = os.path.abspath(os.path.join(os.path.curdir,
"results/models/"+ model_name))
checkpoint_dir = os.path.abspath(os.path.join(out_dir, "checkpoints"))
checkpoint_prefix = os.path.join(checkpoint_dir, "model")
if not os.path.exists("results/generated/"):
os.makedirs("results/generated/")
saver = tf.train.Saver(tf.global_variables())
sess.run(tf.global_variables_initializer())
saver.restore(sess,tf.train.latest_checkpoint(checkpoint_dir))
#calculate approximated global sensitivity
for idx in range(num_batches_per_epoch):
X_mb, Y_mb = next_batch(mb_size, x_train, y_train)
enc_zero = np.zeros([mb_size,z_dim]).astype(np.float32)
if USE_DELTA:
enc_noise = np.random.normal(0.0,1.0,[mb_size,z_dim]).astype(np.float32)
else:
enc_noise = np.random.laplace(0.0,1.0,[mb_size,z_dim]).astype(np.float32)
max_curr, min_curr = sess.run([latent_max,latent_min], feed_dict={
X: X_mb,
Y: Y_mb,
Z_noise: enc_zero,
Z_S: enc_zero})
if idx == 0:
z_max = max_curr
z_min = min_curr
else:
z_max = np.maximum(z_max,max_curr)
z_min = np.minimum(z_min,min_curr)
z_sensitivity = np.abs(np.subtract(z_max,z_min))
#print("Approximated Global Sensitivity:")
#print(z_sensitivity)
z_sensitivity = np.tile(z_sensitivity,(mb_size,1))
x_train = np.append(x_train, X_mb, axis=0)
y_train = np.append(y_train, Y_mb, axis=0)
for i in range(num_batches_per_epoch):
X_mb, Y_mb = next_test_batch(i, mb_size, x_train, y_train)
enc_zero = np.zeros([mb_size,z_dim]).astype(np.float32)
if USE_DELTA:
enc_noise = np.random.normal(0.0,1.0,[mb_size,z_dim]).astype(np.float32)
else:
enc_noise = np.random.laplace(0.0,1.0,[mb_size,z_dim]).astype(np.float32)
G_sample_curr = sess.run(G_sample,
feed_dict={X: X_mb,
Y: Y_mb,
Z_noise: enc_noise,
Z_S: z_sensitivity})
samples_flat = tf.reshape(G_sample_curr,[mb_size,width,height,channels]).eval()
if i == 0:
img_set = samples_flat
label_set = Y_mb
else:
img_set = np.append(img_set, samples_flat, axis=0)
label_set = np.append(label_set, Y_mb, axis=0)
x_generated = img_set[:len_x_train]
y_generated = label_set[:len_x_train]
outfile = "results/generated/{}".format(model_name)
np.savez(outfile, x=x_generated, y=y_generated)
print("dataset: {} model name: {} fin.".format(dataset, model_name))
print("dataset: {} model name: {} fin.".format(dataset, model_name), file = fp)
sess.close()
return x_train, img_set
