【正文】 9。)==2) { m_blc=1:10 萬 。 21 m=atoi((5))。 x=fkh((m1),12)+1。 y=ykh((m1),12)+1。 (D%03d%03d,x,y)。 } else if(chrn(ch,39。(39。)==0amp。amp。chrn(ch,39。39。)==3) { m_blc=1:5 萬 。 str=(5)。 int w0。 w0=(39。39。)。 m=atoi((w0))。 n=ch[strlen(ch)1]64。 x=fkh((m1),12)*2+fkh((n1),2)+1。 y=ykh((m1),12)*2+ykh((n1),2)+1。 (E%03d%03d,x,y)。 } else if(chrn(ch,39。39。)==4amp。amp。chrn(ch,39。(39。)==0) { m_blc=1: 萬 。 str=(5)。 int w0。 w0=(39。39。)。 m=atoi((w0))。 p=atoi((1))。 n=ch[strlen(ch)3]64。 x=fkh((m1),12)*4+fkh((n1),2)*2+fkh((p1),2)+1。 y=ykh((m1),12)*4+ykh((n1),2)*2+ykh((p1),2)+1。 (F%03d%03d,x,y)。 } else if(chrn(ch,39。39。)==3) { m_blc=1:1 萬 。 str=(5)。 int w0。 w0=(39。39。)。 m=atoi((w0))。 int first,last。 first= (39。(39。)。 last= (39。)39。)。 n = atoi((first+1,lastfirst1))。 x=fkh((m1),12)*8+fkh((n1),8)+1。 22 y=ykh((m1),12)*8+ykh((n1),8)+1。 (G%03d%03d,x,y)。 } else { m_blc=1:5000。 str=(5)。 int w0。 w0=(39。39。)。 m=atoi((w0))。 int first,last。 first= (39。(39。)。 last= (39。)39。)。 n = atoi((first+1,lastfirst1))。 p=ch[strlen(ch)1]64。 x=fkh((m1),12)*16+fkh((n1),8)*2+fkh((p1),2)+1。 y=ykh((m1),12)*16+ykh((n1),8)*2+ykh((p1),2)+1。 (H%03d%03d,x,y)。 } m_tf2+=str。 } UpdateData(false)。 } 23 第四 章 使用說明 編寫目的 地圖的分幅編號，在地圖的生產(chǎn)、管理和使用方面都有重要意義。而地形圖圖幅號快速的查詢，以及新舊圖幅號之間的轉換，對于生產(chǎn)單位的企業(yè)、用戶來說就顯得尤為重要。 本手冊可以讓用戶快速的了解和實現(xiàn)圖幅號查詢系統(tǒng)的各種功能，能讓用戶了解在處理過程中遇到問題的解決方法。《使用說明》適用與從事地形圖管理的企業(yè)或者對地形圖圖幅號查詢有需求的所有用戶。 系統(tǒng) 功能 本軟件參考 GB/T 1398992 國家基本比例尺地形圖分幅和編號及舊地形圖分幅和編號規(guī)則編制。 主要功能如下： （ 1） 根據(jù)經(jīng)緯度或高斯坐標生成圖幅號，包括 1:100 萬、 1:50 萬、 1:25 萬、 1:10 萬、1:5 萬、 1: 萬、 1:1 萬、 1:5000 比例尺的圖幅號； （ 2） 由圖幅號查詢圖幅經(jīng)緯度范圍，圖幅理論面積等信息； （ 3）新舊圖幅號之間的相互轉換； （ 4）對每一模塊的查詢結果導出 txt 文件。 運行說明 啟動方法 雙擊打開圖幅號查詢系統(tǒng)文件夾， 3 個 DLL 文件為動態(tài)鏈接庫，雙擊文件名為“ tfhcx”的 exe 文件即可打開圖幅號查詢系統(tǒng)。 如（ 圖 41） 圖 41 24 查詢 （一）由經(jīng)緯度查詢新圖幅號 啟動程序后，在主界面上，單擊“由經(jīng)緯度查詢圖幅號”復選框，選中后，輸入目標地的經(jīng)度和緯度，格式： 90 度 7 分 34 秒寫成 ， 然后選擇 右邊的國家基本比例尺，點擊生成圖幅號，右下方的圖幅號編輯框會出現(xiàn)生成的新圖幅號。例如， 山西省太原市 市 中心位于北緯 37 度 52 分，東經(jīng) 112 度 33 分，比例尺選擇 1:1 萬，生成的圖幅號為 J49G052073。 如（ 圖 42） 圖 42 （二）由高斯坐標查詢圖幅號 在主界面上，單擊“由高斯坐標查詢圖幅號”復選框，選中后，選擇參考坐標系（如北京 54 坐標系，或西安 80坐標系等），輸入目標地的高斯投影 x,y 坐標 （ y坐標加過500Km） ，并輸入當?shù)刂醒胱游缇€經(jīng)度，選擇比例尺，點擊生成圖幅號。例如， 北京市天安門中心在 BJ54 坐標系中 x=， y=， 中央子午線經(jīng)度為 117 度，比例尺選擇 1:10 萬，生成的圖幅號為 J50D001005。 如（ 圖 43） 25 圖 43 （三）查詢圖幅范圍 點擊查詢圖幅范圍標簽，出現(xiàn)查詢圖幅范圍界面，在左邊輸入所要查詢的新圖幅號，圖幅面積處選擇一個參考系，點擊查詢圖幅范圍按鈕，右邊會出現(xiàn)圖幅范圍，圖幅面積等信息。如輸入上述圖幅號 J50D001005，結果如 （ 圖 44） 。 圖 44 （四）新舊圖幅號轉換 點擊新舊圖幅號轉換標簽，出現(xiàn)圖幅號轉換界面，在單選按鈕處選擇轉換方式，輸入原圖幅號，點擊轉換圖幅號按鈕，右邊的目標圖幅號編輯框會出現(xiàn)轉換后的圖幅號。如選擇轉換方式為 —— 舊圖號，輸入上述圖幅號 J50D001005，轉換的舊圖幅號為 J505。26 如（ 圖 45） 圖 45 清除數(shù)據(jù) 在每個功能頁面上查詢完一次結果或者輸入錯誤想重新查詢，則點擊每個頁面上的清除數(shù)據(jù)按鈕即可。 保存 本圖幅號查詢系統(tǒng)保存數(shù)據(jù)為 txt 文件，方法為點擊每個頁面上的保存按鈕，彈出的對話框選擇保存位置，輸入文件名稱，則會在相應位置生成一個 txt 文件。 （四）中圖幅號轉換 的保存結果如 （ 圖 46） 圖 46 27 退出 點擊每個頁面上的退出按鈕即可退出該系統(tǒng)，或者點擊界面右上角的“ X” 也可關閉該系統(tǒng)。 異常說明 由于某些點處在地形圖分幅的邊緣，可能在兩幅圖的公共邊上，也可能為四幅圖的公共點，此時查詢出的圖幅號就不止一個，會出現(xiàn)多個，如經(jīng)度 114 度，緯度 34 度，比例尺為 1:5 萬的圖幅號為 I49E012024， I50E012021， I49E013024， I50E013001，目標點在四幅圖的公共點上。 如（ 圖 47） 圖 47 28 參考文獻 [1] 王永成 ,楊保山 ,侯來福 . 淺談地形圖分幅與編號的計算方法 [J]. 西部探礦工程 . 2021(06) [2] 劉宏林 . 國家基本比例尺地形圖新舊圖幅編號變換公式及其應用 [J]. 測繪通報 . 1998(08) [3] 項仲貞 ,許承權 ,方子巖 . 中大比例尺地形圖標準分幅與編號查詢系統(tǒng)的實現(xiàn) [J]. 鐵道勘察 . 2021(01) [4] 張保綱北京市大比例尺地形圖圖幅號查詢信息系統(tǒng)的研制 [J] 解放軍測繪學院學報 19970320 29 致 謝 光陰似箭，轉眼間大學生活已接近尾聲。在大學的最后階段，畢業(yè)論文是我的重中之重，在論文完成的過程中，除了我自己的努力之外，也凝聚了很多人的心血。所以我要對幫助我完成論文的所有人表示感謝。 首先，我要對我的指導老師，表示誠摯的感謝。感謝我的老師，他以朋友的身份告訴我怎樣完成這篇論文，怎樣搜集資料。他在忙碌的工作中擠出時間來審查、修改我的論文，嚴格把關，循循善誘，在此我表示衷心感謝。 其次，我的同學也給予了熱情的幫助，感謝他們?yōu)槲姨岢鲇幸娴囊庖姾徒ㄗh，感謝他們對我無私的幫助。 最后，我要感謝我的父母，他們一直是我努力和學習的動力，感謝他們?yōu)槲宜冻龅囊磺小? 在寫論文的過程中，盡管我付出了很大的心血，但由于知識和經(jīng)驗的不足，不當之處懇請各位老師給予指正。 30 外文原文 （出自： Wei Z, Yang B S and Li Q Q. 2021. Automated extraction of building footprints from mobile LIDAR point clouds. Journal of Remote Sensing, 16(2): 286–296） Automated extraction of building footprints from mobile LIDAR point clouds WEI Zheng, YANG Bisheng, LI Qingquan Abstract This paper presents a novel method for automated extraction of building footprints from mobile LIDAR point rst generate the georeferenced feature image of mobile LIDAR point clouds using an interpolation method, and adopt imagesegmentation and contour extraction and tracing to extract building boundaries in the georeferenced feature image as the coarselevel of building footprints in Twodimensional imagery space. Then, the coarse level of building footprints is further rened byapplying planar segmentation on the extracted point clouds in the building boundaries. Finally, the triangulated irregular work (TIN) is used to achieve the ne level of building footprints. Dataset of residential areas captured by Optech’s LYNX mobile mapping system was tested to check the validities of the proposed method. Experimental results show that the proposed method provides a promising and valid solution for automatically 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。 Shibasaki, 2021) of Center for Spatial Information Science, The Uni