Pythonで世界地図-6

グラフで座標を取得する方法がわかったので、
#####座標から経度・緯度を求める
方法を探しました。

実行例
./map1.py 東京
（全ての都市で使えるわけではなさそうです。ではなくて時間が掛かって失敗しているのか？何回かやれば表示されました。ローマ字だとすぐに表示されるようです。）

地図上をクリックすると、緯度・経度が表示されます。

間違っていたので修正しました。

map1.py

#!/usr/bin/python3
# coding: UTF-8

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
from pygeocoder import Geocoder
import sys

font = {'family':'IPAGothic'} #日本語Fontを指定

args = sys.argv
argc = len(args) # 引数の個数

if (argc != 2):   # 引数がない場合
    print ('./map1.py 目的地の名称')
    quit()

area_name = args[1]

try:
    result = Geocoder.geocode(area_name)
    lon1 = float(result.longitude)
    lat1 = float(result.latitude)
except:
    quit()

map = Basemap(llcrnrlon=lon1-2,llcrnrlat=lat1-2,urcrnrlon=lon1+2,urcrnrlat=lat1+2,resolution='h',projection='cyl')

map.shadedrelief() #標高図
map.drawcoastlines()

map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))

try:
  while True:
    a=plt.ginput(n=1)[0]  #クリックして座標を取得
    b=map(a[0], a[1], inverse=True) #座標を緯度・経度に変換
    plt.title("東経{}°　北緯{}°".format(b[0],b[1]), **font)
except:
  pass

plt.show()

今、気が付きましたが、これを指定するとグラフの右下に表示される座標が、緯度・経度に代わっています。

#####2点間距離を求める
地図上を2箇所クリックする。

変更点

from mpl_toolkits.basemap import pyproj

try:
  while True:
    aa=plt.ginput(n=1)[0]  #クリックして座標を取得
    bb=map(aa[0], aa[1], inverse=True) #座標を緯度・経度に変換
    cc=plt.ginput(n=1)[0]  
    dd=map(cc[0], cc[1], inverse=True) 
    map.drawgreatcircle(bb[0],bb[1],dd[0],dd[1],linewidth=2,color='b')  #大圏コース
    gc = pyproj.Geod(a=map.rmajor,b=map.rminor) # ポイント間の距離計算
    az1, az2, dist12 = gc.inv(bb[0],bb[1],dd[0],dd[1])
    plt.title("2点間の距離＝ {} km".format(dist12/1000), **font)
except:
  pass

#####ウインドウの中で地球を動かすことができるか？
私が調べた範囲ではできないような気がする。

plt.show()を終了させることが、できなかった。
無限ループでウインドウの右上のXボタンで終了させて、次の表示まで１０秒程度表示なしで繰り返すことくらいしか出来なかったので、アニメgifにしました。
白黒の地球で数分程度？、なんか寂しいので標高図（山の高さにより色が違う）にすると、ファイルサイズは3〜５倍程度だが、無限ループで360枚pngファイルを作成しようと思ったら、メモリ不足？作成時間はだんだん長くなるわ、エラーで止まるわ４時間程度かかった。
最後２０枚程度は１枚１枚作成したら、それほど時間がかからなかった。メモリを開放しなければいけなかったのか？。

色の付いたアニメgifは出来ましたが、サイズオーバーでアップできませんでした。