1 / 72

Relationale Wachstumsgrammatiken und GroIMP: Spezifikation multiskalierter Struktur-Funktions-Modelle von Pflanzen mit 3

Relationale Wachstumsgrammatiken und GroIMP: Spezifikation multiskalierter Struktur-Funktions-Modelle von Pflanzen mit 3D-Grafik-Komponente http://www.grogra.de. Pflanzenmodelle: bisherige Ansätze in der Biologie. Prozessmodelle Einheiten: Stoffpools in Kompartimenten

gilda
Download Presentation

Relationale Wachstumsgrammatiken und GroIMP: Spezifikation multiskalierter Struktur-Funktions-Modelle von Pflanzen mit 3

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Relationale Wachstumsgrammatiken und GroIMP: Spezifikation multiskalierter Struktur-Funktions-Modelle von Pflanzen mit 3D-Grafik-Komponente http://www.grogra.de

  2. Pflanzenmodelle: bisherige Ansätze in der Biologie • Prozessmodelle • Einheiten: Stoffpools in Kompartimenten • Flüsse zwischen den Pools • Umweltfaktoren Beispiel: Modell der Ozon-wirkung (Martin et al. 2000)

  3. Beispiel: TREEDYN (Bossel 1994, Jansen et al. 1995)

  4. Prozessmodelle – Nachteile: • mangelhaft repräsentierte räumliche Struktur • Parameterfitting, bei großen Modellen zu viele Freiheitsgrade • keine Repräsentation der Objekte, mit denen der Pflanzenzüchter / • Botaniker / Agronom arbeitet • (b) Strukturmodelle • Einheiten: Organe / Module der Pflanze • (botanisch sinnvolle und visualisierbare Einheiten!) • Konkurrenzeffekte (Licht, Raum, Bodenressourcen) ergeben sich (z.T.) aus dem Modell • Parameter überschaubar, messbar

  5. Strukturmodelle – Beispiele: L-System-basierte Modelle (siehe später) Kurth 1998, 1999 Prusinkiewicz & Lindenmayer 1990

  6. Lane & Prusinkiewicz 2002

  7. AMAP (de Reffye et al., seit 1976, CIRAD, Montpellier): Atelier de la Modélisation de l'Architecture des Plantes • - modelliert die Aktivität von Meristemen • entwickelt in enger Zusammenarbeit von Informatikern und Biologen • weiterentwickelt über die letzten 25 Jahre: Aufnahme neuer Ideen (Markoffketten-Modelle, Referenzachse, Kohlenstoff-Allokations-Modelle...) • verschiedene, divergierende Software-Systeme für verschiedene Zwecke

  8. XFROG (Deussen & Lintermann): • interaktives System zur Pflanzenmodellierung, das sich am Szenengraphen-Konzept orientiert • botanisch nur schwach fundiert, keine Wachstumsmodellierung • aber grafisch zufriedenstellende Ergebnisse • hohe Form-Flexibilität

  9. Strukturmodelle – Nachteile: • keine (oder wenig) Berücksichtigung der funktionalen Aspekte von Pflanzen • kein Stoffhaushalt, keine Kopplung mit niedrigen Skalenebenen

  10. Logistische Sigmoidfunktion mit 3 Parametern (c) aggregierte Modelle (statistische Modelle) Beschreibung eines Kollektivs (räumlich / zeitlich) mit Mitteln der Statistik (Momente, Korrelation, Regressionsrechnung etc.) Beispiel: Modell der Blattstreckung bei Gerste (Buck-Sorlin 2004)

  11. "Modell-Dreieck":

  12. Struktur-Funktions-Modelle von Pflanzen (Functional-structural plant models, FSPM): Verbindung der beiden unteren Ecken des Dreiecks Vorgehensweise: Verteilung der Prozesse auf die Module Hilfsmittel: OOP aus Breckling (1996)

  13. Beispiel ALMIS (Eschenbach 2000)

  14. Beispiel LIGNUM (Perttunen et al. 1996, 1998, Dzierzon 2000; METLA 2006)

  15. Nachteile bisheriger ad hoc - FSPMs: • Einzellösungen, oft stark spezialisiert • umfangreicher, komplexer Quellcode, • der technische Einzelheiten gemischt mit • fundamentalen Modellstrukturen enthält

  16. geringe Kompatibilität der Modelle

  17. ?

  18. • in bisherigen Struktur-Funktions-Modellen wird oft nur 1 Konkurrenzfaktor, bzw. 1 funktionaler Aspekt behandelt • Ziel: "Metasprache" für alle Typen von Pflanzenmodellen – Wachstumsgrammatik • Verbesserung der Transparenz, • der Vergleichbarkeit und Kombinierbarkeit von Modellen, • der Reproduzierbarkeit der Ergebnisse • zusätzlich: Kompatibilität mit leistungsfähigen Grafik-Tools

  19. klassisches Simulationsprogramm

  20. klassisches Simulationsprogramm Interpreter für Wachstumsgrammatiken

  21. erster Ansatz: L-Systeme (Lindenmayer-Systeme) dafür verwenden analog zu Grammatiken für natürliche Sprachen in jedem Ableitungsschritt parallele Ersetzung aller Zeichen, auf die eine Regel anwendbar ist von A. Lindenmayer (Botaniker) 1968 zur Modellierung des Wachstums von fadenförmigen Algen eingeführt

  22. mathematisch: • Ein L-System ist ein Tripel (, , R); darin ist: • eine Menge von Zeichen, das Alphabet, • eine Zeichenkette mit Zeichen aus , das Startwort (auch "Axiom"), R eine Menge von Regeln der Form Zeichen  Zeichenkette; darin sind das Zeichen auf der linken Regelseite und die Zeichenkette aus  entnommen.

  23. Ein Ableitungsschritt (rewriting) einer Zeichenkette  besteht aus der Ersetzung aller Zeichen in , die in linken Regelseiten von R vorkommen, durch die entsprechenden rechten Regelseiten. • Man vereinbart: Zeichen, auf die keine Regeln anwenbar sind, werden unverändert übernommen. • Ergebnis zunächst nur: • Ableitungskette von Wörtern, die sich durch wiederholte Anwendung des rewriting-Vorgangs aus dem Startwort ergeben. • 123 ....

  24. Beispiel: Alphabet {A, B}, Startwort A Regelmenge R: A  B B  AB Ableitungskette: A  B  AB  BAB  ABBAB  BABABBAB  ABBABBABABBAB  BABABBABABBABBABABBAB  ... wie lang ist die n-te Zeichenkette in dieser Ableitung?

  25. was für die Modellierung von räumlichen Strukturen noch fehlt: eine geometrische Interpretation Füge also zur Def. eines L-Systems hinzu: eine Abbildung, die jeder Zeichenkette mit Zeichen aus  eine Teilmenge des 3-dimensionalen Raumes zuordnet

  26. was für die Modellierung von räumlichen Strukturen noch fehlt: eine geometrische Interpretation Füge also zur Def. eines L-Systems hinzu: eine Abbildung, die jeder Zeichenkette mit Zeichen aus  eine Teilmenge des 3-dimensionalen Raumes zuordnet dann: "interpretierte" L-System-Abarbeitung 123 ....  S1S2S3 .... S1, S2, S3, ... können als Entwicklungs- oder Entwurfsstufen eines Objekts interpretiert werden.

  27. Als Interpretationsabbildung wird meistens gewählt: Turtle geometry ("Schildkrötengeometrie") befehlsgesteuertes, lokales Navigieren im 2D- oder 3D-Raum (Abelson & diSessa 1982; vgl. Programmier-sprache "LOGO")

  28. Als Interpretationsabbildung wird meistens gewählt: Turtle geometry ("Schildkrötengeometrie") befehlsgesteuertes, lokales Navigieren im 2D- oder 3D-Raum (Abelson & diSessa 1982; vgl. Programmier-sprache "LOGO") "Turtle": Zeichen- oder Konstruktionsgerät (virtuell) - speichert (grafische und nicht-grafische) Informationen - mit einem Zustandsspeicher assoziiert (wichtig für Verzweigungen) - aktueller Zustand der Turtle enthält z.B. Information über aktuelle Liniendicke, Schrittweite, Farbe, weitere Eigenschaften des als nächstes zu konstruierenden Objekts

  29. Der Turtle-Befehlsvorrat wird zu einer Untermenge der Zeichenmenge  des L-Systems. Symbole, die nicht Turtle-Befehle sind, werden von der Turtle ignoriert. Befehle (Auswahl): F0 "Forward", mit Konstruktion eines Elements (Linienstück, Segment, Gebäudetrakt...), benutzt wird die aktuelle Schrittweite für die Länge (die Null steht für "keine explizite Längenfestlegung") M0 forward ohne Konstruktion (Move-Befehl) L(x) ändere die aktuelle Schrittweite (Länge) zu x LAdd(x) inkrementiere die aktuelle Schrittweite um x LMul(x) multipliziere die aktuelle Schrittweite mit x D(x), DAdd(x), DMul(x) analog für die aktuelle Dicke

  30. RU(45) Drehung der turtle um die "up"-Achse um 45° RL(...), RH(...) analog um "left" und "head"-Achse up-, left- und head-Achse bilden ein rechtwinkliges, räumliches Koordinatensystem, das von der turtle mit-geführt wird RV(x) Rotation "nach unten" mit durch x vorgegebener Stärke was ist das Ergebnis der Interpretation der Zeichenkette L(10) F0 RU(45) F0 RU(45) LMul(0.5) F0 M0 F0?

  31. Wiederholung von Abschnitten der Zeichenkette möglich mit dem Schlüsselwort "for" z.B. for ((1:3)) ( A B C ) liefert A B C A B C A B C was ist das Ergebnis der Interpretation von L(10) for ((1:6)) ( F0 RU(90) LMul(0.8) ) ?

  32. Verzweigungen: Realisierung mit Speicher-Befehlen [ lege aktuellen Zustand auf Speicher ("Ablage") ] nimm obersten Zustand von der Ablage und mache diesen zumaktuellen Zustand (damit: Ende der Verzweigung)

  33. Beispiel: Regeln a  F0 [ RU(45) b ] a ; b  F0 b ; Startwort L(10) a (a und b werden normalerweise nicht geometrisch interpretiert.) (in der Abbildung wurde F statt F0 geschrieben)

  34. was für eine Struktur liefert das L-System A  [ LMul(0.25) RU(-45) F0 ] F0 B; B  [ LMul(0.25) RU(45) F0 ] F0 A; mit StartwortL(10)A ?

  35. was für eine Struktur liefert das L-System A  [ LMul(0.25) RU(-45) F0 ] F0 B; B  [ LMul(0.25) RU(45) F0 ] F0 A; mit Startwort L(10)A ? äquivalente Regel: A  [ LMul(0.25) RU(-45) F0 ] F0 RH(180) A;

  36. Weitere Beispiele: Koch'sche Kurve:  L(50) RU(90) A F0; A  A LMul(0.3333); /* Skalierung */ F0  F0 RU(-60) F0 RU(120) F0 RU(-60) F0; jedes Linienstück wird durch 4 neue Linienstücke ersetzt (3. Regel); Skalierung durch Hilfssymbol A, welches sich in jedem Schritt reproduziert und dabei jeweils einen zusätzlichen Faktor 1/3 erzeugt (2. Regel). Das Startwort ist hier "  ".

  37. Ausgabe nach 6 Schritten:

  38. Sierpinski-Dreieck (Realisierung als geschlossene Kurve, Verwendung von Hilfssymbol X für Insertion des inneren Dreiecks):  L(50) RU(90) B F0 X F0 RU(-120) F0 F0 RU(-120) F0 F0; F0  F0 F0; X  RU(-120) F0 X F0 RU(120) F0 X F0 RU(120) F0 X F0 RU(-120); B  B LMul(0.5);

  39. Verzweigungsbeispiel: F0 F0 [ RU(25.7) F0 ] F0 [ RU(-25.7) F0 ] F0 ; Ergebnis nach 7 Schritten: (Startwort L(10) F0)

  40. Verzweigung, alternierende Zweigstellung und Verkürzung:  L(10) F0 A ; A  LMul(0.5) [ RU(90) F0 ] F0 RH(180) A ;

  41. welche Struktur liefert  F(10) A ; A  [ RU(-60) F(6) RH(180) A Sphere(3) ] [ RU(40) F(10) RH(180) A Sphere(3) ]; Sphere  Z; ? (F(n) liefert Linie der vorgegebenen Länge n, Sphere(n)eine Kugel mit Radius n)

  42. Erweiterung des Konzepts: • Lasse reellwertige Parameter nicht nur bei Turtle-Kommandos wie "RU(45)" und "F(3)" zu, sondern bei allen Zeichen • parametrische L-Systeme beliebig lange, endliche Parameterlisten • Parameter werden bei Regel-Matching mit Werten belegt • Beispiel: • Regel A(x, y)  F(7*x+10) B(y/2) • vorliegendes Zeichen z.B.: A(2, 6) • nach der Regelanwendung: F(24) B(3) • Parameter können in Bedingungen abgeprüft werden • (Bedingungen mit Java-Syntax): • A(x, y) (x >= 17 && y != 0)  ....

  43. Welche Struktur wird von folgendem L-System erzeugt?  [ RU(90) M(1) RU(90) A(1) ] A(1); A(n)  F(n) RU(90) A(n+1);

  44. Welche Struktur wird von folgendem L-System erzeugt?  [ RU(90) M(1) RU(90) A(1) ] A(1); A(n)  F(n) RU(90) A(n+1); Variante: in der zweiten Regel "RU(90)" etwa durch "RU(92)" ersetzen.

  45. Nachteile von L-Systemen: • in L-Systemen mit Verzweigungen (über Turtle-Kommandos) nur 2 mögliche Relationen zwischen Objekten: "direkter Nachfolger" und "Verzweigung" Erweiterungen: • Zulassen weiterer Relationstypen (beliebig wählbar) • Zulassen von Zyklen (Graph-Grammatik)

  46. • Grammatik modifiziert direkt den Graphen, Umweg über String-Codierung entfällt (bzw. wird nur noch für Regel-Input gebraucht) "relationale Wachstumsgrammatik"

  47. außerdem Nachteil der Turtle-Interpretation von L-Systemen:Segmente sind nur Zylinder, keine Objekte im Sinne der OOP  Erweiterungen: • Knoten des Graphen können beliebige Objekte sein (auch Grafikobjekte) • Einbettung von Code einer höheren, imperativen oder objektorientierten Programmiersprache in die Regeln (für uns: Java)

  48. Relationale Wachstumsgrammatiken (RGG) Aufbau einer Regel einer RGG:

More Related