冷繼全
(貴州省第一測繪院,貴州貴陽550025)
淺析基本比例尺新舊圖幅號批量轉換方法
冷繼全
(貴州省第一測繪院,貴州貴陽550025)
簡要介紹國家基本比例尺新舊圖幅號的編號方法及換算關系,采用VB、Excel VBA等高級程序語言,編程實現新舊圖幅號之間的批量轉換。結合在第二次土地調查項目中的具體應用,提高工作效益。
新舊圖幅號;編號方法;換算關系;工作效益
在第二次全國土地調查工作之中,經常會查找而第一次土地詳查的歷史資料,而第一次土地詳查的基礎資料編號為舊圖幅號,第二次全國土地調查則采用新圖幅號,兩者之間存在一定的對應關系。在時間緊、任務重的情況下,對資料的充分利用,顯得尤為重要。本文通過對國家基本比例尺新舊圖幅號的編號方法進行分析,利用VB、Excel VBA等編程語言,對基本比例尺新舊圖號的換算關系進行程序化,達到新舊圖號的批量轉換功能,以提高工作效率。
1.舊圖號的編號方法
基本比例尺是在1∶1 00 000的基礎上,按不同的經差和緯差,形成不同比例尺的基本圖幅號。首先,根據大地坐標(B、L),利用如下公式式中,B為緯度,以度為單位;L為經度,以度為單位;a表示緯度數字碼;b表示經度數字碼;[]表示取整。
其次,再根據不同基本比例尺分幅之間經差ΔL和緯差ΔB,利用如下公式
式中,ΔB表示相應基本比例尺圖幅的緯差;ΔL表示相應基本比例尺圖幅的經差;[]表示商數取整; ()表式商數取余;c表示基本比例尺圖幅的行號(數字碼);d表示基本比例尺圖幅的列號(數字碼)。
在1∶1 00 000的基礎上,通過式(1)求出字符碼和數字碼,用式(2)求出行列數字碼,將所得的行列數字碼對應為相應的字符碼,結合基本比例尺之間的分幅對應關系,就可得到相應的基本圖號。
2.新圖號的編號方法
在新圖幅號的編號中,各基本比例尺地圖的比例尺代碼用英文字母“A、B、C、D、E、F、G、H”表示,具體如表1所示。
表1
通過圖廓西南角點的大地坐標(B、L)值,再根據不同基本比例尺分幅之間經差ΔL和緯差ΔB,按式(1)、式(2),求出1∶1 000 000字符碼、數字碼和行列數字碼,新圖幅號采用10位編號組成(行號+列號+比例尺代碼+行號+列號)。
通過對上述基本比例尺分幅編號規(guī)則的了解,在1∶1 000 000的基礎上,分出1∶500 000、1∶200 000、1∶100 000基本比例尺;再以1∶100 000為基礎,分出1∶50 000、1∶10 000基本比例尺。根據各比例尺之間的行列關系,有如下各基本比例尺新舊圖幅號換算公式。
1.1 ∶500 000舊到新換算公式
由于1幅1∶1 000 000圖幅,可分為4幅1∶500000圖幅,按2行2列組成甲、乙、丙、丁,于是有下列公式
式中,H50為1∶500 000新圖號的行號;L50為1∶500 000新圖號的列號;X50為1∶500 000舊圖幅編號中相應比例尺地形圖的圖幅代碼值,甲、乙、丙、丁的代碼值為1、2、3、4;[]表示小數取整;()表示小數取余。
例如:一幅1∶500 000基本比例尺圖號為G—48—丙,由式(3)可知
則
取3位,不足3位前面用“0”補足,所以新圖號為G48B002001。
2.1 ∶250 000舊到新換算公式
同樣 1幅1∶1 000 000圖幅,可分為 16幅1∶250 000圖幅,按4行4列組成[1]、[2]、…、[16]編號,于是有下列公式
式中,H25為 1∶250 000新圖號的行號;L25為1∶250 000新圖號的列號;X25為1∶250 000舊圖幅編號中相應比例尺地形圖的圖幅代碼值;[]表示小數取整;()表示小數取余。
例如:一幅1∶250 000基本比例尺圖號為G—48—[10],由式(4)可知
則
取3位,不足3位前面用“0”補足,所以新圖號為G48C003002。
3.1 ∶100 000舊到新換算公式
同理,1幅1∶1 000 000圖幅可分為 144幅 1∶100 000圖幅,按12行12列組成1、2、…、144編號,于是有下列公式
式中,H10為 1∶100 000新圖號的行號;L10為1∶100 000新圖號的列號;X10為1∶100 000舊圖幅編號中相應比例尺地形圖的圖幅代碼值;[]表示小數取整;()表示小數取余。
例如:一幅1∶100 000基本比例尺圖號為G—48—90,由式(5)可知
則
取3位,不足3位前面用“0”補足,所以新圖號為G48D008006。
4.1 ∶50 000舊到新換算公式
1∶50 000是在1∶100 000的基礎上進行編號,1幅1∶100 000圖幅可分為4幅1∶50 000圖幅,按2行2列組成甲、乙、丙、丁編號,于是有下列公式
式中,H10為 1∶100 000新圖號的行號;L10為1∶100 000新圖號的列號;H5為1∶50 000新圖號的行號;L5為1∶50 000新圖號的列號;X5為1∶50 000舊圖幅編號中相應比例尺地形圖的圖幅代碼值,甲、乙、丙、丁的代碼值為1、2、3、4;[]表示小數取整;()表示小數取余。
例如:一幅1∶50 000基本比例尺圖號為G—48—90—乙,由式(6)可知
則
取3位,不足3位前面用“0”補足,所以新圖號為G48E015012。
5.1 ∶10 000舊到新換算公式
1∶10 000是在1∶100 000的基礎上進行編號,1幅1∶100 000圖幅可分為64幅1∶10 000圖幅,按8行8列組成(1)、(2)、…、(64)編號,于是有下列公式
式中,H10為 1∶100 000新圖號的行號;L10為1∶100 000新圖號的列號;H1為1∶10 000新圖號的行號;L1為 1∶10 000萬新圖號的列號;X1為1∶10 000舊圖幅編號中相應比例尺地形圖的圖幅代碼值;[]表示小數取整;()表示小數取余。
例如:一幅1∶10 000基本比例尺圖號為G—48—90—(8),由式(7)可知
則
取3位,不足3位前面用“0”補足,所以新圖號為G48G057048。
綜上所述,以上公式,通過多項式平移,可得到對應的新圖號到舊圖號的換算公式,計算方法同上,不再敘述。
通過對新舊分幅編號規(guī)則及換算關系的分析,利用VB、Excel VBA等編程語言,在VB中進行界面設計,在VB工程中引入Excelctl type library類型庫,編寫Excel VBA代碼,調入外部數據。在Command控件中對上述公式進行代碼編譯、調試,即可達到預期的結果,以下為計算a緯度、b經度數字碼部分程序代碼
因篇幅有限,其上述公式代碼在此省略。
本文所研究的轉換方法,無論是在第二次土地調查工作中,還是在高速公路、高速鐵路、水利水電站等項目的圖號資料快速查詢中,都得到了廣泛的應用和驗證。操作簡捷,界面人性化,具體操作時點擊“圖號輸入”按鈕,調入不同的外部數據,點擊“舊圖號→新圖號”或“新圖號→舊圖號”按鈕選擇換算方式,再點擊“圖號換算”即可完成運算結果,圖1為某縣第二次全國土地調查中舊圖幅號到新圖幅號的換算結果。
圖1
基本比例尺在地圖的編制方面經常使用,經過大量的工程實踐數據驗證,程序計算結果正確,運算速度快,可靠性良好,操作簡捷,界面人性化。在國家基本比例尺圖件縮編、土地資源調查、土地勘界、礦權核查等基礎資料分析工作中應用廣泛,可實現對不同基本比例尺新舊圖號的批量轉換和快速查找新舊圖號的目的,縮短了工程建設項目中基本圖件資料的查詢換算時間,能夠產生事半功倍的效果。
[1] 國家技術監(jiān)督局.GB/T 13989—1992國家基本比例尺地形圖分幅和編號[S].北京:中國標準出版社,1993.
[2] 劉復剛、張立人.關于我國國家基本比例尺地形圖分幅編號問題[J].齊齊哈爾師范學院學報:自然科學版,1994(2):49-51.
[3] 贠小蘇.第二次全國土地調查培訓教材[M].北京:中國農業(yè)出版社,2007(1):13-17.
[4] 楊莉.Visual Basic程序設計教程[M].北京:中國水利水電出版社,2002:339-343.
[5] 高春艷,劉彬彬,王斌.Visual Basic開發(fā)技術大全[M].人民郵電出版社,2007(1):134-166.
[6] [美]WALKENBACH J.Excel 2003高級VBA編程寶典[M].蓋江南,王勇,等譯.北京:電子工業(yè)出版社,2005:599-602.
A Tentative Analysis of the Mass Transformation between New and Old Sheet Numbers of Basic Scale
LENG Jiquan
0494-0911(2011)06-0051-03
P284
B
2010-07-23
冷繼全(1974—),男,貴州甕安人,工程師,主要從事數字化測量數據挖掘與處理研究。