WCSTools programs

sub1 Examples

1. Search with results to standard output
2. Search with results to a file
3. Find closest star
4. List selected stars from catalog
5. Match stars in input catalog
6. Search over the web with results to standard output
7. Return objects found on 4 or more plates

Example 1. Search with results to standard output

Find the 10 brightest (in R2) USNO-B1.0 stars in a 1-degree square box centered on 10:00:00 +30:00:00 J2000 and sort them by right ascension.
sub1 -hn 10 -s m4 -r -1800 10:00:00 +30:00:00 J2000
produces the following, where PM is the probability that there is a proper motion from 0=none to 9=probable and ID is the number of plates, out of 5, on which the star was found (0 means it is bright and from the Tycho-2 catalog):
sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
   USNO-B1.0 10:00:00.000 +30:00:00.00 J2000 +- +00:34:38.5 +00:30:00.0 at epoch 2000.0000
10 / 6370 USNO-B1.0 Sources (MagR2 brighter than 9.66)
USNO_B1_number  RA2000       Dec2000     MagB1  MagR1  MagB2  MagR2  MagN  PM NI SG  Arcsec
1196.0179331 09:59:36.214 +29:38:42.83    6.99   5.16   6.30   5.10   4.61 10  0 12  1314.15
1195.0176556 09:58:29.831 +29:32:03.95    9.21   7.47   8.55   7.41   6.94 10  0 12  2046.36
1204.0178692 09:59:07.128 +30:25:18.94    9.78   7.51   8.91   7.44   6.76 10  0 12  1666.40
1198.0180130 09:59:15.755 +29:52:46.95   10.70   8.50   9.85   8.43   7.79 10  0 12   719.92
1199.0177782 09:58:19.171 +29:55:34.88   10.25   8.70   9.69   8.64   8.21 10  0 12  1336.84
1202.0177596 10:00:52.275 +30:12:16.50   10.72   9.14  10.14   9.09   8.65 10  0 12  1001.31
1200.0179599 10:01:49.216 +30:00:35.48   11.38   9.50  10.66   9.44   8.93 10  0 12  1419.13
1196.0179473 10:00:51.334 +29:41:03.72   10.45   9.51  10.16   9.47   9.19 10  0 12  1318.04
1200.0179274 09:59:19.887 +30:00:27.47   10.36   9.62  10.13   9.60   9.37 10  0 12   521.78
1204.0178926 10:00:24.300 +30:25:51.43   11.32   9.71  10.73   9.66   9.22 10  0 12  1583.08

Example 2. Search with results to a file

Find all USNO A stars within 3 arcminutes of M51 and write RA-sorted results to a file
sub1 -ts r -wr 180 -o m51 13:27:48.311 +47:27:28.98 B1950
produces file m51.ub1:
catalog USNO-B1.0
ra      13:27:48.311
dec     +47:27:28.98
radecsys        fk4
equinox 1950.0000
epoch   1950.0000
rpmunit mas/year
dpmunit mas/year
radsec  180.00
catsort ra
program sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
usnob1_id       ra1950          dec1950         magb1   magr1   magb2   magr2   magn    pm      ni      sg      pmra    pmdec   arcsec
------------    ------------    ------------    -----   -----   -----   -----   -----   --      --      --      ------  ------  ------
1371.0282379    13:27:31.970    +47:26:32.41    16.32   99.99   99.99   99.99   17.74    0       2      10         0.0     0.0  175.14
1371.0282384    13:27:33.248    +47:25:59.82    19.19   18.00   99.99   99.99   18.73    7       3       7        39.3    61.9  176.91
1371.0282382    13:27:33.337    +47:26:01.32    17.40   99.99   99.99   99.99   18.14    0       2      10         0.0     0.0  175.37
1372.0290336    13:27:33.715    +47:28:25.46    99.99   16.92   99.99   99.99   17.54    0       2      10         0.0     0.0  158.42
1371.0282386    13:27:34.111    +47:27:19.48    19.70   99.99   99.99   99.99   18.29    0       2      10         0.0     0.0  144.34
 . . .

1372.0290378    13:28:02.568    +47:29:03.42    18.15   17.74   99.99   99.99   99.99    0       2       1         0.0     0.0  172.67
1371.0282503    13:28:03.011    +47:25:57.51    19.18   19.65   17.83   18.33   99.99    0       4       2         0.0     0.0  174.94
1371.0282504    13:28:03.173    +47:25:59.04    19.85   99.99   19.41   99.99   99.99    0       2       1         0.0     0.0  175.56
1372.0290380    13:28:04.408    +47:27:48.71    99.99   19.25   99.99   99.99   18.72    0       2      11         0.0     0.0  164.43
1372.0290381    13:28:05.181    +47:27:42.70    99.99   99.99   99.99   19.75   18.51    0       2      13         0.0     0.0  171.63
(56 stars in file)

