• arcengine坐标转换(ZZ)


    现在好多人在使用ArcObject的时候都可能需要作经纬度和平面坐标的相互转换。由于经纬度是球面坐标,平面坐标是X-Y的笛卡尔坐标系统,所以这是一个看起来比较难的问题。

      
    好多人一上来就搬出地图学、地图投影学或者测绘学中的投影公式(如高斯投影 或 墨卡托投影),又是基准坐标又是角度的搞的人头都大了。

      
    实际上要想实现这个功能非常easy。 ARC Engine中的IPoint就可以进行投影和反投影运算了。

    投影过程(C#):

    ///   flatref
    投影的坐标系统,这里的54013是世界投影   ,世界投影所有经纬度都可以转换为平面坐标,但是由于投影面积大失真也会比较大(相当于把整个地球劈成一片片的,然后拉伸最后贴到平面上,失真当然大了)。当然也可以选择精度更高的平面如:esriSRProjCS_Beijing1954GK_23N    对应数值21483   仅仅把北京附近的地球平面拉伸铺在平面上,由于投影面积变小,所以投影经度提高。但是由于面积变小,所以有些经度纬度不能转换,比如所美国的精度纬度用北京投影就投不了。

           flatref = pfactory.CreateProjectedCoordinateSystem(54013);



    //
    没什么说的,标准大地经纬度,可以将X-Y逆投影为经度和纬度

       earthref = pfactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_NAD1983);


           /// 将经纬度点转换为平面坐标。
      
           private IPoint GetProject(double x, double y)
           {

             IPoint pt = new PointClass();
             
             pt.PutCoords(x, y);

             IGeometry geo = (IGeometry)pt;
             geo.SpatialReference = earthref;
             geo.Project(flatref);

             return pt;
           }



              /// 将平面坐标转换为经纬度。
          private IPoint GetGeo(double x, double y)
           {

             IPoint pt = new PointClass();

             pt.PutCoords(x, y);

             IGeometry geo = (IGeometry)pt;
             geo.SpatialReference = flatref;
             geo.Project(earthref);
             double xx = pt.X;
             return pt;
           }

    实际上IPoint的投影和任何地图都没什么大关系,完全可以不用地图,直接调用IPoint进行投影的转换

  • 相关阅读:
    rsync--数据镜像备份_转
    netcat
    tcpdump抓包
    find命令应用exec及xargs
    traceroute/tracert--获取网络路由路径
    TCP/IP各层协议数据格式
    (转)mq经验总结-转
    (转)WebSphere MQ基础命令
    MQ通道搭建以及连通性检查
    (转)java并发之CountDownLatch、Semaphore和CyclicBarrier
  • 原文地址:https://www.cnblogs.com/zhangjun1130/p/1750578.html
Copyright © 2020-2023  润新知