最近在學(xué)習(xí)一種理論，真TMD的難啊，又沒有多少圖可以看，說到的以前的分鐘線數(shù)據(jù)，在網(wǎng)上找了很久，都不容易看到。

不過利用通達(dá)信可以幫我們做到。通達(dá)信可以下載每天的分筆數(shù)據(jù)，就利用這個(gè)文本文件去產(chǎn)生以前的分時(shí)圖！

用到的工具是著名的python，也是我最喜歡的語言。

需要說明的是，我安裝的是國信，路徑是e:\cwork\guosen ， 為了方便看 歷史分時(shí)線我 我是將產(chǎn)生的5分鐘數(shù)據(jù)文件和1分鐘數(shù)據(jù)文件拷貝到英大的目錄下看的。由于我不知道通達(dá)信的1分鐘數(shù)據(jù)文件到底是如何命名的，所以我也把它取成lc5的文件名，真正看1分鐘圖時(shí)切換到5分鐘圖下看即可。當(dāng)然可能與真正的1分鐘圖、5分鐘圖有些許差異，但這并無大礙，況且同花順的1分鐘圖就和通達(dá)信的也存在細(xì)微的差異！

以下是CODE(注意：PYTHON 對(duì)縮進(jìn)式敏感的,空格不要隨意刪除！另外我用的python 2.5 的 2.6應(yīng)該可以運(yùn)行，3.0及以上應(yīng)該不可以！ )

#!/usr/bin/python
#-*- encoding: gbk -*-
from __future__ import division
from   struct   import   *
import os,time ,datetime,string,sys,math,re,shutil,glob
import zipfile,StringIO,getopt
#from readths2 import *
# 2010-09-02 by wanghp

basedir = r'e:\cwork\guosen' #如果你的安裝路徑不同,請(qǐng)改這里

exp_dir    = basedir + r'\T0002\export'
#exp_dir    = basedir + r'\T0002\export_back'
lc5_dir_sh = basedir + r'\Vipdoc\sh\fzline'
#lc5_dir_sh =  r'D:\2965\ydzqwsjy\Vipdoc\sh\fzline'
lc5_dir_sz = basedir + r'\Vipdoc\sz\fzline'
day_dir_sh = basedir + r'\Vipdoc\sh\lday'
day_dir_sz = basedir + r'\Vipdoc\sz\lday'

stkdict = {} #存儲(chǔ)股票ID和上海市、深圳市的對(duì)照

#############################################################
# read 通達(dá)信分筆數(shù)據(jù)
# example readfbtxt(readlines(),'20100831-600000.TXT')
# 返回的data格式為
# (stkid,datetime,price,amount,vol(股數(shù)),筆數(shù),buy or sale) 的list
#############################################################
def readfbtxt(p_lines,p_name):
    """讀通達(dá)信分筆數(shù)據(jù) """
    data = []
    shortname = os.path.split(p_name)[1]
    shortname = os.path.splitext(shortname)[0]
    sDay,stkID = shortname.split('-')
    if len(sDay) != 8 : return data
    stky = int(sDay[0:4])
    stkm = int(sDay[4:6])
    stkd = int(sDay[6:8])
    line_no = 0

    for l in p_lines:
        line_no += 1
        if line_no <=3: continue
        l = l.strip()
        t = re.split('\s+',l)
        k =  datetime.datetime(stky,stkm,stkd,int(t[0][0:2]),int(t[0][3:5]))
        p = float(t[1])      #price
        vol = int(t[2])*100  #股數(shù)
        amt = p * vol        #成交量
        bscnt = 0            #筆數(shù)
        bstag = ''           #buy or sale
        try:
            bscnt = int(t[3])    #筆數(shù)
            bstag = t[4]         #buy or sale
        except IndexError,e:
            pass
        data.append((stkID,k,p,amt,vol,bscnt,bstag))
    return data

