process_bam.py 文件源码

def write_psites(tpsites,psites_number,filename):
    with h5py.File(filename,"w") as fout:
        ds = h5py.special_dtype(vlen=str)
        dt = h5py.special_dtype(vlen=np.dtype("int32"))
        fout.create_dataset("transcript_ids",data=tpsites.keys(),dtype=ds)
        fout.create_dataset("p_sites",data=tpsites.values(),dtype=dt, compression="gzip")
        fout.create_dataset("psites_number",data=psites_number,dtype="int32")
    return None