Example 3. Find closest star

Find the closest USNO B-1.0 star to 10:00:00 +50:00:00 J2000
sub1 -a 10:00:00 +50:00:00 J2000
produces the following:
1399.0209660 10:00:00.086 +49:59:45.73   21.05  19.91  99.99  99.99  99.99  0  2   14.29

Example 4. List selected stars from catalog

List USNO-B1.0 Catalog stars 1350.00000001 through 1350.00000010
sub1 -s n -h 1350.0000001-1350.0000010x0.0000001
produces the following:
sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
USNO_B1_number  RA2000       Dec2000     MagB1  MagR1  MagB2  MagR2  MagN  PM NI SG
1350.0000001 00:00:00.001 +45:02:44.55   20.41  99.99  21.06  20.30  19.04  0  4 10
1350.0000002 00:00:00.063 +45:05:28.45   19.78  18.67  19.30  18.11  17.80  0  5  7
1350.0000003 00:00:00.112 +45:05:07.60   20.41  18.46  20.60  18.46  17.83  9  5  8
1350.0000004 00:00:00.223 +45:01:04.77   99.99  20.03  21.12  19.62  18.50  0  4  7
1350.0000005 00:00:00.350 +45:04:10.92   99.99  20.01  21.12  20.30  99.99  9  3  3
1350.0000006 00:00:00.469 +45:03:57.53   99.99  99.99  21.09  20.46  99.99  0  2  6
1350.0000007 00:00:00.531 +45:02:32.74   19.48  16.68  18.26  16.53  15.22  9  5  8
1350.0000008 00:00:00.563 +45:00:16.87   20.41  99.99  21.01  20.22  99.99  8  3  5
1350.0000009 00:00:00.757 +45:01:10.75   99.99  99.99  20.97  20.05  18.96  8  3  7
1350.0000010 00:00:01.044 +45:02:23.97   99.99  19.80  21.00  19.86  18.89  9  4  7
where the columns are USNO-B1.0 region.number, RA, Dec, B1 magnitude, R1 magnitude, B2 magnitude, R2 magnitude, N magnitude, proper motion probability (0-9, 10=Tycho), the number of plate ID's (0=Tycho-2, 2-5), and the closeness to a star PSF (1-11, Tycho-2=12).

Example 5. Match stars in input catalog

Match the stars in the input catalog test.tab to the USNO-B1.0 Catalog and list them one match per line.
sub1 -tl @test.tab
produces the following:
catalog USNO-B1.0
search  test.tab
radecsys        fk5
equinox 2000.0000
epoch   2000.0000
radsec  120.00
rpmunit mas/yr
dpmunit mas/yr
program sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
ppm_id  srch_ra         srch_dec        s_mag   s_TYPE  usnob1_id       ra              dec             magb1   magr1   magb2   magr2   magn    pmra    pmdec   pm      ni      sg       n       dra     ddec    drad   
------- ------------    ------------    -----   ----    ------------    ------------    ------------    -----   -----   -----   -----   -----   ------  ------  --      --      --       --      ---     ------- ------- -------
240020  22:15:54.555    -11:36:56.13     9.80   __      0783.0692289    22:15:54.555    -11:36:56.23    10.87    9.79   10.52    9.76    9.44     12.0     0.0  10       0      12       69        -0.00   -0.10    0.10
267313  17:47:51.210    -22:58:17.35     9.10   A0      0670.0635263    17:47:51.194    -22:58:17.30     9.22    8.70    9.05    8.68    8.53     -2.0   188.0  10       0      12       1038      -0.22    0.05    0.23
196784  13:42:07.204    -09:06:40.06     8.60   K0      0808.0255760    13:42:07.200    -09:06:40.06     9.52    8.08    9.01    8.02    7.62     10.0    24.0  10       0      12       70        -0.06    0.00    0.06
127373  10:24:43.528    +11:21:21.15    10.20   K0      1013.0184462    10:24:43.507    +11:21:21.50    10.70    8.77    9.96    8.70    8.18     -8.0     0.0  10       0      12       33        -0.31    0.35    0.47
095438  06:07:46.792    +23:28:44.76    10.30   K0      1134.0106037    06:07:46.837    +23:28:44.94    10.34    8.73    9.75    8.68    8.24     -4.0   140.0  10       0      12       157        0.61    0.18    0.64
where the columns are input catalog star name, RA, Dec, magnitude, spectral type, B1.0 region.number, RA, Dec, B1 magnitude, R1 magnitude, B2 magnitude, R2 magnitude, N magnitude, B1.0 RA proper motion, B1.0 Dec proper motion, proper motion probability (0-9, 10=Tycho), number of plate ID's (0=Tycho-2, 2-5), closeness to a star PSF (1-11, Tycho-2=12), number of stars in the 120-arcsecond radius default search region, and catalog - search dra, ddec, and dr (in arcseconds).

Example 6. Search over the web with results to standard output

Find the 10 closest USNO A-2.0 stars in a 3-arcminute circle centered on 10:00:00 +30:00:00 J2000 and sort them by right ascension.
setenv UB1_PATH server_URL/cgi-bin/scat
sub1 -hs p -vn 10 -r 180 10:00:00 +30:00:00 J2000
produces the following:
sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
Center:  10:00:00.000   +30:00:00.00 J2000
Radius: +00:03:00.00
   USNO-B1.0 10:00:00.000 +30:00:00.00 J2000 r= +00:03:00.0 at epoch 2000.0000
http://tdc-www.harvard.edu/cgi-bin/scat?catalog=ub1&ra=150.0000000&dec=30.0000000&system=J2000&rad=180.000&sort=distance&epoch=2000.00000&nstar=10
SearchLim: RA: 09:59:46.144 - 10:00:13.856  Dec: +29:57:00.00 - +30:03:00.00
RefLim: RA: 09:59:46.137 - 10:00:13.863  Dec: +29:57:00.00 - +30:03:00.00
TABREAD: Catalog http://tdc-www.harvard.edu/cgi-bin/scat : 10 / 11 / 10 found
10 USNO-B1.0 Sources
USNO_B1_number  RA2000       Dec2000     MagB1  MagR1  MagB2  MagR2  MagN  PM NI SG  Arcsec
1199.0178018 09:59:57.821 +29:59:42.02   18.52  19.10  19.07  19.08  17.89  0  5  0  33.53
1199.0178017 09:59:57.361 +29:59:59.15   99.99  18.46  19.89  18.77  17.51  0  4  0  34.29
1200.0179369 10:00:03.012 +30:00:08.37   18.56  20.01  19.85  19.67  99.99  9  4  0  40.01
1200.0179371 10:00:04.152 +30:00:05.71   99.99  19.96  20.24  99.99  99.99  0  2  0  54.24
1200.0179370 10:00:03.634 +30:00:32.17   99.99  20.56  21.52  20.09  99.99  8  3  0  57.12
1199.0178024 10:00:01.693 +29:58:54.63   23.38  99.99  21.53  99.99  99.99  0  2  0  68.97
1199.0178027 10:00:04.420 +29:59:18.96   24.29  19.81  99.99  99.99  99.99  0  2  0  70.58
1200.0179355 09:59:55.628 +30:00:56.79   99.99  20.54  21.50  99.99  99.99  0  2  0  80.31
1199.0178021 09:59:59.035 +29:58:39.72   99.99  19.41  20.63  18.87  18.33  9  4  0  81.25
1200.0179363 10:00:00.151 +30:01:21.71   99.99  20.30  20.93  20.35  99.99  9  3  0  81.73

Example 7. Return objects found on 4 or more plates

Find objects within 2 arcminutes of 10:00:0 +10:00:00 J2000 which appear on at least 4 of the surveyed plates and sort them by the second epoch R magnitude.
sub1 -s m4 -hr 120 10:00 10:00 J2000 minid=4
produces the following:
sub1 WCSTools 3.6.9, 6 July 2007, Doug Mink SAO
   USNO-B1.0 10:00:00.000 +10:00:00.00 J2000 r= +00:02:00.0 at epoch 2000.0000
USNO_B1_number  RA2000       Dec2000     MagB1  MagR1  MagB2  MagR2  MagN  PM NI SG  Arcsec
0999.0188143 10:00:03.992 +09:59:00.61   15.60  14.03  15.12  14.08  13.47  9  5  0   83.69
1000.0183182 10:00:04.883 +10:01:22.79   20.18  17.91  19.48  17.62  17.08  9  5  0  109.80
1000.0183184 10:00:05.283 +10:01:19.44   18.31  17.66  17.53  17.73  17.80  0  5  0  111.36
0999.0188141 10:00:03.517 +09:58:42.26   19.64  18.51  18.64  18.26  18.18  0  5  0   93.50
0999.0188137 10:00:01.906 +09:58:24.17   99.99  19.34  20.01  18.97  18.19  0  4  0   99.88
1000.0183168 09:59:57.108 +10:00:42.33   99.99  19.81  21.41  19.45  18.65  0  4  0   60.14
1000.0183153 09:59:54.415 +10:01:12.16   20.59  19.91  20.88  19.65  99.99  0  4  0  109.61
0999.0188146 10:00:05.085 +09:58:42.16   99.99  19.63  21.01  19.78  18.98  9  4  0  108.18


Telescope Data Center

Last updated 9 July 2007 by Doug Mink

[WCSTools] [WCSTools programs] [scat] [sub1] [sua2] [sua2 examples]