期刊图片复现|Python绘制Nature和弦图

    #绘制扇区和特征标签    for sector in circos.sectors:  #遍历每一个扇区        #获取颜色        fill_color = color_dict.get(sector.name, "#CCCCCC")        #绘制半径93到100的主颜色环，带黑色细边框        sector.rect(r_lim=(93, 100)        #绘制外圈阴影圆环        sector.rect(r_lim=(101, 108),                    facecolor="#F0F0F0",                    edgecolor="#D0D0D0",                    linewidth=0)        label_text = sector.name  #获取当前扇区的名称作为特征标签文本        #绘制特征标签文本        sector.text(label_text,                    r=104,                    orientation="vertical",                    size=14,                    color="black")

    #绘制连线    for i in range(N):  #遍历行            #只有当数值大于0时才绘制连线            if value > 0:                source_label = labels[i]  #获取源节点标签                target_label = labels[j]  #获取目标节点标签                src_start = df.iloc[i, :j].sum()  #源节点在扇区内的起始弧度位置                color = color_dict.get(source_label, "#888888")                #绘制连线                circos.link(                    (source_label, src_start, src_end),  #源点坐标元组（标签，起点，终点）                    (target_label, tgt_start, tgt_end),  #终点坐标元组                    color=color,  #连线颜色                    alpha=0.6,  #连线透明度                    direction=0  #连线方向参数                )

    circos.plotfig()  #生成并绘制最终的图形对象    #设置标题    plt.title("Circos Graph",              fontsize=18,              pad=80)    #保存    plt.savefig(fr"chord_diagram_{scheme_index}.png", dpi=300, bbox_inches="tight")    plt.savefig(fr"chord_diagram_{scheme_index}.pdf", bbox_inches="tight")

# =========================================================================================# ======================================4.执行部分=======================================# =========================================================================================if __name__ == '__main__':    data = r"data.xlsx"  #输入数据    #读取数据    df= pd.read_excel(data, index_col=0)    labels = df.index.tolist()  #获取数据框的索引列表，作为节点的标签    draw_chord_diagram(df,labels,scheme_id=scheme_index)  #调用绘图函数

scheme_index =20  #选择配色方案

plt.savefig(fr"chord_diagram_{scheme_index}.png", dpi=300, bbox_inches="tight")plt.savefig(fr"chord_diagram_{scheme_index}.pdf", bbox_inches="tight")

data = r"data.xlsx"  #输入数据