程序化分形作为栅格化保真度评估的多尺度靶标

《Journal of Imaging》：Procedural Fractals as Multi-Scale Targets for Evaluating Rasterization Fidelity

　　

程序化分形（Procedural Fractals）提供了可复现的多尺度结构，可作为受控靶标用于评估栅格化（Rasterization）过程中引入的与分辨率相关的图像变化。本研究利用十五个程序化生成的PostScript分形图像，分别在150、300、600和1200 DPI下进行栅格化处理，并以1200 DPI作为操作性高分辨率参考，系统考察了栅格化分辨率对图像几何、频谱、高频及结构特性的影响。栅格化保真度通过分形维数偏差（ΔD）、功率谱密度斜率偏差（Δβ）、高频能量比（EHF）以及多尺度平均结构相似度（QMS）四项指标进行量化。结果表明，随着分辨率提高，ΔD与Δβ显著降低（p < 0.001），EHF在四种分辨率下均系统性下降（p < 0.001）。QMS在150 DPI时已达到较高水平（0.9667 ± 0.0188），并在600 DPI时升至0.9931 ± 0.0058（p < 0.001）。在十五个分形图案中，ΔD、Δβ、EHF和QMS分别有12/15、14/15、15/15和15/15个图案表现出预期的单调变化轨迹。上述发现表明，高结构相似度可以与尺度敏感的几何与频谱特性中可测量的分辨率相关差异共存，从而支持程序化分形作为可复现的候选多尺度靶标，用于评估栅格化保真度。

论文解读：程序化分形作为栅格化保真度评估的多尺度靶标

一、研究背景与问题提出

数字图像并非视觉场景的直接表征，而是经过一系列物理与计算过程将场景信息转化为离散像素表示的结果。计算成像将图像形成、采样与计算处理整合于统一框架之中，而栅格化作为将分辨率无关的图形描述转换为离散栅格表示的基本操作，在可视化、数字出版、显示渲染及计算机图形学中具有基础性地位。尽管栅格化常被视为技术性转换步骤，但其本质上涉及空间采样过程。依据经典采样理论，采样密度不足将限制高频空间信息的准确表达，并可能引入混叠及其他与分辨率相关的失真。

传统成像系统评估通常依赖标准化靶标，如ISO 12233:2024所规定的分辨率靶标，此类靶标主要由预设边缘、线条、楔形或周期结构组成，虽能提供成熟的成像性能测量手段，但并未专门用于表征连续分布于宽范围空间尺度上的确定性结构复杂性。分形几何为描述具有尺度依赖复杂性与自相似组织的结构提供了数学框架，程序化生成的分形结构可在多空间尺度上包含几何信息，同时保持确定性与可复现性，因而成为研究采样与分辨率依赖图像形成的受控合成靶标的有力候选。

然而，从数字图像估计分形特性面临特定方法论挑战。分形维数（FD）的估计受有限空间分辨率、可用尺度区间、像素量化及估计器数值实现的影响。已有研究表明采样分辨率与量化可显著影响数字图像FD估计结果，降采样像素分辨率亦会影响数学定义分形与实验观测分形的广义维数估计。此外，盒计数估计中盒尺寸与空间尺度范围的恰当选择同样至关重要。频率域分析方面，功率谱密度（PSD）分析可表征图像能量在空间频率上的分布，而采样间隔与空间分辨率会影响PSD派生描述符的敏感性。结构图像质量度量方面，结构相似度指数（SSIM）及其多尺度扩展提供了结构层面的质量评估视角。

尽管上述领域均有大量独立研究，但以往工作主要分别考察分辨率对分形性质测量的影响、采样对频谱描述符的影响或图像处理操作对结构质量的影响，鲜有研究从互补视角出发，在统一实验框架内综合运用几何、频谱、高频能量与结构度量，以程序化定义的分形结构作为受控多尺度靶标，定量表征栅格化过程本身。本研究正是针对这一空白，系统考察栅格化分辨率对图像多域特性的影响。

