USNO-B1 implemented for astrometry

mommermi · Aug 30, 2019 · 95e2ba4 · 95e2ba4
1 parent a311b8b
commit 95e2ba4
Show file tree

Hide file tree

Showing 3 changed files with 40 additions and 2 deletions.
diff --git a/_pp_conf.py b/_pp_conf.py
@@ -140,7 +140,7 @@
 
 # list of available catalogs
 allcatalogs = ['URAT-1', '2MASS', 'SDSS-R9', 'SDSS-R13', 'APASS9', 'GAIA',
-               'PANSTARRS', 'TGAS', 'SkyMapper']
+               'PANSTARRS', 'TGAS', 'SkyMapper', 'USNO-B1']
 
 # catalog magnitude systems
 allcatalogs_magsys = {'URAT-1': 'Vega',
@@ -151,7 +151,8 @@
                       'GAIA': 'Vega',
                       'PANSTARRS': 'AB',
                       'TGAS': 'Vega',
-                      'SkyMapper': 'AB'}
+                      'SkyMapper': 'AB',
+                      'USNO-B1': 'Vega'}
 
 # ---- pipeline preferences
 # (if you don't know what you're doing, better don't mess around here)

diff --git a/catalog.py b/catalog.py
@@ -371,6 +371,41 @@ def download_catalog(self, ra_deg, dec_deg, rad_deg,
 
         # --------------------------------------------------------------------
 
+        elif self.catalogname == 'USNO-B1':
+            # astrometric and photometric catalog
+            vquery = Vizier(columns=['USNO-B1.0', 'RAJ2000', 'DEJ2000',
+                                     'e_RAJ2000', 'e_DEJ2000',
+                                     'R2mag'],
+                            column_filters={"R2mag":
+                                            ("<{:f}".format(max_mag))},
+                            row_limit=max_sources,
+                            timeout=300)
+
+            try:
+                self.data = vquery.query_region(field,
+                                                radius=rad_deg*u.deg,
+                                                catalog="I/284/out",
+                                                cache=False)[0]
+            except IndexError:
+                if self.display:
+                    print('no data available from {:s}'.format(
+                        self.catalogname))
+                logging.error('no data available from {:s}'.format(
+                    self.catalogname))
+                return 0
+
+            # rename column names using PP conventions
+            self.data.rename_column('USNO-B1.0', 'ident')
+            self.data.rename_column('RAJ2000', 'ra_deg')
+            self.data.rename_column('DEJ2000', 'dec_deg')
+            self.data.rename_column('e_RAJ2000', 'e_ra_deg')
+            self.data['e_ra_deg'].convert_unit_to(u.deg)
+            self.data.rename_column('e_DEJ2000', 'e_dec_deg')
+            self.data['e_dec_deg'].convert_unit_to(u.deg)
+            self.data['mag'] = self.data['R2mag']  # required for scamp
+
+        # --------------------------------------------------------------------
+
         elif self.catalogname == 'TGAS':
             # astrometric catalog
             vquery = Vizier(columns=['Source', 'RA_ICRS', 'DE_ICRS',

diff --git a/doc/supported.rst b/doc/supported.rst
@@ -133,6 +133,8 @@ photometric calibration:
 +------------------------+--------------------------+---------------+--------------------------+------------------------------------------------------------+
 | Gaia/TGAS (`TGAS`_)    | astrometric              | yes           | G                        | all-sky catalog, excellent astrometry, G < 12 mag          |
 +------------------------+--------------------------+---------------+--------------------------+------------------------------------------------------------+
+| USNO-B1.0              | astrometric              | yes           | R2                       | good astrometry, V < 21 mag                                |
++------------------------+--------------------------+---------------+--------------------------+------------------------------------------------------------+
 | 2MASS (`2MASS`_)       | astrometric/photometric  | yes           | J, H, Ks, K* (Vega)      | all-sky NIR catalog, good astrometry                       |
 +------------------------+--------------------------+---------------+--------------------------+------------------------------------------------------------+
 | URAT-1 (`URAT-1`_)     | astrometric/photometric  | yes (SCAMP    | g, r, i (SDSS AB);       | good coverage over the Northern hemisphere, photometry from|