Неявная поверхность

Неявная поверхность: Тор (R=40, a=15).

Неявная поверхность рода 2.

Неявная неалгебраическая поверхность (рюмка).

Неявная поверхность — это поверхность в евклидовом пространстве, определённая уравнением

$F(x,y,z)=0.$

Неявная поверхность является множеством нулей функции трёх переменных. Термин неявная здесь означает, что уравнение не решено относительно любой из переменных, x, y или z.

График функции обычно описывается уравнением $z=f(x,y)$ и такое представление называется явным. Третьим важным способом описания поверхности является параметрическое представление — $(x(s,t),y(s,t),z(s,t))$, где координаты x, y и z точек поверхности представлены тремя функциями $x(s,t)\,,y(s,t)\,,z(s,t)$, зависящими от общих параметров $s,t$. Обычно изменение представления поверхности осуществляется просто только в случае, если задано явное представление $z=f(x,y)$. Тогда другими двумя представлениями будут $z-f(x,y)=0$ (неявное) и $(s,t,f(s,t))$ (параметрическое).

Примеры:

1. плоскость $x+2y-3z+1=0.$
2. сфера $x^{2}+y^{2}+z^{2}-4=0.$
3. тор $(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+y^{2})=0.$
4. Поверхность рода 2: $2y(y^{2}-3x^{2})(1-z^{2})+(x^{2}+y^{2})^{2}-(9z^{2}-1)(1-z^{2})=0$ (см. рисунок).
5. Поверхность вращения $x^{2}+y^{2}-(\ln(z+3.2))^{2}-0{,}02=0$ (см. рисунок рюмка).

Для плоскости, сферы и тора имеется простое параметрическое представление, что неверно для четвёртого примера.

Теорема о неявной функции описывает условия, при которых уравнение $F(x,y,z)=0$ может быть решено (по меньшей мере неявно) относительно x, y или z. Но в общем случае явного решения может и не существовать. Эта теорема является ключевой для вычисления важных геометрических свойств поверхности, таких как касательные плоскости, нормали к поверхности, кривизны (см. ниже). Однако эти поверхности имеют существенный недостаток — их визуализация затруднительна.

Если $F(x,y,z)$ является многочленом от x, y и z, поверхность называется алгебраической. Пример 5 не является алгебраической поверхностью.

Несмотря на трудность визуализации неявные поверхности дают относительно простые техники для их теоретической генерации (например, Поверхность Штайнера) и поверхности, интересные для практических целей (см. ниже).

Формулы

При следующих соглашениях неявная поверхность представлена уравнением $F(x,y,z)=0$, где функция $F$ удовлетворяет необходимым условиям дифференцируемости. Ниже будем обозначать частные производные функции $F$ как $F_{x},F_{y},F_{z},F_{xx},\ldots$.

Касательная плоскость и нормальный вектор

Точка поверхности $(x_{0},y_{0},z_{0})$ называется регулярной тогда и только тогда, когда градиент функции $F$ в точке $(x_{0},y_{0},z_{0})$ не равен нулевому вектору $(0,0,0)$, что означает

$(F_{x}(x_{0},y_{0},z_{0}),F_{y}(x_{0},y_{0},z_{0}),F_{z}(x_{0},y_{0},z_{0}))\neq (0,0,0)$.

Если точка поверхности $(x_{0},y_{0},z_{0})$ не регулярна, она называется особой (употребляется также термин сингулярная точка).

Уравнение касательной плоскости в регулярной точке $(x_{0},y_{0},z_{0})$

$F_{x}(x_{0},y_{0},z_{0})(x-x_{0})+F_{y}(x_{0},y_{0},z_{0})(y-y_{0})+F_{z}(x_{0},y_{0},z_{0})(z-z_{0})=0,$

а уравнение нормального вектора

$\mathbf {n} (x_{0},y_{0},z_{0})=(F_{x}(x_{0},y_{0},z_{0}),F_{y}(x_{0},y_{0},z_{0}),F_{z}(x_{0},y_{0},z_{0}))^{T}.$

Нормальная кривизна

Чтобы сделать формулы проще, аргументы $(x_{0},y_{0},z_{0})$ в формуле ниже опущены. Тогда

$\kappa _{n}={\frac {\mathbf {v} ^{\top }H_{F}\mathbf {v} }{\|\operatorname {grad} F\|}}$

является нормальной кривизной поверхности в регулярной точке для единичного касательного вектора направления $\mathbf {v}$. $H_{F}$ является гессианом функции $F$ (матрицей вторых производных).

Доказательство этой формулы опирается (как и в случае неявной кривой) на теорему о неявной функции и формулу нормальной кривизны параметрической поверхности.

Приложения неявных поверхностей

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

Эквипотенциальная поверхность четырёх точечных зарядов

Эквипотенциальная поверхность двух точечных зарядов

Точечный заряд $q_{i}$ в точке $\mathbf {p} _{i}=(x_{i},y_{i},z_{i})$ образует потенциал в точке $\mathbf {p} =(x,y,z)$ (физические константы опущены)

$F_{i}(x,y,z)={\frac {q_{i}}{\|\mathbf {p} -\mathbf {p} _{i}\|}}.$

Эквипотенциальная поверхность для значения потенциала $c$ является неявной поверхностью $F_{i}(x,y,z)-c=0$, что есть сфера с центром в точке $\mathbf {p} _{i}$.

Потенциал четырёх точечных зарядов вычисляется пo формуле

$F(x,y,z)={\frac {q_{1}}{\|\mathbf {p} -\mathbf {p} _{1}\|}}+{\frac {q_{2}}{\|\mathbf {p} -\mathbf {p} _{2}\|}}+{\frac {q_{3}}{\|\mathbf {p} -\mathbf {p} _{3}\|}}+{\frac {q_{4}}{\|\mathbf {p} -\mathbf {p} _{4}\|}}.$