#############################################################
# 將分筆數(shù)據(jù)轉(zhuǎn)化為分筆數(shù)據(jù)
# p_data:傳入?yún)?shù) 為readfbtxt所返回
# data:  返回的數(shù)據(jù)格式為
# [stkid,datetime,open,high,low,close,amt,vol(股)]
#############################################################
def fbtxt2lc0(p_data):
    """分筆數(shù)據(jù)轉(zhuǎn)化為1分鐘數(shù)據(jù)"""
    data = []
    for i in p_data:
        t = i[1]     #datetime
        p = i[2]     #price
        data.append([i[0],t,p,p,p,p,i[3],i[4]])
    return data

#############################################################
# 將分筆數(shù)據(jù)轉(zhuǎn)化為1分鐘數(shù)據(jù)
# p_data:傳入?yún)?shù) 為readfbtxt所返回
# data:  返回的數(shù)據(jù)格式為
# [stkid,datetime,open,high,low,close,amt,vol(股)]
#############################################################
def fbtxt2lc1(p_data):
    """分筆數(shù)據(jù)轉(zhuǎn)化為1分鐘數(shù)據(jù)"""
    data = []
    for i in p_data:
        t = i[1]     #datetime
        p = i[2]     #price
        lend = len(data)
        j = lend - 1
        while j >= 0:
            if data[j][1] == t:break
            j -= 1
        if j < 0:  #沒有找到該時(shí)間
            data.append([i[0],t,p,p,p,p,i[3],i[4]])
        else:          #找到該時(shí)間
            if p > data[j][3]:  #high
                data[j][3] = p
            if p < data[j][4]:  #low
                data[j][4] = p
            data[j][5] = p      #close
            data[j][6] += i[3]  #amout
            data[j][7] += i[4]  #vol
    #data.sort(key = lambda x:x[1])  #以datetime 排序
    return data

#############################################################
# 一個(gè)時(shí)間對(duì)應(yīng)的5分鐘區(qū)間段
# dt 傳入?yún)?shù) 為一個(gè)datetime.datetime or datetime.time
# 返回datetime 或time
#############################################################
def which5min(dt):
    """5 分鐘時(shí)間劃分 """
    if type(dt) != datetime.datetime and  type(dt) != datetime.time:
        return None
    t = dt
    ret = None
    if type(dt) == datetime.datetime:
        t = datetime.time(dt.hour,dt.minute,dt.second)

    if t < datetime.time(9,30) : return None
    if   t < datetime.time(9,35): ret = datetime.time(9,35)
    elif t < datetime.time(9,40): ret = datetime.time(9,40)
    elif t < datetime.time(9,45): ret = datetime.time(9,45)
    elif t < datetime.time(9,50): ret = datetime.time(9,50)
    elif t < datetime.time(9,55): ret = datetime.time(9,55)
    elif t < datetime.time(10,0): ret = datetime.time(10,0)
    elif t < datetime.time(10,5): ret = datetime.time(10,5)
    elif t < datetime.time(10,10): ret = datetime.time(10,10)
    elif t < datetime.time(10,15): ret = datetime.time(10,15)
    elif t < datetime.time(10,20): ret = datetime.time(10,20)
    elif t < datetime.time(10,25): ret = datetime.time(10,25)
    elif t < datetime.time(10,30): ret = datetime.time(10,30)
    elif t < datetime.time(10,35): ret = datetime.time(10,35)
    elif t < datetime.time(10,40): ret = datetime.time(10,40)
    elif t < datetime.time(10,45): ret = datetime.time(10,45)
    elif t < datetime.time(10,50): ret = datetime.time(10,50)
    elif t < datetime.time(10,55): ret = datetime.time(10,55)
    elif t < datetime.time(11,0): ret = datetime.time(11,0)
    elif t < datetime.time(11,5): ret = datetime.time(11,5)
    elif t < datetime.time(11,10): ret = datetime.time(11,10)
    elif t < datetime.time(11,15): ret = datetime.time(11,15)
    elif t < datetime.time(11,20): ret = datetime.time(11,20)
    elif t < datetime.time(11,25): ret = datetime.time(11,25)
    elif t <= datetime.time(11,30): ret = datetime.time(11,30)
