国家基本比例尺系列地图是新时期自然资源系统测绘地理信息重要基础数据，广泛应用于地震救灾、滑坡泥石流监测、应急保障等领域^[1-3]。在信息化时代，快速准确获取野外任意点位坐标，并获得其在国家基本比例尺地图中的图幅号，仍具有较大现实意义。随着时代发展，求取图幅号的方式方法出现了不同的形式^[4-7]。目前在外业测量中，往往仍需人工计算、查阅纸质图幅结合表或者回到内业在PC机上查询获取图幅号，效率低，时效性差、容易出错，给野外测量工作带来一定程度的不便。

随着移动通讯进入5G时代，高网速、大容量和低延时的移动通讯网络已成为国家信息高速公路基础设置。智能手机早已成为野外测量人员日常携带的必需设备。Android系统是由Google公司于2007年发布的智能手机平台，是一个开放、自由的终端平台^[8-9]。为在野外测量现场实时获取各种基本比例尺下的图幅号信息提供了可能，本文依据Android智能手机操作系统的开放性，内置GPS信号接收模块和手机基站定位等硬件基础，研制了一款可方便快捷的动态实时获取多比例尺下图幅号的手机软件，可在一定程度上提高野外作业效率。

为了便于测绘、拼接、管理和使用地形图，需要将广大地区的地形图划分为若干单幅地形图。地形图分幅与编号的方法有两大类，一类是按经纬线分幅，即按照经线和纬线来划分，左右以经线为界，上下以纬线为界，图形近似梯形；另外一种是按坐标格网分幅的矩形分幅法，即按照平面直角坐标的纵、横坐标线来划分矩形或者正方形^[10]。前者用于8种国家基本比例尺地形图(1：100万-1：5000)，后者用于工程建设大比例尺地形图。《国家基本比例尺地形图分幅与编号》(GB/T 13989-2012)是目前我国进行地形图分幅与编号的标准依据^[11-12]。

我国1：100万地形图的分幅与编号采用国际1：100万地图分幅与编号标准。按照经差6°、纬差4°，从赤道起算，每纬度4°为一行，至南、北纬88°各分为22行，依次用大写拉丁字母A、B、C、……V表示其相应行号；从180°经线起算，自西向东每经差6°为一列，全球分为60列，依次用数字表示其相应列号。国际1：100万地图编号第一位表示南、北半球，由于我国全部位于北半球，N字省略。各种基本比例尺地形图图幅大小数量关系如表1所示。

表1 不同比例尺地形图经差、纬差和行列数量关系
Table 1 The relationship between longitude difference，latitude difference and the number of rows and columns in different scale maps

在1：50万-1：5000地形图的编号中，第三位是比例尺编码，见表2所示。新图幅号编码中的后6位，分别为行号和列号编码。其规则为从上到下用三位阿拉伯数字来编号，列号编码规则为按照从左往右的顺序用3位阿拉伯数字来编号，前面不足3位的添0来补充。

表2 不同比例尺地形图对应的代码
Table 2 Codes corresponding to different scale topographic maps

核心思路是：根据点的坐标，先依据下列地图学公式，计算点所在的1：100万图幅的编号。

式(1)、(2)中，a为1：100万图幅行号所对应的数字码，b为1：100万图幅列号，^[]表示取整数。

然后再利用下式计算该地在目标比例尺地形图下的图幅行列号

(3)

(4)

式(3)、(4)中，^[]表示取整数，()表示取余数，c为行号，d为列号，φ为该点纬度值，λ为某地经度值，Δφ为某比例尺地图下的纬度差，Δλ为某比例尺地图下的经度差，各基本比例尺地形图经差纬差情况见表1，不同比例尺地形图对应的代码见表2。

本文介绍了地形图分幅与编号的基本原理，给出了以Android操作系统为平台的智能手机动态实时获取图幅号的方法。由于在国家新图幅编号系统中1：500至1：2000大比例地图图幅号采用整公里(或百米)平面坐标进行分幅或采用任意分幅，不直接由经纬度计算，本文暂未讨论由经纬度计算1：2000以上大比例尺地图图幅号的方法。目前用该方法实现的软件已广泛应用于单位野外像控测量、控制点普查等基础测绘项目实践，很好地解决了外业测量中难以方便快捷获取图幅号信息，计算费事费力、效率低、易出错的弊端。