1、經(jīng)緯度坐標(biāo)與高斯坐標(biāo)的轉(zhuǎn)換代碼/功能說明： 將絕對高斯坐標(biāo)(y,x)轉(zhuǎn)換成絕對的地理坐標(biāo)(wd,jd)。 / double y; 輸入?yún)?shù): 高斯坐標(biāo)的橫坐標(biāo)，以米為單位 / double x; 輸入?yún)?shù): 高斯坐標(biāo)的縱坐標(biāo)，以米為單位/ short DH; 輸入?yún)?shù): 帶號，表示上述高斯坐標(biāo)是哪個帶的/ double *L; 輸出參數(shù): 指向經(jīng)度坐標(biāo)的指針，其中經(jīng)度坐標(biāo)以秒為單位/ double *B; 輸出參數(shù): 指向緯度坐標(biāo)的指針，其中緯度坐標(biāo)以秒為單位void GaussToGeo(double y, double x, short DH, double *L, double *B,

2、double LP) double l0; / 經(jīng)差 double tf; / tf = tg(Bf0),注意要將Bf轉(zhuǎn)換成以弧度為單位 double nf ; / n = y * sqrt( 1 + etf * 2) / c, 其中etf = e'*2 * cos(Bf0) * 2 double t_l0; / l0，經(jīng)差，以度為單位 double t_B0; / B0，緯度，以度為單位 double Bf0; / Bf0 double etf; / etf,其中etf = e'*2 * cos(Bf0) * 2 double X_3 ; double PI=3.141592

3、65358979; double b_e2=0.0067385254147; double b_c=6399698.90178271; X_3 = x / 1000000.00 - 3 ; / 以兆米（1000000）為單位 / 對于克拉索夫斯基橢球，計(jì)算Bf0 Bf0 = 27.11115372595 + 9.02468257083 * X_3 - 0.00579740442 * pow(X_3,2) - 0.00043532572 * pow(X_3,3) + 0.00004857285 * pow(X_3,4) + 0.00000215727 * pow(X_3,5) - 0.00000

4、019399 * pow(X_3,6) ; tf = tan(Bf0*PI/180); / tf = tg(Bf),注意這里將Bf轉(zhuǎn)換成以弧度為單位 etf = b_e2 * pow(cos(Bf0*PI/180),2); / etf = e'*2 * cos(Bf) * 2 nf = y * sqrt( 1 + etf ) / b_c; / n = y * sqrt( 1 + etf * 2) / c / 計(jì)算緯度，注意這里計(jì)算出來的結(jié)果是以度為單位的 t_B0 = Bf0 - (1.0+etf) * tf / PI * (90.0 * pow(nf,2) - 7.5 * (5.0

5、+ 3 * pow(tf,2) + etf - 9 * etf * pow(tf,2) * pow(nf,4) + 0.25 * (61 + 90 * pow(tf,2) + 45 * pow(tf,4) * pow(nf,6) ; / 計(jì)算經(jīng)差，注意這里計(jì)算出來的結(jié)果是以度為單位的 t_l0 = (180 * nf - 30 * ( 1 + 2 * pow(tf,2) + etf ) * pow(nf,3) + 1.5 * (5 + 28 * pow(tf,2) + 24 * pow(tf,4) * pow(nf,5) / ( PI * cos(Bf0*PI/180) ) ; l0 = (t

6、_l0 * 3600.0); / 將經(jīng)差轉(zhuǎn)成秒 if (LP = -1000) *L = (double)(DH * 6 - 3) * 3600.0 + l0); / 根據(jù)帶號計(jì)算出以秒為單位的絕對經(jīng)度，返回指針 else *L = LP * 3600.0 + l0; / 根據(jù)帶號計(jì)算出以秒為單位的絕對經(jīng)度，返回指針 /- *B = (double)(t_B0 * 3600.0) ; / 將緯差轉(zhuǎn)成秒，并返回指針/ 功能說明： （1）將地理坐標(biāo)(wd,jd)轉(zhuǎn)換成絕對的高斯坐標(biāo)(y,x) （2）本函數(shù)支持基于六度帶（或三度帶）、克拉索夫斯基橢球進(jìn)行轉(zhuǎn)換 / 適用范圍： 本函數(shù)適用于將地球東半球

7、中北半球（即東經(jīng)0度到東經(jīng)180度，北緯0度至90度）范圍 內(nèi)所有地理坐標(biāo)到高斯坐標(biāo)的轉(zhuǎn)換 / 使用說明： 調(diào)用本函數(shù)后返回的結(jié)果應(yīng)在滿足精度的條件下進(jìn)行四舍五入 / double jd; 輸入?yún)?shù): 地理坐標(biāo)的經(jīng)度，以秒為單位/ double wd; 輸入?yún)?shù): 地理坐標(biāo)的緯度，以秒為單位/ short DH; 輸入?yún)?shù): 三度帶或六度帶的帶號/ 六度帶(三度帶)的帶號是這樣得到的：從東經(jīng)0度到東經(jīng)180度自西向東按每6度(3度)順序編號 (編號從1開始)，這個順序編號就稱為六度帶(三度帶)的帶號。因此，六度帶的帶號的范圍是1-30， 三度帶的帶號的范圍是1-60。 如果一個點(diǎn)在圖號為TH的

