48137 (608476), страница 2

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 2)

Принцип работы алгоритма.

1. Из проверяемой точки строится луч, направленный на источник света.

2. Производится поиск пересечений этого луча с примитивами сцены между проверяемой точкой и источником.

3. Если найдено хотя бы одно пересечение, то проверяемая точка находится в тени. При расчете ее цвета источник, для которого проводилась проверка, не учитывается.

4. Если пересечений не найдено, точка не в тени. При расчете ее цвета учитываем проверяемый источник.

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

РЕАЛИСТИЧНЫЕ ТЕНИ

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

Алгоритм работает следующим образом.

1. Из проверяемой точки строится луч, направленный на источник света.

2. Находятся ВСЕ пересечения этого луча с примитивами сцены между проверяемой точкой и источником и сохраняются в массив пересечений. Первым элементом в массив заносится 0 (сама точка, для которой истроится тень). Одновременно с этим формируется массив, в котором хранятся данные о том, какому объекту сцены принадлежит найденная точка пересечения.

3. Массив пересечений сортируется по возрастанию расстояния от проверяемой точки до точки пересечения. При этом абсолютно аналогичные операции проводятся и со вторым массивом.

4. Если в массив не была добавлена ни одна точка кроме начальной, то проверяемая точка не затенена. Берется полная интенсивность источника. В противном случае начинаем проход массива со второй точки по всем точкам.

5. Если объект, которому принадлежит текущая точка пересечения, абсолютно непрозрачный, то имеет место сплошная тень, цикл прерываем. Иначе проверяем, входит луч в объект или выходит из него. Если входит – переход в следующей точке.

6. Если выходит, то вычисляем расстояние от предыдущей точки до текущей. Это будет расстояние, пройденное светом внутри объекта. В зависимости этого расстояния ослабляем интенсивность источника. Переходим к следующей точке.

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

Математические и физические предпосылки алгоритма обратной трассировки лучей

Освещение

Согласно модели Уиттеда цвет некоторой точки объекта определяется суммарной интенсивностью

I(l)=Ka*Ia(l)C(l) + Kd*Id(l)*C(l) +Ks*Is(l) + Kr*Ir(l) +Kt*It(l),

где l – длина волны, C(l) – заданный исходный цвет точки объекта, Ka, Kd, Ks, Kr и Kt – коэффициенты, учитывающие свойства конкретного объекта параметрами фоновой интенсивности, диффузного рассеивания, зеркальности, отражения и прозрачности.

Ia – интенсивность фоновой подсветки,

Id – интенсивность, учитываемая для диффузного рассеивания,

Is – интенсивность, учитываемая для зеркальности,

Ir – интенсивность излучения, пришедшая по отраженному лучу,

It – интенсивность излучения, пришедшая по преломленному лучу.

Интенсивность фоновой подсветки (Ia) для некоторого объекта обычно константа. Запишем формулы для остальных интенсивностей.

Для диффузного отражения:

Id=∑I(l)cosӨ,

где I(l) – интенсивность излучения i-го источника света, Ө - угол между нормалью к поверхности объекта и направлением на i-й источник света.

Для зеркальности:

Id=∑I(l)cosα,

где p – показатель степени от единицы до нескольких сотен (согласно модели Фонга), α – угол между отраженным лучом (обратной трассировки) и направлением на i-й источник света.

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

Для вычисления коэффициента при интенсивности отраженного луча используется следующая формула:

r = 0.5(cos α - Ncos β)² /(cos α + Ncos β)² + 0.5(Ncos α – cos β)² /(Ncos α + cos β)²,

где α – угол падения луча, β – угол преломления луча, N – относительный показатель преломления двух сред.

Коэффициент при интенсивности преломленного луча вычисляется по формуле:

t = 1 – r.

Схема рассчета интенсивности. параметры, задающие свойства тел

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

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

Интенсивность фоновой подсветки (I_A) задается некоторой константой.

Интенсивность диффузно отраженного света (I_D) вычисляется по классическому «закону косинуса».

I_D = I_L cos α,

где I_L – интенсивность источника света, α – угол между нормалью к поверхности и направлением на источник.

