Moskva, Moscow, Russian Federation
The article deals with modeling the spatial-angular distribution of luminance (SADL) in lighting scenes based on the solution of the global lighting equation by local estimates of the Monte Carlo method. Based on the SADL, a new criterion for lighting quality is formulated, which, in contrast to the generally accepted UGR to date, does not take into account the point glare sources, but directly the continuous distribution of luminance. The proposed criterion is applied in the analysis of lighting at some stations of the Moscow metro, which shows its correlation with the feelings of real observers.
spatial-angular distribution of luminance, quality of illumination, global illumination
1. Введение
Одним из ключевых направлений светотехники является проектирование осветительных установок (ОУ). При проектировании ОУ определяются типы световых приборов, их количество, расположение, направление и т.п. При этом проектировщик руководствуется нормативными документами, определяющими качественные и количественные показатели освещения.
В современных нормативных документах для неспециальных осветительных установок (офисные, производственные, торговые и прочие) в качестве количественной характеристики нормируется освещенность и различные параметры, производные от нее (отношение минимальной освещенности к максимальной и т.п.). Как правило, все расчеты ведутся для освещенности на полу помещения или на мнимой рабочей плоскости, расположенной на высоте стола. Однако освещенность – это интегральная характеристика падающего света, тогда как глаз человека реагирует на отраженный от поверхности свет. То есть, если взять абсолютно черную поверхность с коэффициентом отражения равным нулю, то формально на ней можно получить требуемую освещенность, при этом визуально мы ничего не увидим, так как отражаться от поверхности ничего не будет. С точки зрения органа зрения человека нормировать необходимо яркость. Сложившаяся ситуация вполне объяснима, еще совсем недавно расчет и измерения яркости представляли собой крайне сложную задачу.
Частично устраняет эту проблему введение среди множества качественных показателей для неспециальных установок объединенного показателя дискомфорта
где Li – яркость блеского источника, кд/м2, 𝜔i – угловой размер блеского источника, стер, pi– индекс позиции блеского источника относительно линии зрения, Lа – яркость адаптации, кд/м2.
Таким образом, UGR позволил выразить всего одним числом качество освещения и вошел в нормативные документы. Таким образом, сегодня проектировщик при проектировании не специальных ОУ, таких как торговые, офисные, общественные, производственные помещения и многие другие, руководствуется освещенностью в качестве количественной характеристики и UGR как качественной оценкой освещения. Где UGR отвечает на вопрос – насколько комфортно человеку будет находиться в пределах осветительной установки.
Однако формула справедлива лишь для малоугловых равномерных блеских источников. Т.е. протяженные неравномерные блики она учесть не может в принципе. Дело осложняется еще и тем, что широко используемые в светотехническом проектировании программы моделирования DIALux и Relux построены на основе метода конечных элементов и решают не уравнение глобального освещения относительно яркости, а уравнение излучательности в диффузном приближении. Что очевидно приводит к тому, что вторичные блики не могут быть учтены в принципе. Существенным шагом вперед в светотехническом проектировании становится внедрение DIALuxEvo основанного на методе фотонных карт. Однако методика расчета UGR в нем не претерпела изменений.
Таким образом, ОУ проектируются, лишь отдаленно оценивая, насколько комфортно будет человеку в ней и нормируют невидимую характеристику – освещенность. Однако еще в начале прошлого века было высказано предположение, что ключевую роль в вопросе комфорта играет пространственно-угловое распределение яркости [1].
2. Пространственно-угловое распределение яркости
В основе моделирования осветительных установок лежит хорошо известное в компьютерной графике уравнение глобального освещения, впервые полученные James T Kajiya в 1986 году [2]
,
где – яркость в точке r по направлению
,
– двунаправленная функция отражения (отражения или пропускания), L0 – прямая компонента яркости, непосредственно от источников,
– нормаль в точке r к элементу поверхности сцены.
Уравнение записано относительно точки r, находящейся на поверхности сцены, однако в задаче оценки качества освещения, наблюдатель находится в объеме сцены. После ряда преобразований, можно получить уравнение, записанное относительно точки в объеме сцены [3]
где – ядро уравнения.
Уравнение (3) описывает пространственно-угловое распределение яркости (ПУРЯ) в каждой точке пространства сцены. Что позволяет подойти к вопросу оценки качества освещения не на основе оценки отдельных бликов как в UGR, а на основе анализа непрерывного пространственно-углового распределения яркости.
Рис. 1. Положения точек в пространстве сцены: слева – изображена точка на поверхности сцены, на рисунке справа – точка в пространстве сцены, а которых требуется определить яркость в заданном направлении
В 1997 году была опубликована работа «Instant Radiosity» [4]. Эта публикация положила начало новому подходу в решении уравнения ГО, которое в конечном счете в компьютерной графике не получило широкого распространения. Автор сформулировал алгоритм нового метода в феноменологическом подходе, не дав при этом полного его математического обоснования, используя лишь разрозненные формулы для описания отдельных частей процесса моделирования. Насколько мы видим, на сегодняшний день полной формулировки алгоритма так и нет.
В нашей работе мы предлагаем применить локальные оценки метода Монте-Карло к решению уравнения глобального освещения.
Решение уравнения глобального освещения (2) может быть разложено в ряд Неймана [5]
После ряда преобразование выражение (4) может быть интерпретирована, как Марковская цепь с вероятностью перехода, определяемой ядром уравнения
В результате построения Марковской цепи мы можем проводить оценку яркости в заданной точке по заданному направлению на поверхности сцены. Такая оценка может быть названа локальной оценкой метода Монте-Карло, так как позволяет вычислять непосредственно яркость в заданной точке поверхности сцены в заданном направлении.
Построить аналогичную схему для уравнения глобального освещения для точки находящейся в объеме сцены (3) мы не можем. Однако, качество освещения планируется определять именно на основе пространственно-углового распределения яркости – то есть необходимо иметь возможность получать угловое распределение яркости для произвольной точки в пространстве сцены. Построить аналогичную схему для уравнения глобального освещения для точки находящейся в объеме сцены (3) мы не можем. Однако, качество освещения планируется определять именно на основе пространственно-углового распределения яркости – то есть необходимо иметь возможность получать угловое распределение яркости для произвольной точки в пространстве сцены. В уравнении появляются две дополнительные δ-функции и
, которые зависят от направления Ȋ. Это делает невозможным непосредственное моделирование уравнения. При моделировании яркости в заданной точке и направлении пространства сцены
, мы не можем с узлов траектории марковской цепи попасть в требуемое нам направление Ȋ. Для этого нам требуется еще один дополнительный узел, фиксирующий промежуточною точку
уравнения (3). Такой подход получил название двойной локальной оценки [5]. Двойная локальная оценка позволяет проводить моделирования уравнения глобального освещения для точки находящейся в пространстве сцены и тем самым получать пространственно-угловое распределение яркости.
3. Качество освещения
3.1. Критерий качества на основе ПУРЯ
На сегодняшний день при проектировании неспециальных осветительных установок, таких как офисные помещения, общественные места, магазины, торговые центры и др., фактически используется один единственный критерий, описывающий качество установки – объединенный показатель дискомфорта UGR. На дискомфорт, в пространственно-угловом распределении яркости, влияет не только абсолютная величина наблюдаемой яркости, но и отношение перепада яркости источник-фон к яркости фона (яркости адаптации) – контраст [6]. Отношение контраста к пороговому контрасту и может служить критерием качества освещения. В случае непрерывного пространственно-углового распределения яркости по сцене освещения естественным обобщением контраста является отношение градиента распределения яркости по полю наблюдения к средней по полю яркости [7][8]. С увеличением значения градиента граница между блеским источником и фоном будет становиться более очерченной, а качество освещения, соответственно, уменьшаться. То есть чем больше источник и чем выше градиент яркости вокруг яркого источника, тем больший вклад в дискомфорт вносит этот источник. Отметим, что в реальной жизни протяженные блики как являются источниками дискомфорта, так и вносят свой вклад в яркость адаптации ПУРЯ. Обобщенный контраст в точке сцены можно определить [9]:
, (4)
где
. (5)
– координаты точки на проекции сцены,
– яркость данной точки в направлении
наблюдения, – средняя по полю зрения яркость, p(x,y) – некоторая весовая функция, учитывающая различный вклад в реакцию глаз точек, расположенных в центре поля зрения и на периферии. В формуле критерия p несет тот же физический смысл, что и индекс позиции в формуле UGR.
Таким образом, можем сформулировать критерий качества освещения Q как средне взвешенный по полю контраст K(x,y), отнесенный к некоторому порогу:
(6)
где Kпор – пороговое значение контраста.
Предполагается, что пороговый контраст будет определяться зрительной задачей.
3.2. Определение пороговой яркости
Очевидно, что в предложенном виде любое малозаметное изменение яркости будет давать вклад в критерий качества, так как будет изменение градиента яркости. При этом очевидно, что изменения контрастов в яркостях ниже некоторого предела вносить реальный вклад не будут. Так, например, очевидно, что при наличии в поле зрения непосредственно источника света или блика от него в помещении, контрасты в темном углу скорее всего роли играть в восприятии качества освещения не будут. Но если формально проводить вычисление, то и они дадут вклад. Таким образом, контрасты ниже некоторого порога Lпор фактически не должны быть учтены.
В 1946 году Blackwell провел огромную исследовательскую работу по установлению пороговых контрастов при решении задачи обнаружения [10]. В работе экспериментально была установлена взаимосвязь между пороговым контрастом и яркостью адаптации при различных угловых размерах цели. На рисунке 2 представлены графики этой зависимости.
Рис. 2. Зависимость логарифма порогового контраста от логарифма яркости адаптации для 5 угловых размеров: 121.0, 55.2, 18.2, 9.68, 3.60 угловых минут
Опираясь на это исследование в качестве Lпор можно рассматривать некоторое количество порогов. То есть, задав минимальный размер элемента, который требуется обнаружить – например это может быть размер символа при чтении текста или вывески с определенного расстояния. Тогда зная яркость адаптации, которая может быть взята как средняя яркость по полю, можно определить по результатам исследований Blackwell’а – пороговое изменение яркости для решения задачи обнаружения. Тогда, задав некоторое число превышений порогов, можно определить уже пороговую яркость как
(7)
где N – некоторое число, зависящее от отличия реально решаемой задачи от пороговой (как правило, для надпороговых задач N лежит в диапазоне 5-10).
То есть, в качестве пороговой яркости отсечения контрастов, не влияющих на критерий качества, мы берем некоторое число порогов, например 10ΔLпор.
Тогда выражение для контраста может быть записано как
(8)
Частью нашей работы стал эксперимент по определению качества освещения станций московского метрополитена. В эксперименте наблюдатели визуально определяли «качество» освещения по 10 бальной шкале при выполнении стандартной зрительной задачи – чтения вывески. Также параллельно с помощью яркомера и широкодиапазонного фотоаппарата проводились замеры яркости станций с последующей обработкой и расчетом критерия качества. На рисунке 3 представлена карта рассеяния оценки наблюдателей и критерия качества освещения. Коэффициент корреляции при этом составил 0.61 при числе порогов 50. Что является сильной корреляцией по шкале Чеддока.
Рис. 3. Карта рассеяний средней оценки наблюдателей от критерия качества освещения для 21 станции московского метрополитена
На рисунке 4 представлена зависимость коэффициента корреляции от числа порогов N. Наибольшая корреляция достигается в интервале десятков порогов обнаружения Блеквелла.
Рис. 4. Зависимость коэффициента корреляции от числа порогов
4. Выводы
Из исследований ещё 1915 года известно, что пространственно-угловое распределение яркости влияет в существенной степени на воспринимаемое наблюдателем качестве освещения. На основе ПУРЯ возможно производить оценку качественных характеристик освещения, в том числе и в основном, нормируемого на сегодняшний день, обобщённого показателя дискомфорта UGR. Однако, до недавнего времени, не было ни математических методов, ни вычислительных возможностей моделировать именно ПУРЯ. Таким образом, существующая нормативная база и инструменты инженеров проектировщиком осветительных установок построены вокруг моделирования излучательности.
Однако на сегодняшний день мы имеем ситуацию, когда с точки зрения теории, математических методов решения, алгоритмов, вычислительных возможностей, физических инструментов измерения светотехническое проектирование готово к переходу на анализ подлинно воспринимаемой глазом человека характеристике – яркости. И как следствие этого процесса возможна смена парадигмы проектирования от проектирования на заданные количественные характеристики, к проектированию с учетом качества освещения на основе ПУРЯ.
1. Ferree, C. Rand, G. A Preliminary Study of the Deficiencies of the Method of Flicker for the Photometry of Lights of Different Color // Psychological Review, 1915. - V.22, № 2. - pp. 110-162.
2. Kajiya, J.T. The rendering equation // Computer Graphics (Proc. SIGGRAPH’86), 1986. V.20, N4. pp.143-150.
3. Budak V. Relation of instant radiosity method with local estimations of Monte Carlo method / V. Budak , V. Zheltov, R. Notfulin, V Chembaev V. // Journal of WSCG, 2016. - pp. 189-196
4. Keller, A. Instant radiosity // SIGGRAPH '97 Proceedings of the 24th annual conference on Computer graphics and interactive techniques. - pp. 49-56.
5. Marchuk, G. I. Metod Monte-Karlo v atmosfernoj optike // Novosibirsk, Nauka, 1976. - pp. 283.
6. Scheir, G. H. Calculation of the Unified Glare Rating based on luminance maps for uniform and non-uniform light sources / G. H. Scheir, P. Hanselaer, P. Bracke, G. Deconinck and W. R. Ryckaert // Building and Environment, 2015. - № 84. - pp. 60-67.
7. Budak, V.P. Fast and accurate algorithm for the numerical simulation of radiative transfer in turbid media / V.P. Budak, V.S. Zheltov, A.V. Lubenchenko, K.S. Freidlin, O.V. Shagalov // Optikaatmosferyiokeana, 2016. - Vol. 26. - pp. 127-134.
8. Budak, V.P. Evaluation of illumination quality based on spatial-angular luminance distribution / V.P. Budak, V.S. Zheltov, T.V. Meshkova, R.Sh. Notfullin // Light & Engineering, 2017. - T 25, № 4. - pp. 24-31.
9. Budak, V.P. ChembaevExperemental study of the new criterion of lighting quality based on analysis of luminance distribution at Moscow metro stations / Vladimir P. Budak, V.S. Zheltov, T.V. Meshkova// Light & Engineering, 2020. - Vol. 28, №. 3.- pp. 98-105.
10. Richard, H. Blackwell Contrast Thresholds of the Human Eye // Josa, 1946. - Volume 36, Number 11. - pp. 624-643.