private void OnValidate(object sender, EventArgs e)
{
try
{
double lat = 0.0, lon = 0.0, x = 0.0, y = 0.0, x1 = 0.0, y1 = 0.0, azi = 0.0, rk = 0.0;
AzimuthalEquidistant azimuthal = new AzimuthalEquidistant(m_geodesic);
azimuthal = new AzimuthalEquidistant();
azimuthal.Forward(32.0, -86.0, 33.0, -87.0, out x, out y, out azi, out rk);
azimuthal.Forward(32.0, -86.0, 33.0, -87.0, out x1, out y1);
if (x != x1 || y != y1)
{
throw new Exception("Error in AzimuthalEquidistant.Forward");
}
azimuthal.Reverse(32.0, -86.0, x, y, out lat, out lon, out azi, out rk);
azimuthal.Reverse(32.0, -86.0, x, y, out x1, out y1);
if (x1 != lat || y1 != lon)
{
throw new Exception("Error in AzimuthalEquidistant.Reverse");
}
CassiniSoldner cassini = new CassiniSoldner(32.0, -86.0, m_geodesic);
cassini = new CassiniSoldner(32.0, -86.0);
cassini.Reset(31.0, -87.0);
cassini.Forward(32.0, -86.0, out x, out y, out azi, out rk);
cassini.Forward(32.0, -86.0, out x1, out y1);
if (x != x1 || y != y1)
{
throw new Exception("Error in CassiniSoldner.Forward");
}
cassini.Reverse(x, y, out lat, out lon, out azi, out rk);
cassini.Reverse(x, y, out x1, out y1);
if (x1 != lat || y1 != lon)
{
throw new Exception("Error in CassiniSoldner.Reverse");
}
Gnomonic gnomonic = new Gnomonic(m_geodesic);
gnomonic = new Gnomonic();
gnomonic.Forward(32.0, -86.0, 31.0, -87.0, out x, out y, out azi, out rk);
gnomonic.Forward(32.0, -86.0, 31.0, -87.0, out x1, out y1);
if (x != x1 || y != y1)
{
throw new Exception("Error in Gnomonic.Forward");
}
gnomonic.Reverse(32.0, -86.0, x, y, out lat, out lon, out azi, out rk);
gnomonic.Reverse(32.0, -86.0, x, y, out x1, out y1);
if (x1 != lat || y1 != lon)
{
throw new Exception("Error in Gnomonic.Reverse");
}
}
catch (Exception xcpt)
{
MessageBox.Show(xcpt.Message, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
MessageBox.Show("No errors detected", "OK", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
static void Main(string[] args)
{
try {
Geodesic geod = new Geodesic(); // WGS84
const double lat0 = 48 + 50 / 60.0, lon0 = 2 + 20 / 60.0; // Paris
AzimuthalEquidistant proj = new AzimuthalEquidistant(geod);
{
// Sample forward calculation
double lat = 50.9, lon = 1.8; // Calais
double x, y;
proj.Forward(lat0, lon0, lat, lon, out x, out y);
Console.WriteLine(String.Format("X: {0} Y: {1}", x, y));
}
{
// Sample reverse calculation
double x = -38e3, y = 230e3;
double lat, lon;
proj.Reverse(lat0, lon0, x, y, out lat, out lon);
Console.WriteLine(String.Format("Latitude: {0} Longitude: {1}", lat, lon));
}
}
catch (GeographicErr e) {
Console.WriteLine(String.Format("Caught exception: {0}", e.Message));
}
}
static void Main(string[] args)
{
try {
Geodesic geod = new Geodesic(); // WGS84
const double lat0 = 48 + 50/60.0, lon0 = 2 + 20/60.0; // Paris
AzimuthalEquidistant proj = new AzimuthalEquidistant(geod);
{
// Sample forward calculation
double lat = 50.9, lon = 1.8; // Calais
double x, y;
proj.Forward(lat0, lon0, lat, lon, out x, out y);
Console.WriteLine( String.Format("X: {0} Y: {1}", x, y ) );
}
{
// Sample reverse calculation
double x = -38e3, y = 230e3;
double lat, lon;
proj.Reverse(lat0, lon0, x, y, out lat, out lon);
Console.WriteLine( String.Format("Latitude: {0} Longitude: {1}", lat, lon ) );
}
}
catch (GeographicErr e) {
Console.WriteLine( String.Format( "Caught exception: {0}", e.Message ) );
}
}
private void ConvertAzimuthalEquidistant()
{
double lat0 = Double.Parse(m_lat0TextBox.Text);
double lon0 = Double.Parse(m_lon0TextBox.Text);
double lat = 0.0, lon = 0.0, x = 0.0, y = 0.0, azi = 0.0, rk = 0.0;
switch (m_functionComboBox.SelectedIndex)
{
case 0:
lat = Double.Parse(m_latitudeTextBox.Text);
lon = Double.Parse(m_longitudeTextBox.Text);
m_azimuthal.Forward(lat0, lon0, lat, lon, out x, out y, out azi, out rk);
m_xTextBox.Text = x.ToString();
m_yTextBox.Text = y.ToString();
break;
case 1:
x = Double.Parse(m_xTextBox.Text);
y = Double.Parse(m_yTextBox.Text);
m_azimuthal.Reverse(lat0, lon0, x, y, out lat, out lon, out azi, out rk);
m_latitudeTextBox.Text = lat.ToString();
m_longitudeTextBox.Text = lon.ToString();
break;
}
m_azimuthTextBox.Text = azi.ToString();
m_scaleTextBox.Text = rk.ToString();
}
private void OnValidate(object sender, EventArgs e)
{
try
{
double lat = 0.0, lon = 0.0, x = 0.0, y = 0.0, x1 = 0.0, y1 = 0.0, azi = 0.0, rk = 0.0;
AzimuthalEquidistant azimuthal = new AzimuthalEquidistant(m_geodesic);
azimuthal = new AzimuthalEquidistant();
azimuthal.Forward(32.0, -86.0, 33.0, -87.0, out x, out y, out azi, out rk);
azimuthal.Forward(32.0, -86.0, 33.0, -87.0, out x1, out y1);
if (x != x1 || y != y1)
throw new Exception("Error in AzimuthalEquidistant.Forward");
azimuthal.Reverse(32.0, -86.0, x, y, out lat, out lon, out azi, out rk);
azimuthal.Reverse(32.0, -86.0, x, y, out x1, out y1);
if ( x1 != lat || y1 != lon )
throw new Exception("Error in AzimuthalEquidistant.Reverse");
CassiniSoldner cassini = new CassiniSoldner(32.0, -86.0, m_geodesic);
cassini = new CassiniSoldner(32.0, -86.0);
cassini.Reset(31.0, -87.0);
cassini.Forward(32.0, -86.0, out x, out y, out azi, out rk);
cassini.Forward(32.0, -86.0, out x1, out y1);
if (x != x1 || y != y1)
throw new Exception("Error in CassiniSoldner.Forward");
cassini.Reverse(x, y, out lat, out lon, out azi, out rk);
cassini.Reverse(x, y, out x1, out y1);
if (x1 != lat || y1 != lon)
throw new Exception("Error in CassiniSoldner.Reverse");
Gnomonic gnomonic = new Gnomonic(m_geodesic);
gnomonic = new Gnomonic();
gnomonic.Forward(32.0, -86.0, 31.0, -87.0, out x, out y, out azi, out rk);
gnomonic.Forward(32.0, -86.0, 31.0, -87.0, out x1, out y1);
if (x != x1 || y != y1)
throw new Exception("Error in Gnomonic.Forward");
gnomonic.Reverse(32.0, -86.0, x, y, out lat, out lon, out azi, out rk);
gnomonic.Reverse(32.0, -86.0, x, y, out x1, out y1);
if (x1 != lat || y1 != lon)
throw new Exception("Error in Gnomonic.Reverse");
}
catch (Exception xcpt)
{
MessageBox.Show(xcpt.Message, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
MessageBox.Show("No errors detected", "OK", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
