第六章地球椭球与椭球计算理论.pptVIP

5、地心纬度坐标系 此两式指明了法线Pn在赤道两侧的长度。 令: pn=N（卯酉圈曲率半径） 由图看出： 与前式 相比得： 于是有 由图看出 上两式相比得： 显然有： 大地纬度B、归化纬度u、 地心纬度φ之间的关系 B和u之间的关系 归化纬度坐标系同子午平面坐标系的关系： 子午平面坐标系同 大地坐标系的关系： 8 、站心地平坐标系大地站心地平坐标系是以测站法线和子午线方向为依据的坐标系。在描述两点间关系时，为方便直观，一般采用站心坐标系。根据坐标表示方法，又可分为站心左手地平直角坐标系和站心地平极坐标系，见图 1）以测站P为原点，测站上P点的垂线(法线)为Z轴（U轴），指向天顶为正；子午线方向为x轴（N轴），指向参考短半轴，向北为正；y轴（E轴）与x，z轴平面垂直（向东为正）构成左手坐标系就称为垂线(或法线)站心直角坐标系。或称为站心天文坐标系。在站心直角坐标系下点的X（N） 、Y（E）、Z（U）坐标为该点在三个坐标轴上的投影长度。 站心极坐标系：以P点为中心的站心极坐标系定义如下： ( 1 ) X(N）PY（E） 平面为基准面； ( 2 ）极点位于P； ( 3 ）极轴为X（N） 轴。 Q点在站心极坐标系下的坐标用极距（ d 为 由极点到该点的距离）、方位角（ α为在 基准面上，以极点为顶点，由极轴顺时针 方向量测PQ在基准面上投影的角度）、 高度角（EL为极点与该点连线与基准面间 的夹角）表示。 空间任意一点Q相对于P的位 置可通过地面观测值——斜距d、 天文方位角α和天顶距z来确定 进行 GPS 观测时，常常采用 GPS 卫星相对于测站的高度角、方位角来描述其在空间中的方位。实际上，如果再加上测站到卫星的距离，就是一个完整的站心坐标。 卯酉圈曲率半径的特点:由公式N=r/cosB得： N与纬度B有关．它随B的增大而增大。卯酉圈曲率半径N与经度L无关，即同一平行圈上所有点的卯酉圈曲率半径N相同。 在极点子午圈曲率半径 和卯酉圈曲率半径一致 子午圈曲率半径M及卯酉圈曲率半径N是两个互相垂直的法截弧的曲率半径，统称为主曲率半径。 常引用的两个符号：（1）=ρ″/M （2）=ρ″/N 主曲率半径的计算 以上讨论的子午圈曲率半径M及卯酉圈曲率半径N，是两个互相垂直的法截弧的曲率半径，这在微分几何中统称为主曲率半径。 （1）主曲率半径是第一纬度函数的表达式 任意法截弧的曲率半径 由欧拉公式知： 在实际应用中，用平均曲率半径R代替N,得： 任意法截弧的曲率半径的变化规律: ＲＡ不仅与点的纬度B有关，而且还与过该点的法截弧的方位角A有关。 当Ａ＝０°或180°时，RA值为最小，公式变为计算子午圈曲率半径的，即Ｒ０＝Ｍ； 当A＝90°或270°时，RA值最大，公式变为计算卯酉圈曲率半径，即Ｒ９０＝Ｎ。主曲率半径M及N分别是ＲＡ的极小值和极大值。 当A由0°→90°时，ＲＡ之值由Ｍ→Ｎ，当A由90°→180°时，ＲＡ值由N→Ｍ，可见ＲＡ值的变化是以90°为周期且与子午圈和卯酉圈对称的。 椭球面上任意一点的平均曲率半径 R 等于该点主曲率半径（子午圈曲率半径M和卯酉圈曲率半径N）的几何平均值。 M，N，R的关系 一般情况下 NRM 在极点上都等于极曲率半径ｃ 对于克拉索夫斯基椭球 如果以B＝90°代入，则得子午椭圆在一个象限内的弧长约为10 002 137m。旋转椭球的子午圈的整个弧长约为40 008 549.995m。即一象限子午线弧长约为10 000km，地球周长约为40 000km。 为求子午线上两个纬度B１及Ｂ２间的弧长，只需按(4-101)式分别算出相应的X１及X２，而后取差：ΔＸ＝Ｘ２－Ｘ１，该ΔＸ即为所求的弧长。 当弧长甚短(例如X≤40km，计算精度到0.001m)，可视子午弧为圆弧，而圆的半径为该圆弧上平均纬度点的子午圈的曲率半径Mｍ X=Mm.（B2-B1)″/ρ″=ΔB″/（1）m 2、由子午弧长求大地纬度 底点纬度计算 在高斯投影反算时，已知高斯平面直角坐标（X，Y）反求其大地坐标（L，B）。首先X当作中央子午线上弧长，反求其纬度，此时的纬度称为底点纬度或垂直纬度。计算底点纬度的公式可以采用迭代解法和直接解法。 （1）迭代解法: 重复迭代直至 　　　　　为止。a0为子午圈曲率半径展开的首项系数。 例如： 4.5 大地线 两点间的最短距离，在平面上是两点间的直线，在