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]