理论基础:标准地图分幅基本比例
我国现行地图中有8种基本比例尺:1:100万,1:50万,1:25万,1:10万,1:5万,1:2.5万,1:1万,1:5000。根据1913年巴黎地图会议指定的国际统一的标准,对于整个地球,经度一共有360°,纬度从赤道开始起算,向南向北分别有90°。根据经差6°,纬差4°,可以将地球沿着经线一刀一刀竖切,像切西瓜一样切成60块。纵向横着从赤道开始可以切成19块。多出的南北两极,单独作为一块。需要注意的是,经度切分,是从西经180度开始,并不是从0°开始切。另外,由于由于子午线向南北极收敛,因此,整个图幅呈梯形。其划分的方法和编号,随比例尺不同而不同。1:100万地图的标准分幅是经差6°和纬差4°,由于随着纬度增高地图面积迅速缩小,所以规定在纬度60°到76°之间双幅合并,即每幅图包括经差12°,纬差4°。在纬度76°到88°之间由四幅合并,即每幅图包括经差24°,纬差4°;纬度88°以上单独为一幅。
我国1:100万地形图的分幅按这个国际标准。如北京在1:100万图上处于第J列第50行,故编号为J50。其它基本比例尺地图分幅编号是基于1:100万地图,按照规定的经度、纬度差采用加密的方法进一步分幅。
实践 : 确定1:10000 蚌埠市的市辖区地形图的数量及编号
需要获得蚌埠市市辖区的shp文件(这里使用的是2015年全国县级矢量数据)
通过在ArcGIS软件中查看蚌埠市市辖区的数据属性,可以了解蚌埠市市辖区的经纬度范围。
根据这个经纬度范围,我们可以确定1:10000的比例尺网格左下角的经纬度。可以确定,蚌埠所处的1比100万地形图编号为I50。可以先把I50这张100万地形图范围先画出来。I50的范围经度为114-120°,纬度为32-36°。使用创建渔网工具,模板范围选择经度:114-120,纬度32-36。注意此时数据框属性中单位的是以度的形式表示。行列数输入96(考虑一下为什么是96)。几何类型选择POLYGON。
这时可以再用渔网工具重新绘制覆盖蚌埠市辖区的1:10000的网格,见下图,从上图中可以确定蓝色网格左下角的经纬度为117°11.25′,32°47.5′,右上角的经纬度为117°33.75′,33°2.5′,用此数据对渔网工具进行设置,就可以绘出覆盖蚌埠市辖区的最小范围的1比1万图幅,当然这是一个矩形,我们只需要市辖区覆盖到的网格就行了。
接下来计算上图左下角这个网格在I50中的行列号(当然也可以把整幅I50都标上1:10000图幅的编号)。需要用到I50的经纬度范围,1:10000地图的经纬度差以及上图左下角经纬度数据,几组数据。利用下面所给的公式计算,结果显示左下角网格在I50中的第52列,第77行。因此,按照地形图编号的规则,该网格的编号为 I50G077052。如何将编号标注在图中,如果你掌握了arcgis软件的基本操作,应该能想到处理方法,如果你还没有这个能力,我在这里给一点提示。
其中根据某地的经纬度或图幅西南图廓点的经纬度计算出它所在的地图编号公式为:
行号 = 460/2.5-int((!纬度!%4)60/2.5)
列号 = int((!经度!%6)*60/3.75)+1
提示:绘制的网格实际上是若干个面对象,可以在其属性列表中看到,可以对这些面对象增加一个属性,这个属性值就可以设置为这个面的地图编号,那么只要将这个属性标注出来就行了。当然这只是一个提示,你可以照着这个思路去做,这个过程中可能还会碰到一些问题,你可以尝试去解决,如果还是有问题可以在评论区留言!