agpy 0.1 documentation

Source code for agpy.reg_gal2cel

try:
    import coords
    import pyregion
except ImportError:
    print "reg_gal2cel requires coords & pyregion"
import posang

[docs]def gal2cel(regfile): """ Converts a region file from galactic to celestial coordinates including position angle reference from the center of the box (right now only works on box regions) Requires pyregion with the ShapeList.write() function implemented... not clear if that exists in 1.0 """ reg = pyregion.open(regfile) for R in reg: if R.name == 'box': x,y,dx,dy,angle = R.coord_list posn = coords.Position([x,y],system='galactic') ra,dec = posn.j2000() newang = posang.posang(x-dx,y,x+dx,y,system='galactic') coord_list = [ra,dec,dx,dy,angle-newang-90] R.coord_format = 'fk5' R.coord_list = coord_list R.params = coord_list reg.write(regfile[:-4] + "_fk5.reg")