<div dir="ltr"><div>Hello,</div><div><br></div><div>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).¬†</div><div><br></div><div>Please find a test code at the end of this email.</div><div><br></div><div>Best wishes,</div><div>Andrzej Szary</div><div><br></div><div><br></div><div>import numpy as np</div><div>from astropy.modeling import models, fitting</div><div>import matplotlib.pyplot as pl</div><div><br></div><div># data</div><div>num = 100</div><div>x, y = np.meshgrid(np.linspace(-5., 5., num), np.linspace(-5, 5, num))</div><div>e0 = models.Ellipse2D(amplitude=1., x_0=0., y_0=0., a=2, b=1, theta=0.)</div><div>z0 = e0(x, y)</div><div>print 'DATA:\n', e0, '\n\n'</div><div><br></div><div># initial model</div><div>ei = models.Ellipse2D(amplitude=1., x_0=0.1, y_0=0.1, a=3, b=2, theta=0.2)</div><div>fi = fitting.LevMarLSQFitter()</div><div><br></div><div># fitted model?</div><div>e1 = fi(ei, x, y, z0)</div><div>z1 = e1(x, y)</div><div>print 'MODEL:\n', e1, '\n\n'</div><div><br></div><div>pl.imshow(z0, extent=[-5, 5, -5, 5], alpha=0.5)</div><div>pl.imshow(z1, extent=[-5, 5, -5, 5], alpha=0.2)</div><div>pl.show()</div><div><br></div></div>