林輝 唐可平 王苗根 厲增新
(1國(guó)家林業(yè)局華東林業(yè)調(diào)查規(guī)劃院 浙江金華 3210012;2 浙江省上虞市豐惠鎮(zhèn)農(nóng)技信息中心 312362;3上虞市丁宅鄉(xiāng)農(nóng)技信息中心 312373)
地形圖是林業(yè)工作重要的基礎(chǔ)資料,因不同空間定位精度需求有多種比例尺的地形圖。地形圖是分圖幅進(jìn)行管理的,每幅地形圖都有圖幅編號(hào),圖幅編號(hào)是地形圖組織管理的標(biāo)識(shí)號(hào),通過(guò)圖幅號(hào)可以快速查找地形圖,也可以通過(guò)圖幅號(hào)查詢(xún)相鄰的地形圖。新舊圖幅號(hào)編碼規(guī)則不同,但可以通過(guò)公式進(jìn)行換算。
地形圖的分幅方法有兩種:一種是經(jīng)緯網(wǎng)梯形分幅法或國(guó)際分幅法;另一種是坐標(biāo)格網(wǎng)正方形或矩形分幅法。前者用于國(guó)家基本比例尺地形圖,后者用于工程建設(shè)大比例尺地形圖。
目前林業(yè)生產(chǎn)常用地形圖是 1∶25 0000、1∶10 0000、1∶50000、1∶10000,其分幅采用經(jīng)緯網(wǎng)梯形分幅法。圖幅編號(hào)有兩種,其一是基于北京54坐標(biāo)系的地形圖圖幅編號(hào),其二是基于西安80坐標(biāo)系的地形圖圖幅編號(hào)。前者稱(chēng)為舊圖幅號(hào),后者稱(chēng)為新圖幅號(hào)。新舊圖幅的分幅原則基本一致,但因采用的坐標(biāo)系不同,橢球體不同,在實(shí)際位置上略有差異。
1.1 基于北京54坐標(biāo)系的地形圖分幅及編號(hào)
1∶100 0000地形圖分幅和編號(hào)是采用國(guó)際標(biāo)準(zhǔn)分幅的經(jīng)差6°、緯差4°為一幅圖。從赤道起向北或向南至緯度88°止,按緯差每4°為一劃分共22個(gè)橫行,依次用A、B、……、V表示;從經(jīng)度180°起向東按經(jīng)差每6°為一劃分,全球共劃分為 60縱列,依次用 1、2、……、60表示。每幅圖的編號(hào)由該圖幅所在的“列號(hào)——行號(hào)”組成。如:J-50。
1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖都是在1∶100 0000地形圖的基礎(chǔ)上進(jìn)行分幅編號(hào)。
一幅1∶100 0000的圖以經(jīng)差3°、緯差2°劃分為4幅1∶50 0000的圖,分別以A、B、C、D表示。將1∶100 0000圖幅的編號(hào)加上代碼,即為該圖幅的編號(hào),如:J-50-A。
一幅1∶100 0000的圖以經(jīng)差1°30′、緯差1°劃分為16幅1∶250000的圖,分別用[1]、[2]、……、[16]表示。將1∶100 0000圖幅的編號(hào)加上代碼,即為該圖幅的編號(hào),如:J-50-[1]。
一幅1∶100 0000的圖以經(jīng)差30′、緯差20′劃分為144幅1∶10 0000的圖,分別用1、2、……、144表示。將1∶100 0000圖幅的編號(hào)加上代碼,即為該圖幅的編號(hào),如:J-50-1。
1∶50000、1∶10000比例尺地形圖都是在1∶10 0000地形圖的基礎(chǔ)上進(jìn)行分幅編號(hào)。
一幅1∶10 0000的圖以經(jīng)差15′、緯差10′劃分為4幅1∶50000的圖,分別以甲、乙、丙、丁表示。將1∶100000圖幅的編號(hào)加上代碼,即為該圖幅的編號(hào),如:J-50-1-甲。
一幅1∶10 0000的圖以經(jīng)差3′45″、緯差2′30″劃分為64幅1∶10000的圖,分別以(1)、(2)、……、(64)表示。將1∶10 0000圖幅的編號(hào)加上代碼,即為該圖幅的編號(hào),如:J-50-1-(1)。
1.2 基于西安80坐標(biāo)系的地形圖分幅及編號(hào)
分幅方法:1∶100 0000地形圖的分幅標(biāo)準(zhǔn)仍按國(guó)際分幅法進(jìn)行(舊分幅相同)。其余比例尺的分幅均以 1∶100 0000地形圖為基礎(chǔ),按照各比例尺地形圖相應(yīng)經(jīng)差、緯差劃分圖幅(經(jīng)差、緯差與舊分幅相同)。
相應(yīng)比例尺經(jīng)緯差表
1∶100 0000圖幅的編號(hào),由圖幅所在的“行號(hào)列號(hào)”組成。與國(guó)際編號(hào)基本相同,但行與列的稱(chēng)謂相反。如北京所在1∶100萬(wàn)圖幅編號(hào)為J50。
1∶50 0000至1∶5000圖幅的編號(hào),由圖幅所在的1∶100 0000圖,行號(hào)(字符碼)1位,列號(hào)(數(shù)字碼)2位,比例尺代碼1位,該圖幅行號(hào)(數(shù)字碼)3位,列號(hào)(數(shù)字碼)3位,共10位代碼組成。
地形圖分幅是建立在經(jīng)緯度坐標(biāo)之上的,所以地形圖圖幅號(hào)與經(jīng)緯度存在著必然的關(guān)聯(lián)性,可以通過(guò)公式進(jìn)行描述。
2.1 舊圖幅號(hào)及經(jīng)緯度換算
2.1.1 已知圖幅內(nèi)某點(diǎn)的經(jīng)緯度求圖幅號(hào)
1∶1 00 0000地形圖圖幅編號(hào):
行號(hào):chr(int(緯度 / 4 )+ asc(“A”))
列號(hào):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地形圖圖幅號(hào)后編號(hào):
行號(hào):4/緯差 - int((緯度 mod 4) /緯差)
列號(hào):int((經(jīng)度 mod 6) / 經(jīng)差)+ 1
編號(hào):(行號(hào)-1 )* 列圖幅數(shù) + 列號(hào)
Mod:取余。
列圖幅數(shù):每幅1∶100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
1∶50000、1∶25000、1∶10000比例尺地形圖在1∶10 0000地形圖圖幅號(hào)后的編號(hào):
行號(hào):(4/緯差 - int((緯度 mod 4)/緯差)-1)mod 列圖幅數(shù) + 1
列號(hào):(int((經(jīng)度 mod 6)/經(jīng)差) mod列圖幅數(shù) + 1
編號(hào):(行號(hào)-1)*列圖幅數(shù) + 列號(hào)
列圖幅數(shù):每幅1∶10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
2.1.2 已知圖幅號(hào)求西南角點(diǎn)的經(jīng)緯度
1∶100 0000、1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖經(jīng)緯度:
經(jīng)度:(經(jīng)度帶 - 31)* 6 +((編號(hào) - 1) mod 列圖幅數(shù))* 經(jīng)差
緯度:(asc(緯度帶)- asc(“A”)-1 )* 4 +
(4 / 緯差–int((編號(hào) - 1 ) / 列圖幅數(shù)))* 緯差
列圖幅數(shù):每幅1∶100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
1∶50000、1∶25000、1∶10000比例尺地形圖的經(jīng)緯度:
經(jīng)度:(經(jīng)度帶 – 31)* 6 + (10萬(wàn)編號(hào) mod 12 -1 )* 10萬(wàn)經(jīng)差 +((編號(hào) - 1) mod列圖幅數(shù))* 經(jīng)差
緯度:(asc(緯度帶)- asc(“A”)- 1 )* 4 + (4 / 10萬(wàn)緯差–int((10萬(wàn)編號(hào) - 1 ) / 12) ) * 10萬(wàn)緯差 +(10萬(wàn)緯差 / 緯差–int((編號(hào) – 1) / 列圖幅數(shù)) * 緯差
列圖幅數(shù):每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
2.2 新圖幅號(hào)及經(jīng)緯度換算
2.2.1 已知圖幅內(nèi)某點(diǎn)的經(jīng)緯度求圖幅號(hào)
1:100 0000地形圖圖幅編號(hào):
行號(hào):chr(int(緯度 / 4 )+ asc(“A”))
列號(hào):int(經(jīng)度 / 6 ) + 31
采用VB函數(shù),chr:數(shù)字轉(zhuǎn)換為字符,asc:字符轉(zhuǎn)換ASCII碼,int:取整。
其它比例尺地形圖在1:100 0000地形圖圖幅號(hào)后的行、列號(hào):
行號(hào):4/緯差 - int((緯度mod4) / 緯差)
列號(hào):int((經(jīng)度 mod 6) / 經(jīng)差)+ 1 Mod:取余。
2.2.2 已知圖幅號(hào)求西南圖廓點(diǎn)的經(jīng)緯度
經(jīng)度:(經(jīng)度帶–31)* 6 +(列號(hào)-1)* 經(jīng)差
緯度:(asc(緯度帶)- asc(“A”)- 1 )* 4 + (4 / 緯差–行號(hào) ) * 緯差
1∶10 00000地形圖圖幅編號(hào)在新舊體系中是一致的,不需∶換算。
3.1 新圖幅號(hào)轉(zhuǎn)換為舊圖幅號(hào)
1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖新圖幅號(hào)轉(zhuǎn)換為舊圖幅號(hào):
(行號(hào)–1 ) * 列圖幅數(shù) + 列號(hào)
列圖幅數(shù):每幅1:100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
1∶50000、1∶25000、1∶10000比例尺地形圖新圖幅號(hào)轉(zhuǎn)換為舊圖幅號(hào):
所在 1∶10 0000比例尺地形圖編號(hào):
int((行號(hào) – 1)/ 列圖幅數(shù)) * 12 + int((列號(hào)-1) / 列圖幅數(shù))
相應(yīng)比例尺地形圖編號(hào):
((行號(hào) – 1)mod 列圖幅數(shù)) *列圖幅數(shù)+ (列號(hào)-1) mod 列圖幅數(shù) + 1
列圖幅數(shù):每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
3.2 舊圖幅號(hào)轉(zhuǎn)換為新圖幅號(hào)
1∶50 0000、1∶25 0000、1∶10 0000比例尺地形圖舊圖幅號(hào)轉(zhuǎn)換為新圖幅號(hào):
行號(hào):int((編號(hào)- 1 ) / 列圖幅數(shù))+ 1
列號(hào):(編號(hào) - 1) mod 列圖幅數(shù) + 1
列圖幅數(shù):每幅1:100 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
1∶50000、1∶25000、1∶10000比例尺地形圖舊圖幅號(hào)轉(zhuǎn)換為新圖幅號(hào):
行號(hào):int((10萬(wàn)編號(hào)-1)/ 12 )* 列圖幅數(shù) +int((編號(hào) - 1 ) / 列圖幅數(shù) )+ 1
列號(hào):((10萬(wàn)編號(hào)-1) mod 12 ) * 列圖幅數(shù) + (編號(hào) - 1 ) mod 列圖幅數(shù) + 1
列圖幅數(shù):每幅1:10 0000地形圖包含相應(yīng)比例尺地形圖圖幅數(shù)量的開(kāi)方。
行號(hào):int((行號(hào) - 1 ) / 列圖幅數(shù))+ 1
列號(hào):int((列號(hào) - 1 ) / 列圖幅數(shù))+ 1
列圖幅數(shù): 小比例尺地形圖包含大比例尺地形圖數(shù)量的開(kāi)方。
我國(guó)幅員遼闊,疆域?qū)拸V,牽涉到的地形圖圖幅眾多,了解地形圖分幅原則,對(duì)實(shí)際工作具有重要意義。新舊圖幅號(hào)換算可以有多種方法,以上所列僅為直接換算的方法,也可以通過(guò)間接求經(jīng)緯度法換算。經(jīng)緯度及新舊圖幅號(hào)之間都有一定的相關(guān)關(guān)系,可以通過(guò)計(jì)算機(jī)編程實(shí)現(xiàn)相互之間的換算,對(duì)實(shí)際工作非常有用,也可以在地理信息系統(tǒng)中通過(guò)編程生成地形圖分幅圖,以備查用。
[1] 國(guó)家基本比例尺地形圖分幅和編號(hào)(GB/T 13989一92).