0

ポイントのコレクションの平均と2番目のポイントの間に大円をプロットしようとしています。ただし、平均を使用すると大円経路は正しくありませんが、平均と同じ事前定義された変数を使用すると正しいです。以下に例を示します。

from mpl_toolkits.basemap import Basemap
import pyproj
from pylab import *

# setup basemap.
lat_1 = 49.5
lat_2 = 55.5
lat_0 = 52.5
lon_1 = -16.5
lon_2 = -1.0
lon_0 = -5.25


m = Basemap(llcrnrlon=lon_1,llcrnrlat=lat_1,urcrnrlon=lon_2,urcrnrlat=lat_2,
        resolution='i',projection='lcc',lon_0=lon_0,lat_0=lat_0)        
m.drawcoastlines()

lon1 = -8.2395881636363644
lat1 = 54.876441136363631

lons = [-8.3771698, -8.177818, -8.2758546, -8.2197834, -8.3409037, -8.247925, -8.3208079, -8.3288984, -8.1521259, -8.0139167, -8.1802664]
lats = [54.9846137, 54.8639587, 54.8730255, 54.89845, 54.9096928, 54.9390132, 54.7620135, 54.8300078, 54.8883947, 54.8743167, 54.8173659]

clon = mean(lons)
clat = mean(lats)

g = pyproj.Geod(ellps='WGS84')
az = 261.0
lon2, lat2, baz1 = g.fwd(lon1, lat1, az, 1000000)

m.drawgreatcircle(lon1, lat1, lon2, lat2, del_s=0.1)
m.drawgreatcircle(clon, clat, lon2, lat2, del_s=0.1)

clon==lon1
clat==lat1

誰が問題が何であるか知っていますか??

4

1 に答える 1

2 
>>> type(clon)
<type 'numpy.float64'>
>>> type(lon1)
<type 'float'>

おそらくdrawgreatcircle()どうすればいいのかわからないnumpy.float64ので、 で通常の float に変換します

>>> clon = float(clon)
>>> type(clon)
<type 'float'>
于 2013-03-12T22:08:48.440 に答える