基于Cesium的BIM與實景三維模型集成可視化研究

方孟元羅年學許 毅祁平利

1武漢大學測繪學院，湖北 武漢，430079

2中鐵第五勘察設計院集團有限公司，北京，102600

建 筑 信 息 模 型（building information model，BIM）概念起源于20世紀70年代，是對建筑的物理和功能特性的數字化表達，貫穿于建筑從概念、設計、施工、運營維護及拆除的整個生命周期中，提供信息支撐和決策依據。隨著信息技術的高速發展，BIM成為了規劃設計領域的研究熱點，在水立方、上海中心大廈等大型項目中得到廣泛應用［1，2］。

在建筑設計、施工等各個環節都離不開周邊環境信息做支撐，GIS作為地理信息的數字化表達載體，正為BIM提供了更大尺度的信息支持。GIS的各類空間分析功能可以為BIM提供決策參考，真正將建筑與環境融為一體［3-6］。吳志強等［7］進一步提出了城市信息模型（city information modeling，CIM）的概念，即CIM=BIM+GIS+物聯網（internet of things，IoT），將微觀的BIM、宏觀的GIS與IoT數據進行集成統一，形成綜合數據處理計算平臺，為智慧城市建設提供了技術支撐。

因此，BIM與GIS的組合越來越多的受到國內外學術界和工程界的廣泛關注。文獻［8-13］分別從語義關聯、數據融合、集成方式和工程應用等方面對BIM與GIS的融合展開了探討。實景三維模型是借助傾斜攝影測量技術獲取的地表三維模型，通過多視角航空攝影、空三加密、密集匹配、紋理映射等環節自動計算生成，可真實反映地物地貌特征、數據精度高、自動化程度高，是目前主流的三維GIS數據源之一［4］。Cesium是一款開源的Web端三維GIS開發包，基于WebGL技術和3D-Tiles數據標準實現海量二三維空間數據的跨平臺繪制展示。

本文以某磁浮鐵路項目為例，介紹BIM與GIS數據集成可視化的實現思路及基于Cesium的實現方法；進一步地介紹BIM與實景三維模型的動態卷簾功能實現思路及效果。

1 集成可視化實現思路

BIM與傾斜攝影模型的集成可視化流程如圖1所示。

圖1 BIM與傾斜攝影模型集成可視化流程Fig.1 Integrated Visualization Process of BIM and Oblique Photography Model

由于各大機構和軟件廠商分別定義了各自的數據標準，不同數據標準對應了不同的數據文件格式、要素組織方法及空間索引結構：BIM標準包括RVT、DGN、IFC標準等；實景三維模型標準包括OSGB、SLPK、S3M標準等。與此同時，需分別對BIM模型和實景三維模型進行分層、分塊處理，以保證Web端加載性能。3D-Tiles結構是Cesium原生支持的數據標準，渲染時可根據視錐體范圍和瓦片索引進行動態調度。因此，加載前需將BIM與實景三維模型分別轉換為Cesium支持的3D-Tiles數據標準［14-15］。

實景三維模型由于引入了影像外方位元素、地面控制點，具有絕對坐標基準；而BIM主要用于建筑尺度的信息管理，未引入絕對位置參考，坐標往往采用局部坐標系表達。模型疊加可通過手工配準的方式實現，而顧及到GIS空間計算、分析的精度和準確性需要對BIM模型的位置、姿態及高程進行嚴密計算，故需通過坐標轉換的方式將BIM與實景三維模型的坐標基準進行統一。

BIM模型表現的是建筑設計方案，實景三維模型還原了影像采集時刻的地形地貌特征而未考慮到建筑施工所帶來的地形變化。直接將模型根據坐標進行疊加勢必會出現重疊和遮擋，兩者無法完美套合。因此需考慮到建筑施工所造成的地形變化，對模型銜接區域進行處理。

綜上所述，通過數據標準統一、坐標基準統一和模型銜接處理，即可將BIM與實景三維模型進行集成展示。本文重點介紹坐標基準統一方式、模型銜接處理的實現思路，而數據標準統一可通過FME軟件及開源程序實現，在此不做贅述。

2 模型坐標基準統一

坐標基準統一需通過對BIM的局部坐標根據坐標原點所在絕對坐標及坐標軸指向進行坐標轉換計算，轉換至地球中心坐標系（earth-centered earthfixed，ECEF），流程如圖2所示。

圖2 BIM坐標轉換流程Fig.2 BIM Coordinate Conversion Process

首先將局部坐標系旋轉至ENU（east-north-Up）坐標系（其中x軸指向南方向，y軸指向東方向，z軸指向天頂方向），再建立ENU坐標系到地球橢球ECEF坐標系之間的轉換關系，最后將ENU坐標進行旋轉和平移，得到ECEF世界坐標［16］。

