林輝 唐可平 王苗根 厲增新

（1國家林業(yè)局華東林業(yè)調(diào)查規(guī)劃院 浙江金華 3210012；2 浙江省上虞市豐惠鎮(zhèn)農(nóng)技信息中心 312362；3上虞市丁宅鄉(xiāng)農(nóng)技信息中心 312373）

地形圖是林業(yè)工作重要的基礎(chǔ)資料，因不同空間定位精度需求有多種比例尺的地形圖。地形圖是分圖幅進(jìn)行管理的，每幅地形圖都有圖幅編號，圖幅編號是地形圖組織管理的標(biāo)識號，通過圖幅號可以快速查找地形圖，也可以通過圖幅號查詢相鄰的地形圖。新舊圖幅號編碼規(guī)則不同，但可以通過公式進(jìn)行換算。

1 地形圖分幅

地形圖的分幅方法有兩種：一種是經(jīng)緯網(wǎng)梯形分幅法或國際分幅法；另一種是坐標(biāo)格網(wǎng)正方形或矩形分幅法。前者用于國家基本比例尺地形圖，后者用于工程建設(shè)大比例尺地形圖。

目前林業(yè)生產(chǎn)常用地形圖是 1∶25 0000、1∶10 0000、1∶50000、1∶10000，其分幅采用經(jīng)緯網(wǎng)梯形分幅法。圖幅編號有兩種，其一是基于北京54坐標(biāo)系的地形圖圖幅編號，其二是基于西安80坐標(biāo)系的地形圖圖幅編號。前者稱為舊圖幅號，后者稱為新圖幅號。新舊圖幅的分幅原則基本一致，但因采用的坐標(biāo)系不同，橢球體不同，在實(shí)際位置上略有差異。

1.1 基于北京54坐標(biāo)系的地形圖分幅及編號

1∶100 0000地形圖分幅和編號是采用國際標(biāo)準(zhǔn)分幅的經(jīng)差6°、緯差4°為一幅圖。從赤道起向北或向南至緯度88°止，按緯差每4°為一劃分共22個橫行，依次用A、B、……、V表示；從經(jīng)度180°起向東按經(jīng)差每6°為一劃分，全球共劃分為 60縱列，依次用 1、2、……、60表示。每幅圖的編號由該圖幅所在的“列號——行號”組成。如：J-50。

1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖都是在1∶100 0000地形圖的基礎(chǔ)上進(jìn)行分幅編號。

一幅1∶100 0000的圖以經(jīng)差3°、緯差2°劃分為4幅1∶50 0000的圖，分別以A、B、C、D表示。將1∶100 0000圖幅的編號加上代碼，即為該圖幅的編號，如：J-50-A。

一幅1∶100 0000的圖以經(jīng)差1°30′、緯差1°劃分為16幅1∶250000的圖，分別用[1]、[2]、……、[16]表示。將1∶100 0000圖幅的編號加上代碼，即為該圖幅的編號，如：J-50-[1]。

一幅1∶100 0000的圖以經(jīng)差30′、緯差20′劃分為144幅1∶10 0000的圖，分別用1、2、……、144表示。將1∶100 0000圖幅的編號加上代碼，即為該圖幅的編號，如：J-50-1。

1∶50000、1∶10000比例尺地形圖都是在1∶10 0000地形圖的基礎(chǔ)上進(jìn)行分幅編號。

一幅1∶10 0000的圖以經(jīng)差15′、緯差10′劃分為4幅1∶50000的圖，分別以甲、乙、丙、丁表示。將1∶100000圖幅的編號加上代碼，即為該圖幅的編號，如：J-50-1-甲。

一幅1∶10 0000的圖以經(jīng)差3′45″、緯差2′30″劃分為64幅1∶10000的圖，分別以（1）、（2）、……、（64）表示。將1∶10 0000圖幅的編號加上代碼，即為該圖幅的編號，如：J-50-1-（1）。

