GeographicTransform.cs
上传用户:sex100000
上传日期:2013-11-09
资源大小:1377k
文件大小:5k
源码类别:

GIS编程

开发平台:

C#

  1. // Copyright 2005, 2006 - Morten Nielsen (www.iter.dk)
  2. //
  3. // This file is part of SharpMap.
  4. // SharpMap is free software; you can redistribute it and/or modify
  5. // it under the terms of the GNU Lesser General Public License as published by
  6. // the Free Software Foundation; either version 2 of the License, or
  7. // (at your option) any later version.
  8. // 
  9. // SharpMap is distributed in the hope that it will be useful,
  10. // but WITHOUT ANY WARRANTY; without even the implied warranty of
  11. // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12. // GNU Lesser General Public License for more details.
  13. // You should have received a copy of the GNU Lesser General Public License
  14. // along with SharpMap; if not, write to the Free Software
  15. // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA 
  16. using System;
  17. using System.Collections.Generic;
  18. using System.Collections.ObjectModel;
  19. using System.Text;
  20. namespace SharpMap.CoordinateSystems.Transformations
  21. {
  22. /// <summary>
  23. /// The GeographicTransform class is implemented on geographic transformation objects and
  24. /// implements datum transformations between geographic coordinate systems.
  25. /// </summary>
  26. public class GeographicTransform : MathTransform
  27. {
  28. internal GeographicTransform(IGeographicCoordinateSystem sourceGCS, IGeographicCoordinateSystem targetGCS)
  29. {
  30. _SourceGCS = sourceGCS;
  31. _TargetGCS = targetGCS;
  32. }
  33. #region IGeographicTransform Members
  34. private IGeographicCoordinateSystem _SourceGCS;
  35. /// <summary>
  36. /// Gets or sets the source geographic coordinate system for the transformation.
  37. /// </summary>
  38. public IGeographicCoordinateSystem SourceGCS
  39. {
  40. get { return _SourceGCS; }
  41. set { _SourceGCS = value; }
  42. }
  43. private IGeographicCoordinateSystem _TargetGCS;
  44. /// <summary>
  45. /// Gets or sets the target geographic coordinate system for the transformation.
  46. /// </summary>
  47. public IGeographicCoordinateSystem TargetGCS
  48. {
  49. get { return _TargetGCS; }
  50. set { _TargetGCS = value; }
  51. }
  52. /// <summary>
  53. /// Returns the Well-known text for this object
  54. /// as defined in the simple features specification. [NOT IMPLEMENTED].
  55. /// </summary>
  56. public override string WKT
  57. {
  58. get
  59. {
  60. throw new NotImplementedException();
  61. }
  62. }
  63. /// <summary>
  64. /// Gets an XML representation of this object [NOT IMPLEMENTED].
  65. /// </summary>
  66. public override string XML
  67. {
  68. get
  69. {
  70. throw new NotImplementedException();
  71. }
  72. }
  73. #endregion
  74. /// <summary>
  75. /// Creates the inverse transform of this object.
  76. /// </summary>
  77. /// <remarks>This method may fail if the transform is not one to one. However, all cartographic projections should succeed.</remarks>
  78. /// <returns></returns>
  79. public override IMathTransform Inverse()
  80. {
  81. throw new Exception("The method or operation is not implemented.");
  82. }
  83. /// <summary>
  84. /// Transforms a coordinate point. The passed parameter point should not be modified.
  85. /// </summary>
  86. /// <param name="point"></param>
  87. /// <returns></returns>
  88. public override SharpMap.Geometries.Point Transform(SharpMap.Geometries.Point point)
  89. {
  90. SharpMap.Geometries.Point pOut = point.Clone();
  91. pOut.X /= SourceGCS.AngularUnit.RadiansPerUnit;
  92. pOut.X -= SourceGCS.PrimeMeridian.Longitude / SourceGCS.PrimeMeridian.AngularUnit.RadiansPerUnit;
  93. pOut.X += TargetGCS.PrimeMeridian.Longitude / TargetGCS.PrimeMeridian.AngularUnit.RadiansPerUnit;
  94. pOut.X *= SourceGCS.AngularUnit.RadiansPerUnit;
  95. return pOut;
  96. }
  97. /// <summary>
  98. /// Transforms a list of coordinate point ordinal values.
  99. /// </summary>
  100. /// <remarks>
  101. /// This method is provided for efficiently transforming many points. The supplied array 
  102. /// of ordinal values will contain packed ordinal values. For example, if the source 
  103. /// dimension is 3, then the ordinals will be packed in this order (x0,y0,z0,x1,y1,z1 ...).
  104. /// The size of the passed array must be an integer multiple of DimSource. The returned 
  105. /// ordinal values are packed in a similar way. In some DCPs. the ordinals may be 
  106. /// transformed in-place, and the returned array may be the same as the passed array.
  107. /// So any client code should not attempt to reuse the passed ordinal values (although
  108. /// they can certainly reuse the passed array). If there is any problem then the server
  109. /// implementation will throw an exception. If this happens then the client should not
  110. /// make any assumptions about the state of the ordinal values.
  111. /// </remarks>
  112. /// <param name="points"></param>
  113. /// <returns></returns>
  114. public override Collection<SharpMap.Geometries.Point> TransformList(Collection<SharpMap.Geometries.Point> points)
  115. {
  116.             //Collection<SharpMap.Geometries.Point> trans = new Collection<SharpMap.Geometries.Point>(points.Count);
  117.             Collection<SharpMap.Geometries.Point> trans = new Collection<SharpMap.Geometries.Point>();
  118. foreach (SharpMap.Geometries.Point p in points)
  119. trans.Add(Transform(p));
  120. return trans;
  121. }
  122. /// <summary>
  123. /// Reverses the transformation
  124. /// </summary>
  125. public override void Invert()
  126. {
  127. throw new Exception("The method or operation is not implemented.");
  128. }
  129. }
  130. }