python绘制弦图-科研作图

一、背景

弦图以其直观、精美的展示方式受到越来越多人的关注,它不仅能够有效展示两个变量之间的联系,还能同时展现多个变量间的复杂互动,本文将通过Python语言中的pycirclize库,带你深入了解如何绘制弦图。

弦图是一种圆形可视化方法,通常用于展示两组对象之间的相互联系和权重,弧线代表对象,弦代表对象之间的连接,弦的宽度表示关系的强弱,弦图直观地反映了多个对象之间的复杂联系,使其非常适用于社交网络、基因表达、系统架构等领域的数据展示。

二、代码实例。

from pycirclize import Circos
import pandas as pd

#定义行名称与列名称
rows_name = ['S1', 'S2', 'S3']
cols_name = ['E1', 'E2', 'E3', 'E4', 'E5', 'E6']

#定义矩阵数据
matrix_data = [
    [4,14,13,17,5,2],
    [5,12,9,15,3,7],
    [6,8,11,19,1,4]
]

#创建数据框,用于表示矩阵数据,行名称与列名称分别为rows_name与cols_name
matrix_df = pd.DataFrame(matrix_data, index=rows_name, columns=cols_name)
matrix_df
#创建弦图对象
circos = Circos.initialize_from_matrix(
    matrix_df,#矩阵数据
    start=-265,#起始角度
    end=95,#结束角度
    r_lim=(93,100),#环形范围
    cmap='tab10',#颜色映射
    label_kws=dict(
        r=94,#标签半径
        size=12,#标签大小
        color='white'#标签颜色
    ),
    link_kws=dict(
        ec='black',#弦的颜色
        lw=0.5#弦的宽度
    )
)
print(matrix_df)
#绘制弦图
fig=circos.plotfig()
fig.savefig('弦图1.pdf',format='pdf',bbox_inches='tight',dpi=1200)
from pycirclize import Circos
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
# 生成一个10x10的随机整数矩阵
row_names = [f"S{i+1}" for i in range(10)]  # 行名称,S1 到 S10
col_names = [f"E{i+1}" for i in range(10)]  # 列名称,E1 到 E10
matrix_data = np.random.randint(1, 20, size=(10, 10))  # 生成随机1到20之间的整数
# 将矩阵数据转换为 DataFrame,以方便后续操作
matrix_df = pd.DataFrame(matrix_data, index=row_names, columns=col_names)
# 初始化 Circos 图
circos = Circos.initialize_from_matrix(    matrix_df,             # 传入矩阵数据    
                            space=5,               # 每个单元之间的间隔设为5度    
                            r_lim=(93, 100),       # 设置圆形图的半径范围为93到100    
                            cmap='tab10',          # 使用颜色映射方案 'tab10',可以生成多种颜色    
                            ticks_interval=30,     # 设置刻度间隔,每隔20个单位标记一个刻度    
                            label_kws=dict(        
                                        r=94,              # 标签的半径距离为94        
                                        size=12,           # 标签的字体大小设为12        
                                        color='white'      # 标签的字体颜色为白色    
                            ),    
                            link_kws=dict(        
                                        ec='black',        # 弦的边框颜色为黑色        
                                        lw=0.5             # 弦的线宽设为0.5    
                            )
                        )
# 绘制 Circos 图
fig = circos.plotfig()
print(matrix_df)
fig.savefig("弦图2.pdf", format='pdf', bbox_inches='tight', dpi=1200)
from pycirclize.parser import Matrix
import pandas as pd
from pycirclize import Circos
# 创建一个包含 'from'、'to' 和 'value' 列的 DataFrame# 该表表示从某个节点到另一个节点的联系及其权重
fromto_table_df = pd.DataFrame([ 
                        ["A", "B", 10],  # A到B的权重为10    
                        ["A", "C", 5],   # A到C的权重为5    
                        ["A", "D", 15],  # A到D的权重为15    
                        ["A", "E", 20],  # A到E的权重为20    
                        ["A", "F", 3],   # A到F的权重为3    
                        ["B", "A", 3],   # B到A的权重为3    
                        ["B", "G", 15],  # B到G的权重为15    
                        ["F", "D", 13],  # F到D的权重为13    
                        ["F", "E", 13],  # F到E的权重为13    
                        ["E", "A", 20],  # E到A的权重为20    
                        ["E", "D", 6]    # E到D的权重为6
                      ], 
                        columns=["from", "to", "value"])  # 设置列名为 'from','to','value'
# 使用 Matrix 类的 parse_fromto_table 方法将 'from-to' 表解析为矩阵
matrix = Matrix.parse_fromto_table(fromto_table_df)
# 使用解析后的矩阵初始化 Circos 图
circos = Circos.initialize_from_matrix( matrix,              # 传入解析后的矩阵    
                                         space=3,             # 每个单元之间的间隔设置为3度    
                                         cmap="viridis",      # 使用 'viridis' 颜色映射表    
                                         ticks_interval=5,    # 设置刻度间隔为5    
                                         label_kws=dict(      # 设置标签的关键字参数        
                                         size=12,         # 标签字体大小为12       
                                         r=110            # 标签距离圆心的半径为110    
                                         ),    
                                link_kws=dict(       # 设置弦的样式        
                                         direction=1,     # 设定弦的方向        
                                         ec="black",      # 弦的边框颜色为黑色       
                                         lw=0.5           # 弦的线宽设置为0.5    
                                )
)
# 打印 from-to 表的内容,方便调试或查看
print(fromto_table_df.to_string(index=False))# 绘制 Circos 图
fig = circos.plotfig()
fig.savefig("弦图3.pdf", format='pdf', bbox_inches='tight', dpi=1200)
相关推荐
亚图跨际15 分钟前
Python和R基因组及蛋白质组学和代谢组学
python·r语言·生物医学
程序员与背包客_CoderZ17 分钟前
C++设计模式——Singleton单例模式
c语言·开发语言·c++·单例模式·设计模式
fanyamin22 分钟前
编程语言的局限
开发语言·python
小猿_0024 分钟前
C语言蓝桥杯组题目
c语言·开发语言·蓝桥杯
落笔映浮华丶30 分钟前
C++(进阶) 第1章 继承
开发语言·c++
百香果果ccc40 分钟前
Java中的File和IO流
java·开发语言·log4j
努力更新中1 小时前
Python浪漫之随机绘制不同颜色的气球
开发语言·python
委婉待续1 小时前
java抽奖系统(一)
java·开发语言
__lost1 小时前
Python 将彩色视频转换为黑白视频(MP4-格式可选)
python·opencv·音视频
Evand J1 小时前
【MATLAB蓝牙定位代码】三维平面定位设计,通过N个蓝牙锚点实现对未知位置的精准定位
开发语言·matlab·平面