Xbim获取圆形IIfcBuildingElementProxy的BoundingBox尺寸偏差10mm问题求助

问题背景 #

使用VB.NET结合Xbim库获取IIfcBuildingElementProxy尺寸时，矩形对象的尺寸显示正常，但圆形对象的尺寸偏差约10mm，存在异常。

核心原因 #

1. 曲线几何体细分精度不足：Xbim将圆形等曲线几何体转换为多边形网格时，默认细分精度较低，导致边界盒计算出现误差
2. 单位转换遗漏：IFC文件可能使用米/厘米作为单位，未转换为毫米会导致数值偏差
3. 边界盒的局限性：圆形的外接矩形边界盒本身是近似值，若能直接读取IFC原生属性会更准确

修正方案 #

1. 提高几何体细分精度

创建几何上下文时，设置更小的细分长度参数，减少曲线转多边形的误差：

' 初始化几何上下文时配置细分精度
Dim contextOptions As New XbimGeometryOptions()
contextOptions.MaximumTessellationLength = 1 ' 数值越小，细分越精细，精度越高
Dim context As New Xbim3DModelContext(model, contextOptions)
context.CreateContext()

2. 强制单位转换为毫米

检查IFC文件的长度单位，统一转换为毫米：

' 获取模型单位转换因子（默认转换为毫米）
Dim conversionFactor As Double = 1
Dim unitAssignment = model.Instances.FirstOrDefault(Of IIfcUnitAssignment)()
If unitAssignment IsNot Nothing Then
    Dim lengthUnit = unitAssignment.Units.OfType(Of IIfcLengthUnit).FirstOrDefault()
    If lengthUnit IsNot Nothing Then
        conversionFactor = lengthUnit.UnitType switch
            Case IfcUnitEnum.METRE : 1000
            Case IfcUnitEnum.CENTIMETRE : 10
            Case Else : 1
        End Select
    End If
End If

3. 优先读取IFC原生属性

如果圆形对象的尺寸在IFC属性集中定义，直接读取属性值比计算边界盒更准确：

' 遍历属性集，尝试读取直径/半径属性
For Each propSetRel In BEproxy.IsDefinedBy.OfType(Of IIfcRelDefinesByProperties)()
    Dim propSet = propSetRel.RelatingPropertyDefinition
    If TypeOf propSet Is IIfcPropertySet Then
        Dim ps = CType(propSet, IIfcPropertySet)
        ' 匹配属性名（根据实际IFC文件调整，比如"Diameter"或"Radius"）
        Dim diameterProp = ps.HasProperties.OfType(Of IIfcPropertySingleValue)().FirstOrDefault(Function(p) p.Name.Equals("Diameter", StringComparison.OrdinalIgnoreCase))
        If diameterProp IsNot Nothing Then
            Dim diameter = CDbl(diameterProp.NominalValue.Value) * conversionFactor
            MsgBox($"BuildingElementProxy: {BEproxy.Name}")
            MsgBox($"直径：{Math.Round(diameter,2)}mm")
            ' 读取到原生属性后跳过边界盒计算
            hasNativeSize = True
            Exit For
        End If
    End If
Next

修正后的完整代码 #

Sub Main1()
    ' 加载IFC文件
    Dim ifcFilePath As String = "C:\IFC\2 sparingen.ifc"

    ' 初始化模型
    Using model As IfcStore = IfcStore.Open(ifcFilePath)
        ' 配置几何细分精度
        Dim contextOptions As New XbimGeometryOptions()
        contextOptions.MaximumTessellationLength = 1
        Dim context As New Xbim3DModelContext(model, contextOptions)
        context.CreateContext()

        ' 获取所有IIfcBuildingElementProxy对象
        Dim buildingElmProxy = model.Instances.OfType(Of IIfcBuildingElementProxy)()

        ' 计算单位转换因子（转毫米）
        Dim conversionFactor As Double = 1
        Dim unitAssignment = model.Instances.FirstOrDefault(Of IIfcUnitAssignment)()
        If unitAssignment IsNot Nothing Then
            Dim lengthUnit = unitAssignment.Units.OfType(Of IIfcLengthUnit).FirstOrDefault()
            If lengthUnit IsNot Nothing Then
                conversionFactor = lengthUnit.UnitType switch
                    Case IfcUnitEnum.METRE : 1000
                    Case IfcUnitEnum.CENTIMETRE : 10
                    Case Else : 1
                End Select
            End If
        End If

        ' 遍历每个BuildingElementProxy
        For Each BEproxy In buildingElmProxy
            Dim hasNativeSize As Boolean = False

            ' 优先读取原生属性尺寸
            For Each propSetRel In BEproxy.IsDefinedBy.OfType(Of IIfcRelDefinesByProperties)()
                Dim propSet = propSetRel.RelatingPropertyDefinition
                If TypeOf propSet Is IIfcPropertySet Then
                    Dim ps = CType(propSet, IIfcPropertySet)
                    Dim diameterProp = ps.HasProperties.OfType(Of IIfcPropertySingleValue)().FirstOrDefault(Function(p) p.Name.Equals("Diameter", StringComparison.OrdinalIgnoreCase))
                    If diameterProp IsNot Nothing Then
                        Dim diameter = CDbl(diameterProp.NominalValue.Value) * conversionFactor
                        MsgBox($"BuildingElementProxy: {BEproxy.Name}")
                        MsgBox($"直径：{Math.Round(diameter,2)}mm")
                        hasNativeSize = True
                        Exit For
                    End If
                End If
            Next

            If hasNativeSize Then Continue For

            ' 通过边界盒计算尺寸
            Dim geom = context.ShapeInstances.Where(Function(si) si.IfcProductLabel = BEproxy.EntityLabel)
            Dim bbox As XbimRect3D = XbimRect3D.Empty
            For Each shape In geom
                bbox.Union(shape.BoundingBox)
            Next

            ' 转换单位并输出
            Dim width As Double = Math.Round(bbox.SizeX * conversionFactor, 2)
            Dim height As Double = Math.Round(bbox.SizeY * conversionFactor, 2)
            Dim depth As Double = Math.Round(bbox.SizeZ * conversionFactor, 2)

            MsgBox($"BuildingElementProxy: {BEproxy.Name}")
            MsgBox($"Width: {width}mm, Height: {height}mm, Depth: {depth}mm")
            Console.WriteLine()
        Next
    End Using
End Sub

额外检查项 #

• 验证IFC文件中的圆形对象是否为标准参数化圆形，还是由小平面拼接的近似图形
• 更新Xbim库到最新版本，修复已知的几何计算bug

内容的提问来源于stack exchange，提问作者Roy TImmers