[讨论]：大地坐标与cad世界坐标的转化

panfanggui

已知cad两点坐标（x1,y1),(x2,y2)及对应的大地坐标（E1，N1），（E2，N2），怎样将任一大地坐标点（E，N）转化为cad对应的坐标及反之将cad内任一坐标（x,y）转化为对应的大地坐标。

wxsong

坐标的转换主要是坐标的平移、旋转、及缩放。
通过对两点长度比较可知缩放的比例。
通过两直线的夹角可计算出旋转的角度。
通过两点间的坐标差可得出平移的坐标。
得出这些转换矩阵后就可以进行转换和逆转换

panfanggui

谢谢回复，我已编程实现。源码如下：
void CTestDlg::OnButtonJs4()
{
        // TODO: Add your control notification handler code here
           double pi = 3.1415926535;
           double n,e,x,y,jj;
           float ax,ae;
// m_x1,m_y1 m_x2,m_y2,m_e1,m_n1,m_e2,m_n2 为两点的cad坐标及对应的大地坐标               
           n=(m_n2-m_n1);
           e=(m_e2-m_e1);
           x=(m_x2-m_x1);
           y=(m_y2-m_y1);

       m_k=sqrt((n*n+e*e)/(x*x+y*y));

           ax=atan2(y,x);
           ae=atan2(n,e);
           jj=(ax-ae);
           m_jj=(ax-ae)*180.0/pi;   

       m_ydx=m_x1-(cos(jj)*m_e1-sin(jj)*m_n1)/m_k;
       m_ydy=m_y1-(cos(jj)*m_n1+sin(jj)*m_e1)/m_k;
//以下计算坐标互换          
       m_x=m_ydx+(cos(jj)*m_e-sin(jj)*m_n)/m_k;
       m_y=m_ydy+(cos(jj)*m_n+sin(jj)*m_e)/m_k;

   m_e=m_k*(cos(jj)*(m_x-m_ydx)+sin(jj)*(m_y-m_ydy));
       m_n=m_k*(-sin(jj)*(m_x-m_ydx)+cos(jj)*(m_y-m_ydy));

           UpdateData(FALSE);

yiyuewanwan

我一直被这个问题捆饶，请问：具体怎么操作？不用list．

spring_1980

还有没简单的方法啊，我现在急需要解决啊，

ymd

大地坐标转换应使用高斯公式，不是简单的平移旋转！