Cartopy で Aster GDEM データを使用しようとしていますが、うまくいきません。たとえば、地理的な範囲がわかっている numpy 配列があるとします。cartopy.io.img_tiles の GoogleTile などを使用してレイヤーとしてプロットできるはずです。
サンプルコード:
proj = ccrs.Mercator()
extent = [32, 33, 29, 30]
data = np.random.rand(10,10)
fig = plt.figure(figsize=(4,4))
ax = fig.add_subplot(111, projection=proj)
ax.set_extent(extent)
ax.imshow(data, extent=extent, origin='upper',
transform=proj, alpha=0.5)
gg_tiles = GoogleTiles()
ax.add_image(gg_tiles, 10, alpha=0.5)
ax.coastlines('10m')
ax.set_title('data with Google Tile')
gl = ax.gridlines(draw_labels=True)
gl.xlabels_top = None
生成:
データが表示されません!
コメントax.set_extent(extent)
すると、データは表示されますが、範囲が間違っています:
助言がありますか?