国家地形图经纬度与图幅号的彼此转换

国家基本比例尺地形图新旧图幅号转换计算方法简介在地图绘制和测量过程中，图幅号是非常重要的标识和查询方式。

由于生产的地图不断更新和更换，地图的图幅号也在不断变动。

其中，国家基本比例尺地形图的图幅号更新速度较快。

为了更好地管理和使用地图，需要对新旧图幅进行转换。

本文将介绍国家基本比例尺地形图新旧图幅号转换的计算方法。

背景国家基本比例尺地形图是我国重要的地图资源，包含了我国全境的地形数据、地势信息、河流、湖泊等重要信息。

为了更好地管理和使用这些地图资源，需要对其图幅号进行统一管理，同时对新旧图幅进行转换。

随着地图的不断更新，图幅号也在不断变化，因此需要掌握转换计算方法，以确保使用最新且准确的地图信息。

计算方法国家基本比例尺地形图的图幅号包括旧版和新版两部分。

其中，旧版图幅号格式为“50W200S01”，新版图幅号格式为“50W200S.BMP”。

需要注意的是，新版图幅号中的数字“50W200S”与旧版图幅号的“50W200S01”表示的是同一区域。

旧版图幅号转新版图幅号如果需要将旧版图幅号转换成新版图幅号，可以按照以下方法进行计算：1.将旧版图幅号中的末尾“01”去掉，如“50W200S01”转化为“50W200S”；2.将去掉“01”的旧版图幅号前两位（“50W”）替换为“50”（即新版图幅号前两位），得到新版图幅号的前缀部分；3.将新版图幅号的前缀部分与旧版图幅号去掉“01”后的后缀部分（即“200S”）拼接，得到新版图幅号。

举个例子，如果需要将旧版图幅号“50W200S01”转换成新版图幅号，可以按照以下步骤进行计算：1.将旧版图幅号中的末尾“01”去掉，得到“50W200S”；2.将去掉“01”的旧版图幅号前两位（“50W”）替换为“50”，得到新版图幅号的前缀部分为“50S”；3.将新版图幅号的前缀部分“50S”与旧版图幅号去掉“01”后的后缀部分“200S”拼接，得到新版图幅号“50S200S.BMP”。

地形图的应用（国家基本比例尺地形图分幅与编号）大比例尺地形图是建筑工程规划设计和施工中的重要地形资料。

特别是在规划设计阶段，不仅要以地形图为底图，进行总平面的布设，而且还要根据需要，在地形图上进行一定的量算工作，以便因地制宜地进行合理的规划和设计。

1.按经纬线分幅的梯形分幅法(又称为国际分幅)2.按坐标格网分幅的矩形分幅法。

一．梯形分幅与编号1、1：100万比例尺图的分幅与编号按国际上的规定，1：100万的世界地图实行统一的分幅和编号。

即自赤道向北或向南分别按纬差4°分成横列，各列依次用A、B...V表示。

自经度180°开始起算，自西向东按经差6°分成纵行，各行依次用l、2 (60)表示。

每一幅图的编号由其所在的“横列一纵行”的代号组成。

例如北京某地的经度为东经118°24′20″，纬度为39°56′30″，则所在的1：100万比例尺图的图号为J—50。

2、1：10万比例尺图的分幅和编号将一幅1：100万的图，按经差30′，纬差20′分为144幅1：10万的图。

3、这三种比例尺图的分幅编号都是以1：10万比例尺因为基础的。

每幅1：10万的图，划分成4幅1：5万的图，分别在1：10万的图号后写上各自的代号A、B、C、D。

每幅1：5万的图又可分为4幅1：2．5万的图，分别以1、2、3、4编号。

每幅1：10万图分为64幅1：1万的图，分别以(1)、(2)、……(64)表示。

1：5000和1：2000比例尺图的分幅编号是在1：10000图的基础上进行的。

每幅1：10000 的图分为4幅1：5000的图，分别在1：10000的图号后面写上各自的代号a、b、c、d。

每幅1：5000的因又分成9幅1：2000的图，分别以l、2、……9表示，图幅的大小及编二、矩形分幅与编号大比例尺地形图大多采用矩形分幅法，它是按统一的直角坐标格网划分的。

采用矩形分幅时，大比例尺地形图的编号，一般采用图幅西南角坐标公里数编号法。

1:1万和1:5000比例尺地形图分幅最近就第一次全国水利普查工作，很多水利单位的朋友来向我询问1:1万新旧图幅号对应关系。

正好之前研究过，下面这篇新旧图幅号转换的方法当时好像是在黑龙江某个单位的网站上看到的，时间久了，具体记不清是黑龙江的那个单位了。

无论是否记得清，在此都表示衷心的感谢！或许这次能给咱们第一次全国水利普查带来很大的帮忙。

