public static void Main(string[] args)
{
try
{
/* -------------------------------------------------------------------- */
/* Initialize srs */
/* -------------------------------------------------------------------- */
SpatialReference src = new SpatialReference("");
src.ImportFromProj4("+proj=latlong +datum=WGS84 +no_defs");
Console.WriteLine( "SOURCE IsGeographic:" + src.IsGeographic() + " IsProjected:" + src.IsProjected() );
SpatialReference dst = new SpatialReference("");
dst.ImportFromProj4("+proj=somerc +lat_0=47.14439372222222 +lon_0=19.04857177777778 +x_0=650000 +y_0=200000 +ellps=GRS67 +units=m +no_defs");
Console.WriteLine( "DEST IsGeographic:" + dst.IsGeographic() + " IsProjected:" + dst.IsProjected() );
/* -------------------------------------------------------------------- */
/* making the transform */
/* -------------------------------------------------------------------- */
CoordinateTransformation ct = new CoordinateTransformation(src, dst);
double[] p = new double[3];
p[0] = 19; p[1] = 47; p[2] = 0;
ct.TransformPoint(p);
Console.WriteLine("x:" + p[0] + " y:" + p[1] + " z:" + p[2]);
ct.TransformPoint(p, 19.2, 47.5, 0);
Console.WriteLine("x:" + p[0] + " y:" + p[1] + " z:" + p[2]);
}
catch (Exception e)
{
Console.WriteLine("Error occurred: " + e.Message);
System.Environment.Exit(-1);
}
}
public static string ToCoordinateSystem(string srs, bool isName)
{
if (string.IsNullOrEmpty(srs))
return srs;
string result = null;
if (_srCache.TryGetValue(srs, out result))
return result;
using (SpatialReference sr = new SpatialReference(null))
{
if (isName)
{
int srsId;
if (int.TryParse(srs, out srsId))
sr.ImportFromEPSG(srsId);
else
{
// TODO
}
}
else
{
if (srs.StartsWith("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0"))
{ // +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs
//sr.ImportFromProj4(srs);
sr.ImportFromEPSG(900913);
}
else
{
sr.ImportFromProj4(srs);
}
}
sr.ExportToWkt(out result);
if (result != null)
_srCache.TryAdd(srs, result);
return result;
}
}
private static NetTopologySuite.Geometries.Point Wgs84ToWgs84UpsNorth(Point location)
{
if (location.SRID != 4326)
{
throw new Exception("Unsupported coordinate system: " + location.SRID);
}
OSGeo.OSR.SpatialReference wgs84Src = new OSGeo.OSR.SpatialReference("");
wgs84Src.ImportFromProj4("+proj=longlat +datum=WGS84 +no_defs");
OSGeo.OSR.SpatialReference stereoNorthPoleDest = new OSGeo.OSR.SpatialReference("");
stereoNorthPoleDest.ImportFromProj4("+proj=stere +lat_0=90 +lat_ts=90 +lon_0=0 +k=0.994 +x_0=2000000 +y_0=2000000 +datum=WGS84 +units=m +no_defs");
OSGeo.OSR.CoordinateTransformation ct = new OSGeo.OSR.CoordinateTransformation(wgs84Src, stereoNorthPoleDest);
double[] point = new double[3];
point[0] = location.X;
point[1] = location.Y;
point[2] = location.Z;
ct.TransformPoint(point);
return(new Point(point[0], point[1]));
}