import matplotlib.pyplot as plt
import numpy as np
from spectractor.extractor.psf import Moffat, Order0
from spectractor.extractor.images import Image, find_target
im = Image('tests/data/reduc_20170605_028.fits', target_label="PNG321.0+3.9")
im.plot_image()
guess = [820, 580]
parameters.VERBOSE = True
parameters.DEBUG = True
x0, y0 = find_target(im, guess)
p = [1,40,50,1,1e20]
psf = Order0(target=im.target, p=p)
yy, xx = np.mgrid[:80, :100]
out = psf.evaluate(pixels=np.array([xx, yy]), values=p)
fig = plt.figure(figsize=(5,5))
plt.imshow(out, origin="lower")
plt.xlabel("X [pixels]")
plt.ylabel("Y [pixels]")
plt.grid()
plt.show()