地方独立坐标系与WGS_84坐标系转换方法及应用

[2]　.[M].:,1996[3]　.[J].,2009(6):26-27[4]　,.[M].:,2001[5]　,.[J].:,1990[6]　,,.[M].:,1989[7]　.[J].,2009,29(3):128-130[8]　.[J].,2009,35(5):8-9:2010-05-27:(1986—),,2008,。:1672-7479(2010)04-0008-03WGS-84蒋小军　宋占峰　吴清华(,　410004)ConversionMethodbetweenLocalIndependentCoordinateSystemandWGS-84CoordinateSystemasWellasItsApplicationsJiangXiaojun　SongZhanfeng　WuQinghua　　　　根据n(n≥3)个点的地方独立坐标及对应的WGS-84坐标,结合平面坐标转换模型、布尔莎模型和三维坐标差转换模型,完成了地方独立坐标系与WGS-84坐标系转换参数的计算。基于该转换模型,利用GoogleEarthCOMAPI、LibKML等开发接口,完成了AutoCAD数字地形图数据到GoogleEarth上的三维表示。　地方独立坐标系　WGS-84坐标系　坐标转换　GoogleEarth:U412.3　　　:A1　、,,。,,,,,、。:①,;②,,;③,、。,、GIS、,,GoogleEarth,。GoogleEarth。,,、、。,。GoogleEarth,WGS-84,54,54,54WGS-84,,。2　(1)WGS-848　　　20104a=(6378137±2)m;()GM=(3986005×108±0.6×108)m3/s2;　C2.0=-484.16685×10±0.6×10-6;ω=(7292115×10-11±0.15×10-11)rad/s;WGS-84f=1/298.257223563。(2)54()a=6378245m;f=1/298.3;,、54,Δa,1。1　54　3　2。3.1　同一参考系下的平面坐标转换模型XS=ΔX+k×(XTcosα-YTsinα)YS=ΔY+k×(XTsinα+YTcosα)(1)　ΔX、ΔY———xy;k、α———。3.2　布尔莎(Bursa)模型xyzS=ΔxΔyΔz+(1+k)xyzT+0rz-ry-rz0rxry-rx0xyzT(2),Δx,Δy,Δz,k,rx,ry,rz。(32,C++15～16),2　WGS-84　,,。3.3　三维坐标差转换模型,n(n≥3),i(0i≤n,i≠j)(xi,yi,zi),j(0j≤n,j≠i)(xj,yj,zj),Δxij=xj-xi,Δyij=yj-yi,Δzij=zj-zi,ΔxijΔyijΔzijS=(1+k)ΔxijΔyijΔzijT+0rz-ry-rz0rxry-rx0ΔxijΔyijΔzijTΔxijΔyijΔzijS-ΔxijΔyijΔzijT=Δxij0-ΔzijΔyijΔyijΔzij0-ΔxijΔzij-ΔyijΔxij0Tkrxryrz(3)　　b=Ax9WGS-84:　　,,Gauss-Jordan、。,:①,,,,。②,。,。,x。,,。。,Gauss-Jordan、。,。,,。:,。4　。,,112°50′00.0000″E。1。1　WorldGeodeticSystem1984X/mY/mZ()/m(B)(L)(H)769393985.159050581.050072.227028°09′04.6877″112°59′45.3821″64.5884JT0293864.182049763.000081.037228°09′00.5381″112°59′15.4378″73.3985TK0194799.524049909.745985.570528°09′30.9607″112°59′20.5316″77.9322TK0294750.957449614.1457108.458928°09′29.3034″112°59′09.7120″100.8205TK0395368.305349366.258373.496428°09′49.2898″112°59′00.4378″65.8583TK0495406.716249422.603673.839128°09′50.5528″112°59′02.4913″66.20104.1　四参数计算54x:Δx=3021287.8529my:Δy=465391.0823m:k=1.0×10-6:α=-1975.2243″GPS,。4.2　七参数计算54,,54,(2)。2　ΔxΔyΔzkrxryrzWGS-84※54-3424061.78188904563.30314615438.4672-1.66125472912°11′54.0913″-35°11′42.7900″-21°30′28.6601″54※WGS-84-3928813.43218582189.42014824084.8990-1.661561748-14°09′36.9082″33°58′04.8015″22°22′38.6847″4.3　实验结果,AutoCAD,GoogleEarth,GoogleEarth,GoogleEarth,(3)。,GoogleEarth(4)。5　WGS-84,:①,,。,54,。②,,,,。(下转第13页)10　　　20104。,,,(2)。2　5　,CPⅠCPⅡ10mm/km25mm/km,,,,,。、、CPI。[1]　,,.[M].:,2005[2]　,.[M].:,1997[3]　.(TB10601—2009)[S].:,2010[4]　.(TB10101—2009)[S].:,2010[5]　　.[M].:,2008[6]　　,.MicroStationV8[M].:,2002[7]　.[M].:,1999[8]　.(GB50026—2007)[S].:,20084　GoogleEarth　　　(上接第10页)3　AutoCAD　[1]　.[J].,2004(3):16-17[2]　,,,.WGS-84BJ-54[J].,2007,30(5):33-42[3]　,.[J].,2007,27(3):25-38[4]　,,,.BURSABJ54WGS84[J].,2006,18(3):32-35[5]　　,.GPS[M].:,2005[6]　,.GPS[M].:,2005[7]　,,.[M].:,2001[8]　(),.,.C++[M].:,200513: