>>> import numpy as np

（4）定义一个指标数组

>>> import numpy as np

>>> a=np.arange(1,11,1)

>>> b=np.arange(1,11,1).reshape(2,5)

3. 计算
（1）傅里叶变换

>>> import numpy as np

>>> a=np.ones(100)

>>> b=np.fft.fftn(a)

>>> c=np.fft.ifftn(b).real

4. 画图
（1）积分强度图

# in intensity_map.py:

import pyfits
import numpy as np
import matplotlib.pyplot as plt

hdulist = pyfits.open('xxx.fits')
image = hdulist[0].data
nz = hdulist[0].header['naxis3']
for i in range(0,nz):
 Z = Z+image[i,:,:]

ax = plt.subplot(111)
im = plt.imshow(Z, cmap=cm.gist_yarg,
               origin='lower', aspect='equal',
               extent=[max(x),min(x),min(y),max(y)])
plt.show()

（2）简单曲线图

# in simple_plot.py:

from pylab import *

t=arange(0.0, 2.0, 0.01)

s=sin(2*pi*t)

plot(t,s)

xlabel('time (s)')

ylabel('voltage (mV)')

savefig("test.eps")

show()

5. 查看帮助