pgLatLon
view latlon--0.3.sql @ 9:dc8da756b761
Changed version to 0.3
| author | jbe |
|---|---|
| date | Fri Sep 02 14:04:04 2016 +0200 (2016-09-02) |
| parents | latlon--0.2.sql@8f79d4a1cdc1 |
| children | 684a78d2f9f0 |
line source
2 ----------------------------------------
3 -- forward declarations (shell types) --
4 ----------------------------------------
6 CREATE TYPE epoint;
7 CREATE TYPE ebox;
8 CREATE TYPE ecircle;
9 CREATE TYPE ecluster;
12 ------------------------------------------------------------
13 -- dummy input/output functions for dummy index key types --
14 ------------------------------------------------------------
16 CREATE FUNCTION ekey_point_in_dummy(cstring)
17 RETURNS ekey_point
18 LANGUAGE C IMMUTABLE STRICT
19 AS '$libdir/latlon-v0003', 'pgl_notimpl';
21 CREATE FUNCTION ekey_point_out_dummy(ekey_point)
22 RETURNS cstring
23 LANGUAGE C IMMUTABLE STRICT
24 AS '$libdir/latlon-v0003', 'pgl_notimpl';
26 CREATE FUNCTION ekey_area_in_dummy(cstring)
27 RETURNS ekey_area
28 LANGUAGE C IMMUTABLE STRICT
29 AS '$libdir/latlon-v0003', 'pgl_notimpl';
31 CREATE FUNCTION ekey_area_out_dummy(ekey_area)
32 RETURNS cstring
33 LANGUAGE C IMMUTABLE STRICT
34 AS '$libdir/latlon-v0003', 'pgl_notimpl';
37 --------------------------
38 -- text input functions --
39 --------------------------
41 CREATE FUNCTION epoint_in(cstring)
42 RETURNS epoint
43 LANGUAGE C IMMUTABLE STRICT
44 AS '$libdir/latlon-v0003', 'pgl_epoint_in';
46 CREATE FUNCTION ebox_in(cstring)
47 RETURNS ebox
48 LANGUAGE C IMMUTABLE STRICT
49 AS '$libdir/latlon-v0003', 'pgl_ebox_in';
51 CREATE FUNCTION ecircle_in(cstring)
52 RETURNS ecircle
53 LANGUAGE C IMMUTABLE STRICT
54 AS '$libdir/latlon-v0003', 'pgl_ecircle_in';
56 CREATE FUNCTION ecluster_in(cstring)
57 RETURNS ecluster
58 LANGUAGE C IMMUTABLE STRICT
59 AS '$libdir/latlon-v0003', 'pgl_ecluster_in';
62 ---------------------------
63 -- text output functions --
64 ---------------------------
66 CREATE FUNCTION epoint_out(epoint)
67 RETURNS cstring
68 LANGUAGE C IMMUTABLE STRICT
69 AS '$libdir/latlon-v0003', 'pgl_epoint_out';
71 CREATE FUNCTION ebox_out(ebox)
72 RETURNS cstring
73 LANGUAGE C IMMUTABLE STRICT
74 AS '$libdir/latlon-v0003', 'pgl_ebox_out';
76 CREATE FUNCTION ecircle_out(ecircle)
77 RETURNS cstring
78 LANGUAGE C IMMUTABLE STRICT
79 AS '$libdir/latlon-v0003', 'pgl_ecircle_out';
81 CREATE FUNCTION ecluster_out(ecluster)
82 RETURNS cstring
83 LANGUAGE C IMMUTABLE STRICT
84 AS '$libdir/latlon-v0003', 'pgl_ecluster_out';
87 --------------------------
88 -- binary I/O functions --
89 --------------------------
91 CREATE FUNCTION epoint_recv(internal)
92 RETURNS epoint
93 LANGUAGE C IMMUTABLE STRICT
94 AS '$libdir/latlon-v0003', 'pgl_epoint_recv';
96 CREATE FUNCTION ebox_recv(internal)
97 RETURNS ebox
98 LANGUAGE C IMMUTABLE STRICT
99 AS '$libdir/latlon-v0003', 'pgl_ebox_recv';
101 CREATE FUNCTION ecircle_recv(internal)
102 RETURNS ecircle
103 LANGUAGE C IMMUTABLE STRICT
104 AS '$libdir/latlon-v0003', 'pgl_ecircle_recv';
106 CREATE FUNCTION epoint_send(epoint)
107 RETURNS bytea
108 LANGUAGE C IMMUTABLE STRICT
109 AS '$libdir/latlon-v0003', 'pgl_epoint_send';
111 CREATE FUNCTION ebox_send(ebox)
112 RETURNS bytea
113 LANGUAGE C IMMUTABLE STRICT
114 AS '$libdir/latlon-v0003', 'pgl_ebox_send';
116 CREATE FUNCTION ecircle_send(ecircle)
117 RETURNS bytea
118 LANGUAGE C IMMUTABLE STRICT
119 AS '$libdir/latlon-v0003', 'pgl_ecircle_send';
122 -----------------------------------------------
123 -- type definitions of dummy index key types --
124 -----------------------------------------------
126 CREATE TYPE ekey_point (
127 internallength = 8,
128 input = ekey_point_in_dummy,
129 output = ekey_point_out_dummy,
130 alignment = char );
132 CREATE TYPE ekey_area (
133 internallength = 9,
134 input = ekey_area_in_dummy,
135 output = ekey_area_out_dummy,
136 alignment = char );
139 ------------------------------------------
140 -- definitions of geographic data types --
141 ------------------------------------------
143 CREATE TYPE epoint (
144 internallength = 16,
145 input = epoint_in,
146 output = epoint_out,
147 receive = epoint_recv,
148 send = epoint_send,
149 alignment = double );
151 CREATE TYPE ebox (
152 internallength = 32,
153 input = ebox_in,
154 output = ebox_out,
155 receive = ebox_recv,
156 send = ebox_send,
157 alignment = double );
159 CREATE TYPE ecircle (
160 internallength = 24,
161 input = ecircle_in,
162 output = ecircle_out,
163 receive = ecircle_recv,
164 send = ecircle_send,
165 alignment = double );
167 CREATE TYPE ecluster (
168 internallength = VARIABLE,
169 input = ecluster_in,
170 output = ecluster_out,
171 alignment = double,
172 storage = external );
175 --------------------
176 -- B-tree support --
177 --------------------
179 -- begin of B-tree support for epoint
181 CREATE FUNCTION epoint_btree_lt(epoint, epoint)
182 RETURNS boolean
183 LANGUAGE C IMMUTABLE STRICT
184 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_lt';
186 CREATE FUNCTION epoint_btree_le(epoint, epoint)
187 RETURNS boolean
188 LANGUAGE C IMMUTABLE STRICT
189 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_le';
191 CREATE FUNCTION epoint_btree_eq(epoint, epoint)
192 RETURNS boolean
193 LANGUAGE C IMMUTABLE STRICT
194 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_eq';
196 CREATE FUNCTION epoint_btree_ne(epoint, epoint)
197 RETURNS boolean
198 LANGUAGE C IMMUTABLE STRICT
199 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_ne';
201 CREATE FUNCTION epoint_btree_ge(epoint, epoint)
202 RETURNS boolean
203 LANGUAGE C IMMUTABLE STRICT
204 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_ge';
206 CREATE FUNCTION epoint_btree_gt(epoint, epoint)
207 RETURNS boolean
208 LANGUAGE C IMMUTABLE STRICT
209 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_gt';
211 CREATE OPERATOR <<< (
212 leftarg = epoint,
213 rightarg = epoint,
214 procedure = epoint_btree_lt,
215 commutator = >>>,
216 negator = >>>=,
217 restrict = scalarltsel,
218 join = scalarltjoinsel
219 );
221 CREATE OPERATOR <<<= (
222 leftarg = epoint,
223 rightarg = epoint,
224 procedure = epoint_btree_le,
225 commutator = >>>=,
226 negator = >>>,
227 restrict = scalarltsel,
228 join = scalarltjoinsel
229 );
231 CREATE OPERATOR = (
232 leftarg = epoint,
233 rightarg = epoint,
234 procedure = epoint_btree_eq,
235 commutator = =,
236 negator = <>,
237 restrict = eqsel,
238 join = eqjoinsel,
239 merges
240 );
242 CREATE OPERATOR <> (
243 leftarg = epoint,
244 rightarg = epoint,
245 procedure = epoint_btree_eq,
246 commutator = <>,
247 negator = =,
248 restrict = neqsel,
249 join = neqjoinsel
250 );
252 CREATE OPERATOR >>>= (
253 leftarg = epoint,
254 rightarg = epoint,
255 procedure = epoint_btree_ge,
256 commutator = <<<=,
257 negator = <<<,
258 restrict = scalargtsel,
259 join = scalargtjoinsel
260 );
262 CREATE OPERATOR >>> (
263 leftarg = epoint,
264 rightarg = epoint,
265 procedure = epoint_btree_gt,
266 commutator = <<<,
267 negator = <<<=,
268 restrict = scalargtsel,
269 join = scalargtjoinsel
270 );
272 CREATE FUNCTION epoint_btree_cmp(epoint, epoint)
273 RETURNS int4
274 LANGUAGE C IMMUTABLE STRICT
275 AS '$libdir/latlon-v0003', 'pgl_btree_epoint_cmp';
277 CREATE OPERATOR CLASS epoint_btree_ops
278 DEFAULT FOR TYPE epoint USING btree AS
279 OPERATOR 1 <<< ,
280 OPERATOR 2 <<<= ,
281 OPERATOR 3 = ,
282 OPERATOR 4 >>>= ,
283 OPERATOR 5 >>> ,
284 FUNCTION 1 epoint_btree_cmp(epoint, epoint);
286 -- end of B-tree support for epoint
288 -- begin of B-tree support for ebox
290 CREATE FUNCTION ebox_btree_lt(ebox, ebox)
291 RETURNS boolean
292 LANGUAGE C IMMUTABLE STRICT
293 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_lt';
295 CREATE FUNCTION ebox_btree_le(ebox, ebox)
296 RETURNS boolean
297 LANGUAGE C IMMUTABLE STRICT
298 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_le';
300 CREATE FUNCTION ebox_btree_eq(ebox, ebox)
301 RETURNS boolean
302 LANGUAGE C IMMUTABLE STRICT
303 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_eq';
305 CREATE FUNCTION ebox_btree_ne(ebox, ebox)
306 RETURNS boolean
307 LANGUAGE C IMMUTABLE STRICT
308 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_ne';
310 CREATE FUNCTION ebox_btree_ge(ebox, ebox)
311 RETURNS boolean
312 LANGUAGE C IMMUTABLE STRICT
313 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_ge';
315 CREATE FUNCTION ebox_btree_gt(ebox, ebox)
316 RETURNS boolean
317 LANGUAGE C IMMUTABLE STRICT
318 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_gt';
320 CREATE OPERATOR <<< (
321 leftarg = ebox,
322 rightarg = ebox,
323 procedure = ebox_btree_lt,
324 commutator = >>>,
325 negator = >>>=,
326 restrict = scalarltsel,
327 join = scalarltjoinsel
328 );
330 CREATE OPERATOR <<<= (
331 leftarg = ebox,
332 rightarg = ebox,
333 procedure = ebox_btree_le,
334 commutator = >>>=,
335 negator = >>>,
336 restrict = scalarltsel,
337 join = scalarltjoinsel
338 );
340 CREATE OPERATOR = (
341 leftarg = ebox,
342 rightarg = ebox,
343 procedure = ebox_btree_eq,
344 commutator = =,
345 negator = <>,
346 restrict = eqsel,
347 join = eqjoinsel,
348 merges
349 );
351 CREATE OPERATOR <> (
352 leftarg = ebox,
353 rightarg = ebox,
354 procedure = ebox_btree_eq,
355 commutator = <>,
356 negator = =,
357 restrict = neqsel,
358 join = neqjoinsel
359 );
361 CREATE OPERATOR >>>= (
362 leftarg = ebox,
363 rightarg = ebox,
364 procedure = ebox_btree_ge,
365 commutator = <<<=,
366 negator = <<<,
367 restrict = scalargtsel,
368 join = scalargtjoinsel
369 );
371 CREATE OPERATOR >>> (
372 leftarg = ebox,
373 rightarg = ebox,
374 procedure = ebox_btree_gt,
375 commutator = <<<,
376 negator = <<<=,
377 restrict = scalargtsel,
378 join = scalargtjoinsel
379 );
381 CREATE FUNCTION ebox_btree_cmp(ebox, ebox)
382 RETURNS int4
383 LANGUAGE C IMMUTABLE STRICT
384 AS '$libdir/latlon-v0003', 'pgl_btree_ebox_cmp';
386 CREATE OPERATOR CLASS ebox_btree_ops
387 DEFAULT FOR TYPE ebox USING btree AS
388 OPERATOR 1 <<< ,
389 OPERATOR 2 <<<= ,
390 OPERATOR 3 = ,
391 OPERATOR 4 >>>= ,
392 OPERATOR 5 >>> ,
393 FUNCTION 1 ebox_btree_cmp(ebox, ebox);
395 -- end of B-tree support for ebox
397 -- begin of B-tree support for ecircle
399 CREATE FUNCTION ecircle_btree_lt(ecircle, ecircle)
400 RETURNS boolean
401 LANGUAGE C IMMUTABLE STRICT
402 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_lt';
404 CREATE FUNCTION ecircle_btree_le(ecircle, ecircle)
405 RETURNS boolean
406 LANGUAGE C IMMUTABLE STRICT
407 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_le';
409 CREATE FUNCTION ecircle_btree_eq(ecircle, ecircle)
410 RETURNS boolean
411 LANGUAGE C IMMUTABLE STRICT
412 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_eq';
414 CREATE FUNCTION ecircle_btree_ne(ecircle, ecircle)
415 RETURNS boolean
416 LANGUAGE C IMMUTABLE STRICT
417 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_ne';
419 CREATE FUNCTION ecircle_btree_ge(ecircle, ecircle)
420 RETURNS boolean
421 LANGUAGE C IMMUTABLE STRICT
422 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_ge';
424 CREATE FUNCTION ecircle_btree_gt(ecircle, ecircle)
425 RETURNS boolean
426 LANGUAGE C IMMUTABLE STRICT
427 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_gt';
429 CREATE OPERATOR <<< (
430 leftarg = ecircle,
431 rightarg = ecircle,
432 procedure = ecircle_btree_lt,
433 commutator = >>>,
434 negator = >>>=,
435 restrict = scalarltsel,
436 join = scalarltjoinsel
437 );
439 CREATE OPERATOR <<<= (
440 leftarg = ecircle,
441 rightarg = ecircle,
442 procedure = ecircle_btree_le,
443 commutator = >>>=,
444 negator = >>>,
445 restrict = scalarltsel,
446 join = scalarltjoinsel
447 );
449 CREATE OPERATOR = (
450 leftarg = ecircle,
451 rightarg = ecircle,
452 procedure = ecircle_btree_eq,
453 commutator = =,
454 negator = <>,
455 restrict = eqsel,
456 join = eqjoinsel,
457 merges
458 );
460 CREATE OPERATOR <> (
461 leftarg = ecircle,
462 rightarg = ecircle,
463 procedure = ecircle_btree_eq,
464 commutator = <>,
465 negator = =,
466 restrict = neqsel,
467 join = neqjoinsel
468 );
470 CREATE OPERATOR >>>= (
471 leftarg = ecircle,
472 rightarg = ecircle,
473 procedure = ecircle_btree_ge,
474 commutator = <<<=,
475 negator = <<<,
476 restrict = scalargtsel,
477 join = scalargtjoinsel
478 );
480 CREATE OPERATOR >>> (
481 leftarg = ecircle,
482 rightarg = ecircle,
483 procedure = ecircle_btree_gt,
484 commutator = <<<,
485 negator = <<<=,
486 restrict = scalargtsel,
487 join = scalargtjoinsel
488 );
490 CREATE FUNCTION ecircle_btree_cmp(ecircle, ecircle)
491 RETURNS int4
492 LANGUAGE C IMMUTABLE STRICT
493 AS '$libdir/latlon-v0003', 'pgl_btree_ecircle_cmp';
495 CREATE OPERATOR CLASS ecircle_btree_ops
496 DEFAULT FOR TYPE ecircle USING btree AS
497 OPERATOR 1 <<< ,
498 OPERATOR 2 <<<= ,
499 OPERATOR 3 = ,
500 OPERATOR 4 >>>= ,
501 OPERATOR 5 >>> ,
502 FUNCTION 1 ecircle_btree_cmp(ecircle, ecircle);
504 -- end of B-tree support for ecircle
507 ----------------
508 -- type casts --
509 ----------------
511 CREATE FUNCTION cast_epoint_to_ebox(epoint)
512 RETURNS ebox
513 LANGUAGE C IMMUTABLE STRICT
514 AS '$libdir/latlon-v0003', 'pgl_epoint_to_ebox';
516 CREATE CAST (epoint AS ebox) WITH FUNCTION cast_epoint_to_ebox(epoint);
518 CREATE FUNCTION cast_epoint_to_ecircle(epoint)
519 RETURNS ecircle
520 LANGUAGE C IMMUTABLE STRICT
521 AS '$libdir/latlon-v0003', 'pgl_epoint_to_ecircle';
523 CREATE CAST (epoint AS ecircle) WITH FUNCTION cast_epoint_to_ecircle(epoint);
525 CREATE FUNCTION cast_epoint_to_ecluster(epoint)
526 RETURNS ecluster
527 LANGUAGE C IMMUTABLE STRICT
528 AS '$libdir/latlon-v0003', 'pgl_epoint_to_ecluster';
530 CREATE CAST (epoint AS ecluster) WITH FUNCTION cast_epoint_to_ecluster(epoint);
532 CREATE FUNCTION cast_ebox_to_ecluster(ebox)
533 RETURNS ecluster
534 LANGUAGE C IMMUTABLE STRICT
535 AS '$libdir/latlon-v0003', 'pgl_ebox_to_ecluster';
537 CREATE CAST (ebox AS ecluster) WITH FUNCTION cast_ebox_to_ecluster(ebox);
540 ---------------------------
541 -- constructor functions --
542 ---------------------------
544 CREATE FUNCTION epoint(float8, float8)
545 RETURNS epoint
546 LANGUAGE C IMMUTABLE STRICT
547 AS '$libdir/latlon-v0003', 'pgl_create_epoint';
549 CREATE FUNCTION epoint_latlon(float8, float8)
550 RETURNS epoint
551 LANGUAGE SQL IMMUTABLE STRICT AS $$
552 SELECT epoint($1, $2)
553 $$;
555 CREATE FUNCTION epoint_lonlat(float8, float8)
556 RETURNS epoint
557 LANGUAGE SQL IMMUTABLE STRICT AS $$
558 SELECT epoint($2, $1)
559 $$;
561 CREATE FUNCTION empty_ebox()
562 RETURNS ebox
563 LANGUAGE C IMMUTABLE STRICT
564 AS '$libdir/latlon-v0003', 'pgl_create_empty_ebox';
566 CREATE FUNCTION ebox(float8, float8, float8, float8)
567 RETURNS ebox
568 LANGUAGE C IMMUTABLE STRICT
569 AS '$libdir/latlon-v0003', 'pgl_create_ebox';
571 CREATE FUNCTION ebox(epoint, epoint)
572 RETURNS ebox
573 LANGUAGE C IMMUTABLE STRICT
574 AS '$libdir/latlon-v0003', 'pgl_create_ebox_from_epoints';
576 CREATE FUNCTION ecircle(float8, float8, float8)
577 RETURNS ecircle
578 LANGUAGE C IMMUTABLE STRICT
579 AS '$libdir/latlon-v0003', 'pgl_create_ecircle';
581 CREATE FUNCTION ecircle(epoint, float8)
582 RETURNS ecircle
583 LANGUAGE C IMMUTABLE STRICT
584 AS '$libdir/latlon-v0003', 'pgl_create_ecircle_from_epoint';
586 CREATE FUNCTION ecluster_concat(ecluster[])
587 RETURNS ecluster
588 LANGUAGE sql IMMUTABLE STRICT AS $$
589 SELECT array_to_string($1, ' ')::ecluster
590 $$;
592 CREATE FUNCTION ecluster_concat(ecluster, ecluster)
593 RETURNS ecluster
594 LANGUAGE sql IMMUTABLE STRICT AS $$
595 SELECT ($1::text || ' ' || $2::text)::ecluster
596 $$;
598 CREATE FUNCTION ecluster_create_multipoint(epoint[])
599 RETURNS ecluster
600 LANGUAGE sql IMMUTABLE STRICT AS $$
601 SELECT
602 array_to_string(array_agg('point (' || unnest || ')'), ' ')::ecluster
603 FROM unnest($1)
604 $$;
606 CREATE FUNCTION ecluster_create_path(epoint[])
607 RETURNS ecluster
608 LANGUAGE sql IMMUTABLE STRICT AS $$
609 SELECT CASE WHEN "str" = '' THEN 'empty'::ecluster ELSE
610 ('path (' || array_to_string($1, ' ') || ')')::ecluster
611 END
612 FROM array_to_string($1, ' ') AS "str"
613 $$;
615 CREATE FUNCTION ecluster_create_outline(epoint[])
616 RETURNS ecluster
617 LANGUAGE sql IMMUTABLE STRICT AS $$
618 SELECT CASE WHEN "str" = '' THEN 'empty'::ecluster ELSE
619 ('outline (' || array_to_string($1, ' ') || ')')::ecluster
620 END
621 FROM array_to_string($1, ' ') AS "str"
622 $$;
624 CREATE FUNCTION ecluster_create_polygon(epoint[])
625 RETURNS ecluster
626 LANGUAGE sql IMMUTABLE STRICT AS $$
627 SELECT CASE WHEN "str" = '' THEN 'empty'::ecluster ELSE
628 ('polygon (' || array_to_string($1, ' ') || ')')::ecluster
629 END
630 FROM array_to_string($1, ' ') AS "str"
631 $$;
634 ----------------------
635 -- getter functions --
636 ----------------------
638 CREATE FUNCTION latitude(epoint)
639 RETURNS float8
640 LANGUAGE C IMMUTABLE STRICT
641 AS '$libdir/latlon-v0003', 'pgl_epoint_lat';
643 CREATE FUNCTION longitude(epoint)
644 RETURNS float8
645 LANGUAGE C IMMUTABLE STRICT
646 AS '$libdir/latlon-v0003', 'pgl_epoint_lon';
648 CREATE FUNCTION min_latitude(ebox)
649 RETURNS float8
650 LANGUAGE C IMMUTABLE STRICT
651 AS '$libdir/latlon-v0003', 'pgl_ebox_lat_min';
653 CREATE FUNCTION max_latitude(ebox)
654 RETURNS float8
655 LANGUAGE C IMMUTABLE STRICT
656 AS '$libdir/latlon-v0003', 'pgl_ebox_lat_max';
658 CREATE FUNCTION min_longitude(ebox)
659 RETURNS float8
660 LANGUAGE C IMMUTABLE STRICT
661 AS '$libdir/latlon-v0003', 'pgl_ebox_lon_min';
663 CREATE FUNCTION max_longitude(ebox)
664 RETURNS float8
665 LANGUAGE C IMMUTABLE STRICT
666 AS '$libdir/latlon-v0003', 'pgl_ebox_lon_max';
668 CREATE FUNCTION center(ecircle)
669 RETURNS epoint
670 LANGUAGE C IMMUTABLE STRICT
671 AS '$libdir/latlon-v0003', 'pgl_ecircle_center';
673 CREATE FUNCTION radius(ecircle)
674 RETURNS float8
675 LANGUAGE C IMMUTABLE STRICT
676 AS '$libdir/latlon-v0003', 'pgl_ecircle_radius';
678 CREATE FUNCTION ecluster_extract_points(ecluster)
679 RETURNS SETOF epoint
680 LANGUAGE sql IMMUTABLE STRICT AS $$
681 SELECT "match"[2]::epoint
682 FROM regexp_matches($1::text, e'(^| )point \\(([^)]+)\\)', 'g') AS "match"
683 $$;
685 CREATE FUNCTION ecluster_extract_paths(ecluster)
686 RETURNS SETOF epoint[]
687 LANGUAGE sql IMMUTABLE STRICT AS $$
688 SELECT (
689 SELECT array_agg("m2"[1]::epoint)
690 FROM regexp_matches("m1"[2], e'[^ ]+ [^ ]+', 'g') AS "m2"
691 )
692 FROM regexp_matches($1::text, e'(^| )path \\(([^)]+)\\)', 'g') AS "m1"
693 $$;
695 CREATE FUNCTION ecluster_extract_outlines(ecluster)
696 RETURNS SETOF epoint[]
697 LANGUAGE sql IMMUTABLE STRICT AS $$
698 SELECT (
699 SELECT array_agg("m2"[1]::epoint)
700 FROM regexp_matches("m1"[2], e'[^ ]+ [^ ]+', 'g') AS "m2"
701 )
702 FROM regexp_matches($1::text, e'(^| )outline \\(([^)]+)\\)', 'g') AS "m1"
703 $$;
705 CREATE FUNCTION ecluster_extract_polygons(ecluster)
706 RETURNS SETOF epoint[]
707 LANGUAGE sql IMMUTABLE STRICT AS $$
708 SELECT (
709 SELECT array_agg("m2"[1]::epoint)
710 FROM regexp_matches("m1"[2], e'[^ ]+ [^ ]+', 'g') AS "m2"
711 )
712 FROM regexp_matches($1::text, e'(^| )polygon \\(([^)]+)\\)', 'g') AS "m1"
713 $$;
716 ---------------
717 -- operators --
718 ---------------
720 CREATE FUNCTION epoint_ebox_overlap_proc(epoint, ebox)
721 RETURNS boolean
722 LANGUAGE C IMMUTABLE STRICT
723 AS '$libdir/latlon-v0003', 'pgl_epoint_ebox_overlap';
725 CREATE FUNCTION epoint_ecircle_overlap_proc(epoint, ecircle)
726 RETURNS boolean
727 LANGUAGE C IMMUTABLE STRICT
728 AS '$libdir/latlon-v0003', 'pgl_epoint_ecircle_overlap';
730 CREATE FUNCTION epoint_ecluster_overlap_proc(epoint, ecluster)
731 RETURNS boolean
732 LANGUAGE C IMMUTABLE STRICT
733 AS '$libdir/latlon-v0003', 'pgl_epoint_ecluster_overlap';
735 CREATE FUNCTION ebox_overlap_proc(ebox, ebox)
736 RETURNS boolean
737 LANGUAGE C IMMUTABLE STRICT
738 AS '$libdir/latlon-v0003', 'pgl_ebox_overlap';
740 CREATE FUNCTION ecircle_overlap_proc(ecircle, ecircle)
741 RETURNS boolean
742 LANGUAGE C IMMUTABLE STRICT
743 AS '$libdir/latlon-v0003', 'pgl_ecircle_overlap';
745 CREATE FUNCTION ecircle_ecluster_overlap_proc(ecircle, ecluster)
746 RETURNS boolean
747 LANGUAGE C IMMUTABLE STRICT
748 AS '$libdir/latlon-v0003', 'pgl_ecircle_ecluster_overlap';
750 CREATE FUNCTION epoint_distance_proc(epoint, epoint)
751 RETURNS float8
752 LANGUAGE C IMMUTABLE STRICT
753 AS '$libdir/latlon-v0003', 'pgl_epoint_distance';
755 CREATE FUNCTION epoint_ecircle_distance_proc(epoint, ecircle)
756 RETURNS float8
757 LANGUAGE C IMMUTABLE STRICT
758 AS '$libdir/latlon-v0003', 'pgl_epoint_ecircle_distance';
760 CREATE FUNCTION epoint_ecluster_distance_proc(epoint, ecluster)
761 RETURNS float8
762 LANGUAGE C IMMUTABLE STRICT
763 AS '$libdir/latlon-v0003', 'pgl_epoint_ecluster_distance';
765 CREATE FUNCTION ecircle_distance_proc(ecircle, ecircle)
766 RETURNS float8
767 LANGUAGE C IMMUTABLE STRICT
768 AS '$libdir/latlon-v0003', 'pgl_ecircle_distance';
770 CREATE FUNCTION ecircle_ecluster_distance_proc(ecircle, ecluster)
771 RETURNS float8
772 LANGUAGE C IMMUTABLE STRICT
773 AS '$libdir/latlon-v0003', 'pgl_ecircle_ecluster_distance';
775 CREATE OPERATOR && (
776 leftarg = epoint,
777 rightarg = ebox,
778 procedure = epoint_ebox_overlap_proc,
779 commutator = &&,
780 restrict = areasel,
781 join = areajoinsel
782 );
784 CREATE FUNCTION epoint_ebox_overlap_commutator(ebox, epoint)
785 RETURNS boolean
786 LANGUAGE sql IMMUTABLE AS 'SELECT $2 && $1';
788 CREATE OPERATOR && (
789 leftarg = ebox,
790 rightarg = epoint,
791 procedure = epoint_ebox_overlap_commutator,
792 commutator = &&,
793 restrict = areasel,
794 join = areajoinsel
795 );
797 CREATE OPERATOR && (
798 leftarg = epoint,
799 rightarg = ecircle,
800 procedure = epoint_ecircle_overlap_proc,
801 commutator = &&,
802 restrict = areasel,
803 join = areajoinsel
804 );
806 CREATE FUNCTION epoint_ecircle_overlap_commutator(ecircle, epoint)
807 RETURNS boolean
808 LANGUAGE sql IMMUTABLE AS 'SELECT $2 && $1';
810 CREATE OPERATOR && (
811 leftarg = ecircle,
812 rightarg = epoint,
813 procedure = epoint_ecircle_overlap_commutator,
814 commutator = &&,
815 restrict = areasel,
816 join = areajoinsel
817 );
819 CREATE OPERATOR && (
820 leftarg = epoint,
821 rightarg = ecluster,
822 procedure = epoint_ecluster_overlap_proc,
823 commutator = &&,
824 restrict = areasel,
825 join = areajoinsel
826 );
828 CREATE FUNCTION epoint_ecluster_overlap_commutator(ecluster, epoint)
829 RETURNS boolean
830 LANGUAGE sql IMMUTABLE AS 'SELECT $2 && $1';
832 CREATE OPERATOR && (
833 leftarg = ecluster,
834 rightarg = epoint,
835 procedure = epoint_ecluster_overlap_commutator,
836 commutator = &&,
837 restrict = areasel,
838 join = areajoinsel
839 );
841 CREATE OPERATOR && (
842 leftarg = ebox,
843 rightarg = ebox,
844 procedure = ebox_overlap_proc,
845 commutator = &&,
846 restrict = areasel,
847 join = areajoinsel
848 );
850 CREATE OPERATOR && (
851 leftarg = ecircle,
852 rightarg = ecircle,
853 procedure = ecircle_overlap_proc,
854 commutator = &&,
855 restrict = areasel,
856 join = areajoinsel
857 );
859 CREATE OPERATOR && (
860 leftarg = ecircle,
861 rightarg = ecluster,
862 procedure = ecircle_ecluster_overlap_proc,
863 commutator = &&,
864 restrict = areasel,
865 join = areajoinsel
866 );
868 CREATE FUNCTION ecircle_ecluster_overlap_commutator(ecluster, ecircle)
869 RETURNS boolean
870 LANGUAGE sql IMMUTABLE AS 'SELECT $2 && $1';
872 CREATE OPERATOR && (
873 leftarg = ecluster,
874 rightarg = ecircle,
875 procedure = ecircle_ecluster_overlap_commutator,
876 commutator = &&,
877 restrict = areasel,
878 join = areajoinsel
879 );
881 CREATE OPERATOR <-> (
882 leftarg = epoint,
883 rightarg = epoint,
884 procedure = epoint_distance_proc,
885 commutator = <->
886 );
888 CREATE OPERATOR <-> (
889 leftarg = epoint,
890 rightarg = ecircle,
891 procedure = epoint_ecircle_distance_proc,
892 commutator = <->
893 );
895 CREATE FUNCTION epoint_ecircle_distance_commutator(ecircle, epoint)
896 RETURNS float8
897 LANGUAGE sql IMMUTABLE AS 'SELECT $2 <-> $1';
899 CREATE OPERATOR <-> (
900 leftarg = ecircle,
901 rightarg = epoint,
902 procedure = epoint_ecircle_distance_commutator,
903 commutator = <->
904 );
906 CREATE OPERATOR <-> (
907 leftarg = epoint,
908 rightarg = ecluster,
909 procedure = epoint_ecluster_distance_proc,
910 commutator = <->
911 );
913 CREATE FUNCTION epoint_ecluster_distance_commutator(ecluster, epoint)
914 RETURNS float8
915 LANGUAGE sql IMMUTABLE AS 'SELECT $2 <-> $1';
917 CREATE OPERATOR <-> (
918 leftarg = ecluster,
919 rightarg = epoint,
920 procedure = epoint_ecluster_distance_commutator,
921 commutator = <->
922 );
924 CREATE OPERATOR <-> (
925 leftarg = ecircle,
926 rightarg = ecircle,
927 procedure = ecircle_distance_proc,
928 commutator = <->
929 );
931 CREATE OPERATOR <-> (
932 leftarg = ecircle,
933 rightarg = ecluster,
934 procedure = ecircle_ecluster_distance_proc,
935 commutator = <->
936 );
938 CREATE FUNCTION ecircle_ecluster_distance_commutator(ecluster, ecircle)
939 RETURNS float8
940 LANGUAGE sql IMMUTABLE AS 'SELECT $2 <-> $1';
942 CREATE OPERATOR <-> (
943 leftarg = ecluster,
944 rightarg = ecircle,
945 procedure = ecircle_ecluster_distance_commutator,
946 commutator = <->
947 );
950 ----------------
951 -- GiST index --
952 ----------------
954 CREATE FUNCTION pgl_gist_consistent(internal, internal, smallint, oid, internal)
955 RETURNS boolean
956 LANGUAGE C STRICT
957 AS '$libdir/latlon-v0003', 'pgl_gist_consistent';
959 CREATE FUNCTION pgl_gist_union(internal, internal)
960 RETURNS internal
961 LANGUAGE C STRICT
962 AS '$libdir/latlon-v0003', 'pgl_gist_union';
964 CREATE FUNCTION pgl_gist_compress_epoint(internal)
965 RETURNS internal
966 LANGUAGE C STRICT
967 AS '$libdir/latlon-v0003', 'pgl_gist_compress_epoint';
969 CREATE FUNCTION pgl_gist_compress_ecircle(internal)
970 RETURNS internal
971 LANGUAGE C STRICT
972 AS '$libdir/latlon-v0003', 'pgl_gist_compress_ecircle';
974 CREATE FUNCTION pgl_gist_compress_ecluster(internal)
975 RETURNS internal
976 LANGUAGE C STRICT
977 AS '$libdir/latlon-v0003', 'pgl_gist_compress_ecluster';
979 CREATE FUNCTION pgl_gist_decompress(internal)
980 RETURNS internal
981 LANGUAGE C STRICT
982 AS '$libdir/latlon-v0003', 'pgl_gist_decompress';
984 CREATE FUNCTION pgl_gist_penalty(internal, internal, internal)
985 RETURNS internal
986 LANGUAGE C STRICT
987 AS '$libdir/latlon-v0003', 'pgl_gist_penalty';
989 CREATE FUNCTION pgl_gist_picksplit(internal, internal)
990 RETURNS internal
991 LANGUAGE C STRICT
992 AS '$libdir/latlon-v0003', 'pgl_gist_picksplit';
994 CREATE FUNCTION pgl_gist_same(internal, internal, internal)
995 RETURNS internal
996 LANGUAGE C STRICT
997 AS '$libdir/latlon-v0003', 'pgl_gist_same';
999 CREATE FUNCTION pgl_gist_distance(internal, internal, smallint, oid)
1000 RETURNS internal
1001 LANGUAGE C STRICT
1002 AS '$libdir/latlon-v0003', 'pgl_gist_distance';
1004 CREATE OPERATOR CLASS epoint_ops
1005 DEFAULT FOR TYPE epoint USING gist AS
1006 OPERATOR 11 = ,
1007 OPERATOR 22 && (epoint, ebox),
1008 OPERATOR 23 && (epoint, ecircle),
1009 OPERATOR 24 && (epoint, ecluster),
1010 OPERATOR 31 <-> (epoint, epoint) FOR ORDER BY float_ops,
1011 OPERATOR 33 <-> (epoint, ecircle) FOR ORDER BY float_ops,
1012 OPERATOR 34 <-> (epoint, ecluster) FOR ORDER BY float_ops,
1013 FUNCTION 1 pgl_gist_consistent(internal, internal, smallint, oid, internal),
1014 FUNCTION 2 pgl_gist_union(internal, internal),
1015 FUNCTION 3 pgl_gist_compress_epoint(internal),
1016 FUNCTION 4 pgl_gist_decompress(internal),
1017 FUNCTION 5 pgl_gist_penalty(internal, internal, internal),
1018 FUNCTION 6 pgl_gist_picksplit(internal, internal),
1019 FUNCTION 7 pgl_gist_same(internal, internal, internal),
1020 FUNCTION 8 pgl_gist_distance(internal, internal, smallint, oid),
1021 STORAGE ekey_point;
1023 CREATE OPERATOR CLASS ecircle_ops
1024 DEFAULT FOR TYPE ecircle USING gist AS
1025 OPERATOR 13 = ,
1026 OPERATOR 21 && (ecircle, epoint),
1027 OPERATOR 23 && (ecircle, ecircle),
1028 OPERATOR 24 && (ecircle, ecluster),
1029 OPERATOR 31 <-> (ecircle, epoint) FOR ORDER BY float_ops,
1030 OPERATOR 33 <-> (ecircle, ecircle) FOR ORDER BY float_ops,
1031 OPERATOR 34 <-> (ecircle, ecluster) FOR ORDER BY float_ops,
1032 FUNCTION 1 pgl_gist_consistent(internal, internal, smallint, oid, internal),
1033 FUNCTION 2 pgl_gist_union(internal, internal),
1034 FUNCTION 3 pgl_gist_compress_ecircle(internal),
1035 FUNCTION 4 pgl_gist_decompress(internal),
1036 FUNCTION 5 pgl_gist_penalty(internal, internal, internal),
1037 FUNCTION 6 pgl_gist_picksplit(internal, internal),
1038 FUNCTION 7 pgl_gist_same(internal, internal, internal),
1039 FUNCTION 8 pgl_gist_distance(internal, internal, smallint, oid),
1040 STORAGE ekey_area;
1042 CREATE OPERATOR CLASS ecluster_ops
1043 DEFAULT FOR TYPE ecluster USING gist AS
1044 OPERATOR 21 && (ecluster, epoint),
1045 FUNCTION 1 pgl_gist_consistent(internal, internal, smallint, oid, internal),
1046 FUNCTION 2 pgl_gist_union(internal, internal),
1047 FUNCTION 3 pgl_gist_compress_ecluster(internal),
1048 FUNCTION 4 pgl_gist_decompress(internal),
1049 FUNCTION 5 pgl_gist_penalty(internal, internal, internal),
1050 FUNCTION 6 pgl_gist_picksplit(internal, internal),
1051 FUNCTION 7 pgl_gist_same(internal, internal, internal),
1052 FUNCTION 8 pgl_gist_distance(internal, internal, smallint, oid),
1053 STORAGE ekey_area;
1056 ---------------------
1057 -- alias functions --
1058 ---------------------
1060 CREATE FUNCTION distance(epoint, epoint)
1061 RETURNS float8
1062 LANGUAGE sql IMMUTABLE AS 'SELECT $1 <-> $2';
1064 CREATE FUNCTION distance(ecluster, epoint)
1065 RETURNS float8
1066 LANGUAGE sql IMMUTABLE AS 'SELECT $1 <-> $2';
1068 CREATE FUNCTION distance_within(epoint, epoint, float8)
1069 RETURNS boolean
1070 LANGUAGE sql IMMUTABLE AS 'SELECT $1 && ecircle($2, $3)';
1072 CREATE FUNCTION distance_within(ecluster, epoint, float8)
1073 RETURNS boolean
1074 LANGUAGE sql IMMUTABLE AS 'SELECT $1 && ecircle($2, $3)';
1077 --------------------------------
1078 -- other data storage formats --
1079 --------------------------------
1081 CREATE FUNCTION coords_to_epoint(float8, float8, text = 'epoint_lonlat')
1082 RETURNS epoint
1083 LANGUAGE plpgsql IMMUTABLE STRICT AS $$
1084 DECLARE
1085 "result" epoint;
1086 BEGIN
1087 IF $3 = 'epoint_lonlat' THEN
1088 -- avoid dynamic command execution for better performance
1089 RETURN epoint($2, $1);
1090 END IF;
1091 IF $3 = 'epoint' OR $3 = 'epoint_latlon' THEN
1092 -- avoid dynamic command execution for better performance
1093 RETURN epoint($1, $2);
1094 END IF;
1095 EXECUTE 'SELECT ' || $3 || '($1, $2)' INTO STRICT "result" USING $1, $2;
1096 RETURN "result";
1097 END;
1098 $$;
1100 CREATE FUNCTION GeoJSON_to_epoint(jsonb, text = 'epoint_lonlat')
1101 RETURNS epoint
1102 LANGUAGE sql IMMUTABLE STRICT AS $$
1103 SELECT CASE
1104 WHEN $1->>'type' = 'Point' THEN
1105 coords_to_epoint(
1106 ($1->'coordinates'->>1)::float8,
1107 ($1->'coordinates'->>0)::float8,
1108 $2
1109 )
1110 WHEN $1->>'type' = 'Feature' THEN
1111 GeoJSON_to_epoint($1->'geometry', $2)
1112 ELSE
1113 NULL
1114 END
1115 $$;
1117 CREATE FUNCTION GeoJSON_to_ecluster(jsonb, text = 'epoint_lonlat')
1118 RETURNS ecluster
1119 LANGUAGE sql IMMUTABLE STRICT AS $$
1120 SELECT CASE $1->>'type'
1121 WHEN 'Point' THEN
1122 coords_to_epoint(
1123 ($1->'coordinates'->>1)::float8,
1124 ($1->'coordinates'->>0)::float8,
1125 $2
1126 )::ecluster
1127 WHEN 'MultiPoint' THEN
1128 ( SELECT ecluster_create_multipoint(array_agg(
1129 coords_to_epoint(
1130 ("coord"->>1)::float8,
1131 ("coord"->>0)::float8,
1132 $2
1133 )
1134 ))
1135 FROM jsonb_array_elements($1->'coordinates') AS "coord"
1136 )
1137 WHEN 'LineString' THEN
1138 ( SELECT ecluster_create_path(array_agg(
1139 coords_to_epoint(
1140 ("coord"->>1)::float8,
1141 ("coord"->>0)::float8,
1142 $2
1143 )
1144 ))
1145 FROM jsonb_array_elements($1->'coordinates') AS "coord"
1146 )
1147 WHEN 'MultiLineString' THEN
1148 ( SELECT ecluster_concat(array_agg(
1149 ( SELECT ecluster_create_path(array_agg(
1150 coords_to_epoint(
1151 ("coord"->>1)::float8,
1152 ("coord"->>0)::float8,
1153 $2
1154 )
1155 ))
1156 FROM jsonb_array_elements("coord_array") AS "coord"
1157 )
1158 ))
1159 FROM jsonb_array_elements($1->'coordinates') AS "coord_array"
1160 )
1161 WHEN 'Polygon' THEN
1162 ( SELECT ecluster_concat(array_agg(
1163 ( SELECT ecluster_create_polygon(array_agg(
1164 coords_to_epoint(
1165 ("coord"->>1)::float8,
1166 ("coord"->>0)::float8,
1167 $2
1168 )
1169 ))
1170 FROM jsonb_array_elements("coord_array") AS "coord"
1171 )
1172 ))
1173 FROM jsonb_array_elements($1->'coordinates') AS "coord_array"
1174 )
1175 WHEN 'MultiPolygon' THEN
1176 ( SELECT ecluster_concat(array_agg(
1177 ( SELECT ecluster_concat(array_agg(
1178 ( SELECT ecluster_create_polygon(array_agg(
1179 coords_to_epoint(
1180 ("coord"->>1)::float8,
1181 ("coord"->>0)::float8,
1182 $2
1183 )
1184 ))
1185 FROM jsonb_array_elements("coord_array") AS "coord"
1186 )
1187 ))
1188 FROM jsonb_array_elements("coord_array_array") AS "coord_array"
1189 )
1190 ))
1191 FROM jsonb_array_elements($1->'coordinates') AS "coord_array_array"
1192 )
1193 WHEN 'Feature' THEN
1194 GeoJSON_to_ecluster($1->'geometry', $2)
1195 WHEN 'FeatureCollection' THEN
1196 ( SELECT ecluster_concat(array_agg(
1197 GeoJSON_to_ecluster("feature", $2)
1198 ))
1199 FROM jsonb_array_elements($1->'features') AS "feature"
1200 )
1201 ELSE
1202 NULL
1203 END
1204 $$;