ezdxf实现DXF中bulge转ARC及闭合实体面积周长计算问题

报错分析与修复 #

调用ezdxf.math.bulge_to_arc()和bulge_center()时出现TypeError: object of type 'float' has no len()，原因是参数传递错误：这两个函数要求传入(bulge, 起点坐标元组, 终点坐标元组)，但代码误将单个坐标值（x/y）、线宽等无关参数传入，导致Vec2初始化失败。

修正后的bulge转换代码 #

import ezdxf
from ezdxf.math import bulge_to_arc, bulge_center

dwg = ezdxf.readfile("example.dxf")
for num, pline in enumerate(dwg.query('LWPOLYLINE'), start=1):
    print(f"\n{num}. POLYLINE:")
    vertices = list(pline.points())
    vertex_count = len(vertices)
    # 闭合多段线需循环到最后一个顶点（连接回起点）
    loop_range = vertex_count if pline.closed else vertex_count - 1
    
    for i in range(loop_range):
        x1, y1, _, _, bulge = vertices[i]
        # 获取当前线段的终点
        if i == vertex_count - 1:
            x2, y2, _, _, _ = vertices[0]  # 闭合多段线连接回起点
        else:
            x2, y2, _, _, _ = vertices[i+1]
        
        start_point = (x1, y1)
        end_point = (x2, y2)
        
        if bulge != 0:
            midpoint, radius, start_angle, end_angle = bulge_to_arc(bulge, start_point, end_point)
            center = bulge_center(bulge, start_point, end_point)
            print(f"\t弧段: 起点{start_point} -> 终点{end_point}")
            print(f"\t\t圆心: {center}, 半径: {radius:.2f}, 角度范围: {start_angle:.1f}°~{end_angle:.1f}°")
        else:
            print(f"\t直线段: 起点{start_point} -> 终点{end_point}")

两类闭合实体的面积与周长计算 #

1. 闭合LWPOLYLINE（含直线、圆弧、样条） #

推荐两种方法，根据实体类型选择：

方法一：用virtual_entities()分解为LINE/ARC计算

适合仅含直线和圆弧的LWPOLYLINE，直接分解为基础实体计算：

import math
from ezdxf import path

def calc_closed_lwpolyline(pline):
    perimeter = 0.0
    # 计算周长
    for entity in pline.virtual_entities():
        if entity.dxftype() == 'LINE':
            perimeter += entity.dxf.start.distance(entity.dxf.end)
        elif entity.dxftype() == 'ARC':
            angle_diff = abs(entity.dxf.end_angle - entity.dxf.start_angle)
            angle_diff = angle_diff if angle_diff <= 180 else 360 - angle_diff
            perimeter += (angle_diff / 360) * 2 * math.pi * entity.dxf.radius
    
    # 计算面积：用Path对象直接获取闭合区域面积
    p = path.make_path(pline)
    area = p.area
    return round(area, 2), round(perimeter, 2)

# 使用示例
dwg = ezdxf.readfile("example.dxf")
for pline in dwg.query('LWPOLYLINE'):
    if pline.closed:
        area, perimeter = calc_closed_lwpolyline(pline)
        print(f"闭合LWPOLYLINE - 面积: {area}, 周长: {perimeter}")

方法二：转为Path对象近似处理

适合包含样条曲线的LWPOLYLINE，通过展平路径为折线近似计算：

from ezdxf import path

def calc_lwpolyline_approx(pline, tolerance=0.01):
    p = path.make_path(pline)
    vertices = list(p.flattening(tolerance))
    # 计算周长
    perimeter = sum(vertices[i].distance(vertices[i+1]) for i in range(len(vertices)-1))
    # 鞋带公式计算多边形面积
    area = 0.0
    n = len(vertices)
    for i in range(n):
        x1, y1 = vertices[i]
        x2, y2 = vertices[(i+1)%n]
        area += (x1 * y2) - (x2 * y1)
    area = abs(area) / 2
    return round(area, 2), round(perimeter, 2)

# 使用示例
for pline in dwg.query('LWPOLYLINE'):
    if pline.closed:
        area, perimeter = calc_lwpolyline_approx(pline)
        print(f"闭合LWPOLYLINE(近似) - 面积: {area}, 周长: {perimeter}")

2. 首尾相连的分解轮廓（LINE/ARC组成） #

需先将首尾相连的实体分组为闭合轮廓，再计算：

import math
from collections import defaultdict
from ezdxf import path

def group_closed_contours(msp):
    start_map = defaultdict(list)
    end_map = defaultdict(list)
    all_entities = []
    
    # 构建实体的起点/终点映射
    for entity in msp.query('LINE, ARC'):
        all_entities.append(entity)
        start = tuple(entity.dxf.start)
        end = tuple(entity.dxf.end)
        start_map[start].append(entity)
        end_map[end].append(entity)
    
    visited = set()
    contours = []
    for entity in all_entities:
        if entity in visited:
            continue
        contour = [entity]
        visited.add(entity)
        current_end = tuple(entity.dxf.end)
        
        while True:
            # 查找下一个相连的实体
            next_entities = [e for e in start_map.get(current_end, []) if e not in visited]
            if not next_entities:
                # 检查是否闭合
                first_start = tuple(contour[0].dxf.start)
                if current_end == first_start:
                    contours.append(contour)
                break
            next_entity = next_entities[0]
            contour.append(next_entity)
            visited.add(next_entity)
            current_end = tuple(next_entity.dxf.end)
    return contours

def calc_contour_area_perimeter(contour):
    perimeter = 0.0
    p = path.Path()
    # 构建Path并计算周长
    for entity in contour:
        if entity.dxftype() == 'LINE':
            perimeter += entity.dxf.start.distance(entity.dxf.end)
            p.line_to(entity.dxf.end)
        elif entity.dxftype() == 'ARC':
            angle_diff = abs(entity.dxf.end_angle - entity.dxf.start_angle)
            angle_diff = angle_diff if angle_diff <= 180 else 360 - angle_diff
            perimeter += (angle_diff / 360) * 2 * math.pi * entity.dxf.radius
            p.arc_to(entity.dxf.end, radius=entity.dxf.radius, start_angle=entity.dxf.start_angle)
    p.close()
    area = p.area
    return round(area, 2), round(perimeter, 2)

# 使用示例
dwg = ezdxf.readfile("example.dxf")
msp = dwg.modelspace()
contours = group_closed_contours(msp)
for idx, contour in enumerate(contours, start=1):
    area, perimeter = calc_contour_area_perimeter(contour)
    print(f"闭合分解轮廓{idx} - 面积: {area}, 周长: {perimeter}")

方案总结 #

• 方案1（手动匹配起止坐标）：逻辑繁琐，易出错，不推荐。
• 方案2（利用ezdxf内置工具）：借助virtual_entities()、Path等工具简化计算，兼容性强，是最优选择。

内容的提问来源于stack exchange，提问作者Micul Geala