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);
}
