def evaluate(self, position):
"""Evaluate model at a given position.
Parameters
----------
position : `~astropy.coordinates.SkyCoord`
Position on the sky.
"""
glon, glat = position.galactic.l, position.galactic.b
width = self.width(glon)
amplitude = self.peak_brightness(glon)
mean = self.peak_latitude(glon)
return Gaussian1D.evaluate(glat, amplitude=amplitude, mean=mean, stddev=width)
def evaluate(self, position):
"""Evaluate model at a given position.
Parameters
----------
position : `~astropy.coordinates.SkyCoord`
Position on the sky.
"""
glon, glat = position.galactic.l, position.galactic.b
width = self.width(glon)
amplitude = self.peak_brightness(glon)
mean = self.peak_latitude(glon)
return Gaussian1D.evaluate(glat, amplitude=amplitude, mean=mean, stddev=width)
def _evaluate_y(self, y, pars):
"""Evaluate Gaussian model at a given ``y`` position.
"""
return Gaussian1D.evaluate(y, **pars)
def _evaluate_y(self, y, pars):
"""Evaluate Gaussian model at a given ``y`` position.
"""
return Gaussian1D.evaluate(y, **pars)
def evaluate(self, x, y, sigma, c0, c1, c2, c3):
self.polynomial.coef = np.array([c0, c1, c2, c3])
poly_eval = self.polynomial(x)
poly_eval = (poly_eval + (self.window[0] + 1)) * self.window_delta / 2.
return Gaussian1D.evaluate(y, 1., poly_eval, sigma)
