Возьмите линейку, карандаш или ручку и начертите на бумаге линию. У вас выйдет как‑то так:

Вы скажете, что сложного в том, чтобы нарисовать линию? А теперь попробуйте нарисовать с помощью того же инструмента линию шириной в 5–10 раз больше ширины кегля. Уже сложнее не так ли?:

Ладно, скажете вы, у меня есть корректор, который гораздо толще ручки! Да, вы упростите себе задачу, но если нужно будет нарисовать линию толще? Будет снова не просто:

Вы придёте к выводу, что тонкая линия лишь кажется тонкой, если её приблизить, можно увидеть ширину и погрешность, что и при рисовании линии большей ширины. А всё дело в том, что линии это не геометрические объекты (т. е. не 1d объекты).

Все линии, которые мы можем создавать имеют толщину, а линии с толщиной это уже определённые плоскости (опустим тот факт, что в нашем мире всё имеет не нулевые размеры во всех трёх плоскостях). Но поскольку, чаще всего нам надо иметь однородную и симметричную линию, то линию можно спокойно представить как прямоугольник. А это значит, что задача нарисовать линию, сводится к задаче нарисовать прямоугольник.

В качестве визуализации процесса будем использовать OpenGL, а именно геометрические шейдеры на языке GLSL. Всё что они делают — это берут на вход точки фигуры, а далее строят из них фигуры. Весь их код будет прокомментирован, поэтому знания GLSL не обязательны для понимания процесса.

Для построения этого прямоугольника нам известны координаты центров его сторон, расположенных друг напротив друга. А так же ширина этих сторон. Найти же нам надо, координаты четырёх углов прямоугольника.

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

Точки P0 и P2 очень похожи на Ц0 и Ц1, однако имеют смещения по осям. Поэтому их можно выразить напрямую через эти смещения и координаты точек Ц0 и Ц1:

Координаты же точек P1 и P2 получить ещё проще — через симметрию относительно Ц0 и Ц1:

В результате нам надо найти лишь смещения x и y и через него можно будет получить координаты остальных точек. А сделать это можно как минимум двумя разными способами. Первый — требует минимальных знаний геометрии, второй — знаний работы с векторами. Какой из них проще для понимания решать вам, поэтому опишу их оба.

Метод на основе длины линий

Начнём с метода, для которого нужны только знания теоремы Пифагора. В чём суть данного метода. Мы можем получить расстояние до точки P0 двумя разными способами, а из этого получить два уравнения, решив которые, мы получим её координаты. Вычтя из них координаты точки Ц0, мы и получим данное смещение.

Для начала нарисуем рисунок:

Нас интересуют две прямые L и M. Их длину через исходные данные можно получить вот так:

{L} ^ {2} = {( {x}_ {1} - {x} _ {0} )} ^ {2} + {( {y} _ {1} - {y} _ {0} )} ^ {2}{M} ^ {2} = {L} ^ {2} + {S} ^ {2} = {( {x} _ {1} - {x} _ {0} )} ^ {2} + {( {y} _ {1} - {y} _ {0} )} ^ {2} + {S} ^ {2}

Или через интересующею нас точку вот так:

\begin{cases}{M}^{2}={\left({x}_{1}-{\left({x}_{}+{x}_{0}\right)}_{}\right)}^{2}{+\left({y}_{1}-{{\left({y}_{}+{y}_{0}\right)}_{}}_{}\right)}^{2}\\  {S}^{2}={\left({x}_{0}-{\left({x}_{}+{x}_{0}\right)}_{}\right)}^{2}+{\left({y}_{0}-{{\left({y}_{}+{y}_{0}\right)}_{}}_{}\right)}^{2}\end{cases}

Два уравнения, две переменных значит, решение точно есть. Осталось лишь его найти. Для строгости, приведу полной этап их решения, тем более он тривиален. Раскроем скобки:

\begin{cases}{M}^{2}={\left({x}_{1}-{x}_{0}{-x}_{}\right)}^{2}{+\left({y}_{1}-{y}_{0}-y\right)}^{2}\\  {S}^{2}={\left(-x\right)}^{2}+{\left(-y\right)}^{2}\end{cases}

И ещё одни:

\begin{cases}{M}^{2}={{x}_{1}}^{2}-{2x}_{1}{x}_{0}-{2x}_{1}x+{2x}_{0}x+{{x}_{0}}^{2}+{{x}_{}}^{2}+ \\ + {{y}_{1}}^{2}-{2y}_{1}{y}_{0}-{2y}_{1}y+{2y}_{0}y+{{y}_{0}}^{2}+{{y}_{}}^{2}\\  {S}^{2}={{x}_{}}^{2}+ {y}^{2}\end{cases}

Заметим, что можно выразить квадрат длины расстояния между исходными точками:

\begin{cases}{{x}_{1}}^{2}-{2x}_{1}{x}_{0}+{{x}_{0}}^{2}={\left({x}_{1}-{x}_{0}{}_{}\right)}^{2}\\  {{y}_{1}}^{2}-{2y}_{1}{y}_{0}+{{y}_{0}}^{2}={\left({y}_{1}-{y}_{0}{}_{}\right)}^{2}\end{cases}

В итоге получим:

\begin{cases}{M}^{2}-{L}^{2}=-{2x}_{1}x+{2x}_{0}x+{{x}_{}}^{2}-{2y}_{1}y+{2y}_{0}y+{{y}_{}}^{2}\\  {S}^{2}={{x}_{}}^{2}+{y}^{2}\end{cases}

А благодаря этому можно немного упростить выражение:

\begin{cases}{S}^{2}=-{2x}_{1}x+{2x}_{0}x+{{x}_{}}^{2}-{2y}_{1}y+{2y}_{0}y+{{y}_{}}^{2}\\  {S}^{2}={{x}_{}}^{2}+{y}^{2}\end{cases}

