Voir si lat / long des chutes à l'intérieur d'un polygone à l'aide de mysql

Si vous ne pouvez pas modifier dbs à celui qui a spatiales les opérateurs mis en œuvre correctement, comme PostgreSQL de l'extension PostGIS http://postgis.refractions.net/ , vous pouvez résoudre ce problème en utilisant une approche à deux volets.

Abord, permettez-MySQL vous donner un cadre de la boîte de pré-filtrage des résultats de la boîte englobante (qu'est ce qu'il fait par défaut) à l'aide de leur coupe de l'opérateur (http://dev.mysql.com/doc/refman/5.1/en/functions-that-test-spatial-relationships-between-geometries.html#function_intersects).
Si les requêtes sont lents, assurez-vous que vous avez un indice sur votre géométrie premier champ.

Puis hydrater la géométrie d'origine que vous avez utilisé dans votre requête dans une géométrie de l'objet de SIG de la géométrie de la bibliothèque comme GEOS (http://trac.osgeo.org/geos/) (C++ en fonction, bien qu'il ait également des liaisons pour les différents langages comme Python), Galbées (http://trac.gispython.org/lab/wiki/Shapely), OGR ( ou le Java Topologie de Suite (JTS) http://www.vividsolutions.com/jts/jtshome.htm).

Essai de chacune des géométries que vous revenez de votre résultat de la requête à l'aide de l'opérateur approprié comme dans ou coupe. L'une de ces bibliothèques vous donnera un résultat booléen.

Personnellement, je regarde les échantillons pour OGR car il a une grosse communauté qui est prêt à aider.

Ah oui, et désolé de mettre les liens comme ça... je suppose que depuis que je suis "nouveau", je ne peux que poster un lien (?)

OriginalL'auteur rburhum

J'ai réécrit la fonction qui a été donné dans le post précédent par @danherd, donc il peut travailler avec de vrais multipolygons qui se composent de plus qu'un polygone. Pour ceux d'entre vous qui
toujours garder l'aide d'une vieille version de MySql, il devrait aider.

Ici, il est:

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 ; 

OriginalL'auteur user15