PureBoard

Global Dim Bild.i(1000,1000)
Global maxs.d,mins.d


UseJPEGImageDecoder()

Procedure.d GetMax(a.d, b.d, c.d)
Protected Result.d

Result = a
If b > Result
Result = b
EndIf
If c > Result
Result = c
EndIf

ProcedureReturn Result
EndProcedure

Procedure.d GetMin(a.d, b.d, c.d)
Protected ...

Hast du mal ein Bild, wo man die Vignettierung deutlicher sehen kann als beim Bild vom Lilienstein

Ich habe mal etwas experimentiert, vielleicht kannst Du etwas damit anfangen:
Global Dim Bild.i(1000,1000)
Global maxs.d,mins.d


UsePNGImageDecoder()

Procedure.d GetMax(a.d, b.d, c.d)
Protected Result.d

Result = a
If b > Result
Result = b
EndIf
If c > Result
Result = c
EndIf ...

Noch ein paar (ungeordnete) Bemerkungen:
1. alle genannten Splinetypen lassen sich ineinander umrechnen.
2. bei rationalen (B-)Splines haben alle Kontrollpunkte ein "Gewicht" - je höher das Gewicht, desto näher wird die "Kurve" an den Kontrollpunkt "gezogen"
3. mit rationalen Splines lassen sich ...

Also da ich nur die Geometrie und Darstellungsattribute aus dem STEP einlesen wollte, habe ich mich auf folgende Elemente konzentriert
CARTESIAN_POINT
VERTEX_POINT
DIRECTION
AXIS2_PLACEMENT_3D
AXIS2_PLACEMENT_2D
AXIS1_PLACEMENT
VECTOR
BOUNDED_CURVE
LINE
CIRCLE
ELLIPSE
POLYLINE
POLY_LOOP ...

ich habe für ein CAM-System eine STEP-Import und -Exportschnittstelle gemacht, allerdings in C. Hilft dir also nicht so viel. Bei den ganzen mathematischen Sachen (B-Splines, Triangulation von Flächen etc.) kann ich sicher helfen.

Das ist ein normales Matrix Problem:

Du berechnest eine Richtung R mit einer Darstellungsmatrix m ( aus sin() und cos() der Drehwinkel ) aus einer Ausgangsrichtung A über :
A*m=R
Ax
Ay
Az
m11 m12 m13 Rx
m21 m22 m23 Ry
m31 m32 m33 Rz

um den Ausgangsvektor wieder herzustellen musst du die ...

Dann musst du ja deine 3D Darstellung irgendwie auf die 2D Darstellung umrechnen. Der Codeabschnitt wäre dann interessant.

Den Quelltext kann ich mir nicht anschauen, da mein Browser die Datei nicht runterladen mag.
Schau mal in der Hilfe unter PointPick() nach, da kann man die Richtung mit PickX() etc. auslesen.

hier mein Testcode:

EnableExplicit

Procedure.i PunktInKonvexemPolyeder ( Anzahl.i, Array Indizes.i(2), Array Punkte.f(2), Px.f, Py.f, Pz.f)
Protected i.i
Protected Vx1.d,Vx2.d,Vy1.d,Vy2.d,Vz1.d,Vz2.d,Nx.d,Ny.d,Nz.d,L.d

For i=0 To Anzahl-1
Vx1 = Punkte(Indizes(i,1),0)-Punkte(Indizes(i,0),0 ...