Теперь мы можем приравнять обе части:

{{x}_{}}^{2}+{y}^{2}=-{2x}_{1}x+{2x}_{0}x+{{x}_{}}^{2}-{2y}_{1}y+{2y}_{0}y+{{y}_{}}^{2}

Уберём повторы, и объединим одинаковые переменные, заодно и знаки поменяем:

0={x}_{1}x-{x}_{0}x+{y}_{1}y-{y}_{0}yx\left({x}_{1}-{x}_{0}\right)+y\left({y}_{1}-{y}_{0}\right)=0

Выразим x и y через друг-друга:

\begin{cases}x=\frac{-y\left({y}_{1}-{y}_{0}\right)}{\left({x}_{1}-{x}_{0}\right)}\\  y=\frac{-x\left({x}_{1}-{x}_{0}\right)}{\left({y}_{1}-{y}_{0}\right)}\end{cases}

Подставим в исходное уравнение:

{S}^{2}={{\left(\frac{-y\left({y}_{1}-{y}_{0}\right)}{\left({x}_{1}-{x}_{0}\right)}\right)}_{}}^{2}+{y}^{2}

Приведём к нулю:

\left(1+{\left(\frac{{y}_{1}-{y}_{0}}{{x}_{1}-{x}_{0}}\right)}^{2}\right){y}^{2}-{S}^{2}=0

Решим его, для x:

x=\pm S\sqrt{\frac{1}{1+{\left(\frac{{x}_{1}-{x}_{0}}{{y}_{1}-{y}_{0}}\right)}^{2}}}

Аналогично, для y:

y=\pm S\sqrt{\frac{1}{1+{\left(\frac{{y}_{1}-{y}_{0}}{{x}_{1}-{x}_{0}}\right)}^{2}}}

Их тоже можно упростить, для этого раскроем скобки:

x=\pm S\sqrt{\frac{{\left({y}_{1}-{y}_{0}\right)}^{2}}{{\left({x}_{1}-{x}_{0}\right)}^{2}+{\left({y}_{1}-{y}_{0}\right)}^{2}}}y=\pm S\sqrt{\frac{{\left({x}_{1}-{x}_{0}\right)}^{2}}{{\left({x}_{1}-{x}_{0}\right)}^{2}+{\left({y}_{1}-{y}_{0}\right)}^{2}}}

Заметим, что корень из квадрата это модуль, а делитель это расстояние между исходными точками:

x=\pm S\frac{\left\vert {y}_{1}-{y}_{0}\right\vert }{L}y=\pm S\frac{\left\vert {x}_{1}-{x}_{0}\right\vert }{{L}^{}}

В целом, так как у каждого выражения есть два решения, отличных только по знаку,а S и L всегда положительные, модуль можно и убрать:

x=\pm S\frac{{y}_{1}-{y}_{0}}{L}y=\pm S\frac{{x}_{1}-{x}_{0}}{{L}^{}}

В итоге, мы получили два уравнения, остаётся лишь вопрос, какое из двух решений каждого уравнения соответствует нашему сдвигу относительно изначальных точек.

Для этого рассмотрим четыре положения прямоугольника‑линии в пространстве, относительно центра. Так как точки симметричны, что было показано ранее, покажем влияние лишь на точку P0 :

Так как на знак, влияет лишь разница между координатами исходных точек, то  проанализируем их знаки во всех возможных положениях точки P0:

Как можно заметить из таблицы, смещение для y в точности советует знаку разницы координат, значит уравнение для смещение y будет иметь положительный знак для точки P0 (а для точки Pв силу симметрии отрицательный знак).

Для x ситуация ровно обратная, знак смещения всегда иной, чем знак разницы, значит уравнение смещения x будет с отрицательным знаком для точки P0 (а для точки Pв силу симметрии положительный знак).

В итоге окончательно получим смещение:

\begin{cases}x=-\frac S {{L}}( {y}_{1}-{y}_{0})\\ y=\frac S {{L}}( {x}_{1}-{x}_{0}) \end{cases}

Следующий метод приведёт нас к такому же результату, но совсем другим путём...

Метод на основе ортогональности

Точки на плоскости можно представить как вектора, характеризующийся величиной (длина вектора) и направлением от начала координат. Тогда, пусть точки P0 ,Ц0 и Ц1 это вектора P0,Ц0 и Ц1 направленные от начала координат.

Этот метод основан на том, что P0 и Ц0 ортогональны (т. е. между ними угол 90 градусов). Значит, мы можем повернуть вектор полученный из разницы векторов Ц0 и Ц1 на 90 градусов, и получить вектор Ц2 , такой же длины как расстояние от Ц0 до Ц1. А так же мы знаем расстоянии от P0 до Ц0 . Поэтому для определения вектора P0 достаточно найти соотношения длин векторов, и умножить на него координаты Ц2 . А уже через координаты P0 найти смещение.

Для начала найдём вектор разницы Ц0 и Ц1:

Для полученного вектора Ц2 используем формулы поворота вектора на угол на плоскости :

\begin{cases}{x}^{\text{'}}=\left({x}_{1}-{x}_{0}\right)\ast \cos а-\left({y}_{1}-{y}_{0}\right)\ast \sin а\\{y}^{\text{'}}=\left({x}_{1}-{x}_{0}\right)\ast \sin а+\left({y}_{1}-{y}_{0}\right)\ast \cos а\end{cases}

Подставим в формулу 90 градусов для того чтобы из ортогональному вектору P0 , получить сонаправленный вектор:

\begin{cases}{x}^{\text{'}}=-\left({y}_{1}-{y}_{0}\right)\\{y}^{\text{'}}=\left({x}_{1}-{x}_{0}\right)\end{cases}

