读FITS文件头的参数值和注释

import pyfits

hdulist=pyfits.open(fitsfilename)

hdu0=hdulist[0]

header0=hdu0.header

headercomments0=hdu0.header.comments

print header0['BITPIX'],headercomments0['BITPIX']

删除字符串中的空格

"   xyz   ".strip()# returns "xyz"  

"   xyz   ".lstrip()# returns "xyz   "  

"   xyz   ".rstrip()# returns "   xyz"  

"  x y z  ".replace(' ','')# returns "xyz"

读文件

for line in open(path):
       line=line.replace('n','').split(" ")
       print line[0]

曲线拟合

import lmfit

获取画图范围的参数

ax.get_xlim()

列表（list）拼接和numpy数组拼接

---------------------------

================

np.append(a,b)

---------------------------

np.concatenate((a,b),axis=0)

查看数组形状

import numpy as np

print np.shape(array)

创建三维图形

import matplotlib.pyplot as plt

from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure()

ax = Axes3D(fig)

画三维散点图

ax.scatter(xs, ys, zs)

画二维散点的密度图（二维直方图）

from matplotlib.colors import LogNorm  

from pylab import *

hist2d(x, y, bins=40, norm=LogNorm())

画图的颜色