```python """ FreeCAD无界面点云曲线生成脚本 功能：直接读取包含XYZ坐标的文件，根据内容生成折线或B样条曲线 执行方式：在FreeCAD Python控制台运行，或作为宏执行 """ import FreeCAD import FreeCADGui import Part import Draft import os def create_curve_from_file(file_path, curve_type='polyline', closed=False): """ 核心函数：从文件读取点并创建曲线 参数: file_path (str): 点云数据文件路径 curve_type (str): 曲线类型，'polyline' 或 'bspline' closed (bool): 是否创建封闭曲线 返回: Part.Shape: 生成的曲线对象 """ # 1. 读取点云数据 [ref_1] points = [] try: with open(file_path, 'r', encoding='utf-8') as f: for line_num, line in enumerate(f, 1): line = line.strip() # 跳过空行和注释行 if not line or line.startswith('#'): continue # 支持多种分隔符（空格、逗号、制表符） parts = line.replace(',', ' ').split() if len(parts) >= 3: try: x = float(parts[0]) y = float(parts[1]) z = float(parts[2]) points.append(FreeCAD.Vector(x, y, z)) except ValueError as e: print(f"警告：第{line_num}行数据格式错误: {line}") continue else: print(f"警告：第{line_num}行数据不足3列: {line}") except FileNotFoundError: print(f"错误：文件未找到 - {file_path}") return None except Exception as e: print(f"错误：读取文件时发生异常 - {str(e)}") return None # 2. 验证点数据 if len(points) < 2: print(f"错误：有效点数量不足({len(points)})，至少需要2个点") return None print(f"成功读取 {len(points)} 个点") # 3. 创建曲线 [ref_1] try: if curve_type.lower() == 'bspline': # 创建B样条曲线（需要至少4个点） if len(points) < 4: print("警告：B样条需要至少4个点，自动切换到折线模式") curve_type = 'polyline' else: bspline = Part.BSplineCurve() bspline.interpolate(points) # 插值生成B样条 curve_shape = bspline.toShape() if closed: # 创建封闭B样条 points.append(points[0]) # 添加起点作为终点 bspline_closed = Part.BSplineCurve() bspline_closed.interpolate(points) curve_shape = bspline_closed.toShape() print("已创建B样条曲线") return curve_shape if curve_type.lower() == 'polyline': # 创建折线 if closed: points.append(points[0]) # 添加起点作为终点以封闭 curve_shape = Part.makePolygon(points) # [ref_1] print("已创建折线") return curve_shape except Exception as e: print(f"错误：创建曲线时发生异常 - {str(e)}") return None def main(): """主执行函数""" # 配置参数 DAT_FILE = "line.dat" # 点云数据文件 [ref_1] CURVE_TYPE = "polyline" # 曲线类型: 'polyline' 或 'bspline' [ref_1] CLOSED_CURVE = False # 是否封闭曲线 [ref_1] CREATE_IN_NEW_DOC = True # 是否在新文档中创建 # 1. 检查文件是否存在 if not os.path.exists(DAT_FILE): print(f"错误：文件 {DAT_FILE} 不存在") print("请确保数据文件位于以下目录之一:") print(f" 1. 当前工作目录: {os.getcwd()}") print(f" 2. FreeCAD宏目录") return # 2. 创建或获取活动文档 if CREATE_IN_NEW_DOC: doc = FreeCAD.newDocument("PointCloudCurve") print("已创建新文档: PointCloudCurve") else: doc = FreeCAD.ActiveDocument if not doc: doc = FreeCAD.newDocument() print("已创建新文档") else: print(f"使用现有文档: {doc.Name}") # 3. 生成曲线 print(f"开始处理文件: {DAT_FILE}") print(f"曲线类型: {CURVE_TYPE}") print(f"封闭曲线: {CLOSED_CURVE}") curve_shape = create_curve_from_file(DAT_FILE, CURVE_TYPE, CLOSED_CURVE) if curve_shape: # 4. 将曲线添加到文档 curve_obj = doc.addObject("Part::Feature", f"{CURVE_TYPE}_curve") curve_obj.Shape = curve_shape # 5. 设置显示属性 if CURVE_TYPE == 'bspline': curve_obj.ViewObject.LineColor = (0.0, 0.0, 1.0) # 蓝色 else: curve_obj.ViewObject.LineColor = (0.0, 1.0, 0.0) # 绿色 curve_obj.ViewObject.LineWidth = 2.0 curve_obj.ViewObject.PointSize = 3.0 # 6. 更新文档和视图 doc.recompute() # 7. 自动调整视图 FreeCADGui.SendMsgToActiveView("ViewFit") print("✓ 曲线创建成功！") print(f" 文档: {doc.Name}") print(f" 对象: {curve_obj.Name}") print(f" 曲线类型: {CURVE_TYPE}") print(f" 点数量: {len(curve_shape.Vertexes) if hasattr(curve_shape, 'Vertexes') else 'N/A'}") else: print("✗ 曲线创建失败") # 批量处理版本 def batch_process_dat_files(folder_path, curve_type='polyline'): """ 批量处理文件夹中的所有.dat文件 参数: folder_path (str): 文件夹路径 curve_type (str): 曲线类型 """ import glob if not os.path.exists(folder_path): print(f"错误：文件夹不存在 - {folder_path}") return # 获取所有.dat文件 dat_files = glob.glob(os.path.join(folder_path, "*.dat")) if not dat_files: print(f"警告：在 {folder_path} 中未找到.dat文件") return print(f"找到 {len(dat_files)} 个.dat文件，开始批量处理...") # 创建新文档用于批量处理 batch_doc = FreeCAD.newDocument("BatchCurves") results = [] for file_path in dat_files: file_name = os.path.basename(file_path) print(f"\n处理文件: {file_name}") curve_shape = create_curve_from_file(file_path, curve_type) if curve_shape: # 添加到文档 obj_name = f"{curve_type}_{os.path.splitext(file_name)[0]}" curve_obj = batch_doc.addObject("Part::Feature", obj_name) curve_obj.Shape = curve_shape results.append({ 'file': file_name, 'object': obj_name, 'points': len(curve_shape.Vertexes) if hasattr(curve_shape, 'Vertexes') else 0, 'success': True }) print(f" ✓ 成功创建: {obj_name}") else: results.append({ 'file': file_name, 'success': False, 'error': '曲线生成失败' }) print(f" ✗ 失败") # 更新文档 batch_doc.recompute() FreeCADGui.SendMsgToActiveView("ViewFit") # 输出汇总 print(f"\n{'='*50}") print("批量处理完成汇总:") print(f"{'='*50}") success_count = sum(1 for r in results if r['success']) print(f"总文件数: {len(results)}") print(f"成功: {success_count}") print(f"失败: {len(results) - success_count}") return results # 执行主函数 if __name__ == "__main__": print("开始执行无界面点云曲线生成脚本...") main() ``` ## 详细使用指南 ### 1. **数据文件格式要求** | 格式类型 | 示例 | 说明 | |---------|------|------| | 空格分隔 | `0 0 0`<br>`10 5 2`<br>`20 10 4` | 每行三个数值，用空格分隔 [ref_1] | | 逗号分隔 | `0,0,0`<br>`10,5,2`<br>`20,10,4` | 每行三个数值，用逗号分隔 | | 混合格式 | `0 0,0`<br>`10,5 2`<br>`20 10,4` | 脚本会自动处理多种分隔符 | | 带注释 | `# 起始点`<br>`0 0 0`<br>`# 控制点`<br>`10 5 2` | 以#开头的行会被忽略 | ### 2. **脚本配置参数详解** ```python # 主要配置参数说明 CONFIG = { "DAT_FILE": "line.dat", # 数据文件名，支持绝对或相对路径 [ref_1] "CURVE_TYPE": "polyline", # 可选: 'polyline'(折线) 或 'bspline'(B样条) [ref_1] "CLOSED_CURVE": False, # 是否创建封闭曲线 [ref_1] "CREATE_IN_NEW_DOC": True, # 是否在新文档中创建 "ENABLE_PREPROCESS": False, # 是否启用点云预处理 "PREPROCESS_METHOD": "simple" # 预处理方法: 'simple', 'distance', 'grid' } ``` ### 3. **不同执行方式对比** | 执行方式 | 命令/操作 | 适用场景 | 优点 | |---------|----------|---------|------| | FreeCAD控制台 | 复制脚本到Python控制台执行 | 快速测试、调试 | 实时反馈，便于调试 | | 宏执行 | 保存为.FCMacro文件后执行 | 重复使用、自动化 | 一键执行，方便管理 | | 命令行 | `FreeCADCmd -c "exec(open('script.py').read())"` | 批量处理、集成 | 无需GUI，适合自动化 | | 模块导入 | `import script; script.main()` | 复杂项目集成 | 模块化，可重用 | ### 4. **点云预处理功能** ```python def preprocess_point_cloud(points, method='simple'): """ 点云数据预处理函数 参数: points: FreeCAD.Vector点列表 method: 预处理方法 支持的预处理方法对比: """ # 预处理方法对比表 methods_comparison = { 'simple': { '功能': '简单去重', '适用场景': '去除完全重复的点', '处理效果': '保留所有不重复的点', '时间复杂度': 'O(n)' }, 'distance': { '功能': '距离过滤', '适用场景': '去除过于密集的点', '处理效果': '保留间距大于阈值的点', '时间复杂度': 'O(n)' }, 'grid': { '功能': '网格简化', '适用场景': '大规模点云简化', '处理效果': '网格内点合并为平均值', '时间复杂度': 'O(n)' } } # 实际处理逻辑 if method == 'distance': # 距离过滤实现 if len(points) < 2: return points processed_points = [points[0]] min_distance = 1.0 # 最小距离阈值 for i in range(1, len(points)): last_point = processed_points[-1] current_point = points[i] distance = last_point.distanceToPoint(current_point) if distance >= min_distance: processed_points.append(current_point) print(f"距离过滤: {len(points)} -> {len(processed_points)} 点") return processed_points # 其他方法实现... ``` ### 5. **应用场景示例** #### 场景1：CAD设计自动化 ```python def cad_automation_pipeline(input_file, output_format='step'): """ CAD设计自动化流程：从点云数据到3D模型 流程步骤: 1. 读取点云数据 2. 生成曲线 3. 创建曲面 4. 导出为CAD格式 """ # 1. 读取点云数据 [ref_2] points = [] with open(input_file, 'r') as f: for line in f: coords = line.strip().split() if len(coords) >= 3: x, y, z = map(float, coords[:3]) points.append(FreeCAD.Vector(x, y, z)) # 2. 生成曲线 [ref_1] if len(points) >= 4: # 使用B样条生成平滑曲线 bspline = Part.BSplineCurve() bspline.interpolate(points) curve = bspline.toShape() else: # 点较少时使用折线 curve = Part.makePolygon(points) # 3. 创建曲面（示例：旋转曲面） axis = FreeCAD.Vector(0, 0, 1) # Z轴 surface = curve.revolve(axis, 360) # 旋转360度生成曲面 # 4. 导出 if output_format.lower() == 'step': surface.exportStep("output.step") elif output_format.lower() == 'stl': surface.exportStl("output.stl") return surface ``` #### 场景2：科学数据可视化 ```python def scientific_data_visualization(data_file, visualization_type='3d_curve'): """ 科学实验数据可视化 支持的数据类型: - 3D轨迹数据（粒子运动轨迹） - 测量点云（激光扫描数据） - 实验曲线（物理实验数据） """ # 读取科学数据 [ref_2] scientific_data = { 'points': [], 'attributes': [], # 可存储额外属性如温度、压力等 'metadata': {} # 元数据 } with open(data_file, 'r') as f: for line_num, line in enumerate(f, 1): # 支持多种科学数据格式 if line.startswith('#METADATA:'): # 解析元数据 key_value = line.replace('#METADATA:', '').strip().split('=') if len(key_value) == 2: scientific_data['metadata'][key_value[0]] = key_value[1] elif not line.startswith('#'): # 解析数据点 parts = line.strip().split() if len(parts) >= 3: x, y, z = map(float, parts[:3]) scientific_data['points'].append(FreeCAD.Vector(x, y, z)) if len(parts) > 3: scientific_data['attributes'].append(parts[3:]) # 根据数据类型选择可视化方式 if visualization_type == '3d_curve': # 生成3D曲线 curve = create_curve_from_points(scientific_data['points'], 'bspline') # 根据属性设置曲线颜色（示例：根据温度） if scientific_data['attributes']: # 这里可以添加根据属性设置颜色的逻辑 pass elif visualization_type == 'point_cloud': # 生成点云显示 point_cloud = Part.makeCompound([ Part.Vertex(p) for p in scientific_data['points'] ]) return curve ``` ### 6. **错误处理与调试** | 常见错误 | 原因分析 | 解决方案 | |---------|---------|---------| | 文件未找到 | 文件路径错误或文件不存在 | 检查文件路径，使用绝对路径 | | 数据格式错误 | 文件包含非数字字符 | 清理数据文件，移除注释和空行 | | 内存不足 | 点云数据量过大 | 启用预处理，减少点数 | | FreeCAD未运行 | 脚本在FreeCAD外部执行 | 在FreeCAD Python控制台内执行 | | 权限问题 | 文件只读或无写入权限 | 检查文件权限，更改保存目录 | ### 7. **性能优化建议** 1. **大规模数据处理策略**： ```python def process_large_point_cloud(file_path, chunk_size=10000): """分块处理大规模点云数据""" points = [] chunk_count = 0 with open(file_path, 'r') as f: for line in f: # 解析点数据 # ... points.append(point) # 每处理chunk_size个点生成一个曲线段 if len(points) >= chunk_size: curve_segment = Part.makePolygon(points) # 保存或处理曲线段 points = [] # 清空列表 chunk_count += 1 return chunk_count ``` 2. **内存管理技巧**： - 使用生成器逐行读取大文件 - 及时释放不再使用的对象 - 使用`del`