ワイヤーフレームにおける線の非表示

Question

#質問
Pythonでワイヤーフレーム表示をした際に次に示す結果グラフにて約160hz近辺で時間軸上に伸びている線を消去する方法が分からず困っています。

#現状の写真

#コード

import numpy as np
import scipy.io
import matplotlib.pyplot as plt
import makegraph as mg
import analysisfunction as af
from scipy import signal
from scipy import fftpack
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import mathtext
import matplotlib as mpl

def FFTparameter(Fs,N):
    fmax=int(Fs/2.56)#最大周波数
    Line=int(N/2.56)#ライン数
    dt=1/Fs#サンプリング間隔
    df=fmax/Line
    t=np.arange(0,N*dt,dt)#時間軸
    freq=np.arange(0,Line*df,df)#周波数軸
    return fmax,Line,dt,df,t,freq
def wave1_FFT(data,samplerate,df,Fs,overlap,window):
    
    Ts = len(data) / samplerate
    Fc = Fs / samplerate
    x_ol = Fs * (1 - (overlap/100))
    N_ave = int((Ts - (Fc * (overlap/100))) / (Fc * (1-(overlap/100)))) 
    xdata=np.linspace(0,samplerate,N)#周波数データ
    ydata=np.linspace(0,len(data)/Fs,N_ave)#時間データ
    zdata=[]#振幅データ
    for j in range(N_ave):
        array = []
        ps = int(x_ol * j)
        array.append(data[ps:ps+N:1])
        array_window=array*window
        F = fftpack.fft(array_window)
        Famp=np.abs(F)
        F_total=np.sum(Famp,axis=0)
        F_average=F_total /  Fs*2
        acf=1/(sum(window)/Fs)
        F_average=acf*F_average
        zdata.append(F_average)
    return xdata,ydata,zdata

def wave_plot(xmin,xmax,ymin,ymax,zmin,zmax,ver,hor,savename):
    fig, ax = plt.subplots(figsize=(10, 10), subplot_kw={'projection' : '3d'})
    mpl.rcParams["font.size"] = 20
    mpl.rcParams['axes.labelpad'] = 20
    ax.set_xlim(xmin,xmax)
    X[X>xmax]=np.nan
    ax.set_xticks(np.linspace(xmin, xmax, 11))
    ax.set_ylim(ymin,ymax)
    ax.set_yticks(np.linspace(ymin, ymax, 9))
    ax.set_zticks(np.linspace(zmin,zmax,3))
    ax.set_xlabel("Frequency [Hz]")
    ax.set_ylabel("Time [s]")
    ax.set_zlabel("Mag [m/s$^{^2},_{0-p}$]")
    ax.plot_wireframe(X,Y,Z,color='black')
    ax.view_init(elev=ver, azim=hor)
    plt.savefig(savename)
    
Fs=8192
N=8192
overlap=25

#flattop_window = signal.flattop(N)  # フラットトップ窓
hanning_window = np.hanning(N)  # ハン窓
#rec_window=signal.boxcar(N)#レクタンギュラ窓

matdata = scipy.io.loadmat('SweepSignalFs8192.mat')
wavedata=np.array(matdata['Data']).squeeze()
fmax,Line,dt,df,t,freq=FFTparameter(Fs,N)

x1,y1,z1=wave1_FFT(wavedata,N,df,Fs,overlap,hanning_window)
X,Y=np.meshgrid(x1,y1)
Z=np.array(z1)

wave_plot(0.0,200,0.0,80,0.0,0.10,60,270,'wavetest.png')

#考察
本来は最大周波数が8192hzまであるグラフを一部拡大したもので、160Hz近傍の縦線はワイヤーフレームの一部と考えられるので、ワイヤーフレーム表示以外の表示方法があればそれを試すべきかと考えています。
そのため、今回はax.plot_wireframeの部分をax.plot_lineなどいろいろ試してみたがエラーが発生しました。（どうやらこのようなコマンドは定義されていないみたいでした）
以下にグラフ作成データを載せますので、誰か解決できる方がいらっしゃいましたら宜しくお願い致します。
#グラフ用データ
全て公開するのは量的に非現実的なので、一部だけ抜粋します。各値に書いてあるデータ数は本来のデータの個数になります。

X値　データ数153個
0
1.30719
2.61438
3.92157
5.22876
6.53595
7.84314
9.15033
10.4575
11.7647
13.0719

Y値　データ数8192個
0
0.00012207
0.000244141
0.000366211
0.000488281
0.000610352
0.000732422
0.000854492
0.000976562
0.00109863

