Projekt How-to

[freiBär]

Da es sicher einigen so geht, wie es mir gegangen ist, hier mal eine kleine Erklärung zum Projekt. Grundsätzlich stellt sich mal die Frage, was eigentlich zu machen ist. Die etwas schwammige Angabe findet sich in den Folien. Man findet sie inklusive Quellcode auf der LVA Seite: http://www.viskom.oeaw.ac.at/~milos/lecture/
Die Angaben variieren jedes Jahr ein bisschen, dieses Jahr lauten sie wie folgt:
(Very) simple rendering of 3D data sets
Projection in ±X, ±Y, ± Z directions
1. Simple rendering techniques
Sum, max, compositing FTB, BTF, isosurface
2.Submission and exam:
1. Select 6 images you like representing
each category, submit them in a pdf file
2. Submit source of the project

Das bedeutet grundsätzlich mal, dass man 3D Daten hat, die man jeweils in ±X, ±Y, ± Z Richtung projizieren soll. Dafür gibts eben ein paar Techniken, wie z.B. Sum. Projiziert man z.B. in +X Richtung, schießt man quasi Strahlen in X Richtung durch, und summiert die Helligkeitswerte auf. Dadurch erhält man eine Art Röntgen bild (beim Röntgen passiert ja auch nichts anderes). Zum Darstellen sollte man das ganze noch durch die Anzahl der durchlaufenen Pixel dividieren (also ist Sum eigentlich der Durchschnitt), da man ja nur Farbwerte zwischen 0 und 255 zeichnen kann.

So, und wie geht man das ganze jetzt an?
Also erstens braucht man den Quellcode von der Seite. Dann könnt ihr mal versuchen, dass ganze unter Windows laufen zu lassen, und das ganze schnell aufgeben. Weil erstens zu viele Fehler auftreten, und zweitens das Framework unter Windows falsche Ergebnisse liefert. Also ladet ihr euch VM Ware Player runter, und installiert das neue Ubuntu drauf. Geht ganz einfach und schnell. Dann kopiert ihr das Projekt auf den Ubuntu Desktop. Wenn ihr die Konsole öffnet, wechselt ihr mit cd Desktop/f3dProjBasicNoComp11 ins Projektverzeichnis. Gebt einfach "make" ein (ohne irgendwelche Parameter) um das Projekt zu kompilieren. Das ganze wird immer noch einen Fehler ausspucken. Das ganze lässt sich aber beheben, indem ihr euch von http://zlib.net/ zlib compiled DLL runterladet und entpackt. Dann kopiert ihr zlib1.dll, zconf.h und zlib.h (die Header Dateien sind im include Verzeichnis) in das Projektverzeichnis. Nun lässt sich das ganze kompilieren.

Mittels ./project -t 2 data/bigbrain-n.f3d > bigbrain-n.ppm lässt sich nun eine Projektion erstellen. Die lässt sich unter Ubuntu ohne zusätzliche Installation öffnen und betrachten.

So, wie löse ich nun die Aufgabe?
Öffnet project.cpp mit dem Editor. Ihr müsst nur in diesem File arbeiten, sonst nirgends. Ab ca. Zeile 85 fängt der Code an, für den ihr euch interessiert. Hier gibt es 3 for-Schleifen, die das Bild in alle 3 Richtungen durchlaufen. Die innerste Schleife (z-Schleife) ist die, die die Projektionsrichtung angibt. Mittels in->get(x, y, z) wird der aktuelle Pixel ausgelesen. Das heißt, ihr müsst euch damit nur ein bisschen spielen, um die Projektionsrichtung zu ändern. FTB wurde bereits implementiert, der Code für BTF ist auch schon vorhanden (in FTB, aber auskommentiert). Ihr müsst jetzt nur mehr die restlichen Funktionen ergänzen. Average entspricht hierbei Sum. Ihr ergänzt jetzt einfach bei Average den Case um I= I + in->get(x,y,z); und nach der z-Schleife um

switch(type){
case AVERAGE:
{
I= I/in->getNZ();
}
}

Damit habt ihr schon den Durchschnitt. Jetzt noch ein make in der Konsole und ihr könnt das ganze mit
./project -t 0 data/bigbrain-n.f3d > bigbrain-n.ppm testen. Das -t 0 entspricht dabei dem Sum-Modus. Welche Zahl welcher Modus ist, kann man in Zeile 10 sehen (die Reihenfolge der Einträge). Wenn ihr euch das Bild jetzt anschaut, habt ihr ein schönes Röntgenbildähnliches Bild in Z Richtung.

