The best way to do absolute astrometry for GSAOI is (probably) to align the individual images to themselves and ignore an external REFCAT until near the end and then add this absolute astrometric alignment to all of the input images. Effectively this would mean moving addReferenceCatalog() and determineAstrometricSolution() later in the recipe and having a significant rewrite of the latter primitive.
Having done adjustWCSToReference() there are two ways forward:
Produce a final stack in a stream, run detectSources() on this stack, and compute an alignment/distortion map from that, which can be applied to each input. This is memory intensive.
Simply combine the individual OBJCATs (since we use the KDTreeFitter, we might not even need to merge multiple detections of the same source in the different inputs but can leave them as multiple entries, which would make this much simpler).
I guess what we want is to take the universal set of OBJCAT (RA, DEC) coords and the REFCAT coords and convert them both to (x, y) using the reference image's WCS, and then compute the mapping between these pixel coords and add that into the gWCS pipeline of all inputs. The reference should have a WCS like
pixels -> static -T> world
while the other inputs have
pixels -> static -> variable -T> world
and the transform marked "-T>" should be identical in all inputs after adjustWCSToReference() I think. If not, more thought is required. Then we replace "-T> world" with "-> aligned -T> world". I think "-T>" is just a standard imaging projection.
The best way to do absolute astrometry for GSAOI is (probably) to align the individual images to themselves and ignore an external REFCAT until near the end and then add this absolute astrometric alignment to all of the input images. Effectively this would mean moving
addReferenceCatalog()anddetermineAstrometricSolution()later in the recipe and having a significant rewrite of the latter primitive.Having done
adjustWCSToReference()there are two ways forward:I guess what we want is to take the universal set of OBJCAT (RA, DEC) coords and the REFCAT coords and convert them both to (x, y) using the reference image's WCS, and then compute the mapping between these pixel coords and add that into the gWCS pipeline of all inputs. The reference should have a WCS like
pixels -> static -T> world
while the other inputs have
pixels -> static -> variable -T> world
and the transform marked "-T>" should be identical in all inputs after
adjustWCSToReference()I think. If not, more thought is required. Then we replace "-T> world" with "-> aligned -T> world". I think "-T>" is just a standard imaging projection.