Я нашел хранимую процедуру, которая вычисляет, находится ли точка в многоугольнике или вне его. Но у меня небольшие проблемы с адаптацией к моим потребностям.
Хранимая процедура вызывается следующим образом:
SET @point = PointFromText('POINT(5 5)') ;
SET @polygon = PolyFromText('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))');
SELECT myWithin(@point, @polygon) AS result;
Я бы хотел получить широту и долготу из myPointsTable и установить точку @ в результат. Что-то вроде (псевдокод):
SET @point = PointFromText('POINT( SELECT latitude, longitude
FROM `myPointsTable`
WHERE color = 'red')') ;
А затем выполните хранимую процедуру итерации по этим точкам.
Точно так же я хотел бы вытащить полигон из myPolygonTable. Нечто подобное (псевдокод):
SET @polygon = PolyFromText('POLYGON((SELECT polygonPoints FROM 'myPolygonTable' WHERE town = 'Newport'))');
Хранимая процедура:
DROP FUNCTION IF EXISTS myWithin;
DELIMITER $$
CREATE FUNCTION myWithin(p POINT, poly POLYGON) RETURNS INT(1) DETERMINISTIC
BEGIN
DECLARE n INT DEFAULT 0;
DECLARE pX DECIMAL(9,6);
DECLARE pY DECIMAL(9,6);
DECLARE ls LINESTRING;
DECLARE poly1 POINT;
DECLARE poly1X DECIMAL(9,6);
DECLARE poly1Y DECIMAL(9,6);
DECLARE poly2 POINT;
DECLARE poly2X DECIMAL(9,6);
DECLARE poly2Y DECIMAL(9,6);
DECLARE i INT DEFAULT 0;
DECLARE result INT(1) DEFAULT 0;
SET pX = X(p);
SET pY = Y(p);
SET ls = ExteriorRing(poly);
SET poly2 = EndPoint(ls);
SET poly2X = X(poly2);
SET poly2Y = Y(poly2);
SET n = NumPoints(ls);
WHILE i<n DO
SET poly1 = PointN(ls, (i+1));
SET poly1X = X(poly1);
SET poly1Y = Y(poly1);
IF ( ( ( ( poly1X <= pX )
&& ( pX < poly2X ) ) || ( ( poly2X <= pX )
&& ( pX < poly1X ) ) )
&& ( pY > ( poly2Y - poly1Y ) * ( pX - poly1X ) / ( poly2X - poly1X ) + poly1Y ) ) THEN
SET result = !result;
END IF;
SET poly2X = poly1X;
SET poly2Y = poly1Y;
SET i = i + 1;
END WHILE;
RETURN result;
End
$$
DELIMITER ;
Как использовать инструкцию SELECT при установке переменных? Является ли это возможным? Есть ли лучший способ?
Спасибо.
-Laxmidi
UPDATE:
Привет
У меня есть часть проблемы решена. Я понял, как получить широту и долготу.
DROP PROCEDURE IF EXISTS latlongGrabber;
DELIMITER $$
CREATE PROCEDURE latlongGrabber(IN offense_in VARCHAR(255))
BEGIN
DECLARE latitude_val VARCHAR(255);
DECLARE longitude_val VARCHAR(255);
DECLARE no_more_rows BOOLEAN;
DECLARE latlongGrabber_cur CURSOR FOR
SELECT latitude, longitude FROM myTable WHERE offense = offense_in;
DECLARE CONTINUE HANDLER FOR NOT FOUND SET no_more_rows = TRUE;
OPEN latlongGrabber_cur;
the_loop: LOOP
FETCH latlongGrabber_cur
INTO latitude_val, longitude_val;
IF no_more_rows THEN CLOSE latlongGrabber_cur;
LEAVE the_loop;
END IF;
SELECT latitude_val, longitude_val;
END LOOP the_loop;
END
$$
DELIMITER ;
Теперь я должен выяснить, как объединить две функции.
-Laxmidi