В случае освещения сцены несколькими источниками I_d вычисляется для каждого из них и затем суммируются.

I_Di =Σ I_Li cos α_i.

Интенсивность блика от источника (I_S) вычисляется в соответствии с моделью Фонга.

_S = I_L cos^p β,

где I_L – интенсивность источника света (0<=I_L<=1), β – угол между отраженным лучом от источника света и направлением на точку, в которой расположена камера (центр проекции), p – некоторая степень от 1 до 200 –влияет на размытость блика. При маленьких значениях p блик более размытый.

Как и при вычислении I_D в случае освещения сцены несколькими источниками IS вычисляется отдельно для каждого источника, а затем результаты суммируются.

_Si =Σ I_Li cos^p β _i.

Интенсивности зеркально отраженного (I_R) и преломленного (I_T) света рассчитываются для отраженного и преломленного лучей на следующем шаге рекурсии. Если достигнут предел глубины рекурсии, то эти интенсивности берутся нулевыми. От интенсивности I_R берется r процентов, а от I_T – t = 1 – r (см. предыдущий раздел).

Кроме того, вводятся следующие коэффициенты: KD – коэффициент диффузного отражения поверхности, KS – коэффициент блика.

KD – этот коэффициент является характеристикой неровности отражающей поверхности. Чем больше неровность поверхности, тем меньше света отражается от неё зеркально и меньше света она пропускает, и соответственно больше света она отражает диффузно. 0 <= K_D <= 1.

При K_D = 0 - весь свет, падающий на поверхность, отражается и преломляется. K_D = 1 - весь свет отражается диффузно. На этот коэффициент умножаются интенсивность диффузно отраженного света и интенсивность фоновой подсветки. Интенсивности зеркально отраженного и преломленного света умножаются на (1 – K_D).

K_S – этот коэффициент отвечает за яркость блика от источника. 0<=KS<=1. При K_S = 0 - блик не виден, при KS = 1 – яркость блика максимальна.

Таким образом, окончательная формула для расчета интенсивности объекта в какой-либо точке будет выглядеть следующим образом:

I = I_AK_D + Σ(I_LiK_Dcos α_i + I_LiK_Scos^p β _i) + (1 – K_D )(I_Rr + I_T(1 – r)).

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

Для получения цветного изображения необходимо провести расчеты отдельно для красной, зеленой и синей компоненты света. Цвет пиксела изображения будет вычисляться путем умножения каждой компоненты интенсивности на число, определяющее максимальное количество градаций интенсивности изображения. Для 32-битного изображения оно равно 255 на каждый из цветов(R,G,B).

R = 255*I_R,

G = 255*I_G,

B = 255*I_B.

Здесь I_R (не путать с интенсивностью зеркально отраженного света), I_G, I_B – интенсивности трех компонент света в точке, полученная по формуле, указанной выше.

Коэффициенты K_D, K_S, p – это индивидуальные характеристики объекта, отражающие его свойства. Кроме этого имеется еще один коэффициент – абсолютный показатель преломления n. n = c / v, где c – скорость света в вакууме, v – скорость света в среде (внутри объекта). Для абсолютно непрозрачных тел этот коэффициент равен ∞ (т.к. скорость света внутри тела нулевая). В программе для задания абсолютно непрозрачного тела необходимо поставить этот коэффициент >> 1 (порядка 10 000). При этом доля зеркально отраженного света r будет стремиться к единице, а преломленного, соответственно, к нулю.

Вычисление нормалей

В алгоритме трассировки нормали к объектам необходимы для вычисления отраженного и преломленного лучей, а также для определения освещенности согласно модели Фонга.

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

Вычисление нормали к полигону (Треугольнику).

Вычисление нормали к треугольнику сводится к операции векторного умножения. Пусть задан треугольник ABC координатами трех своих вершин: X_A, Y_A, Z_A, X_B, Y_B, Z_B, X_C, Y_C, Z_C.

Вычислим координаты двух векторов, например AB и AC:

X_AB = X_B – X_A,

Y_AB = X_B – X_A,

Z_AB = X_B – X_A,