Ich glaub den Rest schafft ihr jetzt auch allein

[mOfl]

Hey, das ist ne super Sache, vielen Dank! Ich habs mir selber noch gar nicht angeschaut, weil ich mit den anderen Übungen noch beschäftigt bin, aber das freut mich, dass einer sich die Mühe gemacht hat, den Kollegen den Einstieg zu erleichtern, insbesondere weil ich - wie vermutlich die meisten hier mit Windows - erst ewig versucht hätten, das im Visual Studio zum Laufen zu kriegen. Wenn du sagst, dass es unter Windows falsche Ergebnisse gibt, muss ich mir die Mühe erst gar nicht machen, es da zu probieren.

[freiBär]

Habs am Anfang auch mit Visual Studio probiert. Im Forum gibts sogar irgendwo eine Anleitung dafür. Das Problem ist aber, dass dort einige berichten, dass es zu unterschiedlichen Ergebnissen kommen kann. Hier gibts übrigens eine einfache Installationsanleitung für das neue Ubuntu unter VM Ware: http://www.sysprobs.com/install-natt...l-vmware-tools

Geht einfach und relativ schnell (unter 1h). Hab mich diese Woche mit ein paar anderen aus der VU getroffen, und sind dann gemeinsam auf die oben genannten Schritte gekommen. Ist also nicht allein mein Verdienst Ich freu mich wenn ich helfen kann....immerhin kann man schon ziemlich viel Zeit in den Sand setzen, wenn man nicht genau weis was man machen soll...

[freiBär]

Weis wer, wos noch mehr Bilder gibt? Wir sollen ja von jeder Kategorie 6 Bilder zeigen, wir haben aber nur 3 3D Datasets...

[Kuoni]

ist das bei euch aus so, dass wenn man averaging oder auch sum macht, der kopf nicht erkennbar ist? sondern diese komischen querstreifen?

[Newbie0815]

jop hab ich auch @ kuoni ... ich hätte eine andere frage: wann, wie wo muss man den source etc abgeben? das pdf muss man zur prüfung mitnehmen und den soruce auch=?

[freiBär]

jop hab ich auch @ kuoni ... ich hätte eine andere frage: wann, wie wo muss man den source etc abgeben? das pdf muss man zur prüfung mitnehmen und den soruce auch=?

per mail schicken.

[mOfl]

Was soll man denn unter "case SHADED" implementieren? Wie habt ihr das gelöst? Unter shaded hätte ich spontan eine Gouraud- oder Phong-geshadete Isofläche des Datensatzes und dann Ray Casting mit "First"-Rückgabe verstanden, aber ich wage zu behaupten, dass das den "sehr einfachen" Umfang des Übungsteils sprengen würde, oder liege ich da falsch?

[freiBär]

Was soll man denn unter "case SHADED" implementieren? Wie habt ihr das gelöst? Unter shaded hätte ich spontan eine Gouraud- oder Phong-geshadete Isofläche des Datensatzes und dann Ray Casting mit "First"-Rückgabe verstanden, aber ich wage zu behaupten, dass das den "sehr einfachen" Umfang des Übungsteils sprengen würde, oder liege ich da falsch?

du verwechselst da irgendwas. wo hast dubittedie komische Angabe her? Die Angabe lautet doch eindeutig Sum, max, compositing FTB, BTF, isosurface, und das ist doch ziemlich einfach...

[mOfl]

"case SHADED" ist eine Zeile aus dem Code, den wir vervollständigen sollen, eben bei der Unterscheidung der verschiedenen Rendering Types. Der Tatsache, dass der Rendering Type nicht ISOSURFACE sondern SHADED heißt, habe ich entnommen, dass man sich hier auch noch um ein Beleuchtungsmodell kümmern soll.

[freiBär]

"case SHADED" ist eine Zeile aus dem Code, den wir vervollständigen sollen, eben bei der Unterscheidung der verschiedenen Rendering Types. Der Tatsache, dass der Rendering Type nicht ISOSURFACE sondern SHADED heißt, habe ich entnommen, dass man sich hier auch noch um ein Beleuchtungsmodell kümmern soll.

Vervollständigen ist relativ (steht auch nicht in der Angabe). Das ist nämlich seit Jahren der selbe Code, aber mit leicht veränderter Angabe. Aber um deine Frage zu beantworten: Threshold auswählen, sobald der überschritten wird Gradient (central difference estimation) bestimmen - fertig. Also ja, eigentlich hast du dadurch schon Beleuchtung drinnen. Anleitung dafür gibts schon - suche nach isosurface hilft

