犀牛如何导入CAD等高线生成三维地形？

Question

在使用Rhino（犀牛）导入CAD等高线创建三维地形时，一个常见问题是：**等高线因高程属性缺失或Z值未正确赋值，导致生成的地形高度错误或完全平坦**。许多用户从AutoCAD导入DWG等高线后，仅将二维线条带入Rhino，却未将其标高信息转化为实际的Z坐标。这会导致“放样”或“地形网格”工具无法识别高度差异，最终生成错误模型。此外，等高线间距不一致、线条未闭合或存在冗余重叠，也会造成曲面扭曲或重建失败。如何准确提取CAD中的高程数据并为每条等高线赋予对应Z值，是实现精准三维地形建模的关键技术难点。

Answer 1

1. 问题背景与常见现象

在使用Rhino（犀牛）进行三维地形建模时，用户常从AutoCAD中导入DWG格式的等高线数据。然而，一个普遍且关键的问题是：这些二维等高线缺乏Z轴高程信息或未正确映射到三维空间中的Z坐标值。这导致即使使用“放样曲面”或“地形网格生成器”等工具，也无法构建出具有真实高度变化的地形模型，最终输出结果往往是完全平坦或高度错乱的曲面。

更复杂的情况包括：

• CAD文件中等高线标注了高程文本，但Rhino未自动识别；
• 等高线之间间距不一致，影响三角网插值精度；
• 存在重复、断裂或未闭合的曲线，造成拓扑错误；
• 不同图层对应不同高程，但未建立图层-高程映射关系。

2. 数据提取与预处理流程

为解决Z值缺失问题，必须首先从原始CAD数据中提取高程属性。常见的方法如下表所示：

提取方式	适用场景	工具支持	自动化程度
图层命名解析	每条等高线位于按高程命名的图层（如“Contour_100”）	Rhino + Python脚本	高
文本标签匹配	每条线旁有独立的高程文字标注	Rhino + Grasshopper + TextObject分析	中
块属性提取（Block Attribute）	等高线以带属性的块形式存在	DWG直接读取，需支持属性导出	较高
外部CSV/Excel对照表	已有明确的图层名与高程对照清单	Grasshopper数据映射组件	高

3. Z值赋值的技术实现路径

一旦获取高程数据，下一步是将每条等高线提升至对应的Z坐标平面。以下是典型实现步骤：

1. 导入DWG文件后分离所有等高线对象；
2. 遍历每个曲线对象，获取其所属图层名称或邻近文本内容；
3. 通过正则表达式提取数字部分作为候选高程值；
4. 构建图层→高程映射字典；
5. 使用RhinoScript或Grasshopper将每条曲线沿Z方向移动至目标高度；
6. 验证所有曲线是否处于唯一且合理的Z层级；
7. 清理重叠或冗余曲线（可借助布尔运算或距离容差筛选）；
8. 确保所有等高线形成嵌套结构（外环包围内环）；
9. 导出为Polyline集合用于后续TIN（不规则三角网）生成；
10. 调用Mesh.CreateFromContours或第三方插件（如Land for Rhino）生成地形网格。

4. Grasshopper可视化编程解决方案

对于具备一定参数化设计经验的5年以上IT/设计从业者，推荐采用Grasshopper进行自动化处理。以下是一个简化的逻辑流程图：

import re

def extract_elevation_from_layer(layer_name):
    # 示例：从图层名 "Contour_250.0" 提取 250.0
    match = re.search(r'(\d+\.?\d*)', layer_name)
    return float(match.group(1)) if match else 0.0

# 应用于Rhino.Python环境
import rhinoscriptsyntax as rs

objects = rs.GetObjects("Select contour lines", rs.filter.curve)
for obj in objects:
    layer = rs.ObjectLayer(obj)
    z = extract_elevation_from_layer(layer)
    bbox = rs.BoundingBox(obj)
    center = rs.PointDivide(rs.AddPoints(bbox[0], bbox[6]), 2)
    current_z = center[2]
    rs.MoveObject(obj, [0,0,z - current_z])

5. 流程图：完整地形建模工作流

mermaid.initialize({startOnLoad:true});
graph TD
    A[导入DWG等高线] --> B{是否存在高程标注？}
    B -- 是，基于图层命名 --> C[解析图层名获取Z值]
    B -- 是，基于文本标签 --> D[执行文本 proximity 匹配]
    B -- 否 --> E[手动输入高程对照表]
    C --> F[批量设置曲线Z坐标]
    D --> F
    E --> F
    F --> G[检查曲线闭合性与拓扑一致性]
    G --> H[去除重复/交叉线段]
    H --> I[生成等高线层级序列]
    I --> J[创建TIN或NURBS地形曲面]
    J --> K[导出为OBJ/DEM/Revit可读格式]