X_AC = X_C – X_A,

Y_AC = X_C – X_A,

Z_AC = X_C – X_A.

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

X_n = Y_ABZ_AC - Y_ACZ_AB,

Y_n = X_ABZ_AC - X_ACZ_AB,

Z_n = X_ABY_AC - X_ACY_AB.

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

Вычисление нормали к поверхности второго порядка.

Поверхность второго порядка задается в общем случае уравнением вида: Q(x,y,z) = a₁x² + a₂y² + a₃z² + b₁yz + b₂xz + b₃xy + c₁x +c₂y +c₃z + d =0.

Но мы будем использовать другую форму записи. Так уравнение эллипсоида будет выглядеть следующим образом:

(x-x₀)²/A² + (y-y₀)²/B² + (z-z₀)² /C² = 1,

где x₀, y₀, z₀ – координаты центра эллипсоида, A, B, C – длины полуосей эллипсоида. Уравнение параболоида:

(x-x₀)²/A² + (y-y₀)²/B² - (z-z₀)² /C² = 1,

где x₀, y₀, z₀ – координаты центра параболоида, A, B, C – длины полуосей параболоида. Ось параболоида расположена вдоль оси Oz мировой системы координат.

Для вычисления координат вектора нормали необходимо вычислить частные производные по x, y, z.

Координаты вектора нормали эллипсоида:

X_n = 2(x-x₀)/A²,

Y_n = 2(y-y₀)/B²,

Z_n = 2(z-z₀)/С².

Направление вектора не изменится, если все его координаты разделить на 2:

X_n = (x-x₀)/A²,

Y_n = (y-y₀)/B²,

Z_n = (z-z₀)/С².

Координаты вектора нормали параболоида вычисляются аналогично:

X_n = (x-x₀)/A²,

Y_n = (y-y₀)/B²,

Z_n = – (z-z₀)/С².

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

Вычисление отраженного луча

Пусть задан вектор падающего луча S, а также известен вектор нормали N. Требуется найти вектор отраженного луча R.

Рассмотрим единичные векторы R₁, S₁и N₁. Поскольку векторы нормали, падающего луча и отраженного луча находятся в одной плоскости, то можно записать R₁ + S₁ = N`, где N` - это вектор, соответствующий диагонали ромба и совпадающий по направлению с нормалью. Длина вектора N` равна 2cosθ. Так как вектор N` по направлению совпадает с N₁, то

N` = N`2cosθ.

Отсюда найдем единичный вектор отраженного луча:

R₁ = N₁ 2cosθ – S₁ = N/|N| 2cosθ – S/|S|.

Найдем cosθ. Это можно сделать, используя скалярное произведение векторов N и S:

cosθ = N*S/(|N|*|S|).

Полагая, что искомый вектор отраженного луча будет иметь такую же длину, что и вектор падающего луча, то есть R = |S| R₁, получим

R = N 2NS/|N|² – S.

Это решение в векторной форме. Запишем координаты вектора:

x_R = 2x_N(x_Nx_S+y_Ny_S+z_Nz_S)/(x_N²+y_N²+z_N²) – x_S,

y_R = 2y_N(x_Nx_S+y_Ny_S+z_Nz_S)/(x_N²+y_N²+z_N²) – y_S,

z_R = 2z_N(x_Nx_S+y_Ny_S+z_Nz_S)/(x_N²+y_N²+z_N²) – z_S.

Вычисление преломленного луча

Пусть даны два единичных вектора: S₁ – вектор падающего луча, и N₁ – вектор нормали к границе раздела двух сред. Также должны быть известны два коэффициента преломления для данных сред – n₁ и n₂ (или их отношение).

Требуется найти единичный вектор преломленного луча T₁. Для решения выполним некоторые геометрические построения.

Искомый вектор T₁ равен сумме двух векторов:

T₁ = N_T + B.

Найдем вначале вектор NT. Он противоположен по направлению вектору нормали, а его длина равна |T₁| cos α₂ = cos α₂ (поскольку T₁ – единичный). Таким образом, N_T = -N₁ cos α₂. Необходимо определить cos α₂. Запишем закон преломления n₁ sin α₁ = n₂ sin α₂ в виде:

sin α₂ = n sin α₁,

где n = n₁ / n₂.

Воспользуемся тождеством cos²α + sin²α = 1. Значение cos α₁ можно выразить через скалярное произведение единичных векторов S₁ и N₁, то есть cos α₁ = S₁N₁. Тогда мы можем записать такое выражение для вектора N_T:

N_T = -N₁√1+n²((S₁N₁)² – 1).

Осталось найти выражение для вектора B. Он располагается на одной прямой с вектором A, причем A = S₁ – N_S. Учитывая, что N_S равен N₁ cos α₁, то A = S₁ – N₁ cos α₁. Так как cos α₁ = S₁N₁, то A = S₁ – N₁ (S₁N₁).

Поскольку длина вектора A равна sin α₁, а длина вектора B равна sin α₂,

|B|/|A| = sin α₂/ sin α₁ = n₂/n₁ = n,

откуда |B| = n |A|. Учитывая взаимное расположение векторов A и B, получим

B = –nA =n(N₁(S₁N₁) – S₁).

Теперь мы можем записать искомое выражение для единичного вектора луча преломления T₁:

Вычисление точки пересечения с примитивами

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

R = A + Vt,

где R – радиус вектор произвольной точки, принадлежащей лучу, A – радиус- вектор начальной точки луча, V – направляющий вектор луча, t – параметр.Если направляющий вектор V нормализовать, то параметр t будет численно равен расстоянию от начальной точки луча A до точки R.

Можно записать это уравнение в координатном виде:

x = x₁ + at,

y = y₁ + bt,

z = z₁ + ct.

Здесь x1, y1, z1 – координаты начальной точки луча в прямоугольной декартовой мировой системе координат, a,b,c – координаты направляющего вектора луча.

Вычисление точки пересечения луча с поверхностью второго порядка.

Для нахождения точки пересечения луча, заданного уравнениями (2) с поверхностью второго порядка, заданной уравнениями (2.2.2.3) или (2.2.2.4):

(x–x₀)²/A² + (y–y₀)²/B² + (z–z₀)² /C² = 1 (эллипсоид)

(x–x₀)²/A² + (y–y₀)²/B² – (z–z₀)² /C² = 1 (параболоид),

нужно подставить в уравнение поверхности второго порядка вместо x, y и z соответствующие уравнения луча. В результате этого после раскрытия всех скобок и приведения подобных мы получим квадратное уравнение относительно параметра t. Если дискриминант квадратного уравнения меньше нуля, то луч и поверхность второго порядка общих точек пересечения не имеют. В противном случае можно будет вычислить два значения параметра t. Дискриминант может быть равен нулю – это соответствует предельному случаю касания луча поверхности, и мы получим два совпадающих значения параметра t.

Для нахождения координат точек пересечения луча и поверхности достаточно подставить найденные значения параметра t в уравнения луча (2).

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

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

Для этого после нахождения точки пересечения анализируется ее z-координата.

Вычисление точки пересечения луча с полигоном (Треугольником).

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

Уравнение плоскости выглядит следующим образом:

Q(x, y, z) = Ax + By + Cz +D = 0.

Здесь коэффициенты A, B, C совпадают с координатами нормали к этой плоскости. Координаты нормали плоскости совпадают с координатами нормали треугольника, которые мы посчитали на этапе загрузки сцены.

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

D = –Ax –By – Cz.

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

Теперь для нахождения точки пересечения подставим уравнения луча (2) в уравнение плоскости.

A (x₁ + at) + B (y₁ + bt) + C (z₁ + ct) + D = 0

t = – (Ax₁ + By₁ + Cz₁ + D) / (Aa + Bb + Cc)

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

Для нахождения координат точки пересечения надо подставить найденное значение параметра t в уравнения луча (2). Назовем точку пересечения D. Мы получим координаты x_D, y_D, z_D.

Теперь необходимо определить, попала ли точка D внутрь треугольника. Найдем координаты векторов AB, BC, CA (A, B, C – вершины треугольника) и координаты векторов AD, BD, CD. Затем найдем три векторных произведения:

nA = AB x AD,

nB = BC x BD,

nC = CA x CD.

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

Операцию проверки принадлежности точки D треугольнику ABC можно ускорить. Если ортогонально спроецировать треугольник ABC и точку D на одну из плоскостей xOy, yOz или xOz, то попадание проекции точки в проекцию треугольника будет означить попадание самой точки в треугольник (конечно же, если уже известно, что точка D лежит в плоскости, содержащей треугольник ABC). При этом число операций заметно сокращается. Так для поиска координат всех векторов нужно искать по две координаты на каждый вектор, а при поиске векторных произведений нужно искать только одну координату (остальные равны нулю).

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

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

Остается добавить, что для проецирования лучше выбирать ту из плоскостей, площадь проекции треугольника на которую больше. При таком условии исключается случай проецирования треугольника в отрезок (при условии, что проверяемый треугольник не вырожден в отрезок). Кроме того, при увеличении площади проекции уменьшается вероятность ошибки. Для определения такой плоскости проецирования достаточно проверить три координаты нормали треугольника. Если z-координата нормали больше (по абсолютному значению) x и y, то проецировать надо на плоскость xOy. Если y больше чем x и z, то проецируем на xOz. В оставшемся случае – на yOz.

Описание типов данных. Структура программы

Список модулей:

TTex.h - описание структуры TTex

TSurfTex.h - описание структур TPlaneTex и TEllipsoidTex

TPoint.h - описание структур TPoint2d и TPoint3d

TRGBColor.h - описание страктуры TRGBColor

TLamp.h - описание класса TLamp

TCam.h - описание класса TCam

TPrimitive.h - описание класса TPrimitive

TFrstSurface.h - описание класса TFrstSurface

TScndSurface.h - описание класса TScndSurface

TTriangle.h - описание класса TTriangle

TEllipsoid.h - описание класса TEllipsoid

TParaboloid.h - описание класса TParaboloid

TScene.h - описание класса TScene

TTracer.h - описание класса TTracer

Модули реализующие, интерфейс программы:

_AboutUnit.h - модуль формы «О программе»

_ZoomUnit.h - модуль формы «Лупа»

_Options.h - модуль формы «Опции»

_ExtraGlassOptions.h - модуль формы «Свойства стекла»

_ExtraTableOptions.h - модуль формы «Свойства стола»

_ExtraCamOptions.h - модуль формы «Свойства камеры»

_MainUnit.h - модуль главной формы программы

Рис. 2.3.1. Схема связей между модулями программы.

Рис.2.3.2. Схема наследования примитивов

Краткое описание структур и классов программы

struct TPoint3d – структура, описывающая точку в мировой системе координат

struct TPoint2d – структура, описывающая точку на плоскости (в текстуре) с целочисленными координатами

struct TRGBColor – структура, описывающая цвет по трем составляющим (RGB)

struct TTex – структура, описывающая текстуру – содержит адрес массива пикселей и его размеры

struct TPlaneTex – структура, описывающая привязку текстуры к плоскости.

Содержит три точки, к которым привязывается текстура

class TLamp – класс, описывающий источник освещения.

Содержит объект TPoint3d coord с координатами источника и три переменные типа float (Ir, Ig, Ib) для хранения интенсивности трех компонент света.

class TCam – класс, описывающий камеру.

Содержит два угла (a, b), указывающих направление зрения камеры, точку, на которую направлена камера (viewP) и расстояние от камеры до этой точки (r).

class TPrimitive – абстрактный класс примитива. От него наследуются поверхности первого и второго порядка.

class TFrstSurface – абстрактный класс поверхности первого порядка. От него наследуется класс треугольника.

class TScndSurface – абстрактный класс поверхности второго порядка. От него наследуются классы эллипсоида и параболоида.

class TTriangle – класс треугольника. Содержит три вершины треугольника и его нормаль

class TParaboloid – класс параболоида.

class TEllipsoid – класс эллипсоида.

class TScene – класс сцены. Содержит информацию о всех примитивах, источниках и камере.

