本文详解如何解析 admiralty tide api 返回的 geojson 数据，基于给定经纬度计算球面距离，精准定位最近站点并提取其 `id` 字段，适用于潮汐查询、地理服务等实际场景。
在处理地理坐标数据（如 Admiralty 的潮汐站列表）时，简单使用欧氏距离（Pythagorean Theorem）会因地球曲率导致显著误差——尤其当坐标跨纬度较大时。正确做法是采用大圆距离公式（Haversine 或球面余弦定理），它基于地球近似球体模型，单位精度可达海里（nmi）、公里（km）或英里（mi）。
以下是一个完整、健壮、可直接部署的 PHP 教程实现：
✅ 正确解析 GeoJSON 结构
Admiralty API 返回的是标准 GeoJSON FeatureCollection，其核心结构为：
{
"type": "FeatureCollection",
"features": [
{
"type": "Feature",
"geometry": { "type": "Point", "coordinates": [-4.702540, 51.672791] },
"properties": { "Id": "1603", "Name": "TENBY", ... }
}
]
}注意：coordinates 数组顺序为 [longitude, latitude]（先经度后纬度），这与常见 lat/lng 变量命名相反，极易出错——务必确认 $lng = $feature->geometry->coordinates[0]，$lat = $feature->geometry->coordinates[1]。
立即学习“PHP免费学习笔记（深入）”；
✅ 使用球面余弦距离函数（高精度 & 轻量）
以下函数基于球面余弦定理，兼顾精度与性能，支持 nmi（默认）、km、mi 单位：
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)
));
// 按单位换算为实际距离（地球平均半径换算系数）
$conversion = match($unit) {
'km' => 111.13384, // 1° ≈ 111.13 km
'mi' => 69.05482, // 1° ≈ 69.05 mi
'nmi' => 59.97662, // 1° ≈ 59.98 nmi （航海常用）
default => 59.97662
};
return round($degrees * $conversion, $decimals);
}? 提示：PHP 8.0+ 推荐使用 match 表达式；若用旧版 PHP，请替换为 switch。
✅ 主逻辑：单次遍历 + 实时比较（O(n) 时间复杂度）
避免冗余循环和中间数组（如原代码中多层 foreach 嵌套），直接遍历 features，动态维护最小距离及对应站点：
// 获取实时站点数据（生产环境请添加错误处理与缓存）
$json = file_get_contents('https://easytide.admiralty.co.uk/Home/GetStations');
if ($json === false) {
die("❌ Failed to fetch stations from Admiralty API.");
}
$stations = json_decode($json);
// 设置参考位置（例如：朴茨茅斯城堡，用于测试）
$user_location = ['lat' => 50.778423, 'lng' => -1.087805];
$closest = null;
$min_distance = null;
foreach ($stations->features as $index => $feature) {
// ✅ 关键：GeoJSON coordinates = [lng, lat]
$lng = $feature->geometry->coordinates[0];
$lat = $feature->geometry->coordinates[1];
$distance = getDistance(
$user_location['lat'], $user_location['lng'],
$lat, $lng,
'nmi', 4
);
if ($min_distance === null || $distance < $min_distance) {
$min_distance = $distance;
$closest = [
'id' => $feature->properties->Id,
'name' => $feature->properties->Name,
'distance' => $distance,
'unit' => 'nmi',
'index' => $index
];
}
}
// ✅ 输出结果（可用于后续 API 请求）
if ($closest) {
echo "✅ Closest station: {$closest['name']} (ID: {$closest['id']})\n";
echo "? Distance: {$closest['distance']} {$closest['unit']}\n";
echo "? Use ID '{$closest['id']}' for tide data: \n";
echo " file_get_contents('https://easytide.admiralty.co.uk/Home/GetTideTimes?stationId={$closest['id']}')";
} else {
echo "⚠️ No stations found.";
}⚠️ 注意事项与最佳实践
- HTTPS 与错误处理：file_get_contents() 在无响应或超时时会返回 false，务必检查并加入重试或降级逻辑；
- CORS / API 限制：Admiralty 站点接口可能限制请求频率或来源，生产环境建议服务端代理或加缓存（如 Redis 存储 GetStations 结果，TTL 1 小时）；
- 精度取舍：对百公里内应用，Haversine 公式已足够；若需亚米级精度（如测绘），应使用 Vincenty 或 GeographicLib；
- 性能优化：若站点数超万级，可预构建空间索引（如 GeoHash 或 MySQL POINT + ST_Distance_Sphere），但 Admiralty 当前约 300+ 站点，线性扫描完全足够；
- 安全提醒：勿将 file_get_contents() 直接用于用户可控 URL，防止 SSRF；此处 URL 固定，风险可控。
✅ 总结
你真正需要的不是“嵌套循环遍历子数组”，而是：
- 理解 GeoJSON 坐标顺序（[lng, lat]）；
- 选用正确的地理距离算法（拒绝平面直角距离）；
- 一次遍历 + 状态保持（无需存储所有距离，节省内存）；
- 结构化返回结果（含 Id、Name、distance，便于下游调用）。
现在，你已掌握从 Admiralty API 动态获取最近潮汐站 ID 的完整链路——下一步，即可用该 Id 请求详细潮时数据，构建完整的海洋信息服务。
