地理坐标到本地CAD坐标的坐标转换公式20140209.docVIP

地理全局坐标到本地局部CAD坐标的坐标转换公式 1、卫星位置用地球地心惯性坐标系ECI，由赤道面X轴（春分，由地心指向太阳），赤道面y轴（由地心指向春分向东90度方向），北极z轴（由地心指向地球北极），构成右手迪卡尔坐标系； 2、地面定位目标位置在导航电文中，使用WGS84地心地固坐标系ECEF，由赤道面X轴（由地心指向本初子午线，即0度经线方向），赤道面y轴（由地心指向东经90度方向），北极z轴（由地心指向地球北极），构成右手迪卡尔坐标系； 3、本地水平坐标系LTP（ENU），由本地切平面X，轴（由切点指向东，即东经方向），切平面y，轴（由切点指向北，即北纬方向），高度z，轴（由切点指向海拔高度方向），构成右手迪卡尔坐标系； 4、地固坐标系ECEF映射到经纬度坐标系的就是地理坐标系，由东经lon轴（由切点指向东经方向），北纬lat轴（由切点指向北纬方向）； 5、在由本地切平面不大时，切平面的EN坐标系和地理坐标系存在近似线性变换关系，使用线性坐标转换矩阵可满足工程应用精度； 6、在坐标转换精度要求不高时，且坐标转换范围不大时，可忽略地球扁率的影响，用球体近似表达地球表面； 7、EN（x，y）坐标系和地理坐标系（lon，lat）的近似坐标变换公式如下： ?x=r*?lon，其中，x的单位是m，lat的单位是弧度；r=R*COS(lat)，单位是m； ?Y=R*?lat，其中，y的单位是m，lat的单位是弧度，R是地球平均半径，6371000m； 所以有 {?x，?y}T= RC{?lon，?lat }T 其中C为二阶对角线矩阵： c11=cos(lat0) c22= 1 c12= c21=0 且：{?lon，?lat }T= {lon- lon0，lat- lat0 }T， 得到 {?x，?y}T = RC{ lon- lon0，lat- lat0 }T 8、EN（?x，?y）坐标系到CAD坐标系（X，Y）的坐标变换公式如下： {X，Y}T= K{?x，?y }T+S，这里{?x，?y}T = RC{ lon- lon0，lat- lat0 }T，S={ X0，Y0 }T 其中K为二阶旋转矩阵： k11= cos(α-θ) k12= -sin(α-θ) k21= sin(α-θ) k22= cos(α-θ) 通过基线测量得到两个旋转角： tanθ=dy/dx= dlat /(dlon* cos(lat0)) tanα=DY/DX 式中 dlon = lon1- lon0 dlat = lat1- lat0 DX = X 1- X 0 DY = Y 1- Y 0 见坐标变换示意图， 9、综合7和8两步变换，得到地理坐标到CAD坐标系（X，Y）的坐标变换公式如下： {X，Y}T= KRC{ lon- lon0，lat- lat0 }T +{ X0，Y0 }T， 地理坐标到CAD坐标变换示意图 10、考虑到CAD坐标系还需要归一化到切图坐标系下，对应的坐标变换公式如下： {x，y}T= (60/Lm)* ({ X，Y }T- { XC，YC } T)， 即：{x，y}T= (60/Lm)* (KRC{ lon- lon0，lat- lat0 }T +{ X0，Y0 }T- { XC，YC } T) 其中：（XC，YC）为 楼面积外接矩形中心坐标，Lm为楼面积外接矩形的最大边长，地理坐标单位是弧度。其中： Lm=Max{Xmax-Xmin，Ymax-Ymin}; XC=(Xmax-Xmin)/2+ Xmin; YC=(Ymax-Ymin)/2+ Ymin; 最终坐标转换计算公式为： x= (60/Lm)* (6371000* (cos(α-θ) * cos(lat0) * (lon- lon0)- sin(α-θ)* (lat- lat0) )+ X0- XC) y= (60/Lm)* (6371000* (sin(α-θ) * cos(lat0) * (lon- lon0)+ cos(α-θ)* (lat- lat0) )+ Y0- YC) 式中：Lm，XC，YC由楼面图数据计算得到，lon0，lat0，X0，Y0，lon1，lat1，X1，Y1由基线测量得到，α，θ由基线坐标数据计算得到。上式用于从GPS坐标向直角绘图坐标的转换。 每栋楼的位置特征由以下9个独立参量确定：lon0，lat0，X0，Y0，α，θ，Lm，XC，YC。 上式中的地理坐标单位为弧度，若地理坐标以度为单位，则： x= (60/Lm)* (6371000* (cos(α-θ) * cos(lat0) * (lon- lon0)*π/180- sin(α-θ)* (lat- lat0) *π/180 )+ X0- XC)； y= (60/L