python

关注公众号 jb51net

关闭
首页 > 脚本专栏 > python > Python绘制折线图

Python绘制三维填充折线图的示例代码

作者:Code_Verse

在数据可视化领域,三维图形能够以更直观的方式展示数据之间的复杂关系,本文将为大家详细介绍如何使用Python绘制三维填充折线图,需要的小伙伴可以了解下

在数据可视化领域,三维图形能够以更直观的方式展示数据之间的复杂关系。今天,我们将从头开始,一步步学习如何使用Python中的matplotlib库绘制一个漂亮的三维填充折线图。即使你之前没有接触过Python绘图,也不用担心,我们将详细解释每一步。

一、准备工作:安装必要的库

在开始之前,我们需要确保安装了numpymatplotlib这两个Python库。如果你还没有安装它们,可以通过以下命令安装:

pip install numpy matplotlib

这两个库是Python数据科学领域中非常基础且强大的工具,numpy用于高效处理数值数据,而matplotlib则用于数据可视化。

二、代码解析:一步步搭建三维填充折线图

1. 导入必要的模块

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d.art3d import PolyCollection
import matplotlib.cm as cm
import matplotlib

2. 设置中文和负号显示

matplotlib.rcParams['font.sans-serif'] = ['SimHei']
matplotlib.rcParams['axes.unicode_minus'] = False

这两行代码是为了确保在图表中可以正常显示中文字符和负号。如果你的系统中没有SimHei字体,可以替换为其他支持中文的字体。

3. 模拟数据

num_samples = 5
mz = np.linspace(0, 100, 300)
spectra = []
np.random.seed(0)
for i in range(num_samples):
    peak_positions = np.random.uniform(20, 80, 3)
    peak_heights = np.random.uniform(0.8, 1.5, 3)
    spectrum = sum(h * np.exp(-0.5 * ((mz - p)/3)**2) for p, h in zip(peak_positions, peak_heights))
    spectrum += np.random.normal(0, 0.05, len(mz))
    spectra.append(spectrum)

这里我们生成了一些模拟数据,用来展示图表的效果。num_samples表示我们有5个样本,mz是一个从0到100的线性分布的数组,表示质荷比(M/Z)。spectra是一个列表,存储每个样本的光谱数据。

4. 创建图形

fig = plt.figure(figsize=(12, 7), facecolor='white')
ax = fig.add_subplot(111, projection='3d')

5. 设置颜色和样本标签

colors = cm.viridis(np.linspace(0.2, 0.8, num_samples))
sample_labels = ['样本A', '样本B', '样本C', '样本D', '样本E']

6. 构建多边形区域

verts = []
for spectrum in spectra:
    verts.append(list(zip(mz, spectrum)))
poly = PolyCollection(
    verts,
    facecolors=colors,
    edgecolors='k',         # 添加黑色边框
    linewidths=0.6,         # 细边框
    alpha=0.85              # 稍微增强填充透明度
)
ax.add_collection3d(poly, zs=range(num_samples), zdir='y')

7. 设置坐标轴和标题

ax.set_xlabel('质荷比 (M/Z)', labelpad=12, fontsize=12, fontweight='bold')
ax.set_xlim(0, 100)
ax.set_ylabel('样本', labelpad=12, fontsize=12, fontweight='bold')
ax.set_ylim(-0.5, num_samples - 0.5)
ax.set_yticks(range(num_samples))
ax.set_yticklabels(sample_labels, fontsize=11)
ax.set_zlabel('离子强度', labelpad=10, fontsize=12, fontweight='bold')
ax.set_zlim(0, 1.6)
ax.set_title('提取的质谱子集', fontsize=14, fontweight='bold', pad=20)

这里我们设置了坐标轴的标签、范围、刻度标签和标题。labelpad参数用于调整标签的位置,fontsizefontweight分别设置字体大小和粗细。

8. 美化网格线和调整边距

ax.tick_params(axis='both', which='major', labelsize=10)
ax.grid(True, linestyle='--', linewidth=0.5, alpha=0.4)
fig.subplots_adjust(left=0.08, right=0.95, top=0.92, bottom=0.08)

9. 显示图形

plt.show()

最后,调用plt.show()就可以显示我们绘制的三维填充折线图了。

三、运行代码并查看结果

当你运行上述代码时,你会看到一个三维填充折线图。每个折线代表一个样本的光谱数据,填充区域的颜色表示不同的样本,如下图所示。

完整代码如下:

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d.art3d import PolyCollection
import matplotlib.cm as cm
import matplotlib
# 设置中文和负号显示
matplotlib.rcParams['font.sans-serif'] = ['SimHei']
matplotlib.rcParams['axes.unicode_minus'] = False
# 模拟数据
num_samples = 5
mz = np.linspace(0, 100, 300)
spectra = []
np.random.seed(0)
for i in range(num_samples):
    peak_positions = np.random.uniform(20, 80, 3)
    peak_heights = np.random.uniform(0.8, 1.5, 3)
    spectrum = sum(h * np.exp(-0.5 * ((mz - p)/3)**2) for p, h in zip(peak_positions, peak_heights))
    spectrum += np.random.normal(0, 0.05, len(mz))
    spectra.append(spectrum)
# 创建图形
fig = plt.figure(figsize=(12, 7), facecolor='white')
ax = fig.add_subplot(111, projection='3d')
# 使用viridis色图
colors = cm.viridis(np.linspace(0.2, 0.8, num_samples))
sample_labels = ['样本A', '样本B', '样本C', '样本D', '样本E']
# 构建多边形区域
verts = []
for spectrum in spectra:
    verts.append(list(zip(mz, spectrum)))
poly = PolyCollection(
    verts,
    facecolors=colors,
    edgecolors='k',         
    linewidths=0.6,         
    alpha=0.85              
)
ax.add_collection3d(poly, zs=range(num_samples), zdir='y')
# 设置坐标轴
ax.set_xlabel('质荷比 (M/Z)', labelpad=12, fontsize=12, fontweight='bold')
ax.set_xlim(0, 100)
ax.set_ylabel('样本', labelpad=12, fontsize=12, fontweight='bold')
ax.set_ylim(-0.5, num_samples - 0.5)
ax.set_yticks(range(num_samples))
ax.set_yticklabels(sample_labels, fontsize=11)
ax.set_zlabel('离子强度', labelpad=10, fontsize=12, fontweight='bold')
ax.set_zlim(0, 1.6)
# 标题
ax.set_title('提取的质谱子集', fontsize=14, fontweight='bold', pad=20)
# 网格线美化
ax.tick_params(axis='both', which='major', labelsize=10)
ax.grid(True, linestyle='--', linewidth=0.5, alpha=0.4)
# 调整边距
fig.subplots_adjust(left=0.08, right=0.95, top=0.92, bottom=0.08)
# 显示图
plt.show()

四、总结

通过以上步骤,我们用Python和matplotlib库成功绘制了一个三维填充折线图。虽然代码看起来有点长,但其实每一步都很简单。你可以根据自己的需求修改数据、颜色和样式,制作出属于自己的图表。

到此这篇关于Python绘制三维填充折线图的示例代码的文章就介绍到这了,更多相关Python绘制折线图内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

您可能感兴趣的文章:
阅读全文