mcclearyj / xy-from-visit

Code to go from master catalog RA, Dec to catalog with X/Y on a particular mosaic

xy-from-visit

Code to go from master catalog RA, Dec to catalog with X/Y on a particular mosaic. Makes use of convex hull dcode kindly provided by @maximilienfranco and WCS -> xy code provided by @antonkoekemoer

Figure description: DS9 image showing a sample of objects from Marko Shuntov's (@mShuntov) Tile B9 catalog that overlap the resampled visit 64 image created by Max Franco.

To run:

python grab_XY_from_visit_runner.py -config xy_from_visit_config.yaml

Dependencies: numpy, astropy, scipy, jwst, matplotlib

To-do

• Incorporate iteration over bandpasses for matching to different coordinate files
• Include bandpass in visit_num column of output catalog

Longer term: find a cleverer way to sort through the different coordinate files, so that we are not looping over thousands of entries for many hundreds of thousands of galaxies in four bandpasses. I added a break statement to at least stop looping once a match has been found, and also implemented map(), but this problem seems ripe for parallelization, at least (maybe starmap()?)

xy-from-visit's Issues

Weird behaviour with the column checj

OS

Mac OSX Sonoma 14.1.1

Issue:

It raises exception when checking if the extra_column exists in the master catalog even though it is there

fix

for some reason ( I have no idea ) assigning the column names to mcat_cols and then comparing doesn't work however this does:

    for col in cols2save:
        if col not in catalog[config['input_catalog']['hdu']].header.values():
            raise KeyError(f"Requested column {col} not in catalog, exiting...")