第六章地球橢球與橢球計算理論

1、第六章 地球橢球與橢球計算理論,本章提要 6.1 地球橢球的基本幾何參數(shù)及其相互關(guān)系 6.2 橢球面上的常用坐標(biāo)系及其相互關(guān)系 6.3 幾種主要的橢球公式 6.4 將地面觀測值歸算至橢球面 習(xí)題,本章提要,本章講述地球橢球與參考橢球的概念，進而介紹橢球的基本幾何參數(shù)，基本坐標(biāo)系及其相互關(guān)系。同時，講述橢球面同地面之間的關(guān)系，如何將地面觀測元素（水平方向及斜距等）歸算至橢球面上。在對本章的學(xué)習(xí)中，要建立起空間的概念，只有建立了地球橢球的這些基本空間概念后，才能更好地學(xué)習(xí)控制測量的內(nèi)業(yè)數(shù)據(jù)處理等相關(guān)知識。,1地球橢球的定義及其幾何意義； 2常用測量坐標(biāo)系統(tǒng)的建立及其在控制測量中的應(yīng)用； 3各種測量

2、坐標(biāo)系統(tǒng)之間的相互轉(zhuǎn)換； 4橢球面上幾種曲率、弧長、大地線的計算； 5地面測量值（水平方向和邊長）歸算到橢球面的方法。,知識點及學(xué)習(xí)要求,在對本章的學(xué)習(xí)中，有大量的公式推導(dǎo)與應(yīng)用。 各種常用測量坐標(biāo)系統(tǒng)的建立與相互轉(zhuǎn)換； 幾種常用的橢球計算公式； 地面觀測值歸算到橢球面的方法與計算、大地主題解算 。,本章重點：地球橢球幾何性質(zhì)、地面觀測值歸算 本章難點：,6.1 地球橢球的基本幾何參數(shù)及其相互關(guān)系,一. 地球橢球的基本幾何參數(shù),地球橢球：在控制測量中，用來代表地球的橢球，它是地球的數(shù)學(xué)模型。,參考橢球：具有一定幾何參數(shù)、定位及定向的用以代表某一地區(qū)大地水準面的地球橢球。地面上一切觀測元素都應(yīng)歸

3、算到參考橢球面上，并在這個面上進行計算。參考橢球面是大地測量計算的基準面，同時又是研究地球形狀和地圖投影的參考面。,地球橢球的幾何定義：O是橢球中心，為旋轉(zhuǎn)軸，a 為長半軸，b 為短半軸。,子午圈(經(jīng)圈，或子午橢圓）：包含旋轉(zhuǎn)軸的平面與橢球面相截所得的橢圓。如NKAS,平行圈（或緯圈）：垂直于旋轉(zhuǎn)軸的平面與橢球面相截所得的圓。如QKQ 赤道：通過橢球中心的平行圈 。如EAE,地球橢球是選擇的旋轉(zhuǎn)橢球,旋轉(zhuǎn)橢球的形狀和大小常用子午橢圓的五個基本幾何參數(shù)(或稱元素): 地球橢球的五個基本幾何參數(shù)： 橢圓的長半軸 橢圓的短半軸b 橢圓的扁率 橢圓的第一偏心率,橢圓的第二偏心率,其中a、b 稱為長度元

4、素；扁率 反映了橢球體的扁平程度。偏心率 和 是子午橢圓的焦點離開中心的距離與橢圓半徑之比， 它們也反映橢球體的扁平程度，偏心率愈大，橢球愈扁。,二、通常用a ,表示橢球的形狀和大小。,我國建立1954年北京坐標(biāo)系應(yīng)用的是克拉索夫斯基橢球；建立1980年國家大地坐標(biāo)系應(yīng)用的是1975年國際橢球；而全球定位系統(tǒng)(GPS)應(yīng)用的是WGS-84系橢球參數(shù)。,幾種常見的橢球體參數(shù)值,三、相互關(guān)系 1e與 e的關(guān)系,2e與其它關(guān)系,四、引用符號及其相互關(guān)系 1引用符號,2W與V關(guān)系,兩個常用的輔助函數(shù)，W第一基本緯度函數(shù)，V第二基本緯度函數(shù)：,C幾何意義：極點處的幾何曲率半徑。,3. 地球橢球參數(shù)間的相

5、互關(guān)系,其他元素之間的關(guān)系式如下：,式中，W 第一基本緯度函數(shù)，V 第二基本緯度函數(shù)。,6.2 橢球面上的常用坐標(biāo)系及其相互關(guān)系,1 .大地坐標(biāo)系,p 點的子午面NPS 與起始子午面 NGS 所構(gòu)成的二面角L，叫做p 點的大地經(jīng)度，由起始子午面起算，向東為正，叫東經(jīng)（0180），向西為負，叫西經(jīng)（0180）。P 點的法線 與赤道面的夾角B，叫做P點的大地緯度。由赤道面起算，向北為正，叫北緯（090）；向南為負，叫南緯(090)。從地面點P沿橢球法線到橢球面的距離叫大地高。大地坐標(biāo)系是用大地經(jīng)度L、大地緯度B和大地高H表示地面點位的。如果點不在橢球面上，表示點的位置除B,L外，還要附加另一參數(shù)大

6、地高H，它同正常高及正高有如下關(guān)系,2. 空間直角坐標(biāo)系,地心 坐標(biāo)系，原點位于總地球橢球(或參考橢球)質(zhì)心；Z軸與地球平均自轉(zhuǎn)軸相重合，亦即指向某一時刻的平均北極點；X軸指向平均自轉(zhuǎn)軸與平均格林尼治天文臺所決定的子午面與赤道面的交點G；Y軸與此平面垂直，且指向東為正。 地心空間直角系與參心空間直角坐標(biāo)系之分。,3 . 子午面直角坐標(biāo)系,設(shè)點 p 的大地經(jīng)度L為，在過p點的子午面上，以子午圈橢圓中心為原點，建立x,y 平面直角坐標(biāo)系。在該坐標(biāo)系中，p 點的位置用L,x,y 表示。,4 . 大地極坐標(biāo)系,M 為橢球體面上任意一點，MN 為過M 點的子午線，S 為連結(jié)MP的大地線長，A 為大地線在

7、M 點的方位角。以M 為極點，MN 為極軸，S 為極半徑，A為極角，這樣就構(gòu)成大地極坐標(biāo)系。在該坐標(biāo)系中p 點的位置用S,A 表示。,橢球面上點的極坐標(biāo)（S,A）與大地坐標(biāo)(L,B）可以互相換算，這種換算叫做大地主題解算。,5、地心緯度坐標(biāo)系,橢球面上P點的大地經(jīng)度L，在此子午面上以橢圓中心O為原點建立地心緯度坐標(biāo)系。連接OP，則 稱為地心緯度，而OP 稱為P點向徑，在此坐標(biāo)系中，點的位置為：,6、 歸化緯度坐標(biāo)系,設(shè)橢球面上P點的大地經(jīng)度為L，在此子午面上以橢圓中心O為圓心，以橢球長半徑a為半徑作輔助圓，延長P2P與輔助圓相交P1點，則OP1與x軸夾角稱為P點的歸化緯度，用u表示，在此歸化緯

8、度坐標(biāo)系中，P點位置用L，u 表示。,7. 各坐標(biāo)系間的關(guān)系,橢球面上的點位可在各種坐標(biāo)系中表示，由于所用坐標(biāo)系不同，表現(xiàn)出來的坐標(biāo)值也不同。,1）子午面直角坐標(biāo)系同大地坐標(biāo)系的關(guān)系,過p 點作法線 ，它與x 軸之夾角為B，過點P作子午圈的切線TP，它與x 軸的夾角為（90+B）。子午面直角坐標(biāo)x,y 同大地緯度B 的關(guān)系式如下：,以下：推導(dǎo)子午平面坐標(biāo)系同大地坐標(biāo)系的關(guān)系,此兩式指明了法線Pn在赤道兩側(cè)的長度。,令: pn=N（卯酉圈曲率半徑） 由圖看出： 與前式 相比得： 于是有 由圖看出 上兩式相比得： 顯然有：,上兩式即為子午面直角坐標(biāo)x，y同大地緯度B的關(guān)系式。,2）子午平面直角坐標(biāo)

9、系同歸化緯度坐標(biāo)系的關(guān)系 P(x，y)，OP1=a 由圖可知：x=OP2=OP1cos=acos 代入公式： 得：y=bsin,3）空間直角坐標(biāo)系同子午面直角坐標(biāo)系的關(guān)系,空間直角坐標(biāo)系中的P2P相當(dāng)于子午平面直角坐標(biāo)系中的y，前者的OP2相當(dāng)于后者的x，并且二者的經(jīng)度L相同。 OP2=x 二面角P2OP1=L,4）空間直角坐標(biāo)系同大地坐標(biāo)系的關(guān)系,同一地面點在地球空間直角坐標(biāo)系中的坐標(biāo)和在大地坐標(biāo)系中的坐標(biāo)可用如下兩組公式轉(zhuǎn)換,式中：e子午橢圓第一偏心率，可由長短半徑按式 算得。,N法線長度，可由式 算得。,以下推導(dǎo)空間直角坐標(biāo)系同大地坐標(biāo)系的關(guān)系,推導(dǎo)空間直角坐標(biāo)系同大地坐標(biāo)系的關(guān)系 由大

10、地坐標(biāo)計算空間直角坐標(biāo)：,如果P點在橢球面上： （1）子午面直角坐標(biāo)同大地坐標(biāo)關(guān)系， （2）空間直角坐標(biāo)同子午面直角坐標(biāo)關(guān)系,如果P點不在橢球面上，設(shè)大地高為H，P點在橢球面上投影為P0,由空間直角坐標(biāo)計算大地坐標(biāo),大地經(jīng)度,大地緯度,大地高,由空間直角坐標(biāo)計算相應(yīng)大地坐標(biāo),大地緯度B、歸化緯度u、 地心緯度之間的關(guān)系 B和u之間的關(guān)系 歸化緯度坐標(biāo)系同子午平面坐標(biāo)系的關(guān)系： 子午平面坐標(biāo)系同 大地坐標(biāo)系的關(guān)系：,U、之間的關(guān)系 地心緯度同子午平面坐標(biāo)系關(guān)系 歸化緯度U同子午平面坐標(biāo)系：,、之間的關(guān)系,大地緯度、地心緯度、歸化緯度之間的差異很小，經(jīng)過計算，當(dāng)B=45時,8 、站心地平坐標(biāo)系大地

11、站心地平坐標(biāo)系是以測站法線和子午線方向為依據(jù)的坐標(biāo)系。在描述兩點間關(guān)系時，為方便直觀，一般采用站心坐標(biāo)系。根據(jù)坐標(biāo)表示方法，又可分為站心左手地平直角坐標(biāo)系和站心地平極坐標(biāo)系，見圖 1）以測站P為原點，測站上P點的垂線(法線)為Z軸（U軸），指向天頂為正；子午線方向為x軸（N軸），指向參考短半軸，向北為正；y軸（E軸）與x，z軸平面垂直（向東為正）構(gòu)成左手坐標(biāo)系就稱為垂線(或法線)站心直角坐標(biāo)系?；蚍Q為站心天文坐標(biāo)系。在站心直角坐標(biāo)系下點的X（N） 、Y（E）、Z（U）坐標(biāo)為該點在三個坐標(biāo)軸上的投影長度。,站心極坐標(biāo)系：以P點為中心的站心極坐標(biāo)系定義如下： ( 1 ) X(N）PY（E） 平面為

12、基準面； ( 2 ）極點位于P； ( 3 ）極軸為X（N） 軸。 Q點在站心極坐標(biāo)系下的坐標(biāo)用極距（ d 為 由極點到該點的距離）、方位角（ 為在 基準面上，以極點為頂點，由極軸順時針 方向量測PQ在基準面上投影的角度）、 高度角（EL為極點與該點連線與基準面間 的夾角）表示。 空間任意一點Q相對于P的位 置可通過地面觀測值斜距d、 天文方位角和天頂距z來確定 進行 GPS 觀測時，常常采用 GPS 衛(wèi)星相對于測站的高度角、方位角來描述其在空間中的方位。實際上，如果再加上測站到衛(wèi)星的距離，就是一個完整的站心坐標(biāo)。,7、大地測量常用坐標(biāo)系比較,6.3 橢球面上幾種曲率半徑（橢球數(shù)學(xué)性質(zhì)）,過橢球

13、面上任意一點可作一條垂直于橢球面的法線，包含這條法線的平面叫做法截面，法截面同橢球面交線叫法截線（或法截弧）。包含橢球面一點的法線，可作無數(shù)多個法截面，相應(yīng)有無數(shù)多個法截線。 橢球面上的法截線曲率半徑不同于球面上的法截線曲率半徑都等于圓球的半徑，而是不同方向的法截弧的曲率半徑都不相同。,1. 子午圈曲率半徑,子午橢圓的一部分上取一微分弧長DKdS，相應(yīng)地有坐標(biāo)增量dx，點n是微分弧dS的曲率中心，于是線段Dn及Kn便是子午圈曲率半徑。 由任意平面曲線的曲率半徑的定義公式，易知,任意平面曲線的曲率半徑的定義公式為： 從微分三角形DKE可求得： 聯(lián)立上兩式得： 子午面直角坐標(biāo)同大地坐標(biāo)關(guān)系： x坐

14、標(biāo)對大地緯度B取導(dǎo)數(shù)： 又由第一緯度函數(shù)得： 則有： 子午圈曲率半徑公式為： 或 M與緯度B有關(guān)它隨B的增大而增大。,變化規(guī)律如下表所示： 表中極曲率半徑c的幾何意義是橢球體在極點處的曲率半徑。,2.卯酉圈曲率半徑,過橢球面上一點的法線，可作無限個法截面，其中一個與該點子午面相垂直的法截面同橢球面相截形成的閉合的圈稱為卯酉圈。在圖中PEE 即為過點的卯酉圈。卯酉圈的曲率半徑用表示。 為了推導(dǎo)的表達計算式，過點作以O(shè)為中心的平行圈的切線，該切線位于垂直于子午面的平行圈平面內(nèi)。因卯酉圈也垂直于子午面，故也是卯酉圈在點處的切線。即垂直于Pn 。所以是平行圈及卯酉圈PEE在點處的公切線。,麥尼爾定理:

15、 假設(shè)通過曲面上一點引兩條截弧，一為法截弧，一為斜截弧，且在該點上這兩條截弧具有公共切線，這時斜截弧在該點處的曲率半徑等于法截弧的曲率半徑乘以兩截弧平面夾角的余弦。 如圖根據(jù)直角三角形有平行圈半徑r： 根據(jù)子午面直角坐標(biāo)同大地坐標(biāo)關(guān)系： 卯酉圈曲率半徑： 根據(jù)（4-8）式得： 由圖可看出： 也就是說 卯酉圈曲率半徑恰好等于法線介于橢球面和短軸之間的長度，亦即卯酉圈的曲率中心位于橢球的旋轉(zhuǎn)軸上。,卯酉圈曲率半徑的特點:由公式N=r/cosB得： N與緯度B有關(guān)它隨B的增大而增大。卯酉圈曲率半徑N與經(jīng)度L無關(guān)，即同一平行圈上所有點的卯酉圈曲率半徑N相同。 在極點子午圈曲率半徑 和卯酉圈曲率半徑一致

16、 子午圈曲率半徑M及卯酉圈曲率半徑N是兩個互相垂直的法截弧的曲率半徑，統(tǒng)稱為主曲率半徑。 常引用的兩個符號：（1）=/M （2）=/N,主曲率半徑的計算 以上討論的子午圈曲率半徑M及卯酉圈曲率半徑N，是兩個互相垂直的法截弧的曲率半徑，這在微分幾何中統(tǒng)稱為主曲率半徑。 （1）主曲率半徑是第一緯度函數(shù)的表達式,橢球面上幾種曲率半徑,橢球面上幾種曲率半徑,橢球面上幾種曲率半徑,（2）主曲率半徑是第二緯度函數(shù)的表達式,3.任意法截弧的曲率半徑,子午法截弧是南北方向，其方位角為0或180。卯酉法截弧是東西方向，其方位角為90或270?，F(xiàn)在來討論方位角為的任意法截弧的曲率半徑RA的計算公式。 由歐拉公式知

17、，任意方向方位角為的法截弧的曲率半徑的計算公式如下： 以下是公式的推導(dǎo)：,任意法截弧的曲率半徑 由歐拉公式知： 在實際應(yīng)用中，用平均曲率半徑R代替N,得：,橢球面上幾種曲率半徑,任意法截弧的曲率半徑的變化規(guī)律: 不僅與點的緯度B有關(guān)，而且還與過該點的法截弧的方位角A有關(guān)。 當(dāng)或180時，RA值為最小，公式變?yōu)橛嬎阕游缛η拾霃降模矗?當(dāng)A90或270時，RA值最大，公式變?yōu)橛嬎忝先η拾霃?，即。主曲率半徑M及N分別是的極小值和極大值。 當(dāng)A由090時，之值由，當(dāng)A由90180時，值由N，可見值的變化是以90為周期且與子午圈和卯酉圈對稱的。,橢球面上幾種曲率半徑,橢球面上任意一點的平均曲率半

18、徑 R 等于該點主曲率半徑（子午圈曲率半徑M和卯酉圈曲率半徑N）的幾何平均值。,4.平均曲率半徑R是指經(jīng)過曲面任意一點所有可能方向上的法截線半徑RA的算術(shù)平均值。 平均曲率半徑在實際際工程應(yīng)用中，根據(jù)測量工作的精度要求，在一定范圍內(nèi)，把橢球面當(dāng)成具有適當(dāng)半徑的球面。取過地面某點的所有方向RA的平均值來作為這個球體的半徑是合適的。這個球面的半徑平均曲率半徑R：,M，N，R的關(guān)系 一般情況下 NRM 在極點上都等于極曲率半徑,橢球面上幾種曲率半徑,對于克拉索夫斯基橢球,橢球面上幾種曲率半徑,從表看出：同一緯度時， NRM； 不同緯度時，隨著緯度B的增大，所有的曲率半徑都是增大的，且在極點上都等于極

19、曲率半徑,4.4 橢球面上的弧長計算 1. 子午線弧長計算公式 子午橢圓的一半，它的端點與極點相重合；而赤道又把子午線分成對稱的兩部分。因此，只需計算從赤道開始到已知緯度B間的子午線弧長的計算公式就可以了。 如下圖所示，取子午線上某微分弧PP=dx，令點緯度為，P點緯度為B+dB，點的子午圈曲率半徑為，于是有： 從赤道開始到任意緯度的平行圈之間的弧長可由下列積分求出： 式中M可用下式表達（第一緯度函數(shù)表達）：,其中：,經(jīng)積分，進行整理后得子午線弧長計算式（1）：,為求子午線上兩個緯度B1及B2間的弧長，只需按上式分別算出相應(yīng)的X1及X2，而后取差：，該即為所求的弧長。,將克拉索夫斯基橢球元素代

20、入,得 將1975年國際橢球元素代入,得 若將（4-70）公式中的冪函數(shù)展開為正弦n次冪和余弦乘積的形式，進行整理后得子午線弧長計算式（2）（此公式適合計算機計算）： 代入克拉索夫斯基橢球參數(shù)得子午線弧長計算公式： 代入1975年國際橢球參數(shù)得子午線弧長計算公式： 如果主曲率半徑用第二緯度函數(shù)的表達式，則可得子午線弧長計算公式（3）： 或,如果以B90代入，則得子午橢圓在一個象限內(nèi)的弧長約為10 002 137m。旋轉(zhuǎn)橢球的子午圈的整個弧長約為40 008 549.995m。即一象限子午線弧長約為10 000km，地球周長約為40 000km。 為求子午線上兩個緯度B及間的弧長，只需按(4-1

21、01)式分別算出相應(yīng)的X及X，而后取差：，該即為所求的弧長。 當(dāng)弧長甚短(例如X40km，計算精度到0.001m)，可視子午弧為圓弧，而圓的半徑為該圓弧上平均緯度點的子午圈的曲率半徑M X=Mm.（B2-B1)/=B/（1）m,2、由子午弧長求大地緯度 底點緯度計算 在高斯投影反算時，已知高斯平面直角坐標(biāo)（X，Y）反求其大地坐標(biāo)（L，B）。首先X當(dāng)作中央子午線上弧長，反求其緯度，此時的緯度稱為底點緯度或垂直緯度。計算底點緯度的公式可以采用迭代解法和直接解法。 （1）迭代解法: 重復(fù)迭代直至 為止。a0為子午圈曲率半徑展開的首項系數(shù)。 例如：,在克拉索夫斯基橢球上計算時，迭代開始時設(shè),以后每次迭

22、代按下式計算：,重復(fù)迭代直至 為止。,在1975年國際橢球上計算時，也有類似公式。,（2）直接解法,1975年國際橢球：,克拉索夫斯基橢球：,3、平行圈弧長公式 旋轉(zhuǎn)橢球體的平行圈是一個圓，其短半軸 就是圓上任意一點的子午面直角坐標(biāo)。 如果平行圈上有兩點，它們的經(jīng)度差 平行圈弧長公式： 平行圈弧長隨緯度變化的微分公式可近似地寫為 由于 于是 式中,4、子午線弧長和平行圈弧長變化的比較 從表中可以看出， 單位緯差的子午線弧長隨緯度升高而緩慢地增長； 而單位經(jīng)差的平行圈弧長則隨緯度升高而急劇縮短。 同時還可以看出，1的子午弧長約為110km，1約為1.8km，1約為30m。 而平行圈弧長，僅在赤道

23、附近才與子午弧長大體相當(dāng)，隨緯度的升高它們的差值愈來愈大。 5.橢球面梯形圖幅面積的計算,4.5 大地線 兩點間的最短距離，在平面上是兩點間的直線，在球面上是兩點間的大圓弧，那么在橢球面上又是怎樣的一條線呢? 它應(yīng)是大地線。 一、相對法截線 在橢球面上取 不在同一子午 面和平行圈上 的點A和B；過 A、B兩點的法 線分別與短軸 PP1相交于N1和N2。 通過AN1含有B的 法面AN1B與過BN2含有A的法面BN2A不可能重合而只能相交，其交線就是A、B間所連的直線。兩個法面所截得的法截線AaB和BbA也不重合而形成一狹小的二面角。這兩條法截線稱為相對法截線。現(xiàn)在證明na和nb將不重合。,設(shè)在橢

24、球上任取不在同一子午面和平行圈上的兩點A和B，緯度分別為B1和B2，且二者不等，過A、B兩點分別做法線與短軸交于na和nb兩點，與赤道面分別交于Q1和Q2。,如圖,又,則,正、反法截線： A點照準B，照準面同橢球面的交線AaB，叫做A點的正法截線，或者B點的反法截線；同樣，B點照準A，照準面同橢球面的交線BbA，叫做B點的正法截線，或者A點的反法截線。 把AaB和BbA叫做A、B兩點的相對法截線。 由公式 知，當(dāng)B2B1時，卯酉圈曲率半徑N隨緯度B的增大而增大，OnbOna。這就是說：某點的緯度愈高，其法線與短軸的交點愈低，AB方向在不同的象限時，正反法截線關(guān)系如圖：,由緯度B低的點照準緯度高

25、的點，法截線南偏； 由緯度B高的點照準緯度低的點，法截線北偏； 在北半球，如A點位于B點以南，則A點的正法 截線AaB在反法截線BbA的南面，如上圖。當(dāng)A、B兩點在同一子午圈或同一平行圈上時，相對法截線是重合的，夾角為0。 相對法截線通常是不重合的，兩者所夾的小角有下列數(shù)值： 當(dāng) S=15km時，=0.001 S=20km時，=0.002 S=30km時，=0.007 S=50km時，=0.01 由上可知: 相對法截線所夾的角是很小很小的。當(dāng)S=50km時，Bm=45，A=45時，兩法截線分開的最大距離為0.0008m。 從上面的數(shù)值看，當(dāng)距離不超過30公里時，可以不考慮所引起的問題，而認為相

26、對法截線是重合的。,當(dāng)兩點間距離較長時，相對法截線不重合會帶來如下的問題：,如上圖，在A、B、C三個點上，由兩條正法截線表示的角度，并不能構(gòu)成一完整的三角形，這就造成幾何圖形的破裂。因此需要在相對法截線間尋求一條合適的曲線，以代替相對法截線。經(jīng)理論證明，這一合適的曲線就是介于相對法截線之間的大地線。,二、大地線的定義和性質(zhì) 在橢球面上兩點間最短的線是大地線。在微分幾何中，大地線（又稱測地線）另有這樣的定義：“大地線上每點的密切面（無限接近的三個點構(gòu)成的平面）都包含該點的曲面法線”，亦即“大地線上各點的主法線與該點的曲面法線重合”。因曲面法線互不相交，故大地線是一條空間曲面曲線。 假設(shè)橢球面沒有

27、任何磨擦力，在橢球面上兩點間緊拉一條細線，則此細線的位置就和大地線一致。 大地線的形狀如圖所示，大地線是一條曲線，其長度可以用下式計算：,大地線與法截線的長度之差甚微，實際上可以不 必考慮兩者之差。大地經(jīng)緯度決定了地面點在橢 球面上的絕對位置，大地方位角則決定了橢球面 上兩點間大地線的方向，大地線和由它們所構(gòu)成 的球面角是組成橢球面上大地控制網(wǎng)的基本元素。,大地線的性質(zhì): 大地線是介于相對法截線之間的一“S” 形曲線。兩端與正法截線的夾角為。不在同一子午圈或同一平行圈上的兩點的正反法裁線是不重合的，它們之間的夾角；大地線是兩點間惟一最短線，而且位于相對法截線之間，并靠近正法截線，它與正法截線間

28、的夾角 在橢球面上進行測量計算時，應(yīng)當(dāng)以兩點間的大地線為依據(jù)。在地面上測得的方向、距離等，應(yīng)當(dāng)歸算成相應(yīng)大地線的方向、距離。 實際計算中：在一等三角測量中，數(shù)值可達干 分之一二秒，可見在一等或相當(dāng)于一等三角測量 精度的工程三角測量中是不容忽略的。 大地線與法截線長度之差只有百萬分之一毫 米，所以在實際計算中，這種長度差異總是可忽 略不計的。 長度差異可忽略,方向差異需改化。,大地線的微分方程和克萊勞方程 設(shè)P為大地線上任意一點，其經(jīng)度為，緯度為，大地線方位角為。當(dāng)大地線增加d到P1點時，則 上述各量相應(yīng)變化dL，dB及dA。 所謂大地線微分方程，即表示 dL、dB和dA與dS的關(guān)系。 dS在子

29、午圈上的分量 dS在平行圈上的分量,大地線的微分方程,以上三式稱為大地線的微分方程,三角形PP2P1是一微分直角三角形,大地線的微分方程,推導(dǎo)大地線的克萊勞方程,此式就是克萊勞方程，也叫克萊勞定理。定理表明：在旋轉(zhuǎn)橢球面上，大地線各點的平行圈半徑與大地線在該點的大地方位角的正弦的乘積等于常數(shù)。式中常數(shù)C也叫大地線常數(shù) 。大地線的克萊勞方程：表明了大地線在橢球面上的走向。,兩邊積分，易得,當(dāng)大地線穿越赤道時 當(dāng)大地線達極小平行圈時 由此可見，某一大地線常數(shù)等于橢球半徑與該大地線穿越赤道時的大地方位角的正弦的乘積，或者等于該大地線上具有最大緯度的那一點的平行圈半徑。 由克萊勞方程可以寫出 此關(guān)系式

30、可以檢查緯度和方位角計算的正確性。 克萊勞方程還有以下形式：,6.4 將地面觀測值歸算至橢球面 1.概述 參考橢球面是測量計算的基準面。在野外的各種測量都是在地面上進行，觀測的基準線不是各點相應(yīng)的橢球面的法線，而是各點的垂線，各點的垂線與法線存在著垂線偏差。因此不能直接在地面上處理觀測成果，而應(yīng)將地面觀測元素（包括方向和距離等）歸算至橢球面。 1）歸算的意義： 歸算是將地面元素轉(zhuǎn)化為橢球面上的橋梁。 通過歸算，為在橢球面的測量計算提供數(shù)據(jù)。 2）歸算的基本要求： 以橢球面法線為基準線。 地面點沿法線投影到橢球面。 橢球面兩點連線用大地線。 將地面觀測元素加入適當(dāng)?shù)母恼龜?shù)化為橢球面上大地線的相應(yīng)

31、元素。3）地面觀測元素的歸算內(nèi)容： 水平觀測方向、觀測天頂距歸算、地面長度歸算、 天文經(jīng)緯度和方位角歸算等方面,1.將地面觀測的水平方向歸算至橢球面的內(nèi)容： 將水平方向歸算至橢球面上，包括垂線偏差改正、標(biāo)高差改正及截面差改正，習(xí)慣上稱此三項改正為三差改正。 1) 垂線偏差改正u 產(chǎn)生原因：地面上所有水平方向的觀測 都是以垂線為根據(jù)的，而在橢球面上則 要求以該點的法線為依據(jù)。把以垂線為 依據(jù)的地面觀測的水平方向值歸算到以 法線為依據(jù)的方向值而應(yīng)加的改正定義 為垂線偏差改正，以u表示。 如上圖所示，以測站為中心作出單位半徑的輔助球，是垂線偏差，它在子午圈和卯酉圈上的分量分別以、表示，M是地面觀測目

32、標(biāo)m在球面上的投影。 垂線偏差對水平方向的影響是(R-R1)，垂線偏差改正的計算公式是： 式中：，為測站點上的垂線偏差在子午圈及卯酉圈上的分量，它們可在測區(qū)的垂線偏差分量圖中內(nèi)插取得；Am為測站點至照準點的大地方位角；Z1為照準點的天頂距；1為照準點的垂直角。,結(jié)論：垂線偏差改正，不僅與測站的垂線偏差有關(guān)，而且與觀測方向的方位角和垂直角有關(guān)。當(dāng)法線與鉛垂線一致，或者照準點在zz1O面內(nèi)，或者照準點在測站水平面上時，垂線偏差改正為零。 u為0的情況：,（鉛垂線與法線一致）,（照準點在鉛垂線與法線組成的平面內(nèi)）,（照準點在測站水平面內(nèi)）,u最大的情況： 即當(dāng)觀測方向與垂線偏差方向垂直時。 垂線偏差

33、改正適用范圍： 一、二等角測量 三、四等三角測量中，當(dāng) 和較大時,2) 標(biāo)高差改正h 標(biāo)高差改正又稱由照準點高度而引起的改正。不在 同一子午面或同一平行圈上的兩點的法線是不共面 的。當(dāng)進行水平方向觀測時，如果照準點高出橢球 面某一高度，則照準面就不能通過照準點的法線同 橢球面的交點，由此引起的方向偏差的改正叫做標(biāo) 高差改正，以 h表示。 產(chǎn)生原因：由于A、B兩點的法線不在同一平面所 產(chǎn)生（即照準點不在橢球面上） 。 如右圖所示，A為測站點，如果測站點觀測值 已加垂線偏差改正，則可認為垂線同法線一致。這 時測站點在橢球面上或者高出橢球面某一高度，對 水平方向是沒有影響的。這是因為測站點法線不變，

34、 則通過某一照準點只能有一個法截面。 設(shè)照準點高出橢球面的高程為H2，Ana和Bnb分別為 A點及B點的法線，B點法線與橢球面的交點為b。因 為通常Ana和Bnb不在同一平面內(nèi)，所以在A點照準B點得出的法截線是Ab而不是Ab ，因而產(chǎn)生了Ab同Ab方向的差異。按歸算的要求，地面各點都應(yīng)沿自己法線方向投影到橢球面上，即需要的是Ab方向值而不是Ab方向值，因此需加入標(biāo)高差改正數(shù)h，以便將Ab 方向改到Ab方向。,標(biāo)高差改正的計算公式： 公式中：H2為照準點高出橢球面的高程， 它由三部分組成： 其中H常為照準點標(biāo)石中心的正常高，為 高程異常，為照準點的覘標(biāo)高， B2是照準點緯度，M2是相應(yīng)的子午圈曲