В итоге получим следующие определение координат точек:

Далее найдём соотношение длин вектора Ц2 и заданной ширины S:

{K}^{}=\frac{S}{L}=\frac{S}{\sqrt{{\left({x}_{1}-{x}_{0}\right)}^{2}+{\left({y}_{1}-{y}_{0}\right)}^{2}}}

Осталось только изменить длину вектора Ц2 , не меняя его направления, для этого умножим его координаты на соотношение длин:

\begin{cases}{x}^{}=-\left({y}_{1}-{y}_{0}\right)K\\{y}^{}=\left({x}_{1}-{x}_{0}\right)K\end{cases}

Поставим на место К его уравнение:

\begin{cases}x=-\frac S {{L}}( {y}_{1}-{y}_{0})\\ y=\frac S {{L}}( {x}_{1}-{x}_{0}) \end{cases}

Мы снова пришли к тем же формулам что и в прошлый раз!

Визуализация результата

И вот у нас есть всё необходимое для того чтобы нарисовать линию. Теперь визуализируем это, заодно собрав всё вместе. Для этого нам нужен конвейер рендера, который состоит из 6 этапов. Каждый из которых обрабатывает свои данные (подробнее в этом цикле статей) Из них мы можем влиять только на три этапа (выделены синим):

Рисунок взят из https://learnopengl.com/Getting-started/Hello-Triangle

Для нашего случая вершинный шейдер просто перекладывает вершины, перемещая их в вершинный шейдер. Матрица модели нужна для позиционирования всей линии в пространстве(подробнее):

#version 330 core

//координаты точек линии
layout (location = 0) in float position_x;
layout (location = 1) in float position_y;

//матрица модели, задающая положения  объекта в пространстве
uniform mat4 model;

void main(){
   //передача позиции точек в геометрический шейдер
    gl_Position = model *vec4(position_x,position_y,0.0f,1.0f);
}

Фрагментный шейдер ещё проще, он задаёт один цвет всей линии:

#version 330 core

out vec4 color;
uniform  vec4 mat_color;

void main(){
    color=mat_color;//передаём цвет в растеризатор
};

Любые 2d объекты могут быть представлены как группа треугольников (с некоторой погрешностью, конечно). Поэтому мы будем рисовать линию с помощью них. Геометрический шейдер получает две точки (сегмент линии) и на основе них строит прямоугольник из двух треугольников. Так как вывод из геометрического шейдера возможен только в виде triangle_strip (каждая вершина после второй создаёт новый треугольник), вершины будут представлены в таком порядке:

Сам код геометрического шейдера:

#version 400 core
layout (lines) in;//входные данные, две точки (x0;y0) и (x1;y1)

//входные данные, точки прямоугольника
//(прямоугольник рисуется с помощью двух треугольников,
//strip значит, что для построения второго треугольника 
//будут использованны две вершины предыдущего)
layout (triangle_strip, max_vertices = 4) out;

in  float S[];//передаём ширину линии

//Функция определения смещения
vec2 get_general_point_corner(vec2 P0,vec2 P1){
    float L =(distance(P0,P1));

    float K=S[0]/L;

    float x_=-K*(P1.y-P0.y);
    float y_=K*(P1.x-P0.x);

    return vec2(x_,y_);
}