由于来不及一个一个挨着传，更不知还有哪位朋友需要，就只能把这篇新旧图幅号转换方法放于此，希望能帮助正需要的朋友。

但愿最近搞第一次全国水利普查需要此方法的朋友能在此找到。

希望给咱们第一次全国水利普查工作带来方便和帮助！1.1国家基本比例尺地形图为1/100万、1/50万、1/25万、1/10万、1/5万、1/2.5万、1/1万和1/5000。

基本比例尺地形图的分幅均以1/100万地形图为基础，按规定的纬差和经差划分图幅，并采用行列式编号，十位数，编码长度相同，编码系列统一为一个根部，便于计算机管理和检索。

国家基本比例尺及代码表:1.2 1/100万地形图的分幅和编号⑴自赤道走每纬差4为一行（向两极），依次A、B、C、D…V表示行号⑵从180°经线走自西向东每经差6为一列，依次1、2、3、4…60表示列号(见下图)如Ａ点所在1/100万地形图编号为J50，B点为L521.3 1/1万地形图分幅和编号以１/100万地形图为基础，划分96个横行和96个纵列，共9216幅地形图。

⑴行列号自左上至右和下排列⑵纬差6°/96=3′45″经差4°/96=2′30″见下图如上图晕线图幅编号为J50G089007⑶编号由所在1/100万图幅编号、比例尺代码、1/1万图幅所在行号和列号组成（十位数）××××××××××①②③④⑤①1/100万行号(字符码)②1/100万列号(数字码)③比例尺代码④1/1万行号(数字码)⑤1/1万列号(数字码)⑷图幅编号的求算例：某位置B=46°02′L＝129°16′。

国家基本比例尺地形图的分幅与编号及新旧图号的转换限于篇幅的关系，本文只介绍l：100000和l：10000地形图新旧图号的转换，其它比例尺与此基本相似。

1．1：100000地形图新、旧图号的转换(1)旧图号到新图号的转换。

由于l：100000地形图是把l：l000000地形图分为l2行和12列共计144幅而成的。

所以，l：100000的位置号码除以l2，如果能整除，商即为新图号中的行号，列号为l2；否则商加l为行号，余数为列号。

例：把1-49-108转换为新图号，108／12商为9，余数为0。

所以新图号为I49D009012。

(2)新图号到旧图号的转换。

由于l：100000地形图的新图号中的(行号-1)行为整行，所以(行号一1)×12+列号即为位置号码。

例：把I49G009012转换为旧图号，(9一1)×12+12=108，所以旧图号为1-49-108。

2．1：10000地形图新、旧图号的转换(1)旧图号到新图号的转换。

先根据l：100000图号中的位置号码算出行号和列号；然后把l：10000的位置号码中的数字除以8，如果能整除，商即为l：10000图幅在l：100000图幅中的行位置(m)，列位置(n)为8；否则商加l为m，余数为n。

由于一行l：100000的图幅可分为8行l：10000的图幅，一列l：100000的图幅可分为8列l：10000的图幅。

(行号一1)×8+m就是l：10000新图号中的行号；(列号一1)×8+n．就是l：10000新图号中的列号。

例：把I一49—108一(64)转换为新图号，108／12商为9，余数为0，求得行号=9，列号=l；64／8商为8，余数为0，求得m=8，n=8；从而新图号中的行号=(9一1)×8+8=72，新图号中的列号=(12一1)×8+8=96。

所以新图号为I49G072096。

(2)新图号到旧图号的转换例：把I49G072096转换为旧图号，72／8商为9，余数为0，求得行号=9，m=8；96／8商为l2，余数为0，求得列号=12，n=8；从而l：100000图号中的位置号码=(9一1)×12+12=108，1：10000图号中的位冠号码=(8一1)×8+8=64。

1∶1万地形图新、旧图幅号的手工换算方法在分析国家1∶1万地形图新、旧图幅号构成及其含义的基础上，阐述1∶1万地形图从旧图幅号换算成新图幅号、从新图幅号换算成旧图幅号时的手工换算方法与换算关系式，并用实例进行了验算。

关键词：新图幅号；旧图幅号；换算方法；1∶1万地形图国家技术监督局于1992年12月批准发布了《中华人民共和国国家标准国家基本比例尺地形图分幅和编号》（GB/T 13989-92）[1]，次年7月1日施行。

在实际使用中，将1993年以前按地形图分幅编号标准产生的地形图图幅号称为旧图幅号，1993年以后按新的国家基本比例尺地形图分幅和编号标准（即GB/T 13989-92）产生的地形图图幅号称为新图幅号。

现阶段，我国正在使用中的国家基本比例尺地形图，其图幅编号有新、旧之分，这给人们尤其是市县级以下基层生产单位专业技术人员带来了较大的障碍或困难，造成了使用中的不便。