35、率半徑。 A1為測站點至照準點的大地方位角。 標(biāo)高差改正主要與照準點的高程有關(guān)。經(jīng)過此項改正后， 便將地面觀測的水平方向值歸化為橢球面上相應(yīng)的法截弧方向。,h為0的三種情況：,H2=0 照準點在橢球面上 A1=0 ，90 ，180 ，270 照準點在測站點的子午圈或平行圈上 B2= 90 照準點在極點上,適用范圍：,一、二等角測量 三、四等三角測量中，當(dāng)海拔高于700m時,2）標(biāo)高差改正 Correction for skew normals,3）截面差改正g 在橢球面上，緯度不同的兩點由于其法線不共面，所以在對向觀測時相對法截弧不重合，應(yīng)當(dāng)用兩點間的大地線代替相對法截弧。這樣將法截弧方向化為

36、大地線方向應(yīng)加的改正叫截面差改正，用g表示。 產(chǎn)生原因：法截線與大地線不一致。 如圖所示，AaB是A至B的法截弧，它在A點處的 大地方位角為A1，ASB 是AB間的大地線，它在 A點的大地方位角是A1，A1與A1 之差g就是截 面差改正 。 截面差改正的計算公式： 式中S為AB間大地線長度， ，N1為測站 點緯度B1相對應(yīng)的卯酉圈曲率半徑。 g為0的情況： A1=0、90、180、270照準點在測站點的子午圈或平行圈上 適用范圍：一等角測量 天文方位角歸算為大地方位角按（3-170）（前面第三章已講） 天文天頂距歸算為大地天頂距按（3-173）,4、三差改正的計算 現(xiàn)行作業(yè)一般規(guī)定，一等三角測

