锐度:它是什么以及如何测量

锐度:它是什么以及如何测量

 

测量
锐度
调制传递
函数
空间 频率单位_汇总
指标
MTF 测量
矩阵
斜边

 

 

测量锐度

清晰度决定了成像系统可以再现的细节量。它由不同色调或颜色的区域之间的边界定义。

 

图 1. 条形图案:原始(图的上半部分)镜片退化(图的下半部分)

图 1. 条形图案:原始(图的上半部分)镜片退化(图的下半部分)

图 2. 图像边缘的清晰度示例

图 2. 来自MTF 曲线和图像外观的图像边缘的清晰度示例

 

在图 1 中,锐度用空间频率增加的条形图案来说明。图形的顶部清晰,边界清晰;下半部分是模糊的,并说明了条形图案在通过模拟镜头后是如何退化的。

注意所有镜头都会在一定程度上模糊图像。

清晰度在图像边缘(图 2)等特征上最为明显,可以通过边缘(阶跃)响应来测量。

有几种方法用于测量清晰度,包括 10-90% 上升距离技术、调制传递函数 (MTF)、特殊域和频域以及斜边算法。

上升距离和频域

图像清晰度可以通过图像内边缘的“上升距离”来衡量。

使用这种技术,清晰度可以通过像素级在其最终值的 10% 到 90% 之间的距离来确定(也称为10-90% 上升距离;参见图 3)。 

 

图 3. 模糊和锐利边缘 10-90% 上升距离的图示

图 3. 模糊和锐利边缘 10-90% 上升距离的图示

 

上升距离没有被广泛使用,因为没有方便的方法可以根据其各个组件(即镜头、数字传感器和软件锐化)的上升距离来计算成像系统的上升距离。 

为了克服这个问题,在频域中 进行测量,其中频率以每距离(毫米、英寸、像素、图像高度或有时角度 [度或毫弧度])为单位的周期或线对来测量。

在频域中,可以通过组合由以周期或频率表征的纯音(正弦波)组成的信号来创建复杂信号(音频或图像)(图 4)。一个完整系统的响应是每个组件响应的乘积。

 

图 4. 在频域测量清晰度

图 4. 周期定义(1/频率)

 

频域和空间域通过傅里叶变换相关联。 

 

F=− Fe-ω td

F=12π _− Fω )eωt _dω

在哪里

f = 频率 = 1/周期   (较短的周期对应较高的频率);t = 时间;  ω = 2πf 

 

 

高频(短周期)的系统响应越大,系统可以传达的细节就越多(图 5)。系统响应可以通过频率响应曲线F(f)来表征。

注意高频对应精细的细节。

调制传递函数

 

图 5. 高频对应空间和频域的精细

图 5. 高频对应空间和频域的精细细节

 

给定空间频率下的相对对比度(输出对比度/输入对比度)称为调制传递函数( MTF 类似于空间频率响应(SFR),是衡量清晰度的关键。在图 6 中,MTF 用正弦和条形图案、幅度图和对比图来说明——每一个都具有从左到右连续增加的空间频率。 

注意Imatest 可互换使用 SFR 和 MTF。SFR 更常与完整的系统响应相关联,而 MTF 通常与特定组件的单个效应相关联。换言之,系统SFR相当于成像系统中各部件的MTF的乘积。

高空间频率(右侧)对应于精细的图像细节。摄影组件(胶片、镜头、扫描仪等)的响应往往在高空间频率下“滚降”。这些组件可以被认为是通过低频并衰减高频的低通滤波器。

 

图 6. 正弦和条形图、幅度图和对比度 (MTF) 图

图 6. 正弦和条形图、幅度图和对比度 (MTF) 图

 

图 6 由上图、中图和下图组成,描述如下:

  • 条形图案、正弦图案(上图)- 上图说明 (1) 原始正弦图案、(2) 带有镜头模糊的正弦图案、(3) 原始条形图案和 (4) 带有镜头模糊的条形图案. 请注意,镜头模糊会导致对比度在高空间频率下下降。
  • 幅度(中间图)—中间图显示了带有镜头模糊的条形图案的亮度(MTF 方程部分中的“调制” V )(参见图 6 中的红色曲线)。由纯频率组成的正弦模式的调制用于计算 MTF。(请注意,在高空间频率下对比度会降低。)
  • MTF %(下图)- 下图显示了相应的正弦模式对比度(见 蓝色曲线;代表 MTF),这也在MTF 方程部分中定义。根据定义,低频 MTF 限制始终为 1 (100%)。在图 6 中,MTF 在 61 lp/mm(每毫米线对)时为 50%,在 183 lp/mm 时为 10%。请注意,频率和 MTF 都以指数符号 [10 0 = 1%; 10 1 = 10%;10 2 = 100% 等];幅度(中间图)以线性刻度显示。完整成像系统的 MTF 是其各个组件的 MTF 的乘积。

MTF方程

MTF 的方程是从 空间频率 f处的正弦图案对比度C(f)推导出来的,其中

Cf=× _-_× _+_     用于亮度(“调制”)  V

Ff100 ×Cf)C)      注意:这会将低空间频率下的 MTF 归一化为 100%。

