OpenGL：基于已知高程数据生成视点正确的透视地形山脊线的算法求助

方案1：预计算山脊纹理（性能最优，适配静态DTED数据） #

• 适用场景：高程数据为静态DTED文件，支持分块加载的地形网格
• 实现逻辑：
  1. 预处理阶段对DTED高程网格逐点计算山脊特征：对每个高程点，检测8邻域梯度方向，若当前点高程高于梯度垂直方向的两侧邻域点，且低于梯度下降方向的邻域点，判定为山脊点
  2. 将山脊特征值存储为与DTED分辨率匹配的单通道纹理，山脊点值为1，非山脊点为0
  3. 运行时地形渲染的片元着色器直接采样该纹理，对山脊点叠加指定高亮颜色即可
• 优势：无额外渲染Pass，运行时开销可忽略，山脊线精度与DTED分辨率一致，支持视距相关的淡出效果

方案2：单Pass后处理实时提取（适配动态视角，无预处理成本） #

• 适用场景：需要适配动态修改的高程数据，或不想做离线预处理的场景，解决此前仅用深度缓冲只能提取天际线的问题
• 实现逻辑：仅新增1个全屏后处理Pass，无需读回CPU，全程GPU计算
  1. 地形渲染时输出深度缓冲、法向量缓冲两张纹理（前向渲染可将法向量编码到额外渲染目标，无需完整延迟渲染管线）
  2. 后处理着色器中对当前像素采样4邻域的深度与法向量，加权计算深度差与法向量差，超过阈值即判定为山脊线
• 参考GLSL代码片段：

#version 330 core
out vec4 FragColor;
in vec2 TexCoords;

uniform sampler2D sceneTex;
uniform sampler2D depthTex;
uniform sampler2D normalTex;
uniform float ridgeThreshold = 0.15;
uniform vec3 ridgeColor = vec3(0.8, 0.1, 0.1);
uniform vec2 screenSize;

float getDepth(vec2 uv) {
    return texture(depthTex, uv).r;
}

vec3 getNormal(vec2 uv) {
    return normalize(texture(normalTex, uv).rgb * 2.0 - 1.0);
}

void main() {
    vec2 pixelOffset = 1.0 / screenSize;
    vec3 baseCol = texture(sceneTex, TexCoords).rgb;
    
    // 计算深度差拉普拉斯算子
    float dC = getDepth(TexCoords);
    float dL = getDepth(TexCoords + vec2(-pixelOffset.x, 0));
    float dR = getDepth(TexCoords + vec2(pixelOffset.x, 0));
    float dU = getDepth(TexCoords + vec2(0, pixelOffset.y));
    float dD = getDepth(TexCoords + vec2(0, -pixelOffset.y));
    float depthEdge = abs(dL + dR - 2 * dC) + abs(dU + dD - 2 * dC);
    
    // 计算法向差
    vec3 nC = getNormal(TexCoords);
    vec3 nL = getNormal(TexCoords + vec2(-pixelOffset.x, 0));
    vec3 nR = getNormal(TexCoords + vec2(pixelOffset.x, 0));
    float normalEdge = 1.0 - dot(nC, normalize(nL + nR));
    
    // 判定山脊线，加权两种特征
    float isRidge = step(ridgeThreshold, depthEdge * 0.3 + normalEdge * 0.7);
    FragColor = vec4(mix(baseCol, ridgeColor, isRidge * 0.7), 1.0);
}

• 优化技巧：可根据像素深度动态调整pixelOffset大小，远处山脊线用更大的偏移，保证不同距离山脊线宽度一致，避免远处山脊线过细消失

方案3：几何着色器生成矢量山脊线（无锯齿，适配矢量高亮需求） #

• 适用场景：需要山脊线为矢量形式，无分辨率缩放锯齿的场景
• 实现逻辑：
  1. 地形网格绘制时传入相邻三角形的法向量信息，或在几何着色器中查询相邻三角形面数据
  2. 对每个三角形的三条边进行判定：若边两侧的三角形法向量夹角超过阈值，且两个三角形分别朝向/背向视点，判定为山脊边
  3. 将山脊边扩展为指定宽度的四边形，输出高亮颜色即可
• 优势：生成的山脊线为矢量图形，缩放无锯齿，可自定义线宽、线型，运行时开销远低于通视计算或多Pass渲染

内容的提问来源于stack exchange，提问作者Fergus Crawford