このような絵をPythonでプロットしようとしています。
プロット用の 3 つのパラメーターがあります。バツ:
[ 0.03570416 0.05201517 0.05418171 0.01868341 0.07116423 0.07547471]
y:
[-0.32079484 -0.53330218 -1.02866859 -0.94808545 -0.51682506 -0.26788337]
z:
[-0.32079484 -0.53330218 -1.02866859 -0.94808545 -0.51682506 -0.26788337]
したがって、x は x 軸、y は y 軸です。ただし、z はピクセルの強度です。私はこのコードを思いつきます:
z = np.array(reals)
x = np.array(ra)
y = np.array(dec)
nrows, ncols = 10, 10
grid = z.reshape((nrows, ncols))
plt.imshow(grid, extent=(x.min(), x.max(), y.max(), y.min()), interpolation='nearest', cmap=cm.gist_rainbow)
plt.title('This is a phase function')
plt.xlabel('ra')
plt.ylabel('dec')
plt.show()
ただし、次のエラーが表示されます。
grid = z.reshape((nrows, ncols))
ValueError: total size of new array must be unchanged
ra、dec、および reals は、同じサイズの通常の配列です。以前にそれらを計算してから、それらを使用してnumpy配列を作成しました