#    elif t < datetime.time(13,0): ret = datetime.time(13,0)
    elif t < datetime.time(13,5): ret = datetime.time(13,5)
    elif t < datetime.time(13,10): ret = datetime.time(13,10)
    elif t < datetime.time(13,15): ret = datetime.time(13,15)
    elif t < datetime.time(13,20): ret = datetime.time(13,20)
    elif t < datetime.time(13,25): ret = datetime.time(13,25)
    elif t < datetime.time(13,30): ret = datetime.time(13,30)
    elif t < datetime.time(13,35): ret = datetime.time(13,35)
    elif t < datetime.time(13,40): ret = datetime.time(13,40)
    elif t < datetime.time(13,45): ret = datetime.time(13,45)
    elif t < datetime.time(13,50): ret = datetime.time(13,50)
    elif t < datetime.time(13,55): ret = datetime.time(13,55)
    elif t < datetime.time(14,0): ret = datetime.time(14,0)
    elif t < datetime.time(14,5): ret = datetime.time(14,5)
    elif t < datetime.time(14,10): ret = datetime.time(14,10)
    elif t < datetime.time(14,15): ret = datetime.time(14,15)
    elif t < datetime.time(14,20): ret = datetime.time(14,20)
    elif t < datetime.time(14,25): ret = datetime.time(14,25)
    elif t < datetime.time(14,30): ret = datetime.time(14,30)
    elif t < datetime.time(14,35): ret = datetime.time(14,35)
    elif t < datetime.time(14,40): ret = datetime.time(14,40)
    elif t < datetime.time(14,45): ret = datetime.time(14,45)
    elif t < datetime.time(14,50): ret = datetime.time(14,50)
    elif t < datetime.time(14,55): ret = datetime.time(14,55)
    elif t <= datetime.time(15,0):  ret = datetime.time(15,0)
    else : return None
    if type(dt) == datetime.datetime:
        return datetime.datetime(dt.year,dt.month,dt.day,ret.hour,ret.minute,ret.second)
    else: return ret

#############################################################
# 將1分鐘數(shù)據(jù)轉(zhuǎn)為5分鐘數(shù)據(jù)
# p_data:傳入?yún)?shù) 為fbtxt2lc1所返回
# data:  返回的數(shù)據(jù)格式為
# [stkid,datetime,open,high,low,close,amt,vol(股)]
#############################################################
def lc1tolc5(p_data):
    """1分鐘數(shù)據(jù)轉(zhuǎn)化為5分鐘數(shù)據(jù) """
    if len(p_data) <= 0: return None
    data = []
    for i in p_data:
        t = which5min(i[1])   #找對(duì)應(yīng)5分鐘的區(qū)段
        if t == None:
            raise ValueError,'time out of range: %s' % i[1]
        lend = len(data)
        j = lend - 1
        while j >= 0:
            if data[j][1] == t:break
            j -= 1
        if j < 0:  #沒有找到該時(shí)間
            data.append([i[0],t,i[2],i[3],i[4],i[5],i[6],i[7]])
        else:          #找到該時(shí)間
            if i[3] > data[j][3]:  #high
                data[j][3] = i[3]
            if i[4] < data[j][4]:  #low
                data[j][4] = i[4]
            data[j][5] = i[5]      #close
            data[j][6] += i[6]     #amout
            data[j][7] += i[7]     #vol
    #data.sort(key = lambda x:x[1])  #以datetime 排序
    return data

