def main(hparams):
# Set up some stuff accoring to hparams
hparams.n_input = np.prod(hparams.image_shape)
utils.set_num_measurements(hparams)
utils.print_hparams(hparams)
# get inputs
data_dict = model_input(hparams)
estimator = utils.get_estimator(hparams, hparams.model_types[0])
print(estimator)
hparams.checkpoint_dir = utils.setup_checkpointing(hparams)
measurement_losses, l2_losses = utils.load_checkpoints(hparams)
h_hats_dict = {model_type: {} for model_type in hparams.model_types}
for key, x in data_dict.iteritems():
if not hparams.not_lazy:
# If lazy, first check if the image has already been
# saved before by *all* estimators. If yes, then skip this image.
save_paths = utils.get_save_paths(hparams, key)
is_saved = all([
os.path.isfile(save_path) for save_path in save_paths.values()
])
if is_saved:
continue
# Get Rx data
Rx = data_dict[key]['Rx_data']
Tx = data_dict[key]['Tx_data']
H = data_dict[key]['H_data']
Pilot_Rx = utils.get_pilot(Rx)
print('Pilot_shape', Pilot_Rx.shape)
Pilot_Rx = Pilot_Rx[0::2] + Pilot_Rx[1::2] * 1j
Pilot_Tx = utils.get_pilot(Tx)
Pilot_Tx = Pilot_Tx[0::2] + Pilot_Tx[1::2] * 1j
Pilot_complex = Pilot_Rx / Pilot_Tx
Pilot = np.empty((Pilot_complex.size * 2, ), dtype=Pilot_Rx.dtype)
Pilot[0::2] = np.real(Pilot_complex)
Pilot[1::2] = np.imag(Pilot_complex)
Pilot = np.reshape(Pilot, [1, -1]) / 2.5
# Construct estimates using each estimator
h_hat = estimator(Tx, Rx, Pilot, hparams)
# Compute and store measurement and l2 loss
# measurement_losses['dcgan'][key] = utils.get_measurement_loss(h_hat, Tx, Rx)
# l2_losses['dcgan'][key] = utils.get_l2_loss(h_hat, H)
print "Processed upto image {0} / {1}".format(key + 1, len(data_dict))
# Checkpointing
if (hparams.save_images) and ((key + 1) % hparams.checkpoint_iter
== 0):
# utils.checkpoint(key,h_hat, measurement_losses, l2_losses, save_image, hparams)
utils.save_channel_image(key + 1, h_hat, hparams)
utils.save_channel_mat(key + 1, h_hat, hparams)
print '\nProcessed and saved first ', key + 1, 'channels\n'
def main(hparams):
# Set up some stuff accoring to hparams
hparams.n_input = np.prod(hparams.image_shape)
utils.set_num_measurements(hparams)
utils.print_hparams(hparams)
# get inputs
data_dict = model_input(hparams)
estimator = utils.get_estimator(hparams, 'vae')
utils.setup_checkpointing(hparams)
measurement_losses, l2_losses = utils.load_checkpoints(hparams)
h_hats_dict = {model_type: {} for model_type in hparams.model_types}
for key, x in data_dict.iteritems():
if not hparams.not_lazy:
# If lazy, first check if the image has already been
# saved before by *all* estimators. If yes, then skip this image.
save_paths = utils.get_save_paths(hparams, key)
is_saved = all([
os.path.isfile(save_path) for save_path in save_paths.values()
])
if is_saved:
continue
# Get Rx data
Rx = data_dict[key]['Rx_data']
Tx = data_dict[key]['Rx_data']
H = data_dict[key]['H_data']
# Construct estimates using each estimator
h_hat = estimator(Tx, Rx, hparams)
# Save the estimate
h_hats_dict['vae'][key] = h_hat
# Compute and store measurement and l2 loss
measurement_losses['vae'][key] = utils.get_measurement_loss(
h_hat, Tx, Rx)
l2_losses['vae'][key] = utils.get_l2_loss(h_hat, H)
print 'Processed upto image {0} / {1}'.format(key + 1, len(data_dict))
# Checkpointing
if (hparams.save_images) and ((key + 1) % hparams.checkpoint_iter
== 0):
utils.checkpoint(key, h_hat, measurement_losses, l2_losses,
save_image, hparams)
print '\nProcessed and saved first ', key + 1, 'channels\n'
def main(hparams):
# Set up some stuff accoring to hparams
hparams.n_input = np.prod(hparams.image_shape)
utils.set_num_measurements(hparams)
utils.print_hparams(hparams)
# get inputs
xs_dict = model_input(hparams)
estimators = utils.get_estimators(hparams)
utils.setup_checkpointing(hparams)
measurement_losses, l2_losses = utils.load_checkpoints(hparams)
x_hats_dict = {model_type: {} for model_type in hparams.model_types}
x_batch_dict = {}
for key, x in xs_dict.iteritems():
if not hparams.not_lazy:
# If lazy, first check if the image has already been
# saved before by *all* estimators. If yes, then skip this image.
save_paths = utils.get_save_paths(hparams, key)
is_saved = all([
os.path.isfile(save_path) for save_path in save_paths.values()
])
if is_saved:
continue
x_batch_dict[key] = x
if len(x_batch_dict) < hparams.batch_size:
continue
# Reshape input
x_batch_list = [
x.reshape(1, hparams.n_input) for _, x in x_batch_dict.iteritems()
]
x_batch = np.concatenate(x_batch_list)
# Construct noise and measurements
A = utils.get_A(hparams)
noise_batch = hparams.noise_std * np.random.randn(
hparams.batch_size, hparams.num_measurements)
if hparams.measurement_type == 'project':
y_batch = x_batch + noise_batch
else:
y_batch = np.matmul(x_batch, A) + noise_batch
# Construct estimates using each estimator
for model_type in hparams.model_types:
estimator = estimators[model_type]
x_hat_batch = estimator(A, y_batch, hparams)
for i, key in enumerate(x_batch_dict.keys()):
x = xs_dict[key]
y = y_batch[i]
x_hat = x_hat_batch[i]
# Save the estimate
x_hats_dict[model_type][key] = x_hat
# Compute and store measurement and l2 loss
measurement_losses[model_type][
key] = utils.get_measurement_loss(x_hat, A, y)
l2_losses[model_type][key] = utils.get_l2_loss(x_hat, x)
print('Processed upto image {0} / {1}'.format(key + 1, len(xs_dict)))
# Checkpointing
if (hparams.save_images) and ((key + 1) % hparams.checkpoint_iter
== 0):
utils.checkpoint(x_hats_dict, measurement_losses, l2_losses,
save_image, hparams)
x_hats_dict = {
model_type: {}
for model_type in hparams.model_types
}
print('\nProcessed and saved first ', key + 1, 'images\n')
x_batch_dict = {}
# Final checkpoint
if hparams.save_images:
utils.checkpoint(x_hats_dict, measurement_losses, l2_losses,
save_image, hparams)
print('\nProcessed and saved all {0} image(s)\n'.format(len(xs_dict)))
if hparams.print_stats:
for model_type in hparams.model_types:
print(model_type)
mean_m_loss = np.mean(measurement_losses[model_type].values())
mean_l2_loss = np.mean(l2_losses[model_type].values())
print('mean measurement loss = {0}'.format(mean_m_loss))
print('mean l2 loss = {0}'.format(mean_l2_loss))
if hparams.image_matrix > 0:
utils.image_matrix(xs_dict, x_hats_dict, view_image, hparams)
# Warn the user that some things were not processsed
if len(x_batch_dict) > 0:
print(
'\nDid NOT process last {} images because they did not fill up the last batch.'
.format(len(x_batch_dict)))
print('Consider rerunning lazily with a smaller batch size.')
def main(hparams):
# if not hparams.use_gpu:
# os.environ['CUDA_VISIBLE_DEVICES'] = '-1'
# Set up some stuff accoring to hparams
hparams.n_input = np.prod(hparams.image_shape)
#hparams.stdv = 10 #adjust to HPARAM in model_def.py
#hparams.mean = 0 #adjust to HPARAM in model_def.py
utils.set_num_measurements(hparams)
utils.print_hparams(hparams)
hparams.bol = False
# hparams.dict_flag = False
# get inputs
if hparams.input_type == 'dict-input':# or hparams.dict_flag:
hparams_load_key = copy.copy(hparams)
hparams_load_key.input_type = 'full-input'
hparams_load_key.measurement_type = 'project'
hparams_load_key.zprior_weight = 0.0
hparams.key_field = np.load(utils.get_checkpoint_dir(hparams_load_key, hparams.model_types[0])+'candidates.npy').item()
print(hparams.measurement_type)
xs_dict, label_dict = model_input(hparams)
estimators = utils.get_estimators(hparams)
utils.setup_checkpointing(hparams)
sh = utils.SaveHandler()
sh.load_or_init_all(hparams.save_images,hparams.model_types,sh.get_pkl_filepaths(hparams,use_all=True))
if label_dict is None:
print('No labels exist.')
del sh.class_loss
# measurement_losses, l2_losses, emd_losses, x_orig, x_rec, noise_batch = utils.load_checkpoints(hparams)
if hparams.input_type == 'gen-span':
np.save(utils.get_checkpoint_dir(hparams, hparams.model_types[0])+'z.npy',hparams.z_from_gen)
np.save(utils.get_checkpoint_dir(hparams, hparams.model_types[0])+'images.npy',hparams.images_mat)
x_hats_dict = {model_type : {} for model_type in hparams.model_types}
x_batch_dict = {}
x_batch=[]
x_hat_batch=[]
# l2_losses2=np.zeros((len(xs_dict),1))
# distances_arr=[]
image_distance =np.zeros((len(xs_dict),1))
hparams.x = [] # TO REMOVE
for key, x in xs_dict.iteritems(): #//each batch once (x_batch_dict emptied at end)
if not hparams.not_lazy:
# If lazy, first check if the image has already been
# saved before by *all* estimators. If yes, then skip this image.
save_paths = utils.get_save_paths(hparams, key)
is_saved = all([os.path.isfile(save_path) for save_path in save_paths.values()])
if is_saved:
continue
x_batch_dict[key] = x
hparams.x.append(x)#To REMOVE
if len(x_batch_dict) < hparams.batch_size:
continue
# Reshape input
x_batch_list = [x.reshape(1, hparams.n_input) for _, x in x_batch_dict.iteritems()]
x_batch = np.concatenate(x_batch_list)
# x_batch, known_distortion, distances = get_random_distortion(x_batch)
# distances_arr[(key-1)*hparams.batch_size:key*hparams.batch_size] = distances
# xs_dict[(key-1)*hparams.batch_size:key*hparams.batch_size] =x_batch
# Construct noise and measurements
recovered, optim = utils.load_if_optimized(hparams)
if recovered and np.linalg.norm(optim.x_orig-x_batch) < 1e-10:
hparams.optim = optim
hparams.recovered = True
else:
hparams.recovered=False
optim.x_orig = x_batch
hparams.optim = optim
A, noise_batch, y_batch, c_val = utils.load_meas(hparams,sh,x_batch,xs_dict)
hparams.optim.noise_batch = noise_batch
if c_val:
continue
if hparams.measurement_type == 'sample_distribution':
plot_distribution(hparams,x_batch)
# for i in range(z.shape[1]):#range(1):
# plt.hist(z[i,:], facecolor='blue', alpha=0.5)
# directory_distr =
# pl.savefig("abc.png")
elif hparams.measurement_type == 'autoencoder':
plot_reconstruction(hparams,x_batch)
else:
# Construct estimates using each estimator
for model_type in hparams.model_types:
estimator = estimators[model_type]
start = time.time()
tmp = estimator(A, y_batch, hparams)
if isinstance(tmp,tuple):
x_hat_batch = tmp[0]
sh.z_rec = tmp[1]
else:
x_hat_batch = tmp
del sh.z_rec
end = time.time()
duration = end-start
print('The calculation needed {} time'.format(datetime.timedelta(seconds=duration)))
np.save(utils.get_checkpoint_dir(hparams, model_type)+'elapsed_time',duration)
# DEBUGGING = []
for i, key in enumerate(x_batch_dict.keys()):
# x = xs_dict[key]+known_distortion[i]
x = xs_dict[key]
y = y_batch[i]
x_hat = x_hat_batch[i]
# plt.figure()
# plt.imshow(np.reshape(x_hat, [64, 64, 3])*255)#, interpolation="nearest", cmap=plt.cm.gray)
# plt.show()
# Save the estimate
x_hats_dict[model_type][key] = x_hat
# Compute and store measurement and l2 loss
sh.measurement_losses[model_type][key] = utils.get_measurement_loss(x_hat, A, y)
# DEBUGGING.append(np.sum((x_hat.dot(A)-y)**2)/A.shape[1])
sh.l2_losses[model_type][key] = utils.get_l2_loss(x_hat, x)
if hparams.class_bol and label_dict is not None:
try:
sh.class_losses[model_type][key] = utils.get_classifier_loss(hparams,x_hat,label_dict[key])
except:
sh.class_losses[model_type][key] = NaN
warnings.warn('Class loss unsuccessfull, most likely due to corrupted memory. Simply retry.')
if hparams.emd_bol:
try:
_,sh.emd_losses[model_type][key] = utils.get_emd_loss(x_hat, x)
if 'nonneg' not in hparams.tv_or_lasso_mode and 'pca' in model_type:
warnings.warn('EMD requires nonnegative images, for safety insert nonneg into tv_or_lasso_mode')
except ValueError:
warnings.warn('EMD calculation unsuccesfull (most likely due to negative images)')
pass
# if l2_losses[model_type][key]-measurement_losses[model_type][key]!=0:
# print('NO')
# print(y)
# print(x)
# print(np.mean((x-y)**2))
image_distance[i] = np.linalg.norm(x_hat-x)
# l2_losses2[key] = np.mean((x_hat-x)**2)
# print('holla')
# print(l2_losses2[key])
# print(np.linalg.norm(x_hat-x)**2/len(xs_dict[0]))
# print(np.linalg.norm(x_hat-x)/len(xs_dict[0]))
# print(np.linalg.norm(x_hat-x))
print('Processed upto image {0} / {1}'.format(key+1, len(xs_dict)))
sh.x_orig = x_batch
sh.x_rec = x_hat_batch
sh.noise = noise_batch
#ACTIVATE ON DEMAND
#plot_bad_reconstruction(measurement_losses,x_batch)
# Checkpointing
if (hparams.save_images) and ((key+1) % hparams.checkpoint_iter == 0):
utils.checkpoint(x_hats_dict, save_image, sh, hparams)
x_hats_dict = {model_type : {} for model_type in hparams.model_types}
print('\nProcessed and saved first ', key+1, 'images\n')
x_batch_dict = {}
if 'wavelet' in hparams.model_types[0]:
print np.abs(sh.x_rec)
print('The average sparsity is {}'.format(np.sum(np.abs(sh.x_rec)>=0.0001)/float(hparams.batch_size)))
# Final checkpoint
if hparams.save_images:
utils.checkpoint(x_hats_dict, save_image, sh, hparams)
print('\nProcessed and saved all {0} image(s)\n'.format(len(xs_dict)))
if hparams.dataset in ['mnist', 'fashion-mnist']:
if np.array(x_batch).size:
utilsM.save_images(np.reshape(x_batch, [-1, 28, 28]),
[8, 8],utils.get_checkpoint_dir(hparams, hparams.model_types[0])+'original.png')
if np.array(x_hat_batch).size:
utilsM.save_images(np.reshape(x_hat_batch, [-1, 28, 28]),
[8, 8],utils.get_checkpoint_dir(hparams, hparams.model_types[0])+'reconstruction.png')
for model_type in hparams.model_types:
# print(model_type)
mean_m_loss = np.mean(sh.measurement_losses[model_type].values())
mean_l2_loss = np.mean(sh.l2_losses[model_type].values()) #\|XHUT-X\|**2/784/64
if hparams.emd_bol:
mean_emd_loss = np.mean(sh.emd_losses[model_type].values())
if label_dict is not None:
mean_class_loss = np.mean(sh.class_losses[model_type].values())
print('mean class loss = {0}'.format(mean_class_loss))
# print(image_distance)
mean_norm_loss = np.mean(image_distance)#sum_i(\|xhut_i-x_i\|)/64
# mean_rep_error = np.mean(distances_arr)
# mean_opt_meas_error_pixel = np.mean(np.array(l2_losses[model_type].values())-np.array(distances_arr)/xs_dict[0].shape)
# mean_opt_meas_error = np.mean(image_distance-distances_arr)
print('mean measurement loss = {0}'.format(mean_m_loss))
# print np.sum(np.asarray(DEBUGGING))/64
print('mean l2 loss = {0}'.format(mean_l2_loss))
if hparams.emd_bol:
print('mean emd loss = {0}'.format(mean_emd_loss))
print('mean distance = {0}'.format(mean_norm_loss))
print('mean distance pixelwise = {0}'.format(mean_norm_loss/len(xs_dict[xs_dict.keys()[0]])))
# print('mean representation error = {0}'.format(mean_rep_error))
# print('mean optimization plus measurement error = {0}'.format(mean_opt_meas_error))
# print('mean optimization plus measurement error per pixel = {0}'.format(mean_opt_meas_error_pixel))
if hparams.image_matrix > 0:
utils.image_matrix(xs_dict, x_hats_dict, view_image, hparams)
# Warn the user that some things were not processsed
if len(x_batch_dict) > 0:
print('\nDid NOT process last {} images because they did not fill up the last batch.'.format(len(x_batch_dict)))
print('Consider rerunning lazily with a smaller batch size.')
def main(hparams):
# Set up some stuff accoring to hparams
hparams.n_input = np.prod(hparams.image_shape)
utils.set_num_measurements(hparams)
utils.print_hparams(hparams)
if hparams.dataset == 'mnist':
hparams.n_z = latent_dim
elif hparams.dataset == 'celebA':
hparams.z_dim = latent_dim
# get inputs
xs_dict = model_input(hparams)
estimators = utils.get_estimators(hparams)
utils.setup_checkpointing(hparams)
measurement_losses, l2_losses = utils.load_checkpoints(hparams)
image_loss_mnist = []
meas_loss_mnist = []
x_hat_mnist = []
x_hats_dict = {model_type : {} for model_type in hparams.model_types}
x_batch_dict = {}
for key, x in xs_dict.iteritems():
if not hparams.not_lazy:
# If lazy, first check if the image has already been
# saved before by *all* estimators. If yes, then skip this image.
save_paths = utils.get_save_paths(hparams, key)
is_saved = all([os.path.isfile(save_path) for save_path in save_paths.values()])
if is_saved:
continue
x_batch_dict[key] = x
if len(x_batch_dict) < hparams.batch_size:
continue
# Reshape input
x_batch_list = [x.reshape(1, hparams.n_input) for _, x in x_batch_dict.iteritems()]
x_batch = np.concatenate(x_batch_list)
# Construct noise and measurements
A = utils.get_A(hparams)
noise_batch = hparams.noise_std * np.random.randn(hparams.batch_size, hparams.num_measurements)
if hparams.measurement_type == 'project':
y_batch = x_batch + noise_batch
else:
measure = np.matmul(x_batch, A)
y_batch = np.absolute(measure) + noise_batch
# Construct estimates using each estimator
for model_type in hparams.model_types:
x_main_batch = 10000*np.ones_like(x_batch)
for k in range(num_restarts):
print "Restart #", str(k+1)
# Solve deep pr problem with random initial iterate
init_iter = np.random.randn(hparams.batch_size, latent_dim)
# First gradient descent
z_opt_batch = init_iter
estimator = estimators[model_type]
items = estimator(A, y_batch, z_opt_batch, hparams)
x_hat_batch1 = items[0]
z_opt_batch1 = items[1]
losses_val1 = items[2]
x_hat_batch = x_hat_batch1
x_hat_batch = utils.resolve_ambiguity(x_hat_batch, x_batch, hparams.batch_size)
# Use reflection of initial iterate
z_opt_batch2 = -1*init_iter
items = estimator(A, y_batch, z_opt_batch2, hparams)
x_hat_batch2 = items[0]
z_opt_batch2 = items[1]
losses_val2 = items[2]
x_hat_batch2 = utils.resolve_ambiguity(x_hat_batch2, x_batch, hparams.batch_size)
x_hat_batchnew = utils.get_optimal_x_batch(x_hat_batch, x_hat_batch2, x_batch, hparams.batch_size)
x_main_batch = utils.get_optimal_x_batch(x_hat_batchnew, x_main_batch, x_batch, hparams.batch_size)
x_hat_batch = x_main_batch
if hparams.dataset == 'mnist':
utils.print_stats(x_hat_batch, x_batch, hparams.batch_size)
for i, key in enumerate(x_batch_dict.keys()):
x = xs_dict[key]
y = y_batch[i]
x_hat = x_hat_batch[i]
# Save the estimate
x_hats_dict[model_type][key] = x_hat
# Compute and store measurement and l2 loss
measurement_losses[model_type][key] = utils.get_measurement_loss(x_hat, A, y)
meas_loss_mnist.append(utils.get_measurement_loss(x_hat, A, y))
l2_losses[model_type][key] = utils.get_l2_loss(x_hat, x)
image_loss_mnist.append(utils.get_l2_loss(x_hat,x))
print 'Processed upto image {0} / {1}'.format(key+1, len(xs_dict))
# Checkpointing
if (hparams.save_images) and ((key+1) % hparams.checkpoint_iter == 0):
utils.checkpoint(x_hats_dict, measurement_losses, l2_losses, save_image, hparams)
x_hats_dict = {model_type : {} for model_type in hparams.model_types}
print '\nProcessed and saved first ', key+1, 'images\n'
x_batch_dict = {}
# Final checkpoint
if hparams.save_images:
utils.checkpoint(x_hats_dict, measurement_losses, l2_losses, save_image, hparams)
print '\nProcessed and saved all {0} image(s)\n'.format(len(xs_dict))
if hparams.print_stats:
for model_type in hparams.model_types:
mean_m_loss = np.mean(measurement_losses[model_type].values())
mean_l2_loss = np.mean(l2_losses[model_type].values())
print 'mean measurement loss = {0}'.format(mean_m_loss)
print 'mean l2 loss = {0}'.format(mean_l2_loss)
if hparams.image_matrix > 0:
utils.image_matrix(xs_dict, x_hats_dict, view_image, hparams)
# Warn the user that some things were not processsed
if len(x_batch_dict) > 0:
print '\nDid NOT process last {} images because they did not fill up the last batch.'.format(len(x_batch_dict))
print 'Consider rerunning lazily with a smaller batch size.'