求助：Stereonet图中Pole动画显示位置异常的解决方法

我刚帮你排查了代码里的核心问题——你在动画函数里直接把**原始的走向(strike)和倾角(dip)**传给了poles.set_data()，但忽略了一个关键细节：mplstereonet的pole()方法会自动将地质领域的strike/dip坐标转换成立体网的投影坐标，而set_data()需要的是转换后的平面坐标，不是原始的角度值。

具体解决思路 #

• 先预计算所有极点对应的立体网投影坐标，避免在动画中重复做转换逻辑
• 在动画更新函数中，用预计算好的投影坐标来更新极点位置，保证和静态图的坐标体系一致

修正后的完整代码 #

import matplotlib as mpl
mpl.use("TkAgg")
from matplotlib import pyplot as plt
from matplotlib import animation
import numpy as np
import mplstereonet

fig, ax = mplstereonet.subplots()
fig2, ax1 = mplstereonet.subplots()
ax.grid(True)
ax1.grid(True)

# 假设带有随机偏差的走向（strike）和倾角（dip）
# 当前参数应将极点绘制在0°或180°位置
strike, dip = 90, 80
num = 10
strikes = strike + 10 * np.random.randn(num)
dips = dip + 10 * np.random.randn(num)

# 预计算所有极点的立体网投影坐标
# 1. 先把strike/dip转为极点的方位角和极距
pole_azimuths, pole_distances = mplstereonet.pole(strikes, dips)
# 2. 把极坐标转为立体网的平面直角坐标
x, y = mplstereonet.polar2xy(pole_azimuths, pole_distances)

# 初始化空的极点绘图对象（只创建artist，不传入初始数据）
poles, = ax.plot([], [], 'o')

def init():
    poles.set_data([], [])
    return poles,

def animate(i):
    # 使用预计算好的投影坐标更新极点位置
    poles.set_data(x[:i], y[:i])
    return poles,

anim = animation.FuncAnimation(fig, animate, init_func=init, frames=len(strikes), interval=100, blit=True, repeat=False)
# 绘制预期的静态对比图
ax1.pole(strikes, dips, 'o')

plt.show()

关键细节解释 #

1. 坐标转换的必要性：ax.pole(strikes, dips)内部会自动完成两次转换——先将strike/dip转为极点的极坐标（方位角+极距），再转为立体网的平面直角坐标。你之前直接把原始strike/dip传给set_data，相当于把角度值当成了平面坐标系的x/y，自然位置完全偏离预期。
2. 预计算的优势：提前计算所有坐标既可以避免动画过程中重复计算，提升运行效率，还能保证动画的最终帧和静态对比图完全一致。

现在运行修正后的代码，动画里的极点位置就会和你预期的静态图完全匹配了。

内容的提问来源于stack exchange，提问作者Aly Abdelaziz