[AstroPy] Ellipse2D fit
Fabricio Ferrari
fabricio.ferrari at gmail.com
Wed Feb 15 08:07:53 EST 2017
yes of course, I did not paid attention. sorry :-D
I have been using Fitzgibbon, Pilu and Fischer solution for fitting
ellipses, which is states the problem such that it can be written as and
eigenvalue equation (it is not iterative). Please, take a look for example
at
http://nicky.vanforeest.com/misc/fitEllipse/fitEllipse.html
The original paper I think it is:
http://cseweb.ucsd.edu/~mdailey/Face-Coord/ellipse-specific-fitting.pdf
regards.
F
..-. ..-.
Fabricio Ferrari [www.ferrari.pro.br]
IMEF - FURG
Rio Grande, RS
Brasil
2017-02-15 10:47 GMT-02:00 Andrzej Szary <aszary at gmail.com>:
> Dear Fabricio,
>
> Note that the DATA and MODEL parameters are different. For a good fit they
> should be similar. You can get better results using different fitter (*fi
> = fitting.SimplexLSQFitter()*), but it is still not good enough for my
> purpose (I need to find theta parameter).
>
> Best wishes
> Andrzej Szary
>
>
> >> On 15 February 2017 at 13:14, Fabricio Ferrari <
> fabricio.ferrari at gmail.com> wrote:
> >>
> >> Dear Andrzej
> >>
> >> here it seemed to work. Maybe the version you are using.
> >>
> >> Here I have:
> >>
> >>In [2]: np.__version__
> >>Out[2]: '1.11.1'
> >>
> >>In [4]: astropy.__version__
> >>Out[4]: u'1.2.1'
> >>
> >>My output:
> >>
> >>DATA:
> >>Model: Ellipse2D
> >>Inputs: (u'x', u'y')
> >>Outputs: (u'z',)
> >>Model set size: 1
> >>Parameters:
> >> amplitude x_0 y_0 a b theta
> >> --------- --- --- --- --- -----
> >> 1.0 0.0 0.0 2.0 1.0 0.0
> >>
> >>
> >>MODEL:
> >>Model: Ellipse2D
> >>Inputs: (u'x', u'y')
> >>Outputs: (u'z',)
> >>Model set size: 1
> >>Parameters:
> >> amplitude x_0 y_0 a b theta
> >> -------------- --- --- --- --- -----
> >> 0.338211382114 0.1 0.1 3.0 2.0 0.2
> >>
> >>p.s. thank you for letting me know about the Ellipse2D fit, that I
> indeed need and did not >>know about.
> >>
> >>Regards
> >>Fabricio
>
>
> On 8 February 2017 at 16:40, Andrzej Szary <aszary at gmail.com> wrote:
>
>> Hello,
>>
>> I am struggling with fitting an ellipse using Ellipse2D model. The
>> fitting procedure does not seem to have any effect on ellipse parameters
>> (maybe except amplitude).
>>
>> Please find a test code at the end of this email.
>>
>> Best wishes,
>> Andrzej Szary
>>
>>
>> import numpy as np
>> from astropy.modeling import models, fitting
>> import matplotlib.pyplot as pl
>>
>> # data
>> num = 100
>> x, y = np.meshgrid(np.linspace(-5., 5., num), np.linspace(-5, 5, num))
>> e0 = models.Ellipse2D(amplitude=1., x_0=0., y_0=0., a=2, b=1, theta=0.)
>> z0 = e0(x, y)
>> print 'DATA:\n', e0, '\n\n'
>>
>> # initial model
>> ei = models.Ellipse2D(amplitude=1., x_0=0.1, y_0=0.1, a=3, b=2, theta=0.2)
>> fi = fitting.LevMarLSQFitter()
>>
>> # fitted model?
>> e1 = fi(ei, x, y, z0)
>> z1 = e1(x, y)
>> print 'MODEL:\n', e1, '\n\n'
>>
>> pl.imshow(z0, extent=[-5, 5, -5, 5], alpha=0.5)
>> pl.imshow(z1, extent=[-5, 5, -5, 5], alpha=0.2)
>> pl.show()
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20170215/c29cf1e2/attachment.html>
More information about the AstroPy
mailing list