General astronomy scripts
crossmatch.py - designed to crossmatch MSSS with TGSS and NVSS (or vice versa), so you'll need those catalogues, but with work with anything so long as the catalogue columns are in the right format. For crossmatching TGSS with other data it will automatically query TGSS through PyVO, but currently hard coded to get the position from the MSSS mosaic name. Mean astrometric positional offsets are calculated as a weighted mean, taking errors into account, so here is a version of the NVSS catalogue that I made which includes RA and DEC positional errors: https://www.dropbox.com/s/zbrvcb0la0tq3bv/NVSS_cat_alex.fits?dl=0