Ein bisschen lineare Algebra

[LDericher]

Also, eigentlich dürfte das kein Problem sein - habe ich mir gedacht...

Ich wollte mir einen Farb-Ring bauen; später will ich daraus auch ein Control machen. Was das Teil genau können soll, müsst ihr noch nicht wissen, denn es hakt schon deutlich vorher.

Ich will die Farbe jedes Pixels über den Winkel errechnen, und wenn ich den kennen würde, wäre das ganze auch kein Problem. Folgendermaßen gehe ich bis jetzt vor:

1. Ich zeichne einen schwarzen Ring auf weißem Grund. (Canvas.Ellipse, temporäres BMP)
2. Ich durchlaufe die einzelnen Pixel vom Bild (im Moment noch über Canvas.Pixels, aber das werde ich noch rausoptimieren) und schreibe auf das Zielcanvas (das von der Form) den berechneten Farbwert, sofern ich kein weißes Pixel habe.

Den Farbwert will ich berechnen, indem ich den Winkel "alpha" verarbeite (siehe angefügte Skizze).
Natürlich sind N, Z und M bekannt.

Meine Überlegung ist nun so:

Code:

N(aRadius, 0).
M(aRadius, aRadius).
Z(x, y).
n = N - M = (0, -aRadius).
z = Z - M = (x - aRadius, y - aRadius).
alpha = arccos( (n * z) / (|n| * |z|) ).

wobei n * z für das Skalarprodukt und |n| für den Betrag/die Länge steht.

In Quellcode habe ich das ganze nun so versucht umzusetzen:

Delphi-Quellcode:

function ScalarProduct(var pVec1, pVec2:TPoint):integer;
begin
Result:=(pVec1.X * pVec2.X) + (pVec1.Y * pVec2.Y);
end;

function Length(var pVec:TPoint):Extended;
begin
Result:=Sqrt(Sqr(pVec.X)+Sqr(pVec.Y));
end;
{...}
function TColorRing.getColor(X, Y: integer): TColor;
var
  vecN,
  vecZ:TPoint;
  Angle:Extended;
begin
vecN:=Point(0, -aRadius);
vecZ:=Point(X-aRadius, Y-aRadius);
Angle:=ScalarProduct(vecN, vecZ)/(Length(vecN)*Length(vecZ));
Angle:=ArcCos(Angle);
Result:=trunc(Angle*(360/Pi)); // Sollte Werte zwischen 0 und 360 geben.
end;

Aber leider scheint das ganze nicht so funktionieren zu wollen, wie ich das gern hätte, denn Result liegt zwar richtigerweise zwischen 0 und 360, aber die Werte sind "gespiegelt", also der rechte Halbkreis liefert 0..360 genau wie der Linke. Ich hätte aber gern oben die 0/360, rechts die 90, unten die 180 und links die 270; der Rest entsprechend dazwischen...

Wer kann mir helfen, meinen Fehler zu finden?

Gruß,
Euer LDer

[gammatester]

Pi entspricht 180° nicht 360°!

[LDericher]

Zitat:

Zitat von gammatester

Pi entspricht 180° nicht 360°!

Wah!
Okay, danke, jetzt hab ich schonmal die 180 unten, das ist schonmal gut, danke :)
Aber jetzt ist immer noch die 90 links UND rechts, aber die soll nur rechts sein, links muss die 270 hin. Ein schöner Kreis solls werden, nicht zwei Halbkreise... irgendwo hakts da noch :cry:

[Medium]

Das Problem ist, dass diese Art der Winkelbestimmung immer zum kleineren Winkel zwischen 2 Vektoren führt. Schau dir mal ATan2() an!

[jfheins]

Kennst du Polarkooradinaten?

Um den Winkel alpha zu ermitteln: alpha = Math.arctan2(z.Y, z.X)

Der Winkel ist allerdings mathematisch - also gleich null wenn der Punkt ganz rechts liegt und dann im Uhrzeigersinn positiv (da ja die y-Achse nach unten geht)

[gammatester]

arccos bildet [-1..1] auf [0..Pi] ab, also wirst Du niemals Winkel größer als 180° erhalten. Dies Situation muß Deine übergeordnete Logik behandeln. Vielleicht so (ohne Gewähr): wenn Winkel(n,z) > 180, dann 360 - Winkel(z,n), genaueres mußt Du eventuell in Deiner LA-Quelle nachlesen.

[gammatester]

Zitat:

Zitat von jfheins

Kennst du Polarkooradinaten?

Um den Winkel alpha zu ermitteln: alpha = Math.arctan2(z.Y, z.X)

Das ist aber nicht alpha, sondern der Winkel zwischen x-Achse und z!

[jfheins]

Zitat:

Zitat von gammatester

Zitat:

Zitat von jfheins

Kennst du Polarkooradinaten?

Um den Winkel alpha zu ermitteln: alpha = Math.arctan2(z.Y, z.X)

Das ist aber nicht alpha, sondern der Winkel zwischen x-Achse und z!

Ich hab ja auch extra dazugeschriben, dass der Winkel dann null ist, wenn der Punkt waagerecht rechts des Mittelpunkts liegt. Einen konstanten Winkel zu addieren sollte noch machbar sein ...

Aber, okay, hier die korrigierte Formel:

alpha = Math.arctan2(z.Y, z.X) + PI/2; // bzw. 90°
if alpha < 0 then alpha = alpha + 2*PI; // Danach ist alpha nicht mehr negativ

[Medium]

atan2()-pi/2 sollte das fixen, wobei man dann natürlich ein Auge auf alpha<0 haben muss. Falls links/rechts dann noch vertauscht sind, den 2. Parameter negieren.

Edit: Negative Winkel kommen eh vor, stimmt ja. Die muss man also eh schon zurecht zuppeln :)

[LDericher]

Das mit [0..Pi] ist mir dann auch aufgefallen... aber der arctan2 ist genau das was ich gebraucht habe! :)

Also hab ich was gebastelt, das genau das tut, was ich will:

Delphi-Quellcode:

function TColorRing.getAngle(X, Y: integer): integer;
var
  vecZ:TPoint;
  Angle:Extended;
begin
vecZ:=Point(X-aRadius, Y-aRadius);
Angle:=Math.ArcTan2(vecZ.X, vecZ.Y);
Result:=trunc(Angle*(180/Pi));
Result:=180-Result;
end;

Danke für eure Mühen!

Euer LDer