0

以下の関数は、4 つの点から多角形の正方形を作成します。最後の 5 番目の点が正方形を閉じていると仮定します。少なくとも 24 個の点があり、円を構成する必要があります。何か案は?

function getRadius($point="POINT(-29.8368 30.9096)", $radius=2)
{
    $km = 0.009;
    $center = "GeomFromText('$point')";
        $radius = $radius*$km;
        $bbox = "CONCAT('POLYGON((',
                X($center) - $radius, ' ', Y($center) - $radius, ',',
                X($center) + $radius, ' ', Y($center) - $radius, ',',
                X($center) + $radius, ' ', Y($center) + $radius, ',',
                X($center) - $radius, ' ', Y($center) + $radius, ',',
                X($center) - $radius, ' ', Y($center) - $radius, '
        ))')";

    $query = $this->db->query("
        SELECT id, AsText(latLng) AS latLng, (SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )))/0.009 AS distance
        FROM crime_listing
        WHERE Intersects( latLng, GeomFromText($bbox) )
        AND SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )) < $radius
        ORDER BY distance
                ");

    if($query->num_rows()>0){
                return($query->result());
        }else{
                return false;
        }
}

jsバージョンの下で、これは完璧に機能します

var findCirclePolygons = function(point, r)
{   
    var d2r = Math.PI / 180;
    this.circleLatLngs = new Array();
    numPoints = 24;
    var circleLat = r * 0.009;  // Convert degrees into km
    var circleLng = circleLat / Math.cos(point.lat() * d2r);
    for (var i = 0; i < numPoints + 1; i++) {
        var theta = Math.PI * (i / (numPoints / 2));
        var vertexLat = point.lat() + (circleLat * Math.sin(theta));
        var vertexLng = parseFloat(point.lng()) + parseFloat(( circleLng * Math.cos(theta)));
        var vertextLatLng = new google.maps.LatLng(vertexLat, vertexLng);
        this.circleLatLngs.push(vertextLatLng);
    }

    // Set options
    var options = {
        paths: circleLatLngs,
        strokeColor: "#0055ff",
        strokeOpacity: 1,
        strokeWeight: 1,
        fillColor: "#0055ff",
        fillOpacity: 0.35
        };

    // Return
    return options;
};
4

1 に答える 1

3

ループを作成できます (疑似コードで):

for(i=0; i <= 360; i += 360/24)
{
    $extra_point = "POINT(". $radius*cos(i)." ". $radius*sin(i) . ")"
}

編集

  • ただし、境界ボックスだけが必要な場合は、正方形を使用するのが正しい方法です。円に見えるバウンディング ボックスは必要ありません。後で並べ替えで除外します。
  • 距離を測定するために、ABS を使用する必要はありません。二乗はあなたのためにそれを行います(パフォーマンスについて非常にうるさい場合は、平方根を計算する代わりに半径の二乗と比較してください)
  • 空間インデックスを追加することを忘れないでください。そうしないと、パフォーマンスががらくたになります

2番目の編集 やや少ない疑似コードで;)(私はphpを使用していない時代がありました)。

$lon = 42;
$lat = 2;
$radius = 0.01;
$bbox = "POLYGON((";
for(i=0; i <= 360; i += 360/24)
{
    $bbox .= $radius*cos(deg2rad(i)) + $lon." ". $radius*sin(deg2rad(i)) + $lat;
    if(i < 360)"
        $bbox .= ", "
}
$bbox .= "))"

距離については、任意の x に対して pow(abs(x),2) == pow(x, 2) があります。したがって、(簡単にするために)次のように書くことができます:

SQRT(POW( X(latLng) - X({$center}), 2) + POW(Y(latLng) - Y({$center}), 2 )) < $radius

または、平方根の計算を節約することもできます(この場合、何も変更されず、初期の最適化は悪いことなので、好きなようにしてください)

POW( X(latLng) - X({$center}), 2) + POW(Y(latLng) - Y({$center}), 2 ) < POW($radius, 2)
于 2009-12-01T15:46:25.053 に答える