大地坐標(biāo)與直角空間坐標(biāo)轉(zhuǎn)換計(jì)算公式-資料下載頁
2025-06-24 14:47本頁面
【正文】 nt ProjNo。 int ZoneWide。 ////帶寬 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval。double e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI。iPI = 。 ////。 a = 。 f = 。 //54年北京坐標(biāo)系參數(shù) ////a=。 f=1/。 //80年西安坐標(biāo)系參數(shù) ZoneWide = 6。 ////6度帶寬 ProjNo = (int)(X/1000000L) 。 //查找?guī)杔ongitude0 = (ProjNo1) * ZoneWide + ZoneWide / 2。 longitude0 = longitude0 * iPI 。 //中央經(jīng)線X0 = ProjNo*1000000L+500000L。 Y0 = 0。 xval = XX0。 yval = YY0。 //帶內(nèi)大地坐標(biāo)e2 = 2*ff*f。e1 = ((1e2))/(+sqrt(1e2))。ee = e2/(1e2)。M = yval。u = M/(a*(1e2/43*e2*e2/645*e2*e2*e2/256))。fai = u+(3*e1/227*e1*e1*e1/32)*sin(2*u)+(21*e1*e1/1655*e1*e1*e1*e1/32)*sin(4*u)+(151*e1*e1*e1/96)*sin(6*u)+(1097*e1*e1*e1*e1/512)*sin(8*u)。C = ee*cos(fai)*cos(fai)。T = tan(fai)*tan(fai)。NN = a/sqrt(*sin(fai)*sin(fai))。R = a*(1e2)/sqrt((1e2*sin(fai)*sin(fai))*(1e2*sin(fai)*sin(fai))*(1e2*sin(fai)*sin(fai)))。D = xval/NN。//計(jì)算經(jīng)度(Longitude) 緯度(Latitude)longitude1 = longitude0+(D(1+2*T+C)*D*D*D/6+(52*C+28*T3*C*C+8*ee+24*T*T)*D*D*D*D*D/120)/cos(fai)。latitude1 = fai (NN*tan(fai)/R)*(D*D/2(5+3*T+10*C4*C*C9*ee)*D*D*D*D/24+(61+90*T+298*C+45*T*T256*ee3*C*C)*D*D*D*D*D*D/720)。//轉(zhuǎn)換為度 DD*longitude = longitude1 / iPI。 *latitude = latitude1 / iPI。}如果有需要程序的，可以直接跟我聯(lián)系，呵呵附：高斯正反算參數(shù)pi= ※※ //π長半軸a=。 扁率f=。 //54年北京坐標(biāo)系參數(shù) 長半軸a=。 扁率f=1/。 //80年西安坐標(biāo)系參數(shù)長半軸a=6378137m；扁率f=1:。//WGS84坐標(biāo)系
試題試卷相關(guān)推薦
鄂ICP備17016276號-1