using System;

using GeoAPI.CoordinateSystems.Transformations;

using GeoAPI.Geometries;

using ProjNet.CoordinateSystems;

using ProjNet.CoordinateSystems.Transformations;                    

​

public class Program

{

    public static void Main()

    {

        CoordinateSystemFactory csFact = new CoordinateSystemFactory();

        CoordinateTransformationFactory ctFact = new CoordinateTransformationFactory();

        var wgs = csFact.CreateFromWkt("GEOGCS[\"WGS 84\", DATUM[\"WGS_1984\", SPHEROID[\"WGS 84\", 6378137, 298.257223563, AUTHORITY[\"EPSG\", \"7030\"]], AUTHORITY[\"EPSG\", \"6326\"]], PRIMEM[\"Greenwich\", 0, AUTHORITY[\"EPSG\", \"8901\"]], UNIT[\"degree\", 0.0174532925199433, AUTHORITY[\"EPSG\", \"9122\"]], AUTHORITY[\"EPSG\", \"4326\"]]");

        var rdm = csFact.CreateFromWkt("PROJCS[\"Amersfoort / RD New\",GEOGCS[\"Amersfoort\",DATUM[\"Amersfoort\",SPHEROID[\"Bessel 1841\",6377397.155,299.1528128,AUTHORITY[\"EPSG\",\"7004\"]],TOWGS84[565.4171,50.3319,465.5524,-0.398957,0.343988,-1.87740,4.0725],AUTHORITY[\"EPSG\",\"6289\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4289\"]],PROJECTION[\"Oblique_Stereographic\"],PARAMETER[\"latitude_of_origin\",52.15616055555555],PARAMETER[\"central_meridian\",5.38763888888889],PARAMETER[\"scale_factor\",0.9999079],PARAMETER[\"false_easting\",155000],PARAMETER[\"false_northing\",463000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"28992\"]]");

        ICoordinateTransformation trans = ctFact.CreateFromCoordinateSystems(wgs, rdm);

        var wgsCoordinate = new Coordinate(4.3593204, 52.0088095);

        var rdmCoord = trans.MathTransform.Transform(wgsCoordinate);

        Console.WriteLine("wgs ({0:#.####}, {1:#.####}) => rdm ({2:#}, {3:#})", wgsCoordinate.Y, wgsCoordinate.X, rdmCoord.X, rdmCoord.Y);

    }

}