《中华人民共和国国家标准国家基本比例尺地形图分幅和编号》（GB/T 13989-92）只是规范了新的图幅分幅与编号规则，并未给出我国国家基本比例尺地形图新、旧图幅号彼此间的换算关系；为解决新、旧图幅号之间的换算关系，我国的一些科技工作者从不同角度对此进行了探索研究。

笔者通过多渠道检索，查到17篇相关期刊论文[2-18]。

最早的关于地形图新旧图幅编号的换算研究文献发表于1997年，其中半数研究文献发表于近5年的相关科技期刊上，这些研究文献基本上是基于国家基本比例尺地形图的经纬度条件下，地形图分幅与图幅编号的新旧图幅号之间的换算，且多侧重于编程自动换算，以方便于科研或生产项目中批量操作管理，但满足不了基层生产单位专业技术人员在实际工作中遇到的少量或个别的只用手工即可进行的新旧图幅号便捷换算方法。

2009～2010年，笔者有幸参与广西新一轮森林资源规划设计调查（即二类资源调查）的部分县区的外、内业工作，尤其是内业制图工作，在工作中常遇到1∶1万地形图新、旧图幅号需要彼此间换算的问题，经过查阅相关规范、文献资料，反复探索研究，找到了适用于工作中遇到的少量或个别的可手工进行的新旧图幅号便捷换算方法，经验证，结果正确，便捷有效，现将研究成果系统整理出来，供业界同仁共享，方便工作。

大地坐标和经纬度之间的换算(整理实用)地质工作中常要对进行大地坐标转经纬度和经纬度换大地坐标,我写一下一般的过程,希望对大家有点帮助.一、大地座标-----→经纬度(地理坐标)1,输入大地坐标数据，格式为 Y空格X，输入到文本就行如下，原始的大地坐标由一8位的Y和一个7位的X组成，这组坐标数据中的Y的前两位为31，是分带号，一般使用的分带有三分带，六分带，这里的坐标是三分带的，记下Y前的这两位数，在原始数据中去除掉，现在数据变为，Y，6位，X，7位保存这个TXT的文本文件。

2，打开MAPGIS，启动坐标投影变形程序接下来选择投影转换>>>>用户文件投影转换点打开文件，打开刚才的大地坐标的文本文件，设置输入数据的格式，点击用户投影参数，并完成设置。

这里我们的大地座标为3度带的第31带，注意填好，坐标单位为米好了以为设置输出的格式，我们要求输出的是经纬度，点结果转换参数，完成以下设置我们输出的经纬度的单位应该是DDDMMMSS。

SS注意点写到文件，保存就大功告成了，（注意保存的文件一定要写上.TXT的后缀）下面是计算出的结果文件XP为经度，1234234.357就是123度42分34.357秒，YP为纬度，403950.225就是40度39分50.255秒（纬度没有最多90，所以没有三位数）~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 二、经纬度——→大地座标输入文件格式如下，这里面的数据前面的为经度，格式为ＤＤＤＭＭＳＳ，后面的为纬度，格式为ＤＤＭＭＳＳ接下来的转换过程和大地坐标转换一样，只要将刚才的用户转换参数和结果转换参数交换即可，要注意分带号的确定，如果你不知道分带号，就应该先计算分带号，算法是经度／３得到的整数为三度带的分带号经度／６得到的整数为六度带的分带号计算所得的结果格式如下其中的ＸＰ为地图上的Ｙ坐标，记得在前面加上带号，其中的ＹＰ为地图上的Ｘ坐标三、如何计算当地的中央经线及带号？子午线GPS要设置北京54或者西安80坐标系的时候，除了五个Delta参数需要输入外，还有一个中央经线参数需要输入，与Delta参数不同，中央经线是可以自己计算的，计算方法如下：1、六度带中央经线经度的计算：当地中央经线经度＝６°×当地带号－３°，例如：地形图上的横坐标为２０３４５，其所处的六度带的中央经线经度为：６°×２０－３°＝１１７°（适用于１∶２．５万和１∶５万地形图2、三度带中央经线经度的计算：中央经线经度＝３°×当地带号（适用于１∶１万地形图）。

( 1.景德镇市国土资源测绘勘察设计院景德镇330046)333000;2.江西省第二测绘院南昌摘要直接给出了1: 10000 地形图新旧图号的换算公式, 提出了计算要点和注意事项。

关键词图号经纬度换算根据国家基础测绘的分级管理原则, 1∶10000基础地理信息数据生产和数据库建设是省级政府的职责, 同时也是国家基础地理信息系统的重要组成部分, 还是更新1∶100 万、1∶25 万和1∶5 万数据库的主要数据源。