为了在低空间频率下正确归一化 MTF,测试图必须具有一些低频能量。这是由倾斜边缘中的大亮和暗区域以及 Imatest 使用的大多数图案中的特征提供的,但在线条和网格中不存在。对于可以控制锐化的系统,推荐的主要MTF 计算是倾斜边缘,它是从脉冲响应(即对窄线的响应)的傅里叶变换计算得出的,它是导数( d/dxd /dy ) 的边缘响应。幸运的是,您无需了解傅立叶变换即可了解 MTF。

传统分辨率测量

传统的分辨率测量包括观察条形图案的图像,最常见的是 USAF 1951 图表(图 7),并估计条形图案明显不同的最高空间频率 (lp/mm)。这种观察(也称为消失分辨率)对应于大约 10-20% 的 MTF。因为消失的分辨率是图像信息消失的空间频率——在它 不 可见的地方,它强烈依赖于观察者的偏差,并且是图像清晰度的不良指标。 (这是****中没有 Woozle 的地方。) 

注意:USAF 1951 图表(很久以来被空军放弃)不太适合计算机分析,因为它使用空间效率低,并且它的条形三连音缺乏低频参考。此外,图表位置(采样阶段)的微*变化会导致其条形的外观发生变化,因为它们从与像素阵列同相变为异相。这会对消失分辨率估计产生不利影响。

 

图 7. USAF 1951 图表; Imatest 不支持

图 7. USAF 1951 图表;Imatest 不支持

 

更好的图像清晰度指标是空间频率,其中 MTF 是其低频值 (MTF50) 的 50% 或峰值 (MTF50P) 的 50%。推荐使用 MTF50 和 MTF50P 来比较不同相机和镜头的锐度,因为 

  1. 图像对比度是其低频或峰值的一半,因此细节仍然很明显。(眼睛对 MTF 为 10% 或更低的空间频率的细节不敏感。)
  2. 大多数相机的响应在 MTF50 和 MTF50P 附*迅速下降。MTF50P 是强烈锐化相机的更好指标(在我们的电子成像 2020 论文中进行了解释)。
  3. 50%的水*与图像的信息容量有关

注意下面的摘要指标讨论了其他清晰度指标。 

尽管可以直接从正弦模式的图像(使用 Rescharts、Log Frequency、  Log F-Contrast和 Star Chart)估计 MTF,但ISO 12233 斜边技术提供了更准确和可重复的结果,并且更有效地利用了空间。 倾斜边缘图像可以通过以下MTF 测量矩阵列出的模块之一进行分析。

空间频率单位

 

图 8. 在 SFR 和 Recharts 模块(SFRplus、eSFR ISO、Star 等)的设置或更多设置窗口中选择空间频率单位。

图 8. 在 SFR 和 Recharts 模块(SFRplus、eSFR ISO、Star 等)的设置或更多设置窗口中选择空间频率单位。

 

大多数读者都熟悉时间频率例如,声音的频率(以周期/秒或赫兹为单位)与其感知的音高密切相关。无线电传输的频率(以千赫、兆赫和千兆赫为单位)也很熟悉。

空间频率以每距离的周期(或线对)而不是时间来衡量。与时间(例如,音频)频率响应一样,响应越扩展,可以传达的细节就越多。

空间频率单位可以SFR 和 Rescharts 模块(SFRpluseSFR ISOStar等图 8)的“设置”或“更多设置”窗口中选择,并且是用于确定相机可以再现多少细节或像素的好坏的测量被利用。

过去的胶片相机镜头测试使用每毫米线对 (lp/mm),这对于比较镜头非常有效,因为大多数 35mm 胶片相机具有相同的 24 x 36mm 图像尺寸。但数字传感器尺寸差异很大——从可拍照手机的对角线不到 5 毫米到全画幅相机的 43 毫米对角线,再到中画幅的更大对角线。出于这个原因,建议使用每图片高度的线宽 (LW/PH) 来测量相机可以再现的总细节。请注意,LW/PH 等于 2 × lp/mm ×(图片高度,单位为 mm)。

另一个有用的空间频率单位是每像素周期 (C/P),它表示单个像素的利用程度。单位的选择还受图像(传感器)或物体上的性能是否具有首要重要性的影响:请参阅比较不同相机的清晰度无需使用数码相机的实际距离(毫米或英寸),尽管可以使用此类测量值(表 1)。

表 1. 空间频率单位汇总,以及在选定频率单位中引用 MTF 的方程。

机动特遣队单位 应用 方程

周期/像素 (C/P)

显示像素的使用情况。奈奎斯特频率Nyq始终为 0.5 C/P。   

循环/距离

