Leaflet多段线点击定位：如何在最近点两侧识别点击所在线段

<?php

/**
 * 计算两个地理坐标点之间的方位角（Bearing）。
 * 方位角是从第一个点指向第二个点的方向，以度为单位（0-360）。
 * 北方为0度，东方为90度，南方为180度，西方为270度。
 *
 * @param float $lat1 起点纬度 (度)
 * @param float $lon1 起点经度 (度)
 * @param float $lat2 终点纬度 (度)
 * @param float $lon2 终点经度 (度)
 * @return float 方位角 (度, 0-360)
 */
function calculateBearing($lat1, $lon1, $lat2, $lon2) {
    // 将度转换为弧度
    $lat1 = deg2rad($lat1);
    $lon1 = deg2rad($lon1);
    $lat2 = deg2rad($lat2);
    $lon2 = deg2rad($lon2);

    // 计算经度差
    $dLon = $lon2 - $lon1;

    // 使用atan2函数计算方位角
    $y = sin($dLon) * cos($lat2);
    $x = cos($lat1) * sin($lat2) - sin($lat1) * cos($lat2) * cos($dLon);

    $bearing = atan2($y, $x);
    // 将弧度转换为度
    $bearing = rad2deg($bearing);
    // 确保方位角在0到360度之间
    $bearing = fmod(($bearing + 360), 360);

    return $bearing;
}

// 示例用法：
// 假设多段线顶点 B: (51.6898945656, -2.0241979535)
// 假设多段线顶点 C: (51.6909727036, -2.0277718088)
// 假设点击点: (51.690195, -2.025175)
// 假设多段线顶点 A: (51.689000, -2.020000)

$pointB_lat = 51.6898945656; $pointB_lon = -2.0241979535;
$pointC_lat = 51.6909727036; $pointC_lon = -2.0277718088;
$click_lat   = 51.690195;    $click_lon   = -2.025175;
$pointA_lat  = 51.689000;    $pointA_lon  = -2.020000;

// 计算 B 到 C 的方位角
$bearing_B_to_C = calculateBearing($pointB_lat, $pointB_lon, $pointC_lat, $pointC_lon);
echo "B到C的方位角: " . round($bearing_B_to_C, 2) . "度\n"; // 示例输出: 304.18度 (或根据实际坐标)

// 计算点击点到 B 的方位角
$bearing_click_to_B = calculateBearing($click_lat, $click_lon, $pointB_lat, $pointB_lon);
echo "点击点到B的方位角: " . round($bearing_click_to_B, 2) . "度\n"; // 示例输出: 298.62度 (或根据实际坐标)

// 计算 A 到 B 的方位角
$bearing_A_to_B = calculateBearing($pointA_lat, $pointA_lon, $pointB_lat, $pointB_lon);
echo "A到B的方位角: " . round($bearing_A_to_B, 2) . "度\n"; // 示例输出: 303.81度 (或根据实际坐标)

// 判断逻辑 (伪代码):
$tolerance = 5; // 容忍度，例如5度
if (abs($bearing_click_to_B - $bearing_A_to_B) < $tolerance || abs($bearing_click_to_B - ($bearing_A_to_B + 180) % 360) < $tolerance) {
    echo "点击点可能在A-B线段附近。\n";
} elseif (abs($bearing_click_to_B - $bearing_B_to_C) < $tolerance || abs($bearing_click_to_B - ($bearing_B_to_C + 180) % 360) < $tolerance) {
    echo "点击点可能在B-C线段附近。\n";
} else {
    echo "无法明确判断线段，可能点击点离线段较远或在其他复杂区域。\n";
}

?>