使用ProjNet.Silverlight进行ArcGis座标转换(有解决出错问题方法)

使用ProjNet.Silverlight可以将ArcGis中的从标由大地座标转换为经纬座标,代码如下

 

        public List<MapPoint> ConvertTo(List<MapPoint> mps)
        {

            CoordinateTransformationFactory ctfac = new CoordinateTransformationFactory();
            
            Point p = new Point(108.94, 34.28);
            string bj1954Lcc = "PROJCS[\"China_Lambert_Conformal_Conic\",GEOGCS[\"GCS_Beijing_1954\",DATUM[\"D_Beijing_1954\",SPHEROID[\"Krasovsky_1940\",6378245.0,298.3]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Conformal_Conic\"],PARAMETER[\"False_Easting\",0.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",105.0],PARAMETER[\"Standard_Parallel_1\",30.0],PARAMETER[\"Standard_Parallel_2\",62.0],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]";
            

            IProjectedCoordinateSystem fromCS = CoordinateSystemWktReader.Parse(bj1954Lcc) as IProjectedCoordinateSystem;
            GeographicCoordinateSystem toCS = (GeographicCoordinateSystem)fromCS.GeographicCoordinateSystem;
            List<MapPoint> list = new List<MapPoint>();
            foreach (var mp in mps)
            {
                var p1 = new Point(mp.X, mp.Y);
                var p2 = ctfac.CreateFromCoordinateSystems(fromCS, toCS).MathTransform.Transform(p1);
                list.Add(new MapPoint(p2.X, p2.Y));
            }
            return list;
        }


下载ProjNet.Silverlight后在Silverlight中使用是会出一个'p' is not recognized.的错误,找了半天没有解决,最后反编译代码然后跟踪发现原来有一个字符处理时使用的是Encoding是Unicode,在中文的C#中字符串的默认编码是UTF-8,将代码中的Unicode改为UTF-8即可!

需要改的代码在:ProjNet.Converters.WellKnownText.IO.StreamTokenizer类中的

 private TokenType NextTokenAny()中!

  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
如果你不想使用ProjNet库,你可以手动实现高斯投影3°和6°带之间的转换。这里有一个简单的示例代码,使用数学公式来进行转换: ```csharp public class GaussProjectionConverter { public double[] ConvertFrom3DegreeTo6Degree(double longitude, double latitude) { double a = 6378137; // WGS84椭球体长半轴 double f = 1 / 298.257223563; // WGS84椭球体扁率 double L0 = Math.Floor(longitude / 3) * 3; // 3°带的中央经线 double l = longitude - L0; // 经度差 double L = latitude; // 纬度不变 double e2 = 2 * f - f * f; double e12 = e2 / (1 - e2); double n = a / Math.Sqrt(1 - e2 * Math.Sin(L) * Math.Sin(L)); double X = n * Math.Cos(L) * l; double Y = n * (1 - e2) * Math.Sin(L); return new double[] { L0 + X / (6367000 * Math.Cos(L)), Y / 6367000 }; } public double[] ConvertFrom6DegreeTo3Degree(double longitude, double latitude) { double a = 6378137; // WGS84椭球体长半轴 double f = 1 / 298.257223563; // WGS84椭球体扁率 double L0 = Math.Floor(longitude / 6) * 6 + 3; // 6°带的中央经线 double l = longitude - L0; // 经度差 double L = latitude; // 纬度不变 double e2 = 2 * f - f * f; double e12 = e2 / (1 - e2); double n = a / Math.Sqrt(1 - e2 * Math.Sin(L) * Math.Sin(L)); double X = n * Math.Cos(L) * l; double Y = n * (1 - e2) * Math.Sin(L); return new double[] { L0 + X / (6367000 * Math.Cos(L)), Y / 6367000 }; } } ``` 这个示例代码中的`ConvertFrom3DegreeTo6Degree`方法将3°带的经纬度转换为6°带的经纬度,而`ConvertFrom6DegreeTo3Degree`方法将6°带的经纬度转换为3°带的经纬度。这两个方法都接受经度和纬度作为输入,并返回转换后的经纬度。 请注意,这只是一个简单的示例代码,它使用了一些近似值和简化的计算公式。在实际应用中,你可能需要根据你的需求和具体的高斯投影参数进行适当的调整。同时,确保你已经了解高斯投影的相关知识,并在使用之前,仔细检查和验证你的代码。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值