利用Blender和DEM(TIF)数据高效生成三维模型的实战指南

1次阅读

没有评论

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

在实际项目中，我们经常需要将 DEM（数字高程模型）数据转换为三维模型，用于地形可视化、游戏开发或 GIS 应用。传统的处理方法存在几个明显的痛点：

内存占用高 ：大范围的 DEM 数据直接加载会导致内存溢出
处理速度慢 ：手动处理步骤繁琐，耗时较长
模型细节丢失 ：简单的插值算法会导致地形特征模糊
坐标系统混乱 ：不同来源的 DEM 数据坐标系不统一

Blender 提供了强大的 Python API，我们可以通过脚本实现全流程自动化：

数据预处理（坐标系转换、裁剪）
网格生成与优化
材质自动映射
模型导出

坐标系转换 ：使用 pyproj 库进行 WGS84 到 UTM 的转换
数据裁剪 ：根据目标区域范围裁剪 DEM，减少处理量
空值处理 ：识别并填充 NODATA 值，避免模型空洞

采样精度 ：决定地形细节层级
Z 轴缩放系数 ：控制地形起伏程度
光滑处理 ：影响模型表面平滑度

以下是完整的 Python 脚本示例（Blender 3.0+）：

import bpy
import numpy as np
from osgeo import gdal
import pyproj

# 1. DEM 数据加载与预处理
def load_dem(file_path):
    dataset = gdal.Open(file_path)
    band = dataset.GetRasterBand(1)
    elevation = band.ReadAsArray()

    # 坐标系转换示例（WGS84 转 UTM）src_proj = pyproj.Proj(init='epsg:4326')
    dst_proj = pyproj.Proj(init='epsg:32651')
    # ... 坐标转换代码...

    return elevation

# 2. 生成地形网格
def create_terrain(elevation, z_scale=1.0):
    x_dim, y_dim = elevation.shape

    # 创建网格
    mesh = bpy.data.meshes.new('TerrainMesh')
    obj = bpy.data.objects.new('Terrain', mesh)
    bpy.context.collection.objects.link(obj)

    # 构建顶点
    verts = [(x, y, elevation[x][y]*z_scale) 
             for y in range(y_dim) for x in range(x_dim)]

    # 构建面
    faces = []
    for y in range(y_dim-1):
        for x in range(x_dim-1):
            idx = x + y * x_dim
            faces.append((idx, idx+1, idx+x_dim+1, idx+x_dim))

    mesh.from_pydata(verts, [], faces)
    mesh.update()

    return obj

# 3. 自动材质映射
def assign_material(obj):
    mat = bpy.data.materials.new(name='TerrainMaterial')
    obj.data.materials.append(mat)
    # ... 材质设置代码...

# 主流程
if __name__ == '__main__':
    dem_data = load_dem('input.tif')
    terrain = create_terrain(dem_data, z_scale=0.5)
    assign_material(terrain)