Changeset 25 for trunk/UGeometric.pas

Timestamp:

Apr 12, 2015, 12:08:05 AM (9 years ago)

Author:

chronos

Message:

• Added: Draw track lines which go in same position side by side to see them all and be able to grab one of them. This is not yet perfectly implemented as this require to implement track points link properties first.
• Fixed: Do not allow connect end station to grabbed track between two stations.
• Fixed: Index on wrong list in shape distance calculation function.

File:

• trunk/UGeometric.pas (modified) (4 diffs)

trunk/UGeometric.pas

@@ -13 +13 @@
 function Distance(P1, P2: TPoint): Integer;
 function Dot(const P1, P2: TPoint): Double;
+function AddPoint(const P1, P2: TPoint): TPoint;
 function SubPoint(const P1, P2: TPoint): TPoint;
 function PointToLineDistance(const P, V, W: TPoint): Integer;
@@ -18 +19 @@
 function RotatePoint(Center, P: TPoint; Angle: Double): TPoint;
 function RotatePoints(Center: TPoint; P: TPointArray; Angle: Double): TPointArray;
+function LineIntersect(LineAP1, LineAP2, LineBP1, LineBP2: TPoint): TPoint;
+function ArcTan2Point(Point: TPoint): Float;
+function ArcTanPoint(Point: TPoint): Float;
 
 implementation
@@ -29 +33 @@
 begin
   Result := P1.X * P2.X + P1.Y * P2.Y;
+end;
+
+function AddPoint(const P1, P2: TPoint): TPoint;
+begin
+  Result.X := P1.X + P2.X;
+  Result.Y := P1.Y + P2.Y;
 end;
 
@@ -88 +98 @@
 end;
 
+function LineIntersect(LineAP1, LineAP2, LineBP1, LineBP2: TPoint): TPoint;
+Var
+  LDetLineA, LDetLineB, LDetDivInv: Double;
+  LDiffLA, LDiffLB: TPoint;
+begin
+  LDetLineA := LineAP1.X * LineAP2.Y - LineAP1.Y * LineAP2.X;
+  LDetLineB := LineBP1.X * LineBP2.Y - LineBP1.Y * LineBP2.X;
+
+  LDiffLA := SubPoint(LineAP1, LineAP2);
+  LDiffLB := SubPoint(LineBP1, LineBP2);
+
+  LDetDivInv := 1 / ((LDiffLA.X * LDiffLB.Y) - (LDiffLA.Y * LDiffLB.X));
+
+  Result.X := Trunc(((LDetLineA * LDiffLB.X) - (LDiffLA.X * LDetLineB)) * LDetDivInv);
+  Result.Y := Trunc(((LDetLineA * LDiffLB.Y) - (LDiffLA.Y * LDetLineB)) * LDetDivInv);
+end;
+
+function ArcTan2Point(Point: TPoint): Float;
+begin
+  Result := ArcTan2(Point.Y, Point.X);
+end;
+
+function ArcTanPoint(Point: TPoint): Float;
+begin
+  Result := ArcTan(Point.Y / Point.X);
+end;
 
 end.