开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > GIS编程 > 查看源码
TransformTests.aspx.cs
资源名称：SharpMap-20976.zip　[点击查看]
上传用户：sex100000
上传日期：2013-11-09
资源大小：1377k
文件大小：13k
源码类别：

GIS编程

开发平台：

C#

TransformTests.aspx.cs：源码内容
  1. using System;
  2. using System.Data;
  3. using System.Configuration;
  4. using System.Collections;
  5. using System.Web;
  6. using System.Web.Security;
  7. using System.Web.UI;
  8. using System.Web.UI.WebControls;
  9. using System.Web.UI.WebControls.WebParts;
  10. using System.Web.UI.HtmlControls;
  12. using SharpMap.Geometries;
  13. using SharpMap.CoordinateSystems;
  14. using SharpMap.CoordinateSystems.Transformations;
  16. public partial class TransformTests : System.Web.UI.Page
  17. {
  18. protected void Page_Load(object sender, EventArgs e)
  19. {
  20. TestMercator_1SP();
  21. TestMercator_2SP();
  22. TestTransverseMercator();
  23. TestLambertConicConformal_2SP();
  24. TestAlbers();
  25. TestGeocentric();
  26. }
  28. private void TestMercator_1SP()
  29. {
  30. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  32. IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("Bessel 1840", 6377397.155, 299.15281, LinearUnit.Metre);
  34. IHorizontalDatum datum = cFac.CreateHorizontalDatum("Bessel 1840", DatumType.HD_Geocentric, ellipsoid, null);
  35. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("Bessel 1840", AngularUnit.Degrees, datum,
  36. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East),
  37. new AxisInfo("Lat", AxisOrientationEnum.North));
  38.         System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>();
  39. parameters.Add(new ProjectionParameter("latitude_of_origin", 0));
  40. parameters.Add(new ProjectionParameter("central_meridian", 110));
  41. parameters.Add(new ProjectionParameter("scale_factor", 0.997));
  42. parameters.Add(new ProjectionParameter("false_easting", 3900000));
  43. parameters.Add(new ProjectionParameter("false_northing", 900000));
  44. IProjection projection = cFac.CreateProjection("Mercator_1SP", "Mercator_1SP", parameters);
  46. IProjectedCoordinateSystem coordsys = cFac.CreateProjectedCoordinateSystem("Makassar / NEIEZ", gcs, projection, LinearUnit.Metre, new AxisInfo("East", AxisOrientationEnum.East), new AxisInfo("North", AxisOrientationEnum.North));
  48. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, coordsys);
  50. SharpMap.Geometries.Point pGeo = new SharpMap.Geometries.Point(120, -3);
  51. SharpMap.Geometries.Point pUtm = trans.MathTransform.Transform(pGeo);
  52. SharpMap.Geometries.Point pGeo2 = trans.MathTransform.Inverse().Transform(pUtm);
  54. result.Text += PrintResultTable(gcs, coordsys, pGeo, pUtm, new Point(5009726.58, 569150.82), pGeo2, "Mercator_1SP test");
  55. }
  57. private void TestMercator_2SP()
  58. {
  59. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  61. IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("Krassowski 1940", 6378245.0, 298.3, LinearUnit.Metre);
  63. IHorizontalDatum datum = cFac.CreateHorizontalDatum("Krassowski 1940", DatumType.HD_Geocentric, ellipsoid, null);
  64. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("Krassowski 1940", AngularUnit.Degrees, datum, 
  65. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East),
  66. new AxisInfo("Lat", AxisOrientationEnum.North));
  67.         System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>();
  68. parameters.Add(new ProjectionParameter("latitude_of_origin", 42));
  69. parameters.Add(new ProjectionParameter("central_meridian", 51));
  70. parameters.Add(new ProjectionParameter("false_easting", 0));
  71. parameters.Add(new ProjectionParameter("false_northing", 0));
  72. IProjection projection = cFac.CreateProjection("Mercator_2SP", "Mercator_2SP", parameters);
  74. IProjectedCoordinateSystem coordsys = cFac.CreateProjectedCoordinateSystem("Pulkovo 1942 / Mercator Caspian Sea", gcs, projection, LinearUnit.Metre, new AxisInfo("East", AxisOrientationEnum.East), new AxisInfo("North", AxisOrientationEnum.North));
  76. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, coordsys);
  78. SharpMap.Geometries.Point pGeo = new SharpMap.Geometries.Point(53,53);
  79. SharpMap.Geometries.Point pUtm = trans.MathTransform.Transform(pGeo);
  80. SharpMap.Geometries.Point pGeo2 = trans.MathTransform.Inverse().Transform(pUtm);
  82. result.Text += PrintResultTable(gcs, coordsys, pGeo, pUtm, new Point(165704.29, 5171848.07), pGeo2, "Mercator_2SP test");
  83. }
  85. private void TestTransverseMercator()
  86. {
  87. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  89. IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("Airy 1830", 6377563.396, 299.32496, LinearUnit.Metre);
  91. IHorizontalDatum datum = cFac.CreateHorizontalDatum("Airy 1830", DatumType.HD_Geocentric, ellipsoid, null);
  92. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("Airy 1830", AngularUnit.Degrees, datum,
  93. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East),
  94. new AxisInfo("Lat", AxisOrientationEnum.North));
  95.         //System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>(5);
  96.         System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>();
  97. parameters.Add(new ProjectionParameter("latitude_of_origin", 49));
  98. parameters.Add(new ProjectionParameter("central_meridian", -2));
  99. parameters.Add(new ProjectionParameter("scale_factor", 0.9996012717));
  100. parameters.Add(new ProjectionParameter("false_easting", 400000));
  101. parameters.Add(new ProjectionParameter("false_northing", -100000));
  102. IProjection projection = cFac.CreateProjection("Transverse Mercator", "Transverse_Mercator", parameters);
  104. IProjectedCoordinateSystem coordsys = cFac.CreateProjectedCoordinateSystem("OSGB 1936 / British National Grid", gcs, projection, LinearUnit.Metre, new AxisInfo("East", AxisOrientationEnum.East), new AxisInfo("North", AxisOrientationEnum.North));
  106. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, coordsys);
  108. SharpMap.Geometries.Point pGeo = new SharpMap.Geometries.Point(0.5, 50.5);
  109. SharpMap.Geometries.Point pUtm = trans.MathTransform.Transform(pGeo);
  110. SharpMap.Geometries.Point pGeo2 = trans.MathTransform.Inverse().Transform(pUtm);
  112. result.Text += PrintResultTable(gcs, coordsys, pGeo, pUtm, new Point(577274.99, 69740.50), pGeo2, "Transverse Mercator test");
  113. }
  114. private void TestLambertConicConformal_2SP()
  115. {
  116. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  118. IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("Clarke 1866", 20925832.16, 294.97470, LinearUnit.USSurveyFoot);
  120. IHorizontalDatum datum = cFac.CreateHorizontalDatum("Clarke 1866", DatumType.HD_Geocentric, ellipsoid, null);
  121. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("Clarke 1866", AngularUnit.Degrees, datum,
  122. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East),
  123. new AxisInfo("Lat", AxisOrientationEnum.North));
  124.         //System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>(5);
  125.         System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>();
  126. parameters.Add(new ProjectionParameter("latitude_of_origin", 27.833333333));
  127. parameters.Add(new ProjectionParameter("central_meridian", -99));
  128. parameters.Add(new ProjectionParameter("standard_parallel_1", 28.3833333333));
  129. parameters.Add(new ProjectionParameter("standard_parallel_2", 30.2833333333));
  130. parameters.Add(new ProjectionParameter("false_easting", 2000000));
  131. parameters.Add(new ProjectionParameter("false_northing", 0));
  132. IProjection projection = cFac.CreateProjection("Lambert Conic Conformal (2SP)", "lambert_conformal_conic_2sp", parameters);
  134. IProjectedCoordinateSystem coordsys = cFac.CreateProjectedCoordinateSystem("NAD27 / Texas South Central", gcs, projection, LinearUnit.USSurveyFoot, new AxisInfo("East", AxisOrientationEnum.East), new AxisInfo("North", AxisOrientationEnum.North));
  136. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, coordsys);
  138. SharpMap.Geometries.Point pGeo = new SharpMap.Geometries.Point(-96, 28.5);
  139. SharpMap.Geometries.Point pUtm = trans.MathTransform.Transform(pGeo);
  140. SharpMap.Geometries.Point pGeo2 = trans.MathTransform.Inverse().Transform(pUtm);
  142. result.Text += PrintResultTable(gcs, coordsys, pGeo, pUtm, new Point(2963503.91, 254759.80), pGeo2, "Lambert Conic Conformal 2SP test");
  143. }
  145. private void TestAlbers()
  146. {
  147. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  149. IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("Clarke 1866", 6378206.4, 294.9786982138982,LinearUnit.USSurveyFoot);
  151. IHorizontalDatum datum = cFac.CreateHorizontalDatum("Clarke 1866", DatumType.HD_Geocentric, ellipsoid, null);
  152. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("Clarke 1866", AngularUnit.Degrees, datum,
  153. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East),
  154. new AxisInfo("Lat", AxisOrientationEnum.North));
  155.         System.Collections.ObjectModel.Collection<ProjectionParameter> parameters = new System.Collections.ObjectModel.Collection<ProjectionParameter>();
  156. parameters.Add(new ProjectionParameter("central_meridian", -96));
  157. parameters.Add(new ProjectionParameter("latitude_of_origin", 23));
  158. parameters.Add(new ProjectionParameter("standard_parallel_1", 29.5));
  159. parameters.Add(new ProjectionParameter("standard_parallel_2", 45.5));
  160. parameters.Add(new ProjectionParameter("false_easting", 0));
  161. parameters.Add(new ProjectionParameter("false_northing", 0));
  162. IProjection projection = cFac.CreateProjection("Albers Conical Equal Area", "albers", parameters);
  164. IProjectedCoordinateSystem coordsys = cFac.CreateProjectedCoordinateSystem("Albers Conical Equal Area", gcs, projection, LinearUnit.Metre, new AxisInfo("East", AxisOrientationEnum.East), new AxisInfo("North", AxisOrientationEnum.North));
  166. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, coordsys);
  168. SharpMap.Geometries.Point pGeo = new SharpMap.Geometries.Point(-75,35);
  169. SharpMap.Geometries.Point pUtm = trans.MathTransform.Transform(pGeo);
  170. SharpMap.Geometries.Point pGeo2 = trans.MathTransform.Inverse().Transform(pUtm);
  172. result.Text += PrintResultTable(gcs, coordsys, pGeo, pUtm, new Point(1885472.7,1535925), pGeo2, "Albers Conical Equal Area test");
  173. }
  175. private void TestGeocentric()
  176. {
  177. CoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory();
  179. IGeographicCoordinateSystem gcs = cFac.CreateGeographicCoordinateSystem("WGS84", AngularUnit.Degrees, HorizontalDatum.WGS84,
  180. PrimeMeridian.Greenwich, new AxisInfo("Lon", AxisOrientationEnum.East), new AxisInfo("Lat", AxisOrientationEnum.North));
  181. IGeocentricCoordinateSystem geoccs = cFac.CreateGeocentricCoordinateSystem("WGS84 geocentric", gcs.HorizontalDatum, LinearUnit.Metre, PrimeMeridian.Greenwich);
  183. ICoordinateTransformation trans = new CoordinateTransformationFactory().CreateFromCoordinateSystems(gcs, geoccs);
  185. SharpMap.Geometries.Point3D pGeo = new SharpMap.Geometries.Point3D(2.12955, 53.80939444, 73);
  186. SharpMap.Geometries.Point3D pGc = trans.MathTransform.Transform(pGeo) as SharpMap.Geometries.Point3D;
  187. SharpMap.Geometries.Point3D pGeo2 = trans.MathTransform.Inverse().Transform(pGc) as SharpMap.Geometries.Point3D;
  189. result.Text += PrintResultTable(gcs, geoccs, pGeo, pGc, new Point3D(3771793.97, 140253.34, 5124304.35), pGeo2, "Geocentric test");
  191. return;
  192. }
  193. private string PrintResultTable(ICoordinateSystem fromCoordSys, ICoordinateSystem toCoordSys, Point fromPnt, Point toPnt, Point refPnt, Point backPnt, string header)
  194. {
  195. string table = "<table style="border: 1px solid #000; margin: 10px;">";
  196. table += "<tr><td colspan="2"><h3>" + header + "</h3></td></tr>";
  197. table += "<tr><td width="130px" valign="top">Input coordsys:</td><td>" + fromCoordSys.WKT + "</td></tr>";
  198. table += "<tr><td valign="top">Output coordsys:</td><td>" + toCoordSys.WKT + "</td></tr>";
  199. table += "<tr><td>Input coordinate:</td><td>" + fromPnt.ToString() + "</td></tr>";
  200. table += "<tr><td>Ouput coordinate:</td><td>" + toPnt.ToString() + "</td></tr>";
  201. table += "<tr><td>Expected coordinate:</td><td>" + refPnt.ToString() + "</td></tr>";
  202. table += "<tr><td>Difference:</td><td>" + (refPnt-toPnt).ToString() + "</td></tr>";
  203. table += "<tr><td>Reverse transform:</td><td>" + backPnt.ToString() + "</td></tr>";
  204. table += "<tr><td>Difference:</td><td>" + (backPnt - fromPnt).ToString() + "</td></tr>";
  205. table += "</table>";
  206. return table;
  207. }
  208. }