CN1963430A - 调制传递函数的测量仪及测量方法 - Google Patents
调制传递函数的测量仪及测量方法 Download PDFInfo
- 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
Links
- 238000012546 transfer Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 11
- 230000003287 optical effect Effects 0.000 claims abstract description 54
- 238000003384 imaging method Methods 0.000 claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims description 17
- 239000004744 fabric Substances 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 abstract description 9
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000004590 computer program Methods 0.000 abstract description 2
- 238000001444 catalytic combustion detection Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 206010038743 Restlessness Diseases 0.000 description 1
- 239000012491 analyte Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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的调制传递函数:
图像的功率谱密度的具体计算方法如下:
长度为N的离散函数g[n],(n=1,2,…N),根据功率谱密度定义,功率谱密度为函数g[n]的傅立叶变换的模平方,再除以数组长度:
设二元函数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行的功率谱密度:
MTF计算过程
采用随机噪声图作为检测图样是因为,根据随机原理,随机噪声的功率谱密度在各频率处值都为1,图2是本发明采用的随机噪声图的功率谱密度曲线示意图,由图2可以看出,采用Matlab生成的随机噪声图的功率谱密度在“1”值附近微小波动。所以一般情况下可以直接采用成像后输出图像的功率谱密度曲线进行MTF计算,对结果的影响很小,从而达到简化测量的目的。
此时,公式2简化为:
图3是根据公式4绘制的曲线,即随机噪声图通过光学***成像后的功率谱密度示意图,由公式5知,光学***的x轴方向MTF计算公式:
,绘制曲线如图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)进行数据处理:计算成像后图像的功率谱密度和待测光学***调制传递函数。
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)
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 | 歌尔股份有限公司 | 调制传递函数测试装置 |
-
2006
- 2006-12-08 CN CN 200610119331 patent/CN1963430A/zh active Pending
Cited By (21)
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 |