3

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)すると、データは表示されますが、範囲が間違っています:

ここに画像の説明を入力

助言がありますか?

4

1 に答える 1