Makie.jl中surface!函数nan_color参数的正确使用方法

在数据可视化过程中，我们经常需要处理包含NaN(非数值)的数据点。Makie.jl作为Julia生态中强大的可视化工具，提供了nan_color参数来专门控制NaN值的显示颜色。然而，用户在使用surface!函数时可能会遇到nan_color参数似乎不起作用的情况，这实际上是由于对参数使用场景的误解。

问题现象

当用户尝试使用以下代码时，期望将NaN区域显示为红色：

using CairoMakie

lons = -180:180
lats = -90:90
field = [exp(cosd(l)) + 3(y/90) for l in lons, y in lats]

field[1:100,1:100] .= NaN

fig = Figure()
ax = Axis(fig[1,1])
surface!(ax, lons, lats, field; shading = NoShading, nan_color = :red)
fig

结果却发现NaN区域显示为透明而非预期的红色。

原因分析

这个现象的根本原因在于对surface!函数参数的理解偏差。在Makie.jl中：

1. 当NaN值出现在Z坐标(高度值)时，这些点会被视为无效点，无法参与网格生成，因此会直接跳过绘制，表现为透明
2. nan_color参数实际上设计用于控制颜色映射中的NaN值显示，而不是Z坐标中的NaN值

正确使用方法

要实现将特定区域显示为红色的效果，正确的做法是将数据作为颜色参数传递，而不是Z坐标：

surface!(ax, lons, lats, zeros(size(field)); 
    color = field, 
    shading = NoShading, 
    nan_color = :red)

这种用法中：

• Z坐标使用全零平面，确保所有点都有效
• 数据通过color参数传递，此时nan_color参数会正确地将NaN值映射为红色

技术背景

Makie.jl的surface!函数实际上创建的是三维曲面图，它需要：

1. 有效的几何数据(顶点位置)来构建网格
2. 可选的颜色数据来进行着色

当Z坐标包含NaN时，会导致几何数据不完整，无法构建有效的三角网格。而颜色数据中的NaN则可以通过nan_color参数进行特殊处理。

其他注意事项

1. 不同后端表现可能不同：CairoMakie和GLMakie在处理NaN时可能有细微差别
2. 性能考虑：使用全零平面作为Z坐标可能会影响某些光照效果的计算
3. 替代方案：对于复杂场景，可以考虑使用mesh函数，它对NaN的处理方式略有不同

理解这些底层机制可以帮助用户更有效地利用Makie.jl进行科学可视化，特别是在处理不完整或有缺失值的数据集时。