Makie.jl中hexbin绘图对数坐标轴问题的分析与解决方案

问题背景

在使用Makie.jl数据可视化库时，hexbin函数是一个非常有用的工具，它能够有效地展示二维数据的分布密度。然而，当用户尝试在对数坐标轴下使用hexbin绘图时，会遇到一个技术障碍：直接设置yscale=log10会导致绘图失败，并抛出关于负实数对数计算的错误。

问题现象

当用户尝试在hexbin图中使用对数坐标轴时，例如以下代码：

using Random; Random.seed!(0)
using CairoMakie
n = 10000
xs = randn(n)
ys = randn(n) .|> exp
hexbin(
    xs,
    ys;
    axis = (; yscale=log10), # 这行会导致错误
)

系统会抛出DomainError错误，提示在对数计算中传入了负实数参数。这是因为hexbin内部在处理数据分箱时，可能会生成一些中间值，这些值在对数变换后变为无效的负值或零值。

技术分析

hexbin绘图的核心原理是将二维空间划分为六边形区域，并统计每个区域内数据点的数量。当启用对数坐标轴时，Makie会先对数据进行对数变换，然后再进行分箱计算。然而，这个过程中存在两个潜在问题：

1. 数据预处理顺序：对数变换应该在分箱计算之前完成，但当前实现中顺序可能不正确
2. 边界处理不足：对于接近零的数值，对数变换会产生极大的负值，影响分箱计算

解决方案

目前，推荐使用以下两种解决方案：

方案一：手动对数变换（推荐）

log_ys = log10.(ys)
yticks = floor(Int, minimum(log_ys)) : ceil(Int, maximum(log_ys))
hexbin(
    xs,
    log_ys;
    axis = (;
        yticks = (yticks, [rich("10", y|>string|>superscript) for y in yticks]),
)

这种方法完全避免了hexbin内部的对数变换，而是:

1. 先对数据进行手动对数变换
2. 然后使用变换后的数据绘图
3. 最后通过自定义yticks来显示原始数值的对数刻度

方案二：数据过滤

对于包含零或负值的数据集，可以先进行过滤：

filtered_ys = filter(y -> y > 0, ys)
hexbin(
    xs[ys .> 0],
    filtered_ys;
    axis = (; yscale=log10),
)

技术展望

虽然当前版本(v0.22.2)存在这个限制，但未来Makie.jl可能会在以下方面进行改进：

1. 更智能的坐标变换处理：在分箱计算前自动应用坐标变换
2. 更完善的错误处理：对于无效的变换值提供更有指导性的错误信息
3. 内置对数分箱支持：提供专门针对对数尺度优化的分箱算法

结论

在Makie.jl中实现hexbin图的对数坐标轴显示需要特别注意数据处理顺序。目前最可靠的解决方案是预先对数据进行对数变换，然后手动设置坐标轴刻度。这种方法不仅避免了系统错误，还能提供更精确的可视化效果。随着Makie.jl的持续发展，我们期待未来版本能够原生支持这一功能，简化用户的操作流程。