#############################################################
# read 5分鐘數(shù)據(jù)
# example readlc5(r'E:\new_gxzq_v6\Vipdoc\sh\fzline\sh600000.lc5')
#############################################################
def readlc5(p_name):
        """tdx 5min 數(shù)據(jù)
           日期上低16位表示月日，高16位表示分鐘
           這個(gè)結(jié)構(gòu)個(gè)人感覺就不如同花順做的巧妙
               在一個(gè)4字節(jié)中把 年 月 日 時(shí) 分 都記錄下來了
        """
        f = open(p_name,'rb')
        stkID = os.path.split(p_name)[1]
        stkID = os.path.splitext(stkID)[0]
        if string.lower(stkID[0:2]) == 'sh' or string.lower(stkID[0:2]) == 'sz':
            stkID = stkID[2:]
        icnt = 0
        data = []
        while 1:
            raw = f.read(4*8)
            if len(raw) <= 0 : break
            t = unpack('IfffffII',raw)
            mins = (t[0] >> 16) & 0xffff
            mds  = t[0] & 0xffff
            month = int(mds / 100)
            day   = mds % 100
            hour = int(mins / 60)
            minute = mins % 60
            #datet = "d-d d:d" % (month,day,hour,minute)
            data.append((stkID,(month,day,hour,minute),t[1],t[2],t[3],t[4],t[5],t[6],t[7]))
            #print datet,t[1],t[2],t[3],t[4],t[5],t[6],t[7]
            icnt += 1
        ## end while
        f.close()
        return data

#############################################################
#構(gòu)造通達(dá)信5min數(shù)據(jù)文件
# data 結(jié)構(gòu)
#[stkID,(月,日,時(shí),分),open,high,low,close,amt,vol,0]
#############################################################
def writelc5(p_name,data,addwrite = True):
    if addwrite :
        fout = open(p_name,'ab')
    else:
        fout = open(p_name,'wb')
    for i in data:
        t = i[1][0]*100+i[1][1] + ( (i[1][2] * 60 + i[1][3]) << 16)
        raw = pack('IfffffII',t,i[2],i[3],i[4],i[5],i[6],i[7],i[8])
        fout.write(raw)
    ## end for
    fout.close()

#############################################################
# outlist
# 將list 或者 tuple 輸出
# 遞歸垂直輸出，格式不好
#############################################################
def outlist(l):
    if type(l) != list and type(l) != tuple:
        print l
    else:
        for i in l:
            outlist(i)

#############################################################
# outlist2 二層輸出，對(duì)于上文中的data 實(shí)用
# 傳入?yún)?shù) data，一個(gè)形如矩陣的list
#############################################################
def outlist2(p_data):
    for i in p_data:
        for j in i:
            print j,
        print    

#############################################################
# 正數(shù)與二進(jìn)制表示互相表示
#
#############################################################
#------------------------------
#-- i2bin 整數(shù)轉(zhuǎn)為 2進(jìn)制字符串
#------------------------------
def i2bin(x):
        result = ''
        x = int(x)
        while x > 0:
               mod = x & 0x01 # 取2的余數(shù)
               x = x >> 0x01  # 右移一位
               result = str(mod) + result
        return result

#------------------------------
#-- bin2i 2進(jìn)制字符串 轉(zhuǎn)為正數(shù)
#------------------------------
def bin2i(bin):
    result = 0
    for s in bin:
        if s != '0' and s != '1':
            raise ValueError,'bad bin string:'+bin
        result = 2*result + int(s)
    return result

#############################################################
# fill_stkdict 填充全局變量字典 stkdict
#
#############################################################
def fill_stkdict():
    global stkdict
    lsh = os.listdir(day_dir_sh)
    for l in lsh:
        if len(l) <= 4:continue
        l = string.lower(l)
        if l[-3:] != 'day': continue
        n = os.path.splitext(l)[0]
        if n[0:2] == 'sh' or n[0:2] == 'sz':
            n = n[2:]
        stkdict[n] = 'sh'
    lsz = os.listdir(day_dir_sz)
    for l in lsz:
        if len(l) <= 4:continue
        l = string.lower(l)
        if l[-3:] != 'day': continue
        n = os.path.splitext(l)[0]
        if n[0:2] == 'sh' or n[0:2] == 'sz':
            n = n[2:]
        stkdict[n] = 'sz'

