sql >> データベース >  >> RDS >> Mysql

レイキャスティングアルゴリズムのMySQL実装?

    次の関数(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
    

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



    1. Oracleでフィールドの値が変更されたときに行番号をインクリメントする

    2. MySQLトリガーが機能しない、単純な構文、複雑ではない

    3. リレーショナルデータベース設計(MySQL)

    4. MySQLのGROUP_CONCATとLongtextの問題