8、圖幅中，那麼該點(diǎn)所處的六度帶的帶號就可以這樣得到：將該圖號的 第3、4位組成的字符串先轉(zhuǎn)換成數(shù)字，再減去30。例如某點(diǎn)在圖幅06490701中，該點(diǎn)所在的帶號就 是49-30，即19。 如果調(diào)用本函數(shù)去進(jìn)行一般的從地理坐標(biāo)到基于六度帶高斯坐標(biāo)的變換（非鄰帶轉(zhuǎn)換），則參 數(shù)DH的選取按前一段的方法去確定。 如果調(diào)用本函數(shù)去進(jìn)行基于六度帶鄰帶轉(zhuǎn)換，則參數(shù)DH的選取先按上述方法去確定，然后看是 往前一個帶還是后一個帶進(jìn)行鄰帶轉(zhuǎn)換再確定是加1還是減1。 /void GeoToGauss(double jd, double wd, short DH, short DH_width, double *y,

9、 double *x, double LP) double t; / t=tgB double L; / 中央經(jīng)線的經(jīng)度 double l0; / 經(jīng)差 double jd_hd,wd_hd; / 將jd、wd轉(zhuǎn)換成以弧度為單位 double et2; / et2 = (e' * 2) * (cosB * 2) double N; / N = C / sqrt(1 + et2) double X; / 克拉索夫斯基橢球中子午弧長 double m; / m = cosB * PI/180 * l0 double tsin,tcos; / sinB,cosB double PI=3.14

10、159265358979; double b_e2=0.0067385254147; double b_c=6399698.90178271; jd_hd = jd / 3600.0 * PI / 180.0 ; / 將以秒為單位的經(jīng)度轉(zhuǎn)換成弧度 wd_hd = wd / 3600.0 * PI / 180.0 ; / 將以秒為單位的緯度轉(zhuǎn)換成弧度 / 如果不設(shè)中央經(jīng)線（缺省參數(shù): -1000），則計(jì)算中央經(jīng)線， / 否則，使用傳入的中央經(jīng)線，不再使用帶號和帶寬參數(shù) /L = (DH - 0.5) * DH_width ; / 計(jì)算中央經(jīng)線的經(jīng)度 if (LP = -1000) L = (DH

11、 - 0.5) * DH_width ; / 計(jì)算中央經(jīng)線的經(jīng)度 else L = LP ; l0 = jd / 3600.0 - L ; / 計(jì)算經(jīng)差 tsin = sin(wd_hd); / 計(jì)算sinB tcos = cos(wd_hd); / 計(jì)算cosB / 計(jì)算克拉索夫斯基橢球中子午弧長X X = 111134.8611 / 3600.0 * wd - (32005.7799 * tsin + 133.9238 * pow(tsin,3) + 0.6976 * pow(tsin,5) + 0.0039 * pow(tsin,7) ) * tcos; et2 = b_e2 * pow

12、(tcos,2) ; / et2 = (e' * 2) * (cosB * 2) N = b_c / sqrt( 1 + et2 ) ; / N = C / sqrt(1 + et2) t = tan(wd_hd); / t=tgB m = PI/180 * l0 * tcos; / m = cosB * PI/180 * l0 *x = X + N * t * ( 0.5 * pow(m,2) + (5.0 - pow(t,2) + 9.0 * et2 + 4 * pow(et2,2) * pow(m,4)/24.0 + (61.0 - 58.0 * pow(t,2) + pow(t

13、,4) * pow(m,6) / 720.0 ) ; *y = N * ( m + ( 1.0 - pow(t,2) + et2 ) * pow(m,3) / 6.0 + ( 5.0 - 18.0 * pow(t,2) + pow(t,4) + 14.0 * et2 - 58.0 * et2 * pow(t,2) ) * pow(m,5) / 120.0 );3度帶與6度帶 1我國采用6度分帶和3度分帶：12.5萬及15萬的地形圖采用6度分帶投影，即經(jīng)差為6度，從零度子午線開始，自西向東每個經(jīng)差6度為一投影帶，全球共分60個帶，用1，2，3，4，5，表示即東經(jīng)06度為第一帶，其中央經(jīng)線的經(jīng)度為

