CC BY
29
УДК 004.92 Дата подачи статьи: 06.05.20
DOI: 10.15827/0236-235X.132.697-705 2020. Т. 33. № 4. С. 697-705
Оптимизация расчета освещения интерьернъх сцен для. стохастической трассировки лучей
С.В. Ершов 1, к.ф.-м.н., старший научный сотрудник И.В. Валиев 1, научный сотрудник
А.Г. Волобой 1, д.ф.-м.н, ведущий научный сотрудник, [email protected]éldysh.ru
1 Институт прикладной математики им. М.В. Келдыша РАН, г. Москва, 125047, Россия
Компьютерное моделирование освещенности и построение реалистичных изображений повышают эффективность проектирования зданий. Наиболее распространенным методом, используемым для этого моделирования, является двунаправленная стохастическая трассировка лучей.
В статье рассматривается улучшение трассировки лучей методом Монте-Карло, которое оптимизирует излучение лучей от источника естественного света для интерьерных сцен. Предложенные алгоритмы позволяют ускорить сходимость метода, то есть уменьшить шум, остающийся после заданного времени трассировки, и таким образом улучшить итоговое изображение виртуальной сцены. Эти алгоритмы эффективны для сцен внутри помещений, освещаемых снаружи через окна или другие отверстия. Лучи от источников света генерируются так, что они направлены в эти окна, заведомо отсекая лучи, не освещающие интерьер. Другими словами, количество лучей увеличивается для направлений, которые вносят вклад в изображение, формируемое камерой.
Первый из предложенных алгоритмов требует дополнительных спецификаций пользователя для маркировки окон, прозрачных дверей или других отверстий, через которые свет может попасть во внутренние помещения. Второй алгоритм является полностью автоматическим. Он принадлежит к широкому классу алгоритмов нахождения и использования оптимальной PDF (Probability density function) в стохастической трассировке лучей, дает меньшее ускорение трассировки, но не требует действий от пользователя, в частности, исключая ошибку, которая может привести к недооценке освещенности помещения.
Приведено математическое обоснование корректности предлагаемых алгоритмов. На примерах показано, что предложенный подход позволяет в несколько раз быстрее рассчитать изображение с желаемым качеством.
Ключевые слова: компьютерная графика, реалистичная графика, моделирование освещенности, трассировка лучей, метод Монте-Карло, оптимальная функция PDF.
Моделирование освещенности и расчет глобального освещения широко применяются в настоящее время в архитектуре, дизайне новых устройств и новых материалов. Предварительные расчеты освещения и построение реалистичного изображения виртуальной модели повышают эффективность проектирования зданий. С их помощью можно заранее определить и добиться наилучшего воплощения целого ряда факторов, например, эстетичности разрабатываемых архитектором интерьеров помещений и внешнего вида. Важным аспектом строительства любого здания является экономичность (энергосбережение) его системы освещения. На основе точных расчетов освещения помещений дневным светом архитектор получает возможность спроектировать его максимальное использование, варьируя положение и ориентацию здания относительно сторон горизонта и окружающих строений, располо-
жение и размеры оконных проемов. При проектировании офисных помещений освещенности рабочих мест должны соответствовать стандартам, в которых указаны минимально допустимые величины освещенности, доля прямого света от источника в поле видимости человека и другие параметры. Использование предварительных расчетов в процессе проектирования возможно только при применении достаточно точных вычислительных методов.
Наиболее последовательным и физически точным методом, позволяющим учесть многие аспекты распространения света в различных средах и его взаимодействие с объектами, является стохастическая трассировка лучей [1]. Трассировка лучей методом Монте-Карло (MCRT - Monte Carlo Ray Tracing) генерирует или преобразует лучи случайным образом с распределением вероятностей, определенным локально, то есть, например, в точке попадания
луча в случае его рассеяния на поверхности или частицах, заданного с помощью двунаправленной функции распределения рассеяния (ДФР). Этот подход не является в конечном итоге оптимальным, например, ДФР может отправлять лучи в сторону от виртуального наблюдателя или камеры. Более того, в случае параллельного источника света возникает неоднозначность с точкой испускания луча: формально для параллельного света точка испускания луча должна выбираться равномерно в бесконечной области, что технически невозможно. Поэтому на практике точка испускания луча выбирается на проекции сцены на плоскость, ортогональную направлению луча, а излучение света определяется не только свойствами источника света, но и геометрией виртуальной сцены.
Если рассматривать построение реалистичного изображения внутренних интерьеров (например, офиса, построенное разработанной системой (см. http://www.swsys.ru/uploaded/ image/2020-4/2020-4-dop/2.jpg), то естественным продолжением такого подхода является дневное освещение через окна или открытые двери. В этом случае лучи, падающие на непрозрачные внешние стены, не усиливают внутреннее освещение, которое создается только лучами, проходящими через окна. Следовательно, необходимо выбирать источник луча в проекции окон на плоскость, перпендикулярную направлению освещения. Это основная идея оптимизации излучения света для интерь-ерных сцен. Однако должны быть правильно промоделированы и более сложные случаи: во-первых, когда источник дневного освещения задается не параллельным источником света, а небесной полусферой с распределением яркости; во-вторых, в случае с несколькими окнами, которые в проекции могут перекрываться.
При задании виртуальной сцены пользователю предоставляется возможность указать окна, открытые двери и люки. Затем распределение происхождения лучей маскируется 1 внутри или 0 вне их проекции. Однако такой подход может привести к ошибке, например, если в сцене дом окружен снегом, то есть сильно отражающей поверхностью. Солнечные лучи, отраженные снегом, значительно влияют на освещение интерьера, но они будут находиться вне проекции, так как свет излучается только к прозрачным окнам.
Эту ситуацию можно исправить, если сделать маску излучения более сложной, разрешив
задавать ненулевые значения для некоторых непрозрачных областей, и, наоборот, давать значения меньше 1 для некоторых прозрачных областей, если они не открываются вовнутрь. В результате получаем некоторую функцию плотности вероятности источника излучения (PDF - Probability density function). Эта функция не может быть просто постулирована, а должна быть рассчитана на основе результатов моделирования. В этом случае данный метод принадлежит к широкому классу алгоритмов нахождения и использования оптимальной PDF в стохастической трассировке лучей.
Величины, представляющие интерес для оптического моделирования (например, яркость точки пересечения), как правило, являются интегралами произведения, например, точечная яркость - это свертка ДФР с освещением точки. Поэтому оптимальным распределением вероятности направления рассеянного луча в MCRT должно быть это произведение. Проблема в том, что один множитель этого произведения (в приведенном примере - освещение) сам является результатом стохастической трассировки и не известен заранее. Уже давно предпринимаются попытки повысить эффективность Монте-Карло трассировки лучей, делая PDF рассеянных лучей как можно ближе к этому произведению, используя различные приближения и эвристические методы. Скажем, яркость поверхности даже при прямом освещении распределеным источником света является интегралом по этой освещенной поверхности. Для повышения эффективности его расчета авторы работы [2] предложили выделить наиболее быстро меняющиеся члены подынтегрального выражения и использовать их для построения PDF. Исследования [3, 4] посвящены повышению эффективности обратной MCRT за счет рассеяния луча камеры не только по ДФР, но и по ее произведению с освещенностью точки. В [3] фотонные карты использовались для оценки углового распределения освещенности. В [4] это было сделано с 5D^-ревом, покрывающим сцену, которая накапливала освещение, приносимое различными отраженными обратными лучами. Таким образом, удалось избежать отдельной фазы прямой трассировки.
С разработанным методом связаны две основные проблемы. Во-первых, оптимальное распределение теперь должно быть сохранено не в каждой точке попадания, а в каждом вок-селе, что требует много памяти. Для его уменьшения используются различные приближения,
например, довольно грубое пространственное разрешение. И это приводит ко второй проблеме: только довольно грубая оценка распределения входящего излучения. Таким образом, он отклоняется от целевого произведения, и разница должна быть компенсирована масштабом энергии луча. Наконец, можно оказаться в тупике с помощью методов, подобных [4], где освещение берется из тех же лучей трассировки, которые были испущены с помощью накопленной PDF. Например, если через какое-то время в каком-либо угловом конусе не было испущено ни одного луча, предварительно оцененное освещение от него равно 0. А при генерации следующего луча процедура отбора также не отправит луч в эту область, поэтому освещение от него остается равным 0 и далее. Самое простое средство борьбы с такой ошибкой - использовать смешанную стратегию, когда ДФР рассеивает некоторую долю лучей без использования накопленной PDF. Сложность заключается в том, как найти соотношение между лучами, которые необходимо генерировать с учетом PDF, и лучами без этого учета. Если оно слишком мало, то узкие пики освещения не повлияют на PDF и берутся только из небольшой доли лучей, рассеянных ДФР, которые создают сильный шум. Если отношение велико, это убивает само преимущество (ускорение) метода.
Эта основная идея была использована многими исследователями [5-9]. Если в [6] используется полусферическая PDF в форме, довольно похожей на работу Дженсена [3], то в [7] авторы используют приближение по сумме гауссианов. Его преимущество в том, что произведение двух гауссовских колоколов также является гауссовым и может эффективно отбираться. В [9] вместо этого используются вейвлеты. Некоторые авторы используют сферические гармоники, произведение которых также допускает аналитическую обработку. Используются и другие подходы, вплоть до нейронных сетей и машинного обучения, которые являются еще одним способом оптимизации PDF направления луча на основе накопленной статистики [8].
Все вышеперечисленные методы меняют PDF в одной точке, чтобы уменьшить дисперсию. Радикальным решением является метод Metropolis Light Transport [10], где, в отличие от MCRT, которая отслеживает луч сегмент за сегментом, случайным образом выбирается весь путь от источника света до камеры. Хотя проблема трассировки лучей устранена, появ-
ляются другие проблемы, анализ которых выходит за рамки данной статьи.
Менее радикальное решение по-прежнему использует MCRT, но применяет разные подходы (как для генерации направления рассеяния, так и для сбора освещенности) и генерирует сразу несколько траекторий лучей. Затем яркость пикселя берется как взвешенная сумма яркости, которую принесла каждая из траекторий. Обычно это называется стратегией множественной важности (multiple importance strategy). Этот метод не оптимизирует PDF в точках попадания отдельно, но вместо этого он смешивает вклады нескольких различных целых путей лучей. Основные идеи были предложены в [11], где, помимо всего прочего, доказано, что относительно простая эвристическая оценка баланса близка к оптимальной, то есть другие оценки рассматриваемого класса не могут существенно уменьшить дисперсию. Однако это не всегда так. В исследовании [12] авторы проанализировали условия в теореме Вича [11] и обнаружили, что иногда оценка может быть построена значительно лучше, чем эвристический баланс.
Освещение сцены через окна
Генерация лучей от источника света.
Для параллельного источника света точки испускания луча распределены равномерно. С другой стороны, лучи, не попадающие в окно, бесполезны, так как не создают внутреннего освещения. Поэтому начало луча выбирается равномерно в проекции окна на плоскость, ортогональную направлению света.
В случае триангулированной геометрии сцены используются треугольники, которые составляют оконное стекло. Сначала выбираем один из них с вероятностью, пропорциональной потоку энергии через этот треугольник, а затем равномерно выбираем начало луча внутри треугольника. Возможно, что окно открыто или не имеет стекла внутри. В этом случае внутри окна нет треугольников для направления лучей. Тогда пользователь может выбрать треугольники, образующие раму окна. После этого из них строится выпуклая оболочка и лучи направляются на ее треугольники.
В случае непараллельного источника света расчет немного сложнее. Скажем, свет, образованный небесной полусферой, бесконечно далек, но у него в основном гладкое распределение яркости (гониограмма). Здесь сначала выбирается направление излучения v, а затем -
(1)
треугольник t с вероятностью 5,|(у ■ и, )| ' I | (V ■ и, )|, где & - площадь треугольника; t и щ - его нормаль. Начало луча х выбирается равномерно в этом треугольнике и переводится за пределы области сцены. Направление V выбирается с плотностью вероятности
1(V ■ и, )|* (V)
p(v) =
X S,J|(v• n,)|g(v)d2v:
где g - интенсивность гониограммы интенсивности.
Предположим, что небесная гониограмма задается таблицей на прямоугольной сетке (0, ф) и билинейно интерполируется внутри ячеек. Тогда сначала выбирается ячейка гониограммы [$, ^^[фу, ф;+1] с вероятностью
1Г |(v ■ и, )|* ^ 2v
Р = ' "-, (2)
'" IА Г | (V ■ и, )|£ ^ \ а затем в этой ячейке выбирается направление с плотностью вероятности
1| (V■ п,)|8(V)
Р1 ! (v) » 1 . г,
1Г^ ■ и ^О*^ (3)
Выбираем треугольник t с вероятностью (1). Если ячейка мала так, что направления излучения для всех ее четырех вершин находятся в одной и той же стороне от нормали треугольника, интеграл по ячейке равен
Г ■ и,)|8(у)с12v = |(, ■ и,)|,
где Fj = j vg(v)d2v
зависит только от ячейки
гониограммы, но не зависит от треугольника. Поэтому для ячейки гониограммы достаточно вычислить три значения - ^ , ^(]), ^^), и тогда для большинства треугольников поток через них вычисляется как скалярное произведение точек этого вектора с нормалью треугольника. Для тех немногих треугольников, которые освещаются с разных сторон ячейкой гониограммы, поток энергии необходимо рассчитать напрямую.
Следовательно, в вышеупомянутом методе плотность вероятности начала и направления луча (х, V) равна
1Л(Хv) 8(v)
Р( *v) =
jS,X,( *v) g (v)d 2 vd 2 * X,X,( *v) g (v)
X, S, j\(v • n )|g (v)d 2v'
где уа(х, V) = 1, если луч (х, V) пересекает треугольник t, и 0 в противном случае.
Перекрытие проекции окон. Предложенный метод создает правильный поток энергии через поверхность, состоящую из треугольников (или других граней), проекции которых не перекрываются. Если проекции N треугольников пересекаются (рис. 1), числитель в (4) увеличивается для этой области, что вызывает искажение.
Рис. 1. Сцена с несколькими окнами (лучи, показанные толстыми стрелками, пересекают два окна)
Fig. 1. The scene with multiple windows (the rays shown by bold arrows cross two windows)
Очевидное исправление состоит в том, чтобы изменить (4) на выражение
IN "Ч х, ^х, ( Х, V) 8 (V)
Р( * v) =
jX N "Ч *, v)x, ( *, v) g (v)d 2vd2 *
ч~> ч~> л ......., (5)
где N(x, v) = YtXt(x, v) - число треугольников, пересекающих луч (x, v).
С этим корректирующим фактором легко справиться с помощью метода отклонения: сначала выбираются направление и начало луча в соответствии с плотностью вероятности (4), то есть фактически с (1), (2) и (3), затем перебираются все треугольники всех окон и проверяется, пересекает ли их луч. Это делается для всех изолированных треугольников, то есть не имеет значения наличие других треугольников. Затем с вероятностью 1 - N_1(x, v) луч отклоняется, то есть заново выбираются направление луча и начало координат. Это продолжается до тех пор, пока луч не будет принят. Тогда результирующая плотность будет (4), и она обеспечивает правильное освещение.
Расчет оптимальной адаптивной функции PDF
Как видно из (5), вышеупомянутый подход маскирует собственную пространственно-однородную плотность вероятности источника
света маской, которая равна 1, если луч пересекает любой оконный треугольник, и 0 в противном случае. Это простейшее изменение не является оптимальным. Рассмотрим более гибкий подход, когда происхождение и направление испускаемого луча выбираются с произвольной плотностью р(х, V).
После выбора луча источника света V) его дальнейшее распространение определяется случайным образом в каждой точке пересечения. Соответственно, его вклад в яркость пикселя p является случайной функцией Lp(x, V). Естественно, если изменяется распределение лучей, это должно быть компенсировано изменением их вклада, чтобы сохранить математическое ожидание. Следовательно,
Т I \ И I \ Р(x,у) Ь (х, V) =Лр (х, V)—-- ,
р( х, V)
где Лр(г, V) - вклад для правильной плотности
(5).
Предполагая, что траектория луча камеры является детерминированной, шум в p полностью обусловлен прямой Монте-Карло трассировкой лучей и может быть рассчитан как
|( Ь >( х, v)р( х, v)d2 хй2 V -
- (|(Ьр>(х, v)р(х,v)d2хй2v)2,
где (•> — усреднение по всем случайным путям луча, испускаемого источником света с началом x и направлением V. Общая ошибка — это просто сумма по всем пикселям. Вариация плотности вероятности при помощи 8р изменяет эту суммарную ошибку на выражение
-}^(ЛР Хх^)4^р(х^)й2хй2^ (6)
р р (х, v)
Оптимальный р(x, V) минимизирует шум, и в этом экстремуме изменение должно уменьшаться до любого допустимого 8р. Вариация нормализованной плотности, по определению, есть произвольная функция, для которой |8р( х, v)d2 хй2 v = 0. Сравнивая с (6), делаем вывод, что она должна быть (общий вывод в методах Монте-Карло) следующей:
£<Л2 >(X,V) Р^ = const Р Р (x,V)
или р(X, V) = Cp(X, v) /У<Лp >(X, V)
(7)
где — нормировочная константа.
Значение (Л 2 > (х, v), по определению, не зависит от того, как испускаются лучи (к, V). Перед трассировкой лучей создается четырехмер-
ная сетка в пространстве V). Затем, когда испускается луч, находится ячейка, к которой принадлежит V), и добавляется (Л 2 >( х, v) к
значению, накопленному в этой ячейке. По мере трассировки лучей идет накопление все лучшей и лучшей оценки (Л2 >(х, v) в каждой
ячейке и начинается использование ее для оптимизации плотности. Начало и направление луча теперь выбираются с плотностью (7). По мере продолжения трассировки лучей и изменения накопленной оценки (Л2 >(х, v) изменяется и плотность излучения. Другими словами, теперь все зависит от времени. Это, однако, не создает петли обратной связи, потому что накопленное значение не зависит от выбора начала луча. Постепенно, когда накопленная оценка (Л2 >(х, v) улучшается, плотность позиции начала луча сходится к своей предельной оптимальной форме.
Результаты
Преимущество предложенных методов продемонстрируем на двух примерах.
Первая сцена задает внутреннюю комнату без окон, свет в которую может попасть только через щель под дверью. Дверь связывает ее с внешней комнатой, которая освещена дневным светом через дверь и большое окно (рис. 2). Камера находится во внутренней комнате и смотрит правее двери.
На рисунке 3 представлены изображения для сцены «Внутренняя комната», рассчитанные в течение одного и того же времени. Средняя освещенность и среднеквадратичные значения по выделенному прямоугольнику, показанному на рисунке 3, приведены в таблице 1.
Надо заметить, что выбор щели под дверью в качестве источника освещения, хотя и позволяет получить наименее зашумленную картинку, является неправильным. Эта щель находится в интерьере и не освещена снаружи по всей полусфере, как, например, окно в правой
щ
а) 6)
: Ш
в) г) Рис. 3. Сцена «Внутренняя комната»: а) освещение на всю область сцены без оптимизации, б) освещение сцены через окно и дверь внешней комнаты, в) освещение сцены через щель под дверью, г) освещение сцены, построенное с помощью адаптивной PDF Fig. 3. Scene "Inner room": а) lighting throughout the scene without optimization, б) lighting through the window and door of the external room, в) lighting through the gap under the door, г) lighting is constructed using adaptive PDF
комнате, а вторичный свет от стен внешней комнаты также потерян, значит, мы явно занижаем освещенность. В результате средняя освещенность в выбранном прямоугольнике в 5,7 раза меньше, чем при правильном выборе окна. Также сама щель выглядит полностью черной, поскольку лучи попадают только в щель, а пол и вся внешняя комната не освещены. Для трех остальных вариантов можно видеть, что наименьший шум наблюдается во втором методе, при освещении через дверь и окно, выбранные пользователем. Средняя освещенность во всех трех вариантах одинаковая, в пределах точности вычислений.
Таблица 1 Значения освещенности (сцена «Внутренняя комната») (кд/м2)
Table 1
The average illumination (Scene "Inner room") (cd/m2)
Другой пример сцены показан на рисунке 4. Дом освещается дневным светом. Здесь много окон, через которые он может проникать внутрь. Поэтому второй метод (явное указание окон в сцене) неудобен.
Рис. 4. Внешний вид сцены «Дом» Fig. 4. The appearance of the scene "House"
На рисунке 5 приведены изображения одной из комнат на втором этаже, рассчитанные в течение одного и того же временного интервала, без оптимизации излучения света (рис. 5 а) и при использовании адаптивной функции PDF (рис. 5б). Значения средней освещенности и среднеквадратичной ошибки в области, обозначенной выделенным прямоугольником, приведены в таблице 2. Видно, что уровень шума заметно ниже (в 4,25 раза) при использовании адаптивной функции PDF. Визуально такое же качество изображения может было получено примерно в 18 раз быстрее при использовании функции PDF.
Таблица 2
Значения освещенности (внутренняя комната для сцены «Дом») (кд/м2)
Table 2
The average illumination (room for the scene "House") (cd/m2)
Освещение Среднее Среднеквадратичная ошибка
Без оптимизации 266 1 425
Адаптивное 203 335
Заключение
Разработанные авторами методы позволяют ускорить стохастическую трассировку лучей для интерьерных сцен, где интересующая область освещается снаружи через окна или
Освещение Среднее Среднеквадратичное
Без оптимизации 57,1 30,3
Через дверь и окно 57,7 11,7
Через щель 10,0 1,36
Адаптивное 57,1 18,2
а) б)
Рис. 5. Изображение внутренней комнаты для сцены «Дом»: а) излучение света на всю область сцены, б) освещение, оптимизированное адаптивной функцией PDF
Fig. 5. The interior room for the scene "House": а) light emission over the entire area of the scene,
6) lighting is optimized by adaptive PDF
небольшие отверстия. Продемонстрирована математическая корректность предложенных методов. Оба метода снижают уровень шума или, что то же самое, позволяют быстрее достичь своего целевого уровня. Первый метод требует явной маркировки этих окон в геометрии сцены, что создает дополнительную нагрузку на пользователя, однако он дает лучшее шумоподавление. В случе неправильного выбора окна он может недооценивать освещенность, как это видно из первого примера, приведенного выше. Второй метод создает оптимальную функцию PDF излучения света, использование которой минимизирует шум. Это медленнее, но дает правильное освещение без необходимости явно маркировать окна и возможности допустить в этом ошибку. Тем не ме-
нее, второй метод ускоряет генерацию изображения в разы по сравнению с излучением света без оптимизации для всей области.
Оба метода были реализованы и встроены в комплекс моделирования освещенности и реалистичной визуализации Lumicept [13], разработанный в ИПМ им. М.В. Келдыша РАН. По умолчанию для интерьерных сцен всегда строится функция оптимального освещения PDF, которая позволяет значительно ускорить расчеты независимо от сложности виртуальной цены. Однако в простых случаях пользователю дается возможность маркировать окна, дыры или другие прозрачные поверхности. В таком случае расчеты еще ускоряются, но ответственность за правильность маркировки окон ложится на пользователя.
Литература
1. Pharr M., Humphreys G., Wenzel J. Physically Based Rendering, From Theory to Implementation. USA, SF, Morgan Kaufmann Publ., 2016, 1266 p.
2. Shirley P., Wang C., Zimmerman K. Monte Carlo techniques for direct lighting calculations. ACM TOG, 1996, vol. 15, no. 1, pp. 1-36. DOI: 10.1145/226150.226151.
3. Jensen H.W. Importance driven path tracing using the photon map. Proc. Rendering Techniques, 1995, pp. 326-335. DOI: 10.1007/978-3-7091-9430-0_31.
4. Lafortune E.P., Willems Y.D. A 5D tree to reduce the variance of Monte Carlo ray tracing. Proc. Rendering Techniques, 1995, pp. 11-20. DOI: 10.1007/978-3-7091-9430-0_2.
5. Burke D., Ghosh A., Heidrich W. Bidirectional importance sampling for direct illumination. EGSR, 2005, pp. 147-156. DOI: 10.2312/EGWR/EGSR05/147-156.
6. Hey H., Purgathofer W. Importance sampling with hemispherical particle footprints. Proc. 18th SCCG, 2002, pp. 107-114. DOI: 10.1145/584458.584476.
7. Herholz S., Elek O., Vorba J., Lensch H., Krivanek J. Product importance sampling for light transport path guiding. Comp. Graph. Forum, 2016, vol. 35, no. 4, pp. 67-77. DOI: 10.1111/cgf.12950.
8. Vorba J., Karlik O., Sik M., Ritschel T., Krivanek J. On-line learning of parametric mixture models for light transport simulation. ACM Trans. Graph., 2014, vol. 33, no. 4, art. 101, pp. 1-11. DOI: 10.1145/2601097. 2601203.
9. Clarberg P., Jarosz W., Akenine-Möller W., Jensen H.W. Wavelet importance sampling: efficiently evaluating products of complex functions. ACM Trans. Graph., 2005, vol. 24, no. 3, pp. 1166-1175. DOI: 10.1145/1073204.1073328.
10. Hachisuka T., Kaplanyan A.S., Dashsbacher C. Multiplexed metropolis light transport. ACM Trans. Graph., 2014, vol. 33, no. 4, art. 100, pp. 1-10. DOI: 10.1145/2601097.2601138.
11. Veach E., Guibas L.J. Optimally combining sampling techniques for Monte Carlo rendering. Proc. SIGGRAPH, 1995, pp. 419-428. DOI: 10.1145/218380.218498.
12. Sbert M., Havran V., Szirmay-Kalos L. Multiple importance sampling revisited: breaking the bounds. EURASIP J. Adv. Signal Processing, 2018, vol. 2018, art. 15. DOI: 10.1186/s13634-018-0531-2.
13. Жданов Д.Д., Потемин И.С., Галактионов В.А., Барладян Б.Х., Востряков К.А., Шапиро Л.З. Спектральная трассировка лучей в задачах построения фотореалистичных изображений // Программирование. 2011. № 5. C. 13-26.
Software & Systems Received 06.05.20
DOI: 10.15827/0236-235X. 132.697-705 2020, vol. 33, no. 4, pp. 697-705
Optimization of lighting calculation for interior scenes for stochastic ray tracing
S.V. Ershov 1 Ph.D. (Physics and Mathematics), Senior Researcher I. V. Valiev 1 Researcher
A.G. Voloboy 1 Dr.Sc. (Physics and Mathematics), Leading Researcher, [email protected] 1 Keldysh Institute of Applied Mathematics of RAS, Moscow, 125047, Russian Federation
Abstract. Lighting simulation and the creation of realistic images increases the efficiency of building design. The most common method used for this simulation is bidirectional Monte Carlo ray tracing.
The paper discusses the improvement of the Monte Carlo ray tracing which optimizes the emitting of rays from a natural light source for interior scenes. The proposed algorithms can accelerate the convergence of the method, i.e. reduce the noise remaining after a certain simulation time, and thus improve the final image of the virtual scene. These algorithms are effective for indoor scenes illuminated externally through windows or other openings. Rays from light sources are generated so that they are directed into these windows, deliberately cutting off the rays that do not illuminate the interior. In other words, the number of rays increases for directions that contribute to the image formed by the camera.
The first of the proposed algorithms require supplementary user specifications for marking windows, transparent doors, or other openings through which light can enter into interior rooms. The second algorithm is fully automatic. It belongs to a wide class of algorithms for finding and using the optimal PDF (Probability density function) in stochastic ray tracing. The second algorithm provides less acceleration for simulation but does not require action from the user, in particular, eliminating the potential error that can lead to underestimation of the room illumination.
The paper provides the mathematical justification of the correctness of the proposed algorithms. The examples show that the proposed approach allows calculating the image with the desired quality several times faster.
Keywords: computer graphics, realistic graphics, lighting simulation, ray tracing, Monte Carlo method, optimal PDF.
References
1. Pharr M., Humphreys G., Wenzel J. Physically Based Rendering, from Theory to Implementation. USA, SF, Morgan Kaufmann Publ., 2016, 1266 p.
2. Shirley P., Wang C., Zimmerman K. Monte Carlo techniques for direct lighting calculations. ACM TOG, 1996, vol. 15, no. 1, pp. 1-36. DOI: 10.1145/226150.226151.
3. Jensen H.W. Importance driven path tracing using the photon map. Proc. Rendering Techniques, 1995, pp. 326-335. DOI: 10.1007/978-3-7091-9430-0_31
4. Lafortune E.P., Willems Y.D. A 5D Tree to reduce the variance of Monte Carlo ray tracing. Proc. Rendering Techniques, 1995, pp. 11-20. DOI: 10.1007/978-3-7091-9430-0_2.
5. Burke D., Ghosh A., Heidrich W. Bidirectional importance sampling for direct illumination. EGSR, 2005, pp. 147-156. DOI: 10.2312/EGWR/EGSR05/147-156.
6. Hey H., Purgathofer W. Importance sampling with hemispherical particle footprints. Proc. 18th SCCG, 2002, pp. 107-114. DOI: 10.1145/584458.584476.
7. Herholz S., Elek O., Vorba J., Lensch H., Krivanek J. Product importance sampling for light transport path guiding. Comp. Graph. Forum, 2016, vol. 35, no. 4, pp. 67-77. DOI: 10.1111/cgf.12950.
8. Vorba J., Karlik O., Sik M., Ritschel T., Krivanek J. On-line learning of parametric mixture models for light transport simulation. ACM Trans. Graph., 2014, vol. 33, no. 4, art. 101, pp. 1-11. DOI: 10.1145/2601097. 2601203.
9. Clarberg P., Jarosz W., Akenine-Moller W., Jensen H.W. Wavelet importance sampling: efficiently evaluating products of complex functions. ACM Trans. Graph., 2005, vol. 24, no. 3, pp. 1166-1175. DOI: 10.1145/1073204.1073328.
10. Hachisuka T., Kaplanyan A.S., Dashsbacher C. Multiplexed metropolis light transport. ACM Trans. Graph., 2014, vol. 33, no. 4, art. 100, pp. 1-10. DOI: 10.1145/2601097.2601138.
11. Veach E., Guibas L.J. Optimally combining sampling techniques for Monte Carlo rendering. Proc. SIGGRAPH, 1995, pp. 419-428. DOI: 10.1145/218380.218498.
12. Sbert M., Havran V., Szirmay-Kalos L. Multiple importance sampling revisited: breaking the bounds. EURASIP J. Adv. Signal Processing, 2018, vol. 2018, art. 15. DOI: 10.1186/s13634-018-0531-2.
13. Zhdanov D.D., Potemin I.S., Galaktionov V.A., Barladyan B.Kh., Vostryakov K.A., Shapiro L.Z. Spectral ray tracing in problems of photorealistic imagery construction. Progr. and Comp. Soft., 2011, vol. 37, no. 5, pp. 236-244. DOI: 10.1134/S0361768811050069.
Для цитирования
Ершов С.В., Валиев И.В., Волобой А.Г. Оптимизация расчета освещения интерьерных сцен для стохастической трассировки лучей // Программные продукты и системы. 2020. Т. 33. № 4. С. 697-705. DOI: 10.15827/0236-235X.132.697-705.
For citation
Ershov S.V., Valiev I.V., Voloboy A.G. Optimization of lighting calculation for interior scenes for stochastic ray tracing. Software & Systems, 2020, vol. 33, no. 4, pp. 697-705 (in Russ.). DOI: 10.15827/0236-235X. 132.697-705.
