利用EXCEL实现新旧图幅号转换

第 27卷 � 第 4期 2011年 2月 甘肃科技 G ansu Sc ience and T echnology Vol. 27� No. 4 F eb. � 2011 利用 EXCEL实现新旧图幅号转换 于树晖 (甘肃省核工业地质局,甘肃 兰州 730000) 摘 � 要:在使用和制图时常会碰到国家基本比例尺地形图新旧图幅编号之间的变换问题。本文通过对新旧图幅编 号方法和规律的研究, 提出了利用 EXCEL完成对新旧图幅编号之间的变换。 关键词:图幅号; 转换; EXCEL; 函数 中图分类号: P264 � � 1� 1万、1� 2. 5万、1� 5万、1� 10万、1� 20 万、1� 50万和 1� 100万地形图是我国的国家基本 比例尺地形图, 其图幅编号现有两种形式: 一种是 1991年以前地形图分幅编号标准产生的,称为旧图 幅编号;另一种是 1991年后新的国家地形图分幅编 号标准所产生的,称为新图幅编号。 1 � 新旧图幅分幅编号公式 1. 1� 旧图幅分幅编号公式 1)按下式计算 1� 1 000 000地形图图幅编号: a = [ �/4 ] + 1 ( 1) b= [�/6 ] + 31 ( 2) 式中: [ ] ! ! ! 表示商取整; a! ! ! 1� 1 000 000地形图图幅所在纬度带字 符码所对应的数字码; b! ! ! 1� 1 000 000地形图图幅所在经度带的 数字码; �! ! ! 图幅内某点的经度或图幅西南图廓点 的经度; �! ! ! 图幅内某点的纬度或图幅西南图廓点 的纬度。 上式计算出的 a需要转换成 1� 1 000 000地形 图图幅所在纬度带字符码如图 1所示, 按下式计算: A= CHAR( a + 64) 式中: A! ! ! 1� 1 000 000地形图图幅所在纬度带字 符码。 2) 1� 500 000、1� 200 000和 1� 100 000地形 图的分幅以 1� 1 000 000地形图的基础上分幅; 1� 50 000、1� 10 000地形图的分幅以 1� 100 000地形图的基础上分幅; 1� 25000地形图以 1� 50 000地形图的基础上 分幅。具体分幅见表 1。 图 1� 北半球 1� 1 000 000地形图分幅编号 � � 3)按下式迭代计算所求比例尺地形图图号: W = V - [ ( �/∀ �#) /∀ �] ∃b+ [ (�/∀ �#) /∀ �] 式中: ( ) ! ! ! 表示商取余; [ ] ! ! ! 表示商取整; �! ! ! 图幅内某点的经度或图幅西南图廓点 的经度; �! ! ! 图幅内某点的纬度或图幅西南图廓点 的纬度; n! ! ! 所求比例尺地形图可在基础图幅分得的 列数; ∀ �#! ! ! 所求比例尺地形图前基础图幅的经 差; ∀ �#! ! ! 所求比例尺地形图前基础图幅的纬 差; ∀ �! ! ! 所求比例尺地形图分幅的经差; ∀ �! ! ! 所求比例尺地形图分幅的纬差。 表 1� 具体分幅 1. 2� 新图幅分幅计算公式 1) 1� 1 000 000沿用国际分幅标准, 计算公式 和旧图幅分幅一样。 2)各比例尺地形图采用不同的字符作为其比 例尺的代码,见表 2。 1� 500 000~ 1� 5000地形图的图号均有其所 在 1� 1 000 000的图号、比例尺代码和各图幅的行 列号共 10位码组成,如图 2所示。 表 2� 比例尺的代码 比例尺 1� 500 000 1� 250 000 1� 100 000 1� 50 000 1� 25 000 1� 10 000 1� 5 000 代码 B C D E F G H 图 2� 1� 500 000~ 1� 5000地形图的图号示意 � � 3)按下式计算所求比例尺地形图在 1� 1 000 000地形图图号后的行、列号: c= 4 /∀ �! [ (�/4 ) /∀ �] ( 3) d = [ ( �/6 ) /∀ �] + 1 ( 4) 式中: ( ) ! ! ! 表示商取余; [ ] ! ! ! 表示商取整; c! ! ! 所求比例尺地形图在 1� 1 000 000地形 图图号后的行号; � � d! ! ! 所求比例尺地形图在 1� 1 000 000地形 图图号后的列号; �! ! ! 图幅内某点的经度或图幅西南图廓点 的经度; �! ! ! 图幅内某点的纬度或图幅西南图廓点 的纬度; ∀ �! ! ! 所求比例尺地形图分幅的经差; ∀ �! ! ! 所求比例尺地形图分幅的纬差。 2� 利用 EXCEL实现新旧比例尺图幅编码转 换 2. 1� 已知图幅内某点的经纬度,计算其新旧图幅编 号 � � 以北纬 36 52#58%,东经 118 17#40%在 1� 50 000 地形图上为例,如图 3所示。 41第 4期 � � � � � � � � � � � � 于树晖: 利用 EXCEL实现新旧图幅号转换 图 3� 某点在 1� 50 000新旧图幅编号 � � B5 = INT( B2 /6) ∃ 6, B6= INT( B3 /4) ∃ 4是确 定该点所在 1� 1000 000基础图幅西南图廓点经、 纬度,从而确定该点所在 1� 1 000 000基础图幅的 范围。 2. 2� 已知新图幅号计算该图幅点西南图廓点的经、 纬度,并计算其对应的旧图幅编码 � � 将式 ( 1) ~ ( 4)进行整理便可得到以新图幅号 求西南图廓点的经、纬度的公式: B5= ( C4- 31) ∃ 6+ ( F4- 1) ∃ 15 /60 E5= ( B4- 1) ∃ 4+ ( 4 ∃ 60 /10- E4) ∃ 10 /60 然后根据求得的经纬度计算其所在的旧图幅 号。 (注:上式 & 15∋、& 10∋分别是 1� 50 000地形图 的经差和纬差,实际应用中可以根据图幅号不同判 断其比例尺,确定经差、纬差。如图 4所示。 图 4� 由新图幅号变换成旧图幅号 2. 3� 已知旧图幅号计算该图幅点西南图廓点的经、 纬度,并计算其对应的新图幅编码 � � 这里要特 别说明一下 , B 6 = MOD ( D2, 12) , B7 = INT ( D2 /12 )是判断该图幅在 1� 1 000 000图幅上的行列号 ; C 6 = MOD ( CODE ( E 2 ) - 64, 2 ) , C7 = INT ( ( CODE ( E 2 ) - 64 ) / 2)是判断该图幅号在 1� 10000图幅上的行列 号。 B6、C6等于 0时, 表示该图幅基础图幅的 最右面一列; B7、C7小于等于 1时, 表示该图幅 在基础图幅的最上面一列 , 可以利用 if( )函数 分情况实现图幅计算。 3� 结论 利用了 EXCEL函数解决了新旧图幅转化的问 题,既避免了查询接图表的繁琐,又避免了使用编程 软件的专业性,在实际工作中根据不同比例尺随时调 整参数,可以达到任意比例尺新旧图幅之间的转化。 参考文献: [ 1] � GB /T 13989- 92,国家基本比例尺地图国家基本比例 尺地形图分幅和编号 [ S]. [ 2] � 高井祥. 测量学 [ M ]. 徐州: 中国矿业大学出版社, 2006. (上接第 60页 ) 参考文献: [ 1] � 张培红,陈宝智, 刘丽珍. 大型公共建筑物火灾时人员 疏散行为规律研究 [ J]. 中国安全科学学报, 2001, 11 ( 2). [ 2] � 霍然,衷宏永. 性能化建筑防火 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 与 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 [ M ]. 安徽 科学技术出版社, 1997. [ 3] � 阎卫东. 建筑物火灾时人员行为规律及疏散研究 [M ] .东北大学, 2006. [ 4] � GB50084- 2001, 自动喷水灭火系统 设计规范 民用建筑抗震设计规范配电网设计规范10kv变电所设计规范220kv变电站通用竖流式沉淀池设计 [ S]. [ 5] � DGJ08- 88- 2000,民用建筑防排烟技术古城 [ S] . 42 � � � � � � � � � � � 甘 � 肃 � 科 � 技 � � � � � � � � � � � � � � � � 第 27卷