ArcGIS制作并输出1:5 000比例尺地形图分幅及编号(一)


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

参考链接二:https://www.cnblogs.com/yzhyingcool/p/10328390.html

参考链接三:https://www.cnblogs.com/yzhyingcool/p/10132746.html


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM