cannot solve milky way photo

Skip to first unread message

James Sun

May 14, 2024, 10:56:52 AMMay 14
to astrometry

I tried various downsampling options, still cannot resolve this photo:

Should at least be able to see Sagittarius?


Dustin Lang

May 15, 2024, 9:52:42 AMMay 15
to James Sun, astrometry

Unfortunately, is quite bad at very wide-field images!  They should be the easiest, but the way the algorithm works causes it to not work well when there is a lot of distortion.

One thing you can do is crop out the middle section of your image, where the distortion isn't as bad.  For example, here's the middle one of your images:


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
To view this discussion on the web visit

Danko Dnevic

Jun 23, 2024, 2:49:58 AM (3 days ago) Jun 23
to astrometry

it seems that issue is not about lens distortion (AN is quite happy solving very wide fields around the zenith, or with bottom above 30 deg elevation), but with atmospheric refraction (particularly its gradient) below 45 deg elevation.

With narrow fields, above 15 degrees elevation, refraction is almost constant over the field, but as field widens in the vertical sense (especially those portrait oriented shots, like above) gradient might become extreme.
In the above shots, atmospheric refraction is approx.:
- 14 arcmin at 3 deg elevation
- 10 arcmin at 5 deg elevation
- 6 arcmin at 10 deg elevation
- 4 arcmin at 15 deg elevation
- 2 arcmin at 30 deg elevation
- 1 arcmin at 45 deg elevation
That makes vertical position error change, by cot() law, between 10(20) and 4(8) pixels in just about 12 degrees. Again, gradient is problematic, not refraction itself.
This is almost impossible to model by TAN WCS and quite hard to compensate by SIP. To make things worse, any previous camera calibration coeffs are useless.

For modeling, useful article to read is Bennett's: The Calculation of Astronomical Refraction in Marine Navigation doi:10.1017/s0373463300022037

Currently, I'm trying to make compensation, but have a hard time positioning compensator. Despite the fact that atmospheric refraction is extrinsic parameter, pre SIP placement doesn't seems correct, because forcing tweak on shifted coords and (I think) that invalidates SIP portability. Second option (the only one for me), is to try index tweaking, based on source elevation (yielded from catalog Ra/Dec and JD/lat/lon). Dustin, it would be of great help, based on your knowledge of AN internal flow, if you can give suggestion where to position that code (probably just set of vector multiplications)?

With Best Regards,

Reply all
Reply to author
0 new messages