def getMarketByID(id):
    global stkdict
    if len(stkdict) == 0: fill_stkdict()
    return stkdict.setdefault(id,'')

#copy file   
#dstname = os.path.join(lc5_dir_sh,fout)
#shutil.copy(fout,dstname)

def writelcfiles(p_lines,p_name,lctype = 'lc5lc1',addfile = True):
    """
寫內(nèi)容
p_lines 文件行
l_name  短文件名
lctype  要轉(zhuǎn)化生成的分鐘類型 lc5 表示5分鐘 lc1 表示1分鐘 lc0表示分筆的
addfile True表示追加文件 False 表示覆蓋
    """
    data1  = readfbtxt(p_lines,p_name)
    if len(data1) ==0: return
    data2 = fbtxt2lc1(data1)
    data3 = lc1tolc5(data2)

    ## lc5 5分鐘文件
    if 'lc5' in lctype:
        data = []
        for i in data3:
            data.append([i[0],(i[1].month,i[1].day,i[1].hour,i[1].minute),i[2],i[3],i[4],i[5],i[6],i[7],0])
        if len(data) == 0 :
            sys.stderr.write('Error:no data in data\n')
            sys.exit(1)
        stkID = data[0][0]
        mark = getMarketByID(stkID)
        if mark == '':
            sys.stderr.write('不能確定它的市場：%s.請(qǐng)檢查代碼!\n' % stkID)
        else:
            fout = mark + stkID + '.lc5'
            if mark == 'sh':
                fout = os.path.join(lc5_dir_sh,fout)
            else:
                fout = os.path.join(lc5_dir_sz,fout)       

            writelc5(fout,data,addfile)
        #endif
    #endif
    if 'lc1' in lctype :
        data = []
        for i in data2:
            data.append([i[0],(i[1].month,i[1].day,i[1].hour,i[1].minute),i[2],i[3],i[4],i[5],i[6],i[7],0])
        if len(data) == 0 :
            sys.stderr.write('Error:no data in data\n')
            sys.exit(1)
        stkID = data[0][0]
        mark = getMarketByID(stkID)
        if mark == '':
            sys.stderr.write('不能確定它的市場：%s.請(qǐng)檢查代碼!\n' % stkID)
        else:
            fout = mark + stkID + '.lc1'
            if mark == 'sh':
                fout = os.path.join(lc5_dir_sh,fout)
            else:
                fout = os.path.join(lc5_dir_sz,fout)       
   
            writelc5(fout,data,addfile)
        #endif.
    #endif.

    #lc0 分筆的K線文件
    if 'lc0' in lctype:
        data0 = fbtxt2lc0(data1)
        data = []
        for i in data0:
            data.append([i[0],(i[1].month,i[1].day,i[1].hour,i[1].minute),i[2],i[3],i[4],i[5],i[6],i[7],0])
        if len(data) == 0 :
            sys.stderr.write('Error:no data in data\n')
            sys.exit(1)
        stkID = data[0][0]
        mark = getMarketByID(stkID)
        if mark == '':
            sys.stderr.write('不能確定它的市場：%s.請(qǐng)檢查代碼!\n' % stkID)
        else:
            fout = mark + stkID + '.lc0'
            if mark == 'sh':
                fout = os.path.join(lc5_dir_sh,fout)
            else:
                fout = os.path.join(lc5_dir_sz,fout)       
   
            writelc5(fout,data,addfile)
        #endif.
    #endif.

