mysql - guia - qgis manual

Distancia espacial Mysql usando POINT-No funciona

Question

Mi objetivo es utilizar mysql POINT (lat, long) para encontrar entidades cercanas en la base de datos. Estoy tratando de hacer algo como en la parte inferior de este tutorial http://www.scribd.com/doc/2569355/Geo-Distance-Search-with-MySQL . Esto es lo que tengo:

Mesa:

CREATE TABLE mark (
id INT UNSIGNED NOT NULL AUTO_INCREMENT PRIMARY KEY,
name VARCHAR(20) DEFAULT NULL,
loc POINT NOT NULL,
SPATIAL KEY loc (loc)
) ENGINE=MyISAM;

Insertar algunos datos de prueba:

 INSERT INTO mark (loc,name) VALUES (POINT(59.388433,10.415039), ''Somewhere 1'');
INSERT INTO mark (loc,name) VALUES (POINT(63.41972,10.39856), ''Somewhere 2'');

Declarando la función de distancia:

 DELIMITER $$
CREATE FUNCTION `distance`
(a POINT, b POINT)
RETURNS double DETERMINISTIC
BEGIN
RETURN
round(glength(linestringfromwkb(linestring(asbinary(a),
asbinary(b)))));
END $$
DELIMITER; 

Intentando usar la función para buscar ex .:

 SELECT name, distance(mark.loc, GeomFromText( '' POINT(31.5 42.2) '' )) AS cdist
FROM mark
ORDER BY
cdist limit 10;

o:

 SELECT DISTINCT
dest.name,
distance(orig.loc, dest.loc) as sdistance
FROM
mark orig,
mark dest
having sdistance < 10
ORDER BY
sdistance limit 10;

El problema que estoy obteniendo es: ERROR 1367 (22007): valor ''aswkb (a @ 0)'' no geométrico ilegal encontrado durante el análisis, o ERROR 1416 (22003): No se puede obtener el objeto de geometría de los datos que envía al campo GEOMETRÍA

Parece que no puedo entender cómo resolver esto. Lo importante es que la función ''distancia'' se puede usar dinámicamente.

También probé esta solución: encuentra la distancia entre dos puntos en MYSQL. (usando el tipo de datos Point)

Esta es mi versión mysql mysql Ver 14.14 Distrib 5.5.23, para Linux (x86_64) usando readline 5.1

Espero que la experiencia de alguien me pueda ayudar. ¡Aclamaciones!

Answer 1

Así que terminé con esto como una consulta para calcular la distancia, un ejemplo:

 SELECT glength(LineStringFromWKB(LineString(GeomFromText(astext(PointFromWKB(POINT(63.424818,10.402457)))),GeomFromText(astext(PointFromWKB(POINT(663.422238,10.398996)))))))*100
AS distance;

Lo estoy multiplicando por 100 para obtener una aproximación en kilómetros. El resultado no es exacto, sino "bien". Si alguien conoce una mejor manera, no dude en comentar.

Answer 2

definir una función personalizada

CREATE DEFINER=`test`@`%` FUNCTION `geoDistance`(`lon1` DOUBLE, `lat1` DOUBLE, `lon2` DOUBLE, `lat2` DOUBLE)
RETURNS double
LANGUAGE SQL
DETERMINISTIC
NO SQL
SQL SECURITY DEFINER
COMMENT ''''
BEGIN
DECLARE v DOUBLE;
SELECT cos(radians(lat1))
* cos(radians(lat2))
* cos(radians(lon2) - radians(lon1))
+ sin(radians(lat1))
* sin(radians(lat2)) INTO v;
RETURN IF(v > 1, 0, 6371000 * acos(v));
END

luego llame

SELECT geoDistance(X(point1), Y(point1), X(spoint2), Y(point2))

el resultado viene en metros