14

グーグルマップ上の緯度/経度の値とポリゴンのポイントインポリゴンの迅速でかなり正確な方法を理解する必要があります。いくつかの調査の後-mysqlの幾何学的拡張に関するいくつかの投稿に出くわし、それも実装しました-

SELECT id, Contains( PolyFromText( 'POLYGON(".$polygonpath.")' ) , PointFromText( concat( \"POINT(\", latitude, \" \", longitude, \")\" ) ) ) AS
            CONTAINS
FROM tbl_points

ただし、これは、多数のポイントで構成されるポリゴンでは機能しませんでした:(

さらに調査した後、レイキャスティングアルゴリズムと呼ばれる標準アルゴリズムに出くわしましたが、MySQLでそのクエリを開発する前に、誰かがすでにそれを経験しているか、実装方法を示す便利なリンクに出くわした場合は、チャンスをつかみたいと思いましたMySQL/SQLサーバーのアルゴリズム。

だから、それを短くする-質問は：

誰かがレイキャスティングアルゴリズムのMySQL/SQLサーバー実装を提供できますか？

追加の詳細：

  • ポリゴンは、凹面、凸面、または複雑のいずれかです。
  • 100％以上の精度で迅速な実行を目標としています。
4

7 に答える 7

18

次の関数（RaycastingアルゴリズムのMYSQLバージョン）は私の世界を揺るがしました：

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 myWithin(point, polygon) as result;

どこ

 point  = Point(lat,lng) 
 polygon = Polygon(lat1 lng1, lat2 lng2, lat3 lng3, .... latn lngn, lat1 lng1)

ポリゴンを閉じる必要があることに注意してください（通常、標準のkmlまたはgooglemapデータを取得する場合は閉じますが、必ず閉じてください。最後にlat1 lng1セットが繰り返されることに注意してください）。

データベースに幾何学的フィールドとしてのポイントとポリゴンがなかったため、次のような操作を行う必要がありました。

 Select myWithin(PointFromText( concat( "POINT(", latitude, " ", longitude, ")" ) ),PolyFromText( 'POLYGON((lat1 lng1, ..... latn lngn, lat1 lng1))' ) ) as result

これが誰かに役立つことを願っています。

于 2011-10-10T06:17:44.903 に答える
5

CやDelphi、または使用する高レベルのツールでレイキャスティングアルゴリズムを実装するカスタムUDFを作成します。

UDFを作成するためのリンク
球上のポイントを検索するMySQLGIS実装のソースコードを次に示します（MySQLの空間データ型と対話する方法を確認するためのテンプレートとして使用してください）。
http://www.lenzg.net/archives/220-New-UDF-for-MySQL-5.1-provides-GIS-functions-distance_sphere-and-distance_spheroid.html

MySQLマニュアルから：http：
//dev.mysql.com/doc/refman/5.0/en/adding-functions.html

MS VisualC++のUDFチュートリアル
http://rpbouman.blogspot.com/2007/09/creating-mysql-udfs-with-microsoft.html

DelphiでのUDFチュートリアル：Delphi
でのMySQL用のUDFの作成

レイキャスティングアルゴリズムに関するソースコード
擬似コード：http
： //rosettacode.org/wiki/Ray-casting_algorithm drDobbsの記事（記事の上部にあるコードへのリンクに注意してください）：http ：//drdobbs.com / cpp / 184409586
Delphi（実際にはFreePascal）： http: //www.cabiatl.com/mricro/raycast/

于 2011-10-03T20:57:53.983 に答える
2

念のため、入力としてMULTIPOLYGONを受け入れる1つのMySQL関数：http：//forums.mysql.com/read.php？23,286574,286574

DELIMITER $$

CREATE DEFINER=`root`@`localhost` 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
于 2012-11-02T13:52:01.943 に答える
2

ポリゴン内の緯度/経度を検索するためのzarun関数への応答。

緯度/経度の情報を持つプロパティテーブルがありました。そのため、lat/longがポリゴンlats/longs内にあるレコードを取得する必要がありました（Google APIから取得しました）。最初は、Zarun関数の使い方がばかげていました。これがそのソリューションクエリです。

  • 表：プロパティ
  • フィールド： id、緯度、経度、ベッドなど...
  • クエリ：
SELECT id
FROM properties
WHERE myWithin(
    PointFromText(concat( "POINT(", latitude, " ", longitude, ")")),
    PolyFromText('POLYGON((37.628134 -77.458334,37.629867 -77.449021,37.62324 -77.445416,37.622424 -77.457819,37.628134 -77.458334))' )
) = 1 limit 0,50;

それが私のような愚か者のための時間を節約することを願っています;）

于 2013-10-26T17:38:30.027 に答える
1

MySQL5.6.1以降では、SpatialExtensionになりました。ドキュメントのfunction_st-containsを参照してください。

于 2013-02-19T10:35:28.903 に答える
1

上記のmywithinストアドプロシージャをポリゴンのテーブルで使用したかったので、これを実行するコマンドを次に示します。

次のようにogr2ogrを使用して、ポリゴンを含むシェープファイルをmysqlにインポートした後

ogr2ogr -f "mysql" MYSQL:"mydbname,host=localhost,user=root,password=mypassword,port=3306" -nln "mytablename" -a_srs "EPSG:4326" /path/to/shapefile.shp

次に、MBRwithinを使用してテーブルを事前フィルタリングし、mywithinを使用して次のように終了できます。

DROP TEMPORARY TABLE IF EXISTS POSSIBLE_POLYS;
CREATE TEMPORARY TABLE POSSIBLE_POLYS(OGR_FID INT,SHAPE POLYGON);
INSERT INTO POSSIBLE_POLYS (OGR_FID, SHAPE) SELECT mytablename.OGR_FID,mytablename.SHAPE FROM mytablename WHERE MBRwithin(@testpoint,mytablename.SHAPE);

DROP TEMPORARY TABLE IF EXISTS DEFINITE_POLY;
CREATE TEMPORARY TABLE DEFINITE_POLY(OGR_FID INT,SHAPE POLYGON);
INSERT INTO DEFINITE_POLY (OGR_FID, SHAPE) SELECT POSSIBLE_POLYS.OGR_FID,POSSIBLE_POLYS.SHAPE FROM POSSIBLE_POLYS WHERE mywithin(@testpoint,POSSIBLE_POLYS.SHAPE);

@testpointは、たとえば、

SET @longitude=120;
SET @latitude=-30;
SET @testpoint =(PointFromText( concat( "POINT(", @longitude, " ", @latitude, ")" ) ));
于 2012-06-07T21:15:36.390 に答える
1

これは、MULTIPOLYGONS（POLYGONSでのみ機能するZarunのバージョンの適応）で動作するバージョンです。

CREATE FUNCTION GISWithin(p POINT, multipoly MULTIPOLYGON) RETURNS INT(1) DETERMINISTIC 
BEGIN 
DECLARE n,i,m,x INT DEFAULT 0; 
DECLARE pX,pY,poly1X,poly1Y,poly2X,poly2Y DECIMAL(13,10); 
DECLARE ls LINESTRING; 
DECLARE poly MULTIPOLYGON;
DECLARE poly1,poly2 POINT; 
DECLARE result INT(1) DEFAULT 0; 
SET pX = X(p); 
SET pY = Y(p); 
SET m = NumGeometries(multipoly);
WHILE x<m DO 
SET poly = GeometryN(multipoly,x);
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; 
SET x = x + 1;
END WHILE;
RETURN result; 
End;
于 2014-05-08T04:18:45.597 に答える