一、通用计算思路 #

如果不想依赖第三方API，你可以通过以下步骤自行实现：

• 获取海岸线数据集：比如使用OpenStreetMap导出的全球海岸线shapefile数据，或是NOAA提供的高精度海岸线数据集，这类数据能覆盖全球绝大多数海岸区域。
• 空间查询与距离计算：把目标点和海岸线数据导入支持空间计算的数据库（比如PostGIS），利用内置的空间函数（如ST_Distance）直接计算点到最近海岸线要素的最短距离。这种方法适合离线场景或对精度要求极高的需求，但缺点是需要维护庞大的数据集，且对空间数据库的使用有一定技术门槛。

如果追求快速落地，直接调用第三方地理API是更省心的选择，下面重点讲Google Maps生态下的PHP实现方案。

二、Google Maps API的PHP实现方案 #

首先明确：Google Maps没有专门的「计算到最近海岸线距离」的API端点，但我们可以组合现有API来实现需求——核心思路是：搜索目标点附近的海岸线/沿海自然特征点，再计算目标点到这些点的最短距离。

步骤1：准备前置条件 #

你需要先在Google Cloud控制台创建项目，启用Places API和Geocoding API，并生成API密钥（记得给密钥设置使用限制，避免被盗用产生不必要的费用）。

步骤2：完整PHP代码示例 #

下面是一个可直接运行的示例，用cURL调用Places API搜索附近海岸线，再通过Haversine公式计算最短距离：

<?php
// 目标地点的经纬度（示例：北京）
$targetLat = 39.9042;
$targetLng = 116.4074;
$apiKey = '你的Google Maps API密钥';

// 1. 调用Places API搜索附近的海岸线自然特征
$placesUrl = "https://maps.googleapis.com/maps/api/place/nearbysearch/json?location={$targetLat},{$targetLng}&radius=50000&type=natural_feature&keyword=coastline&key={$apiKey}";

$ch = curl_init($placesUrl);
curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
$placesResponse = curl_exec($ch);
curl_close($ch);

$placesData = json_decode($placesResponse, true);

if ($placesData['status'] !== 'OK') {
    die("Places API调用失败: " . $placesData['status']);
}

// 2. 计算目标点到每个搜索结果的距离，取最小值
$minDistance = PHP_INT_MAX;
$closestCoastPoint = null;

// Haversine公式：计算两个经纬度点的直线距离（单位：千米）
function calculateDistance($lat1, $lng1, $lat2, $lng2) {
    $earthRadius = 6371; // 地球平均半径（千米）
    $dLat = deg2rad($lat2 - $lat1);
    $dLng = deg2rad($lng2 - $lng1);
    $a = sin($dLat/2) * sin($dLat/2) +
         cos(deg2rad($lat1)) * cos(deg2rad($lat2)) *
         sin($dLng/2) * sin($dLng/2);
    $c = 2 * atan2(sqrt($a), sqrt(1-$a));
    return $earthRadius * $c;
}

foreach ($placesData['results'] as $place) {
    $coastLat = $place['geometry']['location']['lat'];
    $coastLng = $place['geometry']['location']['lng'];
    $distance = calculateDistance($targetLat, $targetLng, $coastLat, $coastLng);
    
    if ($distance < $minDistance) {
        $minDistance = $distance;
        $closestCoastPoint = [
            'lat' => $coastLat,
            'lng' => $coastLng,
            'name' => $place['name']
        ];
    }
}

// 输出结果
if ($closestCoastPoint) {
    echo "到最近海岸线的距离：" . round($minDistance, 2) . " 千米\n";
    echo "最近海岸线点的经纬度：{$closestCoastPoint['lat']}, {$closestCoastPoint['lng']}\n";
    echo "地点名称：{$closestCoastPoint['name']}\n";
} else {
    echo "未找到附近的海岸线，请尝试扩大搜索半径。";
}
?>

关键细节说明 #

• 搜索参数调整：radius参数可根据目标点位置调整（内陆地区建议设为50000即50千米以上，沿海地区可缩小）；keyword=coastline确保我们精准定位海岸线相关的自然特征。
• 单位转换：如果需要英里为单位，只需把Haversine公式中的$earthRadius改为3956（地球平均半径英里数）即可。
• 精度说明：Places API返回的是海岸线附近的特征点，并非精确的海岸线线段，所以计算结果是近似值。如果需要更高精度，建议结合OpenStreetMap的离线数据进行空间分析。

三、补充：快速判断是否在海边 #

如果目标点本身就在海岸线上，你可以先调用Reverse Geocoding API验证，避免不必要的搜索：

// 调用Reverse Geocoding API检查是否在海岸线上
$geocodeUrl = "https://maps.googleapis.com/maps/api/geocode/json?latlng={$targetLat},{$targetLng}&key={$apiKey}";
$ch = curl_init($geocodeUrl);
curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
$geocodeResponse = curl_exec($ch);
curl_close($ch);
$geocodeData = json_decode($geocodeResponse, true);

foreach ($geocodeData['results'][0]['address_components'] as $component) {
    if (in_array('natural_feature', $component['types']) && $component['long_name'] === 'Coastline') {
        echo "当前地点就在海岸线上，距离为0千米。";
        exit;
    }
}

内容的提问来源于stack exchange，提问作者test1gddg