光学传递函数及像质评价实验
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验十一 光学传递函数测量及像质评价实验
光学成像系统是信息(结构、灰度、色彩)传递系统,从物面到像面,输出图像的质量取决于光学系统的传递特性。
在频域中分析光学系统的成像质量时,可以把光学成像系统看成是一个低通空间滤波器,将输入信息分解成各种空间频率分量。
通过考察这些空间频率分量在通过系统的传递过程中丢失、衰减、相位移动等变化,也就是研究系统的空间频率传递特性即光学传递函数(OTF ,Optical Transfer Function ),来获取成像的空间频谱特性。
光学传递函数的性质主要体现在:它定量反映了光学系统的孔径、光谱成分以及像差大小所引起的综合效果;用它来讨论光学系统时,其可靠性依赖于光学系统对线性和空间不变性的满足程度;用它来分析讨论物像之间的关系时,不受试验物形式的限制;可以用各个不同方位的一维光学传递函数来分析处理光学系统,简化了二维处理;它可以根据设计结果进行计算,也能对已制成的光学系统进行测量。
可见,光学传递函数表征光学系统对物体或图像中不同频率的信息成分的传递特征,可用于光学系统成像质量的评价。
本实验利用非相干面光源、光栅、透镜、CCD (Charge-coupled Device ,电荷耦合元件)图像传感器、数据采集和处理系统,测出光学成像系统的光学传递函数曲线图,并对成像质量作出评价。
一、实验目的
1.了解光学传递函数及其测量方法。
2.掌握传递函数测量和像质评价的近似方法。
3.熟悉抽样、平均和统计算法。
二、实验仪器
面光源、凸透镜、CCD 图像传感器、数据采集及处理系统、计算机、导轨(滑块)、调节支座(支架)、干版架、可调节光阑。
三、实验原理
1. 光学传递函数
一个确定的物分布可看成许多个δ函数的线性组合,每个δ函数在像面上均有对应的脉冲响应。
如果是非相干照明,则物面上任意两个脉冲都是非相干的,它们的脉冲响应在像面上也是非相干叠加,也就是强度叠加。
假设非相干成像系统是强度的线性系统,成像空域不变,则该系统物像关系满足以下卷积积分:
000000ˆˆˆˆˆˆ(,)(,)(,)(,)(,)
i i i I i i g i i I i i I x y K I x
y h x x y y dx dy K I x y h x y ∞
∞
-∞-∞=--=⊗⎰
⎰
(1)
式中(,)g i i I x y 是物体000(,)I x y 理想像的强度分布,(,)i i i I x y 是物体000(,)I x y 通过衍射受限系统后成像的强度分布,(,)I i i h x y 是强度脉冲响应,为点物产生的像斑的强度分布。
I I g
i
图1 物、像光强分布与调制度、光学传递函数关系示意图
对式(1)两边做傅里叶变换,得到频域中的物像关系,
ˆˆˆ(,)(,)(,)i g I
G u v G u v H u v = (2) 式中ˆ(,){(,)}i i i i G u v F I x y =,ˆ(,){(,)}g g i i G u v F I x y =,ˆ(,){(,)}I I i i
H u v F h x y =。
由于i I 、g I 、I h 均为非负实函数,必有零频分量。
对频谱进行归一化,可得到
2()ˆ(,)(,)ˆ(,)ˆ(0,0)
(,)i i j ux vy i i i i i
i i i
i
i i i i
I x y e dx dy G u v u v G I x y dx dy π∞∞
-+-∞-∞
∞
∞-∞-∞Γ==
⎰⎰⎰⎰
2()ˆ(,)(,)ˆ(,)ˆ(0,0)
(,)i i j ux vy g i i i i
g g
g
g i i i i
I x y e dx dy G
u v u v G I x y dx dy π∞
∞-+-∞-∞
∞
∞-∞-∞
Γ==⎰⎰
⎰⎰
2()(,)ˆ(,)ˆ(,)ˆ(0,0)
(,)i i j ux vy I i i i i
I I I
I i i i i
h x y e dx dy H u v u v H h x y dx dy π∞∞-+-∞-∞
∞
∞-∞-∞
Γ==⎰⎰
⎰⎰
因此有 ˆˆˆ(,)(,)(,)i g I u v u v u v Γ=ΓΓ,其中(,)
ˆ(,)(,)i u v I
u v M u v e j Γ=称为非相干成像系统的
光学传递函数(OTF ,Optical Transfer Function ),是一个复函数,模部分(,)M u v 为调制传递函数(MTF ,Modulation Transfer Function ),相位部分(,)u v j 为相位传递函数(PTF ,
Phase Transfer Function )。
MTF 和PTF 分别反映空间频率为(,)u v 的余弦分量通过系统的降低度和横向位移量。
根据调制度的定义max min max min ()/()V I I I I =-+知,物体理想像强度分布g I 和通过成像系统之后的像强度分布i I 的调制度分别为max min max min
g g g g g I I V I I -=
+和max min
max min
i i i i i I I V I I -=
+,且
(,)i g V M u v V =。
为方便理解相关概念,图1给出了物、像的光强分布与调制度、光学传递
函数的关系。
对于一维正弦光栅,其透过率是空间的正弦函数,可表示为0()cos(2)I x a b x πν=+,其中物分布初相位为零,ν为x 方向的空间频率,则物的调制度为/g V b a =,物光经过系统后像光强的频谱为
[][]ˆˆˆ()()()ˆ
()()()(
)2ˆˆ(0)()()()()2
i x g x I x
x x x I x I x x x I x
G G H b a H b aH H νννδνδννδνννδνδννδννν=⎛⎫=+-++ ⎪⎝⎭
=+-++ 对上式进行傅里叶逆变换,得到像面光强分布函数()(0)()cos(2)i I I I x ah b h x νπνφ=+-,其中()()i I I h h e
φ
νν-=。
因此,该光学系统对应于空间频率ν余弦输入信号调制传递函数
为()()/(0)I I M h h νν=。
考虑到正弦光栅分划板难于制造,且达不到精度要求,实验中用矩形光栅代替正弦光栅作为成像物,测出系统的对比传递函数,再求出调制传递函数。
取空间频率为ν的矩形光
栅,进行傅里叶展开,得到112sin(2)()2n n x I x n πνπ∞==+∑。
假设成像时系统的调制传递函数
为()i M ν,则矩形光栅像的强度分布为112sin(2)
()()2i n n x I x M n n
πννπ∞==+∑
,像的调制度为411
()(3)(5)35i V M M M νννπ⎡⎤
=
-+-⎢⎥⎣⎦。
假设矩形光栅的调制度为1,则该成像系统的对比传递函数411
()()(3)(5)35
C M M M ννννπ⎡⎤
=
-+-⎢⎥⎣⎦
,从而得到成像系统对应空
间频率ν的调制传递函数11
()()(3)(5)435M M M M πνννν⎡
⎤
=
+-+⎢⎥⎣⎦。
考虑到成像系统都存在截止频率off ν,故实际使用的项数有限。
低频情况下,一般取前三项就达到精度要求。
大于三分之一截止频率时,系统的对比传递函数均为零,此时调制传递函数取
()4
M π
ν。
实验中由于衍射和光学像差的影响,使得实际光学系统的像不再是矩形光栅。
通过对矩形光栅形成的像进行抽样统计,作出直方图,两端极大值h M 和极小值l M 之差,可表示在相应空间频率下调制传递函数MTF 值,如图2所示。
强度
10
x
强度
1
x
(a ) (b )
1强度
概率
M M h
l
(c )
图2 (a )调制度为1的矩形光栅目标函数;(b )对矩形光栅的不完善像进行抽样;(c )直方统计图
2. 利用调制传递函数(MTF )曲线进行像质评价
对数码照片成像素质影响最大的是镜头的分辨率和反差,一般通过拍摄正弦光栅来进行测试。
分辨率的单位是线对/毫米(lp/mm ),相邻的黑白两条线称为一个线对,每毫米能够分辨出的线对数就是分辨率。
正弦光栅的疏密程度被称为空间频率,单位也是lp/mm ,表示单位长度(每毫米)的亮度按照正弦变化的图形的周期数。
在对镜头的反差和分辨率进行测试时,将正弦光栅置于镜头前方,测量镜头成像处的调制度。
这时由于镜头像差的影响,会出现以下情况。
当空间频率很低时,测量出的调制度M 几乎等于正弦光栅的调制度;当所拍
摄的正弦光栅空间频率提高时,镜头成像的调制度逐渐下降。
镜头成像的调制度随空间频率变化的函数称为调制度传递函数(MTF )。
对于原来调制度为0V 的正弦光栅,如果经过镜头到达像平面的像的调制度为1V ,则MTF 函数值为10/V V 。
由此可见, MTF 值越接近1,说明镜头的性能越良好。
MTF 值不但可以反映镜头的反差,也可以反映镜头的分辨率。
由于MTF 值排除了成像介质(胶片或感光元件)的影响,因此较为客观。
当空间频率很低时,MTF 趋于1,这时的MTF 值可以反映镜头的反差。
当空间频率提高,也就是正弦光栅的密度提高时,MTF 值逐渐下降,这时的MTF 曲线可以反映镜头的分辨率。
由于人眼的分辨能力有限,一般取MTF 值为0.03时的空间频率作为镜头的目视分辨率极限。
空间频率低于这个值时,镜头成像素质的变化人眼难以察觉。
a b c
MTF 值
空间频率
1.0
0.03
图3 MTF 曲线示意图 图4 Cabon 公司公布的MTF 曲线图
典型的MTF 曲线图见图3,其中曲线a 所代表的镜头在低频段反差适中,但随着空间频率的提高,它的衰减过程比较慢,说明其像质不错;曲线b 所代表的镜头在低频表现较好,说明镜头的反差较好,但随着空间频率的提高,它的曲线衰减很快,说明镜头的分辨率不算很好;曲线c 所代表的镜头在低频时就出现衰减,综合像质较低。
专业绘制MTF 曲线时,需要固定空间频率和光圈,其中固定空间频率低频(10线对/mm )曲线代表镜头反差特性,这条曲线越高,镜头反差越大。
而固定高频(30线对/mm )曲线代表镜头分辨率特性,这条曲线越高,镜头分辨率越高。
虽然纵坐标还是MTF 值,但横坐标改为了像场中心到测量点的距离,见图4。
镜头是以光轴为中心的对称结构,中心向各方向的成像素质变化规律是相同的。
由于像差等因素的影响,像场中某点与像场中心的距离越远,其MTF 值一般呈下降的趋势。
因此以像场中心到像场边缘的距离为横坐标,可以反映镜头边缘的成像素质。
另外,在偏离像场中心的位置,由沿切线方向的线条与沿径向方向的线条的正弦光栅所测得的MTF 值并不一样。
将平行于直径的线条产生的MTF
曲线称为弧矢曲线,标
为S(Sagittal),而将平行于切线的线条产生的MTF曲线称为子午曲线,标为M(Meridional)。
图4中实线为S方向测出,虚线为M方向测出,黑色曲线对应最大光圈得,蓝色曲线对应F8光圈,粗线是空间频率为10 线对/毫米时测得,细线是空间频率为30线对/毫米时测得。
因此,在做像质评价时,需要注意以下几点:(1)MTF 曲线与横轴、纵轴所围成的空间,面积越大越好。
MTF 曲线越平直越好,平直性说明镜头边缘和中心部分的成像均匀性;
(3)低频(10LP/mm)曲线代表镜头的反差特性,高频(30LP/mm)(2)S曲线与M 曲线越接近越好;
曲线代表镜头的分辨率特性;(4)不要将不同焦距段、不同档次、不同规格以及定焦、变焦等不同镜头的MTF图进行比较,因为它们的特性受设计规格、光学特性、像差以及成本、用料的影响很大,不具有可比性。
只有同档次、同规格镜头的MTF图才具有比较意义;(5)在不同测试环境下,镜头MTF曲线图会不一致。
四、实验内容
透镜
图4 成像系统示意图
1. 成像系统调整
(1)如图4,在导轨上调整面光源、目标板、透镜、CCD(白屏)共轴。
(2)让光源滑座与光栅滑座紧贴,在目标板上放置25线对的光栅,调节光栅、透镜、CCD(白屏)之间的距离,在CCD(白屏)上得到光栅的清晰像。
2. MTF测定和数据处理
(1)取出25线对成像用光栅,调节光源强度,采集峰值,取峰值大小为200左右。
(2)放置镀铬玻璃测出暗电背景,得到灰度值下限。
装上K9透明水晶玻璃测峰值,得到灰度值上限。
(3)依次横向放置10、25、50、80线对光栅,在采集系统“设置参数”处改变各通道模式和线对数,采集相应数据。
采集过程中,保证CCD有一定的响应时间(AES电子快门选择OFF,背光补偿选择OFF),确保数据稳定。
设定采样步长时,线对数越大,步长应越
小。
设定CCD数据处理有效区域时,应去除CCD感光边缘数据。
(4)依次纵向放置10、25、50、80线对光栅,同时在采集系统“设置参数”处改变各通道模式和线对数,采集相应数据。
3. 像质评价
根据获取的MTF图评定透镜像质,包括低频清晰度、高频分辨率、不同模式MTF的比较。
五、实验提示
1. 测定MTF曲线图时,应在暗环境下进行,以避免杂散光的影响。
2. 光栅成像时,实际观测到的图像与光栅偏差较大,这与成像透镜孔径有限(高频信息丢失)、透镜本身像差、光路调整有关。
六、预习题
1. 什么是光学成像?为什么会产生成像失真?
2. 解释25lp/mm 的正弦光栅。
七、作业题
1. 光学传递函数、调制传递函数在光学镜头评价中有何作用?
2. 光学传递函数与目标物的空间频率有何关系?结合实验加以分析。