用EXCEL进行高斯投影换算
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、用ＥXCＥL进行高斯投影换算
从经纬度BＬ换算到高斯平面直角坐标XY（高斯投影正算）,或从XＹ换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差，大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。笔者发现,用ＥXCＥL可以很直观、方便地完成坐标换算工作，不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。下面以5４系为例,介绍具体的计算方法。
COS(B)
K２
＝0．0０6７385254１５*Ｊ2＊Ｊ2
L2
=Ｉ2*I２
M2
=1+K2
N2
=6３996９8.９01８/SQRＴ(M2)
O2
=Ｈ2*Ｈ2＊Ｊ2*J2
P2
=Ｉ2*J2
Q２
=P２*Ｐ２
R2
=(３20０5.78006＋Q2*(１33.９２133+Q２*0.7０31)）
S2
=６367558．496８6*E２/５7．２95779513０8-P2*J2*R2+(((（L2-5８）＊L2+6１)*
起算数据Ｌ0
B2
=INT(A2）＋(INＴ(A2*100)－INT(A2)*100)／60+(A2＊１０0００-IＮＴ(A2*１00)＊１0０)/３６00
把L0化成度
C2
以度小数形式输入纬度值,如38°１４′2０〃则输பைடு நூலகம்３８．1420
起算数据B
D2
以度小数形式输入经度值
起算数据Ｌ
E2
=ＩNT(C2)＋(IＮＴ（C２*100)-ＩNT(C２)*１０0）/６０+（C2＊1０00０－INT（C2*１00)＊100)/３６00
完成经纬度ＢＬ到平面直角坐标XY的换算,在EXＣEL中大约需要占用２1列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。在ＥＸCＥL中，输入公式的起始单元格不同，则反映出来的公式不同，以公式从第2行第1列（A２格）为起始单元格为例，各单元格的公式如下:
单元格
单元格内容
说明
A2
输入中央子午线,以度.分秒形式输入，如115度30分则输入115.3０
O2/３0+（4*K２+5）*Ｍ2-L２)*O2/1２+1)*Ｎ2*Ｉ2*O2/２
计算结果Ｘ
T２
＝((((Ｌ2－18）＊L2－（58*L2-14）＊K2+5）*O2/20+Ｍ2-L2）*Ｏ2／６+1)*Ｎ２*（H2*Ｊ2)
计算结果Y
表中公式的来源及EXCＥＬ软件的操作方法，请参阅有关资料，这里不再赘述。按上面表格中的公式输入到相应单元格后，就可方便地由经纬度求得平面直角坐标。当输入完所有的经纬度后,用鼠标下拉即可得到所有的计算结果。表中的许多单元格公式为中间过程，可以用ＥXCEL的列隐藏功能把这些没有必要显示的列隐藏起来，表面上形成标准的计算报表,使整个计算表简单明了。从理论上讲，可计算的数据量是无限的,当第一次输入公式后，相当于自己完成了一软件的编制,可另存起来供今后重复使用，一劳永逸。
把Ｂ化成度
F2
=INT(D2)＋(INT（Ｄ2＊1０0）-ＩNＴ(D2)*10０)／60+(D2＊1000０-ＩＮT(D2*１00)＊10０)/３６00
把Ｌ化成度
Ｇ2
=F2-B2
Ｌ-L0
H2
=G２／5７.29577９5130８23
化作弧度
Ｉ2
=ＴAN(ＲADIAＮＳ(Ｅ2))
Taｎ（B）
J2
＝CＯＳ(RＡＤIANＳ(E2))
从经纬度BＬ换算到高斯平面直角坐标XY（高斯投影正算）,或从XＹ换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差，大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。笔者发现,用ＥXCＥL可以很直观、方便地完成坐标换算工作，不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。下面以5４系为例,介绍具体的计算方法。
COS(B)
K２
＝0．0０6７385254１５*Ｊ2＊Ｊ2
L2
=Ｉ2*I２
M2
=1+K2
N2
=6３996９8.９01８/SQRＴ(M2)
O2
=Ｈ2*Ｈ2＊Ｊ2*J2
P2
=Ｉ2*J2
Q２
=P２*Ｐ２
R2
=(３20０5.78006＋Q2*(１33.９２133+Q２*0.7０31)）
S2
=６367558．496８6*E２/５7．２95779513０8-P2*J2*R2+(((（L2-5８）＊L2+6１)*
起算数据Ｌ0
B2
=INT(A2）＋(INＴ(A2*100)－INT(A2)*100)／60+(A2＊１０0００-IＮＴ(A2*１00)＊１0０)/３６00
把L0化成度
C2
以度小数形式输入纬度值,如38°１４′2０〃则输பைடு நூலகம்３８．1420
起算数据B
D2
以度小数形式输入经度值
起算数据Ｌ
E2
=ＩNT(C2)＋(IＮＴ（C２*100)-ＩNT(C２)*１０0）/６０+（C2＊1０00０－INT（C2*１00)＊100)/３６00
完成经纬度ＢＬ到平面直角坐标XY的换算,在EXＣEL中大约需要占用２1列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。在ＥＸCＥL中，输入公式的起始单元格不同，则反映出来的公式不同，以公式从第2行第1列（A２格）为起始单元格为例，各单元格的公式如下:
单元格
单元格内容
说明
A2
输入中央子午线,以度.分秒形式输入，如115度30分则输入115.3０
O2/３0+（4*K２+5）*Ｍ2-L２)*O2/1２+1)*Ｎ2*Ｉ2*O2/２
计算结果Ｘ
T２
＝((((Ｌ2－18）＊L2－（58*L2-14）＊K2+5）*O2/20+Ｍ2-L2）*Ｏ2／６+1)*Ｎ２*（H2*Ｊ2)
计算结果Y
表中公式的来源及EXCＥＬ软件的操作方法，请参阅有关资料，这里不再赘述。按上面表格中的公式输入到相应单元格后，就可方便地由经纬度求得平面直角坐标。当输入完所有的经纬度后,用鼠标下拉即可得到所有的计算结果。表中的许多单元格公式为中间过程，可以用ＥXCEL的列隐藏功能把这些没有必要显示的列隐藏起来，表面上形成标准的计算报表,使整个计算表简单明了。从理论上讲，可计算的数据量是无限的,当第一次输入公式后，相当于自己完成了一软件的编制,可另存起来供今后重复使用，一劳永逸。
把Ｂ化成度
F2
=INT(D2)＋(INT（Ｄ2＊1０0）-ＩNＴ(D2)*10０)／60+(D2＊1000０-ＩＮT(D2*１00)＊10０)/３６00
把Ｌ化成度
Ｇ2
=F2-B2
Ｌ-L0
H2
=G２／5７.29577９5130８23
化作弧度
Ｉ2
=ＴAN(ＲADIAＮＳ(Ｅ2))
Taｎ（B）
J2
＝CＯＳ(RＡＤIANＳ(E2))