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