[mOfl]

Gradient (central difference estimation) bestimmen

Danke, ich glaube, das war das richtige Stichwort. Ich hatte die Schätzung des Gradienten nicht mit den benötigten Normalen für die Beleuchtung in Verbindung gebracht

[maxputz]

Vielleicht übersehe ich was, aber -

wo gibts die angaben zur UE? Laut freibär war die Angabe letztes Jahr auf den Folien. Das Framework findet sich online auf der LVA-Website. Aber wo genau findet man eine Angabe zur UE? Hat jdm. eine Idee?

[mOfl]

Puh da gräbst einen alten Schinken aus Die Angabe ist sehr vage, aber da du ohnehin nur grob 20 Zeilen ändern/schreiben musst, genügt das wohl auch. Die wenigen Informationen, die es dazu gibt, sollten sich in http://cg.tuwien.ac.at/courses/MedVi...uction2012.pdf finden lassen.

[freiBär]

Die Angabe findet sich in den Folien...ist ungefähr ein bis 2 Zeilen lang ;-)

Gesendet von meinem Nexus S mit Tapatalk

[maxputz]

Die Angabe findet sich in den Folien...ist ungefähr ein bis 2 Zeilen lang ;-)

Gesendet von meinem Nexus S mit Tapatalk

das ist doch alles ein Witz!

[freiBär]

Na was glaubst, warum ich das Projekt-How-To geschrieben hab

edit: Das ist die Angabe:

(Very) simple rendering of 3D data sets Projection in ±X, ±Y, ± Z directions
1. Simple rendering techniques
Sum, max, compositing FTB, BTF, isosurface
2.Submission and exam:
1. Select 6 images you like representing
each category, submit them in a pdf file
2. Submit source of the project

Sogar die selbe wie letztes mal.

[maxputz]

Na was glaubst, warum ich das Projekt-How-To geschrieben hab

hehe, danke nochmal!!!

[freiBär]

Also was wirklich ein Witz war: Die Aufgabe zu lösen, ohne mein Projekt-How-to zu haben

[maxputz]

Sie dürften dieses Semester den Source Code überarbeitet haben (damit man den oben beschriebenen Workaround mit den zlib Sachen nicht braucht, nehme ich an - bei funktionierts auch so)
Bei mir ist die Source ein Fragmenthaufen:
- als Parameter für die direction gibts nur +x und die direction (dir) selbst wir einmal auf MX gesetzt und kommt nie wieder vor.
- Sollen wir jetzt SHADED oder ISOSURFACE / SUM oder AVERAGE implementieren?
- Gabs dazu Infos in der Vorlesung?

[maxputz]

Achja nochwas: Wie verwendet man das Argument +x beim Programmaufruf richtig?

[Thomas G]

Erstmal riesigen Dank an die obige Anleitung und Hilfe!
Leider bin ich trotzdem nicht allzuweit damit gekommen. Ich scheitere an dem ausführen von "make" und zwar bekomme ich als Fehler

g++ -O -I. -o project.o -c project.cpp
make: g++: Command not found
make: *** [project.o] Error 127

Nun steht einmal oben noch was dazu
"Das ganze lässt sich aber beheben, indem ihr euch von http://zlib.net/ zlib compiled DLL runterladet und entpackt. Dann kopiert ihr zlib1.dll, zconf.h und zlib.h (die Header Dateien sind im include Verzeichnis) in das Projektverzeichnis. Nun lässt sich das ganze kompilieren."

Nur kann ich die entsprechenden files nicht finden? bzw weiß schonmal gar nicht was ich auf http://zlib.net/ runterladen muss.
Und maxputz meinte ja das sie dies in der jetzigen ss2012 version nicht mehr nötig ist.

Nur wie bringe ich das dann zum laufen?

[maxputz]

Erstmal riesigen Dank an die obige Anleitung und Hilfe!
Leider bin ich trotzdem nicht allzuweit damit gekommen. Ich scheitere an dem ausführen von "make" und zwar bekomme ich als Fehler

g++ -O -I. -o project.o -c project.cpp
make: g++: Command not found
make: *** [project.o] Error 127

Nun steht einmal oben noch was dazu
"Das ganze lässt sich aber beheben, indem ihr euch von http://zlib.net/ zlib compiled DLL runterladet und entpackt. Dann kopiert ihr zlib1.dll, zconf.h und zlib.h (die Header Dateien sind im include Verzeichnis) in das Projektverzeichnis. Nun lässt sich das ganze kompilieren."

