Sehen Sie, wenn lat / long, falls innerhalb eines Polygons mit mysql

Wenn Sie nicht ändern können dbs eine, die räumliche Operatoren implementiert korrekt, wie PostgreSQL mit der PostGIS-Erweiterung http://postgis.refractions.net/ lösen Sie dieses problem mit einem zweiteiligen Ansatz.

Ersten lassen MySQL geben Sie eine bounding-box pre-filtering Ergebnis basiert auf der bounding-box (das ist das, was Sie standardmäßig tut) mit Ihrem intersects-operator (http://dev.mysql.com/doc/refman/5.1/en/functions-that-test-spatial-relationships-between-geometries.html#function_intersects).
If-Abfragen sind langsam, machen Sie sicher, dass Sie über einen index auf Ihre geometrie-Feld zuerst.

Dann Hydrat der ursprünglichen geometrie, die Sie in Ihre Abfrage in ein geometry-Objekt von GIS geometrie-Bibliothek GEOS (http://trac.osgeo.org/geos/) (C++ basiert, obwohl es auch bindings für andere Sprachen wie Python), Wohlgeformten (http://trac.gispython.org/lab/wiki/Shapely)OGR ( oder der Java Topology Suite (JTS) http://www.vividsolutions.com/jts/jtshome.htm).

Test jeweils die Geometrien, die Sie wieder von Ihrem Abfrage-Ergebnis bei Verwendung der entsprechenden Operatoren wie innerhalb oder schneidet. Jede dieser Bibliotheken wird Ihnen ein Boolesches Ergebnis.

Persönlich würde ich schauen bei den Proben für OGR, da es eine große community, die bereit ist zu helfen.

Ach ja, und sorry für das setzen des links so... ich denke, da bin ich "neu" kann ich nur post einen link (?)

InformationsquelleAutor der Antwort rburhum

Habe ich umgeschrieben Funktion gegeben wurde, dass im vorherigen post von @danherd, so kann es die Arbeit mit echten Multipolygone, die bestehen aus mehr als einem polygon. Für diejenigen von Euch, die
immer noch mit der alten MySql-version sollte es helfen.

Hier ist es:

DELIMITER //

CREATE FUNCTION GISWithin(pt POINT, mp MULTIPOLYGON) RETURNS INT(1) DETERMINISTIC

BEGIN

DECLARE str_big, str, xy LONGTEXT;
DECLARE x, y, p1x, p1y, p2x, p2y, m, xinters DECIMAL(16, 13) DEFAULT 0;
DECLARE counter INT DEFAULT 0;
DECLARE p, pb, pe, sb, se, ct DECIMAL(16, 0) 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_big;
SET str_big = REPLACE(str_big, 'MULTIPOLYGON(((','');
SET str_big = REPLACE(str_big, ')))', '');
SET str_big = REPLACE(str_big, ')),((', '|');
SET str_big = CONCAT(str_big, '|');

SET sb = 1;
SET se = LOCATE('|', str_big);
SET str = SUBSTRING(str_big, sb, se - sb);

WHILE se > 0 DO
SET ct = ct + 1;
SET str = SUBSTRING(str_big, sb, se - sb);

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; 

SET sb = se + 1;
SET se = LOCATE('|', str_big, sb);

END WHILE;

RETURN counter % 2;

END

DELIMITER ; 

InformationsquelleAutor der Antwort user15