二、研究设计与主要技术方法

研究人员构建了一个包含十五个程序化生成分形图像的实验数据集，每个分形均由独立的PostScript程序定义，描述确定性的生成过程，包含迭代几何变换、随机坐标生成及重复路径构建操作。所有分形通过固定伪随机初始化保证完全可复现。研究人员使用GSview 5.0的File → Convert界面，通过Ghostscript 10.06.0的png16m输出设备，将各分形程序分别在150、300、600和1200 DPI四种空间采样密度下进行栅格化。输出图像尺寸分别为1275×1650、2550×3300、5100×6600和10200×13200像素，对应1:2:4:8的采样关系，保持相同物理视场。1200 DPI栅格化作为操作性高分辨率参考。

图像预处理方面，RGB图像经亮度变换转换为灰度，并线性归一化至[0,1]范围，未施加滤波、去噪、直方图均衡或锐化等增强操作。分形维数分析中，使用OpenCV Canny边缘检测器（固定阈值50和150）提取二值边缘图，再以经典盒计数法估计FD，回归拟合在双对数域进行。PSD斜率分析中，对去均值灰度图像计算二维离散傅里叶变换，经径向平均后，在双对数坐标下线性回归估计频谱标度指数β。高频能量比（E_HF）定义为预定义高频区域（径向频率阈值0.35 r_max）内光谱功率占总光谱功率的比例。多尺度平均结构相似度（Q_MS）为四个空间尺度（s∈{1,2,4,8}）下常规SSIM值的未加权算术平均值，计算前将1200 DPI参考图像经面积插值下采样至对应测试图像尺寸以匹配空间维度。统计分析采用Friedman检验（重复测量非参数检验）、Kendall一致性系数（W）、Holm校正的Wilcoxon符号秩检验及描述性最小二乘线性回归分析。

三、研究结果

3.1 高频频谱保持

高频能量比（E_HF）随栅格化分辨率提高而系统性下降，均值从150 DPI时的0.0355 ± 0.0175降至300 DPI时的0.0226 ± 0.0133、600 DPI时的0.0141 ± 0.0106及1200 DPI时的0.0084 ± 0.0079。全部十五个分形均呈现非递增轨迹（15/15，100%）。Friedman检验确认分辨率对E_HF具有统计学显著影响（p < 0.001），Kendall W = 1.00，表明十五个分形在分辨率依赖排序上完全一致。

3.2 分形几何保持

分形维数偏差（ΔD）随分辨率提高而总体下降，均值从150 DPI时的0.0647 ± 0.0365降至300 DPI时的0.0529 ± 0.0226和600 DPI时的0.0237 ± 0.0100。十二个分形（12/15，80%）在三个非参考分辨率上呈现非递增轨迹。Friedman检验显示分辨率对ΔD具有统计学显著影响（p < 0.001），Kendall W = 0.72，表明存在实质性一致性。

3.3 频谱斜率保持

PSD斜率偏差（Δβ）随分辨率提高而逐步下降，均值从150 DPI时的0.2608 ± 0.0926降至300 DPI时的0.1688 ± 0.0565和600 DPI时的0.0820 ± 0.0272。十四个分形（14/15，93.3%）呈现非递增轨迹。Friedman检验确认分辨率对Δβ具有统计学显著影响（p < 0.001），Kendall W = 0.81，表明强一致性。

3.4 结构相似度

多尺度平均结构相似度（Q_MS）随分辨率提高而系统性增加，均值从150 DPI时的0.9667 ± 0.0188升至300 DPI时的0.9813 ± 0.0131和600 DPI时的0.9931 ± 0.0058。全部十五个分形（15/15，100%）在三个非参考分辨率上呈现非递减轨迹。Friedman检验确认分辨率对Q_MS具有统计学显著影响（p < 0.001），Kendall W = 1.00，表明完全一致性。

3.5 整体分辨率依赖行为

