Python 3D 直方体のプロットができない

Question

python, Matplotlibを使って以下の図の様な直方体のグラフが出力ができないかとトライしております。

現状はいろいろとネットの情報を参照しながら試行錯誤しておりますが、以下のコードの出力結果にあります通り現状の出力結果は目指している形とは大分乖離している状態です。

是非ヒントやアドバイス頂けますと幸いです。
よろしくお願いいたします。

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

z = 2.4
x = np.arange(0,90,1)
y = np.arange(0,400000,500)

X, Y = np.meshgrid(x, y)

def func_disk(x, y, z):
    return (x*400000)*z

Z = func_disk(X, Y, z)

fig = plt.figure()
ax = Axes3D(fig)

ax.set_xlim(0, 90)
ax.set_ylim(0, 400000)
ax.set_zlim(0, 2000)
ax.plot_wireframe(X, Y, Z, color='Blue')
plt.show()

実行結果

Answer 1

ただ単に上の直方体を書くんであれば、

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np

fig = plt.figure()
ax = Axes3D(fig)
ax.set_xlim([0,150])
ax.set_ylim([0,800])
ax.set_zlim([0,7.0])
ax.bar3d( 0, 0,  0, 90 , 400, 2.4, color ="blue")
plt.show()

ちゃんとデータをとって書くのであれば、

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np

x_list = [i for i in range(0,91)]
y_list = [i for i in range(0,401)]
Z = np.zeros([91, 401]) + 2.4

fig = plt.figure()
ax = Axes3D(fig)#, rect=(0.1,0.1,0.8,0.8))
X,Y = np.meshgrid(x_list, y_list)
ax.set_xlim([0,150])
ax.set_ylim([0,800])
ax.set_zlim([0,7.0])
ax.bar3d( X.ravel(), Y.ravel(), 0,  1 , 1, Z.ravel(), color ="blue")
plt.show()

Comments

[minhouse10]

Quantum_walker様、ご回答頂き本当にありがとうございます！
実際にシェア頂いたコードを自分の環境でも実行してみました。とてもきれいにグラフが出て感動しました。
初の３Dグラフ挑戦でしたが、これから頑張って自分の力で書けるように精進致します。とても勉強になりました！改めてありがとうとうございました！

[Fallout_18]

こちらこそ！
追記ですが、上の一つ目のコードは、一つの値で直方体を描いており、
2つ目のコードは400*90個の高さ2.4の棒グラフの集合で直方体を描いています。
この意味はおそらくpythonをこれからやっていくと理解していくでしょう。

僕も今年の3月からpythonを始めたばかりでこのサイトには何度も助けられたので、お互い頑張りましょう！投げやりな質問をすると炎上します。一度teratailのルールは読んどいたほうが良いです。

[minhouse10]

再度アドバイス頂き、ありがとうとうございます！
はい、僕も質問が丸投げにならないように、気を付けたいと思います！
確かに2つ目のコードは実行結果を得るまでに結構2，3分かかりました。ｘとｙの値をループで回している分時間がかかっているのですね。

[Fallout_18]

というよりは、プロットして図を書いているのに時間がかかっていると思います。その解釈が正解なのかもしれませんが。

[minhouse10]

ありがとうございます！確かに仰る通りPythonとは言えこの２つのforループだけで２，３分は少しかかりすぎな感じですね。（当然私のPCのスペックもあるかとは思いますが）