private void ComputeWcs()
{
if (header.ContainsKey("CROTA2"))
{
rotation = Convert.ToDouble(header["CROTA2"].Trim());
hasRotation = true;
}
if (header.ContainsKey("CDELT1"))
{
scaleX = Convert.ToDouble(header["CDELT1"].Trim());
if (header.ContainsKey("CDELT2"))
{
scaleY = Convert.ToDouble(header["CDELT2"].Trim());
hasScale = true;
}
}
if (header.ContainsKey("CRPIX1"))
{
referenceX = Convert.ToDouble(header["CRPIX1"].Trim()) - 1;
if (header.ContainsKey("CRPIX2"))
{
referenceY = Convert.ToDouble(header["CRPIX2"].Trim()) - 1;
hasPixel = true;
}
}
bool galactic = false;
bool tan = false;
if (header.ContainsKey("CTYPE1"))
{
if (header["CTYPE1"].Contains("GLON-"))
{
galactic = true;
tan = true;
}
if (header["CTYPE2"].Contains("GLAT-"))
{
galactic = true;
tan = true;
}
if (header["CTYPE1"].Contains("-TAN"))
{
tan = true;
}
if (header["CTYPE1"].Contains("-SIN"))
{
tan = true;
}
}
if (!tan)
{
throw new System.Exception("Only TAN projected images are supported: ");
}
hasSize = true;
if (header.ContainsKey("CRVAL1"))
{
centerX = Convert.ToDouble(header["CRVAL1"].Trim());
if (header.ContainsKey("CRVAL2"))
{
centerY = Convert.ToDouble(header["CRVAL2"].Trim());
hasLocation = true;
}
}
if (galactic)
{
double[] result = Earth3d.GalactictoJ2000(centerX, centerY);
centerX = result[0];
centerY = result[1];
}
if (header.ContainsKey("CD1_1") && header.ContainsKey("CD1_2") &&
header.ContainsKey("CD2_1") && header.ContainsKey("CD2_2"))
{
cd1_1 = Convert.ToDouble(header["CD1_1"].Trim());
cd1_2 = Convert.ToDouble(header["CD1_2"].Trim());
cd2_1 = Convert.ToDouble(header["CD2_1"].Trim());
cd2_2 = Convert.ToDouble(header["CD2_2"].Trim());
if (!hasRotation)
{
CalculateRotationFromCD();
}
if (!hasScale)
{
CalculateScaleFromCD();
}
hasScale = true;
hasRotation = true;
}
ValidWcs = hasScale && hasRotation && hasPixel && hasLocation;
}
private void GoToRADec_Click(object sender, EventArgs e)
{
int index = coordinateType.SelectedIndex;
double ra = 0;
double dec = 0;
bool raValid = false;
bool decValid = false;
switch (index)
{
case 0: // Equitorial
{
ra = Coordinates.ParseRA(raText.Text, false);
dec = Coordinates.ParseDec(decText.Text);
raValid = Coordinates.ValidateRA(raText.Text);
decValid = Coordinates.ValidateDec(decText.Text);
}
break;
case 2: // Galactic
{
double l = Coordinates.Parse(raText.Text);
double b = Coordinates.ParseDec(decText.Text);
raValid = Coordinates.Validate(raText.Text);
decValid = Coordinates.ValidateDec(decText.Text);
if (raValid && decValid)
{
double[] result = Earth3d.GalactictoJ2000(l, b);
ra = result[0] / 15;
dec = result[1];
}
}
break;
case 3: // Ecliptic
{
double l = Coordinates.Parse(raText.Text);
double b = Coordinates.ParseDec(decText.Text);
raValid = Coordinates.Validate(raText.Text);
decValid = Coordinates.ValidateDec(decText.Text);
if (raValid && decValid)
{
AstroCalc.AstroRaDec radec = AstroCalc.AstroCalc.EclipticToJ2000(l, b, SpaceTimeController.JNow);
ra = radec.RA;
dec = radec.Dec;
}
}
break;
case 4: // Geo
{
ra = -Coordinates.Parse(raText.Text) / 15;
dec = Coordinates.ParseDec(decText.Text);
raValid = Coordinates.Validate(raText.Text);
decValid = Coordinates.ValidateDec(decText.Text);
}
break;
case 1: // alt/az
{
double az = Coordinates.Parse(raText.Text);
double alt = Coordinates.ParseDec(decText.Text);
raValid = Coordinates.Validate(raText.Text);
decValid = Coordinates.ValidateDec(decText.Text);
Coordinates radec = Coordinates.HorizonToEquitorial(Coordinates.FromLatLng(alt, az), SpaceTimeController.Location, SpaceTimeController.Now);
ra = radec.RA;
dec = radec.Dec;
}
break;
}
if (raValid && decValid)
{
if (Earth3d.MainWindow.SolarSystemMode)
{
Vector3d pnt = Coordinates.GeoTo3dDouble(dec, Coordinates.Parse(raText.Text));
pnt = Vector3d.TransformCoordinate(pnt, Planets.EarthMatrix);
pnt.Normalize();
Vector2d radec = Coordinates.CartesianToLatLng(pnt);
Earth3d.MainWindow.TargetLat = radec.Y;
Earth3d.MainWindow.TargetLong = radec.X - 90;
}
else
{
Earth3d.MainWindow.GotoTargetRADec(ra, dec, true, false);
}
}
}