imshow и pcolormesh показывают два разных вывода с одинаковым вводом - PullRequest
0 голосов
/ 28 мая 2019

Когда я пытаюсь понять разницу между imshow и pcolormesh, я сталкиваюсь с такой проблемой, что не могу понять, почему.

Например, у меня синусоида с частотой 10 кГц

import numpy as np
import matplotlib.pyplot as plt 
from scipy import fft

f = 10000  # Frequency, in cycles per second, or Hertz
fs = 44100  # Sampling rate, or number of measurements per second

t = np.linspace(0, 5, 5 * fs, endpoint=False)
X = np.sin(f * 2 * np.pi * t)

Применение fft к этому сигналу дает спектрограмму Xs

window_size = 2048  # 2048-sample fourier windows
stride = 512        # 512 samples between windows
wps = fs/float(512) # ~86 windows/second
Xs = np.empty([int(2*wps),2048])

for i in range(Xs.shape[0]):
    Xs[i] = np.abs(fft(X[i*stride:i*stride+window_size]))

Теперь я хотел бы построить спектрограмму Xs с imshow и pcolormesh.

Imshow

fig = plt.figure(figsize=(20,7))
plt.imshow(Xs.T,aspect='auto',cmap='seismic')
plt.gca().invert_yaxis()
plt.colorbar()

enter image description here

pcolormesh

fig = plt.figure(figsize=(20,7))
plt.pcolormesh(Xs.T,cmap='seismic')
plt.colorbar()

enter image description here

Чем они отличаются?

1 Ответ

0 голосов
/ 31 мая 2019

Интересно, что добавление аргумента dpi решает проблемы. Даже установив dpi=100, который является значением по умолчанию, если аргумент dpi не задан, проблема все равно будет решена.

Imshow

fig = plt.figure(figsize=(20,7), dpi=100)
plt.imshow(Xs.T,aspect='auto',cmap='seismic')
plt.gca().invert_yaxis()
plt.colorbar()

enter image description here

Хотя аргумент dpi решает проблему, мне все еще неясно, почему это происходит. Любые дополнительные предложения приветствуются.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...