[AstroPy] Gaussian with constant background

[Erik Tollerud]

Depending on how you think, you might also find pymodelfit
(http://packages.python.org/PyModelFit/ - for full disclosure, I'm the
author) useful - it has gaussian models for both 1D and 2D datasets
that you can just fit like:

gmod = pymodelfit.GaussianModel()
gmod.fitData(x, y)

And if you want a constant background as you're saying here, you could
just add in an offset:

gmod2 = pymodelfit.offset_model(pymodelfit.GaussianModel, 'bkgd')
gmod2.bkgd = whatever your constant background should be
gmod2.fitData(x, y, fixedpars='bkgd')


This is all scipy.optimize functions under the hood, so it should
perform pretty much the same - just rather easier to use.



On Wed, Jun 27, 2012 at 1:10 PM, Vivienne Baldassare
<vivienne.baldassare at gmail.com> wrote:
> Hello,
>
> I've been trying to figure out how to fit a gaussian with a constant
> (non-zero) background.  Aside from just subtracting off the background, I
> haven't been able to do it, and googling around was not particularly
> helpful.  I'd appreciate any tips!
>
> Best,
> Vivienne Baldassare
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at scipy.org
> http://mail.scipy.org/mailman/listinfo/astropy
>



-- 
Erik Tollerud