使用软件：Arcgis10.7 坐标系：CGCS2000 需了解的基础工具与知识：渔网工具、空间连接工具、利用python进行字段计算、model模型工具箱

1、模板数据 模板数据来源于全国地理信息资源目录服务系统 1:100w全国基础地理数据库分幅数据提取得到的边界，该网站公布了1:100w地理数据库，按标准图幅分为77幅，只能单独下载。拼接该数据的BOUA图层后新建一字段“country”，令其全为China，利用该字段融合所有离散面，结果如下图：

2、制作流程 首先利用渔网工具制作中国区域1:100w基本比例尺图幅，中国区域范围：维度0° ~ 56°、经度72° ~ 138，每一幅1:100w分幅的维差4°，经差6°。如下图，像元宽度为每幅图的经差，故为6，计算得到列数为11；像元高度为每幅图的维差，故为4，计算得到行数为14；点击取消标注点，选择输出几何类型为面，输出要素类命名为A，即对应1:100w比例尺编码。故填写参数如下： 如下图，打开生成的属性表，新建双精度字段:lat_center、lon_center，利用“计算几何”得到每个分幅中心点经纬度。 但发现出现偏差，具体原因鄙人不太清楚，故只能选择以上述模板数据作为模板范围生成1:100w图幅，查看其属性表后发现其结果正确无误。

接下来通过“空间连接”工具将“country”字段连接"A图层"（为了后续筛选中国区域分幅，当然也可以直接在下图取消“保留所有目标要素”），具体参数如下图：

完成“country”字段连接后，删除"A"图层，更改"A_SpatialJoin"图层为"A"图层。 下图为经过"空间连接"后，该图层属性表，拥有“country”字段，质心坐标正确。

进行下一步：添加一文本类型、长度为10的字段，命名为“图幅编号A”，然后利用字段计算器根据分幅图的质心经纬度计算得到各分幅的图幅编号，计算公式、结果与python脚本函数如下：

1:100万图幅的行号和列号计算公式 1:100行号 = [纬度/4] +1    []表示取整,1为A，2为B… 1:100列号 = [经度/6]+31    []表示向下取整

以“图幅编号A”字段标注图层，结果如下：

1、分幅过程 创建渔网，1:50w分幅维差为2、经差为3，行列数相比1:100w翻倍，输出要素类命名为B，即对应1:50w比例尺编码，具体参数如下:

2、编码过程 进行“空间连接”操作，谈一下这个操作的意义：比如对于J50分幅，将其图幅编号赋予对应的4个1:50w分幅，同时将质心经纬度、“country”字段保留，后续不用新建直接使用，填写参数如下图（注意匹配选项）：

如上图，打开修改图层名称后的“B”图层的属性表，删除和关闭不需要的字段，更新质心经纬度，添加一文本类型、长度为10的字段，命名为“图幅编号B”，然后利用字段计算器根据分幅图的质心经纬度计算得到各分幅的图幅编号，计算公式、结果与python脚本函数如下：

非1:100万图幅的行号和列号计算公式 行号 = 4º/纬差 - [（纬度/4º）/纬差]    []为取整，（）为取余 列号 = [（经度/4º）/经差]+1       []为取整，（）为取余 1:50万图幅，经差为：3°，纬差为：2°

以“图幅编号B”字段标注图层，叠加1:100w分幅，结果如下：

其余基本比例尺分幅编码步骤与上述1:50w完全一致

按照上表经纬差进行其余比例尺分幅出图，下方展示其中一些过程数据

如下图：1:25w基本比例尺分幅图"创建渔网"参数设置

如下图：J50D012012分幅内部1:5w、1:2.5w分幅及编号

如下图：1:2.5w分幅出图结果，J50下共计48*48=2304个1:2.5w分幅。

如下图：J50F048048分幅内部1:1w、1:5000分幅及编号

以下为一个小流域所占据的1:1w比例尺地形图的图幅与编号，写这篇文章也是因为我想快速知道一些区域它位于哪些比例尺地形图上，编码是多少，具体有哪些，这是我的目的所在。 下图则是点击提取了该小流域左上方区域一点度分秒格式经纬度，并在我个人制作的安卓工具集中使用图幅编号工具计算得到该位置1:1w地形图编号，结果与之相同，没有差错。

通过上述流程，我们获得了CGCS2000大地坐标系下，中国区四至坐标范围内，1:5000到1:100w八个基本比例尺分幅图及其编号，统一以其比例尺等级命名，如1:100w命名为A、1:10000命名为G。 不过上述数据是在四至范围内部，是一个矩形数据，为了精简数据，我们可以筛选中国区范围内的分幅数据，之前空间连接的“country”字段便是用来筛选的，在这里我使用Arcgis模型工具进行遍历筛选及导出，过程如下：

上述模型流程是遍历地理数据库要素类，对8个要素类进行筛选、再导出，过程及结果如下： 1650724963702.png

到这里并没有结束，当我打开八个图层属性表的时候，我发现除了A图层，其他图层的图幅编号字段都没有了… 好一番找发现是这个原因，如下： 选择工作空间后就直接确定了“字段映射”，而且还不随迭代要素的变化而变化，所以最后8个要素类都只有这些字段。

最后实在没找到解决办法，python脚本写起来没必要，所以直接进行批处理，浏览填充一波走起，速度也不慢，如下：

#（二）小结 以上，便是本篇博文主要流程。 同时声明，此教程仅为个人依据分幅编号原理所制，可能与标准数据有所差异，仅供参考。