BC/NW 2007, №1, (10) :16.9

ПОСТРОЕНИЕ РЕАЛИСТИЧНЫХ ЛАНДШАФТОВ С ПОМОЩЬЮ СТОХАСТИЧЕСКИХ ФРАКТАЛОВ С ПРИМЕНЕНИЕМ БЫСТРОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ

А.С. Целищев, В.А. Логинов

(Москва, Московский энергетический институт (технический университет), Россия)

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

Известно, что в общем случае амплитуда гармоник, составляющих случайный сигнал, равна A_f = k / f ^γ . При абсолютно случайном сигнале (белый шум) γ = 0, и амплитуда всех гармоник, как высших, так и низших одинакова, и соседние значения сигнала не зависят друг от друга. Если γ стремится к бесконечности, амплитуды всех гармоник равны нулю, а сигнал – константа. При генерации ландшафта необходимо принять за высоту каждой точки поверхности значения случайного двумерного сигнала с 0<γ<∞. Предложен алгоритм получения такого сигнала для последующего построения ландшафта:

1. Генерация матрицы D_ст, n*n элементов. D_ст[i,j] - случайное целое от 0 до r, где r – расстояние между низшей и высшей точками желаемого ландшафта.

2. Проведение двумерного БПФ над D_ст. F_ст – матрица коэффициентов БПФ размерностью n*n .

3. Изменение коэффициентов двумерного БПФ для получения амплитудного спектра A_f=k / f ^γ при сохранении значений фазового спектра (что обеспечит уникальность ландшафта, обусловленную изначальной случайностью элементов D_{ст [i,j]}). В силу специфики двумерного БПФ коэффициенты низших гармоник будут находиться ближе к угловым элементам матрицы F_ст., а высших – ближе к центральным. Таким образом, новые коэффициенты необходимо считать в четырёх циклах (либо в одном объединённом) для каждой четверти матрицы F_ст. Так, например, для коэффициентов F_ст[i,j]: i=0..n/2, j=0..n/2 (верхняя левая часть матрицы F_ст) новые коэффициенты предложено рассчитывать следующим образом:

F_{ст [i,j]}=a_ст+b_ст*i – коэффициент двумерного БПФ после шага 2.

F_{нов [i,j]}=a_нов+b_нов*i – коэффициент двумерного БПФ, посчитанный на шаге 3, где

a_нов=A_[i,j]*cos(φ_[i,j]) , b_нов=A_[i,j]*sin(φ_[i,j]) , где A_[i,j] = k/ (i+j) ^γ – новый модуль коэффициента F_{ст [i,j]},

φ_[i,j] = arctg(b_ст/a_ст) – аргумент F_{ст [i,j].}

4. Проведение обратного БПФ, получение D_нов.

D_нов – это стохастический фрактал, представленный картой высот [0,r]. Построив его для всех значений D_нов и наложив на поверхность текстуры, можно получить искомый ландшафт. Варьируя γ и k, можно менять характер ландшафта, варьируя элементы D_ст[i,j] – его внешний вид. Пример ландшафта, полученного с помощью написанной программы, представлен на рис. 1.