經(jīng)緯度坐標(biāo)及高斯坐標(biāo)的轉(zhuǎn)換代碼

..經(jīng)緯度坐標(biāo)與高斯坐標(biāo)的轉(zhuǎn)換代碼/功能說(shuō)明：將絕對(duì)高斯坐標(biāo)(y,x)轉(zhuǎn)換成絕對(duì)的地理坐標(biāo)(wd,jd)。///doubley;輸入?yún)?shù):高斯坐標(biāo)的橫坐標(biāo)，以米為單位//doublex;輸入?yún)?shù):高斯坐標(biāo)的縱坐標(biāo)，以米為單位//shortDH;輸入?yún)?shù):帶號(hào)，表示上述高斯坐標(biāo)是哪個(gè)帶的//double*L;輸出參數(shù):指向經(jīng)度坐標(biāo)的指針，其中經(jīng)度坐標(biāo)以秒為單位//double*B;輸出參數(shù):指向緯度坐標(biāo)的指針，其中緯度坐標(biāo)以秒為單位voidGaussToGeo(doubley,doublex,shortDH,double*L,double*B,doubleLP){doublel0;//經(jīng)差doubletf;//tf=tg(Bf0),注意要將Bf轉(zhuǎn)換成以弧度為單位doublenf;//n=y*sqrt(1+etf**2)/c,其中etf=e'**2*cos(Bf0)**2doublet_l0;//l0，經(jīng)差，以度為單位doublet_B0;//B0，緯度，以度為單位doubleBf0;//Bf0doubleetf;//etf,其中etf=e'**2*cos(Bf0)**2doubleX_3;doublePI=3.979;doubleb_e2=0.47;doubleb_c=6399698.90178271;X_3=x/1000000.00-3;//以兆米〔1000000〕為單位//對(duì)于克拉索夫斯基橢球，計(jì)算Bf0Bf0=27.+9.*X_3-0.*pow(X_3,2)-0.*pow(X_3,3)+0.*pow(X_3,4)+0.*pow(X_3,5)-0.*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)**2nf=y*sqrt(1+etf)/b_c;//n=y*sqrt(1+etf**2)/c//計(jì)算緯度，注意這里計(jì)算出來(lái)的結(jié)果是以度為單位的t_B0=Bf0-(1.0+etf)*tf/PI*(90.0*pow(nf,2)-7.5*(5.0+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ì)算出來(lái)的結(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_l0*3600.0);//將經(jīng)差轉(zhuǎn)成秒if(LP==-1000){*L=(double)((DH*6-3)*3600.0+l0);//根據(jù)帶號(hào)計(jì)算出以秒為單位的絕對(duì)經(jīng)度，返回指針}else{*L=LP*3600.0+l0;//根據(jù)帶號(hào)計(jì)算出以秒為單位的絕對(duì)經(jīng)度，返回指針}//----------------------------------*B=(double)(t_B0*3600.0);//將緯差轉(zhuǎn)成秒，并返回指針}/功能說(shuō)明：〔1〕將地理坐標(biāo)(wd,jd)轉(zhuǎn)換成絕對(duì)的高斯坐標(biāo)(y,x)〔2〕本函數(shù)支持基于六度帶〔或三度帶〕、克拉索夫斯基橢球進(jìn)展轉(zhuǎn)換//適用圍：本函數(shù)適用于將地球東半球中北半球〔即東經(jīng)0度到東經(jīng)180度，北緯0度至90度〕圍所有地理坐標(biāo)到高斯坐標(biāo)的轉(zhuǎn)換//使用說(shuō)明：調(diào)用本函數(shù)后返回的結(jié)果應(yīng)在滿足精度的條件下進(jìn)展四舍五入///doublejd;輸入?yún)?shù):地理坐標(biāo)的經(jīng)度，以秒為單位//doublewd;輸入?yún)?shù):地理坐標(biāo)的緯度，以秒為單位//shortDH;輸入?yún)?shù):三度帶或六度帶的帶號(hào)/六度帶(三度帶)的帶號(hào)是這樣得到的：從東經(jīng)0度到東經(jīng)180度自西向東按每6度(3度)順序編號(hào)(編號(hào)從1開(kāi)場(chǎng))，這個(gè)順序編號(hào)就稱為六度帶(三度帶)的帶號(hào)。因此，六度帶的帶號(hào)的圍是1-30，三度帶的帶號(hào)的圍是1-60。如果一個(gè)點(diǎn)在圖號(hào)為T(mén)H的圖幅中，那麼該點(diǎn)所處的六度帶的帶號(hào)就可以這樣得到：將該圖號(hào)的第3、4位組成的字符串先轉(zhuǎn)換成數(shù)字，再減去30。例如某點(diǎn)在圖幅06490701中，該點(diǎn)所在的帶號(hào)就是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的選取先按上述方法去確定，然后看是往前一個(gè)帶還是后一個(gè)帶進(jìn)展鄰帶轉(zhuǎn)換再確定是加1還是減1。/voidGeoToGauss(doublejd,doublewd,shortDH,shortDH_width,double*y,double*x,doubleLP){doublet;//t=tgBdoubleL;//中央經(jīng)線的經(jīng)度doublel0;//經(jīng)差doublejd_hd,wd_hd;//將jd、wd轉(zhuǎn)換成以弧度為單位doubleet2;//et2=(e'**2)*(cosB**2)doubleN;//N=C/sqrt(1+et2)doubleX;//克拉索夫斯基橢球中子午弧長(zhǎng)doublem;//m=cosB*PI/180*l0doubletsin,tcos;//sinB,cosBdoublePI=3.979;doubleb_e2=0.47;doubleb_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)線，不再使用帶號(hào)和帶寬參數(shù)//L=(DH-0.5)*DH_width;//計(jì)算中央經(jīng)線的經(jīng)度if(LP==-1000){L=(DH-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ì)算sinBtcos=cos(wd_hd);//計(jì)算cosB//計(jì)算克拉索夫斯基橢球中子午弧長(zhǎng)XX=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(tcos,2);//et2=(e'**2)*(cosB**2)N=b_c/sqrt(1+et2);//N=C/sqrt(1+et2)t=tan(wd_hd);//t=tgBm=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,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度分帶：1∶2.5萬(wàn)及1∶5萬(wàn)的地形圖采用6度分帶投影，即經(jīng)差為6度，從零度子午線開(kāi)場(chǎng)，自西向東每個(gè)經(jīng)差6度為一投影帶，全球共分60個(gè)帶，用1，2，3，4，5，……表示．即東經(jīng)0～6度為第一帶，其中央經(jīng)線的經(jīng)度為東經(jīng)3度，東經(jīng)6～12度為第二帶，其中央經(jīng)線的經(jīng)度為9度。1∶1萬(wàn)的地形圖采用3度分帶，從東經(jīng)1.5度的經(jīng)線開(kāi)場(chǎng)，每隔3度為一帶，用1，2，3，……表示，全球共劃分120個(gè)投影帶，即東經(jīng)1.5～4.5度為第1帶，其中央經(jīng)線的經(jīng)度為東經(jīng)3度，東經(jīng)4.5～7.5度為第2帶，其中央經(jīng)線的經(jīng)度為東經(jīng)6度．我省位于東經(jīng)113度－東經(jīng)120度之間，跨第38、39、40共計(jì)3個(gè)帶，其中東經(jīng)115.5度以西為第38帶，其中央經(jīng)線為東經(jīng)114度；東經(jīng)115.5～118.5度為39帶，其中央經(jīng)線為東經(jīng)117度；東經(jīng)118.5度以東到山海關(guān)為40帶，其中央經(jīng)線為東經(jīng)120度。地形圖上公里網(wǎng)橫坐標(biāo)前2位就是帶號(hào)，例如：1∶5萬(wàn)地形圖上的橫坐標(biāo)為20345486，其中20即為帶號(hào)，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ù)貛?hào)－3°，例如：地形圖上的橫坐標(biāo)為20345，其所處的六度帶的中央經(jīng)線經(jīng)度為：6°×20－3°＝117°〔適用于1∶2．5萬(wàn)和1∶5萬(wàn)地形圖〕。三度帶中央經(jīng)線經(jīng)度的計(jì)算：中央經(jīng)線經(jīng)度＝3°×當(dāng)?shù)貛?hào)〔適用于1∶1萬(wàn)地形圖〕。3、如何計(jì)算當(dāng)?shù)氐闹醒胱游缇€？當(dāng)?shù)刂醒胱游缇€決定于當(dāng)?shù)氐闹苯亲鴺?biāo)系統(tǒng)，首先確定您的直角坐標(biāo)系統(tǒng)是3度帶還是6度帶投影公式推算:6度帶中央子午線計(jì)算公式：當(dāng)?shù)亟?jīng)度/6=N；中央子午線L=6*N〔帶號(hào)〕當(dāng)沒(méi)有除盡，N有余數(shù)時(shí)，中央子午線L=6*N-33度帶中央子午線計(jì)算公式：當(dāng)?shù)亟?jīng)度/3=N；中央子午線L=3XN我國(guó)的經(jīng)度圍西起73°東至135°，可分成六度帶十一個(gè)〔13號(hào)帶—23號(hào)帶〕，各帶中央經(jīng)線依次為〔75°、81°、……123°、129°、135°〕；三度帶二十二個(gè)〔24號(hào)帶—45號(hào)帶〕。各帶中央經(jīng)線依次為〔72°、75°、……132°、135°〕；六度帶可用于中小比例尺〔如1：250000〕測(cè)圖，三度帶可用于大比例尺〔如1：10000〕