キョリの計算を rgeo に変更しようと思ったが、精度が出せないのでやめて、さらに既存の式をヒュベニの近似式に変更