web-dev-qa-db-ja.com

Python Matplotlib Basemapは、マッププロットに小さな画像をオーバーレイします

航空機のデータを地図上にプロットしています。この75pxx 29pxの飛行機のPNG画像を、プロット上の最新のデータポイントの座標に挿入したいと思います。

airplane

私が知っていて読んだ限りでは、pyplot.imshow()がこれを達成するための最良の方法です。しかし、ステップ1でハングアップし、画像を表示することもできます。ベースマップの代わりに通常のプロットを使用すると、imshowを使用して画像を表示するのは簡単ですが、ベースマップを使用すると、画像をまったく表示できません。サンプルコードを参照してください。

画像を地図上に表示できる場合は、試行錯誤により、imshow()のextent属性とプロット座標を使用して、その位置と適切な寸法を設定できると想定しています。マップ座標x,y = m(lons,lats)から変換されます。

これがサンプルコードです（試してみるには、上の飛行機の画像をダウンロードすることをお勧めします）。

_from matplotlib import pyplot as plt
from mpl_toolkits.basemap import Basemap
import Image
from numpy import arange

lats = arange(26,29,0.5)
lons = arange(-90,-87,0.5)

m = Basemap(projection='cyl',llcrnrlon=min(lons)-2,llcrnrlat=min(lats)-2,
            urcrnrlon=max(lons)+2,urcrnrlat=max(lats)+2,resolution='i')

x,y = m(lons,lats)
u,v, = arange(0,51,10),arange(0,51,10)
barbs = m.barbs(x,y,u,v)
m.drawcoastlines(); m.drawcountries(); m.drawstates()

img = Image.open('c130j_75px.png')
im = plt.imshow(img, extent=(x[-1],x[-1]+50000,y[-1],y[-1]+50000))
plt.show()
_

これが結果の画像で、飛行機の痕跡は含まれていません。 extentを使用していくつかの異なるサイズを試しましたが、小さすぎたのではないかと思いましたが、成功しませんでした。 _zorder=10_も設定してみましたが、うまくいきませんでした。どんな助けでもいただければ幸いです。

result

更新：前者はマップのaxesインスタンスを渡すため、少なくとも_m.imshow_の代わりに_plt.imshow_を使用して画像を表示できるようになりましたが、extent引数はのディメンションに影響を与えないようです画像は、ゼロに設定した場合でも、extentディメンションをいくら小さくしても、常にプロット全体を埋めるためです。飛行機の画像を適切に拡大縮小し、最後のデータポイントの近くに配置するにはどうすればよいですか？

im = m.imshow(img, extent=(x[-1],x[-1]+5,y[-1],y[-1]+2))

result2

pythonimagematplotliboverlaymatplotlib-basemap
27
Levi Cowan

実際、これには、matplotlibの文書化されていない機能を使用する必要があります。matplotlib.offsetboxモジュール。ここに例があります： http://matplotlib.sourceforge.net/trunk-docs/examples/pylab_examples/demo_annotation_box.html

あなたの場合、あなたはこのようなことをするでしょう：

import matplotlib.pyplot as plt
import numpy as np
import Image

from mpl_toolkits.basemap import Basemap
from matplotlib.offsetbox import OffsetImage, AnnotationBbox

# Set up the basemap and plot the markers.
lats = np.arange(26, 29, 0.5)
lons = np.arange(-90, -87, 0.5)

m = Basemap(projection='cyl',
            llcrnrlon=min(lons) - 2, llcrnrlat=min(lats) - 2,
            urcrnrlon=max(lons) + 2, urcrnrlat=max(lats) + 2,
            resolution='i')

x,y = m(lons,lats)
u,v, = np.arange(0,51,10), np.arange(0,51,10)
barbs = m.barbs(x,y,u,v)

m.drawcoastlines()
m.drawcountries()
m.drawstates()

# Add the plane marker at the last point.
plane = np.array(Image.open('plane.jpg'))
im = OffsetImage(plane, zoom=1)
ab = AnnotationBbox(im, (x[-1],y[-1]), xycoords='data', frameon=False)

# Get the axes object from the basemap and add the AnnotationBbox artist
m._check_ax().add_artist(ab)

plt.show()

enter image description here

これの利点は、平面が軸座標にあり、ズームインしたときに図形のサイズに対して同じサイズのままになることです。

26
Joe Kington

ベースマップでは、最初にマップインスタンスを使用して座標を変換する場合、通常は通常のピプロットスタイルのコマンドを使用できます。この場合、次の方法でエクステントをuv座標に変換できます。

x0, y0 = m(x[-1], y[-1])
x1, y1 = m(x[-1] + 0.5, y[-1] + 0.5)

その後、次のことができるようになります。

im = plt.imshow(img, extent=(x0, x1, y0, y1))

これに対する私の完全な解決策は次のようになります：

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


lats = np.arange(26, 29, 0.5)
lons = np.arange(-90, -87, 0.5)

m = Basemap(projection='cyl', llcrnrlon=min(lons)-2, llcrnrlat=min(lats)-2,
            urcrnrlon=max(lons)+2, urcrnrlat=max(lats)+2, resolution='h')

x, y = m(lons,lats)
u, v = np.arange(0, 51, 10), np.arange(0, 51, 10)
barbs = m.barbs(x, y, u, v)

m.drawcoastlines()
m.fillcontinents()

x_size, y_size = 0.8, 0.4
x0, y0 = m(x[-1] - x_size/2., y[-1] - y_size/2.)
x1, y1 = m(x[-1] + x_size/2., y[-1] + y_size/2.)
im = plt.imshow(plt.imread('mslr86.png'), extent=(x0, x1, y0, y1))

plt.show()

のような画像を生成しますenter image description here

更新：ズームに関係なく、画像を固定サイズのままにしたい場合は、ジョーの回答を参照してください。

14
pelson