def main():
parser = argparse.ArgumentParser()
parser.add_argument('input_all_iterations_Xy_path')
parser.add_argument('output_dir')
parser.add_argument('--visualise-progress',
action='store_true',
default=False)
parser.add_argument('--ftol', type=float, default=1e-8)
parser.add_argument('--xtol', type=float, default=1e-8)
parser.add_argument('--maxfev', type=int, default=0)
parser.add_argument('--epsilon', type=float, default=1e-6)
args = parser.parse_args()
ensure_output_path = partial(vis.ensure_path, args.output_dir)
all_iterations_Xy, orig_args = np.load(args.input_all_iterations_Xy_path)
print '<-', orig_args['input_path']
I = plt.imread(orig_args['input_path']).astype(np.float64)[..., :3]
if orig_args['base'] == 'white':
J0 = np.ones_like(I)
elif orig_args['base'] == 'black':
J0 = np.zeros_like(I)
else:
head, tail = os.path.split(orig_args['base'])
root, ext = os.path.splitext(tail)
if ext == '.dat':
J0 = np.load(orig_args['base'])
else:
J0 = plt.imread(orig_args['base']).astype(np.float64)[..., :3]
Xs0, ys0 = zip(*map(itemgetter(-1), all_iterations_Xy))
print 'Solving with `fit_and_colour_shards` ...'
np.seterr(over='ignore')
t0 = time()
(Xs, ys, all_Xs_ys), (exit_code, E0, E1, J, J1) = fit_and_colour_shards(
I, J0, orig_args['alpha'],
Xs0, ys0,
k=orig_args['k'],
epsilon=args.epsilon,
ftol=args.ftol,
xtol=args.xtol,
maxfev=args.maxfev,
return_info=True,
verbose=True)
t1 = time()
np.seterr(over='warn')
print 'E0:', E0
print 'E1:', E1
print 'Exit code: %d' % exit_code
print 'Time taken: %.3fs' % (t1 - t0)
output_path = ensure_output_path('all_Xs_ys.dat')
print '->', output_path
dump(output_path, (all_Xs_ys, args.__dict__), raise_on_failure=False)
output_path = ensure_output_path('J.dat')
print '->', output_path
dump(output_path, (J, args.__dict__), raise_on_failure=False)
output_path = ensure_output_path('J1.dat')
print '->', output_path
dump(output_path, (J1, args.__dict__), raise_on_failure=False)
if args.visualise_progress:
output_path = ensure_output_path('J.png')
print '->', output_path
f, ax = vis.make_image_figure(J)
vis.save_image_figure(output_path, f, J.shape)
output_path = ensure_output_path('J1.png')
print '->', output_path
f, ax = vis.make_image_figure(J1)
vis.save_image_figure(output_path, f, J1.shape)
def main():
parser = argparse.ArgumentParser()
parser.add_argument('input_path')
parser.add_argument('output_dir')
parser.add_argument('base', help='{white, black, BASE}')
parser.add_argument('k', type=float)
parser.add_argument('alpha', type=float)
parser.add_argument('shards', type=int)
parser.add_argument('--n', type=int, default=3)
parser.add_argument('--xtol', type=float, default=1e-4)
parser.add_argument('--epsilon', type=float, default=1e-6)
parser.add_argument('--visualise-progress',
action='store_true',
default=False)
parser.add_argument('--continue-from-previous',
action='store_true',
default=False)
parser.add_argument('--num-restarts-per-shard',
type=int,
default=1)
parser.add_argument('--dont-solve-joint',
dest='solve_joint',
default=True,
action='store_false')
args = parser.parse_args()
ensure_output_path = partial(vis.ensure_path, args.output_dir)
print '<-', args.input_path
I = plt.imread(args.input_path).astype(np.float64)[..., :3]
if args.base == 'white':
J = np.ones_like(I)
elif args.base == 'black':
J = np.zeros_like(I)
else:
head, tail = os.path.split(args.base)
root, ext = os.path.splitext(tail)
if ext == '.dat':
J = np.load(args.base)
else:
J = plt.imread(args.base).astype(np.float64)[..., :3]
sr = ShardReconstructor(I, args.alpha, n=args.n, k=args.k)
if args.continue_from_previous:
for n in xrange(args.shards - 1, -1, -1):
p = os.path.join(args.output_dir, str(n), 'J.dat')
try:
J = np.load(p)
except IOError:
continue
break
n += 1
else:
n = 0
print 'n: [%d, %d)' % (n, args.shards)
t1 = time()
all_iterations_Xy = []
for n in xrange(n, args.shards):
E0 = sr.reconstruction_energy(J)
max_E_diff = 0.0
for i in xrange(args.num_restarts_per_shard):
X, y, all_Xy = sr.candidate_shard(J, args.solve_joint,
verbose=True,
epsilon=args.epsilon,
xtol=args.xtol)
J1 = sr.add_shard_to_reconstruction(J, X, y)
E1 = sr.reconstruction_energy(J1)
E_diff = E0 - E1
if E_diff > max_E_diff:
max_X_y = (X, y, all_Xy)
max_E_diff = E_diff
min_E1 = E1
if max_E_diff == 0.0:
# no shard produced a decrease in energy.
print 'shard %d ignored' % n
continue
X, y, all_Xy = max_X_y
all_iterations_Xy.append(all_Xy)
print 'E:', E0
print 'E(next):', min_E1
J1s = map(lambda t: sr.add_shard_to_reconstruction(J, t[0], t[1]),
all_Xy)
J = J1s[-1]
output_dir = ensure_output_path(n, is_dir=True)
if args.visualise_progress:
vis.make_visualisations_inplace(map(itemgetter(0), all_Xy),
J1s,
output_dir, verbose=True)
vis.make_residual_image(I, J1s[-1],
output_dir, verbose=True)
output_path = ensure_output_path(n, 'all_Xy.dat')
print '->', output_path
dump(output_path, (all_Xy, args.__dict__), raise_on_failure=False)
output_path = ensure_output_path(n, 'all_iterations_Xy.dat')
print '->', output_path
dump(output_path, (all_iterations_Xy, args.__dict__),
raise_on_failure=False)
output_path = ensure_output_path(n, 'J.dat')
print '->', output_path
dump(output_path, J, raise_on_failure=False)
print '%d shards in %.3fs' % (n + 1, time() - t1)