solving a sequence of photos

80 views
Skip to first unread message

Per Ljung

unread,
Aug 17, 2025, 9:02:12 AMAug 17
to astrometry
Hi Dustin - 
I have a noob question on solving a sequence of photos. Using a static wide-angle camera I take a sequence of 1000 photos with 1sec exposures. Using image magick I can create a nice looking composite star streak image. Using celestial navigation algorithms I can compute the camera position for each photo after identifying the stars. My problem is that the solve-field results are rather noisy.

When I (locally) independently solve for each photo (typically >100 stars) this results in separate SIP2 coefficients for each photo, where the SIP coefficients are similar but noisy. The identified photo center is also noisy. Plotting the corr files field values results in noisy ~3 pixel wide star streaks. The corr index star streaks are even noisier. My CelNav algorithms solve for the lat,lon of each solved photo, but the noisy solve results in undesired noise in the lat,lon solution. 

To get less noise, can I re-use the SIP coefficients from the 1st photo for all subsequent photos?  If so then how? Or is there a better way to analyze a sequence of related photos from a static camera?

Here are some images illustrating the noise. 
streak.png
streak_field_index.png
field_index_detail.png


Thx

Sharaf Z. Muhanna

unread,
Dec 2, 2025, 9:41:39 AM (13 days ago) Dec 2
to astrometry
  Hi bro,
I’m working on a hobby project to calculate my position on Earth using star images. Could you share the celestial navigation algorithms with me, please? I’d really appreciate it. Or could you at least explain the general pipeline?
thanks  

Dustin Lang

unread,
Dec 2, 2025, 5:57:13 PM (13 days ago) Dec 2
to Per Ljung, astrometry
Hi,

Thanks for the detailed explanation and plots!

First, a couple of suggestions that *might* help get better solutions -- maybe you have already done this --
- add "--crpix-center" (make the center of the image the center of the projection)
- add "--tweak-order 3" (use 3rd-order polynomial rather than 2nd-order for SIP)
- run solve-field twice; the second time, pass the first-round WCS file as the "--verify IMAGE.wcs" argument.

To answer your original question, you should be able to re-use a SIP solution with the
"--predistort IMAGE.wcs" argument.
This has honestly not been tested a huge amount, so any reports of success or failure would be appreciated!

Best wishes,
dustin





--
You received this message because you are subscribed to the Google Groups "astrometry" group.
To unsubscribe from this group and stop receiving emails from it, send an email to astrometry+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/astrometry/aa460e98-4202-4126-bbee-ddc64b9e73c0n%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages