1
cartopy、matplotlib、およびimshowを使用して、等間隔(緯度/経度)データの正方形グリッドをプロットしようとしています。データはデータラインを横切っており、正しく動作するようにマップを取得することに問題がありました。imshowとcartopyのデータ交差点を含む画像で
ここに私の問題の例です:
import numpy as np
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
lat = np.arange(6000)*0.02 + (-59.99)
lon = np.arange(6000)*0.02 + (85.01)
dat = np.reshape(np.arange(6000*6000),[6000,6000])
tran = ccrs.PlateCarree()
proj = tran
plt.figure(figsize=(8,8))
ax = plt.axes(projection=proj)
print([lon[0],lon[-1],lat[0],lat[-1]])
ax.imshow(dat, extent=[lon[0],lon[-1],lat[0],lat[-1]],transform=tran,interpolation='nearest')
ax.coastlines(resolution='50m', color='black', linewidth=2)
ax.gridlines(crs=proj,draw_labels=True)
plt.show()
tran = ccrs.PlateCarree(central_longitude=180)
proj = tran
plt.figure(figsize=(8,8))
ax = plt.axes(projection=proj)
print([lon[0]-180,lon[-1]-180,lat[0],lat[-1]])
ax.imshow(dat, extent=[lon[0]-180,lon[-1]-180,lat[0],lat[-1]],transform=tran,interpolation='nearest')
ax.coastlines(resolution='50m', color='black', linewidth=2)
ax.gridlines(crs=tran,draw_labels=True)
plt.show()
は、最初のプロットは180Eで降りチョッピング、この画像が得られます。
二フィックスマップの問題を、しかし、グリッドのダニは今間違っています:
私は再投射を試みましたが、私は(tran!= proj)と思っていますが、見かけ上ハングアップしすぎていました。
私は基本的には下の画像が必要ですが、適切なラベルが付いています。私は、より多くの地理的に離れたデータを重ねて表示するつもりだから、今はちょうどハックのように見えるのではなく、正確にそれをやりたいのです。
はcartopyを持っていない、しかし、あなたはちょうどあなたが 'ax.set_xticklabels'で欲しいものにラベルを変更することはできませんか? –
私の答えに関するコメント/質問はありますか? – swatchai