Re: Fitting Circles [message #18085] |
Wed, 01 December 1999 00:00  |
Vince Hradil
Messages: 574 Registered: December 1999
|
Senior Member |
|
|
How about doing a transformation and solving a simpler problem:
What you want to fit is r0 = [(x-x0)^2+(y-y0)^2]^(0.5)
So minimize:
r0 - [(x-x0)^2+(y-y0)^2]^(0.5)
for all pairs of (x,y). This can be done using your favorite fitting
routine (I like Craig markwardt's MPFIT... routines,
http://cow.physics.wisc.edu/~craigm/idl/idl.html)
Good luck
F.N.Hatfield@Leeds.ac.uk wrote:
>
> Dear Colleagues,
>
> I am trying to fit a circle to a set of points (x,y), and from this
> determine the centre point (x0,y0) and radius r0. It is also very
> important that I obtain the standard deviations or errors in the x0,y0
> and r0.
>
> So far, I have looked at the idl routine, curvefit.pro.
> This is a useful program but doesn't seem to allow you to pass 2
> independent variables x and y.
> I was wondering if someone has written something similar in idl,
> or could suggest a way to solve this problem.
>
> Cheers
>
> Fraser Hatfield
>
> University of Leeds, UK.
>
> Sent via Deja.com http://www.deja.com/
> Before you buy.
|
|
|