基于Cass高程点生成三维地形的完整指南：从数据准备到可视化实现

1次阅读

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

Cass 格式是我国测绘行业常用的数据交换格式，其高程点数据通常以文本文件存储，包含点号、X/ Y 坐标、高程值等信息。这类数据的特点包括：

结构化存储 ：每行记录代表一个高程点，字段间用逗号或空格分隔
高密度采样 ：在地形复杂区域（如山脊、河谷）会有更密集的布点
工程友好 ：可直接用于土木工程设计、洪水模拟等专业场景

典型应用案例：

水利工程中的库容计算
城市规划的土方量估算
军事模拟的战场环境构建

算法类型	优点	缺点	适用场景
IDW（反距离加权）	计算简单，保留原始数据特征	易产生 ” 牛眼 ” 效应	数据分布均匀的平坦区域
Kriging（克里金）	可评估插值误差，统计最优	计算复杂度高	存在空间自相关的复杂地形
TIN（不规则三角网）	精确还原地形特征	需要后处理平滑	工程设计等高精度要求场景

# 数据读取与预处理
def load_cass_data(file_path):
    """解析 Cass 格式高程点文件"""
    points = []
    with open(file_path) as f:
        for line in f:
            # 示例数据格式：点号,X,Y,H
            parts = line.strip().split(',')
            if len(parts) >= 4:
                points.append((float(parts[1]), float(parts[2]), float(parts[3])))
    return np.array(points)

# 使用 PyQGIS 创建栅格表面
def create_raster_layer(points, cell_size=10, method='idw'):
    """生成 DEM 栅格数据"""
    extent = QgsRectangle(min(points[:,0]), min(points[:,1]),
        max(points[:,0]), max(points[:,1]))

    # 设置插值参数
    if method.lower() == 'idw':
        interpolation = QgsIdwInterpolator()
    else:
        interpolation = QgsTINInterpolator()

    # 执行插值计算
    raster = QgsGridFileWriter(
        interpolation, "/output/dem.tif",
        extent, cell_size, cell_size)
    raster.writeFile()

数据准备阶段
检查数据完整性（缺失值、异常值）
坐标系统一转换为 EPSG:4326 或当地投影坐标系
去除明显错误点（如高程值超过珠峰高度）
插值计算阶段
根据地形复杂度选择网格大小（平坦地形可用 20m 网格，山区建议 5 -10m）
对大面积区域建议分块处理
保存中间结果避免重复计算

可视化实现

# 使用 Plotly 生成 3D 曲面
import plotly.graph_objects as go

def plot_3d_surface(z_values):
    fig = go.Figure(data=[go.Surface(z=z_values)])
    fig.update_layout(
        title='三维地形模型',
        scene=dict(
            xaxis_title='经度',
            yaxis_title='纬度',
            zaxis_title='高程 (m)'))
    fig.show()