本文详解如何解析 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 请求详细潮时数据，构建完整的海洋信息服务。