imported>Nnamllit11: /* growthexperiments-addlink-summary-summary:1|0|0 */

2025-10-02T22:07:15Z

growthexperiments-addlink-summary-summary:1|0|0

Neue Seite

'''Radiosity''' bzw. '''Radiosität''' ist ein Verfahren zur Berechnung der Verteilung von Wärme- oder Lichtstrahlung innerhalb eines virtuellen Modells. In der [[Bildsynthese]] ist Radiosity neben auf [[Raytracing]] basierenden [[Algorithmus|Algorithmen]] eines der beiden wichtigen Verfahren zur Berechnung des Lichteinfalls innerhalb einer [[Szenengraph|Szene]]. Es beruht auf dem [[Energieerhaltungssatz]]: Alles [[Licht]], das auf eine Fläche fällt und von dieser nicht absorbiert wird, wird von ihr zurückgeworfen. Außerdem kann eine Fläche auch selbstleuchtend sein.

Das Radiosity-Verfahren basiert auf der Annahme, dass alle Oberflächen [[Lambertsches Gesetz|ideal diffuse]] Reflektoren bzw. alle Lichtquellen ideal diffuse Strahler sind. Ideal diffus bedeutet dabei, dass Licht in alle Richtungen gleichmäßig reflektiert bzw. abgestrahlt wird.

Im Gegensatz zu Raytracing ist Radiosity nicht vom Blickpunkt abhängig; die Beleuchtung der Flächen wird also für die gesamte Szene unabhängig von der Position des Betrachters berechnet. Die blickpunktabhängige [[Sichtbarkeitsproblem|Verdeckungsberechnung]] hat in einem unabhängigen Schritt zu erfolgen.

== Eigenschaften ==
=== Vorteile ===
[[Datei:Radiosity-no.jpg|mini|300px|Ein mittels [[Diffuses Raytracing|diffusem Raytracing]] gerendertes Bild ohne Simulation der indirekten Beleuchtung diffuser Oberflächen.]]

[[Datei:Radiosity-yes.jpg|mini|300px|Ein mittels Radiosity gerendertes Bild. Man erkennt deutlich, dass das Licht von den Kugeln auf den Boden reflektiert wird. Außerdem erkennt man auch, dass die weiße Wand die Kugeln von hinten indirekt beleuchtet.]]

; Blickpunktunabhängigkeit
: Ein Vorteil des Radiosity-Verfahrens ist, dass die Berechnung vom Standort und Blickwinkel des Betrachters unabhängig erfolgt. Die Lichtverteilung muss so für eine Szene nur einmal berechnet werden. Danach kann die Szene in Echtzeit [[Bildsynthese|gerendert]] werden (meist mittels [[Scanline-Algorithmus|Scanline-Algorithmen]] oder [[Z-Buffer]]ing), was für Anwendungen wie virtuelle Architekturmodelle interessant ist. Diesen Vorteil haben jedoch nicht alle Programme.
; Einfache indirekte, ideal diffuse Lichtreflexion
: [[Lambertsches Gesetz|Ideal diffuse]] Lichtreflexionen werden vom Radiosity-Verfahren auf natürliche Weise unterstützt. Die Helligkeit und Farbe einer Fläche werden nicht allein aufgrund der direkten Beleuchtung einer Lichtquelle, sondern auch durch diffus reflektiertes Licht anderer Flächen bestimmt. Ein Beispiel dafür ist ein Zimmer, das durch einfallendes Sonnenlicht nicht nur an den direkt beschienenen Stellen, sondern auch insgesamt heller wird.

=== Nachteile ===