1.2 基于西安80坐標(biāo)系的地形圖分幅及編號

分幅方法：1∶100 0000地形圖的分幅標(biāo)準(zhǔn)仍按國際分幅法進(jìn)行（舊分幅相同）。其余比例尺的分幅均以 1∶100 0000地形圖為基礎(chǔ)，按照各比例尺地形圖相應(yīng)經(jīng)差、緯差劃分圖幅（經(jīng)差、緯差與舊分幅相同）。

相應(yīng)比例尺經(jīng)緯差表

1∶100 0000圖幅的編號，由圖幅所在的“行號列號”組成。與國際編號基本相同，但行與列的稱謂相反。如北京所在1∶100萬圖幅編號為J50。

1∶50 0000至1∶5000圖幅的編號，由圖幅所在的1∶100 0000圖，行號（字符碼）1位，列號（數(shù)字碼）2位，比例尺代碼1位，該圖幅行號（數(shù)字碼）3位，列號（數(shù)字碼）3位，共10位代碼組成。

2 地形圖圖幅號及經(jīng)緯度換算

地形圖分幅是建立在經(jīng)緯度坐標(biāo)之上的，所以地形圖圖幅號與經(jīng)緯度存在著必然的關(guān)聯(lián)性，可以通過公式進(jìn)行描述。

2.1 舊圖幅號及經(jīng)緯度換算

2.1.1 已知圖幅內(nèi)某點(diǎn)的經(jīng)緯度求圖幅號

1∶1 00 0000地形圖圖幅編號:

行號：chr(int(緯度 / 4 )+ asc(“A”))

列號：int(經(jīng)度 / 6 ) + 31

采用VB函數(shù)，chr：數(shù)字轉(zhuǎn)換為字符，asc：字符轉(zhuǎn)換ASCII碼，int：取整。

1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖在1∶100 0000地形圖圖幅號后編號:

行號：4/緯差 - int((緯度 mod 4) /緯差)

列號：int（(經(jīng)度 mod 6) / 經(jīng)差）+ 1

編號：（行號－1 ）* 列圖幅數(shù) + 列號

Mod：取余。

列圖幅數(shù)：每幅1∶100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

1∶50000、1∶25000、1∶10000比例尺地形圖在1∶10 0000地形圖圖幅號后的編號:

行號：(4/緯差 - int((緯度 mod 4)/緯差)－1)mod 列圖幅數(shù) + 1

列號：（int（(經(jīng)度 mod 6)/經(jīng)差） mod列圖幅數(shù) + 1

編號：（行號-1）*列圖幅數(shù) + 列號

列圖幅數(shù)：每幅1∶10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

2.1.2 已知圖幅號求西南角點(diǎn)的經(jīng)緯度

1∶100 0000、1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖經(jīng)緯度：

經(jīng)度：（經(jīng)度帶 - 31）* 6 +（(編號 - 1) mod 列圖幅數(shù)）* 經(jīng)差

緯度：（asc（緯度帶）- asc（“A”）-1 ）* 4 +

(4 / 緯差–int(（編號 - 1 ) / 列圖幅數(shù)）)* 緯差

列圖幅數(shù)：每幅1∶100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

1∶50000、1∶25000、1∶10000比例尺地形圖的經(jīng)緯度:

經(jīng)度：（經(jīng)度帶 – 31）* 6 + （10萬編號 mod 12 -1 ）* 10萬經(jīng)差 +（(編號 - 1) mod列圖幅數(shù)）* 經(jīng)差

緯度：（asc（緯度帶）- asc（“A”）- 1 ）* 4 + （4 / 10萬緯差–int(（10萬編號 - 1 ) / 12） ) * 10萬緯差 +（10萬緯差 / 緯差–int((編號 – 1) / 列圖幅數(shù)) * 緯差

列圖幅數(shù)：每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

2.2 新圖幅號及經(jīng)緯度換算

