ratt-ru / bullseye

A GPU-based facet imager
GNU General Public License v2.0
1 stars 1 forks source link

FITS Projections #61

Closed bennahugo closed 9 years ago

bennahugo commented 9 years ago

@o-smirnov I've noticed the SIN projection doesn't do very well at higher declinations (this is actually what causes the facets to overlap when you stitch them together. It will be even worse when you start doing w-projected images and it is giving me some headaches comparing these two...

When you create these facets and drag them into tigger you can actually see the overlap very clearly. Not sure what projection you are using in tigger? Is there any other coordinate system that is better suited for larger images at any declination?

One way around this is to make the coordinates in the facets relative to ra=0,dec=0 in the centre facet, stitch them together and then set the centre pixel coordinate. At least that way the w-projected image and the stitched faceted image have the same wrong coordinates at higher declination, but I feel this is very hacky and there ought to be a better solution?

bmerry commented 9 years ago

@bennahugo what type of faceting are you doing i.e. do all the facets lie in a single plane, or is each facet tangent to the celestial sphere at its centre?

o-smirnov commented 9 years ago

I'm surprised -- I thought SIN projection was supposed to behave the same way regardless of declination. Isn't SIN just an orthogonal projection onto the tangent plane, or am I missing something?

bmerry commented 9 years ago

I'm surprised -- I thought SIN projection was supposed to behave the same way regardless of declination.

That's also my understanding of the FITS format: distortion depends on distance away from the reference value, and shouldn't depend on declination.

bennahugo commented 9 years ago

Hmm... I agree that there should be increasing distortion further from field centre. Something must be off with my coordinate generation if it is not the SIN coordinates.

o-smirnov commented 9 years ago

This to FITS coordinate systems is like the mountain was to Moses: http://fits.gsfc.nasa.gov/fits_wcs.html

bmerry commented 9 years ago

@bennahugo, if you're in the SKA offices at some point I can work through this with you. I spent a lot of time reading the WCS specs to figure out how to correctly set up the keywords in my output files, so I'm reasonably familiar with it.

bennahugo commented 9 years ago

@bmerry @o-smirnov Thanks I will take a look and come in from tomorrow