2023-05-28 19:14:02 +03:00
|
|
|
package gx
|
|
|
|
|
|
|
|
import (
|
|
|
|
"math"
|
|
|
|
)
|
|
|
|
|
2023-05-28 19:45:58 +03:00
|
|
|
// The type represents mathematical equation of line and line itself.
|
|
|
|
type Line struct {
|
2023-05-28 19:14:02 +03:00
|
|
|
K, C Float
|
|
|
|
}
|
|
|
|
|
2023-05-28 20:02:07 +03:00
|
|
|
type Liner interface {
|
|
|
|
Line() Line
|
|
|
|
}
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
type LinerPointContainer interface {
|
|
|
|
Liner
|
|
|
|
PointContainer
|
|
|
|
}
|
|
|
|
|
2023-05-28 19:14:02 +03:00
|
|
|
// The type represents a line segment.
|
|
|
|
type LineSegment [2]Point
|
|
|
|
|
2023-06-03 11:25:19 +03:00
|
|
|
// The type represents multiple line segments.
|
|
|
|
type LineSegments []LineSegment
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
type Edge = LineSegment
|
|
|
|
type Edges []Vertex
|
|
|
|
|
|
|
|
|
|
|
|
// Check if two LinerPointContainers do cross and return the
|
|
|
|
// crossing point.
|
|
|
|
func LinersCross(lp1, lp2 LinerPointContainer) (Point, bool) {
|
|
|
|
l1 := lp1.Line()
|
|
|
|
l2 := lp2.Line()
|
2023-05-28 20:02:07 +03:00
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
p, crosses := l1.crossesLine(l2)
|
|
|
|
if !crosses ||
|
|
|
|
!lp1.ContainsPoint(p) ||
|
|
|
|
!lp2.ContainsPoint(p) {
|
|
|
|
return Point{}, false
|
2023-05-28 20:02:07 +03:00
|
|
|
}
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
return p, true
|
|
|
|
}
|
|
|
|
|
|
|
|
// Check whether the liner is parallel to the other liner.
|
|
|
|
func LinersParallel(first, second Liner) bool {
|
|
|
|
l1 := first.Line()
|
|
|
|
l2 := second.Line()
|
|
|
|
|
|
|
|
return l1.K == l2.K
|
2023-05-28 20:02:07 +03:00
|
|
|
}
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
// Returns the line itself. Made to implement the Liner interface.
|
2023-05-28 21:07:05 +03:00
|
|
|
func (l Line) Line() Line {
|
|
|
|
return l
|
|
|
|
}
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
// Returns corresponding to the segment line line.
|
|
|
|
func (l LineSegment) Line() Line {
|
|
|
|
p0 := l[0]
|
|
|
|
p1 := l[1]
|
|
|
|
|
|
|
|
k := (p0.Y - p1.Y) / (p0.X - p1.X)
|
|
|
|
c := p0.Y - p0.X*k
|
|
|
|
|
|
|
|
return Line{k, c}
|
|
|
|
}
|
|
|
|
|
2023-05-28 21:22:52 +03:00
|
|
|
func (l Line) ContainsPoint(p Point) bool {
|
|
|
|
buf := Line{0, p.Y}
|
|
|
|
pc, ok := l.crossesLine(buf)
|
|
|
|
if !ok {
|
|
|
|
return false
|
2023-05-28 21:07:05 +03:00
|
|
|
}
|
|
|
|
|
2023-05-28 21:22:52 +03:00
|
|
|
return pc == p
|
|
|
|
}
|
|
|
|
|
2023-06-03 12:26:31 +03:00
|
|
|
func (l LineSegment) ContainsPoint(p Point) bool {
|
|
|
|
line := l.Line()
|
|
|
|
if !line.ContainsPoint(p) {
|
|
|
|
return false
|
2023-05-28 21:07:05 +03:00
|
|
|
}
|
2023-06-03 12:26:31 +03:00
|
|
|
|
|
|
|
xMax := Max(l[0].X, l[1].X)
|
|
|
|
xMin := Min(l[0].X, l[1].X)
|
|
|
|
|
|
|
|
yMax := Max(l[0].Y, l[1].Y)
|
|
|
|
yMin := Min(l[0].Y, l[1].Y)
|
|
|
|
|
|
|
|
if !(xMin < p.X && p.X < xMax) ||
|
|
|
|
!(yMin < p.Y && p.Y < yMax) {
|
|
|
|
return false
|
|
|
|
}
|
|
|
|
|
|
|
|
return true
|
2023-05-28 21:07:05 +03:00
|
|
|
}
|
|
|
|
|
|
|
|
func (l1 Line) crossesLine(l2 Line) (Point, bool) {
|
2023-06-03 12:26:31 +03:00
|
|
|
if LinersParallel(l1, l2) {
|
2023-05-28 21:22:52 +03:00
|
|
|
return Point{}, false
|
|
|
|
}
|
|
|
|
|
2023-05-28 21:07:05 +03:00
|
|
|
x := (l1.C - l2.C) / (l2.K - l1.K)
|
|
|
|
y := l1.K*x + l1.C
|
|
|
|
return Point{x, y}, true
|
|
|
|
}
|
|
|
|
|
2023-05-30 13:12:11 +03:00
|
|
|
|
2023-05-28 19:14:02 +03:00
|
|
|
// Get square of length of line segment.
|
|
|
|
func (ls LineSegment) LenSqr() Float {
|
|
|
|
return Sqr(ls[0].X - ls[1].X) +
|
|
|
|
Sqr(ls[0].Y - ls[1].Y)
|
|
|
|
}
|
|
|
|
|
|
|
|
// Get length of the line segment.
|
|
|
|
func (ls LineSegment) Len() Float {
|
|
|
|
return math.Sqrt(ls.LenSqr())
|
|
|
|
}
|
|
|
|
|