Радиосити (компьютерная графика)

Сцена визуализирована с помощью RRV (простая реализация рендеринга радиосити на основе OpenGL) 79-я итерация

Корнеллский ящик , визуализированный с помощью BMRT и без него.

В 3D-компьютерной графике излучательность — это применение метода конечных элементов для решения уравнения рендеринга для сцен с поверхностями, которые диффузно отражают свет . В отличие от методов рендеринга , использующих алгоритмы Монте-Карло (такие как трассировка пути ), которые обрабатывают все типы световых путей, типичное излучение учитывает только пути (представленные кодом «LD*E»), которые покидают источник света и отражаются диффузно в некоторых случаях. количество раз (возможно, ноль), прежде чем попадет в глаз. Излучение — это глобального освещения алгоритм в том смысле, что освещение, попадающее на поверхность, исходит не только непосредственно от источников света, но и от других поверхностей, отражающих свет. Радиоситивность не зависит от точки зрения, что увеличивает объем необходимых вычислений, но делает их полезными для всех точек зрения.

Методы радиации были впервые разработаны примерно в 1950 году в области техники теплопередачи . Позже в 1984 году они были усовершенствованы специально для проблемы рендеринга компьютерной графики исследователями Корнелльского университета. и Университет Хиросимы .

Известными коммерческими движками радиосити являются Enlighten от Geomerics (используется в таких играх, как Battlefield 3 и Need for Speed: The Run ); 3дс Макс ; форма•Z ; LightWave 3D и электрическая система анимации изображений .

Визуальные характеристики

Разница между стандартным прямым освещением без теневой полутени и лучистостью с теневой полутенью

Включение вычислений излучательности в процесс рендеринга часто придает дополнительный элемент реализма готовой сцене, поскольку она имитирует явления реального мира. Рассмотрим простую сцену в комнате.

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

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

Обзор алгоритма радиосити

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

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

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

По мере выполнения алгоритма можно увидеть, как свет проникает в сцену, поскольку вычисляются множественные отражения. Отдельные пятна видны в виде квадратов на стенах и полу.

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

формулировка

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

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

Точнее, излучательность B — это энергия на единицу площади, покидающая поверхность пятна за дискретный интервал времени, и представляет собой комбинацию излучаемой и отраженной энергии:

$B(x)\,dA=E(x)\,dA+\rho (x)\,dA\int _{S}B(x'){\frac {1}{\pi r^{2}}}\cos \theta _{x}\cos \theta _{x'}\cdot \mathrm {Vis} (x,x')\,\mathrm {d} A'$

где:

• B(x) _i d A _i — полная энергия, оставляющая небольшую область d A _i вокруг точки x .
• E(x) _i d A _i — излучаемая энергия.
• ρ(x) — это отражательная способность точки, определяющая отраженную энергию на единицу площади путем умножения на падающую энергию на единицу площади (общая энергия, поступающая от других участков).
• S означает, что переменная интегрирования x ' проходит по всем поверхностям сцены.
• r - расстояние между x и x'
• θ _x и θ _x’ — это углы между линией, соединяющей x и x’ , и векторами, нормальными к поверхности в точках x и x’ соответственно.
• Vis( x , x' ) — это функция видимости, равная 1, если две точки x и x' видны друг от друга, и 0, если это не так.

Геометрический форм-фактор (или «проецируемый телесный угол») F _ij .

F _ij можно получить, проецируя элемент A _j на поверхность единичной полусферы, а затем проецируя его, в свою очередь, на единичный круг вокруг интересующей точки в плоскости A _i . Тогда форм-фактор равен доле единичного круга, охватываемой этой проекцией.

Формфакторы подчиняются соотношению взаимности A _i F _ij = A _j F _ji

Если поверхности аппроксимируются конечным числом плоских участков, каждый из которых имеет постоянную излучательность B _i и отражательную способность ρ _i , приведенное выше уравнение дает дискретное уравнение излучательности:

$B_{i}=E_{i}+\rho _{i}\sum _{j=1}^{n}F_{ij}B_{j}$

где Fij _{— коэффициент} геометрической видимости для излучения, выходящего из j и попадающего в пятно i .

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

Методы решения

Уравнение формально можно решить как матричное уравнение, чтобы получить векторное решение:

$B=(I-\rho F)^{-1}E\;$

Это дает полное решение «бесконечного отскока» для B напрямую. Однако количество вычислений для вычисления матричного решения масштабируется в соответствии с n ³ , где n — количество патчей. Это становится неприемлемым для реально больших значений n .

Вместо этого уравнение легче решить итеративно, многократно применяя приведенную выше формулу обновления однократного отскока. Формально это решение матричного уравнения итерацией Якоби . Поскольку коэффициенты отражения ρ _i меньше 1, эта схема быстро сходится, обычно требуя всего несколько итераций для получения разумного решения. Также можно использовать другие стандартные итерационные методы решения матричных уравнений, например метод Гаусса-Зейделя , где обновленные значения для каждого патча используются в расчете сразу после их расчета, а не все обновляются синхронно в конце каждого мести. Решение также можно настроить так, чтобы он поочередно перебирал каждый из отправляющих элементов в главном внешнем цикле для каждого обновления, а не для каждого из принимающих патчей. Это известен как для стрельбы вариант алгоритма , в отличие от варианта для сбора . Используя взаимность коэффициента обзора, A _i F _ij = A _j F _ji , уравнение обновления также можно переписать в терминах коэффициента просмотра. F _ji, видимый каждым отправляющим патчем A _j :

$A_{i}B_{i}=A_{i}E_{i}+\rho _{i}\sum _{j=1}^{n}A_{j}B_{j}F_{ji}$

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

F Сам коэффициент обзора ij _можно рассчитать несколькими способами. В ранних методах использовался полукуб (воображаемый куб с центром на первой поверхности, на который была проецирована вторая поверхность, изобретенный Майклом Ф. Коэном и Дональдом П. Гринбергом в 1985 году). Поверхность полукуба была разделена на пиксельные квадраты, для каждого из которых коэффициент обзора можно легко рассчитать аналитически. Полный форм-фактор затем можно было бы аппроксимировать путем сложения вклада каждого из пиксельных квадратов. Проекция на полукуб, которую можно было адаптировать из стандартных методов определения видимости многоугольников, также решила проблему промежуточных пятен, частично закрывающих те, что позади.

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

Подходы к выборке

Формфакторы F _ij сами по себе на самом деле явно не нужны ни в одном из уравнений обновления; ни для оценки общей интенсивности Σ _j F _ij B _j, мощность A _j B _j собранной из всего изображения, ни для оценки того, как распределяется излучаемая . Вместо этого эти обновления можно оценить с помощью выборочных методов без необходимости явного расчета форм-факторов. С середины 1990-х годов такие подходы к выборке стали наиболее часто используемыми методами практических расчетов радиационности.

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

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

Сокращение времени вычислений

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

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

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

Преимущества

Современный рендер чайника из Юты . В этой сцене для всего рассеянного освещения использовалось излучение.

Одним из преимуществ алгоритма Radiosity является то, что его относительно просто объяснить и реализовать. Это делает его полезным алгоритмом для обучения студентов алгоритмам глобального освещения. Типичный рендерер прямого освещения уже содержит почти все алгоритмы ( преобразования перспективы , наложение текстур , удаление скрытых поверхностей ), необходимые для реализации излучательности. Для понимания или реализации этого алгоритма не требуется глубокого понимания математики. .

Ограничения

Типичные методы измерения излучаемости учитывают только пути света формы LD*E, т.е. пути, которые начинаются от источника света и совершают множественные диффузные отражения, прежде чем достичь глаза. Хотя существует несколько подходов к интеграции других эффектов освещения, таких как зеркальное отражение. и глянцевый отражений, методы, основанные на излучении, обычно не используются для решения полного уравнения рендеринга.

Базовая излучательность также имеет проблемы с разрешением внезапных изменений видимости (например, резких теней), поскольку грубая, регулярная дискретизация на кусочно-постоянные элементы соответствует коробчатому фильтру нижних частот пространственной области. Построение сетки разрывов [1] использует знания о событиях видимости для создания более интеллектуальной дискретизации.

Путаница в терминологии

Radiosity был, пожалуй, первым широко используемым алгоритмом рендеринга, который учитывал рассеянное непрямое освещение. Более ранние алгоритмы рендеринга, такие как в стиле Уиттеда, трассировка лучей были способны вычислять такие эффекты, как отражения, преломления и тени, но, несмотря на то, что эти эффекты были весьма глобальными явлениями, эти эффекты обычно не назывались « глобальным освещением ». Как следствие, термины « диффузное взаимное отражение стали путать с «глобальным освещением» » и «излучение» в просторечии . Однако это три разные концепции.

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

См. также

• Корнелл Бокс
• Фотонное картирование
• Трассировка лучей
• Отслеживание пути

Дальнейшее чтение

• Обзор радиации из HyperGraph компании SIGGRAPH (предоставляет полный матричный алгоритм радиации и алгоритм прогрессивной радиации)
• Radiosity, Хьюго Элиас (также содержит общий обзор алгоритмов освещения и примеры программирования)
• Радиоситивность Аллена Мартина (немного более математическое объяснение радиоситивности)
• ROVER, доктор Тралвекс Йип (Библиотека рефератов и библиографии Radiosity)
• Радиосити: базовые реализации (Базовое исследование радиосити)

Внешние ссылки

• RADical от Парага Чаудхури (реализация варианта стрельбы и сортировки алгоритма прогрессивной радиации с ускорением OpenGL, являющегося продолжением GLUTRAD от Колбека)
• Radiosity Renderer и Visualizer (простая реализация Radiosity Renderer на основе OpenGL )
• Enlighten (Лицензионный программный код, обеспечивающий радиацию в реальном времени для компьютерных игровых приложений. Разработан британской компанией Geomerics )