江西省测绘局依据国家测绘局“基础框架建设”及“数据江西”的规划要求, 正采用“3S”技术和现代数据通信技术进行测绘生产技术改造, 以数字高程模型(DE M)、数字正射影像图(DOM)、数字线划图(DLG)为代表的基础地理信息的生产已作为近段时期的主要任务, 逐步对全省的1∶10000 地形图进行更新。

其中:GA——1∶100万地形图图幅所在纬度带字符码所对应的数字码, 即新图号中的“H”;GB——1∶100地形图图幅所在经度带的数字码, 即新图号中的“50”:GC——1∶10000地形图在1∶100 万地形图图号后的行号, 即新图号中的“078”;GD——1∶10000地形图在1∶100 万地形图图号后的列号, 即新图号中的“034”; LXN——图幅西南图廓点的经度; BXN——图幅西南图廓点的纬度。

(2) 计算1∶10 万地形图在1∶100 万地形图图号后的数字码SU=1+INT((BXN- (GA- 1)×4°)/20′)SV=1+INT((LXN- (GB- 31)×6°) /30′)SW=SV+(12- SU)×12其中SW——1∶10地形图在1∶100 地形图图号后的数字码, 即旧图号中的“113”。

(3) 计算1∶10000 地形图关1∶10 万地形图图号后的数字码TU=(GA- 1)×4°+20′×INT ((B XN- (GA- 1)×4°) / 20′)TV=(GB- 31)×6°+30′×INT((LXN- (GB- 31)×6°) /30′)TF=1+INT((BXN- TU) /2′30″)TG=1+INT((LXN- TV) /3′45″)TW=TG+(8- TF)×8其中:TW——1∶10000 地形图在1∶10 万地形图图号后的数字码, 即旧图号中的“42”。

6种国家基本比例尺地形图新、旧图幅号的直接换算方法唐世斌【摘要】在解决1∶1万地形图新、旧图幅号的手工换算方法基础上,研究了我国其余6种国家基本比例尺(1∶50万～1∶5 000)地形图直接利用其图幅号进行新、旧图幅号间的相互换算问题,给出各比例尺地形图新、旧图幅号间直接换算的方法与关系式,并用实例进行了验算.【期刊名称】《林业调查规划》【年(卷),期】2013(038)004【总页数】9页(P5-13)【关键词】地形图;新图幅号;旧图幅号;换算方法;国家基本比例尺【作者】唐世斌【作者单位】广西大学林学院,广西南宁530004;广西大学行健文理学院工程与设计学部,广西南宁530005【正文语种】中文【中图分类】S711;S758.6;P217从1993年7月1日《中华人民共和国国家标准国家基本比例尺地形图分幅和编号》(GB/T 13989-92)［1］施行起，我国使用中的国家基本比例尺地形图其图幅编号有了新、旧图幅号之分;《国家基本比例尺地形图分幅和编号》(GB/T 13989-92)只是规范了新的图幅分幅与编号规则，而并未给出地形图新、旧图幅号彼此间的换算关系。

为解决新、旧图幅号之间的换算关系，我国的一些科技工作者从不同角度进行了探索研究，但基本上是基于地形图的经纬度条件下新旧图幅号之间的换算，且多侧重于编程自动换算，无法满足实际工作中的手工换算需要［2］。

笔者在《1∶1万地形图新、旧图幅号的手工换算方法》的基础上［2］，扩展思路，对我国其余6种国家基本比例尺地形图直接利用其图幅号进行新、旧图幅号间的相互换算问题进行了探究，找出了各比例尺地形图新、旧图幅号间的直接相互换算方法与关系，经检验便捷有效，现将研究成果整理出来，供业界同仁共享。

1 1∶50万地形图1.1 1∶50万地形图旧图幅号的构成及其含义我国基本比例尺地形图的图符编号都是在1∶100万地形图图幅号的基础上进行的，新、旧图幅号均是如此［3］;其具体的分幅和编号相关知识不是本文应述内容，请读者查阅相关规范、文献资料。

1:1万和1:5000比例尺地形图分幅最近就第一次全国水利普查工作，很多水利单位的朋友来向我询问1:1万新旧图幅号对应关系。

正好之前研究过，下面这篇新旧图幅号转换的方法当时好像是在黑龙江某个单位的网站上看到的，时间久了，具体记不清是黑龙江的那个单位了。

无论是否记得清，在此都表示衷心的感谢！或许这次能给咱们第一次全国水利普查带来很大的帮忙。

由于来不及一个一个挨着传，更不知还有哪位朋友需要，就只能把这篇新旧图幅号转换方法放于此，希望能帮助正需要的朋友。

但愿最近搞第一次全国水利普查需要此方法的朋友能在此找到。

希望给咱们第一次全国水利普查工作带来方便和帮助！1.1国家基本比例尺地形图为1/100万、1/50万、1/25万、1/10万、1/5万、1/2.5万、1/1万和1/5000。

基本比例尺地形图的分幅均以1/100万地形图为基础，按规定的纬差和经差划分图幅，并采用行列式编号，十位数，编码长度相同，编码系列统一为一个根部，便于计算机管理和检索。

国家基本比例尺及代码表:1.2 1/100万地形图的分幅和编号⑴自赤道走每纬差4为一行（向两极），依次A、B、C、D…V表示行号⑵从180°经线走自西向东每经差6为一列，依次1、2、3、4…60表示列号(见下图)如Ａ点所在1/100万地形图编号为J50，B点为L521.3 1/1万地形图分幅和编号以１/100万地形图为基础，划分96个横行和96个纵列，共9216幅地形图。

⑴行列号自左上至右和下排列⑵纬差6°/96=3′45″经差4°/96=2′30″见下图如上图晕线图幅编号为J50G089007⑶编号由所在1/100万图幅编号、比例尺代码、1/1万图幅所在行号和列号组成（十位数）××××××××××①②③④⑤①1/100万行号(字符码)②1/100万列号(数字码)③比例尺代码④1/1万行号(数字码)⑤1/1万列号(数字码)⑷图幅编号的求算例：某位置B=46°02′L＝129°16′。

国家地形图经纬度与图幅号的相互转换地形图的分幅编号及图幅号与经纬度相互转换1绪论1.1国家基本地形图的特点（1）具有统一的大地坐标系统和高程系统，我国国家基本地形图统一采用1980年“中国国家大地坐标系统”和“1985年国家高程基准”。

（2）具有完整的比例尺系列和分幅编号系统。

国家基本地形图按统一规定的经差和纬差进行分幅，每幅图的内图廓都是由经线和纬线构成，在国际百万分之一地图分幅编号基础上，建立了各级比例尺地形图的图幅编号系统。

（3）国家基本地形图是依据国家测绘管理部门统一定制的测量与编绘规范和《地形图图式》完成的，以便能够保证由各地方测绘部门分工完成的地形图，在质量和规格上的完全统一。

1.2国家基本地形图分幅的编号地形图分幅编号的方法很多，1990年底之前的国家基本地形图采用的是行列-自然序列的编号方法，我国称为原国家基本地形图编号法。

其中1:100万地形图采用的是行列式编号法，而1:50万以下七种比例尺地形图则采用行列-自然序列编号法。

1991年之后的国家基本地形图编号一律采用行列式的编号方法，称为新国家基本地形图编号法。

本文所介绍的就是1991年之后的新国家基本地形图编号法。

以1993-07-01实施的国家基本比例尺地形图分幅和编号标准（GB/T 13989-92）为依据。

1.3主要内容我国地形图的分幅以及地形图编号方法；我国地形图编号应用的公式以及图幅号与经纬度相互转换应用举例；2 我国基本比例尺地形图分幅和编号2.1我国基本比例尺地形图的分幅我国基本比例尺地形图均以1:1000 000地形图为基础，按规定的经差和纬差划分图幅。

1:1000 000地形图的分幅采用国际1：1000 000地图分幅标准。

每幅1：1000 000地形图的范围是经差6°、纬差4°；纬度60°~76°之间为经差12°、纬差4°；纬度76°~88°之间为经差24°、纬差4°（由于我国没有纬度60°以上的区域。

比例尺°′1:100万104910-493601：50万1049B 1210-49-23801：25万1049C 2410-49-(8)3801：20万10493610-49-(18)3801：10万1049D 61110-49-713801：5万1049E 112210-49-71-238101：2.5万1049F 224310-49-71-2-338101：1万1049G 448510-49-71-(29)3891:5千1049H 8716910-49-71-(29)-13811新 图 号旧 图 号纬 度 范根据一点的经纬度数据求所在地形图图幅的编号及经纬度范围蓝色数字AB 红色数字对应甲乙绿色数字对应a行号公式：c=4/△φ-[(φ/4)/△φ]，△φ为所在比例尺图幅的纬差列号公式：d=[(λ/6)/△λ]+1，△λ为所在比例尺图幅的经差,()表示取余。

1：100万图幅编号计算公式行号公式：a=[φ/4]+1,[]表示取整，φ为测区内任一点（或图幅西南角）纬度列号公式：b=[λ/6]+31，[]表示取整，λ为测区内任一点（或图幅西南角）经度1:100万图幅号之后的行列号计算色数字对应ABCDE..红色数字对应甲乙丙丁绿色数字对应abcd度范围经度范围″°′″°′″°′″0~400010800~114000~400011100~114000~3900112300~114000~3840011300~114000~3820011300~1133000~38200113150~1133000~38150113150~113223060~381230113150~113184515~381230113150~1131653。

图幅与经纬度之间的换算⼀、基本知识了解/*** 根据图幅计算经纬度* 1:100万 A 经6度，纬4度; 1:50万 B 经3度，纬2度; 1:25万 C 经1度30分，纬1度; 1:10万 D 经30分纬20分;* 1:5万 E 经15分纬10分; 1：2.5万 F 经3分45秒纬2分30秒;1：1万 G 经1分52.5秒纬1分15秒;* 1：5000 H; 1:2000 J; 1:1000 K; 1:500 L; 1:20万 R;* 1:150万 S; 1:250万 T; 1:400万 U; 1:500万 V; 1:600万 W; 1:800万 X;1:1000万 Y* mappableUnit=J50D002002*/public static MappableUnitHelper doComputerByMappableUnit(String mappableUnitStr){// 分步计算// y=h*y1-(h1-1)*y2,x=(l-31)*x1+(l1-1)*x2// ,y代表纬度，y1代表⽐例尺100万地图的纬差，x代表经度，y1代表纬差，x1代表经差,y2代表纬差，x2代表经差 //下代表⽐例尺100万地图的经差。

[]为商取整符号，如9.3则取9，h代表图号⾏号，l代表图号列号，h1代表// h1代表图号的地形图在基础图幅内所位于的⾏号，l1代表图号的地形图在基础图幅内所位于的列号if(StringUtils.isBlank(mappableUnitStr) ||mappableUnitStr.length() != 10 ){("GeotoolsUtil.doComputerByMappableUnit.mappableUnit is not valid");return null ;}//图号⾏号String h = mappableUnitStr.substring(0,1) ;//图号列号double l =Double.valueOf(mappableUnitStr.substring(1,3)) ;//⽐例尺代号String mappableCode =mappableUnitStr.substring(3,4) ;//基础图幅内位于的⾏号double h1 =Double.valueOf(mappableUnitStr.substring(4,7)) ;//基础图幅内位于的列号double l1 =Double.valueOf(mappableUnitStr.substring(7,10)) ;//⽐例尺Object bilichi = MappableUnitHelper.getMappableUnitMap().get(mappableCode) ;if(bilichi == null ){return null ;}String xy = MappableUnitHelper.getxYDvalueMap().get(bilichi).toString() ;double x2 = Double.valueOf(xy.toString().split(",")[0]) ;double y2 = Double.valueOf(xy.toString().split(",")[1]) ;Object doubleLine = MappableUnitHelper.getMappableUnitLineMap().get(h) ;if(doubleLine == null){return null ;}double hDouble = Double.valueOf(doubleLine.toString()) ;double ymax=(hDouble*14400-(h1-1)*y2)/3600 ;double xmin=((l-31)*21600+(l1-1)*x2)/3600;double ymin= ymax-y2/3600;double xmax=xmin+x2/3600;MappableUnitHelper mappableUnit = new MappableUnitHelper() ;mappableUnit.setXmin(xmin);mappableUnit.setXmax(xmax);mappableUnit.setYmin(ymin);mappableUnit.setYmax(ymax);return mappableUnit;}/*** 根据经纬度计算图幅* 1:100万 A 经6度，纬4度; 1:50万 B 经3度，纬2度; 1:25万 C 经1度30分，纬1度; 1:10万 D 经30分纬20分;* 1:5万 E 经15分纬10分; 1：2.5万 F 经3分45秒纬2分30秒;1：1万 G 经1分52.5秒纬1分15秒;* 1：5000 H; 1:2000 J; 1:1000 K; 1:500 L; 1:20万 R;* 1:150万 S; 1:250万 T; 1:400万 U; 1:500万 V; 1:600万 W; 1:800万 X;1:1000万 Y**/public static String changeToMappale(double x,double y,String scale){// 分步计算// h=[y/y1]+1;l=[x/x1]+31;h1=y1/y2-[(y/y1)/y2];l1=[(x/x1)/x2]+1;// ,y代表纬度，y1代表⽐例尺100万地图的纬差，x代表经度，y1代表纬差(1:100万),值为4度，x1代表经差(1:100万)，值为6度,y2代表纬差，x2代表经差 //下代表⽐例尺100万地图的经差。

大地坐标和经纬度之间的换算大地坐标-----→经纬度(地理坐标)1、输入大地坐标数据，格式为Y空格X，输入到文本就行如下，原始的大地坐标由一8位的Y和一个7位的X组成，这组坐标数据中的Y的前两位为31，是分带号，一般使用的分带有三分带，六分带，这里的坐标是三分带的，记下Y前的这两位数，在原始数据中去除掉，现在数据变为，Y，6位，X，7位保存这个TXT的文本文件。

2，打开MAPGIS，启动坐标投影变形程序接下来选择投影转换>>>>用户文件投影转换点打开文件，打开刚才的大地坐标的文本文件，设置输入数据的格式，点击用户投影参数，并完成设置。

这里我们的大地座标为3度带的第31带，注意填好，坐标单位为米好了以为设置输出的格式，我们要求输出的是经纬度，点结果转换参数，完成以下设置这个参数设置里第一行有一个很重要的参数，椭球参数这一项，很多博友没有到正确的结果就是因为该项参数没有设定，中国地图一般只有两种，一种是北京54，还有就是80黄海正确设定好这个就能得到正确的结果啦，笔者用的是80黄海的参考椭球，不过中国大部分是用54坐标系的！！！我们输出的经纬度的单位应该是DDDMMMSS。

SS注意点写到文件，保存就大功告成了，注意保存的文件要写上.TXT的后缀下面是计算出的结果文件XP为经度，1234234。

357就是123度42分34。

357秒，YP为纬度，403950。

225就是40度39分50。

255秒（纬度最多90，所以没有三位数）~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 经纬度——→大地座标输入文件格式如下，这里面的数据前面的为经度，格式为ＤＤＤＭＭＳＳ，后面的为纬度，格式为ＤＤＭＭＳＳ接下来的转换过程和大地坐标转换一样，只要将刚才的用户转换参数和结果转换参数交换即可，要注意分带号的确定，如果你不知道分带号，就应该先计算分带号，算法是经度／３得到的整数为三度带的分带号经度／６得到的整数为六度带的分带号计算所得的结果格式如下其中的ＸＰ为地图上的Ｙ坐标，记得在前面加上带号，其中的ＹＰ为地图上的Ｘ坐标。

地形图的分幅编号及图幅号与经纬度相互转换1绪论1.1国家基本地形图的特点（1）具有统一的大地坐标系统和高程系统，我国国家基本地形图统一采用1980年“中国国家大地坐标系统”和“1985年国家高程基准”。

（2）具有完整的比例尺系列和分幅编号系统。

国家基本地形图按统一规定的经差和纬差进行分幅，每幅图的内图廓都是由经线和纬线构成，在国际百万分之一地图分幅编号基础上，建立了各级比例尺地形图的图幅编号系统。

（3）国家基本地形图是依据国家测绘管理部门统一定制的测量与编绘规范和《地形图图式》完成的，以便能够保证由各地方测绘部门分工完成的地形图，在质量和规格上的完全统一。

1.2国家基本地形图分幅的编号地形图分幅编号的方法很多，1990年底之前的国家基本地形图采用的是行列-自然序列的编号方法，我国称为原国家基本地形图编号法。

其中1:100万地形图采用的是行列式编号法，而1:50万以下七种比例尺地形图则采用行列-自然序列编号法。

1991年之后的国家基本地形图编号一律采用行列式的编号方法，称为新国家基本地形图编号法。

本文所介绍的就是1991年之后的新国家基本地形图编号法。

以1993-07-01实施的国家基本比例尺地形图分幅和编号标准（GB/T 13989-92）为依据。

1.3主要内容我国地形图的分幅以及地形图编号方法；我国地形图编号应用的公式以及图幅号与经纬度相互转换应用举例；2 我国基本比例尺地形图分幅和编号2.1我国基本比例尺地形图的分幅我国基本比例尺地形图均以1:1000 000地形图为基础，按规定的经差和纬差划分图幅。

1:1000 000地形图的分幅采用国际1：1000 000地图分幅标准。

每幅1：1000 000地形图的范围是经差6°、纬差4°；纬度60°~76°之间为经差12°、纬差4°；纬度76°~88°之间为经差24°、纬差4°（由于我国没有纬度60°以上的区域。

地形图分幅及图幅号转换林辉;唐可平;王苗根;厉增新【摘要】地形图是分图幅进行管理的,林业生产常用地形图采用经纬网梯形分幅法,不同比例尺的地形图通过不同经度差、纬度差进行分幅,不同比例尺地形图在分幅时存在相关关系,图幅号之间也存在关联,可以通过数学公式进行换算.分幅方法是建立在经纬度基础上,经纬度与图幅号之间可以换算.存在新旧两套地形图编号体系,分幅方法基本相同,但编号差异较大,新旧分幅号相互之间可以换算.通过对地形图分幅及编号的介绍,归纳总结相互之间关系,采用数学公式进行描述,具有较强的可操作性.【期刊名称】《华东森林经理》【年(卷),期】2011(025)001【总页数】4页(P59-62)【关键词】地形图;图幅号;经纬度;西安80坐标系;北京54坐标系【作者】林辉;唐可平;王苗根;厉增新【作者单位】国家林业局华东林业调查规划院,浙江金华,3210012;国家林业局华东林业调查规划院,浙江金华,3210012;浙江省上虞市丰惠镇农救信息中心,312362;上虞市丁宅乡农技信息中心,312373【正文语种】中文【中图分类】TP391地形图是林业工作重要的基础资料，因不同空间定位精度需求有多种比例尺的地形图。

地形图是分图幅进行管理的，每幅地形图都有图幅编号，图幅编号是地形图组织管理的标识号，通过图幅号可以快速查找地形图，也可以通过图幅号查询相邻的地形图。

新旧图幅号编码规则不同，但可以通过公式进行换算。

1 地形图分幅地形图的分幅方法有两种：一种是经纬网梯形分幅法或国际分幅法；另一种是坐标格网正方形或矩形分幅法。

前者用于国家基本比例尺地形图，后者用于工程建设大比例尺地形图。

目前林业生产常用地形图是1∶25 0000、1∶10 0000、1∶50000、1∶10000，其分幅采用经纬网梯形分幅法。

图幅编号有两种，其一是基于北京54坐标系的地形图图幅编号，其二是基于西安80坐标系的地形图图幅编号。

前者称为旧图幅号，后者称为新图幅号。

比例尺°′1:100万104910-493601：50万1049B 1210-49-23801：25万1049C 2410-49-(8)3801：20万10493610-49-(18)3801：10万1049D 61110-49-713801：5万1049E 112210-49-71-238101：2.5万1049F 224310-49-71-2-338101：1万1049G 448510-49-71-(29)3891:5千1049H 8716910-49-71-(29)-13811新 图 号旧 图 号纬 度 范根据一点的经纬度数据求所在地形图图幅的编号及经纬度范围蓝色数字AB 红色数字对应甲乙绿色数字对应a行号公式：c=4/△φ-[(φ/4)/△φ]，△φ为所在比例尺图幅的纬差列号公式：d=[(λ/6)/△λ]+1，△λ为所在比例尺图幅的经差,()表示取余。

1：100万图幅编号计算公式行号公式：a=[φ/4]+1,[]表示取整，φ为测区内任一点（或图幅西南角）纬度列号公式：b=[λ/6]+31，[]表示取整，λ为测区内任一点（或图幅西南角）经度1:100万图幅号之后的行列号计算色数字对应ABCDE..红色数字对应甲乙丙丁绿色数字对应abcd度范围经度范围″°′″°′″°′″0~400010800~114000~400011100~114000~3900112300~114000~3840011300~114000~3820011300~1133000~38200113150~1133000~38150113150~113223060~381230113150~113184515~381230113150~1131653。

Visual Basic 环境下经纬度与图符号转换程序设计该应用程序包括两大部分功能，即1、已知经纬度而求地形图的图幅号；2、已知地形图图幅号求该图幅西南图廓点的经纬度。

如果将两部分功能都放在一个用户界面上，将会产生过多的控件而使窗口显得凌乱，影响用户界面的美观；所以，为了使界面看起来简洁美观，操作简单快捷，将程序设计为多重窗口（一个主窗口form1，两个子窗口form2、form3）模式，通过主窗口的操作来打开其子窗口。

下面就按照Visual Basic 开发应用程序的步骤来详细介绍该程序设计的过程：1 建立用户界面，设置窗口和控件属性Form1（经纬度与图符号相互转换）的界面设计如下：图3.1Form1窗体和控件的属性设置：Form2（已知经纬度求图幅号）的界面设计如下：图3.2Form2窗体和控件的属性设置：表 3.2Form3（已知图幅号求经纬度）界面设计如下：图3.3Form3窗体和控件的属性设置：表 3.32 编写程序代码2.1设计思路（1）Form1（图3.1）该窗口的主要作用是来选择计算的模式，通过选择Option1或者Option2，单击“确定”按钮来打开相应的窗口进行运算。

单击“返回”按钮就会结束程序的运行。

另外，单击“关于”按钮，会弹出一个消息框，来提示用户本软件的使用范围。

如果没有对模式进行选择就单击“确定”按钮，应该弹出一个提示性质的消息框，提醒用户选择需要的计算模式。

（2）Form2 (图3.2)该窗口是已知经纬度来求图幅号的运算界面，要求用户输入图幅某点或者图幅西南图廓点的经纬度，来计算（转换）求得该图幅的图幅编号。

用户在相应文本框输入经纬度，并且选择比例尺之后，单击“转换”按钮，就会在下面的文本框内显示运算结果，如果没有按照要求输入经纬度（漏输或者超范围输入）或者没有选择比例尺，都应该弹出提示性消息框，提示用户正确输入经纬度，选择比例尺。

单击“输入说明”，会有提示消息框，说明输入数据的范围和格式，从而保证计算结果的正确性。