我想显示图像的幅度谱。 我能够做到这一点使用下面的代码:
import numpy as np
import matplotlib.pyplot as plt
import pylab
pylab.gray()
pic = pylab.imread("C:/pic.png")[::-1,:]
amp_pic = pylab.subplot(1,4,1)
amp_pic.xaxis.set_ticks_position('top')
pylab.imshow(np.abs(np.fft.fftshift(np.fft.fft2(pic))),\
interpolation='nearest')
pylab.show()
但轴没有标记的幅度谱应该被标记的方式。 在1D-功能的情况下,重新标记是比较容易:
import math
import numpy as np
import matplotlib.pyplot as plt
array = np.arange(149)
frequency_scale = np.fft.fftfreq(array.shape[0])
function = np.cos(2*math.pi*array/10)
fft = np.fft.fft(function)
amp_fft = np.abs(fft)
plt.subplot(1,4,1)
plt.plot(fkt,'r')
plt.show()
我要为我的imshow情节的情况下,X轴相同的标签。 这可能吗?