调整Matplotlib极坐标图以显示亚度级运动（即拉伸极坐标图切片）

我完全懂你的困扰——当theta的变化只有0.01度这种亚度级水平时，直接缩窄极坐标的角度范围会把目标区域压成一条几乎看不出细节的细线，根本没法观察数据的变化趋势。咱们可以通过把微小的真实角度区间线性映射到一个更宽的显示角度区间来实现"拉伸放大"，既能保留真实的rho数值，又能让theta的细微变化清晰可见，同时还能让轴刻度对应真实的角度值，不会造成混淆。

核心思路 #

简单来说，就是把你那只有一点点变化的真实theta值，映射到一个你指定的宽角度区间（比如300°到340°），让极坐标轴用这个宽区间来展示原本微小的变化，相当于把细节"拉宽"了。之后再自定义轴刻度，让显示的刻度标签对应真实的角度数值，保证可读性。

修改后的完整代码 #

import numpy as np
import matplotlib.pyplot as plt
import matplotlib

print('matplotlib version : ', matplotlib.__version__)

# 原始观测数据
ra = np.asarray([1.67484,1.67485,1.67485,1.67486,1.67486,1.67488,1.67487,1.67488,1.67487, 1.67487]) # 弧度
dec = np.asarray([-0.92147,-0.92147,-0.92147,-0.92147,-0.92147,-0.92147,-0.92147, -0.92147,-0.92147, -0.92147]) # 弧度

rho = np.sqrt(ra**2 + dec**2) # 计算rho
theta = np.arctan2(dec, ra) # 计算真实theta（弧度）

# --------------------------
# 关键步骤：将真实theta映射到宽显示区间
# --------------------------
# 先算出真实theta的实际范围
theta_min_real = theta.min()
theta_max_real = theta.max()
delta_theta_real = theta_max_real - theta_min_real

# 定义我们想要用来展示的宽角度区间（单位：度，转成弧度）
# 这里选300°到340°，你可以根据自己的需求随意调整
thetamin_display = np.deg2rad(300)
thetamax_display = np.deg2rad(340)
delta_theta_display = thetamax_display - thetamin_display

# 线性映射：把真实的小角度区间，映射到我们定义的宽显示区间
theta_mapped = thetamin_display + (theta - theta_min_real) * delta_theta_display / delta_theta_real

# --------------------------
# 绘图与轴参数设置
# --------------------------
fig = plt.figure(figsize=(8,6))
ax = fig.add_subplot(1,1,1,polar=True)

# 使用映射后的theta绘制数据
ax.plot(theta_mapped, rho, '*-', color='y', markersize=8)

# 保留你原来设置的rho显示范围
ax.set_ylim(1.9114, 1.9117)

# 设置显示的角度范围，对应我们定义的宽区间
ax.set_thetamin(np.rad2deg(thetamin_display))
ax.set_thetamax(np.rad2deg(thetamax_display))

# 自定义刻度：让显示的刻度对应真实的角度值
# 生成几个真实theta的刻度点（转成角度）
real_theta_ticks = np.linspace(theta_min_real, theta_max_real, 5)
real_theta_ticks_deg = np.rad2deg(real_theta_ticks)
# 计算这些真实刻度在映射后的位置
mapped_ticks = thetamin_display + (real_theta_ticks - theta_min_real) * delta_theta_display / delta_theta_real
# 替换轴刻度与标签
ax.set_xticks(mapped_ticks)
ax.set_xticklabels([f"{t:.6f}°" for t in real_theta_ticks_deg])

# 加个标题说明
ax.set_title("Stretched Polar Plot (Real Theta Stretched to 300°-340° Display)", y=1.1)

plt.show()

效果说明 #

这样处理后，原本只有~0.01度的theta变化会被拉伸到40度的显示区间，你能非常清晰地看到每个数据点的theta变化趋势，同时轴上的刻度显示的是真实的角度值，不会造成混淆。你可以根据自己的需求随意调整thetamin_display和thetamax_display的数值，来控制拉伸的程度——比如想要更宽的显示区间，改成270°到350°也是完全可以的。

备注：内容来源于stack exchange，提问作者tronbons