Nur kann ich die entsprechenden files nicht finden? bzw weiß schonmal gar nicht was ich auf http://zlib.net/ runterladen muss.
Und maxputz meinte ja das sie dies in der jetzigen ss2012 version nicht mehr nötig ist.

Nur wie bringe ich das dann zum laufen?

Compiler installiert?

[Thomas G]

Lol. Nein... ;-) Welchen? Woher?

[maxputz]

http://wiki.ubuntuusers.de/GCC

[freiBär]

Sie dürften dieses Semester den Source Code überarbeitet haben (damit man den oben beschriebenen Workaround mit den zlib Sachen nicht braucht, nehme ich an - bei funktionierts auch so)
Bei mir ist die Source ein Fragmenthaufen:
- als Parameter für die direction gibts nur +x und die direction (dir) selbst wir einmal auf MX gesetzt und kommt nie wieder vor.
- Sollen wir jetzt SHADED oder ISOSURFACE / SUM oder AVERAGE implementieren?
- Gabs dazu Infos in der Vorlesung?

ISOSURFACE, nicht shaded. SUM und Average ist das selbe (also ihr müsst average implementieren). SUM könntest du nämlich grafisch gar nicht darstellen (weil die Werte größer als 255) sind - deshalb den Durchschnitt nehmen.

Für die Direction brauchst du ein switch (weis nicht obs schon im Code drinnen steht). Das ganze soll in alle Richtungen funktionieren.