class TTracer – класс, отвечающий за построения изображения. Содержит буфер (buffer) разметом 400x400 пикселей, в котором формируется изображение сцены. Перед генерацией необходимо вызвать функцию

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

Все классы – потомки TPrimitive предоставляют следующие функции:

float getT(TPoint3d p0, TPoint3d viewDir) – возвращает расстояние от точки начала(p0) луча viewDir до ближайшей точки пересечения с примитивом void getTArr(float* arr, int& n, TPoint3d p0, TPoint3d viewDir) – заполняет массив arr расстояниями от точки начала(p0) луча viewDir до ближайшей всех точек пересечения с примитивом void getNormal(TPoint3d& n, const TPoint3d& p) – возвращает координаты вектора нормали к примитиву в точке p void getColor(TRGBColor& c, const TPoint3d& p) - возвращает цвет примитива точке p (с учетом текстуры).

2. Технологическая часть

Выбор языка программирования

При разработке программы был использован язык программирования высокого уровня C++ в составе среды визуального программирования CBuilder6.

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

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

Описание интерфейса программы

Главное меню программы:

- Файл

- Сохранить картинку – записать в файл содержимое окна вывода изображения

- Выход – выход из программы

- Настройки

- Опции – открытие формы «Опции»

- Инструменты

- Лупа – открытие формы «Лупа»

- Генерировать – запуск генерации изображения

- Помощь

- О программе – сведения о разработчике и версия программы

Для генерации изображения нажмите кнопку ГЕНЕРИРОВАТЬ или выберите аналогичный пункт в главном меню.

В нижней строке во время генерации изображения выдается сообщение ИДЕТ ГЕНЕРАЦИЯ ИЗОБРАЖЕНИЯ. По завершении генерации здесь будет отображено время, затраченное на генерацию.

На этой вкладке находятся средства по настройке освещения сцены.

Координаты источника – координаты в мировой системе координат источника света, выбранного в выпадающем списке.

Интенсивность источника – значения трех компонент интенсивности источника света, выбранного в выпадающем списке.

Фоновая интенсивность – значения трех компонент фоновой интенсивности.

Кнопка “+” (рядом с выпадающим списком) – добавление нового источника света.

Кнопка “–” (рядом с выпадающим списком) – удаление источника света, выбранного в выпадающем списке.

ФОРМА «ОПЦИИ». ВКЛАДКА «КАМЕРА».

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

Предосмотр – здесь можно увидеть примерный вид изображения до его генерации.

Навигация – настройки положения камеры.

Дополнительно – при нажатии на эту кнопку появляется форма Свойства камеры с дополнительными параметрами камеры.

ФОРМА «СВОЙСТВА КАМЕРЫ».

Радиус –расстояние от камеры до точки, на которую она направлена.

Шаг изменения радиуса – приращение радиуса камеры при однократном нажатии кнопки “–” на вкладке “Камера” формы “Опции” (или уменьшение при однократном нажатии кнопки “+”).

ФОРМА «ОПЦИИ». ВКЛАДКА «МАТЕРИАЛЫ».

При нажатии на кнопку «Стекло» отображаются свойства материала стакана.

Цвет – цвет материала стакана.

Коэф. диффузного отражения – коэффициент K_D материала стакана (см. раздел 2.2.1). Коэф. преломления – абсолютный коэффициент преломления материала стакана. Дополнительно – при нажатии на эту кнопку появляется форма Свойства стекла с дополнительными параметрами материала стакана.

ФОРМА «СВОЙСТВА СТЕКЛА».

Коэффициент блика – коэффициент K_S материала стакана (см. раздел 2.2.1). Размытость блика – показатель степени p материала стакана. Расстояние затухания – максимальное расстояние, которое может пройти свет в материале до полного затухания. При нажатии на кнопку «Стол» отображаются свойства материала стола.

Цвет – цвет материала стола.

Коэф. диффузного отражения – коэффициент K_D материала стола (см. раздел 2.2.1).

Текстура – если галочка установлена, то на столе будет отображаться текстура.

Выбрать текстуру – выбор файла изображения (*.bmp), который будет использоваться как текстура стола.

Дополнительно – при нажатии на эту кнопку появляется форма Свойства стола с дополнительными параметрами материала стола.

