用EXCEL处理高斯坐标计算问题

《用EXCEL处理高斯坐标计算问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、第22卷第1期安徽地质Vo1．22No．12012年3月GeologyofAnhuiMarch2012文章编号：1005—6157(2012)01—70—03用EXCEL处理高斯坐标计算问题吕建国(安徽省煤田地质局第三勘探队，安徽宿州234000)摘要：本文介绍了使用EXCEL处理高斯坐标的正算、反算问题，着重阐述了计算公式的编制，具有使用方便、实用性较强的特点。关键词：EXCEL；高斯坐标；地理坐标；正算；反算中国分类号：TP317．3文献标识码：A地理坐标的单位一般是由度分秒组成的，而EX．0

2、引言CEL中角度单位是弧度制，所以首先必须将以度分秒在地质勘探、地形图测绘和公路、铁路等工程为单位的坐标转换为以弧度为单位的坐标。在EX施工中，经常涉及地理坐标与高斯坐标的换算，即CELSE作簿中分别建立正算和反算两个工作表，如下高斯投影坐标计算，包括高斯投影正算和高斯投影图1所示，在B2、C2两个单元格中分别输入纬度和反算。无论是正算还是反算，为了确保精度要求，经度坐标，特别注意，输入纬度和经度的坐标是以度常常会用到高阶导数(一般算至六阶)按照公式计为整数位，小数点后两位为十分位和分位，十分位和

3、算，或者使用《高斯、克吕格投影计算表》进行计分位是零的用0补齐，分位后为秒位。为了便于说算，不管是利用公式计算还是查表计算，计算都比明，对建立的每个单元格公式按表1进行说明：较复杂，而且精度较低，容易出错，计算速度慢。2反算随着计算机的广泛应用，利用EXCEL电子表格计算，不仅速度快而且精度高。高斯投影反算是由高斯平面坐标(x，Y)求地理坐标(L，B)。l正算建立反算工作表后，在B2、C2单元格中分别输高斯投影正算是由地理坐标(L，B)求高斯平入x、Y值，象正算一样，按照高斯投影反算公式逐面坐标(

4、X，Y)。个建立单元格公式，如表2所示：图1EXCEL作簿中分别建立正算和反算两个工作表Fig．1SeparateestablishmentofpositiveandreversecalculationworktableinEXCELworkbook收稿日期：2010—08-16作者简介：吕建国(1967一)，男，安徽六安人，工程师，现从事煤田地质测绘工作。第22卷第1期吕建国：用EXCEL处理高斯坐标计算问题71表1高斯正算单元格公式说明Table1NotesOnGausspositivecal

5、culationcellformula·单元格代码公式或输入内容单元格说明A2输入点名或代号B2输入纬度C2输入经度D2=INT(B2)+(INT(B2100、-INT(B2)100)／60+(B210000-1NT(B2IOO)100)／3600角度单位转换，即以度分秒形E2=INT(C2)+(INT(C2100)-INT(C2)100)／60+(C2IO000-INT(C2lOO)100)／3600式转换为以弧度为计算单位F2=INT((E2+6)／6)求出6。分带带号G2=(EZ-F26+3

6、13600／206264．8063求出经差H2=COS(D23．1415926536，l8o1纬度余弦值I2=6399698．902—21562．267*POWER(H2，2)+108．973*POWER(H2，4)一0．612*POWER(H2，6)+0．004POWER(H2，8)J2=32140．404-(135．3302-(O．7092-0．OO4H2H2)H2H2)H2H2K2=f0．25+0．00252H2H2)H2H2-O．04166公式计算处理L2=(0．166"H2H2-0．08

7、4)H2H2M2=f0．3333333+0．001123H2H2)H2H2-0．1666667N2=0．0083一(0．1667-(0．1968+0．004H2H2)H2H2)H2H2O2=6367558．4969D23600／206264．8063-02-(0．5+(K2+L2G202)G2G2)输出高斯坐标X值G2G2"I2)*SIN(D23．1415926536／180)*H2P2=F21000000+500000+(1+(M2+N2G2G2)G2G2)G2I2"H2输出高斯坐标Y值表2高斯

8、反算单元格公式及说明Table2NotesonGaussIleversecalculationcelformula单元格代码公式或输入内容单元格说明A2输入点名或代号B2输入x值C2输入Y值D2=INT(C21000000)E2=(B2／6367558．4969)206264．8063F2=COS((E2／3600)P10／18o)COS((E2／3600)P10／180)G2=E2+(50221746+(293622+(2350+22F2)F2)F2)POWER(IO，-1O)S