2.2.1 已知圖幅內(nèi)某點(diǎn)的經(jīng)緯度求圖幅號

1:100 0000地形圖圖幅編號:

行號：chr(int(緯度 / 4 )+ asc(“A”))

列號：int(經(jīng)度 / 6 ) + 31

采用VB函數(shù)，chr：數(shù)字轉(zhuǎn)換為字符，asc：字符轉(zhuǎn)換ASCII碼，int：取整。

其它比例尺地形圖在1:100 0000地形圖圖幅號后的行、列號:

行號：4/緯差 - int((緯度mod4) / 緯差)

列號：int（(經(jīng)度 mod 6) / 經(jīng)差）+ 1 Mod：取余。

2.2.2 已知圖幅號求西南圖廓點(diǎn)的經(jīng)緯度

經(jīng)度：（經(jīng)度帶–31）* 6 +（列號-1）* 經(jīng)差

緯度：（asc（緯度帶）- asc（“A”）- 1 ）* 4 + (4 / 緯差–行號 ) * 緯差

3 新舊圖幅號換算

1∶10 00000地形圖圖幅編號在新舊體系中是一致的，不需∶換算。

3.1 新圖幅號轉(zhuǎn)換為舊圖幅號

1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖新圖幅號轉(zhuǎn)換為舊圖幅號：

（行號–1 ） * 列圖幅數(shù) + 列號

列圖幅數(shù)：每幅1:100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

1∶50000、1∶25000、1∶10000比例尺地形圖新圖幅號轉(zhuǎn)換為舊圖幅號：

所在 1∶10 0000比例尺地形圖編號：

int(（行號 – 1）/ 列圖幅數(shù)) * 12 + int((列號-1) / 列圖幅數(shù))

相應(yīng)比例尺地形圖編號：

(（行號 – 1）mod 列圖幅數(shù)) *列圖幅數(shù)+ (列號-1) mod 列圖幅數(shù) + 1

列圖幅數(shù)：每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

3.2 舊圖幅號轉(zhuǎn)換為新圖幅號

1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖舊圖幅號轉(zhuǎn)換為新圖幅號：

行號：int(（編號- 1 ) / 列圖幅數(shù)）+ 1

列號：(編號 - 1) mod 列圖幅數(shù) + 1

列圖幅數(shù)：每幅1:100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

1∶50000、1∶25000、1∶10000比例尺地形圖舊圖幅號轉(zhuǎn)換為新圖幅號：

行號：int（（10萬編號-1）/ 12 ）* 列圖幅數(shù) +int(（編號 - 1 ) / 列圖幅數(shù) ）+ 1

列號：(（10萬編號-1） mod 12 ) * 列圖幅數(shù) + （編號 - 1 ) mod 列圖幅數(shù) + 1

列圖幅數(shù)：每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開方。

4 新體系大比例尺地形圖圖幅號求算小比例尺地形圖圖幅號

行號：int(（行號 - 1 ) / 列圖幅數(shù)）+ 1

列號：int(（列號 - 1 ) / 列圖幅數(shù)）+ 1

列圖幅數(shù): 小比例尺地形圖包含大比例尺地形圖數(shù)量的開方。

5 結(jié)束語

我國幅員遼闊，疆域?qū)拸V，牽涉到的地形圖圖幅眾多，了解地形圖分幅原則，對實(shí)際工作具有重要意義。新舊圖幅號換算可以有多種方法，以上所列僅為直接換算的方法，也可以通過間接求經(jīng)緯度法換算。經(jīng)緯度及新舊圖幅號之間都有一定的相關(guān)關(guān)系，可以通過計算機(jī)編程實(shí)現(xiàn)相互之間的換算，對實(shí)際工作非常有用，也可以在地理信息系統(tǒng)中通過編程生成地形圖分幅圖，以備查用。

[1] 國家基本比例尺地形圖分幅和編號（GB/T 13989一92）.