1. 前言
由于项目要求,需要使用ArcGIS对某一测区范围进行1::5000比例尺分幅及编号;查阅了相关博客及《国家基本比例尺地形图分幅和编号》GB/T 13989-2012,现对制作流程以及制作过程中遇到的问题进行记录。
2. 数据介绍
因为不方便详细的介绍测区的经纬度范围,所以笼统的概括一下,本次所处理的测区涉及到1:1000 000比例尺地形图的经度范围为78°E-90°E,纬度范围为36°N-44°N,即涉及到的图幅编号分别为J44、J45、K44、K45。
1:5000比例尺地形图的分幅和编号是以1:1000 000比例尺地形图为基础,按规定的经差和纬差划分图幅,采用行列编号方法进行编号。
1:5000比例尺代码:H
经差:1′52.5″,即1.875′,亦即0.03125°
纬差:1′15″,即1.25′,亦即0.02083333……°
单个1:1000 000比例尺地形图图幅可以划分为192 * 192 的1:5000比例尺图幅。
3. 地理坐标转换
3.1 数据导出
原始测区范围在奥维地图中,因此,首先将测区从奥维中导出,格式选择 shp,坐标类型选择 WGS84(地理坐标),参数如下图所示:
3.2 将数据添加到ArcMap中
简单,不解释......
3.3 WGS84坐标到CGCS2000坐标的转换
在ArcMap的ArcToolbox工具箱中依次找到 数据管理工具-->投影和变换-->要素-->投影(ArcToolbox工具在地理处理菜单栏中),双击 投影工具,在弹出的窗口中进行如下设置:
图中 输出坐标系 选择 地理坐标系-->World-->ITRF2000:
然后,在目录中选择上一步输出的ITRF2000坐标系的shp文件,单击右键-属性,将图层坐标重新定义成GCGS2000地理坐标系:
在弹出的属性窗口中选择 地理坐标系 下的CGCS2000坐标,如下图:
此时,完成WGS84地理坐标向CGCS2000地理坐标的转换。
4. 分幅与编号
4.1 分幅
在ArcMap的ArcToolbox工具箱中依次找到 数据管理工具-->要素类-->创建渔网,双击 创建渔网 工具,在弹出的对话框中进行如下图参数设置(图为J44图幅的示例,据需设置):
上图中还需在 ‘环境...’ 中设置输出坐标系:
根据上图设置,创建渔网后,会生成一个面层(如下图),接下来根据国标对网格进行编号:
4.2 编号
4.2.1 编号规定及形式
《国家基本比例尺地形图分幅和编号》GB/T 13989-2012 中规定,1:5 000地形图的行、列编号是将1:1 000 000 地形图按所含比例尺地形图的经差和纬差划分成若干行和列,横行从上到下、纵列从左到右,按顺序分别用三位阿拉伯数字(数字码)表示,不足三位者前面补零,取行号在前、列号在后的排列形式标记。因此,J44图幅的1:5 000的编号就如下图所示:
4.2.2 属性表编辑与导出
① 首先,在ArcMap的 内容列表 的图层中找到生成的分幅,然后打开其属性表,添加X、Y、Id(长整型)三个字段:
② 然后,分别计算X、Y的质心坐标,并通过字段计算器为Id赋值:
③ 导出属性表,在属性表中选择 表选项-->导出,导出表格:
4.2.3 Excel编行、列号
① 根据GB标准以及4.2.1的示例,可以知道X方向编号应为升序列,Y方向编号应为降序列;
② 首先,将X所在列按照升序排序(注意排序时要选择 扩展选定区域);然后,新添加一列进行编号,函数式如下图,第一行的单元格要空着; 最后,复制计算出来的编号列(图中D列)重新在该列粘贴为数值格式,并在第一行单元格加上表头(见下图):
③ Y所在列按降序排序后,其操作同X列类似,不做过多解释。
4.2.4 连接Excel表
① 选择生成的分幅,右键选择 连接和关联-->连接...,在打开的 连接数据 对话框中进行如下设置,点击确定即可将Excel表挂接到图层的属性表上:
② 打开连接Excel表后的属性表,添加一个Code字段(文本类型),然后使用字段计算器,计算行列号(脚本和结果如下图所示):
5. 地理坐标转换为投影坐标
在ArcMap的ArcToolbox工具箱中依次找到 数据管理工具-->投影和变换-->要素-->投影,双击 投影工具,在弹出的窗口中进行如下图设置(此处输出坐标系选择的是 3度带+中央经度的形式,也可以选择其他形式,详见 参考链接一):
需要注意的是上图的 ‘环境...’ 中需要将 输出坐标系 设置成和 ‘投影’ 对话框中的相同:
6. 将分幅编号转为注记层
6.1 标注要素
首先,需要选择投影后图层,然后点击右键,打开“图层属性”对话框,切换到“(标注)”页,进行如下图设置;之后,同样选择该图层并右键选择 ‘标注要素’,即可完成要素标记。
6.2 标注转注记
标注完成后,继续选择被标注的图层,然后点击右键,选择 ‘将标注转换为注记’,在弹出的对话框中进行如下设置(注意,注记需要输出在数据库中,因此必须在输出位置建立一个数据库):
标注转换为注记后,会作为一个注记图层加载到 内容列表 里,其最终结果如下图所示(转换为投影坐标后已经呈现出梯形形状):
7. 输出为CAD格式
选择图层并右键,依次选择 ‘数据’-->‘导出至CAD...’,在弹出的对话框中进行如下设置即可:
上述内容是J44图幅的分幅流程,其他分幅参照该流程依次生产,最后根据测区范围,保留需要的1:5 000分幅和编号,导出即可。
<上述内容参考其他博客,并结合实际进行更改,在此感谢其他无私奉献的博主,另外如有不足请各位查阅者指出>
附言:如果测区完全包含在 1:1 000 000 图幅中的一个(例如完全在J44图幅中),不涉及跨图幅,则1:5 000比例尺的的分幅操作可以根据4.1进行。如果涉及到其他图幅,则4.1在进行创建渔网时,需要改变参数。因为上述参数中的像元高度使用的单位是度,而这个度的数字是无限循环小数,所以生成的1:5 000比例尺分幅在上下两个图幅(例如J44和K44)的接边处会存在缝隙。
该问题的具体描述及解决方法可以参考“ArcGIS制作并输出1:5 000比例尺地形图分幅及编号(二)”。
(解决方法还没有写,先在此挖个坑,回来填...... •∩•)
参考链接
参考链接一:https://blog.csdn.net/weixin_39541750/article/details/111252510