|
1
|
|
|
<?php |
|
2
|
|
|
|
|
3
|
|
|
declare(strict_types=1); |
|
4
|
|
|
|
|
5
|
|
|
namespace Location\Processor\Polyline; |
|
6
|
|
|
|
|
7
|
|
|
use Location\Bearing\BearingEllipsoidal; |
|
8
|
|
|
use Location\GeometryInterface; |
|
9
|
|
|
use Location\Polygon; |
|
10
|
|
|
use Location\Polyline; |
|
11
|
|
|
use RuntimeException; |
|
12
|
|
|
|
|
13
|
|
|
/** |
|
14
|
|
|
* Simplify Polyline. |
|
15
|
|
|
* |
|
16
|
|
|
* @author Marcus Jaschen <[email protected]> |
|
17
|
|
|
*/ |
|
18
|
|
|
class SimplifyBearing implements SimplifyInterface |
|
19
|
|
|
{ |
|
20
|
|
|
/** |
|
21
|
|
|
* @var float |
|
22
|
|
|
*/ |
|
23
|
|
|
private $bearingAngle; |
|
24
|
|
|
|
|
25
|
|
|
/** |
|
26
|
|
|
* SimplifyBearing constructor. |
|
27
|
|
|
* |
|
28
|
|
|
* @param float $bearingAngle |
|
29
|
|
|
*/ |
|
30
|
|
|
public function __construct(float $bearingAngle) |
|
31
|
|
|
{ |
|
32
|
|
|
$this->bearingAngle = $bearingAngle; |
|
33
|
|
|
} |
|
34
|
|
|
|
|
35
|
|
|
/** |
|
36
|
|
|
* @param Polyline $polyline |
|
37
|
|
|
* |
|
38
|
|
|
* @return Polyline |
|
39
|
|
|
* @throws RuntimeException |
|
40
|
|
|
*/ |
|
41
|
|
View Code Duplication |
public function simplify(Polyline $polyline): Polyline |
|
|
|
|
|
|
42
|
|
|
{ |
|
43
|
|
|
$result = $this->simplifyGeometry($polyline); |
|
44
|
|
|
|
|
45
|
|
|
if (!($result instanceof Polyline)) { |
|
46
|
|
|
throw new RuntimeException('Result is no Polyline', 4231694400); |
|
47
|
|
|
} |
|
48
|
|
|
|
|
49
|
|
|
return $result; |
|
50
|
|
|
} |
|
51
|
|
|
|
|
52
|
|
|
/** |
|
53
|
|
|
* Simplifies the given polyline |
|
54
|
|
|
* |
|
55
|
|
|
* 1. calculate the bearing angle between the first two points p1 and p2: b1 |
|
56
|
|
|
* 2. calculate the bearing angle between the next two points p2 and p3: b2 |
|
57
|
|
|
* 3. calculate the difference between b1 and b2: deltaB; if deltaB is |
|
58
|
|
|
* smaller than the threshold angle, remove the middle point p2 |
|
59
|
|
|
* 4. start again at (1.) as long as the polyline contains more points |
|
60
|
|
|
* |
|
61
|
|
|
* This method will be merged with `simplify()` in the next major release. |
|
62
|
|
|
* |
|
63
|
|
|
* @param GeometryInterface $geometry |
|
64
|
|
|
* |
|
65
|
|
|
* @return GeometryInterface |
|
66
|
|
|
*/ |
|
67
|
|
|
public function simplifyGeometry(GeometryInterface $geometry): GeometryInterface |
|
68
|
|
|
{ |
|
69
|
|
|
if (!($geometry instanceof Polyline) && !($geometry instanceof Polygon)) { |
|
70
|
|
|
return $geometry; |
|
71
|
|
|
} |
|
72
|
|
|
|
|
73
|
|
|
$counterPoints = $geometry->getNumberOfPoints(); |
|
74
|
|
|
|
|
75
|
|
View Code Duplication |
if ($geometry instanceof Polygon) { |
|
|
|
|
|
|
76
|
|
|
if ($counterPoints <= 3) { |
|
77
|
|
|
return clone $geometry; |
|
78
|
|
|
} |
|
79
|
|
|
$result = new Polygon(); |
|
80
|
|
|
} else { |
|
81
|
|
|
if ($counterPoints < 3) { |
|
82
|
|
|
return clone $geometry; |
|
83
|
|
|
} |
|
84
|
|
|
$result = new Polyline(); |
|
85
|
|
|
} |
|
86
|
|
|
|
|
87
|
|
|
$bearingCalc = new BearingEllipsoidal(); |
|
88
|
|
|
|
|
89
|
|
|
$points = $geometry->getPoints(); |
|
90
|
|
|
|
|
91
|
|
|
$index = 0; |
|
92
|
|
|
|
|
93
|
|
|
// add the first point to the resulting polyline |
|
94
|
|
|
$result->addPoint($points[$index]); |
|
95
|
|
|
|
|
96
|
|
|
do { |
|
97
|
|
|
$index++; |
|
98
|
|
|
|
|
99
|
|
|
// preserve the last point of the original polyline |
|
100
|
|
|
if ($index === ($counterPoints - 1)) { |
|
101
|
|
|
$result->addPoint($points[$index]); |
|
102
|
|
|
break; |
|
103
|
|
|
} |
|
104
|
|
|
|
|
105
|
|
|
$bearing1 = $bearingCalc->calculateBearing($points[$index - 1], $points[$index]); |
|
106
|
|
|
$bearing2 = $bearingCalc->calculateBearing($points[$index], $points[$index + 1]); |
|
107
|
|
|
|
|
108
|
|
|
$bearingDifference = min( |
|
109
|
|
|
fmod($bearing1 - $bearing2 + 360, 360), |
|
110
|
|
|
fmod($bearing2 - $bearing1 + 360, 360) |
|
111
|
|
|
); |
|
112
|
|
|
|
|
113
|
|
|
if ($bearingDifference > $this->bearingAngle) { |
|
114
|
|
|
$result->addPoint($points[$index]); |
|
115
|
|
|
} |
|
116
|
|
|
} while ($index < $counterPoints); |
|
117
|
|
|
|
|
118
|
|
|
return $result; |
|
119
|
|
|
} |
|
120
|
|
|
} |
|
121
|
|
|
|
mjaschen /
phpgeo
Loading history...
Duplicated code is one of the most pungent code smells. If you need to duplicate the same code in three or more different places, we strongly encourage you to look into extracting the code into a single class or operation.
You can also find more detailed suggestions in the “Code” section of your repository.