在Ruby on Rails中使用RGeo计算经纬度多边形的正确面积

问题分析 #

• 使用RGeo::Geographic.spherical_factory时，SphericalPolygonImpl类未定义area方法，因为球面几何实现未内置面积计算逻辑。
• 直接用RGeo::Cartesian.factory(srid:4326)计算的是平方度，不是实际地表面积——经纬度是地理坐标系，不同纬度的1度对应的实际距离差异极大，结果必然不准确。

解决方案：转换到平面投影坐标系计算 #

要得到准确的面积，需将地理坐标系（WGS84，EPSG:4326）的多边形转换到平面投影坐标系（如UTM、当地高斯-克吕格坐标系），这类坐标系以米为单位，能正确计算实际地表面积。

方法1：手动投影转换

先构建球面多边形，再通过Proj4转换到目标投影坐标系后计算面积：

require 'rgeo'
require 'rgeo/proj4'

# 1. 创建WGS84地理坐标系factory
geo_factory = RGeo::Geographic.spherical_factory(srid: 4326)

# 2. 构建原始球面多边形
ring = geo_factory.linear_ring(self.get_polygon.map { |lon, lat| geo_factory.point(lon, lat) })
polygon = geo_factory.polygon(ring)

# 3. 定义目标投影（示例：UTM 31N，覆盖经度0°-6°E，单位米）
# 请根据你的地块位置替换为对应UTM带的proj4参数
target_proj = RGeo::Proj4.new('+proj=utm +zone=31 +datum=WGS84 +units=m +no_defs')

# 4. 转换多边形到平面投影
projected_polygon = RGeo::Proj4.transform(polygon, target_proj)

# 5. 计算面积（单位：平方米）
area = projected_polygon.area

# 可选：转换为平方公里
# area_sqkm = area / 1_000_000

方法2：使用投影地理Factory

直接创建支持自动投影的Factory，无需手动转换，内部会处理地理坐标到平面投影的转换：

require 'rgeo'
require 'rgeo/proj4'

# 创建投影地理Factory，自动处理坐标转换
# 需根据地块位置替换zone参数（UTM带号），南半球将'N'改为'S'
factory = RGeo::Geographic.projected_factory(
  srid: 4326,
  projection_srid: 32631, # UTM 31N对应的EPSG代码
  projection_proj4: '+proj=utm +zone=31 +datum=WGS84 +units=m +no_defs',
  uses_projection: true
)

# 构建多边形并直接计算面积（单位：平方米）
ring = factory.linear_ring(self.get_polygon.map { |lon, lat| factory.point(lon, lat) })
polygon = factory.polygon(ring)
area = polygon.area

如何选择UTM带号 #

UTM坐标系将地球分为60个经度带，每个带覆盖6°经度：

1. 计算带号：带号 = ((经度 + 180) / 6).floor + 1
2. 北半球后缀用N，南半球用S

例如：

• 经度116°E：((116+180)/6).floor +1 = (296/6).floor+1=49+1=50，对应UTM 50N

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