CN116915972A - 一种基于刃边像分析的高分辨率相机mtf测试装置及方法 - Google Patents
一种基于刃边像分析的高分辨率相机mtf测试装置及方法 Download PDFInfo
- Publication number
- CN116915972A CN116915972A CN202310750288.3A CN202310750288A CN116915972A CN 116915972 A CN116915972 A CN 116915972A CN 202310750288 A CN202310750288 A CN 202310750288A CN 116915972 A CN116915972 A CN 116915972A
- Authority
- CN
- China
- Prior art keywords
- camera
- edge
- tested
- mtf
- image
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000010191 image analysis Methods 0.000 title claims abstract description 26
- 238000004458 analytical method Methods 0.000 claims abstract description 42
- 238000013519 translation Methods 0.000 claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 21
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 238000010998 test method Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000003708 edge detection Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 8
- 238000009499 grossing Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 210000001747 pupil Anatomy 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 230000004075 alteration Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000005304 optical glass Substances 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 abstract description 6
- 238000012634 optical imaging Methods 0.000 abstract description 2
- 238000012372 quality testing Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 230000001276 controlling effect Effects 0.000 description 5
- 230000000737 periodic effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012546 transfer Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010561 standard procedure Methods 0.000 description 2
- 238000012356 Product development Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 229910052724 xenon Inorganic materials 0.000 description 1
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N17/00—Diagnosis, testing or measuring for television systems or their details
- H04N17/002—Diagnosis, testing or measuring for television systems or their details for television cameras
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及光学成像质量测试装置及方法，具体涉及一种基于刃边像分析的高分辨率相机MTF测试装置及方法，解决了现有MTF测试的精度低、测试重复性差的技术问题。本发明提供的MTF测试装置在待测试相机焦面位置安装刀口靶和平行光管模拟标准刃边目标，采用电控旋转台调整电控旋转台刃边的倾斜角度消除非均匀采样引入的测试误差，采用调焦平移台模拟不同物距的刃边目标满足不同物距成像条件下的MTF测试需求。本发明还提供一种MTF测试方法，通过减去背景图像、选择分析窗口、调整刃边倾斜角度、刃边检测、滤波消除噪声、校正运算误差，提高了MTF测试的精度和稳定性。
Description
技术领域
本发明涉及光学成像质量测试装置及方法，具体涉及一种基于刃边像分析的高分辨率相机MTF测试装置及方法。
背景技术
相机的调制传递函数(MTF)反映了相机对目标物不同空间频率信息的传递能力，能全面、客观地表征相机的成像性能，是评价相机成像质量的重要手段，是相机技术协议签订、研制过程质量控制和产品交付验收最常用的技术参数之一。相机MTF的测试方法主要有：周期矩形条纹法、点源法、狭缝法和刀口法。
其中，周期矩形条纹法具有条纹图像对比度，能够直观反应相机的成像质量、数据处理过程简单的优点，应用最为广泛。但是其缺点也非常明显：(1)测试靶标无法通用，需要根据被测相机定制周期矩形条纹靶；(2)靶标精度要求高，需要达到微米级；(3)一组周期矩形条纹只能测试单一频率下的MTF，测试效率低，一般只用于测试相机奈奎斯特频率下的MTF；(4)测试过程对调整精度要求高，需要使矩形条纹像和探测器相位严格对准。周期矩形条纹板法存在测试调整精度要求高、经济性和测试效率低的问题。
狭缝法能够同时给出相机不同频率下的MTF，但狭缝法对狭缝的宽度、厚度、两边直线度、两边平行度要求较高，同样需要根据被测相机定制合适宽度的狭缝，且测试结果要扣除狭缝宽度的影响，数据处理复杂，还存在欠采样的问题。狭缝法优点不突出，缺点较明显，目前很少用于相机MTF测试。
星点法测试相机MTF的原理为：星点像代表相机的脉冲响应函数，其傅里叶变换即为相机的调制传递函数。星点像光强分布规律决定了相机成像的清晰程度及成像质量，其缺点是信号较弱、易受噪声影响、采样点少、存在欠采样问题，重复性及精度较差。实际上星点总是有一定的物理尺寸的，需要扣除星点像自身对MTF的影响。
相比其他方法，刃边法优势明显，主要表现在：(1)评价更全面，测试效率高，一次测试能够同时得到全频段的MTF；(2)测量成本低，一个靶标能满足不同相机MTF的测试需求；(3)只有一条刃边，靶标加工精度相对较低。以上优点使得刃边法近些年来在民用数码相机、军用摄像机的MTF测试中应用日益广泛。刃边法在国内外被广泛用于星载光学遥感相机、机载光学遥感相机的外场性能鉴定和评价，也被用于阵列式图像传感器(COMS、CCD)的实验室性能测试。
刃边法被ISO国际标准化组织采纳，作为测试相机MTF的标准方法。ISO标准中虽然给出采用刃边法计算MTF的关键过程和步骤，但是该算法相对复杂，影响计算结果的因素较多，如：探测器噪声、响应线性度、均匀性、分析窗口大小、刃边斜率大小、刃边斜率误差、算法自身的参数，而标准并没有考虑这些因素的影响和明确的处理方法，使得在利用该标准方法进行MTF测试时精度低、测试重复性差。这对于产品研制过程质量控制、产品验收、性能鉴定以及科学研究都是不能接受的，因此，必须进一步深入研究，确定解决这些问题办法，制定抑制这些影响因素的措施，明确算法参数，提高刃边法测量MTF精度和重复性。
发明内容
本发明的目的是解决现有MTF测试的精度低、测试重复性差的技术问题，而提供一种基于刃边像分析的高分辨率相机MTF测试装置及方法。
为实现上述目的，本发明采用的技术方案为：
一种基于刃边像分析的高分辨率相机MTF测试装置，其特殊之处在于，包括平行光管、调焦平移台、安装在调焦平移台上的电控旋转台和均匀光源、安装在电控旋转台上的刀口靶，以及数据处理单元；
所述刀口靶、平行光管依次设置在均匀光源出射光的光路上，所述调焦平移台用于承载电控旋转台和均匀光源，并在调节成像目标的物距时沿平行光管的光轴移动，所述电控旋转台用于调整刀口靶的角度；
所述平行光管的出射光路上设置待测试相机，使得平行光管的出射光束充满待测试相机的入瞳；
所述数据处理单元分别与电控旋转台、调焦平移台、待测试相机电连接，用于控制调焦平移台的位移、电控旋转台的旋转、接收待测试相机的图像信息并对图像信息进行分析和计算得到待测试相机的MTF。
进一步地，所述电控旋转台上设置有中心孔，所述刀口靶安装在中心孔内；
所述刀口靶为镀铬刃边图案的光学玻璃平板，所述刃边图案亮暗两侧的对比度大于等于500:1，刃边直线度小于等于1μm；
所述电控旋转台的旋转角度定位误差小于等于0.003°。
进一步地，所述平行光管出射光束口径大于待测试相机的入瞳直径，工作波段覆盖待测试相机的工作谱段，波像差RMS≤λ/14，λ＝623.8nm；
所述平行光管为折射式平行光管、折反式平行光管或全反式平行光管。
进一步地，其特征在于，所述均匀光源的出射光束直径大于刀口靶的直径，均匀光源包括光源和积分球，所述光源设置在积分球内，所述刀口靶设置在积分球的出射光路上。
一种基于刃边像分析的高分辨率相机MTF测试方法，其特殊之处在于，包括以下步骤：
步骤1、检查、调整上述的基于刃边像分析的高分辨率相机MTF测试装置，并使刀口靶成像在待测试相机的中心视场；
步骤2、关闭均匀光源，待测试相机连续采集T帧背景图像Ib(i,j,t)，并将其发送至数据处理单元，按照以下公式计算平均背景图像
其中，i为待测试相机的图像像元行序号，i取1,2,3,…,m，j为待测试相机的图像像元列序号，j取1,2,3,…,n，t为图像帧序号，t取1,2,3,…,T；待测试相机的分辨率为n×m，n、m均为正整数；
步骤3、打开均匀光源电源，待测试相机连续采集T帧刀口图像Is(i,j,t)，并将其发送至数据处理单元，按照以下公式计算T帧减去背景图像后的刀口图像I′s(i,j,t)：
步骤4、数据处理单元对减去背景图像后的刀口图像进行分析计算，得到待测试相机中心视场的MTF；
步骤5、调整待测试相机的姿态，使刀口靶成像在待测试相机的其他视场，按照步骤2-步骤4的方法，获得待测试相机其他视场的MTF。
进一步地，所述步骤4具体为：
4.1、在其中一帧减去背景图像后的刀口图像中选择分析窗口，所述分析窗口的长度包含完整的线扩展函数曲线，宽度至少包含6个相位周期；
4.2、对所述分析窗口区域的图像进行中值滤波；
4.3、用canny边缘检测算法检测经中值滤波后的分析窗口区域图像的刃边，利用最大类间方差法确定canny算子的阈值，获得边缘检测图像；
4.4、根据步骤4.3得到的边缘检测图像，获取刃边位置对应的坐标(xa,ya)，对刃边位置的坐标进行最小二乘直线拟合，获得拟合函数y＝kx+b，进而得到刃边斜率k和截距b；(xa,ya)表示第a个刃边位置的坐标；
4.5、根据刃边斜率对分析窗口的宽度进行裁剪，裁剪后的分析窗口的行数为：
其中，u'为分析窗口裁剪后的行数，u为分析窗口原始行数，k为刃边斜率，||为取绝对值运算，round为四舍五入取整函数，floor为向下取整函数；
4.6、将裁剪后的分析窗口区域图像的每一个像元图像沿刃边进行映射，将u'行v列像元的响应值映射为一行u'×v列响应值，得到原始边缘扩散函数ESF(xl)，其中，xl为原始边缘扩散函数ESF(xl)第l个点的横坐标，l取1,2,3,…,u'×v；
4.7、划分合并区间长度，对原始边缘扩散函数进行合并，获得合并后的边缘扩展函数ESF′(p)，其中，p为第p个合并区间；
4.8、对合并后的边缘扩展函数进行SG平滑滤波，并对滤波后的边缘扩展函数进行隔点差分计算，得到线扩展函数LSF(p)：
LSF(p)＝[ESF′(p+1)-ESF′(p-1)]×0.5；
4.9、对线扩展函数进行汉宁窗滤波，并对滤波后的线扩展函数进行离散傅里叶变换、取模并进行归一化，得到MTF(p)：
MTF(p)＝|FFT(LSF(p))|/|FFT(LSF(p))|p＝0
其中，FFT()表示快速傅里叶变换运算；
4.10、对隔点差分计算引入的误差进行校正，得到由一帧图像计算得到的被测相机中心视场的MTF，具体为：
其中，MTF′(p)为由一帧图像计算得到的被测相机中心视场的MTF；
4.11、重复步骤4.1-步骤4.10，得到分别由其余T-1帧图像计算得到的被测相机中心视场的MTF，取分别由T帧图像计算得到的被测相机中心视场的MTF的平均值，得到待测试相机中心视场的MTF。
进一步地，所述步骤1具体为：
1.1、检查平行光管，调整调焦平移台的位置，使刀口靶刻划面位于平行光管的焦面位置，数据处理单元获取并记录调焦平移台的位置；
1.2、根据待测试相机的工作物距或测试要求的物距LO，按照以下公式计算平行光管1的离焦量：
其中，△x为平行光管模拟物距为LO的目标对应的离焦量，fc为平行光管的焦距；
1.3、根据步骤1.2获得的平行光管的离焦量，数据处理单元控制调焦平移台向靠近平行光管的方向移动，调焦平移台的位移量为△x；
1.4、调整待测试相机的位置及姿态，使其对准平行光管，并且刀口靶成像于待测试相机的中心视场；
1.5、将待测试相机的积分时间和增益分别设置为MTF测试要求的积分时间和增益；
1.6、调整均匀光源的亮度使待测试相机得到的刀口像的亮侧灰度均值位于待测试相机的线性工作区；
1.7、调整电控旋转台，使刀口像刃边与待测试相机探测器列方向夹角θ满足：cot(θ)取4、5、6、7、8、9、10、11、12、13、14、15或16。
进一步地，步骤4.7中，所述合并区间长度为0.25像元，所述合并后的边缘扩展函数ESF′(p)为：
其中，rect()为矩形函数，当|xl-p·0.25|≤0.125时，rect(xl-p·0.25)＝1，否则rect(xl-p·0.25)＝0；
N为使rect(xl-p·0.25)＝1的xl的个数；p为第p个合并区间；
ESF′(p)表示xl落在第p个合并区间内所有点的ESF(xl)的均值。
进一步地，步骤4.8中，所述SG平滑滤波选用7点3阶多项式。
进一步地，步骤1.6中，所述待测试相机的线性工作区为待测试相机饱和灰度值的20％～80％。
与现有技术相比，本发明具有的有益技术效果如下：
1.本发明提供了一种基于刃边像分析的高分辨率全色相机MTF测试装置，该装置用焦面位置安装刀口靶的平行光管模拟标准刃边目标，刀口靶安装在电控旋转台上，可以精确调整刃边的倾斜角度，消除非均匀采样引入的测试误差。刀口靶安装在调焦平移台上，可以沿平行光管光轴前后移动，模拟不同物距的刃边目标，满足不同物距成像条件下的MTF测试需求。
2.本发明提供了一种基于刃边像分析的高分辨率全色相机MTF测试方法，该方法通过减去背景图像、合理选择分析窗口大小、精确控制调整刃边倾斜角度、对分析窗口图像中值滤波、用最大类间方差法确定阈值和canny算法检测刃边、对边缘扩展函数SG滤波、校正差分运算误差等关键步骤，提高了MTF测试的精度和稳定性。
附图说明
图1为本发明提供的基于刃边像分析的高分辨率相机MTF测试装置的结构原理图；
图2为本发明实施例中刀口靶调整角度示意图；
图3为本发明实施例中选择的分析窗口示意图；
图4为本发明实施例中获得的线扩展函数LSF曲线图；
图5为本发明实施例中获得的边缘检测图像示意图；
图6为本发明实施例中获得的边缘扩展函数曲线图；
图7为本发明实施例中获得的待测试相机中心视场的MTF曲线图；
附图标记说明如下：
1-平行光管，2-刀口靶，3-电控旋转台，4-调焦平移台，5-均匀光源，6-数据处理单元，7-待测试相机。
具体实施方式
下面结合附图和具体实施方式对本发明提出的一种基于刃边像分析的高分辨率相机MTF测试装置及方法作进一步详细说明。本领域技术人员应当理解的是，这些实施方式仅仅用来解释本发明的技术原理，目的并不是用来限制本发明的保护范围。
一种基于刃边像分析的高分辨率相机MTF测试装置，如图1所示，包括平行光管1、调焦平移台4、安装在调焦平移台4上的电控旋转台3和均匀光源5、安装在电控旋转台3上的刀口靶2，以及数据处理单元6。
刀口靶2、平行光管1依次设置在均匀光源5出射光的光路上，调焦平移台用于承载电控旋转台3和均匀光源5，并在调节成像目标的物距时沿平行光管1的光轴移动。调焦平移台4可沿平行光管1的光轴前后移动，实现不同物距刀口目标的模拟，适应不同成像物距待测试相机7的MTF的测试需求。
电控旋转台3上设置有中心孔，刀口靶2安装在中心孔内并位于均匀光源5的出射光路上，电控旋转台3用于调整刀口靶2上刃边图案的倾斜角度。刀口靶2为镀铬刃边图案的光学玻璃平板上，刃边图案亮暗两侧的对比度大于等于500:1，刃边直线度小于等于1μm。电控旋转台3的旋转角度定位误差小于等于0.003°。
平行光管1的出射光路上设置待测试相机7，使得平行光管1的出射光束充满待测试相机7的入瞳。平行光管1出射光束口径大于待测试相机7的入瞳直径，工作波段覆盖待测试相机7的工作谱段，波像差RMS≤λ/14，λ＝623.8nm。平行光管1为折射式平行光管、折反式平行光管或全反式平行光管。
均匀光源5的出射光束直径大于刀口靶2的直径，均匀光源5包括光源和积分球，光源设置在积分球的入射光路上，刀口靶2设置在积分球的出射光路上。光源选择具有与太阳相似发光谱线的卤素灯和氙灯组合。
数据处理单元6分别与电控旋转台3、调焦平移台4、待测试相机7电连接，用于控制调焦平移台4的位移、电控旋转台3的旋转、接收待测试相机7的图像信息并对图像信息进行分析和计算得到待测试相机7的MTF。
数据处理单元6装载有电控旋转台3和调焦平移台4的测控软件，可根据测试需求调整电控旋转台3和调焦平移台4，并采集反馈其状态信息。数据处理单元6还装载有待测试相机7的图像采集卡和软件，以及分析刃边像计算相机MTF的程序，可接收待测试相机7的图像信息，对图像信息进行分析和计算得到待测试相机7的MTF。
本实施例提供的MTF测试装置用在焦面位置安装刀口靶2的平行光管1模拟标准刃边目标，刀口靶2安装在电控旋转台3上，可以精确调整刃边的倾斜角度，消除非均匀采样引入的测试误差。刀口靶2安装在调焦平移台上，可以沿平行光管1光轴前后移动，模拟不同物距的刃边目标，满足不同物距成像条件下的MTF测试需求。
本实施例还提供一种基于刃边像分析的高分辨率相机MTF测试方法，包括以下步骤：
步骤1、检查、调整上述的基于刃边像分析的高分辨率相机MTF测试装置，并使刀口靶2成像在待测试相机7的中心视场，具体为：
1.1、检查平行光管1，调整调焦平移台4的位置，使刀口靶2刻划面位于平行光管1的焦面位置，此时平行光管1模拟无穷远的刃边目标，数据处理单元6获取并记录此时调焦平移台4的位置。
1.2、根据待测试相机7的工作物距或测试要求的物距LO，按照以下公式计算平行光管1的离焦量：
其中，△x为平行光管1模拟物距为LO的目标对应的离焦量，fc为平行光管1的焦距。
1.3、根据步骤1.2获得的平行光管1的离焦量，数据处理单元6控制调焦平移台4向靠近平行光管1的方向移动，调焦平移台4的位移量为△x，此时平行光管1模拟物距为LO的刃边目标。
1.4、调整待测试相机7的位置及姿态，使其对准平行光管1，并且刀口靶2成像在待测试相机7的中心视场。
1.5、将待测试相机7的积分时间和增益分别设置为MTF测试要求的积分时间和增益。
1.6、调整均匀光源5的亮度使待测试相机7得到的刀口像的亮侧灰度均值位于待测试相机7的线性工作区；本实施例中，使刀口像的亮侧灰度均值为测试相机7饱和灰度值的80％。待测试相机7的线性工作区为待测试相机7饱和灰度值的20％～80％，以消除待测试相机7响应非线性引入的测试误差。
1.7、调整电控旋转台3使刀口像刃边与待测试相机7探测器列方向夹角θ满足：cot(θ)取4、5、6、7、8、9、10、11、12、13、14、15或16，以消除非均匀采样引入的测试误差。如图2所示，为夹角θ的示意图，夹角θ的调整误差≤0.003°。
步骤2、关闭均匀光源5电源，待测试相机7连续采集10帧背景图像Ib(i,j,t)，并将其发送至数据处理单元6，按照以下公式计算平均背景图像
其中，i为待测试相机7的图像像元行序号，i取1,2,3,…,m，j为待测试相机7的图像像元列序号，j取1,2,3,…,n，t为图像帧序号，t取1,2,3,…,10。
步骤3、打开均匀光源5电源，待测试相机7连续采集10帧刀口图像Is(i,j,t)，并将其发送至数据处理单元6，按照以下公式计算10帧减去背景图像后的刀口图像Is′(i,j,t)：
步骤4、数据处理单元6对减去背景图像后的刀口图像进行分析计算，得到待测试相机7中心视场的MTF，具体为：
4.1、在其中一帧减去背景图像后的刀口图像中选择如图3所示的分析窗口，分析窗口的长度包含完整的LSF(线扩展函数)曲线，宽度至少包含6个相位周期，即窗口宽度应≥6cot(θ)行，LSF曲线如图4所示。分析窗口的选择与LSF曲线的计算是相互迭代的过程。根据初步选择的分析窗口计算LSF曲线，若获得的LSF曲线完整，即LSF曲线的两端平直并趋于零，则确定此时的分析窗口为最终选定的分析窗口；若获得的LSF曲线不完整，则根据LSF曲线重新选择分析窗口，直至获得完整的LSF曲线。
4.2、对分析窗口区域的图像进行中值滤波，以消除脉冲噪声的影响，中值滤波模板取3×3。
4.3、用canny边缘检测算法检测经中值滤波后的分析窗口区域图像的刃边，利用最大类间方差法(otsu算法)确定canny算子阈值，获得如图5所示的边缘检测图像。
4.4、根据步骤4.3得到的边缘检测图像，获取刃边位置对应的坐标(xa,ya)，对刃边位置的坐标进行最小二乘直线拟合，获得拟合函数y＝kx+b，进而得到刃边斜率k和截距b，(xa,ya)表示第a个刃边位置的坐标。刃边位置对应的坐标通过查找边缘检测图像的各行最大灰度值所在列坐标获得。
4.5、根据刃边斜率对分析窗口的宽度进行裁剪，裁剪后的分析窗口的行数保留整数个相位周期，具体为：
其中，u'为分析窗口裁剪后的行数，u为分析窗口原始行数，k为刃边斜率，||为取绝对值运算，round为四舍五入取整函数，floor为向下取整函数。
4.6、将裁剪后的分析窗口区域图像的每一个像元图像沿刃边进行映射，将u'行v列像元的响应值映射为一行u'×v列响应值，对所用像元映射后，得到原始边缘扩散函数ESF(xl)，其中，xl为原始边缘扩散函数ESF(xl)第l个点的横坐标，l取1,2,3,…,u'×v。映射方法具体为：映射后的横坐标为该像元到刃边的距离，映射后的纵坐标为该像元的灰度值。如坐标为(i,j)灰度值为y的像素，映射后的横坐标为：映射后的纵坐标为y。
4.7、划分合并区间长度为0.25像元，对原始边缘扩散函数进行合并，获得如图6所示的合并后的边缘扩展函数ESF′(p)：
其中，rect()为矩形函数，当|xl-p·0.25|≤0.125时，rect(xl-p·0.25)＝1，否则rect(xl-p·0.25)＝0；
N为使rect(xl-p·0.25)＝1的xl的个数；p为第p个合并区间；
ESF′(p)为xl落在第p个合并区间内所有点的ESF(xl)的均值。
通过对原始边缘扩散函数合并，一方面可以提高边缘扩散的信噪比，抑制噪声的影响；另一方面可以抑制原始边缘扩展函数采样点非均匀性对MTF测试的影响。
4.8、对合并后的边缘扩展函数进行SG平滑滤波，并对滤波后的边缘扩展函数进行隔点差分计算，得到线扩展函数LSF(p)：
LSF(p)＝[ESF′(p+1)-ESF′(p-1)]×0.5；
SG平滑滤波选取7点3阶多项式平滑滤波，可进一步抑制随机噪声对MTF测试结果的影响。
4.9、对线扩展函数进行汉宁窗滤波，防止频率泄漏，并对滤波后的线扩展函数进行离散傅里叶变换、取模并进行归一化，得到MTF(p)：
MTF(p)＝|FFT(LSF(p))|/|FFT(LSF(p))|p＝0
其中，FFT()表示快速傅里叶变换运算。
4.10、对隔点差分计算引入的误差进行校正，得到由一帧图像计算得到的被测相机中心视场的MTF，具体为：
其中，MTF′(p)为由一帧图像计算得到的被测相机中心视场的MTF；
4.11、重复步骤4.1-步骤4.10，得到分别由其余9帧图像计算得到的被测相机中心视场的MTF，取分别由10帧图像计算得到的被测相机中心视场的MTF的平均值，得到如图7所示的待测试相机7中心视场的MTF。
步骤5、调整待测试相机7的姿态，使刀口靶2成像在其他视场，重复步骤2-步骤4，获得待测试相机7其他视场的MTF。
本实施例提供的MTF测试方法通过减去背景图像、合理选择分析窗口大小、精确控制调整刃边倾斜角度、对分析窗口图像中值滤波、用otsu算法确定阈值和canny算法检测刃边、对边缘扩展函数SG滤波、校正差分运算误差等关键步骤，提高了MTF测试的精度和稳定性。
Claims (10)
1.一种基于刃边像分析的高分辨率相机MTF测试装置，其特征在于：包括平行光管(1)、调焦平移台(4)、安装在调焦平移台(4)上的电控旋转台(3)和均匀光源(5)、安装在电控旋转台(3)上的刀口靶(2)，以及数据处理单元(6)；
所述刀口靶(2)、平行光管(1)依次设置在均匀光源(5)出射光的光路上；所述调焦平移台(4)用于在调节成像目标的物距时沿平行光管(1)的光轴移动，所述电控旋转台(3)用于调整刀口靶(2)的角度；
所述平行光管(1)的出射光路上设置待测试相机(7)，使得平行光管(1)的出射光束充满待测试相机(7)的入瞳；
所述数据处理单元(6)分别与电控旋转台(3)、调焦平移台(4)、待测试相机(7)电连接，用于控制调焦平移台(4)的位移、电控旋转台(3)的旋转、接收待测试相机(7)的图像信息并对图像信息进行分析和计算得到待测试相机(7)的MTF。
2.根据权利要求1所述的基于刃边像分析的高分辨率相机MTF测试装置，其特征在于：所述电控旋转台(3)上设置有中心孔，所述刀口靶(2)安装在中心孔内；
所述刀口靶(2)为镀铬刃边图案的光学玻璃平板，所述刃边图案亮暗两侧的对比度大于等于500:1，刃边直线度小于等于1μm；
所述电控旋转台(3)的旋转角度定位误差小于等于0.003°。
3.根据权利要求2所述的基于刃边像分析的高分辨率相机MTF测试装置，其特征在于：所述平行光管(1)出射光束口径大于待测试相机(7)的入瞳直径，工作波段覆盖待测试相机(7)的工作谱段，波像差RMS≤λ/14，λ＝632.8nm；
所述平行光管(1)为折射式平行光管、折反式平行光管或全反式平行光管。
4.根据权利要求1-3任一所述的基于刃边像分析的高分辨率相机MTF测试装置，其特征在于：所述均匀光源(5)的出射光束直径大于刀口靶(2)的直径，均匀光源(5)包括光源和积分球，所述光源设置在积分球内，所述刀口靶(2)设置在积分球的出射光路上。
5.一种基于刃边像分析的高分辨率相机MTF测试方法，其特征在于，包括以下步骤：
步骤1、检查、调整权利要求1-4任一所述的基于刃边像分析的高分辨率相机MTF测试装置，并使刀口靶(2)成像在待测试相机(7)的中心视场；
步骤2、关闭均匀光源(5)，待测试相机(7)连续采集T帧背景图像Ib(i,j,t)，并将其发送至数据处理单元(6)，按照以下公式计算平均背景图像
其中，i为待测试相机(7)的图像像元行序号，i取1,2,3,…,m，j为待测试相机(7)的图像像元列序号，j取1,2,3,…,n，t为图像帧序号，t取1,2,3,…,T；
步骤3、打开均匀光源(5)，待测试相机(7)连续采集T帧刀口图像Is(i,j,t)，并将其发送至数据处理单元(6)，按照以下公式计算T帧减去背景图像后的刀口图像I′s(i,j,t)：
步骤4、数据处理单元(6)对减去背景图像后的刀口图像进行分析计算，得到待测试相机(7)中心视场的MTF；
步骤5、调整待测试相机(7)的姿态，使刀口靶(2)成像在待测试相机(7)的其他视场，按照步骤2-步骤4的方法，获得待测试相机(7)其他视场的MTF。
6.根据权利要求5所述的基于刃边像分析的高分辨率相机MTF测试方法，其特征在于，所述步骤4具体为：
4.1、在其中一帧减去背景图像后的刀口图像中选择分析窗口，所述分析窗口的长度包含完整的线扩展函数曲线，宽度至少包含6个相位周期；
4.2、对所述分析窗口区域的图像进行中值滤波；
4.3、用canny边缘检测算法检测经中值滤波后的分析窗口区域图像的刃边，利用最大类间方差法确定canny算子的阈值，获得边缘检测图像；
4.4、根据步骤4.3得到的边缘检测图像，获取刃边位置对应的坐标(xa,ya)，对刃边位置的坐标进行最小二乘直线拟合，获得拟合函数y＝kx+b，进而得到刃边斜率k和截距b；
4.5、根据刃边斜率对分析窗口的宽度进行裁剪，裁剪后的分析窗口的行数为：
其中，u'为分析窗口裁剪后的行数，u为分析窗口原始行数，k为刃边斜率，||为取绝对值运算，round为四舍五入取整函数，floor为向下取整函数；
4.6、将裁剪后的分析窗口区域图像的每一个像元图像沿刃边进行映射，将u'行v列像元的响应值映射为一行u'×v列响应值，得到原始边缘扩散函数ESF(xl)，其中，xl为原始边缘扩散函数ESF(xl)第l个点的横坐标，l取1,2,3,…,u'×v；
4.7、划分合并区间长度，对原始边缘扩散函数进行合并，获得合并后的边缘扩展函数ESF′(p)，其中，p为第p个合并区间；
4.8、对合并后的边缘扩展函数进行SG平滑滤波，并对滤波后的边缘扩展函数进行隔点差分计算，得到线扩展函数LSF(p)：
LSF(p)＝[ESF′(p+1)-ESF′(p-1)]×0.5；
4.9、对线扩展函数进行汉宁窗滤波，并对滤波后的线扩展函数进行离散傅里叶变换、取模并进行归一化，得到MTF(p)：
MTF(p)＝|FFT(LSF(p))|/FFT(LSF(p))p＝0
其中，FFT()表示快速傅里叶变换运算；
4.10、对隔点差分计算引入的误差进行校正，得到由一帧图像计算得到的被测相机中心视场的MTF，具体为：
其中，MTF′(p)为由一帧图像计算得到的被测相机中心视场的MTF；
4.11、重复步骤4.1-步骤4.10，得到分别由其余T-1帧图像计算得到的被测相机(7)中心视场的MTF，取分别由T帧图像计算得到的被测相机中心视场的MTF的平均值，得到待测试相机(7)中心视场的MTF。
7.根据权利要求6所述的基于刃边像分析的高分辨率相机MTF测试方法，其特征在于，所述步骤1具体为：
1.1、检查平行光管(1)，调整调焦平移台(4)的位置，使刀口靶(2)刻划面位于平行光管(1)的焦面位置，数据处理单元(6)获取并记录调焦平移台(4)的位置；
1.2、根据待测试相机(7)的工作物距或测试要求的物距LO，按照以下公式计算平行光管(1)的离焦量：
其中，△x为平行光管(1)模拟物距为LO的目标对应的离焦量，fc为平行光管(1)的焦距；
1.3、根据步骤1.2获得的平行光管(1)的离焦量，数据处理单元(6)控制调焦平移台(4)向靠近平行光管(1)的方向移动，调焦平移台(4)的位移量为△x；
1.4、调整待测试相机(7)的位置及姿态，使其对准平行光管(1)，并且刀口靶(2)成像于待测试相机(7)的中心视场；
1.5、将待测试相机(7)的积分时间和增益分别设置为MTF测试要求的积分时间和增益；
1.6、调整均匀光源(5)的亮度使待测试相机(7)得到的刀口像的亮侧灰度均值位于待测试相机(7)的线性工作区；
1.7、调整电控旋转台(3)，使刀口像刃边与待测试相机(7)探测器列方向夹角θ满足：cot(θ)取4、5、6、7、8、9、10、11、12、13、14、15或16。
8.根据权利要求6-7任一所述的基于刃边像分析的高分辨率相机MTF测试方法，其特征在于，步骤4.7中，所述合并区间长度为0.25像元，所述合并后的边缘扩展函数ESF′(p)为：
其中，rect()为矩形函数，当|xl-p·0.25|≤0.125时，rect(xl-p·0.25)＝1，否则rect(xl-p·0.25)＝0；
N为使rect(xl-p·0.25)＝1的xl的个数；p为第p个合并区间；
ESF′(p)表示xl落在第p个合并区间内所有点的ESF(xl)的均值。
9.根据权利要求8所述的基于刃边像分析的高分辨率相机MTF测试方法，其特征在于：步骤4.8中，所述SG平滑滤波选用7点3阶多项式。
10.根据权利要求9所述的基于刃边像分析的高分辨率相机MTF测试方法，其特征在于：步骤1.6中，所述待测试相机(7)的线性工作区为待测试相机(7)饱和灰度值的20％～80％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310750288.3A CN116915972A (zh) | 2023-06-25 | 2023-06-25 | 一种基于刃边像分析的高分辨率相机mtf测试装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310750288.3A CN116915972A (zh) | 2023-06-25 | 2023-06-25 | 一种基于刃边像分析的高分辨率相机mtf测试装置及方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116915972A true CN116915972A (zh) | 2023-10-20 |
Family
ID=88354027
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310750288.3A Pending CN116915972A (zh) | 2023-06-25 | 2023-06-25 | 一种基于刃边像分析的高分辨率相机mtf测试装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116915972A (zh) |
-
2023
- 2023-06-25 CN CN202310750288.3A patent/CN116915972A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105486489B (zh) | 电视成像系统调制传递函数测试装置及方法 | |
| CN104062098B (zh) | 一种双线阵ccd扫描成像测量光束质量的装置及方法 | |
| CN100395538C (zh) | 一种新型快速椭圆偏振光测量系统 | |
| CN107403177A (zh) | 基于工业相机的亮度测量方法 | |
| CN109767425B (zh) | 机器视觉光源均匀性评估装置及方法 | |
| US8059282B2 (en) | Reflective film thickness measurement method | |
| CN101706951B (zh) | 一种基于特征融合的气动光学图像质量客观评价方法、装置及系统 | |
| CN107144420B (zh) | 光学镜头像差检测装置及方法 | |
| CN109990985A (zh) | 一种品字型线列红外探测器调制传递函数测试方法 | |
| CN101726316A (zh) | 内方位元素及畸变测试仪 | |
| US7518712B2 (en) | Tilted edge for optical-transfer-function measurement | |
| EP2577265B1 (en) | Apparatus and method for locating the centre of a beam profile | |
| CN117928875A (zh) | 一种风洞流场的时间分辨偏振成像装置及方法 | |
| CN105911716A (zh) | 一种传函测试中被测镜头的光路调节装置及其调节方法 | |
| CN210571299U (zh) | 用于测量小视场投影模组光学参数的系统 | |
| CN109655158B (zh) | 基于大气廓线与led的高光谱遥感器在轨光谱定标的方法 | |
| EP1896792B1 (en) | Method and system for measuring the curvature of an optical surface | |
| US3814520A (en) | Method of, and apparatus for gauging, inspecting or measuring physical properties of objects | |
| CN201277864Y (zh) | 内方位元素及畸变测试仪 | |
| CN116915972A (zh) | 一种基于刃边像分析的高分辨率相机mtf测试装置及方法 | |
| Zhang et al. | Modified slanted-edge method to measure the modulation transfer function of camera | |
| CN114235347B (zh) | 镜头质量评估方法及装置 | |
| CN116105875A (zh) | 一种基于背景辐射校正的偏振遥感器短波红外波段偏振定标方法 | |
| CN106969707B (zh) | 一种检测装置及其控制方法 | |
| CN117041521A (zh) | 一种Bayer滤波型彩色相机的MTF测试装置及测试方法 |
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 |