筛选并合并星表

#coding=utf-8
import pandas as pd
from astropy.table import Table, vstack,hstack

def filter_gaia(oldgaia):

    data = Table.read(oldgaia)
    ra = data['ra']
    dec = data['dec']
    # Gmag = data['phot_g_mean_mag']
    Gmag = data['Gmag']
    id = data['id']
    ndata =hstack([ra,dec,id,Gmag])
    print(ndata)
    ls_data.append(ndata)
def combine_gaia(newgaia):
    comdata = vstack(ls_data)
    print(comdata)
    df = pd.DataFrame(
        {'id': comdata['id'],'ra': comdata['ra'], 'dec': comdata['dec'], 'Gmag':comdata['Gmag']})
    df.to_csv(newgaia, index=False)


alist = [line.rstrip() for line in open('gaia.dat')]
step_num = 19

clist = []
for i in range(0, len(alist), step_num):
    print('processing', i)
    blist = []
    for j in range(i, i + step_num):
        if j > len(alist) - 1:
            break
        else:
            print(alist[j])
        blist.append(alist[j])
    clist.append(blist)



for k in clist:
    print('processing', k)
    ls_data = []
    for l in k:
        oldgaia = l
        filter_gaia(oldgaia)
    newgaia = 'new' + k[0][0:-5] + '.csv'
    combine_gaia(newgaia)