1）局部坐標系旋轉到ENU坐標系。BIM在局部 坐 標 系 下 的 坐 標(x0，y0，z0)到ENU坐 標(x，y，z)，需依次繞三軸旋轉εz、εy、εx角度，對應的旋轉矩陣R(εx，εy，εz)為：

2）計算ENU到ECEF坐標系的平移量(X0，Y0，Z0)及旋轉矩陣R(B0，L0，0)分別為：

式中，BIM模型的坐標系原點位于橢球面上的大地坐標為(B0，L0)；N0表示卯酉圈半徑；a為橢球的長半軸；e為橢球第一偏心率。

3）最終轉換關系式為：

式 中，(x0，y0，z0)為BIM頂 點 的 局 部 坐 標；R(εx，εy，εz)為 局 部 坐 標 到ENU坐 標 旋 轉 矩 陣；(X0，Y0，Z0)和R(B0，L0，0)分 別 為ENU到ECEF坐標系下的平移量和旋轉矩陣。

本文根據該項目局部坐標系定義參數，通過Python對以上過程進行編程實現。在可視化時，將計算得到的旋轉矩陣和平移量合成為4階仿射變換矩陣，輸入到Cesium中3DTileset對象的ModelMatrix屬性中，即可在渲染時逐點進行坐標變換。實現效果如圖3所示，其中灰色部分即BIM所表現的立交橋，與實景三維模型所表現的周邊房屋、地面融為一體。相比于人工三維建模，BIM模型具備更加精細的建筑構件結構并可針對構件進行屬性關聯。

圖3 BIM與實景三維模型疊加效果Fig.3 BIM and Real-Scene 3D Model Overlay Effect

3 模型銜接處理

對于模型銜接區域的處理，本文采用對實景三維模型進行局部壓平的方式：在WebGL渲染過程中，對銜接區域的實景三維模型三角網頂點進行降低，將其統一至設計的地面高度。整個壓平處理的流程如圖4所示。

圖4 局部壓平流程示意圖Fig.4 Partial Flattening Process Diagram

通過修改WebGL渲染過程進行實現：事先采集疊加區域多邊形輪廓及目標高度，通過調用Cesium中的Framebuffer對象渲染一張矩形Texture（疊加區域內為灰度為1，多邊形外為0），將Texture與目標高度以WebGL Uniform屬性的形式傳遞到實景三維模型中的Vertex Shader（WebGL的頂點渲染腳本）中，通過判斷實景三維模型頂點坐標是否落在Texture灰度為1的區域內，進而將落在壓平區域內的頂點高度進行降低。

局部壓平的實現效果如圖5所示，對隧道口外的山體進行壓平，即可直觀表現隧道開挖后的實際效果。

圖5 局部壓平模擬開挖Fig.5 Partial Flattening Simulation Excavation

4 模型動態卷簾

為了在同一視圖中直觀對比BIM模型疊加前后效果，可通過動態卷簾的方式，對BIM與實景三維模型的顯示進行控制：對視圖進行左右分屏，在左屏中展示BIM與實景三維模型疊加效果、右屏中展示未疊加BIM的實景三維模型，并通過視圖分隔UI組件拖動對左右屏寬度比例進行調整。

動態卷簾的實現是通過改寫Cesium中3DTileset類的WebGL著色器（Shader）、配合第3節中介紹的頂點壓平的方式來實現：首先，將分屏UI組件所在位置的屏幕坐標換算到WebGL坐標系下，并分別傳遞給BIM和實景三維模型的Vertex Shader，屏幕坐標系和WebGL坐標系之間關系見文獻［13］；其次，重寫BIM模型的Fragment Shader（WebGL的片元渲染腳本），比較頂點gl_Position的水平坐標和分屏參數，如水平坐標大于分屏參數，通過Discard命令取消Fragment渲染；最后，重寫頂點壓平后的實景三維Vertex Shader，比較處于壓平區域內頂點的gl_Position水平坐標和分屏參數，如水平坐標小于分屏參數，則執行頂點壓平，反之則不執行。動態卷簾的實現效果如圖6所示，其中左側為BIM模型疊加在實景三維模型上的效果；右側為原始實景三維模型，分屏比例可拖動調整。

圖6 卷簾切換實現效果Fig.6 Roller Shutter Switching Effect

5 結束語

BIM與GIS的集成，能夠將建筑與環境信息融為一體，對室內外信息進行更加全面的數字化還原，對規劃設計行業中的方案展示、方案比對具有重要的作用。本文以磁浮鐵路項目為例，探討了BIM與GIS的數據集成可視化流程及思路；基于開源的Cesium引擎，對BIM與實景三維模型集成中的坐標基準統一、模型銜接處理等實現方式進行介紹；介紹了用于對BIM與實景三維模型直觀對比的動態卷簾實現方法。通過對BIM與GIS進行集成可視化，后續可實現更加高級的空間分析功能，進而為建筑設計、路橋施工、城市規劃乃至CIM的建設提供平臺支撐。