def test_gauss_center():
"""Check the gaussian center fitting."""
# make a 2d dataset with a gaussian at the center
from astropy.convolution import Gaussian2DKernel
gaussian_2D_kernel = Gaussian2DKernel(10)
plots = Imexamine()
plots.set_data(gaussian_2D_kernel.array)
a, xx, yy, xs, ys = plots.gauss_center(37, 37)
amp = 0.0015915494309189533
xc = 40.0
yc = 40.0
xsig = 10.0
ysig = 10.0
assert_allclose(amp, a, 1e-6)
assert_allclose(xc, xx, 1e-4)
assert_allclose(yc, yy, 1e-4)
assert_allclose(xsig, xs, 0.01)
assert_allclose(ysig, ys, 0.01)
def test_gauss_center():
"""Check the gaussian center fitting."""
# make a 2d dataset with a gaussian at the center
from astropy.convolution import Gaussian2DKernel
gaussian_2D_kernel = Gaussian2DKernel(10)
plots = Imexamine()
plots.set_data(gaussian_2D_kernel.array)
a, xx, yy, xs, ys = plots.gauss_center(37, 37)
amp = 0.0015915494309189533
xc = 40.0
yc = 40.0
xsig = 10.0
ysig = 10.0
assert_allclose(amp, a, 1e-6)
assert_allclose(xc, xx, 1e-4)
assert_allclose(yc, yy, 1e-4)
assert_allclose(xsig, xs, 0.01)
assert_allclose(ysig, ys, 0.01)
def test_gauss_center():
"""Check the gaussian center fitting."""
from astropy.convolution import Gaussian2DKernel
# This creates a 2D normalized gaussian kernal with
# a set amplitude. Guess off-center
amp = 0.0015915494309189533
size = 81.0
sigma = 10.0
gaussian_2D_kernel = Gaussian2DKernel(sigma, x_size=size, y_size=size)
plots = Imexamine()
plots.set_data(gaussian_2D_kernel.array)
a, xx, yy, xs, ys = plots.gauss_center(37, 37)
assert_allclose(amp, a, 1e-6)
assert_allclose(size // 2, xx, 1e-6)
assert_allclose(size // 2, yy, 1e-6)
assert_allclose(sigma, xs, 0.01)
assert_allclose(sigma, ys, 0.01)
def test_gauss_center():
"""Check the gaussian center fitting."""
from astropy.convolution import Gaussian2DKernel
# This creates a 2D normalized gaussian kernal with
# a set amplitude. Guess off-center
amp = 0.0015915494309189533
size = 81.0
sigma = 10.0
gaussian_2D_kernel = Gaussian2DKernel(sigma, x_size=size, y_size=size)
plots = Imexamine()
plots.set_data(gaussian_2D_kernel.array)
a, xx, yy, xs, ys = plots.gauss_center(37, 37)
assert_allclose(amp, a, 1e-6)
assert_allclose(size // 2, xx, 1e-6)
assert_allclose(size // 2, yy, 1e-6)
assert_allclose(sigma, xs, 0.01)
assert_allclose(sigma, ys, 0.01)