本文详解如何解析 admiralty tide api 的 geojson 格式站点数据，使用球面余弦定理精确计算两点间地理距离，并快速定位最近站点的 id、名称等属性。
在处理地理坐标数据（如潮汐站、气象站或地图标记）时，仅用平面直角坐标系（如勾股定理）计算距离会导致显著误差——尤其当跨度超过几十公里或靠近极地时。正确的做法是采用球面三角法（Haversine 或余弦公式），将地球近似为球体，基于经纬度计算大圆距离（Great-circle distance）。
以下是一个完整、健壮、可直接部署的 PHP 教程实现，适用于 https://easytide.admiralty.co.uk/Home/GetStations 这类标准 GeoJSON 数据源：
✅ 正确的距离计算函数（支持海里/公里/英里）
function getDistance($from_lat, $from_lng, $to_lat, $to_lng, $unit = 'nmi', $decimals = 2) {
// 将角度转为弧度，代入球面余弦公式
$lat1 = deg2rad($from_lat);
$lng1 = deg2rad($from_lng);
$lat2 = deg2rad($to_lat);
$lng2 = deg2rad($to_lng);
$degrees = rad2deg(acos(
sin($lat1) * sin($lat2) +
cos($lat1) * cos($lat2) * cos($lng2 - $lng1)
));
// 按单位换算：1度 ≈ 111.134 km（平均地球周长 / 360）
$conversion = match($unit) {
'km' => 111.13384,
'mi' => 69.05482,
'nmi' => 59.97662,
default => 59.97662
};
return round($degrees * $conversion, $decimals);
}⚠️ 注意：PHP 8.0+ 支持 match 表达式；若使用 PHP < 8.0，请替换为 switch 语句（见原答案）。
✅ 解析 GeoJSON 并查找最近站点
// 1. 获取并解码远程 JSON（生产环境建议加超时与错误处理）
$json_url = 'https://easytide.admiralty.co.uk/Home/GetStations';
$json_content = @file_get_contents($json_url, false, stream_context_create([
'http' => ['timeout' => 10]
]));
if (!$json_content) {
die("❌ 无法获取潮汐站数据，请检查网络或 URL");
}
$stations = json_decode($json_content);
if (json_last_error() !== JSON_ERROR_NONE) {
die("❌ JSON 解析失败：" . json_last_error_msg());
}
// 2. 设置参考位置（例如：朴茨茅斯城堡 GPS 坐标）
$base = [
'lat' => 50.77842324616663,
'lng' => -1.087804949548603
];
// 3. 遍历 features 数组，逐个计算距离并追踪最小值
$closest = null;
$min_distance = null;
foreach ($stations->features as $index => $feature) {
// GeoJSON 坐标格式为 [lng, lat] → 注意顺序！
if (!isset($feature->geometry->coordinates[0], $feature->geometry->coordinates[1])) {
continue; // 跳过无效坐标
}
$lng = (float)$feature->geometry->coordinates[0];
$lat = (float)$feature->geometry->coordinates[1];
$distance = getDistance($base['lat'], $base['lng'], $lat, $lng, 'nmi', 4);
if ($min_distance === null || $distance < $min_distance) {
$min_distance = $distance;
$closest = [
'index' => $index,
'distance' => $distance,
'id' => $feature->properties->Id ?? 'N/A',
'name' => $feature->properties->Name ?? 'Unknown',
'country' => $feature->properties->Country ?? 'N/A',
'raw' => $feature // 可选：保留完整对象供后续扩展
];
}
}
// 4. 输出结果（含关键 ID，可用于后续 API 请求）
if ($closest) {
echo "✅ 最近站点：{$closest['name']} (ID: {$closest['id']})\n";
echo "? 距离：{$closest['distance']} 海里\n";
echo "? 国家：{$closest['country']}\n";
// ? 后续可直接用于其他请求，例如：
// $tide_data = file_get_contents("https://easytide.admiralty.co.uk/Home/GetTideData?stationId={$closest['id']}");
} else {
echo "⚠️ 未找到有效站点。\n";
}? 关键要点总结
- 坐标顺序陷阱：GeoJSON 规范中 coordinates 是 [longitude, latitude]，而非 [lat, lng] —— 错位将导致完全错误的结果；
- 空值防护：务必检查 ->geometry->coordinates 和 ->properties->Id 是否存在，避免 Trying to access property on null；
- 性能优化：无需预先存储所有距离，单次遍历 + 实时比较即可，空间复杂度 O(1)；
-
生产就绪建议：
- 使用 cURL 替代 file_get_contents，便于设置 User-Agent、重试、SSL 验证；
- 对远程接口做缓存（如 Redis 或文件缓存），避免频繁请求；
- 添加 try/catch 包裹 json_decode 和网络调用；
- 若需高精度（如航海级），可升级为 Vincenty 公式（但余弦法对绝大多数应用已足够）。
掌握这一模式后，你不仅能对接 Admiralty Tide API，还可轻松适配 OpenStreetMap、GeoNames、任何返回 GeoJSON 的地理服务——真正实现「一次封装，多处复用」。