; Keine analytischen [[Grafisches Primitiv|Primitiven]]
: Durch die zwingende Aufteilung der Szene in Polygone können keine analytisch definierten Primitiven, wie die bei [[Raytracing]] üblichen Kugeln, verwendet werden. Durch die zur Vermeidung von sichtbaren Kanten notwendige feine Aufteilung ist bei komplexer Szenengeometrie rasch eine sehr hohe Zahl von Flächen notwendig. Dies führt häufig zu einer hohen Wartezeit bei der Berechnung der Formfaktoren (siehe weiter unten).
: Es wurde versucht, mittels adaptativer Vereinfachung der Geometrie dieses Problem teilweise zu lösen,<ref>H. E. Rushmeier u. A: ''Geometric simplification for indirect illumination calculations.'' In: ''Proceedings of Graphics Interface ’93.'' Canadian Information Processing Society, Toronto 1993, S. 227–236 ({{Webarchiv | url=https://smartech.gatech.edu/bitstream/handle/1853/3656/92-01.pdf | archive-is=20160807161022 | text=smartech.gatech.edu}} PDF).</ref> was jedoch ein manueller Vorgang ist, der seine Grenzen hat. Daneben sind die aus den Fehlern dieses Verfahrens resultierenden Auswirkungen auf die Lichtberechnung schwer vorhersehbar.
; Hoher Speicher- und Zeitbedarf
: Wenn <math>N</math> die Anzahl der Primitiven einer Szene ist, so weisen die effizientesten Radiosity-Varianten eine asymptotische [[Zeitkomplexität]] von <math>O(N \log N)</math> auf, wie [[Empirie|empirisch]] festgestellt wurde.<ref>M. F. Cohen u. A.: ''Radiosity and realistic image synthesis.'' Academic Press Professional, San Diego 1993, ISBN 0-12-178270-0.</ref> Dagegen benötigt Raytracing nur eine Laufzeit von <math>O(\log N)</math>. Dies begrenzt schnell den praktisch möglichen Verfeinerungsgrad bei Radiosity. Daneben ergeben sich zur Berechnung und Speicherung der Formfaktoren relativ hohe Speicheranforderungen.
; Globale Beleuchtung schwer realisierbar
: Für die möglichst realistische Darstellung einer Szene muss die [[globale Beleuchtung]] simuliert werden, was aber nur in Spezialfällen mit Radiosity effizient möglich ist. In seiner Grundform ist Radiosity nur zur Simulation ideal diffuser Reflexion fähig. Die Berücksichtigung beliebiger [[Beleuchtungsmodell]]e sowie [[Transluzenz|transluzenter]] Flächen ist möglich,<ref>H. Rushmeier, K. Torrance: ''Extending the radiosity method to include specularly reflecting and translucent materials.'' In: ''ACM Transactions on Graphics.'' Band 9, S. 1–27, ACM Press, New York 1990 ([http://www.graphics.cornell.edu/pubs/1990/RT90.pdf graphics.cornell.edu] PDF)</ref><ref>E. Gobbetti u. A.: ''Hierarchical Higher Order Face Cluster Radiosity for Global Illumination Walkthroughs of Complex Non-Diffuse Environments.'' In: [[Eurographics (Vereinigung)|Computer Graphics Forum]], Band 22-3 (9/2003).</ref> hat jedoch keine breite Verwendung gefunden, da derartige Effekte schneller bzw. präziser mittels auf [[Raytracing]] basierter Lösungen möglich sind.

=== Vergleich mit Raytracing ===
Historisch gesehen war Radiosity interessant, da es die Simulation indirekter diffuser Beleuchtung auf einfache Weise erlaubte, was mit Raytracing lange Zeit nicht möglich war. Andererseits war Raytracing gut für spiegelnde und transparente Objekte geeignet, wozu wiederum Radiosity nicht fähig war. Es wurden daher anfangs Vorschläge zur Kombination von Radiosity mit Raytracing gemacht, die jedoch aufwendig waren und sich letztendlich nicht im großen Maße durchsetzen konnten.

Mit dem Aufkommen moderner globaler Beleuchtungsverfahren wie [[Path Tracing]] und [[Photon Mapping]] wurden jedoch die Möglichkeiten von Raytracing beträchtlich erweitert. Weil derartige Algorithmen alle von Radiosity unterstützten Effekte mit weniger Fehlern und auf elegantere Weise simulieren können, ist Radiosity im Bereich der hochwertigen realistischen [[Bildsynthese]] weitgehend aus der Mode gekommen. Kommerzielle Verwendung findet Radiosity vor allem beim Rendern von Architekturmodellen, bei denen eine zeitaufwendige Vorausberechnung vertretbar ist. Auch derartige Anwendungen sind jedoch mit Raytracing-basierten Verfahren ([[Particle Tracing]]) möglich.<ref>B. Walter u. A.: ''[http://www.graphics.cornell.edu/~bjw/detog.html Global Illumination Using Local Linear Density Estimation.]'' In: ''ACM Transactions on Graphics'', Vol. 16-3 (7/1997), S. 217–259, ACM Press, New York 1997</ref>

Weitere Anwendung finden Radiosity-Verfahren in den Bereichen der Klima- bzw. Wärmeforschung, da Wärmeverteilung eher diffus als gerichtet geschieht und Radiosity hier praktikabler als strahlenbasierte Ansätze ist.

== Prinzip ==
Mit der eingangs formulierten Annahme kann die allgemeine [[Rendergleichung]] in die Radiosity-Gleichung überführt werden.
: <math>B(x) = E(x)+\rho(x)\int_{x' \in S}B(x') \frac{1}{\pi r^2} \cos\phi_x\cos\phi_{x'} \cdot V(x,x') \,\mathrm dx'</math>
mit
: B(x) = gesamte vom Punkt x abgestrahlte Energie (Summe aus Eigenstrahlung und Reflexion als [[Leistung (Physik)|Leistung]] pro Flächeneinheit), genannt [[Bestrahlungsstärke|Radiosity]] im Punkt x
: E(x) = im Punkt x abgegebene Eigenstrahlung
: ρ(x) = [[Reflexion (Physik)|Reflexionsfaktor]] im Punkt x
: S = alle Oberflächen der Szene
: r = Abstand der Punkte x und x'
: <math>\phi (x)</math> = Der [[Winkel]] zwischen der [[Normalenvektor|Normale]] im Punkt x und der Verbindungslinie zwischen den Punkten x und x'
: <math>V(x,x')= \begin{cases} 1, & \mbox{wenn }x\mbox{ von }x'\mbox{ aus sichtbar} \\ 0, & \mbox{sonst}\end{cases}</math>.
Die gesuchte Radiosity im Punkt x ergibt sich aus einem Integral in geschlossener Form, was auf direktem Weg nicht berechnet werden kann. Abhilfe schafft die Diskretisierung der Oberfläche S: anstatt alle infinitesimal kleinen Teilflächen δA' zu betrachten, teilt man die Oberfläche S in zusammenhängende Teilflächen (''Facetten'' oder ''Patches'' genannt) Ai auf ([[Finite-Elemente-Methode]]). Für diese Teilflächen gelten weitere Annahmen: jedes Ai ist planar; die jeweilige Radiosity Bi und der Reflexionsfaktor ρi sind über Ai konstant. Dies führt dann zur diskreten Radiosity-Gleichung.
: <math>B_i = E_i + \rho_i \sum_{j=1}^n B_jF_{ij},\quad 1 \le i \le n</math>
mit
: Bi = gesamte von Teilfläche i abgestrahlte Energie (Summe aus Eigenstrahlung und Reflexion als Leistung pro Flächeneinheit), genannt [[Bestrahlungsstärke|Radiosity]] der Teilfläche i
: Ei = von der Teilfläche i abgegebene Eigenstrahlung
: ρi = [[Reflexion (Physik)|Reflexionsfaktor]] der Teilfläche i
: n = Anzahl der Teilflächen
: Fij = Anteil an der von Fläche j abgegebenen Energie, die auf Fläche i auftrifft, genannt [[Sichtfaktor|''Formfaktor'']].

Die Radiosity <math>B_{i}</math> von Teilfläche <math> i </math> ist also gleich der Eigenstrahlung <math>E_{i}</math> der Teilfläche
<math> i </math> plus der mit dem diffusen Reflexionsfaktor <math>p_{i}</math> gewichteten Summe der Radiosity <math>B_{j}</math> von allen anderen Teilflächen <math> j </math>. Wobei hierbei der Formfaktor <math>F_{i,j}</math>, welcher den Anteil an der von Fläche j abgegebenen Energie die auf Fläche i auftrifft festlegt, eingeht. Der Formfaktor berücksichtigt also die Ausrichtung und den Abstand der Teilflächen <math> i,j </math> zueinander. Da man dies für alle Teilflächen <math> i </math> berechnen muss, ergibt sich ein lineares Gleichungssystem mit so vielen Gleichungen und Unbekannten wie es Teilflächen gibt.

=== 1. Unterteilung der Oberflächen ===
Im ersten Schritt erfolgt die Festlegung der Primitiven: Wie und in welche Teilflächen soll eine gegebene kontinuierliche Oberfläche zerlegt werden? Üblich sind Dreieck und Quadrat. Bereits in dieser Phase wird zwischen Qualität und Effizienz entschieden. Je feinmaschiger das Netz, desto genauer sind die Ergebnisse, aber umso aufwendiger die Berechnungen.

In der Praxis benutzt man meist adaptive Verfahren. Ausgehend von einer z. B. triangulierten Oberfläche, werden nach dem hier angegebenen Schema die zugehörigen Radiosity-Werte aller Facetten bestimmt. Mit diesen Daten erfolgen dann weitere ortsabhängige Netzverfeinerungen. Ausschlaggebend dafür können sein: ein hoher Radiosity-Gradient benachbarter Facetten, Diskontinuitäten im Lichtverlauf (z. B. Lichtfleck) oder eine örtlich ungünstige Netzeinteilung (z. B. T-Knoten).
[[Datei:Radiosity.svg|mini|Bestimmung von Radiosity-Werten bei konstanten, linearen und quadratischen Basisfunktionen]]

=== 2. Festlegung der Basisfunktionen ===
Die diskrete Radiosity-Gleichung stellt eine Möglichkeit der Diskretisierung dar. Sie beruht auf der Annahme, dass die Radiosity über einer gegebenen Facette konstant ist, und verwendet daher konstante Basisfunktionen. Die Wahl von Basisfunktionen höheren Grads ist ebenfalls insbesondere mit dem [[Galerkin-Methode|Galerkin-Ansatz]] möglich.

=== 3. Berechnung der Formfaktoren ===
Der aufwendigste Schritt bei Radiosity ist, unabhängig von dem gewählten [[Algorithmus]], die Berechnung der Formfaktoren. Ein Formfaktor gilt immer zwischen zwei Patches und beschreibt die Menge der ausgetauschten Strahlung, liegt also zwischen null (keine Strahlung wird ausgetauscht) und eins (alle Strahlung wird ausgetauscht).

Der Formfaktor ist rein geometrischer Natur und wird durch die Stellung der Patches zueinander bestimmt. Außerdem spielt die Sichtbarkeit der Patches eine Rolle. Die Sichtbarkeitsberechnung braucht bei weitem die meiste Zeit in der Berechnung.

Die Formel für einen Formfaktor lautet:
: <math>F_{s,e} = \frac{1}{A_s}\int_{v\in A_e}\int_{u\in A_s} \frac{1}{\pi r^2} \cos\phi_u\cos\phi_v V(u,v) \,\mathrm dA_s \mathrm dA_e</math>
mit
: <math>F_{s,e}</math> = Formfaktor zwischen dem Sender <math>S</math> und dem Empfänger <math>E</math>
: <math>A_s</math> = Der Fläche des Senders <math>S</math>
: <math>A_e</math> = Der Fläche des Empfängers <math>E</math>
: <math>\phi_u</math> = Der [[Winkel]] zwischen der [[Normalenvektor|Normale]] des Senders <math>S</math> und der Verbindungslinie zwischen Sender <math>S</math> und Empfänger <math>E</math>
: <math>\phi_v</math> = Der [[Winkel]] zwischen der [[Normalenvektor|Normale]] des Empfängers <math>E</math> und der Verbindungslinie zwischen Empfänger <math>E</math> und Sender <math>S</math>
: <math>r</math> = Entfernung zwischen Sender <math>S</math> und Empfänger <math>E</math>
: <math>V(u,v)= \begin{cases} 1, & \mbox{wenn }u\mbox{ von }v\mbox{ aus sichtbar} \\ 0, & \mbox{sonst}\end{cases}</math>

Da das direkte Berechnen dieses Doppelintegrals sehr schwierig ist, werden im Allgemeinen Annäherungen verwendet.

Das simpelste Verfahren ist nur für Flächen korrekt, die relativ klein und relativ weit entfernt sind und zwischen denen keine partielle Verdeckung besteht.
Man berechnet dabei Winkel und Entfernung nur zwischen zwei repräsentativen Punkten, den Mittelpunkten der beiden Flächen.
: <math>F_{s,e} \approx A_e \frac{\cos\phi_u\cos\phi_v d_{u,v}}{\pi r^2} </math>
Wobei <math>u</math> der Mittelpunkt von <math>A_s</math> und <math>v</math> der Mittelpunkt von <math>A_e</math> ist.

==== Methode nach Nusselt ====
[[Datei:FormfaktorenfuerRadiosityBerechnung.PNG|mini|Ermittlung der Formfaktoren nach Nusselt]]
''(auch „Nusselts Analogon“)''

Es wird ein repräsentativer Punkt der Empfängerfläche, der Mittelpunkt, ausgewählt. Die sichtbaren Teile der Senderfläche werden auf die Einheitshalbkugel um diesen Punkt projiziert. Dadurch wird <math>\frac{\cos\phi_v }{r^2}</math> beachtet. Dann wird die Projektion auf der Einheitskugel wiederum in die Fläche, in der <math>A_j</math> liegt, projiziert. Die entstandene Fläche wird durch <math>\pi</math> geteilt (Fläche des Einheitskreises). Durch diesen Schritt wird der Formfaktor
<math>F_{s,e} \approx \int_{v\in A_e} \frac{\cos\phi_u\cos\phi_v d_{u,v}}{\pi r^2} \,\mathrm dA_e </math>
ermittelt.

==== Hemi-Cube-Verfahren ====
[[Datei:Hemicube Radiosity.png|mini|Approximation der Formfaktoren mittels Hemi-Cube]]
Von Cohen et al. stammt das sogenannte Hemi-Cube-Verfahren.<ref>M. F. Cohen, D. P. Greenberg: ''The hemi-cube: a radiosity solution for complex environments.'' In: ''Proceedings of the 12th annual [[SIGGRAPH|conference on Computer graphics and interactive techniques]]'', S. 31–40, ACM Press, New York 1985, ISBN 0-89791-166-0</ref> Die Einheitshalbkugel nach Nusselt wird durch einen Einheitshalbwürfel approximiert, dessen Seitenflächen in ein diskretes Gitter unterteilt sind. Jeder Gitterfläche wird ein Gewichtungsfaktor, der ''Delta-Formfaktor'', zugeordnet, welcher von der Position der Gitterfläche abhängig ist. Ein Delta-Formfaktor ist damit der Formfaktor der Gitterfläche nach Nusselt. Die Summe der Delta-Formfaktoren ist 1.

Für jede der fünf Halbwürfelflächen wird mit modifizierten Rasteralgorithmen (üblicherweise [[Z-Buffer]]) ein ''Item-Buffer'' berechnet. Dieser enthält für jede Gitterfläche die Identität der Objektfläche (''Item'', rotes Dreieck im Bild), die darauf projiziert wurde. Für jedes ''Item'' wird die Summe der Delta-Formfaktoren der überdeckten Gitterflächen berechnet (rote Gitterflächen im Bild). Diese Summe wird als Formfaktor zwischen ''Item'' und betrachteter Fläche aufgefasst.

==== Sillions Verbesserung ====
Statt eines Halbwürfels wird nun nur noch eine einzige Fläche verwendet, welche zentriert über dem differentiell kleinen Flächenstück ('''dA''') platziert wird. Diese Fläche wird ebenfalls in kleine, diskrete Bereiche aufgeteilt. Diesen werden dann genau wie beim Hemicube-Verfahren in Abhängigkeit von der Geometrie zu '''dA''', den sogenannten Delta-Formfaktoren zugeordnet. Der Vorteil bei diesem Verfahren ist, dass ein Patch nur auf diese eine Fläche projiziert werden muss und nicht mehr auf fünf Flächen eines Würfels. Außerdem ist dieses Verfahren legitim, da Patches, die orthogonal auf der Fläche von '''dA''' stehen, keinen großen Beitrag an der Gesamthelligkeit haben. Diese Beobachtung kann man sich über den Cosinus zwischen den Normalen der beiden Flächenstücke klarmachen.

=== 4. Berechnung der Radiosity-Werte ===
Die diskrete Radiosity-Gleichung kann als [[lineares Gleichungssystem]] aufgefasst werden und lässt sich demzufolge nach einigen Umformschritten wie folgt in [[Matrix (Mathematik)|Matrixform]] darstellen.
: <math>
\begin{pmatrix}
B_1 \\
B_2 \\
\vdots \\
B_n
\end{pmatrix}
=
\begin{pmatrix}
E_1 \\
E_2 \\
\vdots \\
E_n
\end{pmatrix}
+
\begin{pmatrix}
\rho_1 F_{11} & \rho_1 F_{12} & \ldots & \rho_1 F_{1n} \\
\rho_2 F_{21} & \rho_2 F_{22} & \ldots & \rho_2 F_{2n} \\
\vdots & \vdots & \ddots & \vdots \\
\rho_n F_{n1} & \rho_n F_{n2} & \ldots & \rho_n F_{nn}
\end{pmatrix}
\cdot
\begin{pmatrix}
B_1 \\
B_2 \\
\vdots \\
B_n
\end{pmatrix}
</math> oder kurz '''B''' = '''E'''+'''T'''•'''B'''.

Um alle gesuchten Radiosity-Werte zu bestimmen, muss das Gleichungssystem gelöst werden. Naheliegend ist die Invertierung der Matrix ('''I'''-'''T''') ([[Gaußsches Eliminationsverfahren|Matrix-Inversion nach Gauß]]), was aber aufgrund des enormen Aufwands unpraktikabel ist.

Im Allgemeinen existieren zwei unterschiedliche [[Iteration|iterative Lösungsstrategien]], die gegen die exakte Lösung [[Grenzwert (Folge)|konvergieren]]. Damit ist ein frei wählbarer Kompromiss zwischen Darstellungsqualität und Rechenzeit möglich.

* Beim '''Gathering''' wird die Radiosity Bi einer Facette durch ''Einsammeln'' aller einflussnehmenden Bj (abhängig vom jeweiligen Formfaktor) gebildet. Ausgangspunkt ist '''B'''0='''E''' und man sieht nur die selbstleuchtenden Flächen. Nach dem ersten Schritt '''B'''1='''E'''+'''T'''•'''B'''0 sind zusätzlich alle direkt beleuchteten Objekte sichtbar. Im nächsten Iterationsschritt '''B'''2='''E'''+'''T'''•'''B'''1 werden die einfachen Reflexionen berücksichtigt usw. Auffällig ist, dass in jedem Schritt alle Formfaktoren benötigt werden. Auf diesem Prinzip beruhende Verfahren sind [[Jacobi-Verfahren|Jacobi-Iteration]] und [[Gauß-Seidel-Verfahren|Gauß-Seidel-Iteration]].
* Beim '''Shooting''' wird die unverteilte Radiosity Bi auf alle als Empfänger in Frage kommenden Facetten (abhängig vom jeweiligen Formfaktor) ''geschossen''. Ausgangspunkt sind wiederum die selbstleuchtenden Flächen, somit gilt '''B'''='''E'''.
für jede Facette i
Radiosity B = Eigenleuchten E
unverteilte Radiosity ΔB = Eigenleuchten E
wiederhole
i = Facette mit maximaler Restenergie ΔBi•Ai
für jede Facette j
rad = ΔBi•ρj•Fji (wie viel der zu verteilenden Radiosity bekommt Facette j)
ΔBj = ΔBj+rad (zu verteilende Radiosity der Facette j wird um diesen Betrag erhöht)
Bj = Bj+rad (Radiosity der Facette j wird um diesen Betrag erhöht)
ΔBi = 0 (zu verteilende Radiosity der Facette i wurde verteilt und ist nun 0)
bis (Abbruchkriterium)

: Als Abbruchkriterien können neben objektiven Werten wie [[Fehler]]metriken auch der subjektive Eindruck benutzt werden. ''Shooting''-Verfahren sind den ''Gathering''-Verfahren in zwei Gesichtspunkten überlegen. Zum einen kann aus dem [[Pseudocode]] abgeleitet werden, dass für jeden Iterationsschritt lediglich j Formfaktoren benötigt werden. Zum anderen wird in jedem Durchlauf die Radiosity der Facette mit der größten Restenergie verteilt, wodurch ''Shooting''-Verfahren wesentlich schneller gegen „schöne“ Bilder konvergieren als ''Gathering''-Verfahren. Auf diesem Prinzip beruhende Algorithmen sind [[Southwell-Iteration]] und [[Progressive-Refinement]].

=== 5. Rendern ===
Den letzten Schritt bildet das [[Bildsynthese|Rendern]] des fertigen Bildes. Die an den ausgesuchten Stellen berechneten Radiosity-Werte werden gemäß den gewählten Basisfunktionen kombiniert. Wurden Radiosity-Werte in den Netzknoten bestimmt ist z. B. auch [[Gouraud Shading]] möglich. Genügt das so erstellte Bild dem Anspruch nicht, so können weitere Iterationen aus (4) folgen. Treten ungewollte grafische Artefakte auf, so sollte mit entsprechenden Änderungen in der Netzstruktur der Algorithmus von vorn gestartet werden.

== Geschichtliche Entwicklung ==
[[Datei:Radiosity - RRV, step 79.png|mini|Mit Radiosity gerendertes Bild]]

1984 wurde das Verfahren erstmals von Goral et al.<ref>C. Goral, K. E. Torrance, D. P. Greenberg und B. Battaile: [http://www.cs.rpi.edu/~cutler/classes/advancedgraphics/S07/lectures/goral.pdf ''Modeling the interaction of light between diffuse surfaces.''] (PDF; 1,1 MB) In: ''Computer Graphics.'' Band 18, Nr. 3.</ref> vorgestellt. Es stammt aus der [[Thermodynamik]], wo es verwendet wurde, um den Austausch von Wärmestrahlung zu berechnen ([[Strahlungsrechnung]]). Zu dieser Zeit wurde die Full-Form-Factor-Matrix-Methode zum Lösen des Gleichungssystems verwendet. Hier wird das Gleichungssystem für alle Formfaktoren zwischen allen Patches (=Flächen) der Welt aufgestellt und dann mittels eines mathematischen Verfahrens (meist dem Gauß-Seidel-Verfahren) gelöst. Prinzipiell entspricht dieses Vorgehen dem Einsammeln von Radiosity. Für jedes Patch wird berechnet, wie viel Licht es von jedem anderen Patch erhält. Das hat den Nachteil, dass die Berechnung der gesamten Matrix extrem viel Zeit benötigt und viel Speicherplatz belegt, was das Verfahren für komplexe Szenen unbrauchbar macht.

1988 wurde dann das '''Progressive-Refinement-Verfahren''' von Cohen et al. vorgestellt.<ref>M. F. Cohen u. A.: ''A progressive refinement approach to fast radiosity image generation.'' In: ''Proceedings of the 15th annual conference on Computer graphics and interactive techniques'', S. 75–84, ACM Press, New York 1988</ref> Hier wird der Prozess umgedreht, und das Licht wird nicht mehr an jedem Patch eingesammelt, sondern von jedem Patch verschossen. So kann man zuerst einmal das Licht von den Patches mit dem größten Radiositywert versenden und sich dann denen mit wenig Radiosity zuwenden. Hier ist es nicht mehr nötig, die gesamte Matrix zu berechnen, sondern es werden in jedem Schritt nur noch die Formfaktoren von einem einzelnen Patch zu allen anderen benötigt. Dadurch sinkt der benötigte Speicherplatz enorm, und man erhält nach jedem Schritt ein brauchbares Bild. Je länger man wartet, desto besser wird das Bild, da immer mehr Indirektionen berechnet werden. Um zu konvergieren, braucht dieses Verfahren allerdings genauso lange wie die Full-Form-Factor-Matrix-Methode.

Im Prinzip entspricht das Verfahren einer Reihenentwicklung der Matrix. Das [[Lineares Gleichungssystem|Lineare Gleichungssystem]], das sich aus der Strahlungsgleichung ergibt ist dann <math>(I-M) B = E</math>, wobei <math>I</math> die [[Einheitsmatrix]], <math>M</math> die Matrix, bestehend aus den Formfaktoren und den Reflexionskoeffizienten und <math>E</math> die Eigenstrahlung beschreibt. Stellt man die Gleichung nun um zu <math>B = (I-M)^{-1} E</math>, so lässt sich der rechte Teil <math>(I-M)^{-1}</math> als Reihe der folgenden Form entwickeln, und man erhält eine inkrementelle Annäherung an die tatsächliche Strahlungsintensität:

: <math>\sum_{i=1}^\infty x^i = \frac{1}{1-x} = (1-x)^{-1}</math>.

1991 wurde von Hanrahan et al. das '''Hierarchische Radiosity''' vorgestellt.<ref>Pat Hanrahan u. A.: ''[http://graphics.stanford.edu/papers/rad/ A Rapid Hierarchical Radiosity Algorithm].'' In: ''Proceedings of the 18th annual conference on Computer graphics and interactive techniques'', S. 197–206, ACM Press, New York 1991</ref> Bei diesem Verfahren wird eine Patchhierarchie verwendet. Wenige große Patches bestehen aus vielen kleinen. Der Lichtaustausch wird jetzt auf unterschiedlichen Stufen vorgenommen. Wenn der Fehler gering ist, wird das Licht auf einer hohen Hierarchieebene ausgetauscht, wenn ein großer Fehler zu erwarten ist (zum Beispiel wenn viel Licht ausgetauscht wird oder die Flächen sehr nahe beieinander liegen), dann wird das Licht auf einer niedrigeren Ebene ausgetauscht. Dadurch reduziert sich die Anzahl der zu berechnenden Formfaktoren substantiell, und die Berechnung wird deutlich beschleunigt.

Zusätzlich zu diesen Verfahren wurden noch viele Erweiterungen ersonnen. Zum Beispiel gibt es die Methode des '''Clustering'''<ref>B. Smits u. A.: ''A Clustering Algorithm for Radiosity in Complex Environments.'' In: ''Proceedings of the 21st annual conference on Computer graphics and interactive techniques.'' ACM Press, New York 1994, S. 435–442 ({{Webarchiv | url=http://www.cs.ucl.ac.uk/staff/j.mortensen/papers/radiosity/Smits_B__A%20Clustering_Algorithm_for_Radiosity_in_Complex_Environments.pdf | wayback=20040310110423 | text=cs.ucl.ac.uk}} PDF)</ref>, welches eine Erweiterung des Hierarchischen Radiosity ist. Hier wird oberhalb der Patchhierarchie eine weitere Hierarchie erzeugt, die Cluster. Licht kann dann auch zwischen ganzen Clustergruppen ausgetauscht werden, abhängig vom zu erwartenden Fehler. Wieder spart man sich die Berechnung vieler Formfaktoren.

== Literatur ==
* Michael F. Cohen, John R. Wallace: ''Radiosity and Realistic Image Synthesis''. Morgan Kaufmann, San Francisco 1993, ISBN 0-12-178270-0.
* François X. Sillion, Claude Puech: ''Radiosity and Global Illumination''. Morgan Kaufmann, San Francisco 1994, ISBN 1-55860-277-1.

== Quellen ==
<references />

{{Lesenswert|3. September 2005|9019771}}

[[Kategorie:Algorithmus (Computergrafik)]]

Radiosity (Computergrafik) - Versionsgeschichte

imported>Nnamllit11: /* growthexperiments-addlink-summary-summary:1|0|0 */