public static EquatorialCoords GalacticToB1950(GalacticCoords coords) { EquatorialCoords coordsout = new EquatorialCoords(); double y = Math.Atan2(Math.Sin(coords.l - deg2rad(123)), Math.Cos(coords.l - deg2rad(123)) * Math.Sin(deg2rad(27.4)) - Math.Tan(coords.b) * Math.Cos(deg2rad(27.4))); coordsout.ra = deg2rad(quad(12.25 + rad2deg(y))); double sind = Math.Sin(coords.b) * Math.Sin(deg2rad(27.4)) + Math.Cos(coords.b) * Math.Cos(deg2rad(27.4)) * Math.Cos(coords.l - deg2rad(123)); coordsout.dec = Math.Asin(sind); return(coordsout); }
public GalacticCoords B1950ToGalactic() { GalacticCoords coordsout = new GalacticCoords(); double x = Math.Atan2(Math.Sin(Utils.deg2rad(192.25) - this.ra), Math.Cos(Utils.deg2rad(192.25) - this.ra) * Math.Sin(Utils.deg2rad(27.4)) - Math.Tan(this.dec) * Math.Cos(Utils.deg2rad(27.4))); coordsout.l = Utils.deg2rad(Utils.quad(303 - Utils.rad2deg(x))); double sinB = Math.Sin(this.dec) * Math.Sin(Utils.deg2rad(27.4)) + Math.Cos(this.dec) * Math.Cos(Utils.deg2rad(27.4)) * Math.Cos(Utils.deg2rad(192.25) - this.ra); coordsout.b = Math.Asin(sinB); return(coordsout); }
public static GalacticCoords B1950ToGalactic(EquatorialCoords coords) { GalacticCoords coordsout = new GalacticCoords(); double x = Math.Atan2(Math.Sin(deg2rad(192.25) - coords.ra), Math.Cos(deg2rad(192.25) - coords.ra) * Math.Sin(deg2rad(27.4)) - Math.Tan(coords.dec) * Math.Cos(deg2rad(27.4))); coordsout.l = deg2rad(quad(303 - rad2deg(x))); double sinB = Math.Sin(coords.dec) * Math.Sin(deg2rad(27.4)) + Math.Cos(coords.dec) * Math.Cos(deg2rad(27.4)) * Math.Cos(deg2rad(192.25) - coords.ra); coordsout.b = Math.Asin(sinB); return(coordsout); }