def test_explore():
I = Image(make_data_int16())
J = Image(make_data_int16())
IM = IconicMatcher(I.array, J.array, I.toworld, J.toworld)
T = np.eye(4)
T[0:3,3] = np.random.rand(3)
simi, params = IM.explore(ux=[-1,0,1],uy=[-1,0,1])
matcher = IconicMatcher(I.get_data(), J.get_data(),
I.get_affine(), J.get_affine()) ## I: source, J: target
# Params
size = 5
nsimu = 1
depth = 10
import pylab as pl
# Simulations
for i in range(nsimu):
# Select random block
x0 = np.random.randint(I.array.shape[0]-size)
y0 = np.random.randint(I.array.shape[1]-size)
z0 = np.random.randint(I.array.shape[2]-size)
matcher.set(corner=[x0,y0,z0], size=[size,size,size])
# Explore neighborhood
tx = I.voxsize[0] * (np.arange(2*depth + 1)-depth)
s, p = matcher.explore(tx=tx)
# Display
pl.plot(p[:,0],(1-s)**(-.5*size**3))
pl.show()
