簡介面向三維點云處理與逆向工程學習者這份35.02MB的開源軟件包提供CloudCompare的完整源碼與工程文件適用于點云可視化、法向量計算與優(yōu)化、泊松構(gòu)網(wǎng)、濾波等典型任務(wù)。包內(nèi)共2000個文件以C/C源碼590個h、397個cpp、34個c為骨架配套118個ui界面文件、32個qrc資源文件以及png/jpg等圖像資源便于二次開發(fā)與界面定制同時包含lyx/tex/txt等文檔可輔助理解算法與編譯流程。已有4783人學習下載說明其在相關(guān)領(lǐng)域有較廣泛的參考價值。從內(nèi)容構(gòu)成看包體兼顧源代碼、工程配置cmake/pro/sln、腳本bat/sh與說明文檔目錄結(jié)構(gòu)清晰既適合研究者分析點云處理實現(xiàn)思路也適合開發(fā)者基于源碼進行功能擴展是一份可直接用于學習與開發(fā)的三維點云處理工具資源。 第一次打開三維點云PointCloud處理軟件CloudCompare.zip這個壓縮包的時候很多人會覺得它不像一個“軟件”解壓出來一堆文件夾沒有安裝向?qū)щp擊exe就彈出一個灰黑色界面左邊一棵樹右邊一片空白。我當年也是從這一步開始一步步把點云從“看一眼”玩到“能處理”。這篇文章想寫給所有剛接觸點云數(shù)據(jù)、或者已經(jīng)在用別的軟件但想切到CloudCompare的朋友。它不只是介紹按鈕在哪而是把我在實際項目里踩過的坑、反復調(diào)參的經(jīng)驗、以及“為什么這樣做”的邏輯一次性講清楚。無論你手里是無人機航測的LAS點云還是地面掃描的PLY數(shù)據(jù)或者只是從網(wǎng)上下了一個練習用的PointCloud壓縮包思路都適用。1. 從zip到可用這個軟件真正干活的部分在哪1.1 解壓即用的綠版別被目錄結(jié)構(gòu)嚇到CloudCompare的Windows版本通常是一個zip壓縮包解壓后能看到主程序、一堆dll文件、plugins文件夾、translations文件夾等。很多新手習慣性去找安裝包其實這個生態(tài)里壓縮包就是安裝方式解壓到固定目錄雙擊CloudCompare.exe就能用。這里有兩個細節(jié)我建議你第一次就做對目錄路徑不要帶中文和空格盡量放在類似D:\Tools\CloudCompare的位置。雖然新版對中文路徑的兼容性好了不少但在處理某些插件和第三方庫時路徑里一旦出現(xiàn)中文偶爾會出一些奇怪的加載失敗問題。不要單獨把exe拖到桌面運行它運行時需要讀取同目錄下的dll和插件單獨運行會提示缺失組件。解壓后建議先看一眼plugins目錄里面那些以q開頭的文件就是功能插件。CloudCompare的主程序其實承擔的是渲染、編輯、測量等基礎(chǔ)能力很多專業(yè)功能比如CSF地面濾波、M3C2變化檢測是靠插件完成的。如果你從別處下載了額外插件放進來后要在菜單欄的“插件”里看到通常需要重啟程序。1.2 正式處理數(shù)據(jù)前先把這幾個全局選項調(diào)好打開軟件后大多數(shù)人急著拖數(shù)據(jù)進去我建議先花兩分鐘進Tools Options工具 選項看一眼下面的設(shè)置會直接影響后續(xù)操作體驗。界面語言新版本菜單里可以切換中文如果你對英文界面不敏感保留英文反而更容易查資料因為社區(qū)里的教程、論壇討論基本都是英文關(guān)鍵詞。顯示設(shè)置把“顯示當前實體數(shù)量”等選項打開方便隨時確認點云包含了多少個點。OpenGL設(shè)置如果電腦顯卡驅(qū)動正常建議讓軟件走獨立顯卡。CloudCompare對OpenGL的要求不算低千萬級別點云旋轉(zhuǎn)縮放時集顯和獨顯的差距會非常明顯。我的習慣是第一次打開后先隨便生成一個正方體或球體File Primitive factory確認旋轉(zhuǎn)、平移、縮放、點選都正常了再加載真實數(shù)據(jù)。這個習慣幫我排除過不少環(huán)境問題。2. 點云加載不是拖進去就完事格式、坐標系和單位要先對齊2.1 支持的格式和區(qū)分“點云”與“網(wǎng)格”的差異CloudCompare對點云格式的寬容度在同類軟件里屬于很不錯的常見的LAS、LAZ、XYZ、ASC、TXT、PLY、OBJ、PCD、E57等都能讀。但首次拖入文件時建議注意一個容易忽略的地方彈窗會詢問你“是否將此文件作為點云加載還是網(wǎng)格加載”尤其對于PLY、OBJ這種同時支持點云和網(wǎng)格的格式如果選錯后續(xù)很多點云專用工具會變成灰色不可用。如果你手里的數(shù)據(jù)是XYZ或TXT這種純文本格式導入時大概率會彈出一個ASCII導入對話框需要手動指定每一列的含義哪幾列是XYZ坐標哪幾列是顏色哪幾列是強度/回波號。這一步務(wù)必仔細核對選錯一列的話點云可能整體錯位或者把高度數(shù)值當成顏色來渲染圖像會非常詭異。2.2 單位、比例和坐標偏移為什么導入后什么都看不見這是新手遇到的問題里發(fā)生率最高的一個沒有之一。做過測繪、掃描項目的人應該都有體會不同設(shè)備導出的點云單位可能是毫米、厘米或米坐標系可能是地理坐標經(jīng)緯度或高斯平面坐標或局部坐標。如果你的點云是地理坐標X和Y值往往是幾十萬到幾百萬的數(shù)字而Z值可能是幾十到幾百這時候你把點云導入場景它會懸在地圖邊界的老遠之外在默認視角下什么都看不到或者只看到一團不知道是什么的色塊。解決辦法很簡單選中DB樹里的實體在“實體屬性”面板里查看其重心坐標和尺寸。如果發(fā)現(xiàn)坐標范圍極大先手動拖動視角朝著坐標原點方向縮小找更快的辦法是用Edit Multiply / Scale編輯 縮放將點云整體縮放到合理范圍或者在加載前用文本編輯器對源數(shù)據(jù)進行坐標平移把重心平移到原點附近。我做的很多項目數(shù)據(jù)本身就是UWB定位或激光掃描的局部坐標不會有這個問題但每當接手別人的數(shù)據(jù)時第一步永遠都是“先看范圍再談處理”。2.3 加載后先做的三個檢查動作數(shù)據(jù)正常顯示之后先別急著干活我固定會做三件事在DB樹里單擊點云看總點數(shù)是否和源數(shù)據(jù)一致。如果數(shù)量差很多多半是導入時抽樣了或格式轉(zhuǎn)換丟了數(shù)據(jù)。按Alt鼠標左鍵或使用工具欄的“點選”工具點擊幾個點讀取一下坐標值確認坐標體系和預期一致。給點云設(shè)置一個臨時顏色比如用高度字段在z方向做一個連續(xù)漸變著色快速掃一眼地形走勢和噪聲點。這一步能讓你在正式操作前就對數(shù)據(jù)質(zhì)量有個直觀判斷。這些檢查最多花兩分鐘卻能在后續(xù)處理里幫你省掉大量排查時間。3. ICP配準為什么我按教程配了十次還是歪的3.1 配準前必須完成的粗配準才是決定成敗的一步“cloudcompare中ICP配準怎么用”是搜索熱度最高的問題。很多教程上來就說選中兩個點云Tools Registration Fine registration工具 配準 精細配準然后填參數(shù)Apply然就對齊了。但真實情況是如果你兩片點云的初始位置差得遠ICP幾乎一定會收斂到局部最優(yōu)解結(jié)果就是兩塊點云歪七扭八地“拼”在一起。我自己的經(jīng)驗是ICP只是精配準它只能修正小誤差不能解決大錯位。在跑ICP之前至少先用Tools Registration Align (point pairs picking)點對配準手動選三組以上同名點比如墻角、標志物的邊緣等建立粗略坐標變換關(guān)系。這個粗配準操作不用追求特別精確只要讓兩片點云大致重合ICP才有機會找到正確的最優(yōu)解。手動選點時要注意同名點分布要盡量分散在空間的三個方向上不要都在一條線上否則算出來的旋轉(zhuǎn)矩陣容易退化。選點完成后軟件會給出一個RMS誤差值如果粗配準的RMS還很大說明點對選得不準先調(diào)整點對再進ICP。3.2 ICP參數(shù)里最值得關(guān)注的幾個選項在CloudCompare的精細配準對話框中你會看到幾個參數(shù)很多人在這一步直接踩坑。“重疊度”approximate overlap這個參數(shù)表示兩片點云重疊區(qū)域所占的比例默認是100%。如果你配準的是一片掃描的局部區(qū)域與整體模型重疊度可能只有60%-80%這時候仍用100%去算算法會強行把所有點都當成“可對應點”結(jié)果自然偏差很大。我的經(jīng)驗是根據(jù)實際掃描范圍先估一個重疊度不放心就從80%試起看結(jié)果不對再降?！白畲蟮螖?shù)”max iterations默認值通常夠用但如果你發(fā)現(xiàn)結(jié)果還沒收斂就停了可以增大這個值。要注意的是迭代次數(shù)不是越大越好配準本質(zhì)是逐步微調(diào)收斂后繼續(xù)迭代并不會讓結(jié)果更好。“最終誤差閾值”final overlap / threshold這個參數(shù)用于判斷是否停止可以理解為“兩片點云對應點之間距離小于多少就認為配準完成”。它的取值與數(shù)據(jù)單位強相關(guān)。如果你的點云單位是米閾值設(shè)置成0.01就是1厘米精度如果單位是毫米閾值設(shè)置成0.01就是0.01毫米幾乎不可能達到。很多人配不準就是因為單位搞錯了閾值設(shè)得完全不合邏輯。另外如果兩片點云之間尺度和分辨率差異很大可以在對話框里打開“調(diào)整尺度”adjust scale選項。大部分場景不需要打開但處理不同設(shè)備來源的數(shù)據(jù)時打開它往往會有驚喜。3.3 怎么判斷配準結(jié)果到底好不好配準結(jié)束后CloudCompare會在控制臺輸出最終的RMS均方根誤差值。這個數(shù)字是衡量配準質(zhì)量的最直接指標但它不是絕對真理RMS小不代表一定沒有局部扭曲RMS大也不代表完全失敗。我更推薦的做法是配準后把兩層點云疊加顯示開啟剖切面或透明度渲染在不同視角下檢查關(guān)鍵特征邊緣是否重合。還有一個小技巧配準完成后立即用Edit Colors Colorize編輯 顏色 著色給兩片點云分別賦予兩種對比鮮明的顏色比如紅色和青色然后疊加顯示。顏色重疊區(qū)域呈現(xiàn)混合色錯位區(qū)域會出現(xiàn)明顯的顏色分界線任何微小偏移都會被直觀放大。這個技巧在我做變形監(jiān)測和測繪成果比對時幾乎每次都用。4. 分割、測量和渲染高頻操作的最佳路徑4.1 用連通組件和套索分割把“想要的”和“不想要的”分開點云處理里最日常的操作就是分割把地面點、建筑、樹木、噪聲一個個分開。CloudCompare里最高頻的兩種分割方式是手動套索分割先選中點云點擊Edit Segment編輯 分割然后在三維視圖里用鼠標畫閉合區(qū)域雙擊確認后新實體出現(xiàn)在DB樹中。這個操作適合小范圍、快速切分缺點是依賴人眼判斷。連通組件分割Connected Components菜單在Tools Segmentation Connected Components工具 分割 連通組件它基于八叉樹把空間上不連續(xù)的點云拆分成多個獨立實體。適合大范圍地物分離比如把一片點云里所有獨立建筑自動拆開。你需要設(shè)置的參數(shù)主要是八叉樹層級Octree level和最少點數(shù)。層級越低分割粒度越粗最少點數(shù)則用來濾除孤立的噪聲小簇。實際操作中我常常把兩種方式組合起來用先用連通組件粗拆再用套索精修效率比單純手動框選高得多。4.2 距離、面積和體積測量別在高斯圖上翻車CloudCompare的測量工具集中在Tools Measure工具 測量下。最常用的是點距離測量和點到點距離。還有一個好用但容易被忽略的功能是“2.5D體積計算”Tools Volume 2.5D Volume它可以基于兩個點云或網(wǎng)格計算堆體體積常用于土方、料堆、采石場測量。具體操作是選擇基準面實體和被測點云設(shè)置步長grid step軟件會把空間網(wǎng)格化后計算頂面和底面的高度差再積分出體積。這里有個非常關(guān)鍵的教訓步長的設(shè)置不能比點云平均間距小太多。步長設(shè)得過小會導致大量空白格被插值填充計算出的體積失真步長設(shè)得過大地形細節(jié)又被抹平。我一般先通過點云屬性里的“平均點間距”cloud statistics拿到一個參考值再以它乘以2~5倍作為步長結(jié)果比較穩(wěn)定。4.3 用標量場和顏色梯度把隱患“看”出來點云本身通常帶強度、高程、回波號等屬性在CloudCompare里這些屬性被稱為“標量場”Scalar Field。把標量場映射到顏色上是點云數(shù)據(jù)分析里最直觀的一招。操作方法不復雜在DB樹里選中點云展開它的標量字段列表右鍵想要顯示的字段選擇“設(shè)置為活動標量場”然后點頂部工具欄的顏色漸變圖標點云就會按照該字段的變化進行著色。配合Tools Display Color Scale調(diào)整色帶范圍可以把你要關(guān)注的區(qū)間壓縮到某個顏色帶里比如把高程范圍限制到設(shè)計值附近偏差超過閾值的點就變成紅色一眼就能鎖定異常區(qū)域。這個用法在質(zhì)量檢查中替代了人工逐個點看坐標的笨辦法效率和準確率都高了一個量級。5. 大批量點云的效率問題從拖拽操作走向命令行5.1 點云太卡先別急著升級電腦試試降采樣和八叉樹見過不少朋友一加載海量點云就卡到屏幕轉(zhuǎn)圈直接判斷“電腦不行”。其實很多時候不是電腦不行而是數(shù)據(jù)沒有被合理組織。CloudCompare對點云顯示做了一種優(yōu)化數(shù)據(jù)加載后會構(gòu)建八叉樹對同一個區(qū)域的不同層級做多分辨率顯示旋轉(zhuǎn)縮放時軟件只渲染當前視角下的合適點數(shù)。這個功能叫“八叉樹LOD”如果發(fā)現(xiàn)點云特別卡第一步檢查是否已經(jīng)啟用八叉樹顯示。如果你只是想觀察、測量并不需要千萬級別的密度可以直接用Edit Subsample編輯 降采樣把點云按固定百分比或固定間距抽稀。做一次降采樣保留10%甚至1%的點在很多預覽和方案匯報場景下視覺差異完全能接受但流暢度提升是肉眼可見的。5.2 用命令行批量轉(zhuǎn)換和自動保存省掉重復勞動當你有幾十個文件要做格式轉(zhuǎn)換或統(tǒng)一抽稀時在圖形界面里一個個操作會崩潰。CloudCompare自帶命令行模式可以在不打開窗口的情況下批量處理。我經(jīng)常用來批量把LAZ轉(zhuǎn)成LAS或者批量降采樣腳本大致長這樣CloudCompare.exe -SILENT -O input1.laz -O input2.laz -C_EXPORT_FMT LAS -AUTO_SAVE true參數(shù)含義不復雜-O指定輸入文件-C_EXPORT_FMT指定導出格式-AUTO_SAVE true表示處理完直接保存結(jié)果。把這類指令寫成一個循環(huán)腳本放到服務(wù)器或工作站上掛著跑就能讓軟件自己把幾十個文件處理完。命令行模式的限制是它不會自動“智能地”完成你所有鼠標操作但對重復性高、規(guī)則明確的任務(wù)非常合適。5.3 插件生態(tài)里值得裝的幾個CloudCompare真正強大之處是插件生態(tài)。我日常比較依賴的有這三個CSF布料模擬濾波專門做地面點與非地面點分離比單純按高度閾值切割智能得多。它的原理是模擬一塊布從上方掉落到翻轉(zhuǎn)后的地形上貼合地形的那部分就是地面點。M3C2用于兩期點云之間的變化檢測能計算兩個點云間的距離并給出置信區(qū)間做邊坡監(jiān)測、采空區(qū)變化分析很合適。RANSAC Shape Detection從點云中自動提取平面、球體、圓柱、圓環(huán)等幾何形狀對工業(yè)構(gòu)件、室內(nèi)結(jié)構(gòu)的建模很有用。這些插件大多默認不啟用。需要到插件列表里勾選或者在啟動時按住某個快捷鍵選擇加載不同版本方式略有差異如果你打開菜單發(fā)現(xiàn)沒有去插件設(shè)置里看一下即可。6. 排錯記錄加載、內(nèi)存和配準的典型問題排查6.1 文件加載閃退、亂碼和報錯的排查順序如果你的CloudCompare在加載文件時閃退或者中文文件名變成了亂碼不要著急重裝軟件先按這個順序排查確認文件路徑和文件名是否為純英文有時甚至文件所在目錄層級很深也會觸發(fā)問題可以先復制到簡單路徑再試。查看插件是否加載正常LAS/LAZ格式依賴對應插件如果插件缺失文件會被識別失敗。如果是大型文件一開始就加載失敗或閃退建議先用別的小工具比如命令行模式轉(zhuǎn)換格式或用文本工具查看文件頭確認文件不是損壞的空殼。步驟可以自己復現(xiàn)同一份數(shù)據(jù)放到不同路徑下試一次就能很快定位到問題是否出在路徑和命名上。實際項目里因為這個原因浪費時間的案例非常多。6.2 bad allocation 不一定是你內(nèi)存不夠也可能是這一步?jīng)]做“bad allocation”或“Failed to allocate memory”是高頻報錯。很多人的第一反應是內(nèi)存條不夠?qū)嶋H上還有一個非常隱蔽的誘因你試圖對原本存儲格式不連續(xù)的點云執(zhí)行需要隨機訪問的操作導致內(nèi)存中突然需要加載大量數(shù)據(jù)副本。我的處理思路是先把大點云降采樣或按區(qū)塊切分比如用連通組件把大面積點云拆開再重新嘗試原來昂貴的操作。另外確認軟件是否運行的64位版本也很重要32位程序在Windows下默認只能使用約2GB內(nèi)存超大點云必掛無疑。6.3 ICP結(jié)果不收斂時的日志與應對跑完ICP之后如果發(fā)現(xiàn)RMS一直不下降或者兩個點云越配越開先從三個方向排查兩片點云的點數(shù)懸殊某些區(qū)域一個點對應多個點導致約束不穩(wěn)定。先做降采樣讓兩點云密度接近再配準。初始位置偏差太大。退回粗配準增加手動同名點數(shù)量或改善點對分布。閾值和單位不匹配。檢查數(shù)據(jù)單位把迭代停止條件和單位對齊。我在做這類排錯時會在每次配準前把參數(shù)截圖保存配準完對比日志這樣就能知道是哪一步改動讓結(jié)果變好或變差后續(xù)也能沉淀成自己的經(jīng)驗表格。另外還有一個平時不容易想起來的小技巧如果配準完發(fā)現(xiàn)兩塊點云在Z方向上有固定偏差可以在配準對話框里勾選“固定某個軸的旋轉(zhuǎn)/平移”具體選項名稱隨版本略有不同強制算法只在一個方向上修正往往能比全自由度的配準更穩(wěn)定。玩點云這幾年我的體會是CloudCompare這類工具的上手曲線其實不在“操作”而在“理解數(shù)據(jù)”單位是什么、坐標系是什么、點的密度和噪聲有多大、處理目標到底是什么。把這些想清楚了界面上那些按鈕自然就能各就各位。希望這篇圍繞解壓、加載、配準、分割到批量處理的經(jīng)驗整理能幫你少走我當年走過的彎路。本文還有配套的精品資源點擊獲取