四项描述符均表现出统计学显著的与分辨率相关的差异，但幅度与一致性因度量而异。ΔD与Δβ随分辨率提高而降低，E_HF在四种分辨率下持续下降，Q_MS则向操作性1200 DPI参考值1趋近。值得注意的是，Q_MS在150 DPI时已接近1，而同一分辨率下ΔD与Δβ的相对变化幅度明显更大，表明高结构相似度可与尺度敏感几何与频谱特性的可测量差异共存。

3.6 事后两两比较

Holm校正后的Wilcoxon符号秩检验显示，ΔD在150对300 DPI（p_Holm = 0.0151）、150对600 DPI（p_Holm < 0.001）及300对600 DPI（p_Holm < 0.001）之间均存在显著差异。Δβ的三个非参考分辨率两两比较均显著（p_Holm < 0.001）。E_HF的六个两两比较（四种分辨率全部配对）均显著（p_Holm < 0.001）。Q_MS的三个非参考分辨率两两比较均显著（p_Holm < 0.001）。所有纳入相应推断分析的两两比较在Holm校正后均保持统计学显著。

3.7 描述性线性趋势分析

以log2(DPI)表示分辨率时，描述性最小二乘回归显示，ΔD每翻倍分辨率约下降0.0205（R² = 0.9438），Δβ每翻倍约下降0.0894（R² = 0.9997），Q_MS每翻倍约上升0.0132（R² = 0.9963），E_HF拟合斜率为每翻倍?0.0090（R² = 0.9693）。四项趋势方向均与重复测量及单调性分析结果一致。

四、讨论与结论

讨论部分指出，四项描述符表征了栅格化保真度的不同方面，不应视为可互换的单一误差度量。分形维数偏差表征尺度依赖几何复杂性的变化，PSD斜率偏差反映频谱标度变化，E_HF量化预定义高频傅里叶区域的相对贡献，Q_MS则在四个空间尺度上度量与操作性高分辨率参考的结构一致性。各描述符的不同响应具有信息价值而非矛盾：高Q_MS值与几何及频谱描述符的较大相对变化形成对比，表明栅格化可保持整体结构组织，同时改变在跨尺度或频域分析中才显现的属性。

程序化分形作为多尺度栅格化靶标具有独特优势：底层场景描述确定且可复现，递归或尺度依赖组织提供多空间尺度结构内容，程序化生成允许在系统变化栅格化分辨率的同时保持底层描述不变。本研究的受控设计有效分离了栅格化相关变化与场景内容变化。研究结果支持程序化分形作为可复现的候选多尺度测试靶标，而非普遍标准化的栅格化靶标。

研究局限性包括：1200 DPI栅格化作为操作性高分辨率参考而非绝对真值；实验基于单一Ghostscript/png16m栅格化管线，不同栅格化引擎、抗锯齿实现或渲染设置可能产生不同定量响应；数据集仅包含十五个程序化分形，未涵盖程序化几何、矢量图形或自然图像内容的全部多样性；E_HF并非混叠的直接估计量，Q_MS为自定义多尺度度量而非规范加权MS-SSIM。

结论部分指出，提高栅格化分辨率与分形维数和PSD斜率相对操作性1200 DPI参考的偏差逐步减小相关，同时伴随相对高频能量比的系统性下降。多尺度平均结构相似度在150 DPI时已较高并随分辨率进一步提高。重要的是，高结构相似度可与尺度敏感几何与频谱特性的可测量差异共存，因此当栅格化图像用于定量图像分析时，结构相似度单独可能不足以完整表征栅格化保真度。描述性趋势分析显示这些分辨率依赖变化在独立评估范围内随分辨率逐次翻倍近似线性。研究结果支持程序化分形作为受控评估栅格化与采样效应的候选多尺度测试靶标，其确定性源描述与固有结构多尺度性可实现可复现比较，互补定量描述符可揭示分辨率依赖图像形成的不同方面。未来研究应扩展至更多程序化结构、栅格化引擎、抗锯齿与重建策略及更高或独立定义的参考表示。