14、東經(jīng)3度，東經(jīng)612度為第二帶，其中央經(jīng)線的經(jīng)度為9度。11萬的地形圖采用3度分帶，從東經(jīng)1.5度的經(jīng)線開始，每隔3度為一帶，用1，2，3，表示，全球共劃分120個投影帶，即東經(jīng)1.5 4.5度為第1帶，其中央經(jīng)線的經(jīng)度為東經(jīng)3度，東經(jīng)4.57.5度為第2帶，其中央經(jīng)線的經(jīng)度為東經(jīng)6度我省位于東經(jīng)113度東經(jīng)120度之間，跨第38、39、40共計(jì)3個帶，其中東經(jīng)115.5度以西為第38帶，其中央經(jīng)線為東經(jīng)114度；東經(jīng)115.5118.5度為39帶，其中央經(jīng)線為東經(jīng)117度；東經(jīng)118.5度以東到山海關(guān)為40帶，其中央經(jīng)線為東經(jīng)120度。地形圖上公里網(wǎng)橫坐標(biāo)前2位就是帶號，例如：15萬地形圖上

15、的橫坐標(biāo)為20345486，其中20即為帶號，345486為橫坐標(biāo)值。2當(dāng)?shù)刂醒虢?jīng)線經(jīng)度的計(jì)算六度帶中央經(jīng)線經(jīng)度的計(jì)算：當(dāng)?shù)刂醒虢?jīng)線經(jīng)度6°×當(dāng)?shù)貛?°，例如：地形圖上的橫坐標(biāo)為20345，其所處的六度帶的中央經(jīng)線經(jīng)度為：6°×203°117°（適用于125萬和15萬地形圖）。三度帶中央經(jīng)線經(jīng)度的計(jì)算：中央經(jīng)線經(jīng)度3°×當(dāng)?shù)貛枺ㄟm用于11萬地形圖）。3、如何計(jì)算當(dāng)?shù)氐闹醒胱游缇€？ 當(dāng)?shù)刂醒胱游缇€決定于當(dāng)?shù)氐闹苯亲鴺?biāo)系統(tǒng)，首先確定您的直角坐標(biāo)系統(tǒng)是3度帶還是6度帶投影公式推算:6度帶中央子午線計(jì)算公式：

16、當(dāng)?shù)亟?jīng)度/6=N；中央子午線L=6 * N （帶號） 當(dāng)沒有除盡，N有余數(shù)時，中央子午線L=6*N - 3 3度帶中央子午線計(jì)算公式： 當(dāng)?shù)亟?jīng)度/3=N；中央子午線L=3 X N 我國的經(jīng)度范圍西起 73°東至135°，可分成六度帶十一個（13號帶23號帶），各帶中央經(jīng)線依次為（75°、81°、123°、129°、135°）；三度帶二十二 個（24號帶45號帶）。各帶中央經(jīng)線依次為（72°、75°、132°、135°）；六度帶可用于中小比例尺（如 1：250000）測圖，三度帶可用于大比例尺（如 1：10000）測圖，城建坐標(biāo)多采用三度帶的高斯投影4、如何判斷投影坐標(biāo)是