在Matplotlib的3D散点图中清晰地标注特征向量方向
本文介绍如何在Matplotlib的3D散点图中，使用箭头直观地显示特征向量的方向。 许多用户在绘制3D数据及其特征向量后，难以直接判断向量方向。本文提供一种解决方案，通过自定义箭头，增强可视化效果。
Matplotlib本身并不直接支持在3D图中绘制指向已有3D线段的箭头。虽然FancyArrowPatch类可以创建箭头，但它仅限于2D绘图。因此，我们需要创建一个自定义类Arrow3D，继承自FancyArrowPatch并重写其draw方法。
Arrow3D类的核心在于将3D坐标投影到2D坐标系。这通过mpl_toolkits.mplot3d.proj3d.proj_transform()函数实现。该函数将3D坐标(x, y, z)转换为2D坐标(x', y')，然后Arrow3D类使用这些2D坐标设置箭头的起点和终点，最终在2D绘图区域绘制出3D向量对应的箭头。
以下代码展示了Arrow3D类的实现及其应用：
import numpy as np
from matplotlib import pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.patches import FancyArrowPatch
from mpl_toolkits.mplot3d import proj3d
class Arrow3D(FancyArrowPatch):
def __init__(self, xs, ys, zs, *args, **kwargs):
FancyArrowPatch.__init__(self, (0,0), (0,0), *args, **kwargs)
self._verts3d = xs, ys, zs
def draw(self, renderer):
xs3d, ys3d, zs3d = self._verts3d
xs, ys, zs = proj3d.proj_transform(xs3d, ys3d, zs3d, renderer.M)
self.set_positions((xs[0],ys[0]),(xs[1],ys[1]))
FancyArrowPatch.draw(self, renderer)
# ... (假设您已计算出特征向量 eig_vec 和数据中心点 mean_x, mean_y, mean_z) ...
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
# ... (绘制您的3D散点图) ...
for v in eig_vec:
a = Arrow3D([mean_x, mean_x + v[0]], [mean_y, mean_y + v[1]],
[mean_z, mean_z + v[2]], mutation_scale=20,
lw=2, arrowstyle="-|>", color="r") #调整lw和mutation_scale控制箭头粗细和大小
ax.add_artist(a)
# ... (其余绘图代码) ...
plt.show()
通过用Arrow3D类替换之前的线段绘制函数，并将其添加到子图中，即可在3D图中显示带有箭头的特征向量，使特征向量方向更加清晰明了。 此方法利用Matplotlib 3D绘图的底层机制，将3D投影到2D进行绘制。 请确保您已正确计算出特征向量eig_vec和数据中心点mean_x, mean_y, mean_z。 lw和mutation_scale参数可以调整箭头的粗细和箭头的大小。