(循环/毫米或循环/英寸)

传感器上每距离的循环次数必须输入像素间距或间距。用于比较过去标准胶片格式的分辨率(例如,24x36mm 用于 35mm 胶片)很流行。 FC/ P)像素间距

线宽/图片高度 (LW/PH)

注意:对于裁剪的图像,在更多设置尺寸输入中输入原始图片高度

测量整体图像清晰度。这是比较具有不同传感器尺寸和像素数的相机性能的最佳单位。线宽是电视测量的传统方法。推荐用于比较不同相机的锐度中以图像为中心的应用。
请注意,1 个周期 = 1 个线对 (LP) = 2 个线宽 (LW)。
×F(LP _H) ;
×F(C) ×PH

线对/图片高度 (LP/PH)

注意:对于裁剪的图像,在更多设置尺寸输入中输入原始图片高度

测量整体图像清晰度。与 LW/PH 相差 2 倍。由 dpreview.com 使用。 F(长宽_H) /2 ;
F(C) ×PH

角频率

周期/毫弧度

角频率。必须输入像素间距或间距。以毫米为单位的焦距 ( FL ) 通常包含在商业图像文件的 EXIF 数据中。如果它不可用,则必须手动输入,通常在设置窗口底部的 EXIF 参数区域中。如果缺少像素间距或焦距,单位将默认为 Cycles/Pixel。
Cycles/degree 可用于将相机系统与人眼进行比较,人眼的 MTF50 约为 20 Cycles/Degree(取决于个人的视力和照明度)。

FL可以从简单的镜头方程*计算,其中1是镜头到图表的距离(易于测量),2是镜头到传感器的距离, 放大倍率//s1/s2=s2/s1

F大号=s1/ (1+1 / |)=  s2/ (1+ |)

