利用图号查询国家基本比例尺地形图图角坐标程序的开发

刘永涛

(河北省测绘产品质量监督检验站，河北石家庄050031)

一、概 述

在国家和省级的基本比例尺地形图测绘工作中，经常会遇到查询指定图号的图幅图角坐标的问题，而通过图幅结合表查图号，再通过图号找图纸查看图角坐标的方法，在大批量的查询工作中效率不高，而我国基本比例尺地形图之间的数学关系简单明了，地形图编号与图角坐标的关系可利用简单的数学关系式进行表达。本文利用国家基本比例尺地形图图号与图角坐标的数学关系，基于Visual Basic 6．0编程语言，开发了1∶1万地形图根据图号查询图角坐标的程序，除1∶100万地形图外其他采用高斯投影的国家基本比例尺地形图图角坐标查询，只需对本程序中的图幅经、纬差及不同分带的中央子午线进行修改即可实现。

二、国家基本比例尺地形图分幅与编号

我国基本比例尺地形图分为：1∶100万、1∶50万、1∶25 万、1∶10 万、1∶5 万、1∶2．5 万、1∶1 万和1∶5000地形图。

我国基本比例尺地形图分幅与编号均以1∶100万地形图为基础，进行不同比例尺地形图的分幅与编号，1∶100万地形图分幅经度差为6°，纬度差为4°。其他六种比例尺地形图分别在1∶100万地形图的基础上进行划分，基本划分规则如下：

1)1∶50万地形图：2行2列，共4幅，经差3°，纬差 2°。

2)1∶25万地形图：4行4列，共16幅，经差1°30'，纬差 1°。

3)1∶10万地形图：12行12列，共144幅，经差30'，纬差 20'。

4)1∶5万地形图：24行24列，共576幅，经差15'，纬差 10'。

5)1∶2．5 万地形图：48 行48 列，共2304 幅，经差 7'30″，纬差 5'。

6)1∶1万地形图：96行96列，共9216幅，经差3'45″，纬差 2'30″。

7)1∶5000地形图：192行192列，共36 864幅，经差 1'52．5″，纬差 1'15″。

1∶50万 ～1∶5000地形图的编号均以1∶100万地形图编号为基础，采用1∶100万地形图编号加比例尺代码和1∶100万图幅中各比例尺地形图的列号与行号的组合。如J50B002002，即表示1∶100万地形图J50中右下角的第四幅1∶50万地形图的编号，其中B为1∶50万比例尺代码，1∶25万～1∶5000的比例尺代码分别为 C、D、E、F、G、H。

三、国家基本比例尺地形图图号与图角坐标的数学关系

由以上我国基本比例尺地形图的分幅和编号方法可以看出，我国的基本比例尺地形图图号与图角大地坐标的关系易于通过数学公式表达。仍以1∶100万地形图J50为例，以J50西北角坐标(40°，114°)为基础列出其图幅内1∶1万比例尺地形图图号与图角坐标的数学关系。

设C为1∶1万地形图在1∶100万地形图内分幅的列数(列号)，R为行数(行号)，ΔB为该比例尺地形图纬差，ΔL为经差。

西南角：

西北角

东北角

东南角

四、程序开发

本文利用 Visual Basic 6．0编写了 J50(1∶100万)图幅范围内的1∶1万地形图图角坐标查询程序，其他比例尺的图角坐标(计算时应注意不同比例尺地形图分带的不同)可根据图号中的比例尺代码及最大行号、列号对程序添加判断后进行编写。

受篇幅所限，本文只列举主程序代码，辅助计算代码不再赘述。

1．程序界面

查询程序界面如图1所示。

图1 部分动画制作界面

2．程序代码

由于程序在计算图角平面坐标时需要进行高斯正算，故应首先定义高斯正算函数，代码如下：

定义函数完成后，按照图号与图角坐标的关系，编写图角坐标计算程序，代码如下：

五、结束语

利用本文提供的数学关系及程序设计思路，很容易编写出全国范围内1∶50万～1∶5000地形图的图角坐标查询程序，为指定图号的地形图图角坐标的快速查询提供了快捷、准确的计算工具。

[1]GB/T 13989—92国家基本比例尺地形图分幅与编号[S]．北京：中国标准出版社，1993．