CAD等高线地形图转三维效果图：从数据处理到渲染的全流程实战

1次阅读

没有评论

共计 2973 个字符，预计需要花费 8 分钟才能阅读完成。

CAD 等高线地形图是工程设计和地理信息系统中常见的数据形式，但将其转换为三维效果图时面临诸多挑战。等高线数据通常以二维多段线形式存储，缺乏高程信息的直接表达，且可能存在以下问题：

数据精度不一致：不同来源的等高线数据精度差异较大
拓扑关系缺失：等高线之间缺乏明确的拓扑连接关系
数据量大：复杂地形区域的等高线数据可能包含数万个顶点
纹理映射困难：如何将航拍影像或卫星图像准确贴合到三维模型表面

这些挑战使得直接从 CAD 等高线生成高质量三维效果图成为一个非平凡的技术问题。

目前主要有三种主流方案可供选择：

商业软件方案（如 ArcGIS）
优点：一站式解决方案，图形界面友好
缺点：成本高，定制化能力有限
开源 GIS 工具链（如 QGIS+GRASS）
优点：免费，社区支持良好
缺点：处理大数据量时性能不足
自定义开发方案（Python+ 开源库）
优点：完全可控，可针对特定需求优化
缺点：开发成本较高

考虑到灵活性和成本因素，我们选择自定义开发方案，基于 Python 构建处理流水线。

我们使用 pyautocad 库与 AutoCAD 进行交互，提取等高线数据：

import pyautocad

def extract_contours(dwg_path):
    acad = pyautocad.Autocad(create_if_not_exists=True)
    acad.ActiveDocument.Application.Documents.Open(dwg_path)

    contours = []
    for entity in acad.iter_objects(['Polyline']):
        if is_contour(entity):  # 自定义函数判断是否为等高线
            points = list(entity.Coordinates)
            elevation = get_elevation(entity)  # 从图层名或扩展数据获取高程
            contours.append({
                'points': points,
                'elevation': elevation
            })
    return contours

使用 GDAL 将离散的等高线点插值为连续的 DEM：

from osgeo import gdal, osr
import numpy as np

def create_dem(contours, cell_size=5, power=2):
    # 1. 准备插值点数据
    xyz = []
    for contour in contours:
        points = contour['points']
        for i in range(0, len(points), 2):
            xyz.append((points[i], points[i+1], contour['elevation']))

    # 2. 使用 IDW 插值算法
    # ...（具体实现代码）# 3. 创建 GeoTIFF 文件
    driver = gdal.GetDriverByName('GTiff')
    ds = driver.Create('dem.tif', width, height, 1, gdal.GDT_Float32)
    ds.SetGeoTransform(geo_transform)
    ds.SetProjection(srs.ExportToWkt())
    ds.GetRasterBand(1).WriteArray(dem_array)
    ds = None

我们提供两种渲染方案：

// Three.js 地形渲染核心代码
const loader = new THREE.TextureLoader();
const texture = loader.load('texture.jpg');

const geometry = new THREE.PlaneGeometry(width, height, widthSegments, heightSegments);

// 将 DEM 高度数据应用到顶点
const vertices = geometry.attributes.position.array;
for (let i = 0; i < vertices.length; i += 3) {vertices[i + 2] = demData[i / 3] * verticalExaggeration;
}

const material = new THREE.MeshStandardMaterial({map: texture});
const terrain = new THREE.Mesh(geometry, material);
scene.add(terrain);

使用 Blender Python API 进行高质量渲染：

import bpy

def create_blender_terrain(dem_path, texture_path):
    # 1. 创建网格
    bpy.ops.mesh.primitive_grid_add(
        x_subdivisions=width_segments,
        y_subdivisions=height_segments,
        size=terrain_size
    )

    # 2. 应用高度图
    # ...（具体实现代码）# 3. 设置材质和光照
    mat = bpy.data.materials.new(name="TerrainMaterial")
    mat.use_nodes = True
    bsdf = mat.node_tree.nodes["Principled BSDF"]
    tex_node = mat.node_tree.nodes.new('ShaderNodeTexImage')
    tex_node.image = bpy.data.images.load(texture_path)
    mat.node_tree.links.new(bsdf.inputs['Base Color'], tex_node.outputs['Color'])

处理大规模地形数据时，需要考虑以下优化策略：