<div dir="ltr"><div>Dear Andrzej</div><div><br></div><div>here it seemed to work. Maybe the version you are using.</div><div><br></div><div>Here I have:</div><div><br></div><div><div>In [2]: np.__version__</div><div>Out[2]: '1.11.1'</div><div><br></div><div>In [4]: astropy.__version__<br></div><div>Out[4]: u'1.2.1'</div></div><div><br></div><div><br></div><div>My output:</div><div><br></div><div><br></div><div><div><font face="monospace, monospace" size="1">DATA:</font></div><div><font face="monospace, monospace" size="1">Model: Ellipse2D</font></div><div><font face="monospace, monospace" size="1">Inputs: (u'x', u'y')</font></div><div><font face="monospace, monospace" size="1">Outputs: (u'z',)</font></div><div><font face="monospace, monospace" size="1">Model set size: 1</font></div><div><font face="monospace, monospace" size="1">Parameters:</font></div><div><font face="monospace, monospace" size="1"> amplitude x_0 y_0 a b theta</font></div><div><font face="monospace, monospace" size="1"> --------- --- --- --- --- -----</font></div><div><font face="monospace, monospace" size="1"> 1.0 0.0 0.0 2.0 1.0 0.0 </font></div><div><font face="monospace, monospace" size="1"><br></font></div><div><font face="monospace, monospace" size="1"><br></font></div><div><font face="monospace, monospace" size="1">MODEL:</font></div><div><font face="monospace, monospace" size="1">Model: Ellipse2D</font></div><div><font face="monospace, monospace" size="1">Inputs: (u'x', u'y')</font></div><div><font face="monospace, monospace" size="1">Outputs: (u'z',)</font></div><div><font face="monospace, monospace" size="1">Model set size: 1</font></div><div><font face="monospace, monospace" size="1">Parameters:</font></div><div><font face="monospace, monospace" size="1"> amplitude x_0 y_0 a b theta</font></div><div><font face="monospace, monospace" size="1"> -------------- --- --- --- --- -----</font></div><div><font face="monospace, monospace" size="1"> 0.338211382114 0.1 0.1 3.0 2.0 0.2 </font></div></div><div><br></div><div><br></div><div><br></div><div><br></div><div>p.s. thank you for letting me know about the Ellipse2D fit, that I indeed need and did not know about.</div><div><br></div><div>Regards</div><div><br></div><div>Fabricio</div><div><br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div>..-. ..-.<br>Fabricio Ferrari [<a href="http://www.ferrari.pro.br/" target="_blank">www.ferrari.pro.br</a>]<br>IMEF - FURG<br>Rio Grande, RS<br>Brasil</div></div></div>
<br><div class="gmail_quote">2017-02-09 10:00 GMT-02:00 <span dir="ltr"><<a href="mailto:astropy-request@scipy.org" target="_blank">astropy-request@scipy.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send AstroPy mailing list submissions to<br>
<a href="mailto:astropy@scipy.org">astropy@scipy.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://mail.scipy.org/mailman/listinfo/astropy" rel="noreferrer" target="_blank">https://mail.scipy.org/<wbr>mailman/listinfo/astropy</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:astropy-request@scipy.org">astropy-request@scipy.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:astropy-owner@scipy.org">astropy-owner@scipy.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of AstroPy digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Ellipse2D fit (Andrzej Szary)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Wed, 8 Feb 2017 16:40:07 +0100<br>
From: Andrzej Szary <<a href="mailto:aszary@gmail.com">aszary@gmail.com</a>><br>
To: <a href="mailto:astropy@scipy.org">astropy@scipy.org</a><br>
Subject: [AstroPy] Ellipse2D fit<br>
Message-ID:<br>
<<a href="mailto:CAOGYrkXrsAhLfdjyLjUPEQk_aYTkwD9Dn6UZOhbZgo4iaMfBAw@mail.gmail.com">CAOGYrkXrsAhLfdjyLjUPEQk_<wbr>aYTkwD9Dn6UZOhbZgo4iaMfBAw@<wbr>mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hello,<br>
<br>
I am struggling with fitting an ellipse using Ellipse2D model. The fitting<br>
procedure does not seem to have any effect on ellipse parameters (maybe<br>
except amplitude).<br>
<br>
Please find a test code at the end of this email.<br>
<br>
Best wishes,<br>
Andrzej Szary<br>
<br>
<br>
import numpy as np<br>
from astropy.modeling import models, fitting<br>
import matplotlib.pyplot as pl<br>
<br>
# data<br>
num = 100<br>
x, y = np.meshgrid(np.linspace(-5., 5., num), np.linspace(-5, 5, num))<br>
e0 = models.Ellipse2D(amplitude=1., x_0=0., y_0=0., a=2, b=1, theta=0.)<br>
z0 = e0(x, y)<br>
print 'DATA:\n', e0, '\n\n'<br>
<br>
# initial model<br>
ei = models.Ellipse2D(amplitude=1., x_0=0.1, y_0=0.1, a=3, b=2, theta=0.2)<br>
fi = fitting.LevMarLSQFitter()<br>
<br>
# fitted model?<br>
e1 = fi(ei, x, y, z0)<br>
z1 = e1(x, y)<br>
print 'MODEL:\n', e1, '\n\n'<br>
<br>
pl.imshow(z0, extent=[-5, 5, -5, 5], alpha=0.5)<br>
pl.imshow(z1, extent=[-5, 5, -5, 5], alpha=0.2)<br>
pl.show()<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://mail.scipy.org/pipermail/astropy/attachments/20170208/dd96824c/attachment-0001.html" rel="noreferrer" target="_blank">https://mail.scipy.org/<wbr>pipermail/astropy/attachments/<wbr>20170208/dd96824c/attachment-<wbr>0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<br>
AstroPy mailing list<br>
<a href="mailto:AstroPy@scipy.org">AstroPy@scipy.org</a><br>
<a href="https://mail.scipy.org/mailman/listinfo/astropy" rel="noreferrer" target="_blank">https://mail.scipy.org/<wbr>mailman/listinfo/astropy</a><br>
<br>
<br>
------------------------------<br>
<br>
End of AstroPy Digest, Vol 125, Issue 4<br>
******************************<wbr>*********<br>
</blockquote></div><br></div>