1

言葉で説明しにくく申し訳ないです。
以下の図のaxbar3dのxy平面上の水色の四角い図形の影を取り除きたいのですが、どのような対処を施せばよいのでしょうか?公式サイトにも自分が調べる限り載っておらず、stackoverflowの海外の方の似たような質問があったのですが、少し違く悩んでおります。
以下に見本のコードを掲載します。

import numpy as np
import matplotlib.pyplot as plt
import math
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.animation as animation
import matplotlib.colors as colors
from matplotlib import cm


n=10 
m=10 
theta = 3*(math.pi)/12
p_spot=[]
P_spot=[]

P = [[np.cos(theta),np.sin(theta)],[0,0]]
Q = [[0,0],[np.sin(theta),-np.cos(theta)]]
x_list=[]
t_list=[]
p_list=[]
s_list=[]
a = 1/math.sqrt(2)
b = 1j/math.sqrt(2)
p_map=[]
pp_map =np.zeros([2*m+1,2*m+1])#,dtype="complex")
R=1/2
L=1/2
X_list=[]
P_list=[]

for j in range(0,2*n+1):
    if j  == n:
        phai = [a ,b]
        pro = 1
    else:
        phai = [0,0]
        pro =0
    p = np.dot(phai,np.conj(phai))

    x_list.append(j)
    X_list.append(j)
    s_list.append(phai)
    p_list.append(p)
    P_list.append(pro)



for t in range(0,2*m+1):
    t_list.append(t)
    if t ==0:
        s_list
        p_list
        P_list
    else:
        next_s_list = [0]*len(s_list)
        next_P_list = [0]*len(P_list) #listと同じ要素の数ですべて0を用意(初期化)
        for i in range(0,2*n+1):
            if i == 0:
                next_s_list[i] = np.dot(P, s_list[i+1])
                next_P_list[i] = P_list[i+1]*L
            elif i == 2*n:
                next_s_list[i] = np.dot(Q, s_list[i-1])
                next_P_list[i] = P_list[i-1]*R
            else:
                next_s_list[i] = np.dot(P, s_list[i+1]) + np.dot(Q, s_list[i-1])
                next_P_list[i] = P_list[i+1]*L + P_list[i-1]*R

            p_list[i] = np.dot(next_s_list[i],np.conj(next_s_list[i]))
            #pp_map[t]=p_list


        s_list = next_s_list
        P_list = next_P_list
    pp_map[t]=np.real(p_list)
    print(t,np.real(pp_map),np.real(p_list))


fig= plt.figure()
ax = Axes3D(fig, rect=(0.1,0.1,0.8,0.8)) #rect=(x0,y0,width,height)
X,Y = np.meshgrid(x_list, t_list)
ax.set_xlabel("Position",labelpad=10,fontsize=24)
ax.set_ylabel("Time",labelpad=20,fontsize=24)
ax.set_zlabel("|φ|^2",labelpad=10,fontsize=18)
ax.set_xlim(2*n,0)
ax.set_ylim(2*n,0)
ax.set_zlim(0,1)
offset = pp_map.ravel() + np.abs(pp_map.min())
fracs = offset.astype(float)/offset.max()
norm = colors.Normalize(fracs.min(), fracs.max())
clrs = cm.cool(norm(fracs))
ax.bar3d(X.ravel(), Y.ravel(), pp_map.ravel() ,0.5, 0.5, -pp_map.ravel(),color =clrs)
ax.w_xaxis.set_pane_color((0, 0, 0, 0))
ax.w_yaxis.set_pane_color((0, 0, 0, 0))
ax.w_zaxis.set_pane_color((0, 0, 0, 1))
ax.grid(color="white")
ax.grid(False)

plt.show()

画像の説明をここに入力

3
  • これが参考になるかも。 pythonprogramming.net/3d-bar-chart-matplotlib-tutorial
    – user28902
    Commented 2018年6月30日 12:29
  • @Rei さん、回答者は Rei さんがお調べになったページがどれなのか判別できないので、手間を減らすためにも、少し違ったページへのリンクを追記して頂けると嬉しいです。
    – nekketsuuu
    Commented 2018年6月30日 13:09
  • 申し訳ないです、以後気をつけます!
    – 0018
    Commented 2018年7月1日 13:18

1 件の回答 1

0

原因

平面に見えている格子状の正方形たちは、データに 0.0 が混じっていることで生じています。影ではなく、0.0 を示すための高さ 0 の棒です。

方法1

0.0 の部分を欠損値として取り除きたいということであれば、欠損値部分を NaN で埋めるという方法があります。NaN 部分の棒は描画されません。

pp_map[pp_map <= 0.0] = np.nan

また、これに合わせて min, maxnp.nanmin, np.nanmax に変える必要もあるでしょう。

方法2

あるいはそもそも高さ 0 の部分を最初からデータより取り除く方法もあります。方法1と同じようにマスクを作ってそれぞれに付ければ良いです。具体的には

mask = pp_map > 0.0

として、X[mask], Y[mask], pp_map[mask] とすれば良いです。色を計算する部分でも offset の計算で pp_map を制限する必要があります。

2
  • @Rei さん、スタック・オーバーフローでは、回答の左にあるチェックマークを押すことで回答を承認することができます。必須ではありませんが、よろしければ押して頂けると嬉しいです :) 過去のご質問も確認してみてください。
    – nekketsuuu
    Commented 2018年7月1日 13:19
  • あ!そうなんですね、チェックしました!!至らないところも多々ありますが、今後ともよろしくお願い致します!
    – 0018
    Commented 2018年7月2日 1:50

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。