如何为等值线图生成平滑颜色映射并增强高亮区域亮度

我有一个包含数千条数据的数据集，分为X、Y、Z、Z误差四列，部分数据如下：

1.2351e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2417e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2483e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2549e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2615e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2681e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2746e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2812e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.2878e+00  -6.3115e-02 1.8186e+01  1.8186e+01
1.2944e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.3010e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.3075e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.3141e+00  -6.3115e-02 0.0000e+00  0.0000e+00
1.3207e+00  -6.3115e-02 0.0000e+00  0.0000e+00
...

我想生成类似参考图的等值线图，当前使用的代码如下：

import numpy as np  
import matplotlib.pyplot as plt
x, y, z, err = np.genfromtxt(r'dataset_gid.txt', unpack=True)

x1, y1 = np.meshgrid(x, y)
cont1 = plt.tricontourf(x, y, z, cmap='seismic')
plt.colorbar(cont1)
plt.show()

得到的图效果不理想，之后尝试了两种方法：

• 使用scat=plt.scatter(x, y, c=z, s=3);plt.colorbar(scat)，出现异常线条
• 使用plt.tricontourf(x, y, z, levels=50)，效果有所改善，但高亮区域亮度不够，达不到参考图的明亮程度

请问如何调整才能让高亮区域像参考图一样明亮？

解决方案 #

1. 调整颜色映射与显示范围
   如果大部分Z值为0，仅少数区域有高值，默认的全范围色阶会压缩高亮区域的对比度。可以手动限定色阶范围，同时换用更适合突出高亮的色图：

# 获取Z的最大值，聚焦有效数值区间
z_max = np.max(z)
cont1 = plt.tricontourf(x, y, z, levels=50, cmap='hot', vmin=0, vmax=z_max)
plt.colorbar(cont1)
plt.show()

'hot'色图天生偏暖，高亮区域会更明亮；限定vmin和vmax能让高值区域的颜色差异更显著。

2. 自定义高亮度色图
   如果默认色图亮度仍不足，可以手动生成从深色到亮黄/白色的色图，强化高光：

from matplotlib.colors import LinearSegmentedColormap

# 自定义色图：黑→深红→橙→亮黄
colors = [(0,0,0), (0.9,0,0), (1,0.6,0), (1,1,0)]
custom_cmap = LinearSegmentedColormap.from_list('highlight_map', colors, N=50)

z_max = np.max(z)
cont1 = plt.tricontourf(x, y, z, levels=50, cmap=custom_cmap, vmin=0, vmax=z_max)
plt.colorbar(cont1)
plt.show()

3. 叠加散点强化高亮
   如果等值线的高光仍不够突出，可以在等值线图上叠加高Z值的散点，用亮色系标记：

z_max = np.max(z)
cont1 = plt.tricontourf(x, y, z, levels=50, cmap='hot', vmin=0, vmax=z_max)
# 筛选Z值大于阈值的点（比如大于10），用亮黄色标记
high_z_mask = z > 10
plt.scatter(x[high_z_mask], y[high_z_mask], c='yellow', s=4, alpha=0.9)
plt.colorbar(cont1)
plt.show()

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