Astrometry.net Index files, 6000-series. These index files, which only cover a narrow range of scales, were built from a GALEX near-UV catalog. The catalog is "GUVCat" by L.Bianchi et al, described in https://arxiv.org/pdf/1704.05903.pdf It is *not* full-sky -- the GALEX satellite had pointing restrictions that prevented it from observing much of the galactic plane and other parts of the sky containing bright objects. About 25,000 square degrees are covered, with ~83 million sources total. ~80 million of these have NUV mags and are not flagged as artifacts (with artifact bits 0x4 or 0x2, as recommended by the authors). The 6000-series is sorted based on the NUV mag. The 6100-series is cut to sources with FUV mags (about 10% of the sources) and sorted based on FUV mag. The making of: - Downloaded from the author's website http://dolomiti.pha.jhu.edu/uvsky/GUVcat/GUVcat_AIS.html - Cut via: fns = glob('/global/cscratch1/sd/dstn/guvcat/*.fits') fns.sort() TT = [] for fn in fns: T = fits_table(fn, columns=['ra','dec','fuv_mag','nuv_mag', 'fuv_artifact', 'nuv_artifact']) T.cut((T.fuv_artifact & 0x6) == 0) T.cut((T.nuv_artifact & 0x6) == 0) T.delete_column('nuv_artifact') T.delete_column('fuv_artifact') TT.append(T) T = merge_tables(TT) T.writeto('/global/cscratch1/sd/dstn/guvcat/guvcat-cut.fits') FUV = T[T.fuv_mag > -500] FUV.writeto('/global/cscratch1/sd/dstn/guvcat/guvcat-fuv.fits') Then hpsplit via hpsplit -o guvcat-hp%02i.fits -n 1 -m 1 -c ra -c dec -c fuv_mag -c nuv_mag guvcat-cut.fits # FUV measured: hpsplit -o guvcat-fuv-hp%02i.fits -n 1 -m 1 -c ra -c dec -c fuv_mag -c nuv_mag guvcat-fuv.fits Indexing for scale in 6 5 4; do S=$(printf %02i $scale) for hp in 0 1 2 3 4 5 6 7 8 9 10 11; do HP=$(printf %02i $hp) build-astrometry-index -i guvcat-hp${HP}.fits -o index-60${S}-${HP}.fits -P $scale -S nuv_mag -H $hp -s 1 -r 1 -I 60${S} & build-astrometry-index -i guvcat-fuv-hp${HP}.fits -o index-61${S}-${HP}.fits -P $scale -S fuv_mag -H $hp -s 1 -r 1 -I 61${S} & done wait done