基于Ruby on Rails的酒店应用：rgeo计算酒店到海岸线距离问题

核心问题分析 #

你遇到的问题本质是坐标系与几何计算模型不匹配：

• simple_mercator_factory是平面投影模型，计算的是平面距离，高纬度区域会严重失真，完全不适合球面距离计算；
• spherical_factory基于球面几何模型，两点距离计算正确，但点到线/多线串的距离需要注意几何类型处理、单位转换，以及确保所有几何对象来自同一factory。

正确实现步骤 #

1. 选择合适的Factory

使用基于WGS84（SRID 4326）的球面Factory，这是GeoJSON默认的坐标系，适配绝大多数地理数据：

factory = RGeo::Geographic.spherical_factory(srid: 4326)

如果需要更高精度的椭球计算（比如专业地理场景），可以替换为椭球Factory：

factory = RGeo::Geographic.geographic_factory(srid: 4326)

2. 正确解析海岸线GeoJSON

必须用指定的Factory解析GeoJSON，避免跨Factory的几何对象计算错误：

require 'rgeo/geo_json'

# 读取GeoJSON文件或直接传入字符串
coastline_geojson = File.read('path/to/your/coastline.geojson')
coastline_geometry = RGeo::GeoJSON.decode(coastline_geojson, factory: factory)

3. 计算点到海岸线的最小距离

海岸线通常是MultiLineString类型，需要遍历所有子线段计算最小距离，同时注意将球面弧度转换为实际距离（地球平均半径取6371000米）：

def distance_to_coast(hotel_point, coastline_geometry, factory)
  earth_radius_m = 6371000
  min_distance = Float::INFINITY

  # 确保酒店点与海岸线使用同一Factory
  hotel_point = factory.point(hotel_point.x, hotel_point.y) unless hotel_point.factory == factory

  case coastline_geometry
  when RGeo::Feature::LineString
    # 直接计算点到单一线串的距离
    distance_rad = coastline_geometry.distance(hotel_point)
    min_distance = distance_rad * earth_radius_m
  when RGeo::Feature::MultiLineString
    # 遍历所有子线串，取最小距离
    coastline_geometry.each do |line_string|
      distance_rad = line_string.distance(hotel_point)
      actual_distance = distance_rad * earth_radius_m
      min_distance = actual_distance if actual_distance < min_distance
    end
  else
    raise "不支持的几何类型: #{coastline_geometry.class.name}"
  end

  min_distance / 1000 # 转换为公里，如需米则去掉此步
end

4. 调用示例

# 创建酒店点（注意GeoJSON是经纬度顺序：lng, lat）
hotel_point = factory.point(121.4737, 31.2304) # 以上海为例

# 计算距离
coast_distance_km = distance_to_coast(hotel_point, coastline_geometry, factory)
puts "酒店到海洋的距离：#{coast_distance_km.round(2)}公里"

关键注意事项 #

• 坐标系一致性：所有几何对象（酒店点、海岸线）必须使用同一Factory创建/解析，否则计算会出错；
• 单位转换：spherical_factory的distance方法返回的是弧度值，必须乘以地球半径才能得到实际米数；
• 数据优化：如果海岸线数据量极大，建议先通过边界框（Bounding Box）过滤掉与酒店距离遥远的线段，减少计算量；
• 验证测试：用已知距离的点（比如沿海已知距离的酒店）验证结果，确保计算逻辑正确。

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