Z値　153*8192の2重リスト
0.000101528	4.43029e-05	1.03922e-05	1.70963e-05	2.80585e-05	1.65718e-05
0.000244378	0.00031357	0.000348129	0.00024573	7.65761e-05	3.39492e-05
0.000317529	0.000335056	0.000186748	0.000205272	5.63488e-05	2.78178e-05
0.000240313	0.000222441	0.000156231	0.000105827	6.57634e-05	8.20557e-05
0.000393149	0.000314569	0.000144609	0.000138492	7.77692e-05	1.81025e-05
9.54087e-05	0.000199687	0.000153193	3.83784e-05	7.28451e-05	7.07395e-05
0.000202067	0.000247239	0.000311838	0.000290076	0.000184408	6.8552e-05
0.000312637	0.000431914	0.000233941	0.000116044	2.35009e-05	3.24477e-05
8.87849e-05	0.000143231	0.000240549	0.000200959	0.000100509	4.50792e-05

Comments

[etherbeg]

import numpy as np
import scipy.io
import matplotlib.pyplot as plt
from scipy import signal
from scipy import fftpack
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import mathtext
import matplotlib as mpl

X = np.array([
0,
1.30719,
2.61438,
3.92157,
5.22876,
6.53595,
7.84314,
9.15033,
10.4575,
11.7647,
13.0719,
])

Y = np.array([
0,
0.00012207,
0.000244141,
0.000366211,
0.000488281,
0.000610352,
0.000732422,
0.000854492,
0.000976562,
0.00109863,
])

Z = np.array([
[0.000101528, 4.43029e-05, 1.03922e-05, 1.70963e-05, 2.80585e-05, 1.65718e-05,],
[0.000244378, 0.00031357, 0.000348129, 0.00024573, 7.65761e-05, 3.39492e-05,],
[0.000317529, 0.000335056, 0.000186748, 0.000205272, 5.63488e-05, 2.78178e-05,],
[0.000240313, 0.000222441, 0.000156231, 0.000105827, 6.57634e-05, 8.20557e-05,],
[0.000393149, 0.000314569, 0.000144609, 0.000138492, 7.77692e-05, 1.81025e-05,],
[9.54087e-05, 0.000199687, 0.000153193, 3.83784e-05, 7.28451e-05, 7.07395e-05,],
[0.000202067, 0.000247239, 0.000311838, 0.000290076, 0.000184408, 6.8552e-05,],
[0.000312637, 0.000431914, 0.000233941, 0.000116044, 2.35009e-05, 3.24477e-05,],
[8.87849e-05, 0.000143231, 0.000240549, 0.000200959, 0.000100509, 4.50792e-05,],
])

def wave_plot(xmin,xmax,ymin,ymax,zmin,zmax,ver,hor,savename):
fig, ax = plt.subplots(figsize=(10, 10), subplot_kw={'projection' : '3d'})
mpl.rcParams["font.size"] = 20
mpl.rcParams['axes.labelpad'] = 20
ax.set_xlim(xmin,xmax)
ax.set_xticks(np.linspace(xmin, xmax, 11))
ax.set_ylim(ymin,ymax)
ax.set_yticks(np.linspace(ymin, ymax, 9))
ax.set_zticks(np.linspace(zmin,zmax,3))
ax.set_xlabel("Frequency [Hz]")
ax.set_ylabel("Time [s]")
ax.set_zlabel("Mag [m/s$^{^2},_{0-p}$]")
ax.plot_wireframe(X,Y,Z,color='black')
ax.view_init(elev=ver, azim=hor)
plt.savefig(savename)

wave_plot(0.0,200,0.0,80,0.0,0.10,60,270,'wavetest.png')

上記のプログラムで試したところ、以下のエラーが出たのですが、どこが間違っているでしょうか。

Traceback (most recent call last):
File "/Users/xxxxx/main.py", line 120, in <module>
wave_plot(0.0,200,0.0,80,0.0,0.10,60,270,'wavetest.png')
File "/Users/xxxxx/main.py", line 98, in wave_plot
ax.plot_wireframe(X,Y,Z,color='black')
File "/Users/xxxxx/.venv/lib/python3.9/site-packages/matplotlib/_api/deprecation.py", line 431, in wrapper
return func(*inner_args, **inner_kwargs)
File "/Users/xxxxx/.venv/lib/python3.9/site-packages/mpl_toolkits/mplot3d/axes3d.py", line 1898, in plot_wireframe
X, Y, Z = np.broadcast_arrays(X, Y, Z)
File "<__array_function__ internals>", line 5, in broadcast_arrays
File "/Users/xxxxx/.venv/lib/python3.9/site-packages/numpy/lib/stride_tricks.py", line 538, in broadcast_arrays
shape = _broadcast_shape(*args)
File "/Users/xxxxx/.venv/lib/python3.9/site-packages/numpy/lib/stride_tricks.py", line 420, in _broadcast_shape
b = np.broadcast(*args[:32])
ValueError: shape mismatch: objects cannot be broadcast to a single shape

Answer 1

質問してから自分でも調べた結果、

ax.plot_wireframe(X,Y,Z,color='black')
                  ↓
ax.plot_wireframe(X,Y,Z,color='black',ccount=1)

と変更することで上手くいきました。
実際には縦線が消えたというより、ccountの値を調節することで縦線の本数を1本にしているようです。
（１本の縦線は0Hzに存在する模様）
ちなみに横線の本数を設定するのはrcountと使うと良いらしいです。