37、量應(yīng)加三差改正，二等三角測量應(yīng)加垂線偏差改正和標(biāo)高差改正，而不加截面差改正；三等和四等三角測量可不加三差改正。但當(dāng) 時或者H2 000m時，則應(yīng)分別考慮加垂線偏差改正和標(biāo)高差改正。在特殊情況下，應(yīng)該根據(jù)測區(qū)的實際情況作具體分析，然后再做出加還是不加改正的規(guī)定。如下表所示：,將地面觀測的長度歸算至橢球面 根據(jù)測邊使用的儀器不同，地面長度的歸算分兩種：一是基線尺量距的歸算，二是電磁波測距的歸算。 1 基線尺量距的歸算 將基線尺量取的長度加上測段傾斜改正后，可以認為它是基線平均高程面上的長度，以表示，現(xiàn)要把它歸算至參考橢球面上的大地線長度S。 1）垂線偏差對長度歸算的影響,2）高程對長度歸算的影響,

38、如果將上式展開級數(shù)，取至二次項， 3）顧及以上兩項，則地面基線長度歸算到橢球面上長度的公式 經(jīng)過以上計算，便得到橢球面上的基線長度。,2. 電磁波測距邊長歸算橢球面 電磁波測距儀測得的長度是連接地面兩點間的直線斜距，也應(yīng)將它歸算到參考橢球面上。 如圖，大地點Q1和Q2的大地高分別為H1和H2。其間用電磁波測距儀測得的斜距為D，現(xiàn)要求大地點在橢球面上沿法線的投影點Q1和Q2間的大地線的長度S。 在工程測量中邊長一般都是幾公里，最長也不過十幾公里，因此，所求的大地線的長度可以認為是半徑 相應(yīng)的圓弧長。 電磁波測距邊長歸算橢球面上的計算公式為： 式中 電磁波測距的歸算公式推導(dǎo)：,將上式按反正弦函數(shù)展

39、開級數(shù)，舍去五次項，得 進一步化簡得：,下式為兩點間的弦長：,電磁波測距邊長歸算的幾何意義： （1）計算公式中右端第二項是由于控制點之高差引起的傾斜改正的主項，經(jīng)過此項改正，測線已變成平距； （2）第三項是由平均測線高出參考橢球面而引起的投影改正，經(jīng)此項改正后，測線已變成弦線； （3）第四項則是由弦長改化為弧長的改正項。,顯然第一項即為經(jīng)高差改正后的平距。,問題 算例見下表，用上述兩個公式計算將電磁波測距邊長歸算至橢球面上。,電磁波測距邊長歸算至橢球面上的計算公式還可用下式表達：,4.7 大地測量主題解算,4.7.1 大地主題解算的一般說明 大地元素:大地經(jīng)度L、大地緯度B、兩點間 的大地線長

40、度S及其正反大地方位角A12、A21。 大地主題解算：如果知道某些大地元素推 求另一些大地元素，這樣的計算問題就叫大地 主題解算，大地主題解算有正解和反解。 大地主題正解：已知Pl點的大地坐標(biāo)(L1，B1)，P1至P2的大地線長S及其大地方位角A12，計算P2點的大地坐標(biāo)(L2，B2)和大地線S在P2點的反方位角A21，這類問題叫做大地主題正解。 大地主題反解： 如果已知P1和P2點的大地坐標(biāo)(L1，B1)和(L2，B2)，計算P1至P2的大地線長S及其正、反方位角A12和A21這類問題叫做大地主題反解。,根據(jù)大地線的長短，主題解算分為: 短距離(400km)， 中距離(1000km)， 長距