switch(dir){
case PZ:{
...
break;
}
case PX:{
...
break;
}
......

Achja nochwas: Wie verwendet man das Argument +x beim Programmaufruf richtig?

Der Programmaufruf sollte in etwa so ausschaun:
./project +x -t 4 data/tot2-n.f3d > pics/ISO/tot2_px_4.ppm
+x ist die Direction, -t 4 ist dann für die Methode.

[Thomas G]

Wah ich peil nur einen Hauch von diesem Zeug.
1.Also mit dem Compiler (danke nochmal) scheint das nun funktioniert zu haben (siehe unten). (Wenn man die Warnings mal ignoriert)
2. Dann habe ich "./project -t 2 data/bigbrain-n.f3d > bigbrain-n.ppm" in die Konsole eingegebn das zur Folge hatte das nun dieses Röntgen Bild mit dem Kopf zu sehen ist im Ordner. Mehr hatte das nichts zu bedeutenoder?
3. Nun gehts weiter indem ich die project.cpp geöffnet habe und vor diesem Quellcode sitze (was ja schonmal gut ist, dass ich den sehe), Die besagten 3 for schleifen habe ich auch gefunden, nur was soll ich hier jetzt machen? Und wie kann ich mich "spielen" was kann ich hier ausprobieren und wie kann ich es ansehen?

Compiler
ordice@lordice-virtual-machine:~/Desktop/f3dProjBasicNoComp14$ make
g++ -O -I. -o project.o -c project.cpp
g++ -O -I. -o f3d.o -c f3d.c
f3d.c: In function ‘FILE* ReadHeader(FILE*, f3dHeader*, int*, int*, int*)’:
f3d.c:1046:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1059:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1069:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1248:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
g++ -O -I. -o f3dGridLite.o -c f3dGridLite.cpp
g++ -o project project.o f3d.o f3dGridLite.o -lm

[freiBär]

Wah ich peil nur einen Hauch von diesem Zeug.
1.Also mit dem Compiler (danke nochmal) scheint das nun funktioniert zu haben (siehe unten). (Wenn man die Warnings mal ignoriert)
2. Dann habe ich "./project -t 2 data/bigbrain-n.f3d > bigbrain-n.ppm" in die Konsole eingegebn das zur Folge hatte das nun dieses Röntgen Bild mit dem Kopf zu sehen ist im Ordner. Mehr hatte das nichts zu bedeutenoder?
3. Nun gehts weiter indem ich die project.cpp geöffnet habe und vor diesem Quellcode sitze (was ja schonmal gut ist, dass ich den sehe), Die besagten 3 for schleifen habe ich auch gefunden, nur was soll ich hier jetzt machen? Und wie kann ich mich "spielen" was kann ich hier ausprobieren und wie kann ich es ansehen?

Compiler
ordice@lordice-virtual-machine:~/Desktop/f3dProjBasicNoComp14$ make
g++ -O -I. -o project.o -c project.cpp
g++ -O -I. -o f3d.o -c f3d.c
f3d.c: In function ‘FILE* ReadHeader(FILE*, f3dHeader*, int*, int*, int*)’:
f3d.c:1046:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1059:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1069:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
f3d.c:1248:27: warning: ignoring return value of ‘char* fgets(char*, int, FILE*)’, declared with attribute warn_unused_result [-Wunused-result]
g++ -O -I. -o f3dGridLite.o -c f3dGridLite.cpp
g++ -o project project.o f3d.o f3dGridLite.o -lm

Zu 2. Ja das soll so sein.
Zu 3. Indem du das ganze so implementierst, das er die Bilder in der jeweiligen Richtung durchläuft. Du kannst dich dann e an die bereits implementierte Methode halten.

[maxputz]

Danke @freibär!

Was mich stört: Was ist dieses SUM, ISOSURFACE, MAX?
Ich hätte zumindest beim Durchblättern keine passenden Folien mit Formeln dazu gefunden, v.a. was ISOSURFACE betrifft.
War jdm. in der letzten Vorlesung - gabs Erklärungen zur Übung?

[freiBär]

SUM ist aufsummieren, was aber wie gesagt dem Average entspricht. Max ist eigentlich auch selbst erklärend (Maximum Intensity Projection). Isosurface steht irgendwo in den Folien (etwas mit Gradienten) bzw. gibts auch hier im Forum: http://www.informatik-forum.at/showt...4079-Projekt-1

[Thomas G]

Sorry: "Zu 3. Indem du das ganze so implementierst, das er die Bilder in der jeweiligen Richtung durchläuft. Du kannst dich dann e an die bereits implementierte Methode halten."
Ich bin wirklich eine totale Pfeife im Programmieren. Ich weiß wirklich nicht wie ich das implementieren soll.

Ich hab die for Schleifen schon gefunden. hab ich ja geschrieben^^
Nur weiß ich nicht wo ich da jetzt ansetzen muss und und wie oder was ich "implementieren muss, das die Bilder in der jeweiligen Richtung durchläuft"
Kann mir das jemand erklären. Sobald ich das verstanden habe sollte sich der rest ergeben.

for (x=0; x<in->getNX(); x++)
for (y=0; y<in->getNY(); y++){
float I(0); // accumulated color
float b(0); // accumulated alpha
for (z=0; z<in->getNZ(); z++) {
float val = in->get(x, y, z);
switch (type){
case COMP_FTB:
{
// voxel opacity
// we assume that maximum is 255
float a = 1-exp(-val/255);
// voxel color - just gray
float C = in->get(x, y, z)*a;
//FTB
I = I + (1-b)*C;
b = b + (1-b)*a;
//BTF
//I = I*(1-a) + C;
}

[freiBär]

Nun ja, die Basics einer Projektion enes 3D Volumens in en 2D Bild sollte man ehrlich gesagt schon drauf haben (ich glaub für die LVA wird auch "VO Visualisierung" empfohlen. Wenn du nicht gut Programmieren kannst, wirdst du dir hier sicher nicht leicht tun. Ich helf dir trotzdem: Du hast 3 in sich geschachtelte Schleifen, innerhalb werden die Pixel ausgelesen. Die 3 Schleifen entsprechen das durchlaufendes Volumens auf der x,y und z Achse. Wird jetzt in der untersten Schleife aufsummiert, und unterhalb der 2. Schleife (selbe Ebene wie die z Schleife) gespeichert, hast du ein 2D Bild, das für jeden Punkt (x,y) die Summe der Pixel entlang der z-Ahse enthält. Du hast also die Sum Projektion (dividiers noch durch die Länge der z- Achse, da du Werte > 255 nicht rendern kannst = Durchschnittsberechnung). Berechnest du nicht die Summe sondern den Maximalwert hast du die Maximum Intensity Projection (MIP)

[frankenstone]

Es wurde heute gesagt, dass die Prüfung am 21. Juni ist. Abgabe des Projekts soll einfach davor per Mail geschehen. Was auch ganz interessant ist, dass es einen Fragenkatalog geben wird. War das letztes Jahr auch schon und gibt es da vielleicht schon eine Ausarbeitung? Ich habe ein GoogleDocs angelegt: https://docs.google.com/document/d/1...mWFqoHdHw/edit falls sich jemand bei der Ausarbeitung beteiligen möchte.

[freiBär]

Nein, hats nicht gegeben.

[TehSackL]

Hi!

Nachdem ich leider beim isosurface nicht durchblicke (auch nicht mit dem anderen thread), frag ich mal nach

Bitte korrigieren wenn falsch aber soweit hab ichs verstanden und s

Zuerst threshold detection also am ersten pixel welcher den thresh übersteigt die weiteren Berechnungen machen.
Dann die gradient estimation (gx,gy,gz) wie in post nummer 24 in http://www.informatik-forum.at/showt...4079-Projekt-1 berechen.

Ab jetzt versteh ichs leider nicht mehr, aus was berechnet sich jetzt genau der pixel wert dann?

Danke im Voraus

[freiBär]

Steht doch e im Thread den du gepostet hast:
http://www.informatik-forum.at/showt...l=1#post589787

[TehSackL]

ja genau das hab ich dann nicht verstanden,

weil vorher wurde geschrieben dot multiply with the viewing direction

und er schreibt dann im post:
val = gx*vx + gy*vx + gz*gz
was für mich keinen sinn ergibt,

dot multiply with the viewing direction wär für mich sowas:
val = gx*vx + gy*vy + gz*vz
und damit schauts finde ich ok aus

aber stimmt jetzt meins oder das im post? und wenn das im post stimmt wie schauts für die anderen richtungen aus?
danke
TehSackL

[freiBär]

Deins stimmt, das andere ist offensichtlich ein Tippfehler. Der ist mir damals nicht mal aufgefallen...

Gesendet von meinem Nexus S mit Tapatalk

[ReinhardBla]

Bin zwar eigentlich schon fertig, aber eine Sache verwirrt mich da schon ein bissl.

Im Programm ist ja schon eine Projektion implementiert (eh klar, sonst könnt man ja auch kein FTB machen, wie es ja auch bereits implementiert ist). Steht ja auch der Argumentbehandlung:

else if (strcmp(*argv,"-x") == 0) {
            dir = MX;
            ++argv;
            argc--;
        }

Die Projektion die aber ursprünglich im Programm implementiert ist, ist jetzt aber nicht die -x Projektion (wie das in der Argumentbehandlung steht), sondern die -z Projektion oder? Weil er nimmt ja jedes (x,y) Paar und schickt für das einen Strahl in die z-Richtung. Richtig?

[TehSackL]

Seh ich auch so das die innerste Schleife die Richtung angiebt, auch freiBär hats ja in seinem ersten Post geschrieben, ich hät aber gedacht es is die +z Richtung die implementiert is ^^, aber mit den richtungen kämpf ich eben auch noch...

[Gil]

Es wurde heute gesagt, dass die Prüfung am 21. Juni ist. Abgabe des Projekts soll einfach davor per Mail geschehen. Was auch ganz interessant ist, dass es einen Fragenkatalog geben wird. War das letztes Jahr auch schon und gibt es da vielleicht schon eine Ausarbeitung? Ich habe ein GoogleDocs angelegt: https://docs.google.com/document/d/1...mWFqoHdHw/edit falls sich jemand bei der Ausarbeitung beteiligen möchte.

Inzwischen ist der Fragenkatalog auf der LVA-Homepage verlinkt.
Oder direkt: http://www.cg.tuwien.ac.at/courses/M...Topics2012.pdf

[deathhero]

hallo,

ich habe noch eine frage zu isosurface, da mir das noch nicht ganz klar ist.

1. threshold
definiere ich diesen wert selbst? zB einfach 100 ?

2. gx, gy gz ausrechnen sollte machbar sein.
dann multiplikation mit den v werten: (zB +x)
val = gx*1 + gy*0 + gz*0;

was mache ich jetzt mit diesem val wert in der innersten schleife? bzw. was passiert beim nächsten schleifendurchlauf?
zusammenzählen?

bitte um kurzes feedback

danke lg

[ReinhardBla]

1. threshold
definiere ich diesen wert selbst? zB einfach 100 ?

Ja, ich hab aber z.B. das Argumenthandling erweitert, damit man den threshold als Parameter mitgeben kann.

2. gx, gy gz ausrechnen sollte machbar sein.
dann multiplikation mit den v werten: (zB +x)
val = gx*1 + gy*0 + gz*0;

was mache ich jetzt mit diesem val wert in der innersten schleife? bzw. was passiert beim nächsten schleifendurchlauf?
zusammenzählen?

Ich glaub so ganz hab ich das auch noch nicht verstanden, was er da jetzt wirklich will. Ich weiß nur eins, wenn du einen Volume Renderer hast, der ein Isosurface simulieren soll, dann definierst du einen Threshold und wenn der erreicht ist, dann brichst du die Loop fürs Ray-Casting ab und berechnest an dieser Stelle das Shading, siehe z.B. hier: http://code.google.com/p/visualizati...svn1118&r=1118

Nachdem ich noch keine Antwort bezüglich Feedback bekomme habe, hier mal ein paar Bilder von meinem Programm. Ja, so falsch schauts eigentlich gar nicht aus, ist ja nur eine diffuse Beleuchtungssituation bei mir...

[deathhero]

Ich glaub so ganz hab ich das auch noch nicht verstanden, was er da jetzt wirklich will. Ich weiß nur eins, wenn du einen Volume Renderer hast, der ein Isosurface simulieren soll, dann definierst du einen Threshold und wenn der erreicht ist, dann brichst du die Loop fürs Ray-Casting ab und berechnest an dieser Stelle das Shading, siehe z.B. hier: http://code.google.com/p/visualizati...svn1118&r=1118
Nachdem ich noch keine Antwort bezüglich Feedback bekomme habe, hier mal ein paar Bilder von meinem Programm. Ja, so falsch schauts eigentlich gar nicht aus, ist ja nur eine diffuse Beleuchtungssituation bei mir...

ok also du rechnest diesen val wert für jeden punkt nur einmal aus und verlässt danach die innerste schleife, richtig?

[ReinhardBla]

Ja, ich mach es derzeit schon so. Sonst würde meiner Meinung nach der Threshold wenig Sinn machen.

Edit: Ich hatte es aber auch schon anders, weil ich versucht habe ein transparentes Isosurface darzustellen. Bei dem Ansatz hab ich dann aufsummiert und ein alpha berechnet, aber das hat bei den Volumina extrem scheisse ausgesehen.

[deathhero]

ok dann werde ich das auch so machen... thx!

[ReinhardBla]

Ja, post mal ein paar Bilder, wenn du fertig bist

[deathhero]

mache ich....morgen

[Stephe]

Der Programmaufruf sollte in etwa so ausschaun:
./project +x -t 4 data/tot2-n.f3d > pics/ISO/tot2_px_4.ppm
+x ist die Direction, -t 4 ist dann für die Methode.

Wenn ich das so eingeb (oder mit +y, -y etc.) , bekomm ich "Incorrect number of arguments". WTF?

[Hadesskywalker]

Es ist am Anfang nur +x implementiert. Schau in die main-Methode und ergänze einfach die restlichen Richtungen.

[Stephe]

Es ist am Anfang nur +x implementiert. Schau in die main-Methode und ergänze einfach die restlichen Richtungen.

Sorry, ja. Ich dachte, die Argumentbehandlung selbst wurde schon implementiert, was aber nicht der Fall ist Jetzt geht's.

[deathhero]

ich bekomme hin und wieder den fehler: "segmentation fault (core dumped)"
sonst zufällig auch jemand?

[ReinhardBla]

Dann greifst du bei in->get(x, y, z); auf einen Voxel zu, den es nicht gibt...

[deathhero]

alles klar, thx. meine iso bilder schauen übrigens wie deine aus...

[TehSackL]

Hallo,

ich bin mir leider noch unsicher mit den projektionsrichtungen, zb bei den oben geposteten bildern ist der teapot ja verdreht ich durchlaufe momentan die bilder so (von der äußersten zur innersten schleife, vorzeichen gibt die richtung an)
PX: (+y,+z,+x)
MX: (-y,+z,-x)
PY: (+x,+z,+y)
MY: (-x,+z,-y)
PZ: (+x,+y,+z)
MZ: (-x,+y,-z),

Bin mir aber unsicher damit, weil für tot-2 find ich die bilder ok von den richtungen her, aber der teapot steht irgendwie wie bei den bildern oben, die innerste variable ist eigentlich klar, die äußersten verdreh ich damit die bilder für zb mip in +x und -x nicht ident sindern, sondern wirklich so als würde man von der anderen seite draufschauen.

Wie habts ihr das gemacht und weiss vielleicht wer wies der prof gern hätte?

[freiBär]

Dem Professor ist das ziemlich egal. Machs einfach in die Richtungen, egal obs am Kopf steht oder nicht. Das Projekt ist einfach nur Voraussetzung für die Prüfung und zählt nicht zur Note.

[maxputz]

Bei mir war das auch kopfstehend. Dem Prof ist das aber komplett egal. Ich glaub sogar, dass er sich den Source code gar ned anschaut, wenn die bilder passen.

[freiBär]

Der schaut bei jedem ungefähr eine halbe Minute ins PDF, ob die Bilder halbwegs passen....das ist den Professoren sowas von egal....

[maxputz]

Weiß jdm vom Vorjahr wie detailiert die Fragen bei der Prüfung gestellt werden? 4 Fragen mit je 10 Unterpunkten oder sind die mit ein paar Stichwörtern beantwortbar? Schwierigkeit? Benotung?

[ReinhardBla]

Hallo,

ich bin mir leider noch unsicher mit den projektionsrichtungen, zb bei den oben geposteten bildern ist der teapot ja verdreht ich durchlaufe momentan die bilder so (von der äußersten zur innersten schleife, vorzeichen gibt die richtung an)
PX: (+y,+z,+x)
MX: (-y,+z,-x)
PY: (+x,+z,+y)
MY: (-x,+z,-y)
PZ: (+x,+y,+z)
MZ: (-x,+y,-z),

Bin mir aber unsicher damit, weil für tot-2 find ich die bilder ok von den richtungen her, aber der teapot steht irgendwie wie bei den bildern oben, die innerste variable ist eigentlich klar, die äußersten verdreh ich damit die bilder für zb mip in +x und -x nicht ident sindern, sondern wirklich so als würde man von der anderen seite draufschauen.

Wie habts ihr das gemacht und weiss vielleicht wer wies der prof gern hätte?

Nö, das passt schon so. Aber in meinem Post oben fehlt ein Bild! Bei der +y stehts dann nicht am Kopf. Du kannst dir das Ganze ja auch mit dem f3dview ansehen, da steht es auch am Kopf



Betreffend Vertauschen von +z und -z Projektion: Imho kommt es da ja drauf an, welches Koordinatensystem du verwendest, ob rechtshändig oder linkshändig. Und nachdem das nicht angegeben ist, würd ich mal sagen, dass es wurscht ist :P

[freiBär]

Weiß jdm vom Vorjahr wie detailiert die Fragen bei der Prüfung gestellt werden? 4 Fragen mit je 10 Unterpunkten oder sind die mit ein paar Stichwörtern beantwortbar? Schwierigkeit? Benotung?

Naja es waren schon ein paar (3-4) Unterpunkte dabei. Wir hatten aber letztes Jahr keinen Fragenkatalog (zudem war das das erste mal, dass die Prüfung schriftlich war)...vl stellen sie ja die Fragen anders.
Die Fragen waren schon ziemlich ins Detail, also wirklich immer alles hinschreiben was in den Folien stand. Nur Formeln waren nicht nötig (vl eventuell noch die Formel fürs CT). Zudem kam gar nichts von Dicom oder GPUs. Wies mit der Benotung läuft weis ich nicht, ich wusste damals alles...kann also nicht sagen wie sich ein Fehler/Nichtwissen auswirkt.

Gekommen sind (soweit ichs noch weis):
CT oder MRT erklären (abhängig von der Gruppe)
Transferfunktion
Marching Cubes
Perspective Rendering Probleme/Lösung

[Big_Show]

Hallo, ich muss diese Übungsbsp noch abgeben und komme nun bei den Projektionsrichtungen leider nicht mehr weiter. Also bis jetzt verstehe ich, dass die innerste for-Schleife die Projektionsrichtung angbit, somit muss man die for-Schleifen Anordnungen ändern um die Projektionsrichtung zu ändern. Jedoch ist mir nicht klar wo der Unterschied zwischen zB +x und -x ist. Was muss man anders im Code machen? Eventuell kann mir jemand einen kleinen Tipp geben. Danke sehr.

[freiBär]

Das ist einfach von den entgegengesetzten Seiten. Lauf die for Schleife einfach in der anderen Richtung durch.

Gesendet von meinem Nexus S mit Tapatalk 2

[Big_Show]

Ich hätte noch eine Frage zur Angabe: "Select 6 images you like representing each category, submit them in a pdf file"
Müssen wir zusätzlich 3 ( + 3 die in der Angabe vorhanden sind) .f3d Bilder irgendwo raussuchen und mit diesen arbeiten?
Danke sehr.

[freiBär]

Früher wars zumindest nicht so.

Gesendet von meinem Nexus 7 mit Tapatalk 2

[maxputz]

Ich hätte noch eine Frage zur Angabe: "Select 6 images you like representing each category, submit them in a pdf file"
Müssen wir zusätzlich 3 ( + 3 die in der Angabe vorhanden sind) .f3d Bilder irgendwo raussuchen und mit diesen arbeiten?
Danke sehr.

Du hast pro Kategorie ("Funktion") drei Bilder (bigbrain, teapot, tot2) mit jeweils 6 Richtungen. Das sind 18 Kombinationen pro Kategorie. Da suchst du dir jeweils 6 raus, die dir gut gefallen. Das heißt insgesamt mindestens 30 Bilder (5x6) von den 90 Möglichen (5x18).

[Big_Show]

Vielen Dank, sehr nett