def test_plot_tests(imsize, upsample_factor,doplot=False):
"""
Test that, when upsampled, z[::upsample_factor] = g
Tested for offsets in the *input gaussian function*
This doesn't imply you know what's going on, just that the behavior is consistent
ONLY tests for OUTSIZE=INSIZE*UPSAMPLE_FACTOR
"""
# always centered
x = np.arange(imsize)
g = gaussian(x - (imsize-1)/2.)
xz,z = zoom.zoom1d(g,upsample_factor,outsize=g.size*upsample_factor,return_xouts=True)
if doplot:
print x,xz
from pylab import *
figure(1); clf()
plot(xz,z)
plot(x,g,label='input')
plot(xz,gaussian(xz-(imsize-1)/2.))
figure(2); clf()
plot(z[::upsample_factor])
plot(g)
figure(3); clf()
plot(gaussian(xz-(imsize-1)/2.)-z)
assert ((gaussian(xz - (imsize-1)/2.)-z)**2).sum() < 4e-5*upsample_factor+0.12*imsize**-1.83
def test_plot_tests(imsize, upsample_factor, doplot=False):
"""
Test that, when upsampled, z[::upsample_factor] = g
Tested for offsets in the *input gaussian function*
This doesn't imply you know what's going on, just that the behavior is consistent
ONLY tests for OUTSIZE=INSIZE*UPSAMPLE_FACTOR
"""
# always centered
x = np.arange(imsize)
g = gaussian(x - (imsize - 1) / 2.)
xz, z = zoom.zoom1d(g,
upsample_factor,
outsize=g.size * upsample_factor,
return_xouts=True)
if doplot:
print x, xz
from pylab import *
figure(1)
clf()
plot(xz, z)
plot(x, g, label='input')
plot(xz, gaussian(xz - (imsize - 1) / 2.))
figure(2)
clf()
plot(z[::upsample_factor])
plot(g)
figure(3)
clf()
plot(gaussian(xz - (imsize - 1) / 2.) - z)
assert ((gaussian(xz - (imsize - 1) / 2.) - z)**
2).sum() < 4e-5 * upsample_factor + 0.12 * imsize**-1.83
def test_gaussian_upsample_oddoffsets(imsize, upsample_factor, offset):
"""
Test that, when upsampled, z[::upsample_factor] = g
Tested for offsets in the *input gaussian function*
This doesn't imply you know what's going on, just that the behavior is consistent
ONLY tests for OUTSIZE=INSIZE*UPSAMPLE_FACTOR
"""
x = np.arange(imsize)
g = gaussian(x-(imsize-1)/2.-offset)
xz,z = zoom.zoom1d(g,upsample_factor,outsize=g.size*upsample_factor,return_xouts=True)
assert ((gaussian(xz - (imsize-1)/2. -offset)-z)**2).sum() < 4e-5*upsample_factor+0.12*imsize**-1.83
def test_gaussian_upsample_oddoffsets(imsize, upsample_factor, offset):
"""
Test that, when upsampled, z[::upsample_factor] = g
Tested for offsets in the *input gaussian function*
This doesn't imply you know what's going on, just that the behavior is consistent
ONLY tests for OUTSIZE=INSIZE*UPSAMPLE_FACTOR
"""
x = np.arange(imsize)
g = gaussian(x - (imsize - 1) / 2. - offset)
xz, z = zoom.zoom1d(g,
upsample_factor,
outsize=g.size * upsample_factor,
return_xouts=True)
assert ((gaussian(xz - (imsize - 1) / 2. - offset) - z)**
2).sum() < 4e-5 * upsample_factor + 0.12 * imsize**-1.83
