Matplotlib设置3D图对数刻度后出现近乎空白的异常图问题

问题原因 #

Matplotlib的3D轴对对数刻度的支持不如2D轴完善，触发空白图的常见情况有两种：

• 数据中存在0或负数：对数函数无法处理非正值，会导致轴范围计算异常
• 自动轴范围适配失败：3D轴在对数模式下，自动计算的刻度范围可能完全偏离数据实际分布

解决方案 #

1. 先检查数据合法性

确保folder_data_real里的x、y、z数据全部为正值，可以用以下代码快速验证：

# 检查某一组数据是否存在非正值
x_data = folder_data_real[it][0]
y_data = folder_data_real[it][1]
z_data = folder_data_real[it][2]
print("x非正值数量:", sum(x <= 0 for x in x_data))
print("y非正值数量:", sum(y <= 0 for y in y_data))
print("z非正值数量:", sum(z <= 0 for z in z_data))

如果有非正值，优先过滤掉这些数据；若必须保留，可给数据加一个极小的正数（但可能影响数据解读）。

2. 手动指定对数轴的范围

如果数据全为正但还是空白，手动设置轴的上下限可以强制Matplotlib显示正确的范围：

# 计算数据的最小/最大值（取对数前的原始值）
x_min, x_max = min(x_data), max(x_data)
y_min, y_max = min(y_data), max(y_data)
z_min, z_max = min(z_data), max(z_data)

# 设置对数轴范围（用原始值，不是对数后的值）
ax.set_xlim(x_min, x_max)
ax.set_ylim(y_min, y_max)
ax.set_zlim(z_min, z_max)

# 再设置对数刻度
ax.set_xscale('log')
ax.set_yscale('log')
ax.set_zscale('log')

3. 替代方案：自定义对数刻度标签

如果上述方法无效，可以绕过set_xscale，手动设置刻度和标签模拟对数效果：

import numpy as np

# 以x轴为例，生成对数刻度位置（比如10^1,10^2,...）
x_ticks = np.logspace(np.log10(x_min), np.log10(x_max), num=5)
ax.set_xticks(x_ticks)
ax.set_xticklabels([f"{tick:.0e}" for tick in x_ticks])

# y、z轴同理
y_ticks = np.logspace(np.log10(y_min), np.log10(y_max), num=5)
ax.set_yticks(y_ticks)
ax.set_yticklabels([f"{tick:.0e}" for tick in y_ticks])

z_ticks = np.logspace(np.log10(z_min), np.log10(z_max), num=5)
ax.set_zticks(z_ticks)
ax.set_zticklabels([f"{tick:.0e}" for tick in z_ticks])

完整修改后的示例代码 #

for it, folder in enumerate(norm_folder_list):
    fig = plt.figure(figsize=(10, 8))
    ax = fig.add_subplot(111, projection='3d')
    
    x_data = folder_data_real[it][0]
    y_data = folder_data_real[it][1]
    z_data = folder_data_real[it][2]
    
    # 过滤非正值数据（可选，根据需求调整）
    mask = (x_data > 0) & (y_data > 0) & (z_data > 0)
    x_data = x_data[mask]
    y_data = y_data[mask]
    z_data = z_data[mask]
    
    ax.scatter(x_data, y_data, z_data)
    
    ax.set_xlabel('Frequency')
    ax.set_ylabel('Time')
    ax.set_zlabel('Permittivity')
    
    # 设置轴范围
    x_min, x_max = min(x_data), max(x_data)
    y_min, y_max = min(y_data), max(y_data)
    z_min, z_max = min(z_data), max(z_data)
    ax.set_xlim(x_min, x_max)
    ax.set_ylim(y_min, y_max)
    ax.set_zlim(z_min, z_max)
    
    # 设置对数刻度
    ax.set_xscale('log')
    ax.set_yscale('log')
    ax.set_zscale('log')
    
    plt.show()

内容的提问来源于stack exchange，提问作者Entangled Being