def convert(p_stkid,p_type = 'txt',filterfunc = None):
    if p_type == 'txt' : #txt file
        txtfiles = glob.glob(os.path.join(exp_dir,'*-'+p_stkid+'.txt'))
        if filterfunc :
            txtfiles = filter(filterfunc,txtfiles)
        txtfiles.sort()
        l_i = 0
        for fname in txtfiles:
            sys.stderr.write('%s\n' % fname)
            try:
                doc_lines = file(fname).readlines()
            except IOError , e :
                sys.stderr.write("Open file %s fail!\n" % fname)
                continue
            if len(doc_lines) <=4:
                sys.stderr.write('No data in %s\n' % fname)
                continue
            if l_i == 0 :
                writelcfiles(doc_lines,fname,'lc5lc1',False) ##over write
            else:
                writelcfiles(doc_lines,fname,'lc5lc1',True)  ##add write
            l_i += 1
        #endfor
           
    else:                #zipfile
        try:
            fzip = zipfile.ZipFile(os.path.join(exp_dir,p_stkid+'.zip'))
        except IOError,e:
            print 'Can not open file!',e
            return
        zipedfiles = fzip.namelist()
        if len(zipedfiles) == 0:
            return
        if filterfunc :
            zipedfiles = filter(filterfunc,zipedfiles)
        zipedfiles.sort()
        l_i = 0
        for fname in zipedfiles:
            sys.stderr.write('%s\n' % fname)
            doc = fzip.read(fname)
            doc_lines = StringIO.StringIO(doc).readlines()       
            if l_i == 0 :
                writelcfiles(doc_lines,fname,'lc5lc1',False) ##over write
            else:
                writelcfiles(doc_lines,fname,'lc5lc1',True)  ##add write
            l_i += 1
        #endfor
    #endif.

#############################################################
# usage 使用說明
#
#############################################################

def usage(p):
    print """
python %s [-t txt|zip] stkid [from] [to]
-t txt 表示從txt files 讀取數(shù)據(jù)，否則從zip file 讀取(這也是默認(rèn)方式)
for example :
python %s 999999 20070101 20070302
python %s -t txt 999999 20070101 20070302
    """ % (p,p,p)

if __name__ == '__main__':
    """
    python readtdxlc5.py 999999 20070101 20070131
    """
    argv = sys.argv[1:]
    try :
        opts, args = getopt.getopt(argv, "ht:", ["help", "type="])
    except getopt.GetoptError:
        usage(sys.argv[0])
        sys.exit(0)
    l_type = 'zip'  #default type is zipfiles!
    for opt, arg in opts:               
        if opt in ("-h", "--help"):
            usage(sys.argv[0])
            sys.exit(1)
        elif opt in ("-t", "--type"):
            l_type = arg
    if len(args) <1:
        print 'You must specified the stock No.!'
        usage(sys.argv[0])
        sys.exit(1)

    stkid = args[0]
    l_from = None
    l_to = None
    try:
        l_from = args[1]
        l_to   = args[2]
    except :
        pass

    #過濾函數(shù)
    def filfunc(x):
        if l_from == None and l_to == None:
            return True
        ymd = os.path.splitext(os.path.split(x)[1])[0].split('-')[0]
        if l_from and l_to:
            return ymd >= l_from and ymd <= l_to
        elif l_from :
            return ymd >= l_from
        else :
            return ymd <= l_to
   
    if l_type == 'txt': # 從一般txt 文件
        convert(stkid,'txt',filfunc)
    else:
        convert(stkid,'zip',filfunc)
       
    mark = getMarketByID(stkid)
    if mark == '':
        sys.stderr.write('不能確定它的市場：%s.請(qǐng)檢查代碼!\n' % stkid)       
    else:
        os.system('copy E:\\cwork\\guosen\\Vipdoc\\'+ mark +'\\fzline\\'+ mark + stkid +'.lc5 E:\\cwork\\my_yd\\Vipdoc\\'+ mark +'\\fzline\\')
        os.system('copy E:\\cwork\\guosen\\Vipdoc\\'+ mark +'\\fzline\\'+ mark + stkid +'.lc1 E:\\cwork\\ydzqwsjy\\Vipdoc\\' + mark + '\\fzline\\' + mark + stkid +'.lc5')

    #data = readlc5(os.path.join(lc5_dir_sh,'sh601398.lc5'))
    #outlist2(data)