На рисунке четыре заряда имеют величину 1 и расположены в точках $(\pm 1,\pm 1,0)$. Показанная поверхность является эквипотенциальной поверхностью (неявной поверхностью) $F(x,y,z)-2{,}8=0$.

Поверхность постоянного произведения расстояний

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

На рисунке метаморфоз верхняя левая поверхность образуется по такому правилу. Эта поверхность является поверхностью уровня функции $F(x,y,z)-1{,}1=0$, где

${\begin{aligned}F(x,y,z)={}&{\Big (}{\sqrt {(x-1)^{2}+y^{2}+z^{2}}}\cdot {\sqrt {(x+1)^{2}+y^{2}+z^{2}}}\\&\qquad \cdot {\sqrt {x^{2}+(y-1)^{2}+z^{2}}}\cdot {\sqrt {x^{2}+(y+1)^{2}+z^{2}}}{\Big )}\end{aligned}}$

Метаморфоз двух неявных поверхностей (преобразование из одной поверхности в другую) — тора и поверхности постоянного произведения расстояний.

Метаморфозы неявных поверхностей

Другой простой метод создания новых неявных поверхностей называется метаморфозом неявных поверхностей:

Для двух неявных поверхностей $F_{1}(x,y,z)=0,F_{2}(x,y,z)=0$ (на рисунке это поверхность постоянного произведения расстояний и тор) определяются новые поверхности с использованием параметра $\mu \in [0,1]$:

$F(x,y,z)=\mu F_{1}(x,y,z)+(1-\mu )F_{2}(x,y,z)=0$

На рисунке представлены поверхности с величинами параметра $\mu =0,\,0{,}33,\,0{,}66,\,1$ .

Аппроксимация трёх торов (параллельная проекция)

Изображение POV-Ray (центральная проекция) приближения трёх торов.

Гладкая аппроксимация некоторых неявных поверхностей

$\Pi$-поверхности могут быть использованы для аппроксимирования любого гладкого и ограниченного объекта в $R^{3}$, поверхность которого определяется многочленом, который равен произведению других многочленов. Другими словами, мы можем создать любой гладкий объект с единой алгебраической поверхностью. Обозначим многочлены как $f_{i}\in \mathbb {R} [x_{1},\ldots ,x_{n}](i=1,\ldots ,k)$. Тогда аппроксимирующий объект определяется многочленом

$F(x,y,z)=\prod _{i}f_{i}(x,y,z)-r$

где $r\in \mathbb {R}$ определяет смешивающий параметр, контролирующий ошибку аппроксимации.

Аналогично гладкой аппроксимации неявных кривых уравнение

$F(x,y,z)=F_{1}(x,y,z)\cdot F_{2}(x,y,z)\cdot F_{3}(x,y,z)-r=0$

представляет для подходящих параметров $c$ гладкие аппроксимации трёх пересекающихся торов уравнениями

${\begin{aligned}F_{1}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+y^{2})=0,\\[3pt]F_{2}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+z^{2})=0,\\[3pt]F_{3}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(y^{2}+z^{2})=0.\end{aligned}}$

(На рисунке параметры равны $R=1,\,a=0{,}2,\,r=0{,}01.$)

POV-Ray изображение: Метаморфозы между сферой и поверхностью постоянного произведения расстояний (от 6 точек).

Визуализация неявных поверхностей

Есть несколько алгоритмов рендеринга неявных поверхностей, включая алгоритм «marching cubes». Фактически, имеются две идеи визуализации неявных поверхностей — одна создаёт сеть многоугольников, которые затем отрисовываются (см. Триангуляризация поверхности), а вторая опирается на трассировку лучей, когда определяются точки пересечения лучей с поверхностью.

См. также

• Неявная кривая

Литература

• Adriano N. Raposo, Abel J.P. Gomes. Pi-surfaces: products of implicit surfaces towards constructive composition of 3D objects // Journal of WSCG. — 2019. — arXiv:1906.06751. (International Conference in Central Europe on Computer Graphics, Visualization and Computer Vision)

• Jules Bloomenthal, Chandrajit Bajaj, Brian Wyvill. Introduction to Implicit Surfaces. — Morgan Kaufmann, 1997. — ISBN 978-1-55860-233-5.
• Ian Stephenson. Production Rendering: Design and Implementation. — Springer Science & Business Media, 2004. — ISBN 978-1-85233-821-3.
• Eric Haines, Tomas Akenine-Moller. Ray Tracing Gems. — Springer, 2019. — ISBN 978-1-4842-4427-2.
• Gomes A., Voiculescu I., Jorge J., Wyvill B., Galbraith C. Implicit Curves and Surfaces: Mathematics, Data Structures and Algorithms. — London: Springer-Verlag, 2009. — ISBN 978-1-84882-405-8.
• Thorpe J.A. Elementary Topics in Differential Geometry. — New York: Springer-Verlag, 1979. — (Undergraduate Texts in Mathematics). — ISBN 0-387-90357-7.
  • Торп Дж. Начальные главы дифференциальной геометрии. — «Мир», 1982. — (Современная математика. Вводные курсы).

Ссылки

• Sultanow: Implizite Flächen Архивная копия от 19 января 2021 на Wayback Machine
• Hartmann: Geometry and Algorithms for COMPUTER AIDED DESIGN Архивная копия от 30 октября 2017 на Wayback Machine
• GEOMVIEW Архивная копия от 11 июля 2000 на Wayback Machine
• K3Dsurf: 3d surface generator Архивная копия от 10 ноября 2020 на Wayback Machine
• SURF: Visualisierung algebraischer Flächen Архивная копия от 8 марта 2021 на Wayback Machine