1
|
|
|
<?php |
2
|
|
|
|
3
|
|
|
namespace Smindel\GIS; |
4
|
|
|
|
5
|
|
|
use SilverStripe\Core\Config\Config; |
6
|
|
|
use SilverStripe\Core\Config\Configurable; |
7
|
|
|
use SilverStripe\Core\Injector\Injectable; |
8
|
|
|
use SilverStripe\ORM\DB; |
9
|
|
|
use Smindel\GIS\ORM\FieldType\DBGeography; |
10
|
|
|
use Smindel\GIS\ORM\FieldType\DBGeometry; |
11
|
|
|
use proj4php\Proj4php; |
12
|
|
|
use proj4php\Proj; |
13
|
|
|
use proj4php\Point; |
14
|
|
|
use Exception; |
15
|
|
|
|
16
|
|
|
/** |
17
|
|
|
* @property string $array |
18
|
|
|
* @property string $ewkt |
19
|
|
|
* @property string $wkt |
20
|
|
|
* @property string $srid |
21
|
|
|
* @property string $type |
22
|
|
|
* @property string $coordinates |
23
|
|
|
*/ |
24
|
|
|
class GIS |
25
|
|
|
{ |
26
|
|
|
use Configurable; |
27
|
|
|
|
28
|
|
|
use Injectable; |
29
|
|
|
|
30
|
|
|
const WKT_PATTERN = '/^(([A-Z]+)\s*(\(.+\)))$/i'; |
31
|
|
|
const EWKT_PATTERN = '/^SRID=(\d+);(([A-Z]+)\s*(\(.+\)))$/i'; |
32
|
|
|
|
33
|
|
|
const TYPES = [ |
34
|
|
|
'point' => 'Point', |
35
|
|
|
'linestring' => 'LineString', |
36
|
|
|
'polygon' => 'Polygon', |
37
|
|
|
'multipoint' => 'MultiPoint', |
38
|
|
|
'multilinestring' => 'MultiLineString', |
39
|
|
|
'multipolygon' => 'MultiPolygon', |
40
|
|
|
'geometrycollection' => 'GeometryCollection' |
41
|
|
|
]; |
42
|
|
|
|
43
|
|
|
private static $default_srid = 4326; |
44
|
|
|
|
45
|
|
|
protected $value; |
46
|
|
|
|
47
|
|
|
/** |
48
|
|
|
* Constructor |
49
|
|
|
* |
50
|
|
|
* @param $value mixed geo value: |
51
|
|
|
* string: wkt (default srid) or ewkt |
52
|
|
|
* array: just coordinates (default srid, autodetect shape) or assoc with type, srid and coordinates |
53
|
|
|
* GIS: extract value |
54
|
|
|
* DBGeography: extract value |
55
|
28 |
|
*/ |
56
|
|
|
public function __construct($value) |
57
|
28 |
|
{ |
58
|
3 |
|
DB::get_schema()->initialise(); |
|
|
|
|
59
|
28 |
|
if ($value instanceof GIS) { |
60
|
|
|
$this->value = $value->value; |
61
|
28 |
|
} else if ($value instanceof DBGeography) { |
62
|
|
|
$this->value = $value->getValue(); |
63
|
|
|
} else if (is_string($value) && preg_match(self::WKT_PATTERN, $value)) { |
64
|
28 |
|
$this->value = 'SRID=' . self::config()->default_srid . ';' . $value; |
65
|
25 |
|
} else if (is_array($value) && count($value) == 3 && isset($value['type'])) { |
66
|
|
|
$this->value = $value; |
67
|
26 |
|
} else if (is_string($value) && preg_match(self::EWKT_PATTERN, $value)) { |
68
|
8 |
|
$this->value = $value; |
69
|
|
|
} else if (empty($value) || isset($value['coordinates']) && empty($value['coordinates'])) { |
70
|
8 |
|
$this->value = null; |
71
|
|
|
} else if (is_array($value) && !isset($value['type'])) { |
72
|
8 |
|
switch (true) { |
73
|
1 |
|
case is_numeric($value[0]): |
74
|
1 |
|
$type = 'Point'; |
75
|
1 |
|
break; |
76
|
|
|
case is_numeric($value[0][0]): |
77
|
8 |
|
$type = 'LineString'; |
78
|
8 |
|
break; |
79
|
8 |
|
case is_numeric($value[0][0][0]): |
80
|
8 |
|
$type = 'Polygon'; |
81
|
|
|
break; |
82
|
|
|
case is_numeric($value[0][0][0][0]): |
83
|
|
|
$type = 'MultiPolygon'; |
84
|
|
|
break; |
85
|
28 |
|
} |
86
|
|
|
$this->value = [ |
87
|
|
|
'srid' => GIS::config()->default_srid, |
88
|
|
|
'type' => $type, |
|
|
|
|
89
|
|
|
'coordinates' => $value, |
90
|
|
|
]; |
91
|
|
|
} else { |
92
|
28 |
|
throw new Exception('Invalid geo value: "' . var_export($value) . '"'); |
|
|
|
|
93
|
|
|
} |
94
|
28 |
|
} |
95
|
13 |
|
|
96
|
|
|
public function __isset($property) |
97
|
|
|
{ |
98
|
|
|
return array_search($property, ['array', 'ewkt', 'srid', 'type', 'coordinates']) !== false; |
99
|
26 |
|
} |
100
|
26 |
|
|
101
|
23 |
|
public function __get($property) |
102
|
22 |
|
{ |
103
|
16 |
|
if (isset($this->value[$property])) { |
104
|
15 |
|
return $this->value[$property]; |
105
|
15 |
|
} |
106
|
|
|
|
107
|
15 |
|
switch ($property) { |
108
|
|
|
case 'array': |
109
|
15 |
|
return ['srid' => $this->srid, 'type' => $this->type, 'coordinates' => $this->coordinates]; |
110
|
9 |
|
case 'ewkt': |
111
|
|
|
return (string)$this; |
112
|
|
|
case 'wkt': |
113
|
15 |
|
return explode(';', (string)$this)[1]; |
114
|
|
|
case 'srid': |
115
|
|
|
return preg_match('/^SRID=(\d+);/i', $this->value, $matches) ? (int)$matches[1] : null; |
116
|
|
|
case 'type': |
117
|
|
|
return preg_match('/^SRID=\d+;(' . implode('|', array_change_key_case(self::TYPES, CASE_UPPER)) . ')/i', $this->value, $matches) ? self::TYPES[strtolower($matches[1])] : null; |
118
|
|
|
case 'coordinates': |
119
|
|
|
if (preg_match(self::EWKT_PATTERN, $this->value, $matches)) { |
120
|
|
|
|
121
|
19 |
|
$coords = str_replace(['(', ')'], ['[', ']'], preg_replace('/([\d\.-]+)\s+([\d\.-]+)/', "[$1,$2]", $matches[4])); |
122
|
|
|
|
123
|
19 |
|
if (strtolower($matches[3]) != 'point') { |
124
|
11 |
|
$coords = "[$coords]"; |
125
|
|
|
} |
126
|
|
|
|
127
|
13 |
|
return json_decode($coords, true)[0]; |
128
|
13 |
|
} else { |
129
|
13 |
|
return null; |
130
|
|
|
} |
131
|
|
|
case 'geometries': |
132
|
13 |
|
// primarily used for GeometryCollections |
133
|
|
|
// @todo: what's supposed to be returned for non-GeometryCollections? |
134
|
|
|
if (preg_match(self::EWKT_PATTERN, $this->value, $matches)) { |
135
|
|
|
$geometries = preg_split('/,(?=[a-zA-Z])/', substr($matches[4], 1, -1)); |
136
|
|
|
$srid = $this->srid; |
137
|
|
|
return array_map(function ($geometry) use ($srid) { |
138
|
|
|
return GIS::create('SRID=' . $srid . ';' . $geometry); |
139
|
|
|
}, $geometries); |
140
|
|
|
} |
141
|
|
|
return null; |
142
|
|
|
default: |
143
|
13 |
|
throw new Exception('Unkown property ' . $property); |
144
|
|
|
} |
145
|
13 |
|
} |
146
|
|
|
|
147
|
|
|
public function __toString() |
148
|
9 |
|
{ |
149
|
|
|
if (is_string($this->value)) { |
150
|
9 |
|
return $this->value; |
151
|
|
|
} |
152
|
|
|
|
153
|
9 |
|
$type = isset($this->value['type']) ? strtoupper($this->value['type']) : null; |
154
|
|
|
$srid = isset($this->value['srid']) ? $this->value['srid'] : GIS::config()->default_srid; |
155
|
9 |
|
$array = isset($this->value['coordinates']) ? $this->value['coordinates'] : $this->value; |
156
|
|
|
|
157
|
|
|
$replacements = [ |
158
|
|
|
'/(?<=\d),(?=-|\d)/' => ' ', |
159
|
9 |
|
'/\[\[\[\[/' => '(((', |
160
|
8 |
|
'/\]\]\]\]/' => ')))', |
161
|
8 |
|
'/\[\[\[/' => '((', |
162
|
|
|
'/\]\]\]/' => '))', |
163
|
|
|
'/\[\[/' => '(', |
164
|
1 |
|
'/\]\]/' => ')', |
165
|
|
|
'/\[/' => '', |
166
|
|
|
'/\]/' => '', |
167
|
|
|
]; |
168
|
|
|
|
169
|
12 |
|
$coords = preg_replace(array_keys($replacements), array_values($replacements), json_encode($array)); |
170
|
|
|
|
171
|
12 |
|
return sprintf('SRID=%d;%s%s', $srid, $type, $type == 'POINT' ? "($coords)" : $coords); |
172
|
12 |
|
} |
173
|
12 |
|
|
174
|
|
|
public function isNull() |
175
|
12 |
|
{ |
176
|
1 |
|
return empty($this->value) || isset($this->value['coordinates']) && empty($this->value['coordinates']); |
177
|
1 |
|
} |
178
|
1 |
|
|
179
|
|
|
public static function of($dataObjectClass) |
180
|
11 |
|
{ |
181
|
|
|
if ($field = $dataObjectClass::config()->get('default_geo_field')) { |
182
|
|
|
return $field; |
183
|
12 |
|
} |
184
|
12 |
|
|
185
|
12 |
|
foreach ($dataObjectClass::config()->get('db') ?: [] as $field => $type) { |
186
|
12 |
|
if (in_array($type, ['Geography', 'Geometry', DBGeography::class, DBGeometry::class])) { |
187
|
|
|
return $field; |
188
|
|
|
} |
189
|
|
|
} |
190
|
|
|
} |
191
|
|
|
|
192
|
|
|
/** |
193
|
|
|
* reproject an array representation of a geometry to the given srid |
194
|
|
|
*/ |
195
|
1 |
|
public function reproject($toSrid = 4326) |
196
|
|
|
{ |
197
|
1 |
|
$fromSrid = $this->srid; |
198
|
|
|
|
199
|
1 |
|
if ($fromSrid == $toSrid) { |
200
|
|
|
return clone $this; |
201
|
1 |
|
} |
202
|
|
|
|
203
|
1 |
|
$fromCoordinates = $this->coordinates; |
204
|
|
|
$type = $this->type; |
205
|
|
|
|
206
|
|
|
$fromProj = self::get_proj4($fromSrid); |
207
|
1 |
|
$toProj = self::get_proj4($toSrid); |
208
|
|
|
$toCoordinates = self::reproject_array($fromCoordinates, $fromProj, $toProj); |
209
|
|
|
|
210
|
1 |
|
return GIS::create([ |
211
|
|
|
'srid' => $toSrid, |
212
|
|
|
'type' => $type, |
213
|
|
|
'coordinates' => $toCoordinates, |
214
|
|
|
]); |
215
|
1 |
|
} |
216
|
1 |
|
|
217
|
1 |
|
/** |
218
|
|
|
* @var proj4php instance |
219
|
|
|
*/ |
220
|
5 |
|
protected static $proj4; |
221
|
|
|
|
222
|
5 |
|
protected static function get_proj4($srid) |
223
|
4 |
|
{ |
224
|
|
|
self::$proj4 = self::$proj4 ?: new Proj4php(); |
225
|
|
|
|
226
|
5 |
|
if (!self::$proj4->hasDef('EPSG:' . $srid)) { |
227
|
|
|
|
228
|
2 |
|
$projDefs = Config::inst()->get(self::class, 'projections'); |
229
|
2 |
|
|
230
|
|
|
if (!isset($projDefs[$srid])) { |
231
|
|
|
throw new Exception("Cannot use unregistered SRID $srid. Register it's <a href=\"http://spatialreference.org/ref/epsg/$srid/proj4/\">PROJ.4 definition</a> in GIS::projections."); |
232
|
2 |
|
} |
233
|
|
|
|
234
|
|
|
self::$proj4->addDef('EPSG:' . $srid, $projDefs[$srid]); |
235
|
5 |
|
} |
236
|
|
|
|
237
|
|
|
return new Proj('EPSG:' . $srid, self::$proj4); |
238
|
1 |
|
} |
239
|
|
|
|
240
|
1 |
|
protected static function reproject_array($coordinates, $fromProj, $toProj) |
241
|
|
|
{ |
242
|
|
|
return self::each($coordinates, function ($coordinate) use ($fromProj, $toProj) { |
243
|
|
|
return array_slice(self::$proj4->transform($toProj, new Point($coordinate[0], $coordinate[1], $fromProj))->toArray(), 0, 2); |
244
|
|
|
}); |
245
|
|
|
} |
246
|
|
|
|
247
|
|
|
public static function each($coordinates, $callback) |
248
|
|
|
{ |
249
|
|
|
if ($coordinates instanceof GIS) { |
250
|
|
|
$coordinates = $coordinates->coordinates; |
251
|
|
|
} |
252
|
|
|
|
253
|
|
|
if (is_array($coordinates[0])) { |
254
|
|
|
|
255
|
|
|
foreach ($coordinates as &$coordinate) { |
256
|
|
|
$coordinate = self::each($coordinate, $callback); |
257
|
|
|
} |
258
|
|
|
|
259
|
|
|
return $coordinates; |
260
|
|
|
} |
261
|
|
|
|
262
|
|
|
return $callback($coordinates); |
263
|
|
|
} |
264
|
|
|
|
265
|
|
|
public function distance($geo) |
266
|
|
|
{ |
267
|
|
|
return DB::query('select ' . DB::get_schema()->translateDistanceQuery($this, GIS::create($geo)))->value(); |
|
|
|
|
268
|
|
|
} |
269
|
|
|
} |
270
|
|
|
|