import os
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
from mpl_toolkits.basemap import Basemap
#read one pass
passid=1
cfodata=np.loadtxt('pass'+str(passid)+'.txt',skiprows=1)
#数据读取
cfolat=cfodata[:,0]
cfolon=cfodata[:,1]
cfotime=cfodata[:,2]
#read multi passes
for passid in range(2,10):
cfodata=np.loadtxt('pass'+str(passid)+'.txt',skiprows=1)
cfolat=np.concatenate((cfolat,cfodata[:,0]))
cfolon=np.concatenate((cfolon,cfodata[:,1]))
cfotime=np.concatenate((cfotime,cfodata[:,2]))
#draw
m=Basemap(llcrnrlon=-180,llcrnrlat=-80,\
urcrnrlon=180,urcrnrlat=80,projection='mill',resolution='l')
m.drawcoastlines(linewidth=0.3)
m.drawparallels(np.arange(-90.,90.,30.),labels=[1,0,1,0],fontsize=15)
m.drawmeridians(np.arange(0.,360.,45.),labels=[0,0,0,1],fontsize=15)
mlon,mlat=m(cfolon,cfolat)
plt.scatter(mlon,mlat,c=cfotime)
cb=plt.colorbar()
plt.title('Time(s) since 2009-1-1')
plt.show()