однобитныйЗдрасте!
Создал таблицу для хранения геометрической фигуры. Помещаю туда точку.
1.
2.
3.
4.
5.
CREATE TABLE IF NOT EXISTS `geom` (
`g` point DEFAULT NULL
) ENGINE=MyISAM DEFAULT CHARSET=utf8;
INSERT INTO geom VALUES (GeomFromText('POINT(1 1)'));
у меня есть многоугольник (полигон), я хочу узнать входит ли это точка в мой многоугольник
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
DELIMITER //
CREATE FUNCTION GISWithin(pt POINT, mp MULTIPOLYGON) RETURNS INT(1) DETERMINISTIC
BEGIN
DECLARE str, xy TEXT;
DECLARE x, y, p1x, p1y, p2x, p2y, m, xinters DECIMAL(16, 13) DEFAULT 0;
DECLARE counter INT DEFAULT 0;
DECLARE p, pb, pe INT DEFAULT 0;
SELECT MBRWithin(pt, mp) INTO p;
IF p != 1 OR ISNULL(p) THEN
RETURN p;
END IF;
SELECT X(pt), Y(pt), ASTEXT(mp) INTO x, y, str;
SET str = REPLACE(str, 'POLYGON((','');
SET str = REPLACE(str, '))', '');
SET str = CONCAT(str, ',');
SET pb = 1;
SET pe = LOCATE(',', str);
SET xy = SUBSTRING(str, pb, pe - pb);
SET p = INSTR(xy, ' ');
SET p1x = SUBSTRING(xy, 1, p - 1);
SET p1y = SUBSTRING(xy, p + 1);
SET str = CONCAT(str, xy, ',');
WHILE pe > 0 DO
SET xy = SUBSTRING(str, pb, pe - pb);
SET p = INSTR(xy, ' ');
SET p2x = SUBSTRING(xy, 1, p - 1);
SET p2y = SUBSTRING(xy, p + 1);
IF p1y < p2y THEN SET m = p1y; ELSE SET m = p2y; END IF;
IF y > m THEN
IF p1y > p2y THEN SET m = p1y; ELSE SET m = p2y; END IF;
IF y <= m THEN
IF p1x > p2x THEN SET m = p1x; ELSE SET m = p2x; END IF;
IF x <= m THEN
IF p1y != p2y THEN
SET xinters = (y - p1y) * (p2x - p1x) / (p2y - p1y) + p1x;
END IF;
IF p1x = p2x OR x <= xinters THEN
SET counter = counter + 1;
END IF;
END IF;
END IF;
END IF;
SET p1x = p2x;
SET p1y = p2y;
SET pb = pe + 1;
SET pe = LOCATE(',', str, pb);
END WHILE;
RETURN counter % 2;
END;
DELIMITER ;
SET @poly = 'Polygon((0 0,0 3,3 0,0 0),(1 1,1 2,2 1,1 1))';
SELECT AsText(g) FROM `geom` WHERE GISWithin(g, @poly);
все вроде хорошо, но делаю так и ругается на синтаксис...
#1064 - You have an error in your SQL syntax; check the manual that corresponds to your MariaDB server version for the right syntax to use near 'SELECT AsText(g) FROM `geom` WHERE GISWithin(g, @poly)' at line 2
не могу понять что он хочет...
если перезапустить сервер и делать просто
1.
'SELECT AsText(g) FROM `geom`
то этот простой запрос работает
никогда не имел дело с подобным, НО
просто зайдя на страницу документации понял
1)функция ждёт mp MULTIPOLIGON
2)такой тип строиться путём MULTIPOLYGON(poly1,poly2)
3)poly1 = POLYGON(ls11,ls12), poly2 = POLYGON(ls21,ls22)
4)ls11 = LINESTRING(pt111,pt112,pt113), ls12.......
5)pt111 = POINT(0 0);.....
Тоесть задание мультиполигона с полигоном типо кдвадрата 00 11 01 10
будет
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
SET @poly =
MULTIPOLYGON(
POLYGON(
LINESTRING(
POINT(0,0),
POINT(0,3),
POINT(3,3),
POINT(3,0)
),
LINESTRING(
POINT(0,0),
POINT(3,0)
)
),
POLYGON(
LINESTRING(
POINT(1,1),
POINT(1,3),
POINT(3,3),
POINT(3,1)
),
LINESTRING(
POINT(1,1),
POINT(3,1)
)
)
);
