CN1963430A - 调制传递函数的测量仪及测量方法 - Google Patents

调制传递函数的测量仪及测量方法 Download PDF

Info

Publication number
CN1963430A
CN1963430A CN 200610119331 CN200610119331A CN1963430A CN 1963430 A CN1963430 A CN 1963430A CN 200610119331 CN200610119331 CN 200610119331 CN 200610119331 A CN200610119331 A CN 200610119331A CN 1963430 A CN1963430 A CN 1963430A
Authority
CN
China
Prior art keywords
optical system
transfer function
modulation transfer
measured
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN 200610119331
Other languages
English (en)
Inventor
赵罡
朱健强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN 200610119331 priority Critical patent/CN1963430A/zh
Publication of CN1963430A publication Critical patent/CN1963430A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

一种用于光学***调制传递函数的调制传递函数的测量仪及测量方法,本测量仪由计算机及其液晶显示屏,支架,二维彩色CCD和基座组成。测量原理是利用计算机产生的随机噪声图作为检测输入,通过分析成像前后图像的功率谱密度,本发明只需一次测量就可快速检测得到光学***的MTF曲线。该仪器通过计算机程序改变检测输入图像,不仅具有快速检测的特点,也可以方便分析影响调制传递函数的各种因素。

Description

调制传递函数的测量仪及测量方法
技术领域
本发明涉及光学***,特别是一种用于光学***调制传递函数(modulationtransfer function,以下称MTF)的测量仪及测量方法,是一种通过检测随机图样进行MTF快速检测的仪器。
背景技术
光学***调制传递函数(MTF)是一种描述光学***成像质量的物理量,通过MTF曲线,可以方便的观察到***在不同分辨率(lp/mm)下对比度的变化情况。目前,不论是在一次性的精密光学元件的集成组装还是大批量元器件生产的场合,都需要实现MTF的在线快速检测,传统的机械式结构的MTF检测方法需要对检测元件进行分区域扫描,检测的效率不高,应用受限。Glenn D.Boreman在1986年提出的随机目标图法,通过分析物像频谱得到了MTF曲线。
发明内容
本发明的目的是提供一种快速检测光学***调制传递函数的测量仪及测量方法。它应测量装置结构简化,操作简便,很容易得到MTF曲线;通过调制不同的测试图样改变输入图样的频谱成分,可以对影响MTF曲线的各种因素进行分析,在保证实验条件下精度也可以达到合理的范围,并且整套装置成本较低。
本发明通过以下技术方案来实现:
一种用于光学***的调制传递函数的测量仪,其构成包括计算机,该计算机的输出端接液晶显示器,该液晶显示器、支架和二维彩色CCD位于一基座的一直线上,所述的二维彩色CCD的输出端接所述的计算机的输入端,待检测的光学***置放在所述的支架上,调整使所述的液晶显示器、待检测的光学***和二维彩色CCD同光轴并满足成像关系,所述的计算机具有随机图样生成软件和相应的数据处理软件。
所述的计算机由随机图样生成计算机和数据处理计算机构成。
利用上述的调制传递函数的测量仪进行光学***调制传递函数的测量方法,包括下列步骤:
①测量前***的准备工作:
所述的调制传递函数的测量仪置于避光处,将待测光学***置于所述的支架上,经调整,使液晶显示屏、待测光学***和二维彩色CCD的中心位置位于同一光轴上并满足成像关系;
②成像后图像测量:
首先通过具有随机图样生成软件的计算机生成一幅随机噪声图,在液晶显示屏上显示,通过待测光学***成像在所述的二维彩色CCD上得到成像后图像,并将成像图像传输到计算机内进行显示,调节光学镜头,待成像后图像清晰成像后捕捉一帧图像,通过数据线将所得到的成像后图像传送到所述的具有相应的数据处理软件的数据处理计算机;
③数据处理计算机进行数据处理:计算成像后图像的功率谱密度和待测光学***调制传递函数。
所述的光学***可以使单透镜,也可以是组合透镜,组合式光学镜头,测量时需保证液晶显示屏与基座垂直,调整待测光学***和二维彩色CCD的光轴,使其垂直于液晶显示屏,调节要使得由计算机观察到的图像清晰。通过计算机程序生成的检测图像可以包含不同频率、不同颜色,所以可以很方便对影响调制传递函数的因素进行分析。
本发明具有以下优点:
1.接收端不需要扫描,因而简化了传统的测量装置;操作简便,通过一次测量就可以得到MTF曲线。
2.通过计算机,调制各种图样(包含不同的频谱,颜色)改变输入图像,可以对影响MTF曲线的各种因素进行分析。
3.在保证实验条件下,测量的精度可以控制在合理的范围,并且整套装置成本较低,是一种具有很大前景的测量仪。
附图说明
图1本发明调制传递函数的测量仪的光路结构示意图。
图中:100具有随机图样生成软件的计算机,101液晶显示屏,102待测光学***,103二维彩色CCD,104基座,105数据处理计算机,106支架
图2为本发明采用的随机噪声图的功率谱密度曲线示意图。
图3为成像后的随机噪声图的功率谱密度曲线示意图。
图4为由随机噪声图功率谱密度计算得到的调制传递函数MTF曲线
图5为纯白图像成像后图像的功率谱密度曲线示意图
图6为本发明检测不同颜色输入的MTF曲线示意图
具体实施方式
下面结合附图对本发明测量仪的工作过程作进一步详细描述。
先请参阅图1,图1本发明调制传递函数的测量仪一个实施例的光路结构示意图。由图可见,本发明调制传递函数的测量仪的构成包括随机图样生成计算机100,该计算机的输出端接液晶显示屏101,该液晶显示屏101、支架106和二维彩色CCD103位于一基座104的一直线上,所述的二维彩色CCD103的输出端接数据处理计算机105的输入端,待测的光学***102置放在所述的支架106上,调整使所述的液晶显示屏101、待测的光学***102和二维彩色CCD103同光轴并满足成像关系,所述的计算机具有随机图样生成软件和相应的数据处理软件。
所述的待测的光学***102可为单透镜、或多透镜组合镜头。
利用本发明调制传递函数的测量仪进行光学***调制传递函数的测量方法,包括下列步骤:
①测量前***的准备工作
将所述的调制传递函数的测量仪置于避光处,将待测光学***102置于所述的支架106上,经调整,使液晶显示屏101、待测光学***102和二维彩色CCD103的中心位置位于同一光轴上并满足成像关系;
②成像后图像测量:
如图1,首先通过具有随机图样生成软件的计算机100由Matlab生成一幅随机噪声图,在液晶显示屏101上显示,通过待测光学***102成像在所述的二维彩色CCD103上得到成像后图像,并将成像图像传输到数据处理计算机105内进行显示,调节待测光学***102,待清晰成像后,捕捉一帧成像后图像传送到数据处理计算机(105);
③数据处理计算机(105)进行数据处理:计算成像后图像的功率谱密度和待测光学***调制传递函数。
光学***调制传递函数MTF,由线性***理论,成像前后图像满足公式:
PSDobject*MTF2 sys=PSDimage    (1)
其中:
PSDobject为检测图像的功率谱密度(power spectral density,以下简称PSD),
PSDimage为成像后图像的功率谱密度,
MTFsys为光学***的调制传递函数,
由(1)式可知光学***102的调制传递函数:
MTF sys = PSD image / PSD object - - - ( 2 )
图像的功率谱密度的具体计算方法如下:
长度为N的离散函数g[n],(n=1,2,…N),根据功率谱密度定义,功率谱密度为函数g[n]的傅立叶变换的模平方,再除以数组长度:
PSD g [ n ] = | FFT ( g [ n ] ) | 2 N - - - ( 3 )
设二元函数f(x,y)代表一幅图像,在计算机内,f(x,y)被离散化为二维数组f(m,n),取x方向第i行数据f(i,n),即f(i,1)、f(i,2)....f(i,N),由公式3知,第i行的功率谱密度:
PSD f [ i , n ] = | FFT ( f [ i , n ] ) | 2 N ,在y方向上作M次平均,即得到图像f(x,y)在x轴方向的平均功率谱密度:
PSD x = 1 M Σ i = 1 M PSD f [ i , n ] , - - - ( 4 )
MTF计算过程
采用随机噪声图作为检测图样是因为,根据随机原理,随机噪声的功率谱密度在各频率处值都为1,图2是本发明采用的随机噪声图的功率谱密度曲线示意图,由图2可以看出,采用Matlab生成的随机噪声图的功率谱密度在“1”值附近微小波动。所以一般情况下可以直接采用成像后输出图像的功率谱密度曲线进行MTF计算,对结果的影响很小,从而达到简化测量的目的。
此时,公式2简化为:
MTF sys = PSD image - - - ( 5 )
图3是根据公式4绘制的曲线,即随机噪声图通过光学***成像后的功率谱密度示意图,由公式5知,光学***的x轴方向MTF计算公式: MTF = PSD x ,绘制曲线如图4所示;采用Matlab自带的函数polyfit对原始数据进行3次多项式拟合后得到的光滑曲线如图4所示。
误差分析
本发明的测量***可以对由背景噪声引起的影响进行测量分析。
因为纯色图像的功率谱不包含任何的高频成分,所以采用一幅纯白图像作为检测输入,经***成像后的图像的功率谱密度曲线如图5所示,与图2对比,可以看出外界的噪声影响对于引入光学***的误差最大不超过3%,所以本发明方法的测量能够保证足够的测量精度。
本发明的测量应用
本发明的测量***最大优点是可以通过改变输入图样,达到观察MTF曲线变化趋势的目的,下面举例分析某一镜头对不同波长的光的响应变化情况。
在可见光内,不同波长的光对应不同的颜色,因此可以采取不同颜色的图像作为输入,采用不同波长的光对于待测光学***的输入。如图6,取三幅分别只包含红、蓝、绿颜色的噪声图作为输入检测图像,经***成像后按照上述的计算方法得到三条对应不同波长的***MTF曲线,如图6所示,由图6可以看出,该***对于不同波长的光的响应是略有不同的,在这里是对绿色的光更为敏感。

