python 中 对比两序列，若字符一样在二维矩阵中将相应位置赋值为1
001、
>>> import openpyxl ## 导入包 >>> import numpy as np >>> import pandas as pd >>> s1=list('ATGATAGCAGTGAAATGGG') ## 将序列储存为列表 >>> s1 ['A', 'T', 'G', 'A', 'T', 'A', 'G', 'C', 'A', 'G', 'T', 'G', 'A', 'A', 'A', 'T', 'G', 'G', 'G'] >>> s2=list('GATAGCAGTGAAACGGGCA') >>> s2 ['G', 'A', 'T', 'A', 'G', 'C', 'A', 'G', 'T', 'G', 'A', 'A', 'A', 'C', 'G', 'G', 'G', 'C', 'A'] >>> len(s1); len(s2) 19 19 >>> compare = [[0]*19 for i in range(19)] ## 生成二维数组 >>> compare [[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]] >>> for i in range(19): ... for j in range(19): ... if s1[i] == s2[j]: ... compare[i][j] = 1 ... >>> compare [[0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0]] >>> a = np.array(compare) ##将数组转换为矩阵形式 >>> a array([[0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0], [1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0]]) >>> data_df = pd.DataFrame(a) ## 转换为数据框 >>> data_df 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 0 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 2 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 3 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 4 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 5 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 6 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 7 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 1 0 8 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 9 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 10 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 11 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 12 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 13 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 14 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 15 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 16 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 17 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 18 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 >>> data_df.columns = s2 ## 重命名行名和列名 >>> data_df.index = s1 >>> data_df G A T A G C A G T G A A A C G G G C A A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 T 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 T 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 C 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 1 0 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 T 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 A 0 1 0 1 0 0 1 0 0 0 1 1 1 0 0 0 0 0 1 T 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0 G 1 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 0
>>> data_df.to_csv("xxx.csv")
(base) root@PC1:/home/test2# cat xxx.csv ,G,A,T,A,G,C,A,G,T,G,A,A,A,C,G,G,G,C,A A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 T,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 T,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 C,0,0,0,0,0,1,0,0,0,0,0,0,0,1,0,0,0,1,0 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 T,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 A,0,1,0,1,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1 T,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0 G,1,0,0,0,1,0,0,1,0,1,0,0,0,0,1,1,1,0,0
参考：https://blog.csdn.net/m0_57099761/article/details/123464340
