######################################################################################## #------------------------------- Figure ----------------------------- ######################################################################################## colors=["m","b","c","r"] markers=["^","o",".","v"] ncolor=0 nmarker=0 nc = 450*3 nr = 450*3 surface = np.zeros((nr,nc)) classique.lecturesurface(filename,surface) # for j in range(0,len(ordonnee)): # surface[int(round(abscisse_1[j]))+(int(1./(2*lc))-1)/2,int(round(abscisse_2[j]))+(int(1./(2*lc))-1)/2]=ordonnee[j] # mask[int(round(abscisse_1[j]))+(int(1./(2*lc))-1)/2,int(round(abscisse_2[j]))+(int(1./(2*lc))-1)/2]=0 #surface[int(round(abscisse_1[j])),int(round(abscisse_2[j]))]=ordonnee[j] #mask[int(round(abscisse_1[j])),int(round(abscisse_2[j]))]=0 #print(abscisse_1[j],abscisse_2[j],int(round(abscisse_1[j])),int(round(abscisse_2[j])),ordonnee[j]) #print(mask) #print(surface) #if maillage==3: # xextent=[lc*min(abscisse_1)-lc*(int(1./(2*lc))-1)/2,lc*max(abscisse_1)+lc*(int(1./(2*lc))-1)/2] # yextent=[lc*min(abscisse_2)-lc*(int(1./(2*lc))-1)/2,lc*max(abscisse_2)+lc*(int(1./(2*lc))-1)/2] #else:
