地理坐标到本地cad坐标的坐标转换公式20140 209.pdfVIP

地理全局坐标到本地局部CAD坐标的坐标转换公式 1、卫星位置用地球地心惯性坐标系ECI，由赤道面X 轴（春分，由地心指向太阳），赤道面 y轴（由地心指向春分向东90度方向），北极z轴（由地心指向地球北极），构成右手迪卡 尔坐标系； 2、地面定位目标位置在导航电文中，使用WGS84 地心地固坐标系ECEF，由赤道面X 轴 （由 地心指向本初子午线， 0 度经线方向），赤道面y 轴（由地心指向东经90 度方向），北极 z 轴（由地心指向地球北极），构成右手迪卡尔坐标系； 3、本地水平坐标系 LT （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； 所以有 T T {∆x ，∆y} = RC{∆lon，∆lat } 其中C 为二阶对角线矩阵： c =cos(lat ) 11 0 c = 1 22 c = c =0 12 21 T T 且：{∆lon，∆lat } = {lon- lon ，lat- lat } ， 0 0 得到 T T {∆x ，∆y} = RC{ lon- lon ，lat- lat } 0 0 8、EN （∆x ，∆y ）坐标系到CAD 坐标系（X ，Y ）的坐标变换公式如下： T T T T T {X，Y} = K{∆x ，∆y } +S，这里{∆x ，∆y} = RC{ lon- lon ，lat- lat } ，S={ X ，Y } 0 0 0 0 其中K 为二阶旋转矩阵： k = cos(α-θ) 11 k = -sin(α-θ) 12 k = sin(α-θ) 21 k = cos(α-θ) 22 通过基线测量得到两个旋转角： tanθ=dy/dx= dlat /(dlon* cos(lat )) 0 tanα=DY/DX 式中 dlon = lon - lon 1 0 dlat = lat - lat 1 0 1 DX = X - X 1 0 DY = - 1 0 见坐标变换示意图， 9、综合7 和8 两步变换，得到地理坐标到CAD 坐标系（X ， ）的坐标变换公式如下：