Claims (4)

1.一种用于光学***调制传递函数检测的调制传递函数的测量仪,特征在于其构成包括计算机,该计算机的输出端接液晶显示器(101),该液晶显示器(101)、支架(106)和二维彩色CCD(103)位于一基座(104)的一直线上,所述的二维彩色CCD(103)的输出端接所述的计算机的输入端,待测的光学***(102)置放在所述的支架(106)上,调整使所述的液晶显示器(101)、待测的光学***(102)和二维彩色CCD(103)同光轴并满足成像关系,所述的计算机具有随机图样生成软件和相应的数据处理软件。
2、根据权利要求1所述的调制传递函数的测量仪,其特征在于所述的由随机图样生成计算机(100)和数据处理计算机(105)构成。
3、根据权利要求1所述的调制传递函数的测量仪,其特征在于所述的待测的光学***(102)为单透镜、或多透镜组合镜头。
4、利用权利要求2所述的调制传递函数的测量仪进行光学***调制传递函数的测量方法,其特征在于包括下列步骤:
①测量前***的准备工作
将权利要求1所述的调制传递函数的测量仪置于避光处,将待测光学***(102)置于所述的支架(106)上,经调整,使液晶显示屏(101)、待测光学***(102)和二维彩色CCD(103)的中心位置位于同一光轴上并满足成像关系;
②成像后图像测量:
首先通过具有随机图样生成软件的计算机(100)生成一幅随机噪声图,在液晶显示屏(101)上显示,通过待测光学***(102)成像在所述的二维彩色CCD(103)上得到成像后图像,并将成像图像传输到数据处理计算机(105)内进行显示,调节待测光学***(102),待清晰成像后捕捉一侦成像后图像传送数据到数据处理计算机(105);
③数据处理计算机(105)进行数据处理:计算成像后图像的功率谱密度和待测光学***调制传递函数。
CN 200610119331 2006-12-08 2006-12-08 调制传递函数的测量仪及测量方法 Pending CN1963430A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200610119331 CN1963430A (zh) 2006-12-08 2006-12-08 调制传递函数的测量仪及测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200610119331 CN1963430A (zh) 2006-12-08 2006-12-08 调制传递函数的测量仪及测量方法