void bild_quard(vec2 A_0, vec2 P0, vec2 P1){
    // 1:bottom-left
    gl_Position =vec4((A_0.x+P0.x), (A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 2:bottom-right
    gl_Position =vec4((-A_0.x+P0.x),  (-A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 3:top-left
    gl_Position =vec4((A_0.x+P1.x), (A_0.y+P1.y), 0.0, 1.0);    
    EmitVertex();
    // 4:top-right
    gl_Position =vec4((-A_0.x+P1.x),  (-A_0.y+P1.y), 0.0, 1.0);  
    EmitVertex();
}

void main() {
    //получаем смещения
    vec2 p_= get_general_point_corner(gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
    //строим прямоугольник-линию
    bild_quard(p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
    EndPrimitive();
}


Как запустить пример

Если вы плохо знаете OpenGl, но хотите сами запустить этот пример, используйте информацию из этой статьи. (пример из неё очень просто адаптировать для данного случая). Если это вызовет проблемы, я дополню статью подробным методом (для строительства линии будем использовать каждую новую вершину для нового сегмента).

В итоге получим следующие:

Проблема стыка двух линий

С одной линией всё нормально, попробуем нарисовать кривую из разных линий:

Как мы видим, на стыке сегментов есть проблема — они не соединены, нужно отдельно реализовать соединение линий.В принципе можно выделить три типа соединений:

  • Косое

  • Прямое

  • Закруглённое

Для реализации каждого из них нам будет нужна информация о вершинах других сегментов. Для этого вместо layout (lines) in  в шейдере используем layout (lines_adjacency) in, что даёт  информацию об вершинах других сегментах.

Косое соединение

В целом, мы можем просто соединить точки P2 и P4 следующим образом:

Для этого мы получаем в геометрическом шейдере смещение для другого прямоугольника, и находим через него точку P4. Можно было бы найти точно, с какой стороны нам рисовать треугольник для плавного стыка, на основе наклона двух линий, но лучше нарисовать лишний треугольник, чем добавлять условие. Для этого найдём точку P5 , и нарисуем прямоугольник P2P4P3P5.

Разбиение на треугольники его вершин будет таким же:

Код шейдера, на основе вышесказанного будет выглядеть вот так:

#version 400 core
//входные данные, четыре точки (x0;y0),(x1;y1),(x2;y2) и (x3;y4)
layout (lines_adjacency) in;

//входные данные, точки прямоугольника
//(прямоугольник рисуется с помощью двух треугольников,
//strip значит, что для построения второго треугольника 
//будут использованны две вершины предыдущего)
layout (triangle_strip, max_vertices = 6) out;

in  float S[];//передаём ширину линии

//Функция определения смешения
vec2 get_general_point_corner(vec2 P0,vec2 P1){
    float L =(distance(P0,P1));

    float K=S[0]/L;

    float x_=-K*(P1.y-P0.y);
    float y_=K*(P1.x-P0.x);

    return vec2(x_,y_);
}


void bild_quard(vec2 A_0, vec2 P0, vec2 P1){
    // 1:bottom-left
    gl_Position =vec4((A_0.x+P0.x), (A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 2:bottom-right
    gl_Position =vec4((-A_0.x+P0.x),  (-A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 3:top-left
    gl_Position =vec4((A_0.x+P1.x), (A_0.y+P1.y), 0.0, 1.0);    
    EmitVertex();
    // 4:top-right
    gl_Position =vec4((-A_0.x+P1.x),  (-A_0.y+P1.y), 0.0, 1.0);  
    EmitVertex();
}

void oblique_connection(vec2 A_1, vec2 P1){
    // верхний треугольник в стыке
    gl_Position =vec4((A_1.x+P1.x), (A_1.y+P1.y), 0.0, 1.0);       
    EmitVertex();
    // нижний треугольник в стыке
    gl_Position =vec4((-A_1.x+P1.x), (-A_1.y+P1.y), 0.0, 1.0);       
    EmitVertex();
}

void main() {
    //получаем смещения
    vec2 p_= get_general_point_corner(gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
    vec2 p_1=get_general_point_corner(gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);
     //строим прямоугольник-линию
    bild_quard(p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
   
    //строим соединение стыка линий
    oblique_connection(p_1,gl_in[1].gl_Position.xy);
    EndPrimitive();
}

В итоге выйдет вот так:

Это самый эффективный способ по числу треугольников и описанию способ. Однако, возможны и иные решения этой проблемы.

Прямое соединение

Для этого типа разрешения стыка, нам необходимо определить точку P6 :

К счастью, сделать это очень просто, так как наши сегменты линии находятся на одной плоскости, то всегда можно найти точку пресечения двух прямых линий или сказать что они параллельны. А значит, для того чтобы найти эту точку, мы должны найти по точкам P0 и P2 первое уравнение прямой линии, и по точкам P4 и P5 второе уравнение, а потом их приравнять.

Для начала запишем уравнение прямой линии через две точки для P0 и P2 :

\frac{{y}_{2}-{y}_{0}}{{x}_{2}-{x}_{0}}\ast x+\frac{{x}_{2}\ast {y}_{0}-{x}_{0}\ast {y}_{2}}{{x}_{0}-{x}_{0}}=y

Перепишем его через смещения и координаты изначальных точек линии (т.е. центров), для этого просто раскроем переменные:

\begin{cases}{x}_{0}={x}^{\text{'}}+{{x}_{0}}^{\text{'}\text{'}}\\  {x}_{2}={x}^{\text{'}}+{{x}_{1}}^{\text{'}\text{'}}\end{cases}\begin{cases}{y}_{0}={y}^{\text{'}}+{{y}_{0}}^{\text{'}\text{'}}\\  {y}_{2}={y}^{\text{'}}+{{y}_{1}}^{\text{'}\text{'}}\end{cases}\begin{cases}{x}_{1}=-S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{0}}^{\text{'}\text{'}}\\  {x}_{2}=-S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{1}}^{\text{'}\text{'}}\end{cases}\begin{cases}{y}_{1}=S\frac{{x}_{1}-{x}_{0}}{{L}^{}}+{{y}_{0}}^{\text{'}\text{'}}\\  {y}_{2}=S\frac{{x}_{1}-{x}_{0}}{{L}^{}}+{{y}_{1}}^{\text{'}\text{'}}\end{cases}

В итоге получается страшное уравнение:

\frac{S\frac{{x}_{1}-{x}_{0}}{{L}^{}}+{{y}_{1}}^{\text{'}\text{'}}-S\frac{{x}_{1}-{x}_{0}}{{L}^{}}-{{y}_{0}}^{\text{'}\text{'}}}{S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{0}}^{\text{'}\text{'}}-S\frac{{y}_{1}-{y}_{0}}{L}-{{x}_{1}}^{\text{'}\text{'}}}\ast x+\\+\frac{\left(-S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{1}}^{\text{'}\text{'}}\right)\ast \left(S\frac{{x}_{1}-{x}_{0}}{{L}^{}}+{{y}_{0}}^{\text{'}\text{'}}\right)-\left(-S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{0}}^{\text{'}\text{'}}\right)\ast \left(S\frac{{x}_{1}-{x}_{0}}{{L}^{}}+{{y}_{1}}^{\text{'}\text{'}}\right)}{S\frac{{y}_{1}-{y}_{0}}{L}+{{x}_{0}}^{\text{'}\text{'}}-S\frac{{y}_{1}-{y}_{0}}{L}-{{x}_{1}}^{\text{'}\text{'}}}=y

Сократим лишние:

\frac{{{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}\ast x+\frac{\left({{x}_{}}^{\text{'}}+{{x}_{1}}^{\text{'}\text{'}}\right)\ast \left({{y}_{}}^{\text{'}}+{{y}_{0}}^{\text{'}\text{'}}\right)+\left({{x}_{}}^{\text{'}}+{{x}_{0}}^{\text{'}\text{'}}\right)\ast \left({{y}_{}}^{\text{'}}+{{y}_{1}}^{\text{'}\text{'}}\right)}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

Раскроем скобки:

\frac{{{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}\ast x+\\+\frac{{{x}_{}}^{\text{'}}{{y}_{}}^{\text{'}}+{{x}_{1}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}+{{x}_{}}^{\text{'}}{{y}_{0}}^{\text{'}\text{'}}+{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{}}^{\text{'}}{{y}_{}}^{\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}-{{x}_{}}^{\text{'}}{{y}_{1}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

И ещё раз упростим:

\frac{{{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}\ast x+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}+{{x}_{}}^{\text{'}}{{y}_{0}}^{\text{'}\text{'}}+{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}-{{x}_{}}^{\text{'}}{{y}_{1}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

Разделим его на три части:

\frac{{{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}\ast x+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}+{{x}_{}}^{\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{}}^{\text{'}}-{{x}_{}}^{\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

Умножим первую и вторую часть на S/L :

\frac{\frac{S}{L}({{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}})}{\frac{S}{L}({{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}})}\ast x+\frac{{{y}_{}}^{\text{'}}\left({{x}_{1}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}\right)\frac{S}{L}+{{x}_{}}^{\text{'}}\left(y_{0}\text{'}\text{'}-{{y}_{1}}^{\text{'}\text{'}}\right)\frac{S}{L}}{\frac{S}{L}\left({{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}\right)}+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

Поменяем знаки в скобках:

\frac{\frac{S}{L}({{y}_{1}}^{\text{'}\text{'}}-{{y}_{0}}^{\text{'}\text{'}})}{\frac{S}{L}({{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}})}\ast x+\frac{{{y}_{}}^{\text{'}}\left({{x}_{1}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}\right)\frac{S}{L}-{{x}_{}}^{\text{'}}\left(y_{1}\text{'}\text{'}-{{y}_{0}}^{\text{'}\text{'}}\right)\frac{S}{L}}{-\frac{S}{L}\left({{x}_{1}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}\right)}+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

И окончательно упростим:

\frac{-{{x}_{}}^{\text{'}}}{{{y}_{}}^{\text{'}}}\ast x-\frac{{{{y}_{}}^{\text{'}}}^{2}+{{{x}_{}}^{\text{'}}}^{2}}{{{y}_{}}^{\text{'}}}+\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}=y

Перепишем части уравнения ка отдельные константы:

\begin{cases}a=\frac{-{{x}_{}}^{\text{'}}}{{{y}_{}}^{\text{'}}}\\  b=\frac{{{{y}_{}}^{\text{'}}}^{2}+{{{x}_{}}^{\text{'}}}^{2}}{{{y}_{}}^{\text{'}}}\\  c=\frac{{{x}_{1}}^{\text{'}\text{'}}{{y}_{0}}^{\text{'}\text{'}}-{{x}_{0}}^{\text{'}\text{'}}{{y}_{1}}^{\text{'}\text{'}}}{{{x}_{0}}^{\text{'}\text{'}}-{{x}_{1}}^{\text{'}\text{'}}}\end{cases}

Тогда уравнение примет вид:

a\ast x-b+c=y

Зачем мы это сделали?

Во‑первых, мы получили два константных компонента. Первый, c показывает смещение относительно центра координат до центра центральной линии(посмотрев на уравнение линии через две точки, мы увидим точно такую же константу, которая будет является единственной, и показывать смещение).

Во‑вторых, b показывает смещение от центральной линии прямоугольника до границ. Знак говорит о том, с какой стороны находится прямая линия заданная уравнением, слева(‑) или справа(+).

Далее, приравниваем уравнение:

{a}_{0}\ast x-{b}_{0}+{c}_{0}={a}_{1}\ast x-{b}_{1}+{c}_{1}

И находим координату x (y можно будет найти через любое другое уравнение):

x=\frac{-{b}_{1}+{c}_{1}-\left(-{b}_{0}+{c}_{0}\right)}{\left({a}_{0}-{a}_{1}\right)}

Как уже упоминалось в методе про косой скос, использовать условные выражения в шейдере может быть менее эффективно, чем добавление дополнительных треугольников. Поэтому мы найдём дополнительно точку пересечения через прямые линии проложенные через две других стороны прямоугольника. Благо, это просто изменение знака у переменной b на обратный, согласно второму факту.

Рисовать треугольники для этого стыка мы будем в таком порядке, из‑за требования каждая точка после второй — новый треугольник. И тут, в любом случае придётся строить дополнительные точки, даже если учитывать направление скоса, поэтому мы и будем строить его зеркально, располагая вершины вот так:

Я думаю код шейдера уже пояснять не нужно, тут тоже всё довольно просто. Из нового лишь функция, определяющая положение точки P6 (P9 ) :

#version 400 core
//входные данные, четыре точки (x0;y0),(x1;y1),(x2;y2) и (x3;y4)
layout (lines_adjacency) in;

//входные данные, точки прямоугольника
//(прямоугольник рисуется с помощью двух треугольников,
//strip значит, что для построения второго треугольника 
//будут использованны две вершины предыдущего)
layout (triangle_strip, max_vertices = 6) out;

in  float S[];//передаём ширину линии

//Функция определения смешения
vec2 get_general_point_corner(vec2 P0,vec2 P1){
    float L =(distance(P0,P1));

    float K=S[0]/L;

    float x_=-K*(P1.y-P0.y);
    float y_=K*(P1.x-P0.x);

    return vec2(x_,y_);
}


void bild_quard(vec2 A_0, vec2 P0, vec2 P1){
    // 1:bottom-left
    gl_Position =vec4((A_0.x+P0.x), (A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 2:bottom-right
    gl_Position =vec4((-A_0.x+P0.x),  (-A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 3:top-left
    gl_Position =vec4((A_0.x+P1.x), (A_0.y+P1.y), 0.0, 1.0);    
    EmitVertex();
    // 4:top-right
    gl_Position =vec4((-A_0.x+P1.x),  (-A_0.y+P1.y), 0.0, 1.0);  
    EmitVertex();
}

vec2 get_point_corner(int side,
                      vec2 A_0,vec2 P0,vec2 P1,
                      vec2 A_1,vec2 P2,vec2 P3){

    float a_0= -A_0.x/A_0.y;
    float a_1= -A_1.x/A_1.y;

    float b_0 = side*(A_0.x*A_0.x+A_0.y*A_0.y)/A_0.y + (P1.x*P0.y-P0.x*P1.y)/(P1.x-P0.x);
    float b_1 = side*(A_1.x*A_1.x+A_1.y*A_1.y)/A_1.y + (P3.x*P2.y-P2.x*P3.y)/(P3.x-P2.x);

    float x=-(b_0-b_1)/(a_0-a_1);
    float y=a_0*x+b_0;

    return vec2(x,y);
}

void direct_connection(vec2 T_0, vec2 T_1,vec2 A_1, vec2 P1){
    gl_Position =vec4(T_0.xy, 0.0, 1.0);
    EmitVertex();
    gl_Position =vec4(T_1.xy, 0.0, 1.0);
    EmitVertex();
    gl_Position =vec4((A_1.x+P1.x), (A_1.y+P1.y), 0.0, 1.0);
    EmitVertex();
    gl_Position =vec4((-A_1.x+P1.x), (-A_1.y+P1.y), 0.0, 1.0);
    EmitVertex();
}

void main() {
    //получаем смещения
    vec2 p_= get_general_point_corner(gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
    vec2 p_1=get_general_point_corner(gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);

    //ищём точку угла
    vec2 t_= get_point_corner( 1,p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy,                 
                                                  p_1,gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);
    vec2 t_1=get_point_corner(-1,p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy,                    
                                      p_1,gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);

     //строим прямоугольник-линию
    bild_quard(p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);

   
    //строим соединение стыка линий
    direct_connection(t_,t_1,p_1,gl_in[1].gl_Position.xy);
    EndPrimitive();
}

В итоге на экране мы увидим:

Уже лучше, чем прошлый метод, и вполне хороший результат. Но можно поступить ещё интереснее, и сделать закругление линии.

Закруглённое соединение

В предыдущем методе мы определили точки для прямого стыка. А в методе косого стыка, мы использовали точки на разных линиях для создания косой линии. В этом методе нам понадобится все три точки для использования квадратной кривой Безье (подробнее в прекрасной статье). Этот метод является промежуточной между обоими стыками способом.

Чтобы нарисовать данную кривую нам понадобится её уравнение:

{B\left(t\right)=\left(t-1\right)}^{2}{P}_{0}+2\ast t\ast \left(t-1\right){P}_{1}+{t}^{2}{P}_{2}

где t в диапазоне 0 до 1.

Достаточно менять t с определённым шагом, который и определит точность изгиба (чем больше шаг, тем меньше точность). Тут так же надо будет отзеркалить точки, для удобного рисования точек:

Как можно заметить, чтобы нарисовать такой парный изгиб, достаточно по‑очереди добавлять точки из обоих линий. На этом и основан шейдер для рендера этого случая:.

#version 400 core
//входные данные, четыре точки (x0;y0),(x1;y1),(x2;y2) и (x3;y4)
layout (lines_adjacency) in;

//входные данные, точки прямоугольника
//(прямоугольник рисуется с помощью двух треугольников,
//strip значит, что для построения второго треугольника 
//будут использованны две вершины предыдущего)
layout (triangle_strip, max_vertices = 48) out;

in  float S[];//передаём ширину линии

//Функция определения смешения
vec2 get_general_point_corner(vec2 P0,vec2 P1){
    float L =(distance(P0,P1));

    float K=S[0]/L;

    float x_=-K*(P1.y-P0.y);
    float y_=K*(P1.x-P0.x);

    return vec2(x_,y_);
}


void bild_quard(vec2 A_0, vec2 P0, vec2 P1){
    // 1:bottom-left
    gl_Position =vec4((A_0.x+P0.x), (A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 2:bottom-right
    gl_Position =vec4((-A_0.x+P0.x),  (-A_0.y+P0.y), 0.0, 1.0);    
    EmitVertex();
    // 3:top-left
    gl_Position =vec4((A_0.x+P1.x), (A_0.y+P1.y), 0.0, 1.0);    
    EmitVertex();
    // 4:top-right
    gl_Position =vec4((-A_0.x+P1.x),  (-A_0.y+P1.y), 0.0, 1.0);  
    EmitVertex();
}

vec2 get_point_corner(int side,
                      vec2 A_0,vec2 P0,vec2 P1,
                      vec2 A_1,vec2 P2,vec2 P3){

    float a_0= -A_0.x/A_0.y;
    float a_1= -A_1.x/A_1.y;

    float b_0 = side*(A_0.x*A_0.x+A_0.y*A_0.y)/A_0.y + (P1.x*P0.y-P0.x*P1.y)/(P1.x-P0.x);
    float b_1 = side*(A_1.x*A_1.x+A_1.y*A_1.y)/A_1.y + (P3.x*P2.y-P2.x*P3.y)/(P3.x-P2.x);

    float x=-(b_0-b_1)/(a_0-a_1);
    float y=a_0*x+b_0;

    return vec2(x,y);
}

void bezier_connection(vec2 T_0, vec2 T_1,vec2 A_0,vec2 A_1, vec2 P1){
    vec2 H0=vec2( A_0.x+P1.x, A_0.y+P1.y);
    vec2 H1=vec2(-A_0.x+P1.x,-A_0.y+P1.y);

    vec2 H2=vec2( A_1.x+P1.x, A_1.y+P1.y);
    vec2 H3=vec2(-A_1.x+P1.x,-A_1.y+P1.y);

    float t=0.0;
    int count=16;
    for(int z = 0; z < count+1; z++){
        vec2 B_0=(1-t)*(1-t)*H0+2*t*(1-t)*T_0+t*t*H2;
        vec2 B_1=(1-t)*(1-t)*H1+2*t*(1-t)*T_1+t*t*H3;

        gl_Position =vec4(B_0.xy, 0.0, 1.0);
        EmitVertex();
        gl_Position =vec4(B_1.xy, 0.0, 1.0);
        EmitVertex();

        t+=1.0/count;
    }
}


void main() {
    //получаем смещения
    vec2 p_= get_general_point_corner(gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
    vec2 p_1=get_general_point_corner(gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);

    //ищём точку угла
    vec2 t_= get_point_corner( 1,p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy,                 
                                                  p_1,gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);
    vec2 t_1=get_point_corner(-1,p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy,                    
                                      p_1,gl_in[1].gl_Position.xy,gl_in[2].gl_Position.xy);

     //строим прямоугольник-линию
    bild_quard(p_,gl_in[0].gl_Position.xy,gl_in[1].gl_Position.xy);
   
   //строим соединение cтыка линий
    bezier_connection(t_,t_1,p_,p_1,gl_in[1].gl_Position.xy);
    EndPrimitive();
}

И в итоге получим:

Как можно заметить, каждый следующий метод стыка, использует также промежуточные расчёты предыдущего, что наглядно показывает почему был выбран именно этот порядок описания методов.

Вместо заключения

В нашей жизни есть много вещей, которые на первый взгляд кажутся простыми, например, линия. Но если присмотреться, то всё становится гораздо интереснее. Даже чтобы провести простую линию, чуть толще, чем инструмент, нам нужно вспомнить основы геометрии.

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

Комментарии (34)


  1. Newpson
    29.04.2024 16:30
    +6

    Это самый эффективный способ по числу треугольников и описанию способ. Однако, возможны и иные решения этой проблемы.

    Есть ещё один вариант. Он описан здесь, я его реализовывал на OpenGL ES 2.0, писал небольшое приложение-рисовалку для Android, которое принимало по USB данные с графического планшета и определённым образом записывало их в буфер, чтобы отрисовку точек и все расчёты можно было делать прямо на GPU почти без участия CPU.

    Давние конспекты по матлогике


  1. turchan
    29.04.2024 16:30
    +16

    Можно ли описанными способами решить классическую задачу? ;-)
    Семь красных линий, все они должны быть строго перпендикулярны, и кроме того, некоторые нужно нарисовать зеленым цветом, а еще некоторые — прозрачным.


    1. martein
      29.04.2024 16:30
      +2

      "Ну если без дураков, то без дураков выполним..." (с) "Фитиль" 1976г


    1. rsashka
      29.04.2024 16:30
      +29

      Вот так?


      1. PsihXMak
        29.04.2024 16:30

        "Просили же 7 красных линий. А у вас их 5. Да ещё и бледные и кривые какие то. Короче, всё переделывай."


        1. whoisking
          29.04.2024 16:30

          И прозрачных там 5, а не 3


    1. Iustinianus
      29.04.2024 16:30
      +4

      1. Семь перпендикулярных линий без проблем рисуются в неевклидовом пространстве.

      2. Можно нарисовать зеленую линию на цветном фоне, которая будет восприниматься красной - особенности человеческого цветовосприятия (точную пару цветов, увы, не вспомню).

      3. Делаем коэффициент отражения чернил комплексным числом - вуаля, можно рисовать красную линию прозрачными чернилами.

        Как видите, ничего невозможного. :)


      1. Wesha
        29.04.2024 16:30
        +7

        Семь перпендикулярных линий без проблем рисуются в неевклидовом пространстве.

        Кто скажет, что это простраство неевклидово, пусть первым бросит в меня котёнка!


      1. Metotron0
        29.04.2024 16:30

        Кроме того, в архитектуре термин "красная линия" имеет конкретное значение. Она вроде как обычно красная, но вдруг её цвет можно переопределить, оставив функциональность "красной линии"? Ведь на чёрно-белой печати она чёрная, но называется всё равно красной.


    1. RuVl
      29.04.2024 16:30

      Я именно за этим и пришёл :)


  1. HOMPAIN
    29.04.2024 16:30
    +9

    Так и не понял какую задачу вы упорно решали. Как нарисовать прямую линию? Как нарисовать симпатичную линию? Нарисовать линию постоянной толщины? Нарисовать линию максимально быстро? Или это урок по OpenGL?


    1. Porohovnik Автор
      29.04.2024 16:30
      +2

      Всё вместе. Как основная цель - показать, что рисовать линию толще одного пикселя/кегля не так просто как кажется на первый взгляд.

      OpenGl был выбран как удобный инструмент, для отображения и демонстрации формул, на урок по нему тут явно не хватает материала


      1. kinall
        29.04.2024 16:30
        +7

        показать, что рисовать линию толще одного пикселя/кегля не так просто как кажется на первый взгляд.

        Думаю, если бы вы начали статью с этих слов, было бы лучше. Я тоже долго скроллил в самый конец, но так и не понял, что хотел сказать автор.


  1. erzi
    29.04.2024 16:30
    +1

    Интересно, а почему в "закруглённом соединении" используется Безье, а не круг ? Ведь круг проще и быстрее


    1. Porohovnik Автор
      29.04.2024 16:30
      +2

      Кривую Безье можно задать одним уравнением с одной переменной, что удобно для объяснения и для построения.

      С кругом есть два варианта:

      1. Позиционирование и масштабирование круга так, чтобы он соединял две линии

      2. Рисовать два отдельных сегмента

      И в обоих способах придётся выполнить куда больше подготовительных вычеслений.

      Поэтому и был выбран способ с Безье - проще объяснить, и всего одно уравнение


      1. pvvv
        29.04.2024 16:30
        +4

        в пиксельном шейдере дорисовать круг в конце каждой линии, диаметром в её ширину


        1. UserSergeyB
          29.04.2024 16:30

          У вас не выйдет плавный переход из одной линии в другую. Чем толще будут линии, тем лучше сильнее это будет заметно.


          1. pvvv
            29.04.2024 16:30
            +1

            что будет заметно?

            или соединение двух линий разной толщины?

            так оно и с безье коряво (и даже хуже) будет выглядеть если толщины совсем разные.


  1. Porohovnik Автор
    29.04.2024 16:30

    del


  1. mynameco
    29.04.2024 16:30
    +8

    Проблемы начнутся, когда текстура на линии будет с альфой. Там где треугольники накладываюся, будет косяк.


    1. Veritaris
      29.04.2024 16:30

      Да, но это легко решается с помощью glAlphaFunc, glBlendEquation с параметрами GL_MIN либо GL_MAX, или же GL_DEPTH_TEST, что, конечно, намного проще, чем описанное в статье


      1. Porohovnik Автор
        29.04.2024 16:30

        Только есть нюанс - в статье об этой проблеме не слова.

        Ниже я кидал скрины полупрозрачной одноцветной линии - нет артефактов. Думаю что при наличии текстуры эффект будет тот же


  1. JerryI
    29.04.2024 16:30
    +3

    Спасибо за статью!

    Мне всегда было интересно, так как простые линии, и кривые Безье можно описать уравнениями. Можно ли тогда обойтись лишь фрагментным шейдером, который на вход берет их параметры и растеризует их попиксельно с помощью там sdf…?

    Да, и почему линий нет на аппаратном уровне?)


    1. Porohovnik Автор
      29.04.2024 16:30

      который на вход берет их параметры и растеризует их попиксельно с помощью там sdf…

      В принципе можно, но это сложнее и дольше. Потому и используется геометрический шейдер - он делает всю подготовительную работу для растеризатора

      Да, и почему линий нет на аппаратном уровне?)

      Есть примитив - линия, но она однопиксельная по ширине. Так как проблему стыка двух сегментов линий нельзя решить однозначно, то и универсального метода нет


      1. pvvv
        29.04.2024 16:30

        void glLineWidth(GLfloat width);

        правда его сломали уже


        1. Porohovnik Автор
          29.04.2024 16:30

          Да, он сломан.
          Во-первых макс ширина явно ограничена, это максимум для линии из статьи
          Во-вторых места соединение её сегментов не обработаны
          В-третьих нельзя использовать как материал ширину линий, а это значит нужно перед каждым вызовом запроса на рисование добавлять эту функцию.( когда через материал можно отрендерить сколько угодно линий за один запрос)(uniform я только для статьи в шейдерах на писал для простоты. По-факту там полноценный SSBO на всё и рендер через glMultiDrawArraysIndirect )


    1. RAX7
      29.04.2024 16:30
      +1

      На shadertoy можно найти подобные способы рисования кривых Безье

      https://www.shadertoy.com/results?query=bezier


  1. fiego
    29.04.2024 16:30
    +3

    А ещё бы теперь для микроконтроллеров с ручной растеризацией. Чётная толщина становится не таким простым вопросом... А потом вспомнить про антиалиасинг...


  1. chnav
    29.04.2024 16:30
    +4

    Смутило вот это:

    >> Метод на основе ортогональности

    >> Точки на плоскости можно представить как вектора, показывающие направления и скорость

    Какая скорость, при чём тут скорость ? Есть вектор, есть длина вектора aka модуль. Классическая терминология.

    И вообще решение задачи реализовано в стиле древней Греции, в лоб через системы трёхэтажных уравнений. Понятно что гимнастика для ума и пр. Второй способ это переизобретение векторной алгебры. Вместо матриц 2x2 используются раздутые развёрнутые формулы. Зачем - непонятно.


    1. Porohovnik Автор
      29.04.2024 16:30
      +4

      Какая скорость, при чём тут скорость ? Есть вектор, есть длина вектора aka модуль. Классическая терминология.

      Спасибо, поправил. Я в основном занимаюсь моделированием ветра, а там вектра это именно что скорость и направление. Вот и написал по-привычке.

      Понятно что гимнастика для ума и пр.

      Я старался ограничится геометрией для 9 класса и не вводить лишних абстракций.

      Вместо матриц 2x2 используются раздутые развёрнутые формулы

      С этим действием спорный момент. Можно было просто сказать что при повороте вектора на плоскости на 90 градусов, координаты полученного вектора зеркальны исходному, а в зависимости от поворота по часовой или против часовой стрелки меняется знак у x-координаты (как поступили вот здесь). Но я решил более подробно осветить этот момент, но оставаться хотелось в рамки школьной геометрии и алгебры. Поэтому я и не использовал понятие матрицы.


  1. MagisterAlexandr
    29.04.2024 16:30
    +1

    С одной линией всё нормально, попробуем нарисовать кривую из разных линий

    А что, если эта кривая должна быть полупрозрачной? :D


    1. Porohovnik Автор
      29.04.2024 16:30
      +1

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


      1. MagisterAlexandr
        29.04.2024 16:30

        Тогда уж стоит упомянуть Accumulation Buffer в OpenGL, и 3dfx T-Buffer.


  1. Classic_Fungus
    29.04.2024 16:30
    +1

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