def tester_image_simple():
f,g=lib.get_images()
lamb=[0,10,50,70,100]
mu=[0,10,50,70,100]
stepini=[0.001]
nitermax=[100000]
for m in mu:
for s in stepini:
for n in nitermax:
for l in lamb:
produire_image("simple",f,g,l,m,n,s)
def tester_image_simple():
f, g = lib.get_images()
lamb = [0, 10, 50, 70, 100]
mu = [0, 10, 50, 70, 100]
stepini = [0.001]
nitermax = [100000]
for m in mu:
for s in stepini:
for n in nitermax:
for l in lamb:
produire_image("simple", f, g, l, m, n, s)
def produire_image(nom, f, g, lamb, mu, nitermax, stepini):
ux, uy, CF, step = lib.RecalageDG_TP(f, g, lamb, mu, nitermax, stepini)
#ux,uy,CF,step=lib.RecalageGN_TP(f,g,lamb,mu,nitermax,stepini,100)
plt.rcParams['text.usetex'] = 'true'
plt.rcParams['text.latex.unicode'] = 'true'
plt.rcParams['font.family'] = 'serif'
plt.rcParams['font.size'] = 9
fig, ax = plt.subplots(2, 3)
ax[0, 0].imshow(f, origin='lower')
ax[0, 0].set_title('original function $f$')
ax[0, 1].imshow(g, origin='lower')
ax[0, 1].set_title('target function $g$')
ax[1, 0].quiver(uy, ux, color='b')
ax[1, 0].set_title('displacement field $u$')
ax[1, 1].imshow(lib.interpol(f, ux, uy), origin='lower')
ax[1, 1].set_title('final function $f \circ (I+u)$')
ax[0, 2].plot(CF)
ax[0, 2].set_title('objective history')
ax[1, 2].plot(np.log(step))
ax[1, 2].set_title('step history (log scale)')
filename = nom + "-lam_" + str(lamb) + "-mu_" + str(
mu) + "-stepini_" + str(stepini) + "-iter_" + str(step.size) + ".png"
fig.savefig(filename, bbox_inches='tight')
def produire_image(nom,f,g,lamb,mu,nitermax,stepini):
ux,uy,CF,step=lib.RecalageDG_TP(f,g,lamb,mu,nitermax,stepini)
#ux,uy,CF,step=lib.RecalageGN_TP(f,g,lamb,mu,nitermax,stepini,100)
plt.rcParams['text.usetex'] = 'true'
plt.rcParams['text.latex.unicode'] = 'true'
plt.rcParams['font.family'] = 'serif'
plt.rcParams['font.size'] = 9
fig, ax = plt.subplots(2,3)
ax[0,0].imshow(f, origin='lower')
ax[0,0].set_title('original function $f$')
ax[0,1].imshow(g, origin='lower')
ax[0,1].set_title('target function $g$')
ax[1,0].quiver(uy,ux, color='b')
ax[1,0].set_title('displacement field $u$')
ax[1,1].imshow(lib.interpol(f,ux,uy), origin='lower')
ax[1,1].set_title('final function $f \circ (I+u)$')
ax[0,2].plot(CF)
ax[0,2].set_title('objective history')
ax[1,2].plot(np.log(step))
ax[1,2].set_title('step history (log scale)')
filename = nom+"-lam_"+str(lamb)+"-mu_"+str(mu)+"-stepini_"+str(stepini)+"-iter_"+str(step.size)+".png"
fig.savefig(filename, bbox_inches='tight')
#from joblib import Parallel, delayed
def tester_image_simple():
f,g=lib.get_images()
lamb=[0,10,50,70,100]
mu=[0,10,50,70,100]
stepini=[0.001]
nitermax=[100000]
for m in mu:
for s in stepini:
for n in nitermax:
for l in lamb:
produire_image("simple",f,g,l,m,n,s)
#Parallel(n_jobs=4)(delayed(produire_image)("simple",f,g,l,m,n,s) for l in lamb)
if __name__ == '__main__':
f,g=lib.get_images()
produire_image("simple",f,g,50,50,1000,0.001)
#tester_image_simple()
#parser = argparse.ArgumentParser(description='Lancer résolution fonction simple')
#parser.add_argument('--dossier',action='store_true',default=''\
# help="Pour la commande --dot, afficher les numéros des noeuds")
#f,g=lib.get_images()
#produire_image("simple",f,g,50,50,1000,0.001)
tester_image_simple()
import argparse
#from joblib import Parallel, delayed
def tester_image_simple():
f, g = lib.get_images()
lamb = [0, 10, 50, 70, 100]
mu = [0, 10, 50, 70, 100]
stepini = [0.001]
nitermax = [100000]
for m in mu:
for s in stepini:
for n in nitermax:
for l in lamb:
produire_image("simple", f, g, l, m, n, s)
#Parallel(n_jobs=4)(delayed(produire_image)("simple",f,g,l,m,n,s) for l in lamb)
if __name__ == '__main__':
f, g = lib.get_images()
produire_image("simple", f, g, 50, 50, 1000, 0.001)
#tester_image_simple()
#parser = argparse.ArgumentParser(description='Lancer résolution fonction simple')
#parser.add_argument('--dossier',action='store_true',default=''\
# help="Pour la commande --dot, afficher les numéros des noeuds")
#f,g=lib.get_images()
#produire_image("simple",f,g,50,50,1000,0.001)
tester_image_simple()
import argparse
