CN106840617B - 变频栅条靶标及其光电成像***动态传递函数测量方法 - Google Patents
变频栅条靶标及其光电成像***动态传递函数测量方法 Download PDFInfo
- Publication number
- CN106840617B CN106840617B CN201710119156.5A CN201710119156A CN106840617B CN 106840617 B CN106840617 B CN 106840617B CN 201710119156 A CN201710119156 A CN 201710119156A CN 106840617 B CN106840617 B CN 106840617B
- Authority
- CN
- China
- Prior art keywords
- grizzly bar
- target
- frequency
- width
- bar
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于光电成像***成像质量性能评价技术领域,具体涉及一种变频栅条靶标及其基于变频靶标的光电成像***动态传递函数测量方法。本发明包括变频栅条靶标,靶标采用矩形栅条靶标图案,图案的整体分布为栅条的空间频率逐渐增大或逐渐减少。本发明通过对单个栅条靶标图案宽度进行有规律地逐级递增或递减设计,结合三杆靶法,应用栅条靶标图像傅里叶变换算法,解算出不同频率的MTF值。这种方法可以实现从一幅靶标采集图像中,解算出光电成像***从低频到高频较宽频率范围的多点MTF值。
Description
一、技术领域
本发明属于光电成像***成像质量性能评价技术领域,具体涉及一种变频栅条靶标及其基于变频靶标的光电成像***动态传递函数测量方法。
二、背景技术:
光学***调制传递函数(MTF)与星点检验、分辨率、几何像差等像质检验方法相比,可以全面反映成像***的成像性能,是被广泛采用的评价光学***成像质量的标准方法之一。根据MTF检测时所采用的目标靶标不同,可将MTF检测分为点光源法、狭缝法、刀口法及栅条法等。除栅条法外,其他方法MTF均是通过光学***点扩散函数傅里叶变换求解。这些方法的前提条件是光学***的***模型满足线性空间不变***。光学***与数字成像传感器组合形成光电成像***。数字成像***图像采集过程对图像二维空间坐标进行采样,最后输出数字图像。理论上,数字成像***不是线性空间不变***。只有当数字成像传感器的像元尺寸小于光学***最小可分辨空间尺寸的二分之一以上时,可以用线性空间不变***模型近似。
动态调制传递函数(MTF)是成像***观测运动目标、或自身载体摇摆振动时、或两种运动同时存在时针对动态场景成像能力的评估参数之一。由于光电成像***数字成像传感器具有一定的时延特性,也就是光电转换需要一定的时间。当目标运动时,光波波前相对于单个光电成像单元的稳定程度较静止目标时有明显下降。因此相同场景下运动目标相当于静目标,图像的对比度降低。从而引起光电成像***的成像分辨率、光电成像***的可探测和可识别能力下降。
栅条法(如三杆靶法)的数据解算基于图像中灰度的最大值与最小值,对***模型没有严格要求,常用于动态MFT的测量。另外基于栅条法的旋转靶标法(如:专利201710004012.5)也是一种检测动态调制传递函数的有效方法。现有的栅条法每组条纹图只能表示一种频率,若需进行不同频率测量时,需制作多组栅条或多个靶标,通过对不同频率栅条图像进行分析,获得MTF曲线上点的值。每次测量所能测量MTF曲线上点个数受栅条组数的限制,对于大频率范围内MTF测量时,需要多个靶标,且需进行靶标切换。
三、发明内容
本发明针对大频率范围内动态MTF测量现有技术创造的缺陷或不足,提供一种变频栅条靶标及其基于变频靶标的光电成像***动态传递函数测量方法,该方法通过对单个栅条靶标图案宽度进行有规律地逐级递增或递减设计,结合三杆靶法,应用栅条靶标图像傅里叶变换算法,解算出不同频率的MTF值。这种方法可以实现从一幅靶标采集图像中,解算出光电成像***从低频到高频较宽频率范围的多点MTF值。
为实现上述目的,本发明采用的技术方案为:一种变频栅条靶标,其特征在于:包括变频栅条靶标,靶标采用矩形栅条靶标图案,图案的整体分布为栅条的空间频率逐渐增大或逐渐减少;
取任意相邻的四条栅条图案,中间两个栅条的宽度相等均为w(2w是所定区域三杆靶栅条的标准周期),左边栅条的宽度为w-Δ,与左边栅条相邻的非栅条区域宽度为w,右边栅条的宽度为w+Δ,与右边栅条相邻的非栅条区域宽度为w+Δ,Δ为栅条周期宽度变化增量;
取任意相邻三条栅条图案,存在频率递减和频率递增两种情况,三杆靶频率递增图案的右边两个栅条的宽度相等均为w,栅条间的非栅条区域的宽度为w,左边栅条的宽度为w-Δ;三杆靶频率递减图案的左边两个栅条的宽度相等,均为w,左边两栅条间的非栅条区域的宽度为w,右边栅条的宽度为w+Δ,右边栅条与中间栅条间的非栅条区域宽度为w+Δ。
一种基于变频靶标的光电成像***动态传递函数测量方法包括如下步骤:
步骤1、根据对象MTF测试要求,设定栅条周期宽度变化增量Δ,控制MTF曲线检测点的位置,通过空间光调制器编程实现或制作实物变频栅条靶标;
步骤2、在靶标沿栅条垂直方向运动情况下,采集栅条靶标图像;
步骤3、从采集图像中应用图像处理算法从左向右或从右向左连续截取形成三杆靶图像;
步骤4、从所截取的三杆靶图像中沿栅条垂直方向选取一组或多组图像灰度值形成图像灰度值数组;
步骤5,对每一组图像灰度值数组进行一维离散傅里叶变换,并对变换结果进行取模运算,获得图像灰度值数据的频域幅值曲线M像;
步骤6,从频率为零点出发,沿频率增大的方向识别图像灰度值数据的频域幅值曲线值为零的第一个点,确定其频率11,记为
步骤7,从图像灰度值数据的频域幅值曲线上确定频率为f处的M像值13,记为M像(f);若步骤4选取的是多组图像灰度值,则对各组图像灰度值数组所求得的M像(f)进行幅值求平均运算;
步骤8,计算所截取的频率逐级变化的所有三杆靶图像的M像;
步骤9,应用公式(M物(f)表示应用上述步骤3至步骤7解算的靶标物方频率为f的物方调制度),计算出所截取的三杆靶图像代表的所有频率点的MTF值。
与现有技术相比,本发明具有如下优点和效果:
本发明通过频率有规律的逐级递增或递减栅条靶标,应用直线运动平台带动靶标往复运动,采集靶标模拟目标在直线运动下的图像,结合MTF测量的三杆靶法,应用栅条靶标图像傅里叶变换算法,从所采集图像中解算出MTF值。这种方法可以实现从一幅靶标采集图像中,解算出光电成像***从低频到高频较宽频率范围内多点MTF值,且MTF曲线各采样点的值可通过靶标栅条周期宽度变化增量Δ进行设计设定。
本发明较之等频率的栅条靶标,易于实现宽频率范围多采样点的MTF测量。
四、附图说明:
图1为变频栅条靶标图案;
图2为任意相邻的4条靶标图案;
图3为相邻三杆靶标(频率递减)图案;
图4为相邻三杆靶标(频率递增)图案;
图5为被检产品MTF测量时的***组成图;
图6为标准三杆靶与三杆靶(频率递增)图像的傅里叶变换幅值曲线。
附图标记说明:1、变频靶标,2、任意相邻的4条靶标图区域,3、相邻三杆靶标(频率递减)图区域,4、相邻三杆靶标(频率递增)图区域,5、积分球,6、变频栅条靶标,7、直线运动平台,8、平行光管,9、被检光电成像***,10、振动摇摆台,11、从频率为零点出发,沿频率增大的方向识别图像灰度值数据的频域幅值曲线值为零的第一个点12、标准三靶标图像灰度数组傅里叶变换曲线(实线),13、M像(f),14、三靶标(频率递增)图像灰度数组傅里叶变换曲线(点划线)。
五、具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以上所述,仅为本发明较佳实施例而已,并非用于限定本发明的保护范围。
本发明一种变频栅条靶标,
靶标1采用矩形栅条靶标图案,从图案的整体分布来看,栅条的空间频率逐渐增大或逐渐减少(参见图1)。
取任意相邻的四条栅条图案2,中间两个栅条的宽度相等均为w。左边栅条的宽度为w-Δ,与左边栅条相邻的非栅条区域宽度为w(2w是所定区域三杆靶栅条的标准周期)。右边栅条的宽度为w+Δ,与右边栅条相邻的非栅条区域宽度为w+Δ。Δ为栅条周期宽度变化增量(参见图2)。
取任意相邻三条栅条图案,可能存在频率递减3和频率递增4两种情况,三杆靶(频率递增)图案的右边两个栅条的宽度相等,均为w,栅条间的非栅条区域的宽度为w,左边栅条的宽度为w-Δ。三杆靶(频率递减)图案的左边两个栅条的宽度相等,均为w,左边两栅条间的非栅条区域的宽度为w,右边栅条的宽度为w+Δ,右边栅条与中间栅条间的非栅条区域宽度为w+Δ(参见图3和图4)。
一种基于变频靶标的光电成像***动态传递函数测量方法包括如下步骤:
步骤1、根据对象MTF测试要求,设定栅条周期宽度变化增量Δ,控制MTF曲线检测点的位置,通过空间光调制器编程实现或制作实物变频栅条靶标。
步骤2、在靶标沿栅条垂直方向运动情况下,采集栅条靶标图像。
步骤3、从采集图像中应用图像处理算法从左向右或从右向左连续截取形成三杆靶图像。
步骤4、从所截取的三杆靶图像中沿栅条垂直方向选取一组(或多组)图像灰度值形成图像灰度值数组。
步骤5,对每一组图像灰度值数组进行一维离散傅里叶变换,并对变换结果进行取模运算,获得图像灰度值数据的频域幅值曲线M像。
步骤6,从频率为零点出发,沿频率增大的方向识别图像灰度值数据的频域幅值曲线值为零的第一个点,确定其频率11,记为
步骤7,从图像灰度值数据的频域幅值曲线上确定频率为f处的M像值13,记为M像(f)。若步骤4选取的是多组图像灰度值,则对各组图像灰度值数组所求得的M像(f)进行幅值求平均运算。
步骤8,计算所截取的频率逐级变化的所有三杆靶图像的M像。
步骤9,应用公式(M物(f)表示应用上述步骤3至步骤7解算的靶标物方频率为f的物方调制度),计算出所截取的三杆靶图像代表的所有频率点的MTF值。
参见图5:靶标安装在直线运动平7的滑块上置于积分球5和平行光管8之间,积分球发出的光透过靶标,穿过靶标射入平行光管,通过分划板与焦平面一致性调整技术,靶标平面位于平行光管焦平面上,并与焦平面高度重合,积分球的出光口径大于靶标平面,同时大于平行光管的有效靶面。被检光电成像***9安装固定在振动摇摆台10上,通过调整其口径对准平行光管8的物镜通光口,当直线运动平台带动靶标做往复运动时,可为被检光电成像***提供无穷远的直线运动目标。
积分球5的出光口具有极高的光照均匀性,一般要求其出光口处光照均匀性不小于98%,确保不同区域靶标物方调制度M物具有极高的一致性。
针对被检产品测量要求,直线运动导轨带动靶标做匀速和匀加速往复运动。匀速往复运动是指在一个周期内,平台带动靶标作加速、匀速减速运动,当靶标经过平行光管靶面时,靶标处于匀速运动状态,匀加速往复运动是指在一个周期内,平台带动靶标作加速减速运动,当靶标经过平行光管靶面时,靶标处于匀加速运动状态。
当靶标运动经过平行光管靶面区域时,被检成像***采集靶标图像,本发明步骤的步骤3至步骤9,解算出MTF曲线上采样点的值。
标准三杆靶图像灰度频谱幅值过零点,可应用下面公式求解:
X表示栅条周期,即X=2w。本发明中截取形成的三杆靶中,存在栅条周期宽度变化增量Δ,三杆靶图像灰度频谱幅值过零点与标准三杆靶位置相比,在标准三杆靶图像灰度频谱幅值过零点附近向低频或高频产生了微小的移动(如图6)。
Claims (1)
1.一种变频栅条靶标,其特征在于:包括变频栅条靶标(1),靶标(1)采用矩形栅条靶标图案,图案的整体分布为栅条的空间频率逐渐增大或逐渐减少;
取任意相邻的四条栅条图案,中间两个栅条的宽度相等均为w,2w是所定区域三杆靶栅条的标准周期,左边栅条的宽度为w-Δ,与左边栅条相邻的非栅条区域宽度为w,右边栅条的宽度为w+Δ,与右边栅条相邻的非栅条区域宽度为w+Δ,Δ为栅条周期宽度变化增量;
取任意相邻三条栅条图案,存在频率递减和频率递增两种情况,三杆靶频率递增图案的右边两个栅条的宽度相等均为w,栅条间的非栅条区域的宽度为w,左边栅条的宽度为w-Δ;三杆靶频率递减图案的左边两个栅条的宽度相等,均为w,左边两栅条间的非栅条区域的宽度为w,右边栅条的宽度为w+Δ,右边栅条与中间栅条间的非栅条区域宽度为w+Δ;
基于变频靶标的光电成像***动态传递函数测量方法包括如下步骤:
步骤1、根据对象MTF测试要求,设定栅条周期宽度变化增量Δ,控制MTF曲线检测点的位置,通过空间光调制器编程实现或制作实物变频栅条靶标;
步骤2、在靶标沿栅条垂直方向运动情况下,采集栅条靶标图像;
步骤3、从采集图像中应用图像处理算法从左向右或从右向左连续截取形成三杆靶图像;
步骤4、从所截取的三杆靶图像中沿栅条垂直方向选取一组或多组图像灰度值形成图像灰度值数组;
步骤5,对每一组图像灰度值数组进行一维离散傅里叶变换,并对变换结果进行取模运算,获得图像灰度值数据的频域幅值曲线M像;
步骤6,从频率为零点出发,沿频率增大的方向识别图像灰度值数据的频域幅值曲线值为零的第一个点,确定其频率11,记为
步骤7,从图像灰度值数据的频域幅值曲线上确定频率为f处的M像值13,记为M像(f);若步骤4选取的是多组图像灰度值,则对各组图像灰度值数组所求得的M像(f)进行幅值求平均运算;
步骤8,计算所截取的频率逐级变化的所有三杆靶图像的M像;
步骤9,应用公式,M物(f)表示应用上述步骤3至步骤7解算的靶标物方频率为f的物方调制度,计算出所截取的三杆靶图像代表的所有频率点的MTF值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710119156.5A CN106840617B (zh) | 2017-03-02 | 2017-03-02 | 变频栅条靶标及其光电成像***动态传递函数测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710119156.5A CN106840617B (zh) | 2017-03-02 | 2017-03-02 | 变频栅条靶标及其光电成像***动态传递函数测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106840617A CN106840617A (zh) | 2017-06-13 |
CN106840617B true CN106840617B (zh) | 2019-01-01 |
Family
ID=59138673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710119156.5A Active CN106840617B (zh) | 2017-03-02 | 2017-03-02 | 变频栅条靶标及其光电成像***动态传递函数测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106840617B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111289224B (zh) * | 2020-01-17 | 2021-08-06 | 中国科学院上海技术物理研究所 | 一种自动识别ctf测试图像并计算数值的方法 |
CN111896230A (zh) * | 2020-08-10 | 2020-11-06 | 西安工业大学 | 一种光电成像***动态mtf传递函数检测*** |
CN116847067B (zh) * | 2023-08-30 | 2023-11-24 | 荣耀终端有限公司 | 成像***评估方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5303023A (en) * | 1992-03-26 | 1994-04-12 | Allergan, Inc. | Apparatus and method for inspecting a test lens, method of making a test lens |
CN1620112A (zh) * | 2003-11-21 | 2005-05-25 | 中国印钞造币总公司 | 多功能靶标 |
CN2781403Y (zh) * | 2004-12-22 | 2006-05-17 | 南开大学 | 长周期光纤光栅应力模板 |
CN1983036A (zh) * | 2005-11-22 | 2007-06-20 | Asml荷兰有限公司 | 作为对准标记的二进制正弦子波长格栅 |
CN102162739A (zh) * | 2010-12-30 | 2011-08-24 | 中国科学院长春光学精密机械与物理研究所 | 空间相机在轨动态传函测试方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6900884B2 (en) * | 2001-10-04 | 2005-05-31 | Lockheed Martin Corporation | Automatic measurement of the modulation transfer function of an optical system |
-
2017
- 2017-03-02 CN CN201710119156.5A patent/CN106840617B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5303023A (en) * | 1992-03-26 | 1994-04-12 | Allergan, Inc. | Apparatus and method for inspecting a test lens, method of making a test lens |
CN1620112A (zh) * | 2003-11-21 | 2005-05-25 | 中国印钞造币总公司 | 多功能靶标 |
CN2781403Y (zh) * | 2004-12-22 | 2006-05-17 | 南开大学 | 长周期光纤光栅应力模板 |
CN1983036A (zh) * | 2005-11-22 | 2007-06-20 | Asml荷兰有限公司 | 作为对准标记的二进制正弦子波长格栅 |
CN102162739A (zh) * | 2010-12-30 | 2011-08-24 | 中国科学院长春光学精密机械与物理研究所 | 空间相机在轨动态传函测试方法及装置 |
Non-Patent Citations (1)
Title |
---|
基于三线靶标法的在轨MTF评价方法;王琦等;《遥感信息》;20110215;第14-17页,附图2 |
Also Published As
Publication number | Publication date |
---|---|
CN106840617A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106840617B (zh) | 变频栅条靶标及其光电成像***动态传递函数测量方法 | |
CN104160241B (zh) | 利用高维亮度信息的条纹图像的相位分布分析方法、装置及其程序 | |
KR20160018744A (ko) | 규칙성 모양에 의한 변위 분포를 위한 측정방법, 장치 및 그 프로그램 | |
CN106802230B (zh) | 一种光电成像***动态传递函数测量方法及靶标 | |
US20180094918A1 (en) | Measurement method, measurement apparatus, measurement program and computer readable recording medium in which measurement program has been recorded | |
KR20100087103A (ko) | 주기적 패턴 조명과 tdi를 구비한 이미지 측정 시스템 | |
De Silva et al. | Assessment of tomographic PIV in wall-bounded turbulence using direct numerical simulation data | |
Lu et al. | Modulation measuring profilometry with cross grating projection and single shot for dynamic 3D shape measurement | |
CN105579809B (zh) | 测量方法、测量装置以及计算机可读记录介质 | |
CN109341548B (zh) | 基于变密度条纹的三维振动视觉测量***和方法 | |
JP3629532B2 (ja) | 連続移動物体のリアルタイム形状計測方法及びシステム | |
CN108050955B (zh) | 基于结构光投影与数字图像相关的高温空气扰动滤除方法 | |
CN105738073B (zh) | 一种在空间频率域进行像素响应函数测量的方法 | |
CN114577444A (zh) | 一种调制传递函数的检测方法及检测装置 | |
Kozul et al. | A scanning particle tracking velocimetry technique for high-Reynolds number turbulent flows | |
JP2013178174A (ja) | 複数の格子を用いた三次元形状計測装置 | |
US20200370947A1 (en) | Devices and methods for vibration analysis | |
Shang et al. | Analysis and reduction of error caused by tested object using fringe projection technique with wavelet transform | |
CN110458785A (zh) | 一种基于图像传感的磁悬浮球悬浮间隙检测方法 | |
Chen et al. | High-speed Fourier transform profilometry for reconstructing objects having arbitrary surface colours | |
JP6013819B2 (ja) | 表面形状検査装置及び表面形状検査方法 | |
Fuchs et al. | Uncertainty quantification of three-dimensional velocimetry techniques for small measurement depths | |
CN105593635B (zh) | 测量方法、测量装置以及计算机可读记录介质 | |
Jin et al. | The online measurement of optical distortion for glass defect based on the grating projection method | |
Sioma et al. | Focus Estimation Methods for Use in Industrial SFF Imaging Systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221208 Address after: 314100 Floor 2, Building 7, No. 1, Weizhong Road, Weitang Street, Jiashan County, Jiaxing City, Zhejiang Province Patentee after: Zhejiang Zhixiang Photoelectric Technology Co.,Ltd. Address before: 710032 No. 2 Xuefu Middle Road, Weiyang District, Xi'an City, Shaanxi Province Patentee before: XI'AN TECHNOLOGICAL University |