ФОРМА «СВОЙСТВА СТОЛА».

Коэффициент блика – коэффициент K_S материала стола (см. раздел 2.2.1).

Размытость блика – показатель степени p материала стола.

Повторения текстуры – сколько раз текстура стола будет повторяться вдоль осей OX и OY.

ФОРМА «ОПЦИИ». ВКЛАДКА «СИСТЕМНЫЕ».

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

Глубина рекурсии – этот параметр устанавливает глубину рекурсии в алгоритме трассировки. При бОльших значениях этого параметра качество сгенерированного изображения улучшается.

ВНИМАНИЕ!

Глубина рекурсии СИЛЬНО влияет на скорость генерации изображения. Не рекомендуется ставить значения этого параметра больше 10.

Анитиалиазинг – включение алгоритма сглаживания изображения.

Тип тени – выбор алгоритма построения теней.

Лупа позволяет рассмотреть сгенерированное изображение в увеличенном виде.

Можно менять масштаб увеличения лупы (1:1, 1:2, 1:4, 1:8, 1:10, 1:16) в выпадающем меню.

3. Исследовательская часть

Исследования проводились на компьютере со следующей конфигурацией:

CPU AMD Athlon 2100+

512 Mb DDR SDRAM

WindowsXP

Зависимость времени генерации от глубины рекурсии

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

T1 – время генерации без тени в секундах.

T2 – время генерации со сплошной тенью в секундах.

T3 – время генерации с реалистичной тенью в секундах.

N – глубина рекурсии.

Зависимость времени генерации от количества источников

Анализ результатов исследований

Из первого исследования видно, что график зависимости времени генерации изображения от глубины рекурсии по виду совпадает с графиком функции 2ⁿ. Это хорошо соответствует теории, т.к. количество лучей растет с увеличением глубины рекурсии как 2ⁿ.

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

Во втором исследовании показано, что зависимость времени генерации от количества источников света линейна. Из полученных значений можно вычислить время, необходимое для расчета одного источника. На машине, на которой проводились исследования, при глубине рекурсии 5 это время примерно равно 0,8 секунды.

Заключение

В этой программе были продемонстрированы результаты роботы алгоритма генерации реалистичных изображений – обратной трассировки лучей.

Данная реализация демонстрирует возможности алгоритма строить изображения близкие к фотореалистичным. Трассировка является одним из самых совершенных алгоритмов генерации реалистичных изображений. Качество получаемого изображения несравнимо лучше, чем качество изображения, полученного с помощью таких алгоритмов, как Z-буфер. Однако требования к вычислительным мощностям, необходимым для генерации одного кадра изображения намного выше, чем в том же Z-буфере. Это делает невозможным применение этого алгоритма для реализации приложений, где необходима генерация изображений в реальном времени (3D-игры, различного рода симуляторы). Существуют методы повышения быстродействия алгоритма. Для увеличения скорости генерации изображения применяют такие средства, как метод порталов и метод оболочек.

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

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

Но даже со всеми усовершенствованиями на сегодняшний день реализация трассировки в реальном времени представляется затруднительной. Известны попытки создания 3D-ускорителей с аппаратной поддержкой алгоритма обратной трассировки, однако широкого распространения они не получили из-за относительно высокой цены и недостаточного быстродействия (удовлетворительная скорость построения изображения была только для изображения с разрешением 640x480).

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

Список литературы

1. Роджерс Д. Алгоритмические основы машинной графики: пер. с англ.— М.: Мир, 1989.— 512 с.: ил.

2. Порев В. Н. Компьютерная графика. – СПб.: БХВ-Петербург, 2002. – 432 с.: ил.

3. Никулин Е. А. Компьютерная геометрия и алгоритмы машинной графики. СПб.: БХВ-Петербург, 2003. – 560с.: ил.

4. Авдеева С.М., Куров А.В. Алгоритмы трехмерной машинной графики: Учебное пособие. – М.: Изд-во МГТУ им. Н.Э. Баумана, 1996. – 60 с.: ил.

39

Характеристики

Список файлов курсовой работы