def test_missing_inputs():
missing_idx = 10
class TimeOutSequence(DummySequence):
def __getitem__(self, item):
if item == missing_idx:
time.sleep(120)
return super(TimeOutSequence, self).__getitem__(item)
enqueuer = GeneratorEnqueuer(create_finite_generator_from_sequence_pcs(
TimeOutSequence([3, 2, 2, 3])),
use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.warns(UserWarning, match='An input could not be retrieved.'):
for _ in range(4 * missing_idx):
next(gen_output)
enqueuer = OrderedEnqueuer(TimeOutSequence([3, 2, 2, 3]),
use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
warning_msg = "The input {} could not be retrieved.".format(missing_idx)
with pytest.warns(UserWarning, match=warning_msg):
for _ in range(11):
next(gen_output)
def test_missing_inputs():
missing_idx = 10
class TimeOutSequence(DummySequence):
def __getitem__(self, item):
if item == missing_idx:
time.sleep(120)
return super(TimeOutSequence, self).__getitem__(item)
enqueuer = GeneratorEnqueuer(create_finite_generator_from_sequence_pcs(
TimeOutSequence([3, 2, 2, 3])), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.warns(UserWarning, match='An input could not be retrieved.'):
for _ in range(4 * missing_idx):
next(gen_output)
enqueuer = OrderedEnqueuer(TimeOutSequence([3, 2, 2, 3]),
use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
warning_msg = "The input {} could not be retrieved.".format(missing_idx)
with pytest.warns(UserWarning, match=warning_msg):
for _ in range(11):
next(gen_output)
class MultiProcessIterableDataIterator(BaseDataIterator):
def __init__(self, data_loader):
super().__init__(data_loader)
generator = _MPAssistantIterator(self.dataset)
self.enqueuer = GeneratorEnqueuer(generator, use_multiprocessing=True)
self.enqueuer.start(workers=self.num_workers,
max_queue_size=self.queue_size)
self.output_generator = self.enqueuer.get()
self.buffer = []
def __next__(self):
data = []
batch_indices = self._next_index()
for _ in batch_indices:
try:
while len(self.buffer) < self.queue_size:
output = next(self.output_generator)
self.buffer.append(output)
except StopIteration:
pass
if self.buffer:
i = random.randrange(len(self.buffer))
data.append(self.buffer.pop(i))
if len(data) == 0 or (self.drop_last
and len(data) < len(batch_indices)):
raise StopIteration
return self.collator.collate_fn(data)
def compute_prior_dist(data_generator, cf):
prior_path = os.path.join(cf.output_path, 'priors')
if not os.path.exists(prior_path): os.makedirs(prior_path)
test_datagen = GeneratorEnqueuer(
data_generator, use_multiprocessing=cf.use_multiprocessing)
test_datagen.start(cf.workers, cf.max_queue_size)
test_generator = test_datagen.get()
a_cum_dist = np.zeros(219)
b_cum_dist = np.zeros(219)
for _ in tqdm.tqdm(range(data_generator.samples)):
batch_real, _ = next(test_generator)
rgb_batch = data_generator.decoder(batch_real)
for image in rgb_batch:
image = rgb2lab(image)
a_cum_dist += np.histogram(image[..., 1],
bins=np.arange(-110., 110.))[0]
b_cum_dist += np.histogram(image[..., 2],
bins=np.arange(-110., 110.))[0]
a_dist = np.log(a_cum_dist / data_generator.samples)
b_dist = np.log(b_cum_dist / data_generator.samples)
a_dist[a_dist == -np.inf] = np.inf
a_dist[a_dist == np.inf] = a_dist.min()
b_dist[b_dist == -np.inf] = np.inf
b_dist[b_dist == np.inf] = b_dist.min()
np.save(os.path.join(prior_path, 'a_dist.npy'), a_dist)
np.save(os.path.join(prior_path, 'b_dist.npy'), b_dist)
def test_generator_enqueuer_fail_processes():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
FaultSequence()), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(StopIteration):
next(gen_output)
def DatasetIterator(dataset_path,
size,
iter_sampler,
use_multiprocessing=False,
workers=1,
max_queue_size=10):
'''Iterate over the dataset
# Arguments
dataset_path: path to `.h5` dataset
size: batch size
iter_sampler: callable. Should return iteration to be sliced
use_multiprocessing: use multiprocessing for workers, Bool
workers: number of workers, int
max_queue_size: maximum queue size, int
# Output
x_batch, y_batch
x_batch: input tensor of size `(size, height, width, 2)`
y_batch: output mask of size `(size, height, width, 1)`
'''
core_gen = minibatch_generator(dataset_path=dataset_path,
size=size,
iter_sampler=iter_sampler)
enqueuer = GeneratorEnqueuer(generator=core_gen,
use_multiprocessing=use_multiprocessing)
enqueuer.start(workers, max_queue_size=max_queue_size)
generator = enqueuer.get()
while True:
try:
yield next(generator)
except StopIteration:
return
def test_generator_enqueuer_fail_processes():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
FaultSequence()), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(IndexError):
next(gen_output)
def _predict(self,
data_generator_function,
steps_per_epoch,
include_datum=True):
data_generator = data_generator_function(include_datum=True)
enqueuer = GeneratorEnqueuer(data_generator)
enqueuer.start(workers=self._WORKERS, max_queue_size=self._MAX_Q_SIZE)
caption_results = []
datum_results = []
for _ in tqdm(range(steps_per_epoch)):
generator_output = None
while enqueuer.is_running():
if not enqueuer.queue.empty():
generator_output = enqueuer.get()
break
else:
sleep(self._WAIT_TIME)
X, y, datum_batch = next(generator_output)
captions_pred_str = self._predict_batch(X, y)
caption_results += captions_pred_str
datum_results += datum_batch
enqueuer.stop()
if include_datum:
return zip(caption_results, datum_results)
else:
return caption_results
def DatasetIterator(dataset_path, size, iter_sampler,
use_multiprocessing=False,
workers=1, max_queue_size=10):
'''Iterate over the dataset
# Arguments
dataset_path: path to `.h5` dataset
size: batch size
iter_sampler: callable. Should return iteration to be sliced
use_multiprocessing: use multiprocessing for workers, Bool
workers: number of workers, int
max_queue_size: maximum queue size, int
# Output
x_batch, y_batch
x_batch: input tensor of size `(size, height, width, 2)`
y_batch: output mask of size `(size, height, width, 1)`
'''
core_gen = minibatch_generator(dataset_path=dataset_path, size=size, iter_sampler=iter_sampler)
enqueuer = GeneratorEnqueuer(generator=core_gen, use_multiprocessing=use_multiprocessing)
enqueuer.start(workers, max_queue_size=max_queue_size)
generator = enqueuer.get()
while True:
try:
yield next(generator)
except StopIteration:
return
def test_generator_enqueuer_fail_threads():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_threads(
FaultSequence()), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(StopIteration):
next(gen_output)
def test_generator_enqueuer_fail_threads():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_threads(
FaultSequence()), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(IndexError):
next(gen_output)
def test_generator_enqueuer_threadsafe():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
DummySequence([3, 10, 10, 3])), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(RuntimeError) as e:
[next(gen_output) for _ in range(10)]
assert 'thread-safe' in str(e.value)
enqueuer.stop()
def test_generator_enqueuer_threadsafe():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
DummySequence([3, 10, 10, 3])), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
with pytest.raises(RuntimeError) as e:
[next(gen_output) for _ in range(10)]
assert 'thread-safe' in str(e.value)
enqueuer.stop()
def test_finite_generator_enqueuer_threads():
enqueuer = GeneratorEnqueuer(create_finite_generator_from_sequence_threads(
DummySequence([3, 200, 200, 3])), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for output in gen_output:
acc.append(int(output[0, 0, 0, 0]))
assert set(acc) == set(range(100)), "Output is not the same"
enqueuer.stop()
def test_finite_generator_enqueuer_processes():
enqueuer = GeneratorEnqueuer(create_finite_generator_from_sequence_pcs(
TestSequence([3, 200, 200, 3])), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for output in gen_output:
acc.append(int(output[0, 0, 0, 0]))
assert acc != list(range(100)), "Order was keep in GeneratorEnqueuer with processes"
enqueuer.stop()
def test_finite_generator_enqueuer_threads():
enqueuer = GeneratorEnqueuer(create_finite_generator_from_sequence_threads(
TestSequence([3, 200, 200, 3])), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for output in gen_output:
acc.append(int(output[0, 0, 0, 0]))
assert len(set(acc) - set(range(100))) == 0, "Output is not the same"
enqueuer.stop()
def test_generator_enqueuer_processes():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
DummySequence([3, 200, 200, 3])), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for i in range(100):
acc.append(int(next(gen_output)[0, 0, 0, 0]))
assert acc != list(range(100)), ('Order was keep in GeneratorEnqueuer '
'with processes')
enqueuer.stop()
def train_gan(gan, data_generator, cf):
check = _checkpoint(cf.output_path)
logs = _logs_manager(cf.output_path, gan)
train_datagen = GeneratorEnqueuer(
data_generator, use_multiprocessing=cf.use_multiprocessing)
train_datagen.start(cf.workers, cf.max_queue_size)
train_generator = train_datagen.get()
for epoch in range(cf.epochs):
print('Epoch %d/%d' % (epoch + 1, cf.epochs))
progbar = generic_utils.Progbar(data_generator.nb_steps)
# Training loop
for step in range(data_generator.nb_steps):
batch_real = next(train_generator)
batch_gray = np.expand_dims(batch_real[:, :, :, 0], -1)
batch_chroma = batch_real[:, :, :, 1:]
batch_fake = gan.combined.predict(batch_gray)[-1]
dis_real_X = np.concatenate((batch_gray, batch_chroma), axis=-1)
dis_fake_X = np.concatenate((batch_gray, batch_fake), axis=-1)
d_loss = []
for d, real_Y, fake_Y in zip(gan.discriminator, gan.dis_real_Y,
gan.dis_fake_Y):
d_real = d.train_on_batch(dis_real_X, real_Y)
d_fake = d.train_on_batch(dis_fake_X, fake_Y)
d_loss.append(0.5 * np.add(d_real, d_fake))
g_loss = gan.combined.train_on_batch(
batch_gray, gan.gen_real_Y + [batch_chroma])
if step % (data_generator.nb_steps //
cf.plots_per_epoch) == 0 and step != 0:
logs.save_plots(epoch, step,
data_generator.decoder(dis_fake_X),
data_generator.decoder(batch_real))
if step % (data_generator.nb_steps //
cf.weights_per_epoch) == 0 and step != 0 and epoch != 0:
check.save_weights(gan, epoch, step)
d_names = ['d_loss_%d' % i for i in range(cf.d_scales)]
g_names = ['g_cgan_%d' % i for i in range(cf.d_scales)]
logs.update(
names=['d_loss'] + d_names + ['g_loss'] + g_names + ['g_l1'],
values=[sum(d_loss)] + d_loss + [g_loss[0]] + g_loss[1:-1] +
[g_loss[-1]],
progbar=progbar,
display_step=cf.display_step)
train_datagen.stop()
def test_generator_enqueuer_processes():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_pcs(
DummySequence([3, 10, 10, 3])), use_multiprocessing=True)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for i in range(100):
acc.append(int(next(gen_output)[0, 0, 0, 0]))
assert acc != list(range(100)), ('Order was keep in GeneratorEnqueuer '
'with processes')
enqueuer.stop()
class ImageGenerator(object):
def __init__(self,
rgb_generator,
input_processing_function=lambda x: x,
label_processing_function=lambda x: None,
use_multiprocessing=False,
wait_time=0.01,
workers=4,
max_queue_size=10):
class BatchGenerator(object):
def __iter__(self):
return self
def __next__(self):
rgb_images = next(rgb_generator)
inputs, labels = [], []
if type(rgb_images) is tuple:
rgb_images, labels = rgb_images
inputs = np.array([
input_processing_function(rgb_image)
for rgb_image in rgb_images
])
if len(labels) == 0:
labels = [
label_processing_function(rgb_image)
for rgb_image in rgb_images
]
labels = [item for item in labels if item is not None]
if len(labels) == 0: return np.array(inputs)
else: return np.array(inputs), np.array(labels)
# Python2 compatibility
next = __next__
self.generator = GeneratorEnqueuer(
generator=BatchGenerator(),
use_multiprocessing=use_multiprocessing,
wait_time=wait_time)
self.generator.start(workers=workers, max_queue_size=max_queue_size)
def __iter__(self):
return self
def __next__(self):
return next(self.generator.get())
# Python2 compatibility
next = __next__
def test_generator_enqueuer_threads():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_threads(
TestSequence([3, 200, 200, 3])), use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for i in range(100):
acc.append(int(next(gen_output)[0, 0, 0, 0]))
"""
Not comparing the order since it is not guarantee.
It may get ordered, but not a lot, one thread can take the GIL before he was supposed to.
"""
assert len(set(acc) - set(range(100))) == 0, "Output is not the same"
enqueuer.stop()
def test_generator_enqueuer_threads():
enqueuer = GeneratorEnqueuer(create_generator_from_sequence_threads(
DummySequence([3, 200, 200, 3])),
use_multiprocessing=False)
enqueuer.start(3, 10)
gen_output = enqueuer.get()
acc = []
for i in range(100):
acc.append(int(next(gen_output)[0, 0, 0, 0]))
"""
Not comparing the order since it is not guaranteed.
It may get ordered, but not a lot, one thread can take the GIL before he was supposed to.
"""
assert len(set(acc) - set(range(100))) == 0, "Output is not the same"
enqueuer.stop()
def auc_eval(self):
if isinstance(self.validation_data, Generator):
assert self.validation_steps is not None, \
'If validation data is a generator, validation steps must be provided'
y_pred = []
y_true = []
enqueuer = GeneratorEnqueuer(self.validation_data,
use_multiprocessing=False,
wait_time=.01)
enqueuer.start(workers=1, max_queue_size=10)
output_generator = enqueuer.get()
for _ in range(self.validation_steps):
generator_output = next(output_generator)
if not hasattr(generator_output, '__len__'):
raise ValueError('Output of generator should be a tuple '
'(x, y, sample_weight) '
'or (x, y). Found: ' +
str(generator_output))
if len(generator_output) == 2:
x, y = generator_output
elif len(generator_output) == 3:
x, y, _ = generator_output
else:
raise ValueError('Output of generator should be a tuple '
'(x, y, sample_weight) '
'or (x, y). Found: ' +
str(generator_output))
outs = self.model.predict_on_batch(x)
y_pred += outs.tolist()
y_true += y.tolist()
enqueuer.stop()
else:
y_pred = self.model.predict(self.validation_data[0])
y_true = self.validation_data[1].astype(np.bool)
roc_auc = roc_auc_score(y_true=y_true, y_score=y_pred)
self.auc.append(roc_auc)
print('AUC Score is %s' % self.auc[-1])
class SegDataIterator(Iterator):
"""
支持voc和coco图像分割
"""
def __init__(self,
dataset,
image_size=None,
batch_size=32,
shuffle=False,
seed=None,
use_multiprocessing=True,
workers=4,
max_queue_size=10,
deprocess_X=None,
deprocess_Y=None):
super(SegDataIterator, self).__init__(len(dataset), batch_size,
shuffle, seed)
self.dataset = dataset
self.image_size = image_size[1], image_size[0]
self.enqueuer = GeneratorEnqueuer(
self, use_multiprocessing=use_multiprocessing)
self.workers = workers
self.max_queue_size = max_queue_size
self.is_start = False
self.deprocess_X = deprocess_X
self.deprocess_Y = deprocess_Y
def _get_batches_of_transformed_samples(self, index_array):
batch_x = np.zeros(
(len(index_array), self.image_size[1], self.image_size[0], 3),
np.float32)
batch_y = np.zeros((len(index_array), self.image_size[1],
self.image_size[0], self.dataset.num_classes),
np.float32)
for i, j in enumerate(index_array):
x = self.dataset.load_image(j, self.image_size)
y = self.dataset.load_segmentation(j, self.image_size)
batch_x[i] = x
batch_y[i] = y
if self.deprocess_X:
batch_x = self.deprocess_X(batch_x)
if self.deprocess_Y:
batch_y = self.deprocess_Y(batch_y)
return batch_x, batch_y
def next(self):
"""For python 2.x.
# Returns
The next batch.
"""
# Keeps under lock only the mechanism which advances
# the indexing of each batch.
with self.lock:
index_array = next(self.index_generator)
# The transformation of images is not under thread lock
# so it can be done in parallel
return self._get_batches_of_transformed_samples(index_array)
def data_gen(self):
if not self.is_start:
self.enqueuer.start(workers=self.workers,
max_queue_size=self.max_queue_size)
self.output_generator = self.enqueuer.get()
self.is_start = True
batch_x, batch_y = next(self.output_generator)
return batch_x, batch_y
class ModelDiagnoser(Callback):
def __init__(self, data_generator, m_batch_size, num_samples, output_dir, normalization_mean):
super().__init__()
self.epoch_index = 0
self.data_generator = data_generator
self.batch_size = m_batch_size
self.num_samples = num_samples
self.tensorboard_writer = TensorBoardWriter(output_dir)
self.normalization_mean = normalization_mean
is_sequence = isinstance(self.data_generator, Sequence)
if is_sequence:
self.enqueuer = OrderedEnqueuer(self.data_generator,
use_multiprocessing=True,
shuffle=False)
else:
self.enqueuer = GeneratorEnqueuer(self.data_generator,
use_multiprocessing=False, # todo: how to 'True' ?
wait_time=0.01)
# todo: integrate the Sequence generator properly
# import multiprocessing
# self.enqueuer.start(workers=multiprocessing.cpu_count(), max_queue_size=4)
self.enqueuer.start(workers=1, max_queue_size=4)
def on_epoch_end(self, epoch, logs=None):
output_generator = self.enqueuer.get()
steps_done = 0
total_steps = int(np.ceil(np.divide(self.num_samples, self.batch_size)))
sample_index = 0
while steps_done < total_steps:
generator_output = next(output_generator)
x, y = generator_output[:2]
x = next(iter(x.values()))
y = next(iter(y.values()))
y_pred = self.model.predict(x)
self.epoch_index += 1
for i in range(0, len(y_pred)):
n = steps_done * self.batch_size + i
if n >= self.num_samples:
return
# rearranging images for visualization
img_x = self.__reformat_img(x, i)
img_y = self.__reformat_img(y, i)
img_p = self.__reformat_img(y_pred, i)
self.tensorboard_writer.save_image("Epoch-{}/{}/x"
.format(self.epoch_index, sample_index), img_x)
self.tensorboard_writer.save_image("Epoch-{}/{}/y"
.format(self.epoch_index, sample_index), img_y)
self.tensorboard_writer.save_image("Epoch-{}/{}/y_pred"
.format(self.epoch_index, sample_index), img_p)
sample_index += 1
steps_done += 1
def __reformat_img(self, img_np_array, img_index):
img = np.squeeze(img_np_array[img_index, :, :, :])
img = 255. * (img + self.normalization_mean) # mean is the training images normalization mean
img = img[:, :, [2, 1, 0]] # reordering of channels
return img
def on_train_end(self, logs=None):
self.enqueuer.stop()
self.tensorboard_writer.close()
3000) #check point: every 10 epoch
step_lr_drop = 5
disc_losses = []
recont_losses = []
gen_losses = []
pre_loss = 9999
lr_current = lr_schedule[epoch]
real_ratio = 1.0
feed_iter = datagenerator.generator()
K.set_value(discriminator.optimizer.lr, lr_current)
K.set_value(dcgan.optimizer.lr, lr_current)
fed = GeneratorEnqueuer(feed_iter, use_multiprocessing=True, wait_time=5)
fed.start(workers=6, max_queue_size=200)
iter_ = fed.get()
zero_target = np.zeros((batch_size))
for X_src, X_tgt, disc_tgt, prob_gt in iter_:
discriminator.trainable = True
X_disc, y_disc = get_disc_batch(X_src,
X_tgt,
generator_train,
0,
label_smoothing=True,
label_flipping=0.2)
disc_loss = discriminator.train_on_batch(X_disc, y_disc)
X_disc, y_disc = get_disc_batch(X_src,
X_tgt,
generator_train,
def test_gan(gan, data_generator, cf):
results_path = os.path.join(cf.output_path, 'results')
if not os.path.exists(results_path): os.makedirs(results_path)
l1_cum = 0
psnr_cum = 0
perc_cum = 0
a_cum_dist = np.zeros(219)
b_cum_dist = np.zeros(219)
perception = perceptual_model()
gan.load_weights(cf.weights_path)
test_datagen = GeneratorEnqueuer(
data_generator, use_multiprocessing=cf.use_multiprocessing)
test_datagen.start(cf.workers, cf.max_queue_size)
test_generator = test_datagen.get()
for _ in tqdm.tqdm(range(data_generator.samples)):
batch_real, _ = next(test_generator)
batch_gray = np.expand_dims(batch_real[..., 0], -1)
batch_pred = gan.predict(batch_gray)
batch_fake = np.concatenate((batch_gray, batch_pred), axis=-1)
# compute color distribution
real_rgb = data_generator.decoder(batch_real)
fake_rgb = data_generator.decoder(batch_fake)
real_lab, fake_lab = [], []
for fake, real in zip(fake_rgb, real_rgb):
fake = rgb2lab(fake)
fake_lab.append(fake)
real_lab.append(rgb2lab(real))
a_cum_dist += np.histogram(fake[..., 1],
bins=np.arange(-110., 110.))[0]
b_cum_dist += np.histogram(fake[..., 2],
bins=np.arange(-110., 110.))[0]
fake_lab = np.array(fake_lab)
real_lab = np.array(real_lab)
# compute l1
l1_cum += np.abs(np.subtract(fake_lab[..., 1:], real_lab[...,
1:])).mean()
# compute psnr
mse = np.square(np.subtract(fake_lab[..., 1:], real_lab[...,
1:])).mean()
psnr_cum += np.sum(20 * np.log10(255. / np.sqrt(mse)))
# compute perceptual loss
real_rgb = resize(real_rgb)
fake_rgb = resize(fake_rgb)
perc_cum += sum(
perception.predict(
[preprocess_input(i) for i in [real_rgb, fake_rgb]]))
l1 = l1_cum / data_generator.samples
psnr = psnr_cum / data_generator.samples
perc_loss = perc_cum / data_generator.samples
with open(results_path + '/logs.txt', 'a') as f:
f.write('experiment, l1, psnr, perc loss\n')
f.write(cf.experiment_name + ', %f, %f, %f\n' % (l1, psnr, perc_loss))
a_dist = np.log(a_cum_dist / data_generator.samples)
b_dist = np.log(b_cum_dist / data_generator.samples)
a_dist[a_dist == -np.inf] = np.inf
a_dist[a_dist == np.inf] = a_dist.min()
b_dist[b_dist == -np.inf] = np.inf
b_dist[b_dist == np.inf] = b_dist.min()
np.save(os.path.join(results_path, 'a_dist.npy'), a_dist)
np.save(os.path.join(results_path, 'b_dist.npy'), b_dist)
def validate(config, model, val_client, validation_steps, metrics_id, epoch):
val_di = val_client.gen()
from keras.utils import GeneratorEnqueuer
val_thre = GeneratorEnqueuer(val_di)
val_thre.start()
model_metrics = []
inhouse_metrics = []
for i in range(validation_steps):
X, GT = next(val_thre.get())
Y = model.predict(X)
model_losses = [(np.sum((gt - y)**2) / gt.shape[0] / 2)
for gt, y in zip(GT, Y)]
mm = sum(model_losses)
if config.paf_layers > 0 and config.heat_layers > 0:
GTL6 = np.concatenate([GT[-2], GT[-1]], axis=3)
YL6 = np.concatenate([Y[-2], Y[-1]], axis=3)
mm6l1 = model_losses[-2]
mm6l2 = model_losses[-1]
elif config.paf_layers == 0 and config.heat_layers > 0:
GTL6 = GT[-1]
YL6 = Y[-1]
mm6l1 = None
mm6l2 = model_losses[-1]
else:
assert False, "Wtf or not implemented"
m = calc_batch_metrics(i, GTL6, YL6,
range(config.heat_start, config.bkg_start))
inhouse_metrics += [m]
model_metrics += [(i, mm, mm6l1, mm6l2, m["MAE"].sum() / GTL6.shape[0],
m["RMSE"].sum() / GTL6.shape[0], m["DIST"].mean())]
print(
"Validating[BATCH: %d] LOSS: %0.4f, S6L1: %0.4f, S6L2: %0.4f, MAE: %0.4f, RMSE: %0.4f, DIST: %0.2f"
% model_metrics[-1])
inhouse_metrics = pd.concat(inhouse_metrics)
inhouse_metrics['epoch'] = epoch
inhouse_metrics.to_csv("logs/val_scores.%s.%04d.txt" % (metrics_id, epoch),
sep="\t")
model_metrics = pd.DataFrame(model_metrics,
columns=("batch", "loss", "stage6l1",
"stage6l2", "mae", "rmse", "dist"))
model_metrics['epoch'] = epoch
del model_metrics['batch']
model_metrics = model_metrics.groupby('epoch').mean()
with open('%s.val.tsv' % metrics_id, 'a') as f:
model_metrics.to_csv(f,
header=(epoch == 1),
sep="\t",
float_format='%.4f')
val_thre.stop()
def validate(config, model, multi_model, val_client, validation_steps,
metrics_id, epoch):
# 得到的X是包含image, confidence mask, paf mask的list,得到的Y是包含6个stage一共12个groundtruth的heapmap
# 网络一共有三个输入（对于训练时的评估，指标为了反映出训练的效果在测试时网络模型就不用考虑对feature map的输出进行mask了，
# 可以对所有区域预测）,即原始image，以及在训练模型过程中评估时，去除没有标记区域的confidence和paf的mask1, mask2
val_di = val_client.gen()
val_thre = GeneratorEnqueuer(
val_di
) # The provided generator can be finite in which case the class will throw
# a `StopIteration` exception. 但是这里实现的gen貌似不存在这种问题。不过这个函数提供了multiprocess的封装
val_thre.start()
model_metrics = []
inhouse_metrics = []
t0 = time()
for i in range(validation_steps
): # 分成很多个batch进行预测估计的，为了减少validation耗时，在计算validation部分数据
# validation_steps ＝ val_samples//batch_size 为了防止内存OOM，所以要分batch预测
# if random.randint(0, 9) < 5: # 只计算20%的数据
# continue
X, GT = next(val_thre.get())
Y = multi_model.predict(X)
model_losses = [(np.sum((gt - y)**2) / gt.shape[0] / 2)
for gt, y in zip(GT, Y)]
# 与模型定义时的loss保持一致，除以2的好处是，平方项的微分会出现2，抵消，可以减少乘法操作
mm = sum(model_losses)
if config.paf_layers > 0 and config.heat_layers > 0:
GTL6 = np.concatenate([GT[-2], GT[-1]], axis=3)
YL6 = np.concatenate([Y[-2], Y[-1]], axis=3)
mm6l1 = model_losses[-2] # NOTICE! 计算的是模型最后一个阶段的预测和groundtruth的距离
mm6l2 = model_losses[-1]
elif config.paf_layers == 0 and config.heat_layers > 0:
GTL6 = GT[-1]
YL6 = Y[-1]
mm6l1 = None
mm6l2 = model_losses[-1]
else:
assert False, "Wtf or not implemented"
m = calc_batch_metrics(i, GTL6, YL6,
range(config.heat_start, config.bkg_start))
inhouse_metrics += [m]
model_metrics += [(i, mm, mm6l1, mm6l2, m["MAE"].sum() / GTL6.shape[0],
m["RMSE"].sum() / GTL6.shape[0], m["DIST"].mean())]
# 以epoch为key，group之后取平均值
print(
"Validating[BATCH: %d] LOSS: %0.4f, S6L1: %0.4f, S6L2: %0.4f, MAE: %0.4f, RMSE: %0.4f, DIST: %0.2f"
% model_metrics[-1])
t1 = time()
print('The CNN prediction time during validation is : ', t1 - t0)
# inhouse_metrics = pd.concat(inhouse_metrics)
# inhouse_metrics['epoch'] = epoch
# inhouse_metrics.to_csv("logs/val_scores.%s.%04d.csv" % (metrics_id, epoch)) # , sep="\t" 默认的不是\t，而是','
# # 保存的是每个层的细节
#
# model_metrics = pd.DataFrame(model_metrics, columns=("batch","loss","stage6l1","stage6l2","mae","rmse","dist") )
# model_metrics['epoch'] = epoch
# del model_metrics['batch']
# model_metrics = model_metrics.groupby('epoch').mean()
# with open('%s.val.tsv' % metrics_id, 'a') as f:
# model_metrics.to_csv(f, header=(epoch==1), float_format='%.4f') # sep="\t",
#
# print(inhouse_metrics[["layer", "epoch", "MAE", "RMSE", "DIST"]].groupby(["layer", "epoch"]).mean())
val_thre.stop()
