def run_msc(target):
print target.shape
source=makeball(target.shape[0]/2, target.shape[1]/2, 2, target.shape)
trans=[]
#trans.append(BinImScale(1))
trans.append(HShift(target.shape[0]/2))
trans.append(VShift(target.shape[1]/2))
ms=MapSeeker(trans, {'maxiter':200, 'kappa':0.2, 'rule':'kappa', 'checkdomain':False, 'silent':2})
print source.shape, target.shape
ms.run(source.astype(float32),target.astype(float32))
return ms
def run_msc(target):
lv=LocalizingVectorField(target.shape)
source=lv.local(array([target.shape[0]/2, target.shape[1]/2]), 2)
trans=[]
#trans.append(BinImScale(1))
trans.append(AShift(target.shape[0]/2, 0, target.shape))
trans.append(AShift(target.shape[1]/2, 1, target.shape))
ms=MapSeeker(trans, {'maxiter':200, 'kappa':0.2, 'rule':'kappa', 'silent':2})
target=lv.fromFull(target)
o=target.full()
ms.run(source,target)
return ms
