comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Least square fitting
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Least square fitting [message #60510 is a reply to message #60326] Tue, 20 May 2008 11:44 Go to previous messageGo to previous message
MichaelT is currently offline  MichaelT
Messages: 52
Registered: May 2006
Member
> I have performed a lot of this type of astrometry, and your problem
> seems very strange to me.
The way I proceed is described here:
http://home.arcor-online.de/axel.mellinger/mwpan_web/mwpan_w eb.html
It uses the true functional form of the conversion (ra, dec) -> (x,y).
So I am not so sure if we are both talking about the same thing here.
I am not after the distortions, yet (These seem to be small. The
calculated x',y' positions deviate by only 0.5 pixels, on average,
from the real ones - when the algorithm converges).

> you really care about them? All you want is a function that maps
> the known coordinates to the calculated ones. Often what one does

What I want is a function that converts (ra, dec) into pixel
coordinates (x, y). So I am not so sure if that is what you describe
in the following. What you describe looks more like (x', y') -> (x, y)
to me.

> is to find a quadratic or cubic function that will linearize the x,y
> coordinates (i.e. so that they line up with RA, Dec)
>
> xp = x + axy + by^2 + c*x^2 + ..
> yp = y + dxy + ey^2 + f*x^2
>
> You then use least squares to determine the a,b,c.. coefficients.
> This is the 'SIP' convention discussed inhttp://ssc.spitzer.caltech.edu/postbcd/doc/shupeADASS.pdfa nd used
> e.g. by the astrometry.net software for handling distortions.

The webpage I posted above also discusses the elimination of the
distortions and uses a cubic function of the form given by you above
and as discussed in the paper.

Can I also use the cubic/quadratic function to map (ra, dec) to (x, y)
or vice versa?

This is my first attempt as an amateur astronomer to deal with these
things :-)
Many thanks for your help.

Michael
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: interact with iimage from the command line?
Next Topic: Re: problem in using function ll_to_utm.pro

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 06:11:05 PDT 2025

Total time taken to generate the page: 0.56097 seconds