Basemap のさまざまな投影法を使用すると、奇妙な動作が発生します。
世界地図にプロットしたい測定グリッドの形状は [181,83] です。つまり、経度 -180° ～ 180° および緯度 -82° ～ 82° の範囲の 2°/2° ポイントごとに値があります。
from mpl_toolkits.basemap import Basemap
import numpy as np
measurementgrid = np.random.random_sample((181,83))
m = Basemap(projection='cyl',llcrnrlon=-180, llcrnrlat=-82, urcrnrlon=180, urcrnrlat=82, resolution='l')
m.drawcoastlines()
m.drawparallels(np.arange(-90,90,30),labels=[1,0,0,0])
m.drawmeridians(np.arange(-180,180,45), labels=[0,0,0,1])
data, x, y = m.transform_scalar(measurementgrid.T, lons=np.arange(-180,182,2), lats=np.arange(-82,84,2), nx = 181, ny = 83, returnxy=True, order=0)
m.imshow(data, origin='lower', interpolation='none')
円筒投影法を使用すると、返されたデータ グリッドは測定グリッドと等しくなり、すべて問題ありません。投影を「ミル」に変更すると、補間されたデータは元のデータとは異なります。
測定グリッドをそのままプロットする方法はありますが、投影の変化に関してですか?