Publications (1)

Publication Number Publication Date
CN1963430A true CN1963430A (zh) 2007-05-16

Family

ID=38082571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200610119331 Pending CN1963430A (zh) 2006-12-08 2006-12-08 调制传递函数的测量仪及测量方法

Country Status (1)

Country Link
CN (1) CN1963430A (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101872059A (zh) * 2010-06-18 2010-10-27 南京理工大学 Otf测试仪的自动调焦***及其自动调焦方法
CN101949769A (zh) * 2010-08-13 2011-01-19 哈尔滨工业大学 动像调制传递函数测量装置
CN101413843B (zh) * 2007-10-19 2011-03-30 鸿富锦精密工业(深圳)有限公司 调制传递函数值的测量装置及其测量方法
CN101625280B (zh) * 2008-11-14 2011-06-29 昆明理工大学 快速测定投影镜头光学传递函数的测量装置
CN102735429A (zh) * 2012-06-13 2012-10-17 中国科学院长春光学精密机械与物理研究所 一种用于ccd调制传递函数测试的设备及其测试方法
CN103217272A (zh) * 2013-03-04 2013-07-24 中国科学院长春光学精密机械与物理研究所 用光学传递函数测量光电平台运动成像质量退化的方法
CN103499432A (zh) * 2013-09-24 2014-01-08 清华大学深圳研究生院 一种成像***的成像质量检测方法及***
CN105987806A (zh) * 2015-02-02 2016-10-05 宁波舜宇光电信息有限公司 一种转折镜头的测试装置及测试方法
CN106153301A (zh) * 2015-03-24 2016-11-23 北京威斯顿亚太光电仪器有限公司 一种用于硬性内窥镜调制传递函数的检测方法
CN106199379A (zh) * 2016-06-28 2016-12-07 南京理工大学 图像融合电路调制传递函数测试***
CN107421722A (zh) * 2017-07-01 2017-12-01 南京理工大学 基于ccd相机像增强器调制传递函数测试方法及***
CN110736556A (zh) * 2019-10-21 2020-01-31 中国科学院上海光学精密机械研究所 多波长光场能量测量方法
CN112229603A (zh) * 2020-09-07 2021-01-15 中国科学院上海光学精密机械研究所 二维光学传递函数测量装置及方法
CN112880975A (zh) * 2021-01-14 2021-06-01 歌尔股份有限公司 调制传递函数测试装置

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101413843B (zh) * 2007-10-19 2011-03-30 鸿富锦精密工业(深圳)有限公司 调制传递函数值的测量装置及其测量方法
CN101625280B (zh) * 2008-11-14 2011-06-29 昆明理工大学 快速测定投影镜头光学传递函数的测量装置
CN101872059A (zh) * 2010-06-18 2010-10-27 南京理工大学 Otf测试仪的自动调焦***及其自动调焦方法
CN101949769A (zh) * 2010-08-13 2011-01-19 哈尔滨工业大学 动像调制传递函数测量装置
CN101949769B (zh) * 2010-08-13 2012-01-11 哈尔滨工业大学 动像调制传递函数测量装置
CN102735429A (zh) * 2012-06-13 2012-10-17 中国科学院长春光学精密机械与物理研究所 一种用于ccd调制传递函数测试的设备及其测试方法
CN102735429B (zh) * 2012-06-13 2014-08-20 中国科学院长春光学精密机械与物理研究所 一种用于ccd调制传递函数测试的设备及其测试方法
CN103217272A (zh) * 2013-03-04 2013-07-24 中国科学院长春光学精密机械与物理研究所 用光学传递函数测量光电平台运动成像质量退化的方法
CN103499432A (zh) * 2013-09-24 2014-01-08 清华大学深圳研究生院 一种成像***的成像质量检测方法及***
CN103499432B (zh) * 2013-09-24 2015-08-12 清华大学深圳研究生院 一种成像***的成像质量检测方法及***
CN105987806A (zh) * 2015-02-02 2016-10-05 宁波舜宇光电信息有限公司 一种转折镜头的测试装置及测试方法
CN105987806B (zh) * 2015-02-02 2019-12-20 宁波舜宇光电信息有限公司 一种转折镜头的测试装置及测试方法
CN106153301A (zh) * 2015-03-24 2016-11-23 北京威斯顿亚太光电仪器有限公司 一种用于硬性内窥镜调制传递函数的检测方法
CN106153301B (zh) * 2015-03-24 2018-07-06 北京威斯顿亚太光电仪器有限公司 一种用于硬性内窥镜调制传递函数的检测方法
CN106199379A (zh) * 2016-06-28 2016-12-07 南京理工大学 图像融合电路调制传递函数测试***
CN107421722A (zh) * 2017-07-01 2017-12-01 南京理工大学 基于ccd相机像增强器调制传递函数测试方法及***
CN110736556A (zh) * 2019-10-21 2020-01-31 中国科学院上海光学精密机械研究所 多波长光场能量测量方法
CN110736556B (zh) * 2019-10-21 2021-01-01 中国科学院上海光学精密机械研究所 多波长光场能量测量方法
CN112229603A (zh) * 2020-09-07 2021-01-15 中国科学院上海光学精密机械研究所 二维光学传递函数测量装置及方法
CN112880975A (zh) * 2021-01-14 2021-06-01 歌尔股份有限公司 调制传递函数测试装置
CN112880975B (zh) * 2021-01-14 2023-01-24 歌尔股份有限公司 调制传递函数测试装置

Similar Documents

Publication Publication Date Title
CN1963430A (zh) 调制传递函数的测量仪及测量方法
CN104634449B (zh) 微光iccd信噪比测试***及测试方法
CN104296968B (zh) 多通道ccd的调制传递函数测试方法
CN107403177A (zh) 基于工业相机的亮度测量方法
CN102486402B (zh) 一种测量脉冲激光能量的方法及***
CN105547342A (zh) 基于液晶面板的工业镜头测试装置及方法
CN103024427B (zh) 相机调制传递函数的测试方法及其测试装置
CN103257029A (zh) 无刀口机构可见光纹影仪***及其测量方法
CN107063638B (zh) 一种基于模组化高精度显微镜***的测试装置
CN201043930Y (zh) 光学***调制传递函数测量仪
CN110726372B (zh) 一种准确处理单缝衍射图像方法
CN101673043A (zh) 广角畸变测试***及方法
CN114577444A (zh) 一种调制传递函数的检测方法及检测装置
JP2004191112A (ja) 欠陥検査方法
CN102759788A (zh) 表面多点对焦***及方法
CN102901465A (zh) 一种测量光学镜片斜率误差的设备和方法
KR960011413A (ko) 표시화소의 광량측정방법 및 표시화면의 검사방법 및 검사장치
CN104359654B (zh) 一种光纤传像束两端面像元偏移量的测量方法
CN204202849U (zh) 一种ccd调制传递函数测量装置
CN114235347B (zh) 镜头质量评估方法及装置
WO2022057710A1 (zh) 激光脉冲整形装置及方法、脉冲整形器、光学***
CN106323197B (zh) 便携式阵列调零高精度激光大工作距自准直装置与方法
JPS63163150A (ja) プリント基板検査装置
CN203101029U (zh) 一种用于检测中继镜的装置
CN104280214A (zh) 一种ccd调制传递函数测量装置及测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication