0

3D グラフにカラーマップを使用して論文の 4 次元データをプロットしていて、問題が発生しました。私が使用しているカラーマップ メソッドは、コーナー ポイントの値を平均し、その値でタイル全体を色付けしているように見えます。解像度が高ければこれで問題ありませんが、現在のデータを取得するためにシミュレーションを実行するのに約 1 か月かかりました。

誰かがこれを変更する方法を提案できますか? できれば、解像度を上げるためにすべてのポイントの独自の線形補間をコーディングしないでください。それはおそらく、この時点での価値よりも多くの努力です (私にとって)。

Driven=np.zeros((5,9))
Driver=np.zeros((5,9))
Compositions=np.zeros((5,9))
Durations=np.zeros((5,9))
N=np.zeros((5,9))
for i in range(0,5):
    for j in range(0,9):
        Driven[i,j]=ReservoirData[i][1][j]
        Compositions[i,j]=ReservoirData[i][2][j]
        Driver[i,j]=float(ReservoirData[i][0][0][:-3])
        Durations[i,j]=ReservoirData[i][3][j]
maxi=Durations.max()
mini=Durations.min()
for i in range(0,5):
    for j in range(0,9):
        N[i,j]=(Durations[i,j]-mini)/(maxi-mini)

fig = plt.figure() 
ax = fig.gca(projection='3d')

#Important Stuff    Start----------------------------------------------------

surf = ax.plot_surface(Driven,Driver,Compositions, facecolors=cm.jet(N), rstride=1, cstride=1, antialiased=True)
m = cm.ScalarMappable(cmap=cm.jet)

#Important Stuff End---------------------------------------------------------

m.set_array(Durations)
cbar=plt.colorbar(m, shrink=0.8)
cbar.set_label('Nominal Duration')
ax.set_ylabel('Driver Pressure, kPa')
ax.set_xlabel('Compositions, %He')
ax.set_zlabel('Driven Pressure, kPa')
plt.title('Three Dimensional Representation of Tailored Conditions for RS at 9.2MPa')
fig.set_size_inches(14,8)
plt.savefig('RS9.2.png')

ここに画像の説明を入力 どんな意見でも大歓迎です、ありがとう!

4

1 に答える 1