def test_targets():
test_data = \
"""#name lon lat radius coord
object_1 354.36 -63.26 1.0 CEL
object_2 19.45 -17.46 1.0 CEL
#object_3 18.94 -41.05 1.0 CEL
"""
with open('targets.txt','w') as f:
f.write(test_data)
parser = ugali.utils.parser.Parser()
parser.add_coords(targets=True)
args = parser.parse_args(['-t','targets.txt'])
np.testing.assert_array_almost_equal(args.coords['lon'],[316.311,156.487],
decimal=3)
np.testing.assert_array_almost_equal(args.coords['lat'],[-51.903,-78.575],
decimal=3)
np.testing.assert_array_almost_equal(args.coords['radius'],[1.0,1.0],
decimal=3)
os.remove('targets.txt')
return args
def test_targets():
test_data = \
"""#name lon lat radius coord
object_1 354.36 -63.26 1.0 CEL
object_2 19.45 -17.46 1.0 CEL
#object_3 18.94 -41.05 1.0 CEL
"""
with open('targets.txt', 'w') as f:
f.write(test_data)
parser = ugali.utils.parser.Parser()
parser.add_coords(targets=True)
args = parser.parse_args(['-t', 'targets.txt'])
np.testing.assert_array_almost_equal(args.coords['lon'],
[316.311, 156.487],
decimal=3)
np.testing.assert_array_almost_equal(args.coords['lat'],
[-51.903, -78.575],
decimal=3)
np.testing.assert_array_almost_equal(args.coords['radius'], [1.0, 1.0],
decimal=3)
os.remove('targets.txt')
return args
#!/usr/bin/env python
import ugali.candidate.associate
import ugali.utils.parser
import numpy as np
from ugali.utils.projector import gal2cel,ang2const,ang2iau
#CATALOGS = ['McConnachie12','Rykoff14', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'WEBDA14']
CATALOGS = ['McConnachie12', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'Bica08', 'WEBDA14', 'ExtraDwarfs','ExtraClusters']
if __name__ == "__main__":
import argparse
description = "python script"
parser = ugali.utils.parser.Parser(description=description)
parser.add_coords(required=True,radius=True,targets=True)
parser.add_argument('-n','--nnearest',default=1,type=int)
opts = parser.parse_args()
catalog = ugali.candidate.associate.SourceCatalog()
for i in CATALOGS:
catalog += ugali.candidate.associate.catalogFactory(i)
for name,(glon,glat,radius) in zip(opts.names, opts.coords):
ra,dec = gal2cel(glon,glat)
iau = ang2iau(glon,glat)
const = ang2const(glon,glat)[0]
if radius <= 0: radius = None
idx1,idx2,sep = catalog.match([glon],[glat],tol=radius,nnearest=opts.nnearest)
match = catalog[idx2]
if len(match) > 0:
import numpy as np
from ugali.utils.projector import gal2cel, ang2const, ang2iau
#CATALOGS = ['McConnachie12','Rykoff14', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'WEBDA14']
CATALOGS = [
'McConnachie12', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85',
'Kharchenko13', 'Bica08', 'WEBDA14', 'ExtraDwarfs', 'ExtraClusters'
]
if __name__ == "__main__":
import argparse
description = "python script"
parser = ugali.utils.parser.Parser(description=description)
parser.add_coords(required=True, radius=True, targets=True)
parser.add_argument('-n', '--nnearest', default=1, type=int)
opts = parser.parse_args()
catalog = ugali.candidate.associate.SourceCatalog()
for i in CATALOGS:
catalog += ugali.candidate.associate.catalogFactory(i)
for name, (glon, glat, radius) in zip(opts.names, opts.coords):
ra, dec = gal2cel(glon, glat)
iau = ang2iau(glon, glat)
const = ang2const(glon, glat)[0]
if radius <= 0: radius = None
idx1, idx2, sep = catalog.match([glon], [glat],
tol=radius,
nnearest=opts.nnearest)
match = catalog[idx2]
