Gnuplot循环生成2D直方图分箱问题求助

核心问题梳理 #

你在循环处理多.dat文件时遇到的语法错误和数据缺失问题，根源在于循环中无法直接复用全局临时数据块$Data/$Bins，且部分命令格式不符合gnuplot 6.0的语法要求：

• stats命令格式错误：替换$Data为单个文件时，无需额外添加[*:*]，正确格式应为stats infile u 1:2 nooutput
• smooth freq生成分箱数据时，不能直接在plot循环中引用原文件列，需先将当前文件数据读入临时数据块再处理
• 转义\$1的问题：循环迭代时引用数据列需转义，但结合临时数据块可规避该问题

修正后的循环代码示例 #

假设你已定义XYtoBinNo()、BinNoToX()、BinNoToY()三个自定义分箱函数，以及分箱参数dx（x方向箱宽）、dy（y方向箱宽），以下是完整的循环处理代码：

# 定义分箱相关函数（基于@theozh原回答调整）
XYtoBinNo(x,y) = floor(x/dx) + floor(y/dy)*nbinsx
BinNoToX(n) = (n % nbinsx)*dx
BinNoToY(n) = floor(n/nbinsx)*dy

# 设置分箱参数（根据你的数据范围调整）
dx = 0.5
dy = 0.5

# 循环处理所有.dat文件
FILES = system("ls *.dat")
do for [infile in FILES] {
    # 1. 将当前文件数据读入临时数据块$Data
    set table $Data
    plot infile u 1:2:3 with table
    unset table

    # 2. 统计数据范围，计算x方向分箱数
    stats $Data u 1:2 nooutput
    nbinsx = ceil((STATS_max_x - STATS_min_x)/dx)

    # 3. 生成分箱统计数据到临时数据块$Bins
    set table $Bins
    plot $Data u (XYtoBinNo($1,$2)):3 smooth freq
    unset table

    # 4. 绘制分箱图并保存
    set output infile[:-4].".png"  # 输出为对应文件名的PNG
    plot $Bins u (BinNoToX($1)+dx/2):(BinNoToY($1)+dy/2):(dx):(dy):2 w boxxy fc palette z notitle
    unset output
}

关键修正点说明 #

• 循环内为每个文件单独生成$Data和$Bins临时数据块，避免跨文件数据污染
• stats命令直接作用于$Data，解决undefined variable: u错误
• 分箱统计用$Data而非原文件，避免column() called from invalid context错误
• 绘图时将箱位置调整为BinNoToX($1)+dx/2，确保箱中心对应分箱区间中点，避免因位置偏移导致的"数据缺失"假象
• 无需手动转义\$1，临时数据块内的列引用不会被循环迭代变量干扰

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