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。平方会为您完成（作为额外的，如果您对性能非常挑剔，请与半径平方进行比较，而不是计算平方根）
  • 不要忘记添加空间索引，否则表演会很糟糕

第二次编辑 伪代码少一些；）（我已经很久没有使用 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 回答