Мой подход заключается в вычислении двух касательных векторов, параллельных осям X и Y соответственно. Затем вычислите векторное произведение, чтобы найти вектор нормали.
Касательный вектор задается линией, которая пересекает среднюю точку на двух ближайших сегментах, как показано на следующем рисунке.
Мне было интересно, есть ли более прямой расчет или менее затратный с точки зрения циклов ЦП.
На самом деле вы можете рассчитать его без перекрестного произведения, используя «метод конечных разностей» (или, по крайней мере, я думаю, что он называется таким образом).
На самом деле это достаточно быстро, чтобы я использовал его для вычисления нормалей на лету в вершинном шейдере.
// # P.xy store the position for which we want to calculate the normals
// # height() here is a function that return the height at a point in the terrain
// read neightbor heights using an arbitrary small offset
vec3 off = vec3(1.0, 1.0, 0.0);
float hL = height(P.xy - off.xz);
float hR = height(P.xy + off.xz);
float hD = height(P.xy - off.zy);
float hU = height(P.xy + off.zy);
// deduce terrain normal
N.x = hL - hR;
N.y = hD - hU;
N.z = 2.0;
N = normalize(N);
