不同坐標(biāo)系的轉(zhuǎn)換及其程序設(shè)計(jì)畢業(yè)

Goodisgood,butbettercarriesit.

精益求精，善益求善。Goodisgood,butbettercarriesit.

精益求精，善益求善。不同坐標(biāo)系的轉(zhuǎn)換及其程序設(shè)計(jì)畢業(yè)第一章緒論第一章緒論摘要隨著空間技術(shù)的發(fā)展，全球一體化的形成，越來越多的要求全球測繪資料形成統(tǒng)一規(guī)范，尤其是坐標(biāo)系統(tǒng)的統(tǒng)一。由于各測量單位工作目的不同，所選擇的橢球參考系也會有所不同，出現(xiàn)了許多不同形式的坐標(biāo)系，例如WGS-84坐標(biāo)系、國家80坐標(biāo)系、北京54坐標(biāo)系、獨(dú)立地方坐標(biāo)及各種坐標(biāo)。在同一坐標(biāo)系下坐標(biāo)的表示方式又有空間直角坐標(biāo)、大地坐標(biāo)、平面坐標(biāo)。根據(jù)不同的測繪需求，需要將不同的坐標(biāo)系下的坐標(biāo)進(jìn)行相互轉(zhuǎn)換，在這些坐標(biāo)轉(zhuǎn)換的過程中既會運(yùn)用到同一坐標(biāo)系下的坐標(biāo)轉(zhuǎn)換模型，又會用到不同參考系下各坐標(biāo)系間的坐標(biāo)轉(zhuǎn)換模型。首先本文介紹大地測量學(xué)坐標(biāo)的相關(guān)知識，接著介紹了與坐標(biāo)轉(zhuǎn)換相關(guān)的知識以及坐標(biāo)轉(zhuǎn)換模型（包括同一坐標(biāo)系下的坐標(biāo)轉(zhuǎn)換模型和不同參考系下各坐標(biāo)系間的坐標(biāo)轉(zhuǎn)換模型），并利用vb語言實(shí)現(xiàn)坐標(biāo)轉(zhuǎn)換的過程。關(guān)鍵詞：地球橢球，坐標(biāo)系，轉(zhuǎn)換模型，坐標(biāo)轉(zhuǎn)換AbstractAlongwiththedevelopmentofspacetechnology,theformationofglobalintegration,moreandmorerequirementsofsurveyingandmappingmaterialformaunifiedglobalstandard,especiallytheunityofthecoordinatesystem.Becauseeachmeasurementunitworkpurposeisdifferent,choosetheframeofreferenceellipsoidwoulddiffer,theemergenceofmanydifferentformsofcoordinatesystem,suchasWGS-84coordinatesystem,thestate80coordinatesystem,Beijing54coordinatesystem,independentlocalcoordinatesystemandvariouskindsofurbanconstructioncoordinates.Inthesamecoordinatesystemofrepresentationandcoordinatespaceright-anglecoordinate,coordinate,coordinatetheearthplane.Accordingtothedifferentneedsofsurveyingandmapping,needdifferentcoordinatetransformationcoordinatesystem,inwhichtheprocessofcoordinatetransformationcanusetothesamecoordinatesystemcoordinatetransformationmodel,andwillusedifferentreferenceframe,thecoordinatetransformationbetweenthecoordinatesystemmodel.Firstofallthispaperintroducesthegeodeticcoordinatesofrelevantknowledge,thenintroducetheknowledgeofcoordinatetransformationandcoordinatetransformationmodel(includingthesamecoordinatesystemcoordinatetransformationmodelanddifferentreferenceframe,thecoordinatetransformationbetweenthecoordinatesystemmodel),andtheuseofvblanguagerealizationofcoordinatetransformationprocess.Keywords:theearthellipsoid,coordinatesystem,transformationmodel,coordinatetransformationTOC\o"1-3"\h\u摘要 1Abstract 21緒論 31.1研究的背景和意義 31.2國內(nèi)外研究現(xiàn)狀 41.3研究的主要內(nèi)容 52相關(guān)理論和知識介紹 52.1地球橢球 62.2基準(zhǔn) 72.3建立大地坐標(biāo)系的基本原理 72.3.1橢球定位、定向的概念 72.3.2參考橢球定位與定向的實(shí)現(xiàn)方法 82.3.3多點(diǎn)定位 83坐標(biāo)系統(tǒng)簡介 93.1測量常用的坐標(biāo)系 93.1.1大地坐標(biāo)系 93.1.2空間直角坐標(biāo)系 103.1.3平面坐標(biāo)系 113.1.4地方獨(dú)立坐標(biāo)系 123.2我國常用的坐標(biāo)系統(tǒng) 123.2.11954年北京坐標(biāo)系 123.2.21980年國家大地坐標(biāo)系 133.2.31954年北京坐標(biāo)系（整體平差轉(zhuǎn)換值） 143.2.4WGS-84世界大地坐標(biāo)系 143.2.52000國家大地坐標(biāo)系 154坐標(biāo)轉(zhuǎn)換理論與程序設(shè)計(jì) 164.1坐標(biāo)轉(zhuǎn)換的基本概念 164.2坐標(biāo)系轉(zhuǎn)換的模型 164.2.1同一參考橢球下大地坐標(biāo)系與空間直角坐標(biāo)系的相互轉(zhuǎn)換 174.3基準(zhǔn)轉(zhuǎn)換的模型 194.3.1不同地球橢球坐標(biāo)系的空間三參數(shù)或七參數(shù)轉(zhuǎn)換 195全文總結(jié) 23參考文獻(xiàn) 24附錄 25致謝 351緒論1.1研究的背景和意義隨著大地測量學(xué)，衛(wèi)星大地測量學(xué)，攝影測量學(xué)的發(fā)展和電子計(jì)算機(jī)的普及，對各種坐標(biāo)系的研究變得越來越重要了。精確地測量，計(jì)算和表示點(diǎn)的坐標(biāo)，為各種比例尺地形圖和大型工程測量提供控制，大地坐標(biāo)系作為大地測量基準(zhǔn)的一部分，一直是大地測量中最基本的問題。按其原點(diǎn)相對地球質(zhì)心的位置，大地坐標(biāo)系可為局部坐標(biāo)系和地心坐標(biāo)系。過去由于科技水平的制約，人類不能精確地確定地心的位置，局部坐標(biāo)系無疑是國家和地區(qū)的惟一選擇。應(yīng)用傳統(tǒng)技術(shù)建立起來的參心坐標(biāo)系逐漸難以滿足測繪及相關(guān)行業(yè)發(fā)展的需求，甚至在有些應(yīng)用中完全失去了意義。單純采用目前參心、二維、低精度、靜態(tài)的大地坐標(biāo)系統(tǒng)和相應(yīng)的基礎(chǔ)設(shè)施作為中國現(xiàn)行應(yīng)用的測繪基準(zhǔn)，必然會帶來越來越多的不協(xié)調(diào)問題新形勢下，測量坐標(biāo)系問題顯得越來越突出，使用地心坐標(biāo)系的要求也越來越迫切。世界許多發(fā)達(dá)國家和地區(qū)都開始采用地心坐標(biāo)。信息時(shí)代的控制測量儀器和測量系統(tǒng)已形成數(shù)字化，智能化和集成化的新發(fā)展態(tài)勢，空間測量和地面測量儀器和測量系統(tǒng)出現(xiàn)互補(bǔ)共榮的新的發(fā)展格局；傳統(tǒng)的大地測量技術(shù)發(fā)生了質(zhì)的變化，傳統(tǒng)的測繪行業(yè)逐漸向地理信息化產(chǎn)業(yè)轉(zhuǎn)換，工作重點(diǎn)已經(jīng)由外業(yè)轉(zhuǎn)為內(nèi)業(yè)處理。在實(shí)際測量中由于經(jīng)濟(jì)條件和環(huán)境條件的的限制，測量工作者在選取坐標(biāo)系的時(shí)候往往會正對實(shí)際情況選著最實(shí)用的坐標(biāo)系，這就對軟件提出了新的要求，需要數(shù)據(jù)處理的時(shí)候進(jìn)行轉(zhuǎn)換。由于各種轉(zhuǎn)換模型的相繼推出，對我們測量工作者來說，了解這些轉(zhuǎn)換的原理和數(shù)據(jù)的處理的過程方法是必要的。1.2國內(nèi)外研究現(xiàn)狀自60年代以來,各國大地測量學(xué)者,經(jīng)過大量研究,提出了多種坐標(biāo)轉(zhuǎn)換模型及多種解算方法,北美1927基準(zhǔn)面(基于克拉克1966橢球體與北美1983基準(zhǔn)面(基于GRS1980橢球體)之間坐標(biāo)轉(zhuǎn)換是根據(jù)研究區(qū)內(nèi)一系列已知點(diǎn)的大地坐標(biāo)或網(wǎng)格坐標(biāo)改正量進(jìn)行插值進(jìn)行的坐標(biāo)系轉(zhuǎn)換；英國采用北向與東向的雙線性網(wǎng)格插值進(jìn)行坐標(biāo)轉(zhuǎn)換；挪威在海岸帶調(diào)查中，采用經(jīng)緯度多項(xiàng)式用于坐標(biāo)系轉(zhuǎn)換這種方法進(jìn)行新（ED87—?dú)W洲1987基準(zhǔn)面）、舊(ED50—?dú)W洲1950基準(zhǔn)面)坐標(biāo)系之間的轉(zhuǎn)換；歐洲石油勘探組織(EPSG)對新、舊坐標(biāo)系采用“雙線性插值”進(jìn)行坐標(biāo)轉(zhuǎn)換。在國內(nèi)空間三維直角坐標(biāo)轉(zhuǎn)換中,通常采用7參數(shù)布爾沙—沃爾夫模型、莫洛金斯基模型和范式模型，并且劉經(jīng)南院士和其同事在對這三種傳統(tǒng)轉(zhuǎn)換模型進(jìn)行分析的基礎(chǔ)上，從理論和實(shí)踐上證明了這三個(gè)模型的等價(jià)性，并在此基礎(chǔ)上他還提出了第4個(gè)等價(jià)模型—“武測模型”，這些模型雖然表示形式上略有差異，但從坐標(biāo)轉(zhuǎn)換的最終結(jié)果而言，他們是等價(jià)的。1.3研究的主要內(nèi)容測量坐標(biāo)轉(zhuǎn)換問題在測量工程中經(jīng)常遇到，其計(jì)算過程比較繁瑣，采用手動(dòng)計(jì)算式相當(dāng)麻煩，國內(nèi)的許多坐標(biāo)軟件都有一個(gè)缺點(diǎn)操作界面過于復(fù)雜，操作起來也很繁瑣，為了提高軟件的交互性和實(shí)用性，我采用VB實(shí)現(xiàn)這一目的。本文首先介紹研究的背景和意義，國內(nèi)外研究現(xiàn)狀，接著對國內(nèi)外有關(guān)空間三維直角坐標(biāo)系做了系統(tǒng)概述，接著介紹了與坐標(biāo)轉(zhuǎn)換相關(guān)的知識以及坐標(biāo)轉(zhuǎn)換模型，并用VB設(shè)計(jì)了相關(guān)程序，最后是全文總結(jié)。2相關(guān)理論和知識介紹由于當(dāng)今世界上有著有許多的參心坐標(biāo)系和地心坐標(biāo)系等，所以對于地球表面上的任一一點(diǎn)P，表述改點(diǎn)坐標(biāo)的方式有很多。因此對于地面上一點(diǎn)，由于所選擇的坐標(biāo)系不同，其表達(dá)方式也會不同。而且，即使使用同一坐標(biāo)系，也會有不同的表達(dá)方式。想要弄清楚它們之間的聯(lián)系，那么就會涉及到坐標(biāo)轉(zhuǎn)換的問題。坐標(biāo)系統(tǒng)之間的坐標(biāo)轉(zhuǎn)換既包括不同的參心坐標(biāo)之間的轉(zhuǎn)換，或者不同的地心坐標(biāo)系之間的轉(zhuǎn)換，也包括參心坐標(biāo)系與地心坐標(biāo)系之間的轉(zhuǎn)換以及同一坐標(biāo)系下的直角坐標(biāo)與大地坐標(biāo)之間的坐標(biāo)轉(zhuǎn)換，還有大地坐標(biāo)與高斯平面坐標(biāo)之間的坐標(biāo)轉(zhuǎn)換等。2.1地球橢球由于地球內(nèi)部質(zhì)量分布不均，導(dǎo)致大地體其實(shí)是一個(gè)不平的似球體，經(jīng)過長期的理論研究和實(shí)踐認(rèn)為，當(dāng)一個(gè)通過南北兩極的子午圈，繞地球南北極旋轉(zhuǎn)一周而形成一個(gè)橢球體，用這個(gè)橢球面來代替大地水準(zhǔn)面，是一個(gè)很理想的計(jì)算基準(zhǔn)面。一個(gè)與大地體符合最好，最接近地球大小和形狀的旋轉(zhuǎn)橢球，稱為總地球橢球體。其具體條件為：1.總地球橢球體的體積與大地體的體積一致，而且其表面與大地水準(zhǔn)面之間的差距的平方和最小。2.總地球橢球體的總質(zhì)量與地球的總質(zhì)量一致，而且其中心與地球重心相重合，總地球橢球的赤道面也應(yīng)該與地球的赤道一致。3.總地球橢球體的旋轉(zhuǎn)角速度與地球的旋轉(zhuǎn)角速度一致。在眾多橢球體中，WGS-84橢球體被認(rèn)為符合上述條件最好的橢球，由于經(jīng)典大地測量技術(shù)存在一定的局限性，大地測量工作者算出一個(gè)涵蓋整個(gè)大陸和海洋的總地球橢球，而是根據(jù)本國的測繪成果推求出一個(gè)最能表達(dá)本國或者是本地區(qū)的地球橢球體，就是所謂的參考橢球。地球表面、大地水準(zhǔn)面和橢球面三者的關(guān)系及偏差如下圖所示圖2-1地球三個(gè)面及其偏差2.2基準(zhǔn)所謂基準(zhǔn)是指為描述空間位置而定義的點(diǎn)線面。而大地測量基準(zhǔn)是指用以描述地球形狀的地球橢球參數(shù)，包含描述地球橢球幾何特征的長短半軸和物理特征的有關(guān)參數(shù)、地球在空間的定位及定向以及描述這些位置所采用的單位長度的定義。不同的坐標(biāo)系統(tǒng)會使用的基準(zhǔn)也不同。在大地測量中，根據(jù)參考橢球所選原點(diǎn)位置不同，可以分為地心坐標(biāo)系和參心坐標(biāo)系。地心坐標(biāo)系是以地球的質(zhì)心為原點(diǎn)，同樣有地心大地坐標(biāo)系和地心空間直角坐標(biāo)系兩種表述方法。地心空間直角坐標(biāo)系的定義為:以地球質(zhì)心為原點(diǎn)，X軸指向格林尼治子午面與地球赤道的交點(diǎn)，Z軸指向北極，Y軸過原點(diǎn)垂直于平面XOZ，構(gòu)成右手空間直角坐標(biāo)系。地心大地坐標(biāo)系定義為：以地球的質(zhì)心作為原點(diǎn)，以地球自轉(zhuǎn)軸作為橢球的短軸，大地緯度B是過地面點(diǎn)的橢球法線與橢球赤道面之間的夾角，大地經(jīng)度L為過地面點(diǎn)的橢球子午面與格林尼治子午面之間的夾角，大地高度H為地面點(diǎn)沿橢球法線到橢球面的最短距離。參心坐標(biāo)系是這樣定義的：選取一個(gè)參考橢球面作為基本的參考面，選一參考點(diǎn)作為大地測量的起算點(diǎn)，并且通過大地的質(zhì)點(diǎn)來進(jìn)行測量，從而確定參考橢球在地球面的位置和方向。這時(shí)參考橢球的原點(diǎn)一般不會和地球質(zhì)心重合，所以稱為參心。參心坐標(biāo)主要用于大地測量中，如測量某一地區(qū)的控制網(wǎng)等，所以又稱局部坐標(biāo)。它同樣具有參心大地坐標(biāo)系和參心直角坐標(biāo)系兩種表述方法，它們的定義與地心坐標(biāo)系的定義相似。2.3建立大地坐標(biāo)系的基本原理2.3.1橢球定位、定向的概念大地坐標(biāo)系是建立在一定的大地基準(zhǔn)上的用于表達(dá)地球表面空間位置及其相對關(guān)系的數(shù)學(xué)參照系，這里所說的大地基準(zhǔn)是指能夠最佳擬合地球形狀的地球橢球的參數(shù)及橢球定位和定向。橢球定位是確定橢球中心的位置,可分為兩類:局部定位和地心定位。局部定位要求在一定范圍內(nèi)橢球面與大地水準(zhǔn)面有最佳的符合,而對橢球的中心位置無特殊要求；地心定位要求在全球范圍內(nèi)橢球面與大地水準(zhǔn)面有最佳的符合,同時(shí)要求橢球中心與地球質(zhì)心一致或最為接近。橢球定向是指確定橢球旋轉(zhuǎn)軸的方向,不論是局部定位還是地心定位,都應(yīng)滿足兩個(gè)平行條件:①橢球短軸平行于地球自轉(zhuǎn)軸;②大地起始子午面平行于天文起始子午面具有確定參數(shù)(長半徑a和扁率α),經(jīng)過局部定位和定向,同某一地區(qū)大地水準(zhǔn)面最佳擬合的地球橢球,叫做參考橢球。除了滿足地心定位和雙平行條件外,在確定橢球參數(shù)時(shí)能使它在全球范圍內(nèi)與大地體最密合的地球橢球,叫做總地球橢球。2.3.2參考橢球定位與定向的實(shí)現(xiàn)方法建立（地球）參心坐標(biāo)系，需進(jìn)行下面幾個(gè)工作：①選擇或求定橢球的幾何參數(shù)（長短半徑）；②確定橢球中心位置（定位）；③確定橢球短軸的指向（定向）；④建立大地原點(diǎn)。參考橢球定位定向方法選定某一適宜的點(diǎn)K作為大地原點(diǎn)，在該點(diǎn)上實(shí)施精密的天文測量和高程測量，由此得到該點(diǎn)的天文經(jīng)度，天文緯度，至某一相鄰點(diǎn)的天文方位角和正高得到K點(diǎn)相應(yīng)的大地經(jīng)度，大地緯度，至某一相鄰點(diǎn)的大地方位角和大地高（2.1）2.3.3多點(diǎn)定位一點(diǎn)定位的結(jié)果在較大范圍內(nèi)往往難以使橢球面與大地水準(zhǔn)面有較好的密合。所以在國家或地區(qū)的天文大地測量工作進(jìn)行到一定的時(shí)候或基本完成后，利用許多拉普拉斯點(diǎn)（即測定了天文經(jīng)度、天文緯度和天文方位角的大地點(diǎn)）的測量成果和已有的橢球參數(shù)，按照廣義弧度測量方程按=最小或=最?。┻@一條件，通過計(jì)算進(jìn)行新的定位和定向，從而建立新的參心大地坐標(biāo)系。按這種方法進(jìn)行參考橢球的定位和定向，由于包含了許多拉普拉斯點(diǎn)，因此通常稱為多點(diǎn)定位法。多點(diǎn)定位的結(jié)果使橢球面在大地原點(diǎn)不再同大地水準(zhǔn)面相切，但在所使用的天文大地網(wǎng)資料的范圍內(nèi)，橢球面與大地水準(zhǔn)面有最佳的密合。3坐標(biāo)系統(tǒng)簡介所謂坐標(biāo)系指的是描述空間位置的表達(dá)形式,即采用什么方法來表示空間位置。人們?yōu)榱嗣枋隹臻g位置,采用了多種方法,從而也產(chǎn)生了不同的坐標(biāo)系。3.1測量常用的坐標(biāo)系3.1.1大地坐標(biāo)系空間大地坐標(biāo)系以大地經(jīng)度L，地緯度B，大地高H來表示空間某一點(diǎn)的位置。地面上P點(diǎn)的大地子午面NPS與起始大地子午面所構(gòu)成的二面角L，叫做P點(diǎn)的大地經(jīng)度，有起始子午面起算，向東為正，成為東經(jīng)，向西為負(fù)，稱為西經(jīng)，P點(diǎn)對于橢球的法線PK與赤道面的夾角B，叫做P點(diǎn)的大地緯度，有赤道面起算，向北為正，稱為北緯，向南為負(fù)，稱為南緯，如3-1圖所示圖3-1大地坐標(biāo)示意圖P點(diǎn)沿法線到橢球面的距離H，叫大地高，從橢球面起算，向外為正，向里為負(fù)。GPS測量出來的高程為大地高，與我們所選要的正常高存在高程異常。大地高H與水準(zhǔn)測量中的正常高或正高有以下關(guān)系H（大地高）=（正常高）+（高程異常）圖3-2大地水準(zhǔn)面的差距3.1.2空間直角坐標(biāo)系空間直角坐標(biāo)系的坐標(biāo)原點(diǎn)與參考橢球的中心重合，Z軸正向指向參考橢球的北極，X軸正向指向起始子午面與赤道的交點(diǎn)，Y軸按右手系與X軸呈90°夾角且位于赤道面上。某點(diǎn)在空間中的坐標(biāo)可用該點(diǎn)在此空間坐標(biāo)系的各個(gè)坐標(biāo)軸上的投影來表示，如圖所示:圖3-3空間直角坐標(biāo)系示意圖3.1.3平面坐標(biāo)系平面直角坐標(biāo)系是利用投影，將空間坐標(biāo)(空間直角坐標(biāo)或空間大地坐標(biāo))通過某種數(shù)學(xué)變換映射到平面上，這種變換稱為投影變換。投影變換的方法有很多，如Lambuda投影，UTM投影等，在我國一般采用的是高斯一克呂格投影，也稱為高斯投影。地形測圖以及許多的測量定位應(yīng)用在現(xiàn)實(shí)中是我們常見的平面直角坐標(biāo)。對于一個(gè)國家或較大區(qū)域，應(yīng)按照一定的數(shù)學(xué)法則將參考橢球面上的各點(diǎn)的大地經(jīng)緯度投影為平面上相對應(yīng)點(diǎn)的平面直角坐標(biāo)。由于地球橢球面是不可展曲面，所以無論采用什么樣的投影都會產(chǎn)生一定變形。投影變形一般分為長度變形，角度變形和面積變形這三種。根據(jù)制測量的任務(wù)和目的應(yīng)當(dāng)采用等角投影(又稱正形投影)。在采用的正形投影時(shí)，還要求長度和面積變形不大，并且能用簡單的公式來計(jì)算這些變形而帶來的改正數(shù)。為了解決這些的矛盾，測量上往往是將一個(gè)大的區(qū)域按照一定規(guī)律分成若干個(gè)小的區(qū)域(或帶)。每個(gè)區(qū)域單獨(dú)投影，并組成自身的直角坐標(biāo)系，然后，在將這些帶用簡單的數(shù)學(xué)方法聯(lián)系起來，從而組成同一系統(tǒng)。目前測量上廣泛采用的是高斯投影，它是一種正形投影，它的特點(diǎn)是：沒有角度變形，在不同點(diǎn)上的長度比隨點(diǎn)位而異，但在同一點(diǎn)上各方向的長度比相同。高斯-克呂格正形投影又稱橫軸橢圓柱投影，即橢圓柱內(nèi)面橫套在地球橢球的外表面，橢圓柱的中心通過橢球的中心，并在某一中央子午線上相切，該中央子午線就是高斯平面直角坐標(biāo)系的X軸，X軸沒有長度變形，赤道在橢圓柱上的投影是高斯平面直角坐標(biāo)系的Y軸，把橢球柱展開，就得到以(X，Y)為坐標(biāo)的高斯平面直角坐標(biāo)系。3.1.4地方獨(dú)立坐標(biāo)系在我國，平面坐標(biāo)主要采用的是高斯投影，在該投影中，除中央子午線外，其它位置上的任何線段，投影后都會產(chǎn)生一定的長度變形，而且變形隨離開中央子午線的距離增加而增加。因此，一般采用分帶投影的辦法，來限制長度變形。我國規(guī)定了采用3度帶或6度帶進(jìn)行分帶投影。在城市、工礦等工程測量中，如果直接在國家分帶坐標(biāo)系中建立控制網(wǎng)，會使地面長度投影的變形較大，當(dāng)長度變形大于2.5cm/km時(shí)，就難以滿足工程上的需要。因此為了滿足大比例尺測圖和進(jìn)行施工放樣時(shí)的需要，必需是基于一個(gè)與當(dāng)?shù)仄骄０胃叱虒?yīng)的參考橢球。而且該橢球的扁率、中心、軸向一般與國家參考橢球體相同，但其長半軸則必有一個(gè)修正量，這個(gè)參考橢球就被稱為地方參考橢球。另一些特殊的測量，比如大橋施工測量，水利水壩測量，滑坡變形監(jiān)測等，采用國家坐標(biāo)系精度達(dá)不到要求，不實(shí)用也不方便，常常會建立適合本地區(qū)的地方獨(dú)立坐標(biāo)系。3.2我國常用的坐標(biāo)系統(tǒng)3.2.11954年北京坐標(biāo)系1954年，總參測繪局在有關(guān)方面的建議與支持下，鑒于當(dāng)時(shí)的歷史條件，采取先將我國一等鎖與前蘇聯(lián)遠(yuǎn)東一等鎖相聯(lián)接，然后以連接處呼瑪，吉拉林，東寧基線網(wǎng)擴(kuò)大邊端點(diǎn)的前蘇聯(lián)1942年普爾科沃坐標(biāo)系的坐標(biāo)為起算數(shù)據(jù)，平差我國東北及東部一等鎖，這樣從蘇聯(lián)傳算來的坐標(biāo)系定名為1954年北京坐標(biāo)系。1954年北京坐標(biāo)系實(shí)際上是前蘇聯(lián)1942年普爾科沃坐標(biāo)系在我國的延伸，但我國坐標(biāo)系的大地點(diǎn)高程（1956年黃海高程系）卻與前蘇聯(lián)坐標(biāo)系的計(jì)算基準(zhǔn)面不同，因此嚴(yán)格意義上來說，二者不是完全相同的大地坐標(biāo)系。特點(diǎn):1954年北京坐標(biāo)系屬于參心坐標(biāo)系；采用克拉索夫斯基橢球參數(shù)；多點(diǎn)定位:垂線偏差由900個(gè)點(diǎn)解得，大地水準(zhǔn)面差距由43個(gè)點(diǎn)解得；大地原點(diǎn)是前蘇聯(lián)的普爾科沃；大地點(diǎn)高程是以1956年青島驗(yàn)潮站求出的黃海平均海水面為基準(zhǔn)；高程異常是以前蘇聯(lián)1955年大地水準(zhǔn)面重新平差結(jié)果為起算值，按我國天文水準(zhǔn)路線推算出來的；提供的大地點(diǎn)成果是局部平差結(jié)果。問題和缺點(diǎn)：克拉索夫斯基橢球比現(xiàn)代精確橢球相差過大；只涉及兩個(gè)幾何性質(zhì)的橢球參數(shù)（a和α），滿足不了當(dāng)今理論研究和實(shí)際工作中所需四個(gè)地球橢球基本參數(shù)的要求；處理重力數(shù)據(jù)時(shí)采用的是赫爾默特1901到1909年正常重力公式，與之相應(yīng)的赫爾默特扁球不是旋轉(zhuǎn)橢球，它與克拉索夫斯基橢球是不一致的；對應(yīng)的參考橢球面與我國大地水準(zhǔn)面存在著自西向東明顯的系統(tǒng)性傾斜，在東部地區(qū)高程異常最大達(dá)到＋65米，全國范圍平均29米；橢球定向不明確，橢球短軸指向既不是CIO,也不是我國的JYD1968.0；起始子午面不是國際時(shí)間局BIH所定義的格林尼治平均天文臺子午面，給坐標(biāo)換算帶來一些不便和誤差；坐標(biāo)系未經(jīng)整體平差而僅是局部平差成果，點(diǎn)位精度不高，也不均勻；名不副實(shí)，容易引起一些誤解。3.2.21980年國家大地坐標(biāo)系特點(diǎn)：1980年國家大地坐標(biāo)系屬參心大地坐標(biāo)系；采用既含幾何參數(shù)又含物理參數(shù)的四個(gè)橢球基本參數(shù)。數(shù)值采用1975年IUGG第16屆大會的推薦值；多點(diǎn)定位；定向明確。地球橢球短軸平行于由地球質(zhì)心指向地極原點(diǎn)JYD1968.0方向，起始大地子午面平行于我國起始天文子午面；大地原點(diǎn)在我國中部：陜西省涇陽縣永樂鎮(zhèn)，簡稱西安原點(diǎn)；大地點(diǎn)高程以1956年青島驗(yàn)潮站求出的黃海平均海水面為基準(zhǔn)；1980年國家大地坐標(biāo)系建立后，進(jìn)行了全國天文大地網(wǎng)整體平差，計(jì)算了5萬余個(gè)點(diǎn)的成果。新問題：原來的各種關(guān)于橢球參數(shù)的用表均要變更低等點(diǎn)要重新平差，編撰新的三角點(diǎn)成果表地形圖圖廓線和方里網(wǎng)線位置發(fā)生變化，并引起地形圖內(nèi)地形、地物相關(guān)位置的改變新形勢下1980年國家大地坐標(biāo)系的地極原點(diǎn)JYD1968.0已不能適應(yīng)當(dāng)代建立高精度天文地球動(dòng)力學(xué)系帶要求。3.2.31954年北京坐標(biāo)系（整體平差轉(zhuǎn)換值）它是在1980年國家大地坐標(biāo)系的基礎(chǔ)上，改變IUGG1975年橢球至克拉索夫斯基橢球，通過在空間三個(gè)坐標(biāo)軸上進(jìn)行平移而來的。因此，其坐標(biāo)值仍體現(xiàn)了整體平差的特點(diǎn)，精度和1980年國家大地坐標(biāo)系相同，克服了1954年北京坐標(biāo)系局部平差的缺點(diǎn)；其坐標(biāo)軸和1980年國家大地坐標(biāo)系坐標(biāo)軸相互平行，所以它的定向明確；它的橢球參數(shù)恢復(fù)為1954年北京坐標(biāo)系的橢球參數(shù)，從而使其坐標(biāo)值和1954年北京坐標(biāo)系局部平差坐標(biāo)值相差較小。特點(diǎn)：屬參心大地坐標(biāo)系；長短軸采用克拉索夫斯基橢球參數(shù)；多點(diǎn)定位，參心雖和1954年北京坐標(biāo)系參心不相一致，但十分接近；定向明確，與1980年國家大地坐標(biāo)系的定向相同；大地原點(diǎn)與1980年國家大地坐標(biāo)系相同，但大地起算數(shù)據(jù)不同；大地點(diǎn)高程基準(zhǔn)是以1956年青島驗(yàn)潮站求出的黃海平均海水面為基準(zhǔn)；提供坐標(biāo)是1980年國家大地坐標(biāo)系整體平差轉(zhuǎn)換值，精度一致；用于測圖坐標(biāo)系，對于1:5萬以下比例尺測圖，新舊圖接邊，不會產(chǎn)生明顯裂痕。3.2.4WGS-84世界大地坐標(biāo)系該坐標(biāo)系是一個(gè)協(xié)議地球參考系CTS（ConventionalTerrestrialSystem）,其原點(diǎn)是地球的質(zhì)心，Z軸指向BIH1984.0定義的協(xié)議地球極CTP（ConventionalTerrestrialPole）方向，X軸指向BIH1984.0零度子午面和CTP赤道的交點(diǎn)，Y軸和Z、X軸構(gòu)成右手坐標(biāo)系。WGS-84橢球采用國際大地測量與地球物理聯(lián)合會第17屆大會大地測量常數(shù)推薦值自1987年1月10日之后，GPS衛(wèi)星星歷均采用WGS-84坐標(biāo)系統(tǒng)。因此GPS網(wǎng)的測站坐標(biāo)及測站之間的坐標(biāo)差均屬于WGS-84系統(tǒng)。為了求得GPS測站點(diǎn)在地面坐標(biāo)系（屬于參心坐標(biāo)系）中的坐標(biāo)，就必須進(jìn)行坐標(biāo)系的轉(zhuǎn)換。3.2.52000國家大地坐標(biāo)系2008年3月，由國土資源部正式上報(bào)國務(wù)院《關(guān)于中國采用2000國家大地坐標(biāo)系的請示》，并于2008年4月獲得國務(wù)院批準(zhǔn)。自2008年7月1日起，中國將全面啟用2000國家大地坐標(biāo)系，國家測繪局受權(quán)組織實(shí)施。2000國家大地坐標(biāo)系是全球地心坐標(biāo)系在我國的具體體現(xiàn)，其原點(diǎn)為包括海洋和大氣的整個(gè)地球的質(zhì)量中心。2000國家大地坐標(biāo)系采用的地球橢球參數(shù)如下：長半軸a=6378137m扁率f=1/298.257222101地心引力常數(shù)GM=3.986004418×1014m3s-2自轉(zhuǎn)角速度ω=7.292l15×10-5rads-1采用2000國家大地坐標(biāo)系具有科學(xué)意義，隨著經(jīng)濟(jì)發(fā)展和社會的進(jìn)步，我國航天、海洋、地震、氣象、水利、建設(shè)、規(guī)劃、地質(zhì)調(diào)查、國土資源管理等領(lǐng)域的科學(xué)研究需要一個(gè)以全球參考基準(zhǔn)為背景的、全國統(tǒng)一的、協(xié)調(diào)一致的坐標(biāo)系統(tǒng)，來處理國家、區(qū)域、海洋與全球化的資源、環(huán)境、社會和信息等問題，需要采用定義更加科學(xué)、原點(diǎn)位于地球質(zhì)量中心的三維國家大地坐標(biāo)系。4坐標(biāo)轉(zhuǎn)換理論與程序設(shè)計(jì)4.1坐標(biāo)轉(zhuǎn)換的基本概念坐標(biāo)轉(zhuǎn)換是測繪實(shí)踐中經(jīng)常遇到的重要問題之一。坐標(biāo)轉(zhuǎn)換通常包含兩層含義：坐標(biāo)系變換和基準(zhǔn)變換。坐標(biāo)系變換：就是在同一地球橢球下，空間點(diǎn)的不同坐標(biāo)表示形式間進(jìn)行變換。包括大地坐標(biāo)系與空間直角坐標(biāo)系的相互轉(zhuǎn)換、空間直角坐標(biāo)系與站心坐標(biāo)系的轉(zhuǎn)換、以及大地坐標(biāo)系與高斯平面坐標(biāo)系的轉(zhuǎn)換（即高斯投影正反算）基準(zhǔn)變換：是指空間點(diǎn)在不同的地球橢球間的坐標(biāo)變換?？捎每臻g的三參數(shù)或七參數(shù)實(shí)現(xiàn)不同橢球間空間直角坐標(biāo)系或不同橢球間大地坐標(biāo)系的轉(zhuǎn)換。4.2坐標(biāo)系轉(zhuǎn)換的模型4.2.1同一參考橢球下大地坐標(biāo)系與空間直角坐標(biāo)系的相互轉(zhuǎn)換（1）大地坐標(biāo)系轉(zhuǎn)換為空間直角坐標(biāo)系（BLH→XYZ）將同一坐標(biāo)系下的大地坐標(biāo)(B,L,H)轉(zhuǎn)換成空間直角坐標(biāo)(X,Y,Z)的轉(zhuǎn)換公式為：（4.1）式中，e為第一偏心率；N為卯酉圈的半徑；b為短半軸；a為參考橢球長半軸；b為短半軸；并且有若點(diǎn)在橢球面上，則大地高H=0，式可簡化為：（4.2）（2）空間直角坐標(biāo)系轉(zhuǎn)換為大地坐標(biāo)系（XYZ→BLH）將同一坐標(biāo)系下的空間直角坐標(biāo)(X,Y,Z}轉(zhuǎn)換為大地坐標(biāo)(B,L,H)的公式為：（4.3）在使用上式進(jìn)行空間直角坐標(biāo)到大地坐標(biāo)的轉(zhuǎn)換時(shí)，因?yàn)橛?jì)算大地緯度B時(shí)需要用到大地高H，而計(jì)算大地高時(shí)又需要用到大地緯度B。因此不能直接計(jì)算出大地坐標(biāo)，而需要采用迭代計(jì)算的方法。具體計(jì)算時(shí)，可先根據(jù)下式求出大地緯度B的初值:然后利用該初值代入公式來求大地高H,N的初值，再利用所求出的大地高H和N的初值代入公式中再次求出B值。再將B代入公式求H和N,如此反復(fù)，直至求出的B,H,N收斂為止。也可以采用下面的算法直接將空間直角坐標(biāo)轉(zhuǎn)換為大地坐標(biāo):（4.4）e為參考橢球的第二偏心率。在主要代碼如下：各橢球參數(shù)：PI=3.14159265358979XAa=6378140:XAb=6356755.28815753:XAc=6399596.65198801:XAe2=0.006694384999588:XAe12=0.006739501819473BJa=6378245:BJb=6356863.01877305:BJc=6399698.90178271:BJe2=0.006693421622966:BJe12=0.006738525414683WGSa=6378137:WGSb=6356752.3142:WGSc=6399593.6258:WGSe2=0.0066943799013:WGSe12=0.00673949674227度分秒轉(zhuǎn)換下弧度：s=Val(Txt)du=Fix(s)fe=Fix((s-du)*100)mi=Round((s-du-fe*0.01)*10000,4)Rad=(du+fe/60+mi/3600)*PI/180程序設(shè)計(jì)界面如下：每一個(gè)國家大地坐標(biāo)系都有自己的基準(zhǔn)參數(shù)，每一個(gè)坐標(biāo)系都有幾種表現(xiàn)形式，如上圖所示，可以用B、L、H表示也可以用X、Y、Z所示，只要知道橢球參數(shù)利用公式即可實(shí)現(xiàn)這轉(zhuǎn)換。利用VB實(shí)現(xiàn)大地坐標(biāo)與空間直角坐標(biāo)的轉(zhuǎn)換程序直觀易懂，操作方便，在日常的測繪工作中這種類似程序并不經(jīng)常，程序的缺陷沒有把經(jīng)緯度轉(zhuǎn)換成度分秒。4.3基準(zhǔn)轉(zhuǎn)換的模型4.3.1不同地球橢球坐標(biāo)系的空間三參數(shù)或七參數(shù)轉(zhuǎn)換不同地球橢球之間的坐標(biāo)系轉(zhuǎn)換實(shí)際上是不同基準(zhǔn)之間的轉(zhuǎn)換。不同基準(zhǔn)之間的轉(zhuǎn)換方法很多，可以通過空間變換的方法實(shí)現(xiàn)，亦可用平面變換方法進(jìn)行。下面介紹七參數(shù)布爾莎模型設(shè)兩不同地球橢球的對應(yīng)的兩個(gè)空間直角坐標(biāo)系間有7個(gè)轉(zhuǎn)換參數(shù)：3個(gè)平移參數(shù)（原點(diǎn)不重合產(chǎn)生）；3個(gè)旋轉(zhuǎn)參數(shù)（坐標(biāo)軸不平行產(chǎn)生）；1個(gè)尺度參數(shù)（兩坐標(biāo)系間的尺度不一致產(chǎn)生）。見下圖設(shè)（XA，YA，ZA）為某點(diǎn)在A空間直角坐標(biāo)系中的三維坐標(biāo)，（XB，YB，ZB）為某點(diǎn)在B空間直角坐標(biāo)系中的三維坐標(biāo)，（△X0，△Y0，△Z0）為某點(diǎn)從A空間直角坐標(biāo)系轉(zhuǎn)換到B空間直角坐標(biāo)系中的三個(gè)平移參數(shù)，（ωX，ωY，ωZ）為某點(diǎn)從A空間直角坐標(biāo)系轉(zhuǎn)換到B空間直角坐標(biāo)系中的三個(gè)旋轉(zhuǎn)參數(shù)，K為某點(diǎn)從A空間直角坐標(biāo)系轉(zhuǎn)換到B空間直角坐標(biāo)系中的三個(gè)尺度參數(shù)。則點(diǎn)從A空間直角坐標(biāo)系轉(zhuǎn)換到B空間直角坐標(biāo)系中的模型為式中：程序設(shè)計(jì)界面如下：5全文總結(jié)橢球面上的點(diǎn)位可在各種坐標(biāo)系中表示，由于所用坐標(biāo)系不同，表現(xiàn)出來的坐標(biāo)值也不同。在實(shí)際應(yīng)用中，要想對不同測量坐標(biāo)系下的坐標(biāo)進(jìn)行綜合處理，使他們之間有所聯(lián)系，必然會涉及到測量坐標(biāo)轉(zhuǎn)換問題。測量坐標(biāo)轉(zhuǎn)換包含測量坐標(biāo)系轉(zhuǎn)換和測量坐標(biāo)基準(zhǔn)轉(zhuǎn)換兩個(gè)方面內(nèi)容。測量坐標(biāo)系轉(zhuǎn)換是指在同一參考橢球下，空間點(diǎn)的不同坐標(biāo)形式間的坐標(biāo)轉(zhuǎn)換。測量坐標(biāo)基準(zhǔn)轉(zhuǎn)換是指不同基準(zhǔn)下測量坐標(biāo)之間的轉(zhuǎn)換，其關(guān)鍵在于確定坐標(biāo)轉(zhuǎn)換參數(shù)。由于本人的知識水平和計(jì)算機(jī)編程語言掌握的有限，論文還存在許多不足之處，本文存在不足建議在以下方面繼續(xù)做工作：（1）本文主要只介紹幾種轉(zhuǎn)換模型，還有其他的幾種模型有待去實(shí)現(xiàn)。（2）對這些模型和公式的具體適用在什么情況以及他們的優(yōu)缺點(diǎn)要有詳細(xì)闡明。（3）對程序界面設(shè)計(jì)功能進(jìn)行完善，以更好的為用戶使用。（4)在編程方面應(yīng)多下功夫以實(shí)現(xiàn)各種坐標(biāo)系統(tǒng)之間的轉(zhuǎn)換。參考文獻(xiàn)[1]孔祥元，郭際明，劉宗泉.大地測量學(xué)基礎(chǔ)[M].武漢：武漢大學(xué)出版社，2009.[2]李征航，黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢：武漢大學(xué)出版社，2005.[3]劉大杰，白征東等.大地坐標(biāo)轉(zhuǎn)換與GPS控制網(wǎng)平差計(jì)算及軟件系統(tǒng)[M].同濟(jì)大學(xué)出版社，1997：16-23.[4]潘正風(fēng)，程效軍，成樞，王騰軍，宋偉東，鄒進(jìn)貴.數(shù)字測圖原理與方法[M].武漢：武漢大學(xué)出版社，2009.[5]沈云中.獨(dú)立坐標(biāo)系中GPS網(wǎng)的坐標(biāo)變換方法.工程勘察，1998(1).[6]朱華統(tǒng).大地坐標(biāo)系的建立[M].北京:測繪出版社,1986.[7]徐仕琪，張曉帆，周可法，趙同陽.關(guān)于利用七參數(shù)法進(jìn)行WGS—84和BJ—54坐標(biāo)轉(zhuǎn)換問題的探討[J].測繪與空間地理信息，2007,30[8]曾文憲,陶本藻.3維坐標(biāo)轉(zhuǎn)換的非線性模型[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2003,28(5):566-568.[9]韓曉冬，陶華學(xué)，董軍.空間直角坐標(biāo)系非線性坐標(biāo)轉(zhuǎn)換模型[J].工程勘察，2002,(5):27-30.[10]朱華統(tǒng)，楊元喜，呂志平.GPS坐標(biāo)系統(tǒng)的變換[M].北京：測繪出版社，1994.[11]朱華統(tǒng)，大地坐標(biāo)系的建立[M]，測繪出版社，2008.5，19-20[12]孫建，北京54坐標(biāo)系中的坐標(biāo)計(jì)算[J]，水運(yùn)工程，2010.2,25[13]寧津生，現(xiàn)代大地測量參考系統(tǒng)[J]，測繪學(xué)報(bào)，2012.5,5-7[14]周崇山、李紅波、劉偉，測量中坐標(biāo)轉(zhuǎn)換計(jì)算程序的開發(fā)及應(yīng)用[J]，江西科學(xué)，2011，29.1-3。[15]洪大永.GPS全球定位系統(tǒng)技術(shù)及其應(yīng)用[M].廈門大學(xué)出版社.2010，32-35.附錄大地坐標(biāo)與空間直角坐標(biāo)相互轉(zhuǎn)換程序PublicXAa,XAb,XAc,XAe2,XAe12AsDoublePublicBJa,BJb,BJc,BJe2,BJe12AsDoublePublicWGSa,WGSb,WGSc,WGSe2,WGSe12AsDoublePublicPIAsDoublePrivateSubCommand1_Click()Dima,e2,N,Ba,Br,La,Lr,H,X,Y,Z,tgB,tgB1AsDoublea=BJae2=BJe2H=Val(Text3.Text)Br=Rad(Text1.Text)Lr=Rad(Text2.Text)N=a/Sqr(1-e2*(Sin(Br))^2)X=(N+H)*Cos(Br)*Cos(Lr)Y=(N+H)*Cos(Br)*Sin(Lr)Z=(N*(1-e2)+H)*Sin(Br)Text4.Text=Round(X,3)Text5.Text=Round(Y,3)Text6.Text=Round(Z,3)EndSubPrivateSubCommand2_Click()DimB(),tgB()AsDoubleDimBa,Br,La,Lr,c,p,k,t0,t1,t2AsDoublea=BJae2=BJe2e12=BJe12X=Val(Text4.Text)Y=Val(Text5.Text)Z=Val(Text6.Text)c=a*Sqr(1+e12)p=(c*e2)/Sqr(X^2+Y^2)k=1+e12t0=Z/Sqr(X^2+Y^2)t1=t0+(p*t0)/Sqr(k+t0^2)DoUntilAbs(t1-t2)=0t1=t2t2=t0+(p*t1)/Sqr(k+t1^2)LoopBr=Atn(t2)Lr=PI+Atn(Y/X)w=Sqr(1-e2*(Sin(Br)^2))N=a/wH=Sqr(X^2+Y^2)/Cos(Br)-NBa=Txt(Br)La=Txt(Lr)Text1.Text=BaText2.Text=LaText3.Text=HEndSubPrivateSubCommand3_Click()Dima,e2,N,Ba,Br,La,Lr,H,X,Y,Z,tgB,tgB1AsDoublea=XAae2=XAe2H=Val(Text3.Text)Br=Rad(Text1.Text)Lr=Rad(Text2.Text)N=a/Sqr(1-e2*(Sin(Br))^2)X=(N+H)*Cos(Br)*Cos(Lr)Y=(N+H)*Cos(Br)*Sin(Lr)Z=(N*(1-e2)+H)*Sin(Br)Text4.Text=Round(X,3)Text5.Text=Round(Y,3)Text6.Text=Round(Z,3)EndSubPrivateSubCommand4_Click()DimB(),tgB()AsDoubleDimBa,Br,La,Lr,c,p,k,t0,t1,t2AsDoublea=XAae2=XAe2e12=XAe12X=Val(Text4.Text)Y=Val(Text5.Text)Z=Val(Text6.Text)c=a*Sqr(1+e12)p=(c*e2)/Sqr(X^2+Y^2)k=1+e12t0=Z/Sqr(X^2+Y^2)t1=t0+(p*t0)/Sqr(k+t0^2)DoUntilAbs(t1-t2)=0t1=t2t2=t0+(p*t1)/Sqr(k+t1^2)LoopBr=Atn(t2)Lr=PI+Atn(Y/X)w=Sqr(1-e2*(Sin(Br)^2))N=a/wH=Sqr(X^2+Y^2)/Cos(Br)-NBa=Txt(Br)La=Txt(Lr)Text1.Text=BaText2.Text=LaText3.Text=HEndSubPrivateSubCommand5_Click()Dima,e2,N,Ba,Br,La,Lr,H,X,Y,Z,tgB,tgB1AsDoublea=WGSae2=WGSe2H=Val(Text3.Text)Br=Rad(Text1.Text)Lr=Rad(Text2.Text)N=a/Sqr(1-e2*(Sin(Br))^2)X=(N+H)*Cos(Br)*Cos(Lr)Y=(N+H)*Cos(Br)*Sin(Lr)Z=(N*(1-e2)+H)*Sin(Br)Text4.Text=Round(X,3)Text5.Text=Round(Y,3)Text6.Text=Round(Z,3)EndSubPrivateSubCommand6_Click()DimB(),tgB()AsDoubleDimBa,Br,La,Lr,c,p,k,t0,t1,t2AsDoublea=WGSae2=WGSe2e12=WGSe12X=Val(Text4.Text)Y=Val(Text5.Text)Z=Val(Text6.Text)c=a*Sqr(1+e12)p=(c*e2)/Sqr(X^2+Y^2)k=1+e12t0=Z/Sqr(X^2+Y^2)t1=t0+(p*t0)/Sqr(k+t0^2)DoUntilAbs(t1-t2)=0t1=t2t2=t0+(p*t1)/Sqr(k+t1^2)LoopBr=Atn(t2)Lr=PI+Atn(Y/X)w=Sqr(1-e2*(Sin(Br)^2))N=a/wH=Sqr(X^2+Y^2)/Cos(Br)-NBa=Txt(Br)La=Txt(Lr)Text1.Text=BaText2.Text=LaText3.Text=HEndSubPrivateSubCommand7_Click()Text1.Text=""Text2.Text=""Text3.Text=""Text4.Text=""Text5.Text=""Text6.Text=""EndSubPrivateSubCommand8_Click()EndEndSubPrivateSubForm_Load()PI=3.14159265358979XAa=6378140:XAb=6356755.28815753:XAc=6399596.65198801:XAe2=0.006694384999588:XAe12=0.006739501819473BJa=6378245:BJb=6356863.01877305:BJc=6399698.90178271:BJe2=0.006693421622966:BJe12=0.006738525414683WGSa=6378137:WGSb=6356752.3142:WGSc=6399593.6258:WGSe2=0.0066943799013:WGSe12=0.00673949674227EndSubFunctionRad(TxtAsString)AsDoubles=Val(Txt)du=Fix(s)fe=Fix((s-du)*100)mi=Round((s-du-fe*0.01)*10000,4)Rad=(du+fe/60+mi/3600)*PI/180EndFunctionFunctionTxt(ra)AsDoubled=ra*180/PIdu=Fix(d)fe=Fix((d-du)*60)mi=Round(((d-du)*60-fe)*60)Txt=du+fe*0.01+mi*0.0001EndFunction三維坐標(biāo)轉(zhuǎn)換程序設(shè)計(jì)代碼PrivateSubCommand1_Click()Dimk#,Ex#,Ey#,Ez#,dx#,dy#,dz#DimX#,z#,Y#,Yy#,Xx#,Zz#k=Val(Text1.Text)Ex=Val(Text2.Te