Вначале хочу сказать про матрицы. Совершенно заморочили голову с этими преобразованиями. Так вот, матрицы нужны для шейдерного движка при расчете сцены. Если мы обсчитываем координаты одной фигуры в неподвижном окружении, то использование матриц не принесет ускорения расчета и понимание кода.
Я не зря написал в комментарии про "Правильно заданный вопрос". Действительно, именно обдуманная и исправленная формулировка вопроса помогла найти корректный ответ (C++):
// исходное направление камеры: вдоль оси Х в сторону +
// Построим перед камерой плоский прямоугольник, немного
// сместив его (x = t) навстречу камере. После "сдвига" координат
// при помощи сложения с вектором ViewTo это смещение
// компенсируется.
double s = 0.01, t = -0.9;
std::vector<glm::vec3> p {
glm::vec3(t,-s,-s),
glm::vec3(t,-s, s),
glm::vec3(t, s, s),
glm::vec3(t, s,-s)
};
// каждую вершину поворачиваем на углы положения камеры и смещаем
// на вектор ViewTo (к точке, на которую направлена камера).
for (size_t i = 0; i < 4; i++)
{
p[i] = glm::rotateZ(p[i], look_zenith);
// азимут отрицательный из-за реализации поворота в
// в библиотеке GLM, не совпавшей с моей моделью вращения.
p[i] = glm::rotateY(p[i], 0 - look_azimuth);
p[i] += ViewTo;
}

На изображении белый прямоугольник и есть плоскость перед камерой, которая всегда неподвижна относительно камеры.