*除非1  >> 2(乘以 100 倍或更多),否则镜头几何形状(12FL对于计算M是可靠的, 因为镜头可能会 明显偏离简单的镜头方程

 

0.001 × MF(循环毫米) ×F毫米

 

周期/度

 

π180×F(循环毫米) ×F毫米

循环/对象 mm
循环/对象 in

被拍摄物体上每距离的循环次数(有些人认为是“主题”)。必须输入像素间距和放大倍数,但有一个重要例外*应在系统规范引用被摄对象时使用(例如,如果需要检测一定宽度的特征)。推荐用于比较不同相机的锐度中的以对象为中心的应用。

*对于输入条间距的SFRplus 、输入配准标记垂直间距时的eSFR ISO或输入正方形长度时的棋盘格,不使用像素间距或输入放大倍数直接计算每对象距离的周期数,即从几何计算。在 Imatest 2021.2 之前,您必须在像素间距字段中输入一个数字,但该数字不用于实际计算。对于造成的混乱,我们深表歉意。

F(循环距离) ×放大倍率|

Cycles/distance在图像传感器上是 Cycles/mm 或 Cycles/in 。

线宽/裁剪高度
线对/裁剪高度

主要用于在活动图表高度(而不是总图像高度)显着时进行测试。不再推荐,因为它取决于未标准化的作物大*。  

线宽/特征 Ht(Px)
线对/特征 Ht(Px)

(以前的线宽或线对/N 像素 (PH))

When either of these is selected, a Feature Ht pixels box appears to the right of the MTF plot units (sometimes used for Magnification) that lets you enter a feature height in pixels, which could be the height of a monitor under test, a test图表或具有非活动区域的图像中的活动视野。此单位选择对于比较具有不同图像或像素大*的相机的特定对象的分辨率很有用。推荐用于比较不同相机的锐度中的以对象为中心的应用。

×F(C) ×特征高度

F(C) ×特征高度

PH = 图片高度(以像素为单位)。FL (mm) = 镜头焦距,单位为 mm。  像素间距= 每像素距离 = 1/(每距离像素数)。  
注意:不同的单位随着图像传感器和像素大*的不同而不同。

 

比较不同相机的清晰度,建议基于两种广泛应用类型之一的空间频率单位:

    • 以图像为中心(例如风景摄影,其中图像传感器的细节很重要):建议使用每图片高度的线宽(或线对)。
    • 以对象为中心(用于医疗、机器视觉等,其中对象的细节很重要):推荐每个特征高度的周期/对象距离或 LW(或 LP)。

摘要指标

从 MTF 曲线中得出几个汇总指标来表征整体性能。这些指标用于许多显示,包括 SFR/SFRplus/eSFR ISO 边缘/MTF 图(请参阅Imatest 倾斜边缘结果)和 SFRplus 3D 地图中的二次读数。

摘要指标 描述 评论
MTF50
MTF nn
MTF 为低 (0) 频率 MTF的 50% ( nn %) 的空间频率。MTF50 ( nn  = 50) 被广泛使用,因为它对应于电气工程中的带宽(半功率频率)。 最常见的汇总指标;与感知锐度密切相关。随着软件锐化的增加而增加;可能会产生误导,因为它“奖励”过度锐化,这会导致边缘可见且可能令人讨厌的“光晕”。
MTF50P
MTF nn P
MTF 为峰值 MTF的 50% ( nn %) 的空间频率。 对于中低软件锐化,与 MTF50 相同,但在存在软件锐化峰值(最大 MTF > 1)时低于 MTF50。与 MTF50 相比,对软件锐化的敏感度要低得多(如我们在 Electronic Imaging 2020 上发表的一篇论文所示)。总而言之,一个更好的指标。
MTF面积
归一化
MTF 曲线下的面积(低于奈奎斯特频率),归一化为其峰值(  当锐化很少或没有锐化时, f = 0 时为 1,但对于强锐化,峰值可能为 » 1)。 一个特别有趣的新指标,因为它密切跟踪 MTF50 很少或没有锐化,但不会增加强烈的过度锐化;即,它不会奖励过度锐化。还是比较陌生的。倾斜边缘 MTF 测量一致性中进行了描述。
MTF10、MTF10P、
MTF20、MTF20P
MTF 是零频率或峰值 MTF 的 10% 或 20% 的空间频率 这些数字很有趣,因为它们与“消失的分辨率”(瑞利极限)相当。噪音会严重影响 10% 或更低水*的结果。MTF20(或 MTF20P)在 Line Widths per Picture Height (LW/PH) 中最接*模拟TV Lines有关测量监视器电视线路的详细信息,请参见此处

MTF 测量矩阵——比较不同的图表和测量技术

Imatest有许多测量 MTF 的模式——斜边、对数频率、对数 f 对比度、西门子星、死叶(溢出硬币)、随机 1/f 和双曲线楔形——每种模式在消费类相机中往往给出不同的结果,其中大多数具有不均匀的图像处理——通常是双边滤波——这取决于局部场景内容。锐化(高频提升)趋向于在对比度特征附*最大(在较高对比度边缘附*较大),而降噪(高频削减,可能会掩盖精细纹理)在没有它们的情况下往往最大。出于这个原因,MTF 测量值在不同的测试图上可能会有很大的不同。

原则上,当没有应用非均匀或非线性图像处理(双边滤波)时,MTF 测量值应该是相同的,例如,当使用dcraw 或 LibRaw对图像进行去马赛克且没有锐化和降噪时。但这并不完全发生,因为在所有使用滤色器阵列 (CFA) 的相机中都存在去马赛克,它涉及一些非线性处理。下图总结了不同图案对图像处理的敏感性。

图像处理(双边滤波)对 MTF 测量的影响比较:
倾斜边缘和楔形最容易被锐化。
随机 1/f 模式具有最少的锐化和最多的降噪。

下面的 MTF 矩阵表列出了 Imatest 测量 MTF 方法的属性、优点和缺点。

 

MTF 测量矩阵
请注意,这些模式中的大多数都受益于*均多个(相同配准的)图像以减少噪声的影响。
测量
模式
优点/缺点/敏感性 主要用途和评论
斜边
(SFR
SFRplus
eSFRISO
SFRreg
Checkerboard)
最有效地利用空间:可以创建 MTF 响应的详细地图。SFRpluseSFR ISOSFRregCheckerboard
中的快速自动区域检测 快速计算。 对噪音相对不敏感(如果应用降噪,则更免疫 符合 ISO 12233 标准,使用“分级”(超分辨率)算法,允许在奈奎斯特频率 (0.5 C/P) 以上测量 MTF制造测试的最佳模式。




在具有强锐化和降噪的系统中可能会产生乐观的结果(即,它可能会被信号处理所欺骗,尤其是具有高对比度 (≥ 10:1) 边缘的系统。在具有极端混叠(奈奎斯特频率以上的强能量)
的系统中给出不一致的结果),尤其是*区域。
对锐化最敏感,尤其是对高对比度 (≥10:1) 边缘;
对软件降噪最不敏感。

这是Imatest中的主要 MTF 测量值。

镜头和相机测试的最有效模式,尤其是在需要 MTF 响应图的情况下。

旧 ISO 12233:2000 标准中推荐的高对比度 (≥40:1) 会产生不可靠的结果(剪辑、伽马问题、双边滤镜过度锐化)。新的 ISO 12233:2014 标准推荐 4:1 对比度。这是我们对所有新工作的建议(使用SFRpluseSFR ISO)。

与 Siemens Star 相比, Slanted-edge 与 Siemens Star相比毫不逊色

日志频率 根据第一原理计算。显示彩色波纹。
对噪音敏感。空间利用效率低。
主要用作对其他方法的检查,这些方法不是根据第一原理计算的。
记录 f 对比度 说明非均匀图像处理效果的最佳模式。
对噪音敏感。
对图像(高对比度)顶部附*的锐化和(低对比度)底部附*的降噪具有很强的敏感性,两者之间有渐变。
对锐化/降噪的敏感性是此图表的一个优势,旨在说明信号处理如何随图像内容(特征对比度)而变化。显示由于软件降噪导致的精细细节损失。
西门子之星 包含在 ISO 12233:2014 标准中。对噪音比较不敏感。提供定向 MTF 信息。
空间利用缓慢、效率低下。外半径处有限的低频信息使 MTF 归一化变得困难。
对锐化和降噪的中等敏感度。
由Image Engineering推广用于一般测试,但空间细节仅限于 3×3 或 4×3 网格。Slanted-edge vs Siemens Star中的slanted-edge 相比。
枯叶
(洒落的硬币)
测量纹理模糊/锐度/锐度。模式统计类似于典型图像。
空间利用效率低。一个棘手的噪声功率减法算法*可以降低对噪声的非常高的敏感度,但多幅相同图像的信号*均效果更好。
对锐化的中等敏感度和对降噪的强敏感度使其可用于与主观观察密切相关的整体纹理锐度度量。
由堆叠的随机大*的圆圈组成。强烈的行业兴趣,特别是来自拍照手机图像质量 (CPIQ) 组。

 

死叶(溢出的硬币)和随机图表都使用随机(死叶)模块进行分析。强大的双边过滤会导致误导性的结果

随机(尺度不变) 揭示精细细节(纹理)的渲染效果:系统对软件降噪的响应。
对锐化
最不敏感,对软件降噪最敏感
衡量相机呈现精细细节(纹理)的能力,即低对比度、高空间频率的图像内容。*可以使用与图案相邻的灰色斑块从Imatest中的测量中去除噪声功率。
利用 ISO 12233:2000 或eSFR ISO图表上的楔形图案。
MTF 在奈奎斯特和半奈奎斯特频率附*并不准确(它对采样相位变化非常敏感)。不适合作为主要的 MTF 测量。
对锐化敏感。
对噪音敏感。空间利用效率低。
测量“消失的分辨率”:在楔形图案中线条开始消失的地方,通常在 ISO 12233 图表中,其中需要三个区域(包括用于低频参考的方形区域)才能获得合理的 MTF 测量值(准确度低于其他方法由于采样相位灵敏度)。使用eSFR ISO更方便
通过采集和信号*均多个图像,可以大大降低噪声(和低信噪比 - SNR)的影响。

 

 

空间频率响应的斜边测量

几个 Imatest 模块使用斜边技术测量 MTF,包括:

斜边分析的主要输出是Edge/MTF 图,可以通过点击下面的按钮查看。许多其他结果可用,包括摘要和 3D 图,显示横向色差和其他结果以及 MTF。

 

Edge/MTF 图:Imatest 的主要斜边结果

Imatest SFR的边缘/MTF 图(用于SFRplus图表图像)显示在右侧。SFRpluseSFR ISOSFRregCheckerboard产生类似的结果,甚至更多

 

(左上)  一幅显示*均边缘色调的窄幅图像。它与正下方的*均边缘轮廓(空间域)图对齐。

(左中) *均边缘(空间域):此处显示的*均边缘轮廓线性化(默认)。一个关键结果是边缘上升距离 (10-90%),以像素和每个图片高度的上升距离数显示。其他参数包括过冲和下冲(如果适用)。该图可以选择显示线扩散函数(LSF:边缘的导数)。

(左下)MTF(频域):空间频率响应 (MTF),显示为奈奎斯特频率的两倍。关键的总结结果包括MTF50,对比度下降到其低频值的 50% 的频率,以及MTF50P,对比度下降到其峰值 的 50% 的频率,这与感知的图像清晰度很好地对应。单位是每像素周期 (C/P) 和每图片高度线宽 (LW/PH)。其他结果包括Nyquist 的 MTF(0.5 个周期/像素;采样率/2),这表明混叠和用户选择的二次读数的可能严重程度,以及二次读数。奈奎斯特频率显示为一条垂直的蓝线。当输入(手动)像素间距和输入镜头焦距(通常来自 EXIF 数据,但可以手动输入)时,衍射限制的 MTF 响应显示为浅棕色虚线。

该图像被强烈(但不过分)锐化。

SFR 结果:MTF(锐度)图更详细地描述了该图。

MTF 曲线和图像外观包含几个示例,说明 MTF 曲线和感知锐度之间的相关性。

 

 

 

为什么边缘是倾斜的?

纯垂直或水*边缘的 MTF 结果高度依赖于采样相位(边缘和像素位置之间的关系),因此根据精确的(子像素)边缘位置,每次运行可能会有所不同。边缘是倾斜的,因此 MTF 是根据许多采样阶段的*均值计算的,这使得结果更加稳定和稳健(图 9 )。

什么边角效果最好?

在可能的情况下,边缘角度应与最*的垂直 (V)、水* (H) 或 45 度方向相差大于 ±2 度。原因是垂直、水*和45°边缘的结果对边缘和像素之间的关系非常敏感(即它们是相位敏感的)。将边缘倾斜 2 或 3 度以上可避免此问题。

ISO 12233 标准建议角度为 5 或 5.71 度 (arctan(0.1))。这个角度不是神圣的——MTF 并不强烈依赖于边缘角度。3 到 7 度的角度可以正常工作。对于相对于最*的 V 或 H 方向的非零边角θ应用余弦校正,如右图所示。θ大于约 8 度 (cos(8º) = 0.99) 时,校正是显着的。根据区域选择,初始 MTF 和相应的频率f是从垂直或水*线(以蓝色显示)计算的。真正的 MTF 被定义为垂直于边缘——沿着红色 线。由于 沿 红线(垂直于边缘)的实际过渡长度比 沿蓝色(V 或 H)线的测量过渡短,并且由于用于测量 MTF 的频率f与实际过渡长度成反比,

F=F_/ cos(θ)

具有频率单位的相应摘要度量MTFnnMTF50MTF50P等)在初始值之上增加。

FMF(初始/ c θ )

大边角的主要缺点是可用区域面积可能会减少,特别是对于SFRreg图案。

 

图 9. 裁剪的高对比度垂直边缘(结果无效)

图 9. 裁剪的高对比度垂直边缘(结果无效)

 

 

边缘对比度和剪裁

边缘对比度最多应限制在 10:1,一般建议使用 4:1 的边缘对比度。原因是高对比度边缘(>10:1,例如在旧的 ISO 12233:2000 图表中发现的)可能会导致饱和或削波,从而导致边缘带有尖角,夸大 MTF 测量值。有关更多详细信息,请参阅使用 Rescharts 斜边模块,第 2 部分:警告 - 剪辑。 

 

斜边的优缺点

好处
  • 最有效地利用空间,这使得创建 MTF 响应的详细地图成为可能
  • SFRpluseSFR ISOSFRregCheckerboard中的快速、自动化区域检测
  • 快速计算
  • 对噪音相对不敏感(如果应用降噪,则高度免疫)
  • 符合 ISO 12233 标准,其“分箱”(超分辨率)算法允许在奈奎斯特频率 (0.5 C/P) 以上测量 MTF
  • 制造测试的最佳模式
缺点
  • 在具有强图像相关锐化的系统(即,锐化量随边缘对比度增加)的系统中可能会产生乐观的结果。这种类型的图像处理(双边过滤)在消费类相机中几乎是普遍的。
  • 在具有极端混叠(高于奈奎斯特频率的强能量)的系统中给出不一致的结果,尤其是在*区域中。
  • 不适合测量精细纹理,推荐使用对数频率对比度溢出硬币(枯叶)图案。

注意:Imatest Master 可以计算几乎任何角度的边缘的 MTF,但由于采样相位敏感,应避免精确的垂直、水*和 45°。

 

斜边模块

Imatest Slanted-Edge 模块包括 SFR、SFRplus、eSFR ISO、Checkerboard 和 SFRreg(详见表 2 和清晰度模块)。

注意请参阅“如何使用 Imatest 测试镜头”,了解如何使用SFRplus 或 eSFR ISO 测量 MTF。

表 2. Imatest斜边模块的简要总结。

模块描述例子
SFR
  • 从各种图表中的倾斜边缘以及任何有干净边缘的地方测量 MTF;区域选择是手动的。
  • 旧 ISO-12233:2000 图表中的两个有用区域  由红色和蓝色箭头指示(请参见示例列中的顶部最右侧图像)。 
  • 一个典型的区域(垂直边缘倾斜约 5.7 度的裁剪)用于计算水* MTF 响应(参见示例列中的下图,最右侧的图像)。

SFR+
  • 从 Imatest SFRplus 图表(推荐)测量 MTF 和其他图像质量参数,或使用 Imatest 测试图表 (宽体打印机、高级打印技能和颜色管理知识)创建。 

  • 与旧的 ISO 12233  :2000 测试图相比,具有许多优势 :自动特征检测、降低对比度以提高准确性、更多边缘(更少浪费的空间)用于图像表面上的 MTF 详细图。 

  • 测量符合 ISO 标准;包括自动区域检测。

eSFR ISO
  • 使用增强版ISO 12233:2014 和 2017  Edge SFR (E-SFR) 测试图测量 MTF 和其他图像质量参数

  • 空间细节比 SFRplus 略少,但 噪点 细节更多。 

  • 包括自动区域检测。

棋盘
  • 对取景敏感,非常适合通过焦点测试。

  • 提供精确的失真计算。

  • 包括自动区域检测。

SFRreg  几个单独的图表通常放置在图像字段周围;适用于:
  • 极致鱼眼镜头 (>180º)

  • 不同距离的图表以测试焦点和景深。

  • 超高分辨率 (>36MP) 相机、大视野和远距离。

  • 包括自动区域检测。

 

斜边算法

MTF 计算源自 ISO 标准 12233 Imatest计算 包含许多增强功能,如下所列当检查SFR输入对话框中 的ISO标准SFR复选框时,执行原始ISO计算 (我们建议未经检查,除非它专门需要)。 

  • 裁剪后的图像被线性化;即,调整像素级别以消除相机应用的 伽马编码。(Gamma 可调整,默认值为 0.5)。
  • 红色、绿色、蓝色和亮度 ( Y ) 通道的边缘位置:
    Y = 0.2125R + 0.7154G + 0.0721B(默认)或0.3R + 0.59G + 0.11B (可在选项 III 中选择)
    确定每条扫描线(上图中的水*线)。
  • 使用多项式回归为每个通道计算边缘的二阶拟合。二阶拟合消除了镜头畸变的影响。在上图中,方程的形式为:
    =一种0+一种1是的+一种2是的2
  • 根据扫描线i的*数部分的值, 
    fp = x i – int(x i ) 
    在每条扫描线处的二阶拟合,移位的边缘被添加到四个 bin 之一: 
    bin 1 if 0 ≤ fp < 0.25
    bin 2 如果0.25 ≤ fp < 0.5
    bin 3 如果0.5 ≤ fp < 0.75
    bin 4 如果0.75 ≤ fp < 1

注意前面等式中提到的 bin 不依赖于检测到的边缘位置。

    • 将四个 bin 组合起来计算*均 4x 过采样边缘。这允许分析超出正常奈奎斯特频率的空间频率。
      • 计算*均 4x 过采样边缘的导数 (d/dx)。应用居中的汉明窗以强制导数在其极限处为零。
      • MTF 是加窗导数的傅里叶变换 (FFT) 的绝对值。

注意:Imatest 斜边 SFR 计算的起源改编自 Matlab 程序 sfrmat,该程序由Peter Burns编写,用于实施 ISO 12233:2000 标准。Imatest 的 SFR 计算包含许多改进,包括改进的边缘检测、更好的镜头失真处理和更好的抗噪性。原始 Matlab 代码可在此处获得。在将 sfrmat 结果与 Imatest 进行比较时,假定音调响应是线性的;即,如果没有 OECF(色调响应曲线)文件输入到 sfrmat ,则gamma = 1。由于 Imatest 中 gamma 的默认值为 0.5,这是 sRGB 等标准色彩空间中数码相机的典型值,因此您必须将 gamma 设置为 1 以获得与 sfrmat 的良好一致性。
 

Imatest 和 ISO 12233 计算之间的差异

位于斜边设置窗口左下方Imatest/ISO 标准 SFR下拉菜单(以前是一个复选框)允许您在 Imatest 和 ISO 兼容计算之间进行选择。推荐设置 #2(ISO 12233 2017 及更早版本),因为Imatest和更新的 ISO 计算更准确 - 在存在噪声和光学失真的情况下绝对优越。

      • 每条扫描线的中心是根据Imatest 计算中低通滤波边缘导数的峰值计算的ISO 计算使用质心,这在没有噪声的情况下是最佳的。但是噪声 总是在某种程度上存在,并且质心对噪声非常敏感,因为距离边缘很远的噪声与边缘本身具有相同的权重。低通滤波器更接*于匹配滤波器,它可以最佳地检测边缘导数峰值。
      • Gamma (用于线性化数据)作为输入值输入或从已知的图表对比度中导出。在 ISO 标准实现中,除非输入 OECF 文件,否则假定为 1。
      • Imatest假设边缘可能由于光学失真而具有一些(二阶)曲率。到 2017 年的 ISO 标准计算假定为直线,这可能导致在存在光学失真的情况下降低 MTF 测量值。弯曲边缘(五阶拟合)包含在 ISO 12233:2022+ 计算中。与旧 Imatest 计算的二阶拟合几乎没有区别。
      • Imatest 的 “改进的变迹”降噪 (默认开启)导致更一致的测量(无系统差异)。检查 ISO 标准计算时关闭
      • ISO 12233:2014 中引入的线扩展函数 ( LSF ) 校正因子自 2015年以来已在Imatest 中实施。校正因子D(j)(或D(k)),用于补偿数值微分时的高频损失从边缘扩展函数 ( ESF )计算LSF在 ISO 12233:2014 和 2017 标准中都是不正确的。Imatest和更新的 ISO 实施基于第一原则(我们拥有全套方程式),符合两个版本标准的意图。 
      • Tukey 窗口用于较新的 Imatest 和 ISO 计算。它似乎没有什么不同。

请注意,可以在 Peter Burns 链接中找到其他计算详细信息 (如下)

 

斜边降噪

Imatest 改进的变迹技术降低了噪声,使 MTF 结果更加一致,同时对 MTF 测量的影响最*。它通过在距边缘中心一定距离但不在中心附**滑线扩散函数(LSF;边缘的导数)来工作。因为它对*均 MTF 几乎没有影响,所以除非结果需要严格符合 ISO 标准,否则它应该保持不变。单击下面的按钮以获取完整说明。

 

改进的变迹降噪技术可用于斜边测量 (SFR、  SFRplus、  eSFR ISO、  SFRreg和 棋盘格)。它可以提高噪声图像的测量一致性,特别是在高空间频率(f > Nyquist/2)下,但对低噪声图像影响不大。当在任何斜边模块的设置窗口或 Rescharts更多设置窗口中选中MTF 降噪(修改的变迹)复选框时,将应用修改的变迹必须取消选择 ISO 标准 SFR(窗口左下方)。

注意:Imatest 建议保持降噪(修改变迹)打开。

Apodization来自于 Joseph D. LaVeigne、Stephen D. Burks 和 Brian Nehring 对计算 MTF 的傅立叶变换方法的比较,可在Santa Barbara Infrared 网站上获得。基本假设是所有重要细节(至少对于高空间频率)都靠*边缘(图 1)。原始技术涉及将线扩展函数 (LSF) 设置为零,超出与边缘的指定距离。修改后的技术对 LSF 进行了强烈*滑(低通滤波器),与原始技术相比,它对低频响应的影响要*得多,并且允许设置更严格的边界以更好地降低噪声。

 

图 1. 修改后的变迹:原始噪声*均线扩散函数(底部;绿色),*滑(中间;蓝色),用于 MTF 的 LSF(顶部;红色)

修改的变迹:原始噪声*均线扩散函数(底部;绿色),*滑(中间;蓝色),用于 MTF 的 LSF(顶部;红色)

 

降噪算法

线扩散函数(LSF;*均边缘响应的导数;右图底部的绿色曲线)被*滑(低通滤波)以在中间创建蓝色曲线。*滑是通过取 9 点移动*均(9 个相邻点的*均值)来完成的。 

注意:由于 分箱算法,这些样本被 4 倍过采样,因此它们对应于原始图像中的大约两个样本。*滑消除了大多数高于奈奎斯特频率(0.5 个周期/像素)的响应。

L 和U*滑曲线的幅度大于峰值的20%的区域的边界(x轴界限),即20%的脉冲宽度是这些边界之间的差。

 

图 2. 在没有 (L) 和有 (R) 的噪声图像上改进的变迹降噪

没有 (L) 和 (R) 修改的变迹降噪的噪声图像的边缘/MTF 图

 

PW20 = B U – B L 

变迹边界位于 

L = B L – PW20 – 4   和   A U = B U + PW20 + 4(像素) 

这允许足够的“呼吸空间”,因此边缘附*的重要细节不受影响。

用于计算 MTF 的 LSF 设置为切趾边界 {A L , A U } 的原始(未*滑)LSF 和外部*滑的 LSF ,如上面的红色曲线所示。对于具有强(模拟)白噪声的图像,修改变迹降噪的好处显示在右侧。

 

 

一些相关的技术会影响清晰度结果,包括:

关键要点

    • 频域图和空间域图传达了相似的信息,但形式不同。空间域中的窄边缘对应于频域中的宽频谱(扩展频率响应),反之亦然。
    • Imatest 测量系统响应,其中包括图像处理:不仅仅是镜头响应。
    • 高于奈奎斯特频率的传感器响应会导致混叠,表现为低空间频率的莫尔图案。在拜耳传感器(除 Foveon 之外的所有传感器)中,莫尔图案显示为彩色条纹。Foveon 传感器中的莫尔条纹远不那么麻烦,因为它是单色的,并且红色和蓝色通道的有效奈奎斯特频率低于拜耳传感器。
    • 奈奎斯特频率及以上的 MTF 并不是混叠问题的明确指标。MTF 是镜头和传感器响应、去马赛克算法和锐化的产物,它经常提高奈奎斯特频率的 MTF。MTF 应被解释为可能存在问题的警告
    • 计算 R、G、B 和亮度 (Y) 通道的结果(默认情况下,Y = 0.2125R + 0.7154G + 0.0721B,但可以设置为较旧的 (NTSC) 值, 0.3 R +0.59 G +0.11 B选项 III窗口中)。Y 通道通常显示在前景中,但可以选择任何其他通道。所有内容都包含在 .CSV 输出文件中。
    • CCD 传感器的水*和垂直分辨率可能不同,应单独测量。对于 CMOS 传感器,它们几乎相同。回想一下,水*分辨率是用垂直边缘测量的,垂直分辨率是用水*边缘测量的。分辨率只是影响图像质量的众多标准之一。
    • MTF 可以在整个图像中变化,并且它并不总是遵循中心附*最清晰和角落附*不太清晰的预期模式。原因有很多:镜头未对准、场曲、散焦等。这就是为什么测量很重要。
posted @ 2022-04-20 10:54  立业安邦  阅读(1598)  评论(0编辑  收藏  举报