汪林,閆偉,陳艷云
(1.天津市測繪院,天津 300381; 2.國家海洋信息中心,天津 300171)
?
基于GeoDatabase的天津地方坐標系至CGCS2000坐標轉(zhuǎn)換實現(xiàn)
汪林1*,閆偉1,陳艷云2
(1.天津市測繪院,天津 300381; 2.國家海洋信息中心,天津 300171)
作為測繪成果的GIS數(shù)據(jù)多采用GeoDatabase數(shù)據(jù)模型進行存儲和管理,為避免重復測繪,需將GeoDatabase數(shù)據(jù)結(jié)構(gòu)的GIS數(shù)據(jù)從天津地方坐標轉(zhuǎn)換至CGCS2000坐標系。本文從GeoDatabase數(shù)據(jù)模型的介紹入手,總結(jié)了坐標轉(zhuǎn)換相關(guān)的理論,并給出了天津地方坐標系轉(zhuǎn)換至CGCS2000坐標系的方法和步驟,實現(xiàn)了GeoDatabase數(shù)據(jù)結(jié)構(gòu)的GIS數(shù)據(jù)從天津地方坐標系轉(zhuǎn)換至CGCS2000坐標系的轉(zhuǎn)換。最后分別從轉(zhuǎn)換前后從圖形一致性、線性方向一致性、屬性一致性、坐標轉(zhuǎn)換精度和坐標轉(zhuǎn)換效率加以論證。
坐標轉(zhuǎn)換;天津地方坐標系;CGCS2000坐標系
根據(jù)《中華人民共和國測繪法》,國家建立全國統(tǒng)一的大地坐標系。因建設(shè)、城市規(guī)劃和科學研究的需要,大城市和國家重大工程項目可建立相對獨立的平面坐標系統(tǒng)。我國于上世紀50年代和80年代,分別建立了國家大地坐標系統(tǒng)—北京54坐標系和西安80坐標系,測制了各種比例尺地形圖,為國民經(jīng)濟和社會發(fā)展提供了基礎(chǔ)的測繪保障。上述兩種坐標系都是參心坐標系。隨著社會的進步,經(jīng)濟建設(shè)、國防建設(shè)、社會發(fā)展和科學研究等都迫切的需要原點位于地球質(zhì)量中心的坐標系統(tǒng)(簡稱地心坐標系)。我國自 2008 年7 月 1 日起正式啟用 2000 國家大地坐標系[1](簡稱CGCS2000坐標系)。采用地心坐標系,有利于采用現(xiàn)代空間技術(shù)對坐標系進行維護和快速更新,測定高精度大地控制點三維坐標,并提高測圖工作效率。
在啟用CGCS2000坐標系前,天津市為避免高斯投影變形帶來的不便,沿用的平面測繪基準為1990年天津市任意直角坐標系(簡稱天津90坐標系)?,F(xiàn)有的天津90坐標系采用的是北京54坐標系的橢球及相關(guān)參數(shù),是一個局部參心坐標系統(tǒng)。歷史原因造成該基準缺乏時空動態(tài)信息,現(xiàn)實性差,難以維持且精度低,作為現(xiàn)代大地測量的起算基準已經(jīng)不能滿足天津市城市建設(shè)的快速發(fā)展對高精度標準的要求。為避免重復測繪,需要將已有的成果從天津90地方坐標系轉(zhuǎn)換成CGCS2000坐標系。而GIS數(shù)據(jù)作為測繪成果一種存儲管理的重要成果形式,多采用GeoDatabase的空間數(shù)據(jù)模型進行存儲,為此本文嘗試從GeoDatabase的模型出發(fā),探討天津地方坐標系轉(zhuǎn)換至CGCS2000坐標系的理論與方法。
2.1 GeoDatabase介紹
GeoDatabase是一個關(guān)系數(shù)據(jù)庫基礎(chǔ)上擴展的“面向?qū)ο髷?shù)據(jù)庫”,又稱對象關(guān)系模型[2]。它將現(xiàn)實世界抽象為由若干對象類組成的數(shù)據(jù)模型,每個對象類有對應(yīng)的屬性、行為和規(guī)則。GeoDatabase可將空間信息和屬性信息有效的組織在一起,能夠處理復雜的對象。
GeoDatabase作為存儲GIS數(shù)據(jù)的容器有三種類型:Personal Geodatabase、File Geodatabase和ArcSDE Geodatabase。Personal Geodatabase用微軟Access文件存儲數(shù)據(jù),是將幾何數(shù)據(jù)存儲在二進制大字段中的MDB資料庫檔案,最大支持 2 GB。File GeoDatabase采用二進制文件格式存儲,每個表限制為 1 TB。ArcSDE GeoDatabase采用C/S架構(gòu)的企業(yè)級GeoDatabase,支持多用戶共享、編輯和使用數(shù)據(jù)庫。
在GeoDatabase中,存儲的數(shù)據(jù)有四種:矢量數(shù)據(jù),柵格數(shù)據(jù),不規(guī)則三角網(wǎng)和地址與定位數(shù)據(jù)。其中矢量數(shù)據(jù)以要素集或要素類的形式實現(xiàn),要素類具有表示幾何形狀和空間位置的特殊字段,名稱為“Shape”,字段類型為“Geometry”。Shape字段可以是以下七種幾何類型的一種:點、多點、折線、多邊形,注記,尺寸注記和多面體。以下主要介紹常見的點、線和面三種。
點是0維的幾何圖形,有一對(x,y)坐標和可選的高程值(Z)、測量值(M)。折線是相交或不相交的路徑有序幾何,折線是由路徑Path構(gòu)成。多邊形用來表示面狀要素的幾何形狀,多邊形是由閉合的環(huán)Ring構(gòu)成的。曲線Segment是路徑和環(huán)的基本組成,而曲線Segment包含四種基本類型:直線Line、圓弧CircularArc、橢圓弧EllipticArc和貝塞爾曲線BezierCurve。
2.2 坐標轉(zhuǎn)換理論
測量坐標系常見的形式分為大地坐標系、空間直角坐標系和平面直角坐標系[3]。同一基準下,大地坐標和空間直角坐標之間的互相換算公式,文獻[4]已做過推導。大地坐標與平面坐標的轉(zhuǎn)換關(guān)系可通過高斯投影正反算公式進行轉(zhuǎn)換,文獻[5]詳細描述并進行應(yīng)用。而空間直角坐標與平面直角坐標之間的轉(zhuǎn)換通常經(jīng)大地坐標過渡而間接進行轉(zhuǎn)換。
針對不同地球橢球基準下的大地坐標系、空間直角坐標系、平面直角坐標系也能進行相互轉(zhuǎn)換,但一般在不同基準下的同一坐標系統(tǒng)下進行轉(zhuǎn)換。不同基準下的空間直角轉(zhuǎn)換模型較常見的主要有布爾沙模型、莫洛金斯基模型和武測模型[6],三種模型均有七個參數(shù)(包括三平移參數(shù)、三旋轉(zhuǎn)參數(shù)和一個尺度參數(shù))。一般在精度要求不高,較小測區(qū)范圍內(nèi)可用三平移參數(shù)代替七參數(shù)進行轉(zhuǎn)換。
不同基準下的大地坐標系轉(zhuǎn)換通常包括廣義大地坐標微分公式轉(zhuǎn)換模型、基于平面直角坐標的坐標轉(zhuǎn)換模型。前者通常是由已知的公共點求其轉(zhuǎn)換參數(shù)或通過已知參數(shù)實現(xiàn)不同大地坐標系之間的直接轉(zhuǎn)換。后者一般是將大地坐標系轉(zhuǎn)換為同一基準下的空間直角坐標,然后將求出的空間直角坐標作為過渡,進行不同空間直角坐標之間的轉(zhuǎn)換,最后再將轉(zhuǎn)換后的新基準下的空間直角坐標轉(zhuǎn)換為大地坐標,從而間接的實現(xiàn)大地坐標系間的轉(zhuǎn)換。
廣義大地坐標微分公式轉(zhuǎn)換模型[7]的轉(zhuǎn)換參數(shù)有九個,除進行不同基準大地坐標轉(zhuǎn)換所需的七參數(shù)外還需增加轉(zhuǎn)換之間兩地球橢球長半軸和扁率的變化值二參數(shù),通常情況下這兩參數(shù)是已知的。廣義大地坐標微分公式轉(zhuǎn)換模型根據(jù)應(yīng)用一般又分為三維七參數(shù)模型、二維七參數(shù)模型[8]。三維七參數(shù)模型適用于全國及省級橢球面3°及以上不同地球橢球基準下的大地坐標系統(tǒng)間控制點坐標轉(zhuǎn)換。二維七參數(shù)模型適用于全國及省級適用于橢球面3°及以上不同地球橢球基準下的大地坐標系統(tǒng)間控制點坐標轉(zhuǎn)換。對于1954年北京坐標系、1980西安坐標系向2000國家大地坐標系的轉(zhuǎn)換,由于這兩個參心系下的大地高的精度較低,建議采用二維七參數(shù)轉(zhuǎn)換模型。
基于平面直角坐標的坐標轉(zhuǎn)換模型又包含相似變換模型、正形轉(zhuǎn)換模型和多項式轉(zhuǎn)換模型。其中相似轉(zhuǎn)換模型(簡稱二維四參數(shù)模型)包含二平移參數(shù)、一旋轉(zhuǎn)參數(shù)和一尺度參數(shù)共四參數(shù),適用于小范圍內(nèi)不同高斯投影平面坐標系間的轉(zhuǎn)換,沒有考慮兩坐標系間的局部變形和誤差積累,轉(zhuǎn)換精度較低。相似轉(zhuǎn)換模型可適用于省級以下小范圍內(nèi)控制點平面坐標轉(zhuǎn)換、相對獨立的平面坐標系統(tǒng)與2000國家大地坐標系的聯(lián)系。正形轉(zhuǎn)換模型[9]依據(jù)的是轉(zhuǎn)換過程中小范圍內(nèi)圖形保持相似,長度與方向脫離關(guān)系,而僅與點的位置有關(guān)。正形轉(zhuǎn)換相比相似轉(zhuǎn)換模型精度要高,適用于大范圍的坐標轉(zhuǎn)換,但是同樣只能用于高斯投影坐標系之間的轉(zhuǎn)換,不能解決投影變形。多項式轉(zhuǎn)換模型[10]是一種采用平面坐標位置作為因變量的二維坐標變換方法,如果轉(zhuǎn)換區(qū)域較大,公共點較多,可選擇更高階參數(shù)的多項式來提高精度。該轉(zhuǎn)換模型坐標轉(zhuǎn)換精度高,但轉(zhuǎn)換后的精度需進行檢核,適用于全國或大范圍的擬合,多用于相對獨立的平面坐標系統(tǒng)轉(zhuǎn)換。
2.3 坐標轉(zhuǎn)換方法
天津地方坐標系向CGCS2000大地坐標轉(zhuǎn)換方法步驟[11]如下:
(1)重合點選取。選用天津地方坐標系和CGCS2000國家大地坐標系坐標的控制點作為重合點。重合點選取的基本原則是:等級高、精度高、分布均勻、覆蓋整個轉(zhuǎn)換區(qū)域、局部變形小。
(2)轉(zhuǎn)換數(shù)學模型確定。天津地方坐標系作為北京54坐標系的局部投影坐標系,其參心坐標系的大地高通常是未知或不能精確獲取的,無法使用三維坐標轉(zhuǎn)換模型。另大地坐標向二維平面進行高斯投影過程中發(fā)生的變形在天津市域范圍內(nèi)不能忽略的,所以適合于二維七參數(shù)模型。
(3)轉(zhuǎn)換參數(shù)計算及剔除粗差點。首先需要將天津地方坐標系的控制點反投影成北京54坐標系所對應(yīng)的大地坐標點,然后將重合點代入模型利用最小二乘法計算轉(zhuǎn)換參數(shù)。用得到的轉(zhuǎn)換參數(shù)計算重合點坐標殘差,剔除殘差大于3倍點位中誤差的重合點,然后再代入模型計算直至重合點坐標殘差均小于3倍點位中誤差。最終用于計算轉(zhuǎn)換參數(shù)的重合點數(shù)量與轉(zhuǎn)換區(qū)域大小有關(guān),但不得少于6個。
(4)外部檢驗。選擇部分重合點作為外部檢核點,不參與轉(zhuǎn)換參數(shù)計算,用轉(zhuǎn)換參數(shù)計算這些點的轉(zhuǎn)換坐標與已知坐標進行比較進行外部檢核。應(yīng)選定至少6個均勻分布的重合點對坐標轉(zhuǎn)換精度進行檢核。
第一次全國國情普查項目[12]中需要將大量的天津90坐標ArcGIS Geodatabase格式的測繪成果轉(zhuǎn)換至國家要求的CGCS2000坐標系。然而本人單位使用的西安坐標轉(zhuǎn)換軟件(簡稱“西安軟件”)不支持ArgGIS 10.X高版本的GDB文件格式,且轉(zhuǎn)換速度很慢。因此開發(fā)了支持ArcGIS GDB文件格式的高效率坐標轉(zhuǎn)換軟件(簡稱“測繪院軟件”)。
本文的實驗是在Vistual Studio 2010環(huán)境下,用C#和ArcObjects10.2.2作為基本的開發(fā)語言和平臺,采用組件開發(fā)技術(shù),實現(xiàn)GeoDatabase結(jié)構(gòu)下的數(shù)據(jù)從天津90坐標系到國家GCCS2000坐標的轉(zhuǎn)換。
GeoDatabase結(jié)構(gòu)下的坐標轉(zhuǎn)換并不是簡單的點到點的轉(zhuǎn)換,雖然通過2.2節(jié)方法最后能過得到點轉(zhuǎn)換點的坐標,但基于GeoDatabase的轉(zhuǎn)換需充分考慮其本身的數(shù)據(jù)結(jié)構(gòu)、要素幾何模型等。基本思路是首先創(chuàng)建一個遍歷GeoDatabase數(shù)據(jù)結(jié)構(gòu)中的要素類,然后創(chuàng)建一個與轉(zhuǎn)換前要素類相同的Schema,但其空間參考需更改為CGCS2000。然后判斷要素類的幾何類型。如果是點則遍歷點要素直接進行轉(zhuǎn)換。如果是折線,則遍歷折線要素獲取要素幾何進而得到路徑Path集合。然后遍歷Path獲取曲線Segment集合,判斷曲線Segment類型并通過轉(zhuǎn)換后的點重新構(gòu)造生成新的Segment,然后逆向最終得到轉(zhuǎn)換后的要素幾何。多邊形的轉(zhuǎn)換需充分考慮多邊形的內(nèi)部Ring和外部Ring。多邊形轉(zhuǎn)換首先應(yīng)遍歷多邊形要素獲取多邊形的幾何包GeometryBag進而通過遞歸獲取幾何包的內(nèi)部Ring和外部Ring,然后遍歷內(nèi)部Ring和外部Ring同樣獲取曲線Segment幾何,判斷曲線Segment類型并通過轉(zhuǎn)換后的點重新構(gòu)造生成新的Segment,只不過如果是外部Ring的時候,在逆向通過Ring構(gòu)造多邊形幾何的時候需要反轉(zhuǎn)外部Ring的方向才能保證幾何一致。后再將重新轉(zhuǎn)換后的幾何和屬性創(chuàng)建要素,然后將新的要素寫入轉(zhuǎn)換后的要素類中,最后就生成了轉(zhuǎn)換為CGCS2000坐標系統(tǒng)的GeoDatabase數(shù)據(jù)。
實驗數(shù)據(jù)為兩組:第一組數(shù)據(jù)為國情普查[12]中2014年底上交國家局的7幅圖中的地表覆蓋和國情要素成果;第二組數(shù)據(jù)為天津全市域地表覆蓋成果中的旱地、高密度低矮房屋建筑區(qū)和低密度低矮房屋建筑區(qū)成果,GDB文件大小為413 M,包含 448 643個圖斑,約 3 952萬個頂點。第二組測試數(shù)據(jù)分布如圖1所示,范圍覆蓋天津全市域。
圖1 第二組測試數(shù)據(jù)分布范圍
第一組測試數(shù)據(jù)格式轉(zhuǎn)換如圖2所示。
圖2 第一組測試數(shù)據(jù)轉(zhuǎn)換流程
第二組測試數(shù)據(jù)格式轉(zhuǎn)換如圖3所示。
圖3 第二組測試數(shù)據(jù)轉(zhuǎn)換流程
將兩組測試原始數(shù)據(jù)分別用圖2、圖3所示流程轉(zhuǎn)換至CGCS2000坐標系GeoDatabase文件,后對比分析測繪院軟件和西安軟件轉(zhuǎn)換的CGCS2000的GeoDatabase文件,驗證測繪院轉(zhuǎn)換軟件的圖形一致性、線性方向一致性、屬性一致性、坐標轉(zhuǎn)換精度。同時對比測繪院軟件與西安軟件轉(zhuǎn)換同一組數(shù)據(jù)所耗時間,評估測繪院轉(zhuǎn)換軟件效率。
轉(zhuǎn)換圖形一致性測試采用第一組數(shù)據(jù)中的覆蓋LCA層、HYDA層、LCTL層和SFCP層,用ArcMap軟件的按位置選擇(菜單“選擇→按屬性選擇→與源圖層要素完全相同”),驗證測繪院軟件與西安軟件轉(zhuǎn)換后的GDB文件的幾何是否一一匹配。
實驗結(jié)果表明:第一組測試數(shù)據(jù)測繪院軟件轉(zhuǎn)換結(jié)果與西安軟件轉(zhuǎn)換結(jié)果圖形的幾何都是一一對應(yīng)。如覆蓋LCA層的一致性測試結(jié)果如圖4示,其他圖層測試結(jié)果一致,不再一一列出。
圖4 覆蓋LCA層轉(zhuǎn)換一致性測試結(jié)果
轉(zhuǎn)換圖形一致性測試采用第一組數(shù)據(jù)中的HYDL層,用ArcMap軟件顯示對比測繪院軟件與西安軟件轉(zhuǎn)換后的Geodatabase文件的線起點到終點的方向是否一致。
實驗結(jié)果表明:第一組測試數(shù)據(jù)測繪院軟件轉(zhuǎn)換結(jié)果與西安軟件轉(zhuǎn)換結(jié)果圖形的線方向一致。如HYDL層的測試結(jié)果如圖5示。
圖5 HYDL層線方向一致性測試結(jié)果
轉(zhuǎn)換屬性一致性測試采用第一組測試數(shù)據(jù)中的LCA層,用DataReviewer軟件中重復幾何校驗工具,然后選擇“校驗屬性”可進行重復幾何中的屬性比較。
實驗結(jié)果表明:第一組測試數(shù)據(jù)測繪院軟件轉(zhuǎn)換結(jié)果與西安軟件轉(zhuǎn)換結(jié)果屬性一致。如LCA層的測試結(jié)果如圖6示。要素總共個數(shù) 16 776,重復幾何且屬性相等的要素個數(shù)也為 16 776個。
圖6 LCA層重復幾何且重復屬性比較
坐標轉(zhuǎn)換精度測試采用了第一組數(shù)據(jù)的HYDA層和第二組數(shù)據(jù)中的LCA層,首先使用ArcGIS數(shù)據(jù)管理工具箱→要素→要素折點轉(zhuǎn)點工具,提取面或線的所有頂點得到面或線對應(yīng)的點層,然后通過數(shù)據(jù)管理工具箱→要素→添加XY坐標,將點層的X和Y坐標分別作為POINT_X和POINT_Y字段的形式添加到圖層屬性,然后采用工具箱→分析工具→疊加分析→空間連接,將兩組點層進行空間連接,生成一一對應(yīng)的一組數(shù)據(jù)。 然后通過字段計算器將點層中的西安軟件生成的XY坐標與測繪院軟件生成的XY坐標做差處理。
對第一組數(shù)據(jù):選取了HYDA層進行精度統(tǒng)計計算,其中測繪院軟件與西安軟件轉(zhuǎn)換結(jié)果x方向(經(jīng)度方向)和y方向(緯度方向)的最大差值分別為9.0e-009°≈0.000 97 m和7.0e-009°≈0.000 76 m。
對第二組數(shù)據(jù):選取低密度低矮房屋建筑區(qū)、高密度低矮房屋建筑區(qū)和旱地成果共 3 952萬個定點進行轉(zhuǎn)換,轉(zhuǎn)換如表1所示。
數(shù)據(jù)比較結(jié)果 表1
第二組數(shù)據(jù)出現(xiàn)了e-008級別上的偏差,進一步發(fā)現(xiàn)這些點對應(yīng)的多邊形存在自相交的問題。
轉(zhuǎn)換效率比較是分別用測繪院軟件和西安軟件對測試數(shù)據(jù)進行坐標轉(zhuǎn)換,通過對比轉(zhuǎn)換時間評估測繪院軟件坐標轉(zhuǎn)換效率。第二組測試數(shù)據(jù)包含 448 643個圖斑,約 3 952萬個頂點,西安軟件完成坐標轉(zhuǎn)換用時約 48 h,測繪院軟件完成轉(zhuǎn)換用時 42 min。
對兩組測試數(shù)據(jù)的轉(zhuǎn)換分析顯示:測繪院軟件轉(zhuǎn)換結(jié)果與西安軟件轉(zhuǎn)換結(jié)果圖形一致,線方向一致,屬性一致,轉(zhuǎn)換誤差在允許范圍內(nèi),測繪院轉(zhuǎn)換效率高于西安軟件轉(zhuǎn)換效率。
本文總結(jié)了坐標轉(zhuǎn)換相關(guān)的理論和方法,結(jié)合天津地方坐標的特點,選用合適的坐標轉(zhuǎn)換方法,并結(jié)合測繪成果所使用的GeoDatabase數(shù)據(jù)模型,實現(xiàn)了天津地方坐標系向CGCS2000坐標轉(zhuǎn)換。本文還以GIS的角度分別從圖形一致性、線性方向一致性、屬性一致性、坐標轉(zhuǎn)換精度和坐標轉(zhuǎn)換效率進行比較,實驗證明本文的轉(zhuǎn)換結(jié)果有效可行。
[1] 立實. 2000國家大地坐標系7月1日啟用[N]. 中國測繪報,2008-06-27(1).
[2] 吳儂,過靜珺. 現(xiàn)代普通測量學[M]. 北京:清華大學出版社,2001.
[3] 黃謨濤,翟國君,管錚等. 空間直角坐標和大地坐標的轉(zhuǎn)換[J]. 解放軍測繪學院學報,1998,15(3):164~167.
[4] 阮興龍,范東明. 高精度高斯投影正反算公式及其應(yīng)用[J]. 城市建設(shè)理論研究·電子版,2011(26).
[5] 于彩霞,黃文騫,樊沛. Bursa的3參數(shù)模型與7參數(shù)模型的適用性研究[J]. 測繪科學,2008,33(2):96~97.
[6] 武繼軍. 不同大地坐標系間坐標轉(zhuǎn)換模型研究[J]. 河南理工大學學報:自然科學版,2006,25(5):383~385.
[7] 馮里濤,鄧云青. 基于二維七參數(shù)轉(zhuǎn)換模型的坐標轉(zhuǎn)換參數(shù)的計算[J]. 城市勘測,2014(5):108~110.
[8] 張龍. 基于改進模型的坐標轉(zhuǎn)換精度分析[J]. 地礦測繪,2013,29(1):17~19.
[9] 鮑建寬,李永利,李秀海. 大地坐標轉(zhuǎn)換模型及其應(yīng)用[J]. 測繪工程,2013,22(3):56~60.
[10] 中國測繪科學研究院. 現(xiàn)有測繪成果轉(zhuǎn)換到2000國家大地坐標系技術(shù)指南[EB/OL].
Transformation of Tianjin Local Coordinate System to CGCS2000 Coordinate System Based on GeoDatabase
Wang Lin1,Yan Wei1,Chen Yanyun2
(1.Tianjin Institute of Surveying and Mapping,Tianjin 300381,China;2.National Marine Data and Information Service,Tianjin 300171,China)
GIS Data which is one of the most important results of surveying and mapping is often stored and managed by the geodatabase data model. In order to avoid duplication of surveying and mapping,We need to transform the GIS data of GeoDatabase structure from the local coordinate system of Tianjin to the CGCS2000 coordinate system.First of All,This article describes the GeoDatabase data model and summarizes the theory of coordinate transformation.Secendly,The article finds out the method and the steps of transforming the local coordinate system of Tianjin to the CGCS2000 coordinate system.Then,We realize the GIS data of GeoDatabase structure transforming from the local coordinate system of Tianjin to the CGCS2000 coordinate system with the experiment.Finally,it is proved effective and feasible from the graphical consistency,the linear direction consistency,the attribute consistency,the coordinate conversion precision and the coordinate conversion efficiency.
coordinate transformation;local coordinate system of Tianjin;the CGCS2000 coordinate system
1672-8262(2016)05-114-06
P226.3
B
2016—04—12
汪林(1985—),男,碩士,工程師,主要從事城市測繪與地理信息系統(tǒng)開發(fā)技術(shù)工作。