Program

Pythonパッケージpyprojを使って２つの座標の距離を求める

最近はPython使って位置情報をごりごり加工してますコンバンハ。

以前はMongoDBを使って半径XXＫｍ以内にある建物を求めるってことをやりました。今回は２つの座標間の距離を求めてみます。すっごく簡単なのでオススメ。

パッケージはpyprojを使います。

pyprojインストール
$ pip install pyproj

コード

[crayon-5d361dd956c12176447196 lang="python"  title="getdistance.py"]
#!/usr/bin/env python
# -*- coding:utf-8 -*-
from pyproj import Geod

class GetDistance:
    def get_distance(self, start, to):
        """
        pyprojを使って距離を求める
        """
        res = u'[ %(from)s ]から[ %(to)s ]まで %(distance_str)s (%(distance)s)'

        q = Geod(ellps='WGS84')
        fa, ba, d = q.inv(start['lon'], start['lat'], to['lon'], to['lat'])

        print res % {
            'from'  : start['name'],
            'to'    : to['name'],
            'distance_str'  : self._cutdown(d),
            'distance'      : d,
        }


    def _cutdown(self, num):
        """
        距離に単位をつけて返す

        1000m以下の時はメートル表記、
        それ以上の時はキロメートル表記で返す。
        """
        val = int(round(num))
        if val < 1000:
            return '%sm' % val
        else:
            km = val * 0.001
            return '%sKm' % round(km, 1)





if __name__ == '__main__':
    # 東京タワー
    tokyo_tw = {
        'name': u'東京タワー',
        'lat' : 35.65861,
        'lon' : 139.745447,
    }

    # 芝公園
    shiba_park = {
        'name': u'芝公園',
        'lat' : 35.654071,
        'lon' : 139.749838,
    }

    # スカイツリー
    skytree = {
        'name': u'スカイツリー',
        'lat' : 35.71007,
        'lon' : 139.80948,
    }

    # 富士山
    mt_fuji = {
        'name': u'フジサーン',
        'lat' : 35.360556,
        'lon' : 138.727778,
    }

    g = GetDistance()
    g.get_distance(tokyo_tw, shiba_park)
    g.get_distance(tokyo_tw, skytree)
    g.get_distance(skytree, mt_fuji)
    g.get_distance(tokyo_tw, mt_fuji)

[crayon-5d361dd956c12176447196 lang="python" title="getdistance.py"]

#!/usr/bin/env python

# -*- coding:utf-8 -*-

from pyproj import Geod

class GetDistance:

def get_distance(self, start, to):

"""

pyprojを使って距離を求める

"""

res = u'[ %(from)s ]から[ %(to)s ]まで %(distance_str)s (%(distance)s)'

q = Geod(ellps='WGS84')

fa, ba, d = q.inv(start['lon'], start['lat'], to['lon'], to['lat'])

print res % {

'from' : start['name'],

'to' : to['name'],

'distance_str' : self._cutdown(d),

'distance' : d,

}

def _cutdown(self, num):

"""

距離に単位をつけて返す

1000m以下の時はメートル表記、

それ以上の時はキロメートル表記で返す。

"""

val = int(round(num))

if val < 1000:

return '%sm' % val

else:

km = val * 0.001

return '%sKm' % round(km, 1)

if __name__ == '__main__':

# 東京タワー

tokyo_tw = {

'name': u'東京タワー',

'lat' : 35.65861,

'lon' : 139.745447,

}

# 芝公園

shiba_park = {

'name': u'芝公園',

'lat' : 35.654071,

'lon' : 139.749838,

}

# スカイツリー

skytree = {

'name': u'スカイツリー',

'lat' : 35.71007,

'lon' : 139.80948,

}

# 富士山

mt_fuji = {

'name': u'フジサーン',

'lat' : 35.360556,

'lon' : 138.727778,

}

g = GetDistance()

g.get_distance(tokyo_tw, shiba_park)

g.get_distance(tokyo_tw, skytree)

g.get_distance(skytree, mt_fuji)

g.get_distance(tokyo_tw, mt_fuji)

[/crayon]

実行結果：

[ 東京タワー ]から[ 芝公園 ]まで 642m (641.660523655)
[ 東京タワー ]から[ スカイツリー ]まで 8.1Km (8136.21631144)
[ スカイツリー ]から[ フジサーン ]まで 105.5Km (105484.921113)
[ 東京タワー ]から[ フジサーン ]まで 98.1Km (98064.787889)