【正文】 ly extracting building footprints from mobile LIDAR point clouds. Key words: mobile LIDAR, footprints extraction, image segmentation, point cloud segmentation, triangulated irregular work 1 INTRODUCTION Mobile Laser Scanning (MLS) provides an efcient solution forcapturing surface data of roads, buildings and trees fast, and dynamic. It has been widely used in many elds, such as digital city,basic surveying and mapping, urban planning, transportation, andenvironment protection. Some MLS systems have been developed,such as VLMS (Zhao amp。在大學(xué)的最后階段，畢業(yè)論文是我的重中之重，在論文完成的過程中，除了我自己的努力之外，也凝聚了很多人的心血。如輸入上述圖幅號(hào) J50D001005，結(jié)果如 （ 圖 44） 。 系統(tǒng) 功能 本軟件參考 GB/T 1398992 國家基本比例尺地形圖分幅和編號(hào)及舊地形圖分幅和編號(hào)規(guī)則編制。 y=ykh((m1),12)*16+ykh((n1),8)*2+ykh((p1),2)+1。(39。 str=(5)。 last= (39。 w0=(39。 x=fkh((m1),12)*4+fkh((n1),2)*2+fkh((p1),2)+1。 str=(5)。 (E%03d%03d,x,y)。 int w0。 } else if(chrn(ch,39。 m_tf2=國家基本比例尺新圖幅號(hào)不包括 1:20 萬 。)。[39。||ch[5]==39。 int m,x,y,n,p。 (%c,p+64)。: m_blc=1:5000。 (%d,m)。 p=ykh((x1),2)*2+ykh((y1),2)+1。F39。 m=fkh((x1),2)*12+fkh((y1),2)+1。: m_blc=1:10 萬 。C39。 switch(ch[3]) { case 39。 } 新舊圖幅號(hào)轉(zhuǎn)換代碼： void prop3::OnButton1() { UpdateData(true)。 while ((BpB)jtoh(0,0,)||(lpl)jtoh(0,0,)) 18 { B=Bp。 return h7。 return h7。 } double GAOSI::funx (double h2,double h3) //計(jì)算 Fx(B,l) { double h4,h5,h6,q2,q4,q6,h7。 a2=*N*sinh*cosh。 } void GAOSI::an (double m) //計(jì)算 a1a6 { double n,v,N,sinh,cosh,t。 } double GAOSI::htoj (double m0) //弧度轉(zhuǎn)角度 { double u。 b2=b01。 系統(tǒng)設(shè)計(jì)路線 圖 31 系統(tǒng)界面 本次圖幅號(hào)查詢系統(tǒng)是用 VC++語言設(shè)計(jì)實(shí)現(xiàn)的， 在主對(duì)話框中添加標(biāo)簽控件，設(shè)立四個(gè)標(biāo)簽：查詢圖幅號(hào)、查詢圖幅范圍、新舊圖幅號(hào)轉(zhuǎn)換、關(guān)于。 [ ]表示商取整； ( )表示商取余； （ 6） 1: 萬地形圖從新圖幅號(hào)換算舊圖幅號(hào) 舊圖幅號(hào)中的第 3 組代碼 —— 1: 萬地形圖所在的 1:10 萬地形圖在 1:100 萬地形圖中的位置代碼“ m”： m=[???14 ]*12+[???14 ]+1； 公式 211 舊圖幅號(hào)中的第 4 組代碼 —— 1: 萬地形圖所在的 1:5 萬地形圖在 1:10 萬地形圖的位置代碼“ n”： n=([??+12 ]?12 )*2+([??+12 ]?12 )+1； 公式 212 舊圖幅號(hào)中的第 5 組代碼 —— 1: 萬地形圖在 1:5 萬地形圖的位置代碼“ p”： p=(???12 )*2+(???12 )+1； 公式 213 計(jì)算結(jié)果， n=1 時(shí)取 A， n=2 時(shí)取 B， n=3 時(shí)取 C， n=4 時(shí)取 D； x 表示 1: 萬地形圖新圖幅號(hào)中的圖幅行號(hào)， y 表示 1: 萬地形圖新圖幅號(hào)中的8 圖幅列 號(hào)。 )/△ L]； 公式 22 d=[(B/6176。由經(jīng)線和緯線所圍成的每一個(gè)小方格為一幅 1： 1 00 萬 地形圖，它們的編號(hào)由該圖所在的行號(hào)和列號(hào)組合而成，例如北京所在的 1： 1 00 萬 地形圖圖號(hào)為 J 50。39。緯差是 539。構(gòu)成以 ?、144 為代號(hào)的 1:10 萬地形圖，如： J5092； (3)1:5 萬、 1: 萬、 1:1 萬 、 1:5000 地形圖分幅編號(hào) 將一幅 1:10 萬地形圖分為 4 幅，構(gòu)成以 A、 B、 C、 D 為代號(hào)的 1:5 萬地形圖， 圖幅經(jīng)差為 15’，緯差為 10’， 如： J5092A； 將一幅 1:5 萬地形圖分為 4 幅，構(gòu)成以 4 為代號(hào)的 1: 萬地形圖， 圖幅經(jīng)差為 7’30’’，緯差為 5’， 如： J5092A1； 將一幅 1:10 萬地形圖分為 8 行、 8 列，共 64 幅 1:1 萬地形圖，分別以 (1)、 (2)、?(64)表示，經(jīng)差為 339。 將一幅 1:100 萬地形圖分為 4 幅 ,圖幅經(jīng)差為 3176。劃作 22 個(gè)橫列，依次用 A、 B、??V 表示；從經(jīng)度 180176。查詢地形圖庫依據(jù)的元素很多，如圖名、檔案號(hào)、圖幅號(hào)、地區(qū)等，但最基本、最便捷的查詢?cè)厥亲鴺?biāo)和經(jīng)緯度，只要知道某點(diǎn)的坐標(biāo)或經(jīng)緯度，系統(tǒng)便可自動(dòng)搜尋該點(diǎn)的圖幅號(hào)，并進(jìn)行生成 ?？傊?，圖幅號(hào)為了地形圖的管理、保存而出現(xiàn)。 地形圖 分幅有兩種方法：一 種是矩形分幅，另一種 是經(jīng)緯線分幅。I 國家基本比例尺地形圖圖幅號(hào)查詢系統(tǒng)的設(shè)計(jì) 摘 要 圖幅 號(hào) 的快速查詢對(duì)于企業(yè)和用戶 有著廣泛的用途 。常用的地圖編號(hào)有行列式、自然序數(shù)式、行列－自然序數(shù)式和西南角圖廓點(diǎn)坐標(biāo)公里數(shù)編號(hào)等。隨著計(jì)算機(jī)技術(shù)的應(yīng)用和發(fā)展，許多企事業(yè)單位在生產(chǎn)管理中會(huì)用到大量的地形圖，這些地形圖以圖號(hào)和圖名來標(biāo)識(shí)，人們通過查找索引表來找到所需要的圖紙。 2 圖幅號(hào)查詢系統(tǒng)設(shè)計(jì)意義 國家基礎(chǔ)地形圖的使用單位大多已建立了數(shù)據(jù)庫管理查詢系統(tǒng)。止，按緯差每 4176。 （ 2） 1:50 萬、 1:25 萬、 1:20 萬、 1:10 萬地形圖分幅編 號(hào) 1:50 萬、 1:25 萬、 1:20 萬、 1:10 萬地形圖均以 1:100 萬圖為基礎(chǔ)分幅、編號(hào)。緯差為 2039。3039。39。從 180 度經(jīng)線算起自西向東每經(jīng)差 6 度為一列，全球分為 60 列，一次用阿拉伯?dāng)?shù)字 ?? 60 表示其相應(yīng)列號(hào)。 /△ L[(L/4176。 [ ]表示商取整； ( )表示商取余； （ 5） 1:5 萬地形圖從新圖幅號(hào)換算舊圖幅號(hào) 舊圖幅號(hào)中的第 3 組代碼 —— 1:5 萬地形圖所在的 1:10 萬地形圖在 1:100 萬地形圖中的位置代碼“ m”： m=[???12 ]*12+[???12 ]+1； 公式 29 舊圖幅號(hào)中的第 4 組代碼 —— 1:5 萬地形圖在 1:10 萬地形圖的位置代碼“ n”： n=(???12 )*2+(???12 )+1； 公式 210 計(jì)算結(jié)果， n=1 時(shí)取 A， n=2 時(shí)取 B， n=3 時(shí)取 C， n=4 時(shí)取 D； x 表示 1:5 萬地形圖新圖幅號(hào)中的圖幅行號(hào)， y 表示 1:5 萬地形圖新圖幅號(hào)中的圖幅列 號(hào)。簡潔要求系統(tǒng)界面友好，用戶輸入數(shù)據(jù)方便，系統(tǒng)代碼簡潔明了，效率高；準(zhǔn)確指 數(shù)據(jù)精度高，數(shù)據(jù)類型采用雙精度，并且數(shù)據(jù)的格式符合國家規(guī)范中的要求，查詢結(jié)果準(zhǔn)確無誤。 b0=*ep*ep+*mi(ep,4)*mi(ep,6)+*mi(ep,8)。 return q。 } return q1。 a1=N*cosh。 return h1。 h7=q2*h3*h3+q4*mi(h3,4)+q6*mi(h3,6)。 h7=q3*mi(h3,3)+q5*mi(h3,5)。 lp=(m2funy(Bp,l))/a1p。 fsjd=htoj(lp+m3)。 y=atoi((3))。 case 39。D39。: m_blc=1:5 萬 。 case 39。 m_tf2=m_tf2++str。 m=fkh((x1),8)*12+fkh((y1),8)+1。H39。 p=ykh((x1),2)*2+ykh((y1),2)+1。 ch=(LPSTR)(LPCTSTR)m_tf1。C39。 } else if (ch[5]==39。]39。) { m_blc=1:20 萬 。 (D%03d%03d,x,y)。 str=(5)。 y=ykh((m1),12)*2+ykh((n1),2)+1。)==0) { m_blc=1: 萬 。 n=ch[strlen(ch)3]64。 int w0。)。 } else { m_blc=1:5000。 first= (39。 x=fkh((m1),12)*16+fkh((n1),8)*2+fkh((p1),2)+1?！妒褂谜f明》適用與從事地形圖管理的企業(yè)或者對(duì)地形圖圖幅號(hào)查詢有需求的所有用戶。 如（ 圖 43） 25 圖 43 （三）查詢圖幅范圍 點(diǎn)擊查詢圖幅范圍標(biāo)簽，出現(xiàn)查詢圖幅范圍界面，在左邊輸入所要查詢的新圖幅號(hào)，圖幅面積處選擇一個(gè)參考系，點(diǎn)擊查詢圖幅范圍按鈕，右邊會(huì)出現(xiàn)圖幅范圍，圖幅面積等信息。 如（ 圖 47） 圖 47 28 參考文獻(xiàn) [1] 王永成 ,楊保山 ,侯來福 . 淺談地形圖分幅與編號(hào)的計(jì)算方法 [J]. 西部探礦工程 . 2021(06) [2] 劉宏林 . 國家基本比例尺地形圖新舊圖幅編號(hào)變換公式及其應(yīng)用 [J]. 測繪通報(bào) . 1998(08) [3] 項(xiàng)仲貞 ,許承權(quán) ,方子巖 . 中大比例尺地形圖標(biāo)準(zhǔn)分幅與編號(hào)查詢系統(tǒng)的實(shí)現(xiàn) [J]. 鐵道勘察 . 2021(01) [4] 張保綱北京市大比例尺地形圖圖幅號(hào)查詢信息系統(tǒng)的研制 [J] 解放軍測繪學(xué)院學(xué)報(bào) 19970320 29 致 謝 光陰似箭，轉(zhuǎn)眼間大學(xué)生活已接近尾聲。 在寫論文的過程中，盡管我付出了很大的心血，但由于知識(shí)和經(jīng)驗(yàn)的不足，不當(dāng)之處懇請(qǐng)各位老師給予指正。感謝我的老師，他以朋友的身份告訴我怎樣完成這篇論文，怎樣搜集資料。26 如（ 圖 45） 圖 45 清除數(shù)據(jù) 在每個(gè)功能頁面上查詢完一次結(jié)果或者輸入錯(cuò)誤想重新查詢，則點(diǎn)擊每個(gè)頁面上的清除數(shù)據(jù)按鈕即可。 如（ 圖 41） 圖 41 24 查詢 （一）由經(jīng)緯度查詢新圖幅號(hào) 啟動(dòng)程序后，在主界面上，單擊“由經(jīng)緯度查詢圖幅號(hào)”復(fù)選框，選中后，輸入目標(biāo)地的經(jīng)度和緯度，格式： 90 度 7 分 34 秒寫成 ， 然后選擇 右邊的國家基本比例尺，點(diǎn)擊生成圖幅號(hào)，右下方的圖幅號(hào)編輯框會(huì)出現(xiàn)生成的新圖幅號(hào)。 } UpdateData(false)。)39。39。 n = atoi((first+1,lastfirst1))。 m=atoi((w0))。 } else if(chrn(ch,39。39。)==4amp。)。amp。)==2) { m_blc=1:10 萬 。 y=ykh((m1),4)+1。 first= (39。 m=ch[5]64。A39。 case 0: m_blc=1:100 萬 。 m_tf2=m_tf2++str。 (%d,n)。 break。 (%d,m)。 n=ykh((x1),2)*2+ykh((y1),2)+1。 19 m_tf2=m_tf2++str。 (%d,m)。