wktcoordsysparsertests.cs

C# 的地图开发例子（sharp map)

using System;
using System.Collections.Generic;
using System.Text;

using SharpMap.Converters.WellKnownText;
using SharpMap.CoordinateSystems;
using NUnit.Framework;

namespace UnitTests.Converters.WKT
{
	[TestFixture]
	public class WKTCoordSysParserTests
	{
		/// <summary>
		/// Parses a coordinate system WKTs
		/// </summary>
		/// <remarks><code>
		/// PROJCS["NAD83(HARN) / Texas Central (ftUS)",
		/// 	GEOGCS[
		/// 		"NAD83(HARN)",
		/// 		DATUM[
		/// 			"NAD83_High_Accuracy_Regional_Network",
		/// 			SPHEROID[
		/// 				"GRS 1980",
		/// 				6378137,
		/// 				298.257222101,
		/// 				AUTHORITY["EPSG","7019"]
		/// 			],
		///				TOWGS84[725,685,536,0,0,0,0],
		/// 			AUTHORITY["EPSG","6152"]
		/// 		],
		/// 		PRIMEM[
		/// 			"Greenwich",
		/// 			0,
		/// 			AUTHORITY["EPSG","8901"]
		/// 		],
		/// 		UNIT[
		/// 			"degree",
		/// 			0.01745329251994328,
		/// 			AUTHORITY["EPSG","9122"]
		/// 		],
		/// 		AUTHORITY["EPSG","4152"]
		/// 	],
		/// 	PROJECTION["Lambert_Conformal_Conic_2SP"],
		/// 	PARAMETER["standard_parallel_1",31.88333333333333],
		/// 	PARAMETER["standard_parallel_2",30.11666666666667],
		/// 	PARAMETER["latitude_of_origin",29.66666666666667],
		/// 	PARAMETER["central_meridian",-100.3333333333333],
		/// 	PARAMETER["false_easting",2296583.333],
		/// 	PARAMETER["false_northing",9842500.000000002],
		/// 	UNIT[
		/// 		"US survey foot",
		/// 		0.3048006096012192,
		/// 		AUTHORITY["EPSG","9003"]
		/// 	],
		/// 	AUTHORITY["EPSG","2918"]
		/// ]
		/// </code></remarks>
		[Test]
		public void ParseCoordSys()
		{
			CoordinateSystemFactory fac = new CoordinateSystemFactory();
			string wkt = "PROJCS[\"NAD83(HARN) / Texas Central (ftUS)\", GEOGCS[\"NAD83(HARN)\", DATUM[\"NAD83_High_Accuracy_Regional_Network\", SPHEROID[\"GRS 1980\", 6378137, 298.257222101, AUTHORITY[\"EPSG\", \"7019\"]], TOWGS84[725, 685, 536, 0, 0, 0, 0], AUTHORITY[\"EPSG\", \"6152\"]], PRIMEM[\"Greenwich\", 0, AUTHORITY[\"EPSG\", \"8901\"]], UNIT[\"degree\", 0.0174532925199433, AUTHORITY[\"EPSG\", \"9122\"]], AUTHORITY[\"EPSG\", \"4152\"]], PROJECTION[\"Lambert_Conformal_Conic_2SP\"], PARAMETER[\"standard_parallel_1\", 31.883333333333], PARAMETER[\"standard_parallel_2\", 30.1166666667], PARAMETER[\"latitude_of_origin\", 29.6666666667], PARAMETER[\"central_meridian\", -100.333333333333], PARAMETER[\"false_easting\", 2296583.333], PARAMETER[\"false_northing\", 9842500], UNIT[\"US survey foot\", 0.304800609601219, AUTHORITY[\"EPSG\", \"9003\"]], AUTHORITY[\"EPSG\", \"2918\"]]";
			ProjectedCoordinateSystem pcs = CoordinateSystemWktReader.Parse(wkt) as ProjectedCoordinateSystem;
			Assert.IsNotNull(pcs, "Could not parse WKT: " + wkt);			
						
			Assert.AreEqual("NAD83(HARN) / Texas Central (ftUS)", pcs.Name);
			Assert.AreEqual("NAD83(HARN)", pcs.GeographicCoordinateSystem.Name);
			Assert.AreEqual("NAD83_High_Accuracy_Regional_Network", pcs.GeographicCoordinateSystem.HorizontalDatum.Name);
			Assert.AreEqual("GRS 1980", pcs.GeographicCoordinateSystem.HorizontalDatum.Ellipsoid.Name);
			Assert.AreEqual(6378137, pcs.GeographicCoordinateSystem.HorizontalDatum.Ellipsoid.SemiMajorAxis);
			Assert.AreEqual(298.257222101, pcs.GeographicCoordinateSystem.HorizontalDatum.Ellipsoid.InverseFlattening);
			Assert.AreEqual("EPSG", pcs.GeographicCoordinateSystem.HorizontalDatum.Ellipsoid.Authority);
			Assert.AreEqual(7019, pcs.GeographicCoordinateSystem.HorizontalDatum.Ellipsoid.AuthorityCode);
			Assert.AreEqual("EPSG", pcs.GeographicCoordinateSystem.HorizontalDatum.Authority);
			Assert.AreEqual(6152, pcs.GeographicCoordinateSystem.HorizontalDatum.AuthorityCode);
			Assert.AreEqual(new Wgs84ConversionInfo(725,685,536,0,0,0,0), pcs.GeographicCoordinateSystem.HorizontalDatum.Wgs84Parameters);
			Assert.AreEqual("Greenwich", pcs.GeographicCoordinateSystem.PrimeMeridian.Name);
			Assert.AreEqual(0, pcs.GeographicCoordinateSystem.PrimeMeridian.Longitude);
			Assert.AreEqual("EPSG", pcs.GeographicCoordinateSystem.PrimeMeridian.Authority);
			Assert.AreEqual(8901, pcs.GeographicCoordinateSystem.PrimeMeridian.AuthorityCode, 8901);
			Assert.AreEqual("degree", pcs.GeographicCoordinateSystem.AngularUnit.Name);
			Assert.AreEqual(0.0174532925199433, pcs.GeographicCoordinateSystem.AngularUnit.RadiansPerUnit);
			Assert.AreEqual("EPSG", pcs.GeographicCoordinateSystem.AngularUnit.Authority);
			Assert.AreEqual(9122, pcs.GeographicCoordinateSystem.AngularUnit.AuthorityCode);
			Assert.AreEqual("EPSG", pcs.GeographicCoordinateSystem.Authority);
			Assert.AreEqual(4152, pcs.GeographicCoordinateSystem.AuthorityCode, 4152);
			Assert.AreEqual("Lambert_Conformal_Conic_2SP", pcs.Projection.ClassName, "Projection Classname");

			ProjectionParameter latitude_of_origin = pcs.Projection.GetParameter("latitude_of_origin");
			Assert.IsNotNull(latitude_of_origin);
			Assert.AreEqual(29.6666666667, latitude_of_origin.Value);
			ProjectionParameter central_meridian = pcs.Projection.GetParameter("central_meridian");
			Assert.IsNotNull(central_meridian);
			Assert.AreEqual(-100.333333333333, central_meridian.Value);
			ProjectionParameter standard_parallel_1 = pcs.Projection.GetParameter("standard_parallel_1");
			Assert.IsNotNull(standard_parallel_1);
			Assert.AreEqual(31.883333333333, standard_parallel_1.Value);
			ProjectionParameter standard_parallel_2 = pcs.Projection.GetParameter("standard_parallel_2");
			Assert.IsNotNull(standard_parallel_2);
			Assert.AreEqual(30.1166666667, standard_parallel_2.Value);
			ProjectionParameter false_easting = pcs.Projection.GetParameter("false_easting");
			Assert.IsNotNull(false_easting);
			Assert.AreEqual(2296583.333, false_easting.Value);
			ProjectionParameter false_northing = pcs.Projection.GetParameter("false_northing");
			Assert.IsNotNull(false_northing);
			Assert.AreEqual(9842500, false_northing.Value);
			
			Assert.AreEqual("US survey foot", pcs.LinearUnit.Name);
			Assert.AreEqual(0.304800609601219, pcs.LinearUnit.MetersPerUnit);
			Assert.AreEqual("EPSG", pcs.LinearUnit.Authority);
			Assert.AreEqual(9003, pcs.LinearUnit.AuthorityCode);
			Assert.AreEqual("EPSG", pcs.Authority);
			Assert.AreEqual(2918, pcs.AuthorityCode);
			Assert.AreEqual(wkt, pcs.WKT);
		}
		[Test]
		public void ParseAllWKTs()
		{
			CoordinateSystemFactory fac = new CoordinateSystemFactory();
			int parsecount = 0;
			System.IO.StreamReader sr = System.IO.File.OpenText(@"..\..\SharpMap\SRID.csv");
			string line = "";
			while (!sr.EndOfStream)
			{
				line = sr.ReadLine();
				int split = line.IndexOf(';');
				if (split > -1)
				{
					string srid = line.Substring(0, split);
					string wkt = line.Substring(split + 1);
					ICoordinateSystem cs = SharpMap.Converters.WellKnownText.CoordinateSystemWktReader.Parse(wkt) as ICoordinateSystem;
					Assert.IsNotNull(cs,"Could not parse WKT: " + wkt);
					parsecount++;
				}
			}
			sr.Close();
			Assert.AreEqual(parsecount, 2671, "Not all WKT was parsed");
		}
	}
}