41、離(1000km以上) 大地主題解算根據(jù)不同理論基礎(chǔ)可分五類： 1以大地線在大地坐標(biāo)系中的微分方程為基礎(chǔ)，直接在地球橢球面上進行積分運算。 2以白塞爾大地投影為基礎(chǔ)。 3利用地圖投影理論解算大地問題。 4對大地線微分方程進行數(shù)值積分的解法。 5依據(jù)大地線外的其他線為基礎(chǔ)。,1.以大地線在大地坐標(biāo)系中的微分方程為基礎(chǔ)，直接在地球橢球面上進行積分運算。 主要特點：解算精度與距離有關(guān)，距離越長，收斂越慢，因此只適用于較短的距離 典型解法：高斯平均引數(shù)法,2.以白塞爾大地投影為基礎(chǔ) 白塞爾大地主題解算的步驟： 1)按橢球面上的已知值計算球面相應(yīng)值，即實現(xiàn)橢球面向球面的過渡； 2)在球面上解算大地問題；

42、 3)按球面上得到的數(shù)值計算橢球面上的相應(yīng)數(shù)值，即實現(xiàn)從圓球向橢球的過渡。 典型解法：白塞爾大地主題解算 特點：解算精度與距離長短無關(guān)，它既適用于短距離解算，也適用于長距離解算?？蛇m應(yīng)20 000km或更長的距離，這對于國際聯(lián)測，精密導(dǎo)航，遠程導(dǎo)彈發(fā)射等都具有重要意義。 3 利用地圖投影理論解算大地問題 如在地圖投影中，采用橢球面對球面的正形投影和等距離投影以及橢球面對平面的正形投影（如高斯投影），它們都可以用于解算大地主題。這類解法受距離的限制，只在某些特定情況下才比較有利。,4 對大地線微分方程進行數(shù)值積分的解法 這種解法直接進行數(shù)值積分計算以解決大地主題的解算。常用的數(shù)值積分算法有高斯法

43、，龍格庫塔法，牛頓法以及契巴雪夫法等。這種算法易于編寫程序，適用于任意長度距離。缺點是隨著距離的增長，計算工作量大，且精度降低，而在近極地區(qū)，這種方法無能為力。 5 依據(jù)大地線外的其他線為基礎(chǔ) 連接橢球面兩點的媒介除大地線之外，當(dāng)然還有其他一些意義的線，比如弦線 、法截線等。利用弦線解決大地主題實質(zhì)是三維大地測量問題，由電磁波測距得到法截線弧長。,4.7.2 勒讓德級數(shù)式在過已知點P1(L1，B1)且在該點處大地方位角為A12的大地線S上任意一點P2的大地坐標(biāo)(L2，B2)及其方位角A21必是大地線長度S的函數(shù)。 S0時，這些函數(shù)值等于P1點的相應(yīng)數(shù)值因此，可在已知點P1點(S0)上，按麥克勞

44、林公式將Pl和P2點的緯度差、經(jīng)度差及方位角之差展開為大地線長度S的冪級數(shù)。為了計算B、L、A的級數(shù)展開式，關(guān)鍵問題是推求各階導(dǎo)數(shù)。,一階導(dǎo)數(shù)： 二階導(dǎo)數(shù)：,大地測量主題解算,三階導(dǎo)數(shù),大地測量主題解算,引用符號：顧及V/c=1/N及第4、第5階導(dǎo)數(shù)，則得勒讓德級數(shù)式：,大地測量主題解算,勒讓德級數(shù)是大地主題解算的一組基本公式，但他僅適用于邊長短于30km的情況。邊長長的話，級數(shù)收斂慢，且計算復(fù)雜。后來學(xué)者博爾茨、赫里斯托夫、史賴伯對級數(shù)系數(shù)進行了改化。高斯對勒讓德級數(shù)也進行了改化，提出：首先把勒讓德級數(shù)在 P點展開改在大地線長度中點M展開，使級數(shù)公式項數(shù)減少，收斂快，精度高；其次，考慮到求定

45、中點 M 的復(fù)雜性，將 M 點用大地線兩端點平均緯度及平均方位角相對應(yīng)的 m 點來代替，并借助迭代計算便可順利地實現(xiàn)大地主題正解。,4.7.3 高斯平均引數(shù)正算 高斯平均引數(shù)正算公式推導(dǎo)的基本思想：首先把勒讓德級數(shù)在P1點展開改在大地線長度中點M展開，以便級數(shù)公式項數(shù)減少，收斂快，精度高；其次，考慮到求定中點M的復(fù)雜性，將M點用大地線兩端點平均緯度及平均方位角相對應(yīng)的m點來代替， 并借助迭代計算，便可順利地實現(xiàn)大地主題正解。 (1)建立級數(shù)展開式: 兩式相減得：,同理可得:,(2)由于大地線中點M處的緯度和大地方位角均為未知，不能直接用來計算，為此用P1和P2點平均緯度和平均方位角相對應(yīng)的m點

46、代替M點。,大地測量主題解算,將（4-202）（4-204）式中的BM,AM為依據(jù)的導(dǎo)數(shù)值改化為Bm、Am為依據(jù)的導(dǎo)數(shù)值：dB/dS是B和A的函數(shù)，則：將上式展開為以Bm、Am為依據(jù)的級數(shù)：,大地測量主題解算,(3)由大地線微分方程依次求偏導(dǎo)數(shù):,高斯平均引數(shù)正算公式：同理可得：,注意： 從公式可知，欲求，及，必先有及。但由于2和21未知，故精確值尚不知，為此須用逐次趨近的迭代方法進行公式的計算。 除此之外，此方法適合與200公里以下的大地問題解算，其計算經(jīng)緯度計算精度可達到0.0001, 方位角計算精度可達到0.001。,4.7.4 高斯平均引數(shù)反算公式 大地主題反算是已知兩端點的經(jīng)、緯度L

47、1，B1及L2，B2，反求兩點間的大地線長度S及正、反大地方位角A12和A21。 這時，由于經(jīng)差L、緯差 B及平均緯度Bm均為已知，故可依正算公式很容易地導(dǎo)出反算公式。 高斯平均引數(shù)反算公式可以依正算公式導(dǎo)出: 上述兩式的主式為: 將上式代入前兩式，并按L和B集項得：,已知： 求得：,4.7.5 白塞爾大地主題解算方法 白塞爾法解算大地主題的基本思想: 以輔助球面為基礎(chǔ),將橢球面三角形轉(zhuǎn)換為輔助球面的相應(yīng)三角形,由三角形對應(yīng)元素關(guān)系,將橢球面上的大地元素按照白塞爾投影條件投影到輔助球面上，然后在球面上進行大地主題解算，最后再將球面上的計算結(jié)果換算到橢球面上。 這種方法的關(guān)鍵問題是找出橢球面上的

