Blender通过SHP生成三维地形的完整解决方案与避坑指南

1次阅读

没有评论

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

在 GIS 数据可视化领域，将二维的 SHP 格式数据转换为三维地形时，开发者常遇到几个典型问题：

坐标系转换错误：源数据采用地理坐标系（如 WGS84）而 Blender 使用笛卡尔坐标系，直接导入会导致严重变形
高程数据丢失：SHP 属性表中的高程值未被正确解析，导致生成的地形平坦化
纹理拉伸问题：常规 UV 展开方式在山地区域会产生严重拉伸，需特殊处理
性能瓶颈：大规模地形数据直接生成会导致内存溢出和视图卡顿

目前主流方案有三种实现路径：

Blender GIS 插件
优点：可视化操作，无需编程
缺点：对复杂 CRS 支持有限，无法处理自定义高程字段
Python 脚本处理
优点：灵活控制每个处理环节，支持自动化流水线
缺点：需要熟悉 Blender Python API
第三方工具链
优点：专业 GIS 软件（如 QGIS）预处理能力强
缺点：数据流转环节多，易产生精度损失

推荐采用 Python 脚本 +GDAL 的组合方案，兼具灵活性和精度控制。

使用 GDAL 进行关键转换：

from osgeo import ogr

def load_shp_elevation(shp_path):
    """提取 SHP 中的高程信息"""
    datasource = ogr.Open(shp_path)
    layer = datasource.GetLayer()

    # 获取高程字段（示例字段名为 'ELEV'）elevations = []
    for feature in layer:
        elevations.append(feature.GetField('ELEV'))
    return np.array(elevations)

通过 bpy 模块创建地形网格：

import bpy
import numpy as np

def create_terrain_mesh(vertices, faces):
    """生成地形网格"""
    mesh = bpy.data.meshes.new('Terrain')
    mesh.from_pydata(vertices, [], faces)

    # 自动计算法线
    mesh.update(calc_edges=True)

    # 创建材质
    mat = bpy.data.materials.new('TerrainMaterial')
    mesh.materials.append(mat)

    return mesh

基于坡度的材质混合：在着色器节点中使用 Geometry 节点的 Normal 信息驱动不同材质混合
投影 UV 映射：采用正交投影方式减少山地纹理拉伸
微表面细节：通过 Bump 节点增加小尺度地形起伏感

# 地形生成主流程（PEP8 规范）import bpy
import numpy as np
from osgeo import ogr

def generate_terrain(shp_path, texture_path):
    """
    从 SHP 生成带纹理的三维地形
    :param shp_path: SHP 文件路径
    :param texture_path: 纹理图片路径
    """
    # 1. 提取高程数据
    elevations = load_shp_elevation(shp_path)

    # 2. 创建网格顶点（示例为简化逻辑）x = np.linspace(0, 100, len(elevations))
    z = elevations
    vertices = [(xi, 0, zi) for xi, zi in zip(x, z)]

    # 3. 生成 Blender 对象
    mesh = create_terrain_mesh(vertices, [])
    obj = bpy.data.objects.new('Terrain', mesh)
    bpy.context.collection.objects.link(obj)

    # 4. 配置材质（关键步骤）setup_terrain_material(obj, texture_path)

    return obj