def test_curve_of_growth():
"""Test the cog functionality."""
from astropy.convolution import Gaussian2DKernel
data = Gaussian2DKernel(1.5, x_size=25, y_size=25)
plots = Imexamine()
plots.set_data(data.array)
x, y = plots.curve_of_growth(12, 12, genplot=False)
rad = [1, 2, 3, 4, 5, 6, 7, 8]
flux = [0.04535423476987057,
0.34605795394960859,
0.70094774639729907,
0.89675381769455764,
0.97242001951216395,
0.99711765053819645,
0.99965838382174854,
0.99998044756724924]
assert_array_equal(rad, x)
assert_allclose(flux, y, 1e-6)
def test_curve_of_growth():
"""Test the cog functionality."""
from astropy.convolution import Gaussian2DKernel
data = Gaussian2DKernel(1.5, x_size=25, y_size=25)
plots = Imexamine()
plots.set_data(data.array)
x, y = plots.curve_of_growth(12, 12, genplot=False)
rad = [1, 2, 3, 4, 5, 6, 7, 8]
flux = [0.04535423476987057,
0.34605795394960859,
0.70094774639729907,
0.89675381769455764,
0.97242001951216395,
0.99711765053819645,
0.99965838382174854,
0.99998044756724924]
assert_array_equal(rad, x)
assert_allclose(flux, y, 1e-6)
def test_curve_of_growth():
"""Test the cog functionality."""
from astropy.convolution import Gaussian2DKernel
data = Gaussian2DKernel(1.5, x_size=25, y_size=25)
plots = Imexamine()
plots.set_data(data.array)
rad_out, flux_out = plots.curve_of_growth(12, 12, genplot=False)
rads = [1, 2, 3, 4, 5, 6, 7, 8]
flux = []
# Run the aperture phot on this to build up the expected fluxes
plots.aper_phot_pars['genplot'][0] = False
plots.aper_phot_pars['subsky'][0] = False
for rad in rads:
plots.aper_phot_pars['radius'][0] = rad
plots.aper_phot(12, 12)
flux.append(plots.total_flux)
assert_array_equal(rads, rad_out)
assert_allclose(flux, flux_out, 1e-6)
def test_curve_of_growth():
"""Test the curve of growth functionality."""
from astropy.convolution import Gaussian2DKernel
data = Gaussian2DKernel(1.5, x_size=25, y_size=25)
plots = Imexamine()
plots.set_data(data.array)
rad_out, flux_out = plots.curve_of_growth(12, 12, genplot=False)
rads = [1, 2, 3, 4, 5, 6, 7, 8]
flux = []
# Run the aperture phot on this to build up the expected fluxes
plots.aper_phot_pars['genplot'][0] = False
plots.aper_phot_pars['subsky'][0] = False
for rad in rads:
plots.aper_phot_pars['radius'][0] = rad
plots.aper_phot(12, 12)
flux.append(plots.total_flux)
assert_array_equal(rads, rad_out)
assert_allclose(flux, flux_out, 1e-6)
