地図に画像を張る¶

地図上に衛星画像を貼ってみよう。ここではNASA EarthData Galleryから 2017年台風第5号Noruの画像を取り込む。Windowsでは，保存する際jpgの前の.以外は_に変換され，Noru_A2017217_0435_2km.jpgのようになることがある。その際は以下のスクリプト中でも_を使ったファイル名を用いること。

In [1]:

import matplotlib.pyplot as plt
import cartopy.crs as ccrs

fname = "Noru.A2017217.0435.2km.jpg"
img = plt.imread(fname)

imread()でJPEG画像を読み，imgに格納した。imgには何が入っているのだろうか。まずtype()型を調べてみる。

In [2]:

type(img)

Out[2]:

numpy.ndarray

Numpyの数値配列ndarrayであることが分かった。形状をshapeで調べる。

In [3]:

img.shape

Out[3]:

(1050, 800, 3)

1050x800x3の3次元配列である。大きさが1050x800画素でRGBのカラー画像だ。最初の画素は何色だろうか。

In [4]:

img[0,0,:]

Out[4]:

array([248, 247, 245], dtype=uint8)

imshow()で表示する。

In [5]:

plt.imshow(img)
plt.show()

画像のページの「more information」を見ると四隅の座標が分かる。img_extentというタプルをimshow()のキーワード引数extentに与えて地図上に貼る。

In [6]:

fig = plt.figure(figsize=(8, 12))

img_extent = (121.7937, 138.4089, 18.5564, 37.4421)
ax = plt.axes(projection=ccrs.PlateCarree())
plt.title('Aqua/MODIS Typhoon Noru\n' '2017-08-05 04:35 UTC')

ax.imshow(img, origin='upper', extent=img_extent, transform=ccrs.PlateCarree())
ax.coastlines(resolution='10m', color='black', linewidth=1)

ax.set_extent([120, 140, 15, 40])

ax.plot(130.557139, 31.596806, 'bo', markersize=7, transform=ccrs.Geodetic())
ax.text(131, 31.6, '鹿児島', transform=ccrs.Geodetic())

plt.show()

In [ ]: