ある緯度経度で指定された地点からある距離と方向に移動した後の緯度経度を計算する

gistfile1.rb

here = Position.new(lat_A, lng_B)

# 90度方向に100m移動したあとの地点
there = here.move(distance: 100, heading: 90)

position.rb

class Position
  attr_accessor :latitude, :longitude

  # 地球の半径
  EARTH_RADIUS = 6378150

  def initialize(latitude, longitude)
    @latitude = latitude
    @longitude = longitude
  end

  def move(distance: 0, heading: 0)
    # 緯線上の移動距離
    latitude_distance = distance * Math.cos(heading * Math::PI / 180)

    # 1mあたりの緯度
    earth_circle = 2 * Math::PI * EARTH_RADIUS
    latitude_per_meter = 360 / earth_circle

    # 緯度の変化量
    latitude_delta = latitude_distance * latitude_per_meter
    new_latitude = @latitude + latitude_delta

    # 経線上の移動距離
    longitude_distance = distance * Math.sin(heading * Math::PI / 180)

    # 1mあたりの経度
    earth_radius_at_longitude = EARTH_RADIUS * Math.cos(new_latitude * Math::PI / 180)
    earth_circle_at_longitude = 2 * Math::PI * earth_radius_at_longitude
    longitude_per_meter = 360 / earth_circle_at_longitude

    # 経度の変化量
    longitude_delta = longitude_distance * longitude_per_meter

    Position.new(new_latitude, @longitude + longitude_delta)
  end
end