#plt.show(block=True)
plt.gcf().savefig('../data/barim2-1.png')
np.save(fixed,im1)
#plt.imshow(im2, interpolation='nearest')
#plt.show(block=True)
#plt.gcf().savefig('../data/barim2-2.png')
#np.save(moving,im2)
# imm transform
imm = np.load(moving)
[Nx,Ny] = np.shape(imm)
center = np.array([[Nx/2.,Ny/2.]])
def T(XY):
phi = 2*np.pi*0.09
R = np.array([[np.cos(phi),-np.sin(phi)],[np.sin(phi),np.cos(phi)]])
res = np.dot(R,XY.T-np.kron(np.ones((1,XY.shape[0])),center.T))+np.kron(np.ones((1,XY.shape[0])),center.T)
return res.T
immT = T
sim = imsim.get(pointsPerAxis, immname=moving, imfname=fixed, immT=immT, border=border, normalize=True, visualize=visualize, order=order, smoothscaleFactor=smoothing, SIGMAF=SIGMAF, h=h, splineS=splineS)
logging.info("initial point configuration: %s",sim['initial'])
(fstate,res) = match.match(sim,sim['SIGMA'],weights,initial=sim['initial'],gradTol=gradTol,order=order,scalegrad=scalegrad,maxIter=maxIter,visualize=visualize, visualizeIterations=visualizeIterations)
#print res
if True: # res.success:
match.genStateData(fstate,sim)