48、大地元素與球面上相應(yīng)元素之間的關(guān)系式,同時也要解決在球面上進行大地主題解算的方法。 在球面上進行大地主題解算 球面上大地主題正算: 已知 求解 球面上大地主題反算: 已知 求解,1、在球面上進行大地主題解算 極球面三角元素間的相互關(guān)系,球面上大地主題正解,球面上大地主題反解方法,2 、橢球面和球面上坐標(biāo)關(guān)系式,在橢球面上與單位球面上的大地線微分方程為:,白塞爾提出如下三個投影條件： 1.橢球面大地線投影到球面上為大圓弧 2.大地線和大圓弧上相應(yīng)點的方位角相等； 3.球面上任意一點緯度等于橢球面上相應(yīng)點的歸化緯度。=,由（4-240）式及第二投影條件（A=）及公式 得：,以上兩式為白塞爾微分方程

49、。表達了橢球面上大地線長度與球面上大圓弧長度，橢球面上經(jīng)差與球面上經(jīng)差的微分關(guān)系。對這組方程進行積分就可以求得S與，L與的關(guān)系式。,3 、白塞爾微分方程的積分,積分得到S與的關(guān)系式： 利用此公式可以計算赤道至大圓弧任意一點Pi的大地線的長度。為計算兩點P1P2間的大地線長度，對這兩點分別使用上式得： 因為S=S2-S1，=2-1，將上兩式相減得:,適合于反算: 適合于正算: 迭代法: 直接法:,現(xiàn)在研究:,將三角函數(shù)冪級數(shù)用倍角函數(shù)代替，合并同類項，積分。截去4倍角項，其值小于0.0001秒。,正算： 反算：,4 白塞爾法大地主題正算步驟,1.計算起點的歸化緯度 2.計算輔助函數(shù)值，解球面三角

50、形可得: 3. 按公式計算相關(guān)系數(shù)A,B,C以及,4.計算球面長度 迭代法: 直接法:,5.計算經(jīng)度差改正數(shù) 6.計算終點大地坐標(biāo)及大地方位角,5 白塞爾法大地主題反算步驟,1.輔助計算,2.用逐次趨近法同時計算起點大地方位角、球面長度及經(jīng)差 ，第一次趨近時，取。,計算下式,重復(fù)上述計算過程2. 3. 計算大地線長度S 4. 計算反方位角,大地主題解算編程實驗,一、實驗?zāi)康?1、提高運用計算機語言編程開發(fā)能力。 2、 加深對大地主題解算計算公式及輔助參數(shù)的理解并掌握其計算步驟。 3 、實現(xiàn)大地主題解算計算機的計算過程提高解算精度。 二、實驗工具 運用自己熟悉的編程開發(fā)語言（C、VC、VB、FO

51、RTRAN、Matlab等） 三、實驗要求提交報告、實驗總結(jié)及編寫的源代碼程序 1、每人獨立完成大地主題解算的計算程序編制，并調(diào)試運行并計算出正確結(jié)果； 2、此次編程實驗應(yīng)上交成果資料（每人一份）： （1）編程思想、編程過程中出現(xiàn)的問題及如何解決問題的； （2）源程序編碼； （3）計算結(jié)果； 四、注意 1、計算所需的變量多。 2、正反算函數(shù)的編寫。 3、函數(shù)的調(diào)用。 4、弧度及角度的轉(zhuǎn)換。 5、輸出的參數(shù)精度正確。,本章小結(jié),1.地球橢球的幾何性質(zhì)。幾個重要概念：法截線、子午圈曲率半徑、卯酉圈曲率半徑、大地線、相對法截弧。 2.地面觀測值歸算到橢球面的原理及過程：方向歸算、長度歸算。 3.大地

52、測量主題解算方法。,習(xí) 題,1試寫出橢球的基本元素及其基本關(guān)系式。,2在控制測量的橢球解算中，常引用下列符號： 、 、 、 ，,3我國解放后主要采用哪兩種參考橢球?其主要參數(shù)是什么? 4繪圖并說明表示橢球面上點位的三種常用坐標(biāo)系統(tǒng)。 5在報紙上經(jīng)?？吹絏 X號輪船在東經(jīng)XXX度，北緯X X度遇險一 6寫出參考橢球體的五個基本元素及相互間的關(guān)系。 7什么叫子午圈？什么叫平行圈？,8參考橢球體扁率的變化，橢球體的形狀發(fā)生怎樣的變形？ 9簡要說明并圖示地面某一點的大地高、正常高以及大地水準面 差距的幾何意義。 10什么是大地測量的基本坐標(biāo)系？有何優(yōu)點？ 11畫圖表示地心緯度坐標(biāo)系和歸化緯度坐標(biāo)系，這

53、兩種坐標(biāo)系 在大地測量中有何意義？ 12用公式表示空間直角坐標(biāo)系和大地坐標(biāo)系之間的關(guān)系。 13何為大地緯度、歸化緯度、地心緯度？三者間有何關(guān)系？,試問它們之間函數(shù)關(guān)系的一個基本共同特點是什么？,類的報導(dǎo)，試問這是指的什么坐標(biāo)系，為什么?,14簡要敘述M、N、R 三種曲率半徑之間的關(guān)系。 15大地坐標(biāo)系和天文坐標(biāo)系各以什么作基準面和基準線？ 16試推證卯酉圈、子午圈曲率半徑的計算公式。 17B00的平行圈是否有可能是法截線？為什么？ 18卯酉圈曲率半徑N與子午圈曲率半徑M何時有最大值？何時有最小值？ 19，為什么說任意方向法截線曲率半徑 隨A的變化是以900為周期的？這 一結(jié)論對橢球問題的解算有什么意義？ 20什么是法線？什么是法截面？它們對橢球解算有什么意義？ 21當(dāng)橢球元素確定之后，橢球面上任意方向法截線曲率半徑的計算值 取決于哪兩個變量？為什么？ 22已知歐拉公式：,23研究平均曲率半徑R對橢球解算有何意義？在我國中緯度地區(qū)R與的 最大差異是多少？試將它對距離化算（用R代替）的影響作一定量分析。,試由橢球基本元素及公式出發(fā)，用兩種方法導(dǎo)出計算任意方 向法截線曲率半徑的公式和平均曲率半徑R的公式。,24在推導(dǎo)計算子午線弧長公式時，為什么要從赤道起算？若欲求緯度B1和B2間的子午線弧長( 0)，如何計算？ 25當(dāng)子午線弧長不超過45km時，則可將其視為圓弧，試論