电视成像系统调制传递函数测试装置及方法
电视成像系统调制传递函数测试装置及方法
【技术领域】
[0001] 本发明设及光学计量与测量领域，具体为一种电视成像系统调制传递函数测试装 置及方法。
【背景技术】
[0002] 调制传递函数(MTF)是成像系统性能最基本的评价参数之一，反映了成像系统对 景物中各空间频率分量的响应特性，体现了成像系统对不同空间频率的目标对比度的传输 能力，因此MTF是评价成像系统成像性能的一项重要指标。
[0003] 可见光电视成像观瞄系统负责对目标进行捜索、跟踪和瞄准，其成像性能是决定 系统作用距离和跟踪精度的关键所在，因此必须对电视成像系统的成像性能进行准确测量 和评估。
[0004] 可见光电视成像观瞄系统通常是将光学系统和CCD成像器件固定在一起，由CCD将 光学系统所成的像转化为电信号，然后对电信号进行电子放大和处理，最终将视频图像信 号传输到显示器上，进行图像显示。由此可见，电视成像系统成像的性能不仅与光学系统成 像性能、CCD成像器件的性能有关，还与光学系统与CCD的匹配精度、电子电路的处理和转化 精度等有关。因此，对于电视成像系统的MTF评价方法既不同于普通光学透镜的像质评价， 也不同于CCD器件的性能评价方法。为了对电视成像系统的成像性能进行评价，就必须对电 视成像系统整机的MTF进行准确测量。
[000引电视成像系统MTF测试装置需要利用待测电视成像系统自身CCD器件采集并获得 祀标的图像，由计算机软件通过对CCD像面祀标的灰度分布进行分析计算得到待测电视成 像系统的MTF值。
[0006] 在实际中，由于电视成像系统MTF测试装置照明光源强度空间分布的不均匀、强度 不稳定、W及电视成像系统自身的特点，如接收器为离散型器件、具有不同灰度动态范围， 器件具有随机噪声等，会影响CCD器件像面图像的灰度分布W及线扩散函数的判读等，最终 影响MTF测量结果。
【发明内容】
[0007] 为了解决W上问题，本发明提供一种电视成像系统调制传递函数测试装置及方 法，可对不同视频格式、不同灰度等级、W及不同噪声水平的各类电视成像系统MTF测量仪 进行高准确度测量。
[0008] 本发明的技术方案为：
[0009] 所述一种电视成像系统调制传递函数测试装置，其特征在于:包括积分球光源系 统（1)、祀标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜巧)和计算机处理系统 (7);
[0010 ]所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统；
[0011]所述积分球光源系统（1)的输出光强可调；
[001引所述祀标组(2)由单狭缝祀标、多狭缝祀标和全黑祀标组成;根据需要选择祀标组 (2)中的祀标固定在积分球光源系统（1)的出口位置;祀标处于离轴抛物面反射镜(5)的焦 点位置；
[0013] 所述滤光片组（3)安装在积分球光源系统（1)的出口位置;将从积分球光源系统 (1)输出的光变为所需波长的光；
[0014] 所述计算机处理系统(7)与积分球光源系统（1)W及待测电视成像系统(8)连接， 待测的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控制积分球 光源系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行图像采集， 计算机处理系统(7)对采集的图像进行处理。
[0015] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 所述积分球光源系统（1)由积分球（1-1)、面鹤灯（1-2)、精密直流程控电源（1-3)和电动可 变光阔（1-4)组成;计算机处理系统（7)控制精密直流程控电源（1-3)和/或电动可变光阔 (1-4)调整积分球光源系统（1)输出光强。
[0016] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 积分球（1-1)出口端面安装有定位销;祀标组(2)的祀标边缘有缺口，祀标边缘缺口与定位 销配合;带有狭缝的祀标上的狭缝方向垂直于待测电视成像系统(8)的CCD祀面扫描方向。
[0017] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 祀标组(2)中，单狭缝祀标的狭缝宽度选取为50μπι、100Μ?或200μπι，单狭缝祀标和多狭缝祀 标的狭缝开口直线性优于0.5%。
[0018] 所述一种利用上述装置进行电视成像系统调制传递函数测试的方法，其特征在 于:包括W下步骤：
[0019] 步骤1:将祀标组(2)中的多狭缝祀标安装在积分球光源系统（1)出口位置，调节待 测电视成像系统(8)，对多狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图像 数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不满 足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝 祀标上的狭缝个数为T;
[0020] 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝祀标采集的图像， 利用该图像并按照W下步骤获得系统放大率β:
[0021 ] 步骤2.1:采用二值化阔值Tbw对图像进行二值化处理，得到二值化图像Ibw;
[0022] 步骤2.2:在二值化图像18冲选取一个矩形区域11?)1，所述矩形区域11?)1高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧；
[0023] 步骤2.3:计算矩形区域Iroi内每一行中灰度值为1的像素个数，形成数据序列 {Nw}：
[0024]
[0025] 并依次计算系统放大率β:
[0026]
[0027] 其中，v，w均为整数，v = a，a+l，a+2, = ，b+l，b+l，···，B; (Xv,yw)为矩形区域 Iroi中像素(v，w)的坐标，（xa，yb)为矩形区域Iroi中的像素起点（a，b)的坐标，（XA，yB)为矩形 区域Iroi中的像素终点(A，B)的坐标;Iroi (xv，yw)为矩形区域中（xv，yw)处的灰度值；Δ X为待 测电视成像系统(8)的像素尺寸，L为多狭缝祀标的单个狭缝宽度;f'为待测电视成像系统 (8)的焦距，时为准直光学系统焦距；
[0028] 步骤3:采用全黑祀标替换多狭缝祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对全黑祀标进行图像采集;计算机处理系统(7)提取待测电视成像 系统(8)对全黑祀标采集的图像，得到背景灰度Ibk;
[0029] 步骤4:采用单狭缝祀标替换全黑祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对单狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图 像数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不 满足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测 电视成像系统(8)对单狭缝祀标再次进行图像采集，计算机处理系统(7)对采集的若干帖图 像进行平均，得到平均灰度图像lAV;
[0030] 步骤5:通过W下步骤测量传递函数：
[0031] 步骤5.1:在平均灰度图像lAv中选取一个矩形区域1/ AV，所述矩形区域1/ AV高度方 向上的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧；
[0032] 步骤5.2:对I%v中的每一行，均采用高斯函数对该行的灰度数据序列化1,^进行最 小二乘法数据拟合，得到拟合后的灰度数据序列{LF 口 1,^，其中表示I%v中第j行的第i 点的灰度值，LFITi,读示I^Av中第j行的第i点的拟合灰度值；由拟合后的灰度数据序列 (LF 口 ij}得到矩形区域I"Av;
[0033] 步骤5.3:捜索矩形区域I"Av中每一行的灰度最大值，得到矩形区域I"Av中每一行 灰度最大值对应像素的坐标，其中矩形区域I"Av中第U行灰度最大值对应像素(max_u，u)的 坐标为(Xmax_u，yu)，得到数据序列{Xmax_u，Yu}:
[0034] Xmax u _ { Xmax_e , Xmax_e+1 , Xmax_e+2 ,......Xmax_E} j
[003引 Yu= {ye，ye+i，ye+2,......ye}；
[0036] 其中u = e，e+l，e+2, ···EiYe为矩形区域I"AV中起始行的行坐标，yE为矩形区域I"AV 中终止行的行坐标；
[0037] 步骤5.4:对于矩形区域I"AV中的第U行，11 = e，e+1，e+2，…E，W该行灰度最大值对 应像素(max_u，u)为中屯、，计算中屯、两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝 对值的均值：
[003引
[0039] 其中亂同时满足条件沁3[/^争諷',,心.>',,）-^枯。,一,"，乂,）]<半、^片纖,,+.^^^ ，凡）<Φ、. Xn?__u+W,, - X。和 - %; Xd为 I" AV 中起始列的列坐柄，XD 为 I" AV 中 终止列的列坐标，Ψ为像素灰度差阔值，Φ为像素灰度最小阔值，Φ取值为背景灰度值Ibk;
[0040] 得到数据序列 Ay={S(ye)，S(ye+i)，S(ye+2)，......，S(yE)};
[0041 ] 步骤5.5 : 捜索数据序列Δ y的所有极小值：4执& )，)，.….，，、）；得 到所有极小值行坐标组成的新的数据序列〇={.Vln,，Λν,，...，丸Kf }，L化~LNf是对应极 小值所在行在矩形区域I"Av中的行号;f是极小值个数;在数据序列Ω中任选两个相邻数据 LNp和LNp+i，得到倾斜角日为：
[0042] "二 arctan[A、'乃少.|、 - 、；)];
[0043] 步骤5.6:对矩形区域I"Av中，处于第LNp行和第LNp+l行之间的像素的灰度数据进行 合并，形成数据序列化S(q)}:
[0044] 当 α〉〇时，LS(q)为：
[004引其中
，into为取整函数，modO为取余函数;q = l，2，…，RNX(D-d)，RN=LNp+广LNp;d为I"A冲起始列的像素编号，D为I"A冲终止列的像素 编号；
[0049] 步骤5.7:采用W下步骤对数据序列化S(q)}进行采样，得到数据序列化SF(t)}:
[0050] 步骤5.7.1:计算采样间隔INTI:
[0051]
[0052] 式中into为取整函数;S为采样率；
[00閲步骤5.7.2:根据采样间隔1饥'1计算版)0。:1孤。=1110(1(9/1饥'1)，1110(10为取余函数；
[0054]步骤5.7.3:在数据序列化S(q)}捜索MODq为0的数据，组成数据序列化SF(t)}，数 据序列化SF(t)}中的数据个数为INT2;
[005引步骤5.8:根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):
[0065] Κ = INT2; 1为单狭缝祀标的狭缝宽度。
[0066] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤及步骤4中比值的设定要求为30 %~80 %。
[0067] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤4中，计算机处理系统(7)对采集的不少于10帖图像进行平均，得到平均灰度图像lAv。
[0068] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.2中，拟合后的灰度数据序列{LF口1, j为：
[0073] 其中Μ为每行灰度数据序列化i,j}中的数据个数。
[0074] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.5中，p = int(f/2)，int〇为取整函数。
[0075] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.8中，先采用高斯函数对数据序列化SF(t)}进行最小二乘法数据拟合，得到拟合后的 数据序列化SFfit(t)};再根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh)
[0081]
[0082] 有益效果
[0083] 本发明的整体技术效果体现在W下几个方面：
[0084] ( - )本发明提供的电视成像系统MTF测试测量方法可对多种不同视频接口的电视 成像系统的视频图像数据进行采集，可满足不同视频格式电视成像系统MTF准确测量，具有 测量准确性高的特点。
[0085] (二)本发明中，采用高均匀、高稳定的积分球光源作为照明光源，极大地提高了光 源的时间稳定性和空间均匀性，从而可在被测电视成像系统像面上形成稳定的目标狭缝 像，且沿垂直狭缝像方向像素灰度分布是非常光滑的高斯分布，运对提高亚像素线扩散函 数灰度计算精度、减小电视成像系统MTF的测量不确定度具有积极意义。
[0086] (Ξ)本发明采用电动可变衰减器W及程控电源对入射至待测电视成像系统祀面 的光斑强度进行衰减，实现了输出光强的自动调节；同时，采用电动可变衰减器W及程控电 源技术实现了输出光强的高精度连续调整，克服了无衰减措施时光强无法调整的缺点，从 而可满足不同增益、不同动态范围W及不同饱和灰度值的电视成像系统MTF的测量。另一方 面，采用该衰减措施将待测电视成像系统祀面处的光斑强度最大值控制为介于祀面接收强 度饱和值的30%~80%之间，有效避免了光源强度过强导致的待测电视成像系统器件饱和 或光源强度过弱导致图像的信噪比过低，从而合理利用了待测电视成像系统器件的动态范 围，进一步提高了电视成像系统MTF的测量重复性和准确性。
[0087] (四）本发明采用斜狭缝技术测量电视成像系统的MTF值，避免了采用垂直狭缝测 量离散器件MTF时，狭缝图像位于电视成像系统像面不同位置处时MTF测量结果不同的缺 点，使测量结果具有重复性和可比较性;另一方面，利用斜狭缝技术将灰度最大值横向变化 一个像素的一组数据进行合并可实现亚像素图像处理，从而使本发明提供的MTF测量方法 可实现频率高于奈奎斯特频率的MTF测量。
【附图说明】
[0088] 图1是本发明电视成像系统MTF测试方法中所用测试装置构成示意图。
[0089] 图2是本发明电视成像系统MTF测试方法中所用测试装置的积分球光源系统组成 不意图。
[0090] 图3是本发明电视成像系统MTF测试方法中所用测试装置的祀标组示意图。
[0091] 图4是本发明电视成像系统MTF测试方法中计算机的工作流程图。
【具体实施方式】
[0092] 下面结合具体实施例描述本发明：
[0093] 如图1所示，本实施例中的电视成像系统调制传递函数测试装置包括积分球光源 系统（1)、祀标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜(5)和计算机处理系 统(7)。
[0094]所述积分球光源系统（1)由积分球（1-1)、面鹤灯（1-2)、精密直流程控电源（1-3) 和电动可变光阔（1-4)组成。从而提供高时间稳定性、高空间均匀性且亮度可调的照明光， 计算机处理系统(7)控制精密直流程控电源（1-3)和/或电动可变光阔（1-4)调整积分球光 源系统（1)输出光强。在本优选实施例中，积分球光源系统1的积分球出口处照度均匀性可 达到5%，局部均匀性可达到2%，从而有效降低了由于光源的时间和空间不稳定或不均匀 而导致的MTF处理误差。积分球(1-1)出口端面安装有定位销，并且带有螺纹孔，用于安放和 固定祀标。
[009引所述祀标组(2)由单狭缝祀标、狭缝宽度约为2mm的多狭缝祀标和全黑祀标组成； 祀标组(2)的祀标边缘有缺口，祀标边缘缺口与定位销匹配使用。祀标通过螺钉固定在积分 球(1-1)出口处，根据需要通过人工更换的方式可选择祀标组(2)中的祀标固定在积分球光 源系统（1)的出口位置，当将该缺口对准定位销后，祀标上的狭缝位于积分球（1-1)出口中 央位置，且带有狭缝的祀标上的狭缝方向垂直于待测电视成像系统（8)的CCD祀面扫描方 向。本实施例中单狭缝祀标的狭缝宽度选取为50μπι、100μπι或200μπι，单狭缝祀标和多狭缝祀 标的狭缝开口直线性优于0.5%，避免由于狭缝自身不同位置处宽度不相等而引入测量误 差。
[0096] 所述滤光片组(3)安装在积分球光源系统（1)的出口位置，是一组中屯、波长不同且 范围在0.4皿~Ι.Ομπι内的窄带滤光片W及2片波长范围分别为0.4皿~0.化m和0.55皿~ 1. Ομπι的宽带滤光片，将从积分球光源系统（1)输出的光变为所需波长的光。
[0097] 所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统，祀标处于离 轴抛物面反射镜(5)的焦点位置，由祀标出射的光经过平面反射镜4和离轴抛物面反射镜5 的反射后变为平行光。在本实施例中，离轴抛物面反射镜5的口径为210mm，焦距为2000mm。
[0098] 此外，本实施例中还包括用于放置待测电视成像系统(8)的工作台（6)，工作台（6) 具有横滚、俯仰、升降和视场调整功能。
[0099] 所述计算机处理系统(7)与积分球光源系统（1)的精密直流程控电源（1-3)和电动 可变光阔（1-4)?及待测电视成像系统（8)连接。计算机处理系统（7)内置存储器、测试软 件，并配有鼠标、键盘和图像显示器，通过USB接口、网口及数字和模拟视频图像采集卡等设 备可与各类不同接口的待测电视成像系统连接，满足不同视频输出格式、不同视频接口的 电视成像系统的图像处理和MTF测量。
[0100] 待测的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控 制积分球光源系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行 图像采集，计算机处理系统(7)对采集的图像进行处理。
[0101] 下面将使用上述装置进行电视成像系统调制传递函数测试，在测试前，需要将工 作台（6)放置于测试装置的准直光路中，调整工作台（6)的位置使其中屯、轴与准直光路的中 屯、轴重合。将待测电视成像系统(8)放置于工作台（6)上，使用相应的数据线将待测电视成 像系统的视频接口与计算机处理系统(7)相应接口相连接。计算机处理系统(7)通过初始化 设置对待测电视成像系统的帖频W及图像采集时间等参数进行设置，从而将待测电视成像 系统的图像显示到计算机处理系统（7)的显示器上。初始化之后，开启电视成像系统MTF测 试装置并使其处于运行状态。调整工作台（6)的横滚、俯仰、升降和视场，使待测电视成像系 统的光轴与离轴抛物面反射镜(5)的光轴 重合。
[0102] 在本实施例中，待测电视成像系统的焦距为200mm，全分辨率输出像素数为NcniDsX Mcm〇s= 1600 X 1200,像原尺寸 Δ χΧ Δ y = 5.5皿X5.5皿。
[0103] 而测试的具体步骤为：
[0104] 步骤1:将祀标组(2)中的多狭缝祀标安装在积分球光源系统（1)出口位置，调节待 测电视成像系统(8)，对多狭缝祀标进行图像采集，且多狭缝祀标图像清晰呈现在计算机处 理系统(7)的显示器上。计算机处理系统(7)计算成像所得图像数据中最大灰度值与待测电 视成像系统(8)中探测器饱和灰度值的比值，若比值不满足设定要求，则调节积分球 光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝祀标上的狭缝个数为T。
[0105] 本实施例中，比值的设定要求为30%~80%，W保证待测电视成像系统的探测器 不会过饱和且又能充分利用其动态范围。当比值不满足设定要求时，计算机处理系统(7)首 先自动调节电动可变光阔（1-4)的大小，若改变电动可变光阔（1-4)达到其极限位置仍然不 能使比值满足设定要求，则计算机处理系统（7)调节精密直流程控电源1-3的输出电流，然 后再调节电动可变光阔（1-4)的大小，直至使比值满足设定要求。本实施例中Ifuii_t = 255。
[0106] 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝祀标采集的图像， 利用该图像并按照W下步骤获得系统放大率β:
[0107] 步骤2.1:采用二值化阔值Tbw对图像进行二值化处理，得到二值化图像Ibw;本实施 例中 Tbw=150;
[0108] 步骤2.2:在二值化图像Ibw中选取一个矩形区域Iroi，所述矩形区域Iroi高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧;本实施例中矩 形区域Iroi要尽量处于狭缝高度方向的中央位置；
[0109] 步骤2.3:计算矩形区域Iroi内每一行中灰度值为1的像素个数，形成数据序列 {Nw}：
[0110]
[0111] 并依次计算系统放大率β:
[0112]
[0113] 其中，v，w均为整数，v = a，a+l，日+2, = 6，b+l，b+l，…，Β; (Xv,yw)为矩形区域 Iroi中像素(v，w)的坐标，（xa，yb)为矩形区域Iroi中的像素起点（a，b)的坐标，（XA，yB)为矩形 区域Iroi中的像素终点(A，B)的坐标;lR〇i(Xv，yw)为矩形区域中（xv，yw)处的灰度值；Αχ为待 测电视成像系统(8)的像素尺寸，L为多狭缝祀标的单个狭缝宽度;f'为待测电视成像系统 (8)的焦距，时为准直光学系统焦距;xv = v Αχ。
[0114] 步骤3 :采用全黑祀标替换多狭缝祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对全黑祀标进行图像采集;计算机处理系统(7)提取待测电视成像 系统(8)对全黑祀标采集的图像，得到背景灰度Ibk;
[0115]
[0116] I(xi，yj)为本步骤采集的图像中，像素坐标(xi，yj)处的灰度值。
[0117] 步骤4:采用单狭缝祀标替换全黑祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对单狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图 像数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不 满足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测 电视成像系统(8)对单狭缝祀标再次进行图像采集，计算机处理系统(7)对采集的若干帖图 像进行平均，得到平均灰度图像lAV。
[0118] 本实施例中，比值设定要求为30%~80%，W保证待测电视成像系统的探测器不 会过饱和且又能充分利用其动态范围。积分球光源系统（1)输出光强的调节方式与步骤1相 同。
[0119] 本实施例中，计算机处理系统(7)对采集的不少于10帖图像进行平均，得到平均灰 度图像lAv。平均灰度图像中某一点的灰度值等于采集的不少于10帖图像对应点处灰度值的 平均值。
[0120] 步骤5:通过W下步骤测量传递函数：
[0121] 步骤5.1:在平均灰度图像Iav中选取一个矩形区域1/ AV，所述矩形区域1/ AV高度方 向上的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧;本实施例中矩形区 域I%V尽量处于狭缝高度方向的中央位置，且矩形区域I^AV高度为狭缝高度的1/4~1/2;
[0122] 步骤5.2:对I%v中的每一行，均采用高斯函数对该行的灰度数据序列化1,^进行最 小二乘法数据拟合，得到拟合后的灰度数据序列{LF 口 1,^，其中表示I%v中第j行的第i 点的灰度值，LFITi, j表示1/ AV中第j行的第i点的拟合灰度值;根据拟合后的灰度数据序列 (LF 口 ij}得到矩形区域I"Av;
[0123] 本实施例中，拟合后的灰度数据序列{LFITi,j}为：
[0128] 其中Μ为每行灰度数据序列化i,j}中的数据个数；
[0129] 步骤5.3:捜索矩形区域I"Av中每一行的灰度最大值，得到矩形区域I"Av中每一行 灰度最大值对应像素的坐标，其中矩形区域I"Av中第U行灰度最大值对应像素(max_u，u)的 坐标为(Xmax_u，yu )，得到数据序列{ Xmax_u，Yu }:
[01 30] Xmax u _ { Xmax_e , Xmax_e+1 , Xmax_e+2 ,......Xmax_E} j
[0131] Yu= {ye,ye+l,ye+2,......ye}；
[0132] 其中u = e，e+l，e+2,…E，ye为矩形区域I"AV中起始行的行坐标，yE为矩形区域I"AV 中终止行的行坐标；
[0133] 步骤5.4:对于矩形区域I"AV中的第U行，11 = e，e+1，e+2，…E，W该行灰度最大值对 应像素(max_u，u)为中屯、，计算中屯、两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝 对值的均值：
[0134]
[013引其中Nu同时满足条件
Xd为Γ AV中起始列的列坐标，XD为I" AV中 终止列的列坐标，Ψ为像素灰度差阔值，Φ为像素灰度最小阔值，Φ取值为背景灰度值Ibk; 本实施例中像素灰度差阔值Ψ取2;
[0136] 得到数据序列 Ay={S(ye)，S(ye+i)，S(ye+2)，......，S(yE)};
[0137] 步骤5.5:捜索数据序列Δ y的所有极小值：糾Λν, )，)，..........，S(打ν,);得 到所有极小值行坐标组成的新的数据序列，，化1%，·…，，.化Nf}，LNi~LNf是对应极 小值所在行在矩形区域I"Av中的行号;f是极小值个数;在数据序列Ω中任选两个相邻数据 LNp和LNp+i，得到倾斜角日为：
[013 引
[0139] 本实施例中，取p = int(f/2)，int〇为取整函数；
[0140] 步骤5.6:对矩形区域I"Av中，处于第LNp行和第LNp+l行之间的像素的灰度数据进行 合并，形成数据序列化S(q)}:
[0141] 当 α〉〇时，LS(q)为：
[014引其中
，int()为取整函数，modO为取余函数；q = l，2，···，RNX(D-d)，RN=LNp+广LNp;d为I"A冲起始列的像素编号，D为I"A冲终止列的像素 编号；
[0146] 步骤5.7:采用W下步骤对数据序列化S(q)}进行采样，得到数据序列化SF(t)}:
[0147] 步骤5.7.1:计算采样间隔INTI:
[014 引
[0149] 式中into为取整函数;S为采样率;S-般为大于1且小于5的整数，在本实施例中， 选取S为4进行计算；
[0150] 步骤5.7.2:根据采样间隔INT1计算MODq: MODq=mod (q/INT1)，mod ()为取余函数；
[0151] 步骤5.7.3:在数据序列化S(q)}捜索MODq为0的数据，组成数据序列化SF(t)}，数 据序列化SF( t)}中的数据个数为INT2;
[0152] 步骤5.8:根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):
[0162] 式中，Fh为空间频率，h为整数，h = 0，l，2,3,……Κ，Κ= INT2; Δ v为采样间隔，d为 I"A冲起始列的像素编号，D为I"Av中终止列的像素编号，峽J系统放大倍率，Αχ为待测电视 成像系统的像素的水平尺寸，1为单狭缝祀标的狭缝宽度。
[0163] 本实施例中，在步骤5.8中，先采用高斯函数对数据序列化SF(t)}进行最小二乘法 数据拟合，得到拟合后的数据序列化SFfit(t)};再根据W下公式计算Fo的h倍空间频率Fh的 调制传递函数MTF(Fh)
[0169] 计算机处理系统(7)完成上述计算后，将MTF测量值及图像在屏幕上显示并将其存 储于内部的存储器中，至此完成了单狭缝宽度下电视成像系统MTF测量。
[0170] 本优选实施例实现了对0.4μπι~1. Ομπι波长范围内电视成像系统MTF准确测量， MTF 测量不确定度达到0.1。
【主权项】
1. 一种电视成像系统调制传递函数测试装置，其特征在于:包括积分球光源系统（I)、 靶标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜(5)和计算机处理系统(7); 所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统； 所述积分球光源系统（1)的输出光强可调； 所述靶标组(2)由单狭缝靶标、多狭缝靶标和全黑靶标组成;根据需要选择靶标组（2) 中的靶标固定在积分球光源系统（1)的出口位置;靶标处于离轴抛物面反射镜(5)的焦点位 置； 所述滤光片组(3)安装在积分球光源系统（1)的出口位置;将从积分球光源系统（1)输 出的光变为所需波长的光； 所述计算机处理系统(7)与积分球光源系统（1)以及待测电视成像系统（8)连接，待测 的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控制积分球光源 系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行图像采集，计算 机处理系统(7)对采集的图像进行处理。2. 根据权利要求1所述一种电视成像系统调制传递函数测试装置，其特征在于:所述积 分球光源系统（1)由积分球（I -1 )、卤钨灯（1 -2 )、精密直流程控电源（1 -3)和电动可变光阑 (1-4)组成;计算机处理系统(7)控制精密直流程控电源（1-3)和/或电动可变光阑（1-4)调 整积分球光源系统（1)输出光强。3. 根据权利要求2所述一种电视成像系统调制传递函数测试装置，其特征在于:积分球 (1-1)出口端面安装有定位销;靶标组(2)的靶标边缘有缺口，靶标边缘缺口与定位销配合； 带有狭缝的靶标上的狭缝方向垂直于待测电视成像系统(8)的CCD靶面扫描方向。4. 根据权利要求1或2或3所述一种电视成像系统调制传递函数测试装置，其特征在于： 靶标组(2)中，单狭缝靶标的狭缝宽度选取为50μπι、100μπι或200μπι，单狭缝靶标和多狭缝靶 标的狭缝开口直线性优于〇. 5%。5. -种利用权利要求1所述装置进行电视成像系统调制传递函数测试的方法，其特征 在于:包括以下步骤： 步骤1:将靶标组(2)中的多狭缝靶标安装在积分球光源系统（1)出口位置，调节待测电 视成像系统(8)，对多狭缝靶标进行图像采集，计算机处理系统(7)计算成像所得图像数据 中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值I full_t的比值，若比值不满足设 定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝靶标 上的狭缝个数为T; 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝靶标采集的图像，利用 该图像并按照以下步骤获得系统放大率β: 步骤2.1:采用二值化阈值Tbw对图像进行二值化处理，得到二值化图像Ibw; 步骤2.2:在二值化图像Ibw中选取一个矩形区域I·，所述矩形区域〖隱高度方向上的两 条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧； 步骤2.3:计算矩形区域Irq1内每一行中灰度值为1的像素个数，形成数据序列{Nw} :并依次计算系统放大率β:其中，v，w均为整数，v = a，a+l，a+2, '"Aj = I^b+!，b+l，…，B; (XvJw)为矩形区域Iroi中 像素(v，w)的坐标，（Xa，yb)为矩形区域Iro冲的像素起点(a，b)的坐标，（ XA，yB)为矩形区域 Iroi中的像素终点(A，B)的坐标;Irqi (Xv，yw)为矩形区域中（Xv，yw)处的灰度值；Δ X为待测电 视成像系统(8)的像素尺寸，L为多狭缝靶标的单个狭缝宽度;f '为待测电视成像系统(8)的 焦距，fo为准直光学系统焦距； 步骤3:采用全黑靶标替换多狭缝靶标安装在积分球光源系统（1)出口位置，调节待测 电视成像系统(8)，对全黑靶标进行图像采集;计算机处理系统(7)提取待测电视成像系统 (8)对全黑靶标采集的图像，得到背景灰度Ibk; 步骤4:采用单狭缝靶标替换全黑靶标安装在积分球光源系统（1)出口位置，调节待测 电视成像系统(8)，对单狭缝靶标进行图像采集，计算机处理系统（7)计算成像所得图像数 据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值I full_t的比值，若比值不满足 设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测电视 成像系统(8)对单狭缝靶标再次进行图像采集，计算机处理系统(7)对采集的若干帧图像进 行平均，得到平均灰度图像Lw; 步骤5:通过以下步骤测量传递函数： 步骤5.1:在平均灰度图像Iav中选取一个矩形区域V AV，所述矩形区域I^v高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧； 步骤5.2:对V AV中的每一行，均采用高斯函数对该行的灰度数据序列{L^}进行最小二 乘法数据拟合，得到拟合后的灰度数据序列ILFIT1,丄其中表示I'AV中第j行的第i点的 灰度值，LFIT i,j表示F AV中第j行的第i点的拟合灰度值；由拟合后的灰度数据序列{LFITi,j} 得到矩形区域I〃av; 步骤5.3:搜索矩形区域I" AV中每一行的灰度最大值，得到矩形区域I〃AV中每一行灰度最 大值对应像素的坐标，其中矩形区域I〃m中第u行灰度最大值对应像素(maX_u，u)的坐标为 (Xmax_u，yu)，得到数据序列{Xmax_u，Yu}:其中u = e，e+l，e+2,…E，ye为矩形区域I〃AV中起始行的行坐标，yE为矩形区域I〃AV中终 止行的行坐标； 步骤5.4:对于矩形区域I" Av中的第u行，u = e，e+1，e+2，…E，以该行灰度最大值对应像 素(maX_u，u)为中心，计算中心两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝对值 的均值：其中Nu同时满足条件Xd为1〃 AV中起始列的列坐标，XD为1〃 AV中 终止列的列坐标，Ψ为像素灰度差阈值，Φ为像素灰度最小阈值，Φ取值为背景灰度值Ibk; 得到数据序列 Ay={5(ye)，5(ye+1)，5(ye+2)，......,δ(γ Ε)}； 步骤5.5:搜索数据序列&7的所有极小值：6〇^1)，0\^ :)，.一.，01>^：>;得到所 有极小值行坐标组成的新的数据序列〇=_!.V1 ，Λ.Ν2，…，Ymf I^LN1-LNf是对应极小值 所在行在矩形区域I" Av中的行号;f是极小值个数;在数据序列ω中任选两个相邻数据υαρ LNP+i，得到倾斜角α为：步骤5.6:对矩形区域I〃AV中，处于第LNp行和第LNP+1行之间的像素的灰度数据进行合 并，形成数据序列{LS(q)}:其中> int()为取整函数，mod()为取余函数;q=l， 2，-"，1^\(0-(1)，1^ = 1^+1-1^;(1为1〃;《中起始列的像素编号，0为1〃;《中终止列的像素编 号； 步骤5.7:采用以下步骤对数据序列{1^((1)}进行采样，得到数据序列{1^?(〇}: 步骤5.7.1:计算采样间隔INTl:式中int()为取整函数;S为采样率； 步骤5.7.2:根据采样间隔INTl计算MODq: M0Dq=mod (q/INTl )，mod 〇为取余函数； 步骤5.7.3:在数据序列{LS(q)}搜索MODq为0的数据，组成数据序列{LSF(t)}，数据序列 {LSF(t)}中的数据个数为INT2; 步骤5.8:根据以下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):其中：K= INT2; 1为单狭缝靶标的狭缝宽度。6. 根据权利要求5所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤1以及步骤4中比值的设定要求为30 %~80 %。7. 根据权利要求5或6所述一种进行电视成像系统调制传递函数测试的方法，其特征在 于：步骤4中，计算机处理系统（7)对采集的不少于10帧图像进行平均，得到平均灰度图像 Iav 〇8. 根据权利要求7所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤5.2中，拟合后的灰度数据序列{LFIT i, j}为：其中M为每行灰度数据序列{Li, j}中的数据个数。9. 根据权利要求8所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤5.5中，P = int(f/2)，int〇为取整函数。10. 根据权利要求9所述一种进行电视成像系统调制传递函数测试的方法，其特征在 于:步骤5.8中，先采用高斯函数对数据序列{LSF(t)}进行最小二乘法数据拟合，得到拟合 后的数据序列{LSF flt(t)};再根据以下公式计算Fo的h倍空间频率Fh的调制传递函数MTF (Fh)其中：
【专利摘要】本发明提出一种电视成像系统调制传递函数测试装置及方法，属于光学测量技术领域。其特点是将待测电视成像系统放置于由发光强度可调、空间均匀性高的积分球光源照明的准直光学系统中，准直光学系统的焦面放置狭缝靶标。狭缝靶标经过准直光学系统和待测电视成像光学系统后，成像到电视成像系统的靶面上，图像记录与处理系统对狭缝图像进行一系列计算和处理后，获得狭缝靶标的线扩散函数分布以及MTF值，由此完成电视成像系统MTF的测量。
【IPC分类】G01M11/02
【公开号】CN105486489
【申请号】CN201510882364
【发明人】王雷, 康登魁, 姜昌录, 杨斌, 马世帮, 杨红
【申请人】西安应用光学研究所
【公开日】2016年4月13日
【申请日】2015年12月4日
【技术领域】
[0001] 本发明设及光学计量与测量领域，具体为一种电视成像系统调制传递函数测试装 置及方法。
【背景技术】
[0002] 调制传递函数(MTF)是成像系统性能最基本的评价参数之一，反映了成像系统对 景物中各空间频率分量的响应特性，体现了成像系统对不同空间频率的目标对比度的传输 能力，因此MTF是评价成像系统成像性能的一项重要指标。
[0003] 可见光电视成像观瞄系统负责对目标进行捜索、跟踪和瞄准，其成像性能是决定 系统作用距离和跟踪精度的关键所在，因此必须对电视成像系统的成像性能进行准确测量 和评估。
[0004] 可见光电视成像观瞄系统通常是将光学系统和CCD成像器件固定在一起，由CCD将 光学系统所成的像转化为电信号，然后对电信号进行电子放大和处理，最终将视频图像信 号传输到显示器上，进行图像显示。由此可见，电视成像系统成像的性能不仅与光学系统成 像性能、CCD成像器件的性能有关，还与光学系统与CCD的匹配精度、电子电路的处理和转化 精度等有关。因此，对于电视成像系统的MTF评价方法既不同于普通光学透镜的像质评价， 也不同于CCD器件的性能评价方法。为了对电视成像系统的成像性能进行评价，就必须对电 视成像系统整机的MTF进行准确测量。
[000引电视成像系统MTF测试装置需要利用待测电视成像系统自身CCD器件采集并获得 祀标的图像，由计算机软件通过对CCD像面祀标的灰度分布进行分析计算得到待测电视成 像系统的MTF值。
[0006] 在实际中，由于电视成像系统MTF测试装置照明光源强度空间分布的不均匀、强度 不稳定、W及电视成像系统自身的特点，如接收器为离散型器件、具有不同灰度动态范围， 器件具有随机噪声等，会影响CCD器件像面图像的灰度分布W及线扩散函数的判读等，最终 影响MTF测量结果。
【发明内容】
[0007] 为了解决W上问题，本发明提供一种电视成像系统调制传递函数测试装置及方 法，可对不同视频格式、不同灰度等级、W及不同噪声水平的各类电视成像系统MTF测量仪 进行高准确度测量。
[0008] 本发明的技术方案为：
[0009] 所述一种电视成像系统调制传递函数测试装置，其特征在于:包括积分球光源系 统（1)、祀标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜巧)和计算机处理系统 (7);
[0010 ]所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统；
[0011]所述积分球光源系统（1)的输出光强可调；
[001引所述祀标组(2)由单狭缝祀标、多狭缝祀标和全黑祀标组成;根据需要选择祀标组 (2)中的祀标固定在积分球光源系统（1)的出口位置;祀标处于离轴抛物面反射镜(5)的焦 点位置；
[0013] 所述滤光片组（3)安装在积分球光源系统（1)的出口位置;将从积分球光源系统 (1)输出的光变为所需波长的光；
[0014] 所述计算机处理系统(7)与积分球光源系统（1)W及待测电视成像系统(8)连接， 待测的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控制积分球 光源系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行图像采集， 计算机处理系统(7)对采集的图像进行处理。
[0015] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 所述积分球光源系统（1)由积分球（1-1)、面鹤灯（1-2)、精密直流程控电源（1-3)和电动可 变光阔（1-4)组成;计算机处理系统（7)控制精密直流程控电源（1-3)和/或电动可变光阔 (1-4)调整积分球光源系统（1)输出光强。
[0016] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 积分球（1-1)出口端面安装有定位销;祀标组(2)的祀标边缘有缺口，祀标边缘缺口与定位 销配合;带有狭缝的祀标上的狭缝方向垂直于待测电视成像系统(8)的CCD祀面扫描方向。
[0017] 进一步的优选方案，所述一种电视成像系统调制传递函数测试装置，其特征在于： 祀标组(2)中，单狭缝祀标的狭缝宽度选取为50μπι、100Μ?或200μπι，单狭缝祀标和多狭缝祀 标的狭缝开口直线性优于0.5%。
[0018] 所述一种利用上述装置进行电视成像系统调制传递函数测试的方法，其特征在 于:包括W下步骤：
[0019] 步骤1:将祀标组(2)中的多狭缝祀标安装在积分球光源系统（1)出口位置，调节待 测电视成像系统(8)，对多狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图像 数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不满 足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝 祀标上的狭缝个数为T;
[0020] 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝祀标采集的图像， 利用该图像并按照W下步骤获得系统放大率β:
[0021 ] 步骤2.1:采用二值化阔值Tbw对图像进行二值化处理，得到二值化图像Ibw;
[0022] 步骤2.2:在二值化图像18冲选取一个矩形区域11?)1，所述矩形区域11?)1高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧；
[0023] 步骤2.3:计算矩形区域Iroi内每一行中灰度值为1的像素个数，形成数据序列 {Nw}：
[0024]
[0025] 并依次计算系统放大率β:
[0026]
[0027] 其中，v，w均为整数，v = a，a+l，a+2, = ，b+l，b+l，···，B; (Xv,yw)为矩形区域 Iroi中像素(v，w)的坐标，（xa，yb)为矩形区域Iroi中的像素起点（a，b)的坐标，（XA，yB)为矩形 区域Iroi中的像素终点(A，B)的坐标;Iroi (xv，yw)为矩形区域中（xv，yw)处的灰度值；Δ X为待 测电视成像系统(8)的像素尺寸，L为多狭缝祀标的单个狭缝宽度;f'为待测电视成像系统 (8)的焦距，时为准直光学系统焦距；
[0028] 步骤3:采用全黑祀标替换多狭缝祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对全黑祀标进行图像采集;计算机处理系统(7)提取待测电视成像 系统(8)对全黑祀标采集的图像，得到背景灰度Ibk;
[0029] 步骤4:采用单狭缝祀标替换全黑祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对单狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图 像数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不 满足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测 电视成像系统(8)对单狭缝祀标再次进行图像采集，计算机处理系统(7)对采集的若干帖图 像进行平均，得到平均灰度图像lAV;
[0030] 步骤5:通过W下步骤测量传递函数：
[0031] 步骤5.1:在平均灰度图像lAv中选取一个矩形区域1/ AV，所述矩形区域1/ AV高度方 向上的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧；
[0032] 步骤5.2:对I%v中的每一行，均采用高斯函数对该行的灰度数据序列化1,^进行最 小二乘法数据拟合，得到拟合后的灰度数据序列{LF 口 1,^，其中表示I%v中第j行的第i 点的灰度值，LFITi,读示I^Av中第j行的第i点的拟合灰度值；由拟合后的灰度数据序列 (LF 口 ij}得到矩形区域I"Av;
[0033] 步骤5.3:捜索矩形区域I"Av中每一行的灰度最大值，得到矩形区域I"Av中每一行 灰度最大值对应像素的坐标，其中矩形区域I"Av中第U行灰度最大值对应像素(max_u，u)的 坐标为(Xmax_u，yu)，得到数据序列{Xmax_u，Yu}:
[0034] Xmax u _ { Xmax_e , Xmax_e+1 , Xmax_e+2 ,......Xmax_E} j
[003引 Yu= {ye，ye+i，ye+2,......ye}；
[0036] 其中u = e，e+l，e+2, ···EiYe为矩形区域I"AV中起始行的行坐标，yE为矩形区域I"AV 中终止行的行坐标；
[0037] 步骤5.4:对于矩形区域I"AV中的第U行，11 = e，e+1，e+2，…E，W该行灰度最大值对 应像素(max_u，u)为中屯、，计算中屯、两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝 对值的均值：
[003引
[0039] 其中亂同时满足条件沁3[/^争諷',,心.>',,）-^枯。,一,"，乂,）]<半、^片纖,,+.^^^ ，凡）<Φ、. Xn?__u+W,, - X。和 - %; Xd为 I" AV 中起始列的列坐柄，XD 为 I" AV 中 终止列的列坐标，Ψ为像素灰度差阔值，Φ为像素灰度最小阔值，Φ取值为背景灰度值Ibk;
[0040] 得到数据序列 Ay={S(ye)，S(ye+i)，S(ye+2)，......，S(yE)};
[0041 ] 步骤5.5 : 捜索数据序列Δ y的所有极小值：4执& )，)，.….，，、）；得 到所有极小值行坐标组成的新的数据序列〇={.Vln,，Λν,，...，丸Kf }，L化~LNf是对应极 小值所在行在矩形区域I"Av中的行号;f是极小值个数;在数据序列Ω中任选两个相邻数据 LNp和LNp+i，得到倾斜角日为：
[0042] "二 arctan[A、'乃少.|、 - 、；)];
[0043] 步骤5.6:对矩形区域I"Av中，处于第LNp行和第LNp+l行之间的像素的灰度数据进行 合并，形成数据序列化S(q)}:
[0044] 当 α〉〇时，LS(q)为：
[004引其中
，into为取整函数，modO为取余函数;q = l，2，…，RNX(D-d)，RN=LNp+广LNp;d为I"A冲起始列的像素编号，D为I"A冲终止列的像素 编号；
[0049] 步骤5.7:采用W下步骤对数据序列化S(q)}进行采样，得到数据序列化SF(t)}:
[0050] 步骤5.7.1:计算采样间隔INTI:
[0051]
[0052] 式中into为取整函数;S为采样率；
[00閲步骤5.7.2:根据采样间隔1饥'1计算版)0。:1孤。=1110(1(9/1饥'1)，1110(10为取余函数；
[0054]步骤5.7.3:在数据序列化S(q)}捜索MODq为0的数据，组成数据序列化SF(t)}，数 据序列化SF(t)}中的数据个数为INT2;
[005引步骤5.8:根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):
[0065] Κ = INT2; 1为单狭缝祀标的狭缝宽度。
[0066] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤及步骤4中比值的设定要求为30 %~80 %。
[0067] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤4中，计算机处理系统(7)对采集的不少于10帖图像进行平均，得到平均灰度图像lAv。
[0068] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.2中，拟合后的灰度数据序列{LF口1, j为：
[0073] 其中Μ为每行灰度数据序列化i,j}中的数据个数。
[0074] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.5中，p = int(f/2)，int〇为取整函数。
[0075] 进一步的优选方案，所述一种电视成像系统调制传递函数测试方法，其特征在于： 步骤5.8中，先采用高斯函数对数据序列化SF(t)}进行最小二乘法数据拟合，得到拟合后的 数据序列化SFfit(t)};再根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh)
[0081]
[0082] 有益效果
[0083] 本发明的整体技术效果体现在W下几个方面：
[0084] ( - )本发明提供的电视成像系统MTF测试测量方法可对多种不同视频接口的电视 成像系统的视频图像数据进行采集，可满足不同视频格式电视成像系统MTF准确测量，具有 测量准确性高的特点。
[0085] (二)本发明中，采用高均匀、高稳定的积分球光源作为照明光源，极大地提高了光 源的时间稳定性和空间均匀性，从而可在被测电视成像系统像面上形成稳定的目标狭缝 像，且沿垂直狭缝像方向像素灰度分布是非常光滑的高斯分布，运对提高亚像素线扩散函 数灰度计算精度、减小电视成像系统MTF的测量不确定度具有积极意义。
[0086] (Ξ)本发明采用电动可变衰减器W及程控电源对入射至待测电视成像系统祀面 的光斑强度进行衰减，实现了输出光强的自动调节；同时，采用电动可变衰减器W及程控电 源技术实现了输出光强的高精度连续调整，克服了无衰减措施时光强无法调整的缺点，从 而可满足不同增益、不同动态范围W及不同饱和灰度值的电视成像系统MTF的测量。另一方 面，采用该衰减措施将待测电视成像系统祀面处的光斑强度最大值控制为介于祀面接收强 度饱和值的30%~80%之间，有效避免了光源强度过强导致的待测电视成像系统器件饱和 或光源强度过弱导致图像的信噪比过低，从而合理利用了待测电视成像系统器件的动态范 围，进一步提高了电视成像系统MTF的测量重复性和准确性。
[0087] (四）本发明采用斜狭缝技术测量电视成像系统的MTF值，避免了采用垂直狭缝测 量离散器件MTF时，狭缝图像位于电视成像系统像面不同位置处时MTF测量结果不同的缺 点，使测量结果具有重复性和可比较性;另一方面，利用斜狭缝技术将灰度最大值横向变化 一个像素的一组数据进行合并可实现亚像素图像处理，从而使本发明提供的MTF测量方法 可实现频率高于奈奎斯特频率的MTF测量。
【附图说明】
[0088] 图1是本发明电视成像系统MTF测试方法中所用测试装置构成示意图。
[0089] 图2是本发明电视成像系统MTF测试方法中所用测试装置的积分球光源系统组成 不意图。
[0090] 图3是本发明电视成像系统MTF测试方法中所用测试装置的祀标组示意图。
[0091] 图4是本发明电视成像系统MTF测试方法中计算机的工作流程图。
【具体实施方式】
[0092] 下面结合具体实施例描述本发明：
[0093] 如图1所示，本实施例中的电视成像系统调制传递函数测试装置包括积分球光源 系统（1)、祀标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜(5)和计算机处理系 统(7)。
[0094]所述积分球光源系统（1)由积分球（1-1)、面鹤灯（1-2)、精密直流程控电源（1-3) 和电动可变光阔（1-4)组成。从而提供高时间稳定性、高空间均匀性且亮度可调的照明光， 计算机处理系统(7)控制精密直流程控电源（1-3)和/或电动可变光阔（1-4)调整积分球光 源系统（1)输出光强。在本优选实施例中，积分球光源系统1的积分球出口处照度均匀性可 达到5%，局部均匀性可达到2%，从而有效降低了由于光源的时间和空间不稳定或不均匀 而导致的MTF处理误差。积分球(1-1)出口端面安装有定位销，并且带有螺纹孔，用于安放和 固定祀标。
[009引所述祀标组(2)由单狭缝祀标、狭缝宽度约为2mm的多狭缝祀标和全黑祀标组成； 祀标组(2)的祀标边缘有缺口，祀标边缘缺口与定位销匹配使用。祀标通过螺钉固定在积分 球(1-1)出口处，根据需要通过人工更换的方式可选择祀标组(2)中的祀标固定在积分球光 源系统（1)的出口位置，当将该缺口对准定位销后，祀标上的狭缝位于积分球（1-1)出口中 央位置，且带有狭缝的祀标上的狭缝方向垂直于待测电视成像系统（8)的CCD祀面扫描方 向。本实施例中单狭缝祀标的狭缝宽度选取为50μπι、100μπι或200μπι，单狭缝祀标和多狭缝祀 标的狭缝开口直线性优于0.5%，避免由于狭缝自身不同位置处宽度不相等而引入测量误 差。
[0096] 所述滤光片组(3)安装在积分球光源系统（1)的出口位置，是一组中屯、波长不同且 范围在0.4皿~Ι.Ομπι内的窄带滤光片W及2片波长范围分别为0.4皿~0.化m和0.55皿~ 1. Ομπι的宽带滤光片，将从积分球光源系统（1)输出的光变为所需波长的光。
[0097] 所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统，祀标处于离 轴抛物面反射镜(5)的焦点位置，由祀标出射的光经过平面反射镜4和离轴抛物面反射镜5 的反射后变为平行光。在本实施例中，离轴抛物面反射镜5的口径为210mm，焦距为2000mm。
[0098] 此外，本实施例中还包括用于放置待测电视成像系统(8)的工作台（6)，工作台（6) 具有横滚、俯仰、升降和视场调整功能。
[0099] 所述计算机处理系统(7)与积分球光源系统（1)的精密直流程控电源（1-3)和电动 可变光阔（1-4)?及待测电视成像系统（8)连接。计算机处理系统（7)内置存储器、测试软 件，并配有鼠标、键盘和图像显示器，通过USB接口、网口及数字和模拟视频图像采集卡等设 备可与各类不同接口的待测电视成像系统连接，满足不同视频输出格式、不同视频接口的 电视成像系统的图像处理和MTF测量。
[0100] 待测的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控 制积分球光源系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行 图像采集，计算机处理系统(7)对采集的图像进行处理。
[0101] 下面将使用上述装置进行电视成像系统调制传递函数测试，在测试前，需要将工 作台（6)放置于测试装置的准直光路中，调整工作台（6)的位置使其中屯、轴与准直光路的中 屯、轴重合。将待测电视成像系统(8)放置于工作台（6)上，使用相应的数据线将待测电视成 像系统的视频接口与计算机处理系统(7)相应接口相连接。计算机处理系统(7)通过初始化 设置对待测电视成像系统的帖频W及图像采集时间等参数进行设置，从而将待测电视成像 系统的图像显示到计算机处理系统（7)的显示器上。初始化之后，开启电视成像系统MTF测 试装置并使其处于运行状态。调整工作台（6)的横滚、俯仰、升降和视场，使待测电视成像系 统的光轴与离轴抛物面反射镜(5)的光轴 重合。
[0102] 在本实施例中，待测电视成像系统的焦距为200mm，全分辨率输出像素数为NcniDsX Mcm〇s= 1600 X 1200,像原尺寸 Δ χΧ Δ y = 5.5皿X5.5皿。
[0103] 而测试的具体步骤为：
[0104] 步骤1:将祀标组(2)中的多狭缝祀标安装在积分球光源系统（1)出口位置，调节待 测电视成像系统(8)，对多狭缝祀标进行图像采集，且多狭缝祀标图像清晰呈现在计算机处 理系统(7)的显示器上。计算机处理系统(7)计算成像所得图像数据中最大灰度值与待测电 视成像系统(8)中探测器饱和灰度值的比值，若比值不满足设定要求，则调节积分球 光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝祀标上的狭缝个数为T。
[0105] 本实施例中，比值的设定要求为30%~80%，W保证待测电视成像系统的探测器 不会过饱和且又能充分利用其动态范围。当比值不满足设定要求时，计算机处理系统(7)首 先自动调节电动可变光阔（1-4)的大小，若改变电动可变光阔（1-4)达到其极限位置仍然不 能使比值满足设定要求，则计算机处理系统（7)调节精密直流程控电源1-3的输出电流，然 后再调节电动可变光阔（1-4)的大小，直至使比值满足设定要求。本实施例中Ifuii_t = 255。
[0106] 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝祀标采集的图像， 利用该图像并按照W下步骤获得系统放大率β:
[0107] 步骤2.1:采用二值化阔值Tbw对图像进行二值化处理，得到二值化图像Ibw;本实施 例中 Tbw=150;
[0108] 步骤2.2:在二值化图像Ibw中选取一个矩形区域Iroi，所述矩形区域Iroi高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧;本实施例中矩 形区域Iroi要尽量处于狭缝高度方向的中央位置；
[0109] 步骤2.3:计算矩形区域Iroi内每一行中灰度值为1的像素个数，形成数据序列 {Nw}：
[0110]
[0111] 并依次计算系统放大率β:
[0112]
[0113] 其中，v，w均为整数，v = a，a+l，日+2, = 6，b+l，b+l，…，Β; (Xv,yw)为矩形区域 Iroi中像素(v，w)的坐标，（xa，yb)为矩形区域Iroi中的像素起点（a，b)的坐标，（XA，yB)为矩形 区域Iroi中的像素终点(A，B)的坐标;lR〇i(Xv，yw)为矩形区域中（xv，yw)处的灰度值；Αχ为待 测电视成像系统(8)的像素尺寸，L为多狭缝祀标的单个狭缝宽度;f'为待测电视成像系统 (8)的焦距，时为准直光学系统焦距;xv = v Αχ。
[0114] 步骤3 :采用全黑祀标替换多狭缝祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对全黑祀标进行图像采集;计算机处理系统(7)提取待测电视成像 系统(8)对全黑祀标采集的图像，得到背景灰度Ibk;
[0115]
[0116] I(xi，yj)为本步骤采集的图像中，像素坐标(xi，yj)处的灰度值。
[0117] 步骤4:采用单狭缝祀标替换全黑祀标安装在积分球光源系统（1)出口位置，调节 待测电视成像系统(8)，对单狭缝祀标进行图像采集，计算机处理系统(7)计算成像所得图 像数据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值IfuiLt的比值，若比值不 满足设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测 电视成像系统(8)对单狭缝祀标再次进行图像采集，计算机处理系统(7)对采集的若干帖图 像进行平均，得到平均灰度图像lAV。
[0118] 本实施例中，比值设定要求为30%~80%，W保证待测电视成像系统的探测器不 会过饱和且又能充分利用其动态范围。积分球光源系统（1)输出光强的调节方式与步骤1相 同。
[0119] 本实施例中，计算机处理系统(7)对采集的不少于10帖图像进行平均，得到平均灰 度图像lAv。平均灰度图像中某一点的灰度值等于采集的不少于10帖图像对应点处灰度值的 平均值。
[0120] 步骤5:通过W下步骤测量传递函数：
[0121] 步骤5.1:在平均灰度图像Iav中选取一个矩形区域1/ AV，所述矩形区域1/ AV高度方 向上的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧;本实施例中矩形区 域I%V尽量处于狭缝高度方向的中央位置，且矩形区域I^AV高度为狭缝高度的1/4~1/2;
[0122] 步骤5.2:对I%v中的每一行，均采用高斯函数对该行的灰度数据序列化1,^进行最 小二乘法数据拟合，得到拟合后的灰度数据序列{LF 口 1,^，其中表示I%v中第j行的第i 点的灰度值，LFITi, j表示1/ AV中第j行的第i点的拟合灰度值;根据拟合后的灰度数据序列 (LF 口 ij}得到矩形区域I"Av;
[0123] 本实施例中，拟合后的灰度数据序列{LFITi,j}为：
[0128] 其中Μ为每行灰度数据序列化i,j}中的数据个数；
[0129] 步骤5.3:捜索矩形区域I"Av中每一行的灰度最大值，得到矩形区域I"Av中每一行 灰度最大值对应像素的坐标，其中矩形区域I"Av中第U行灰度最大值对应像素(max_u，u)的 坐标为(Xmax_u，yu )，得到数据序列{ Xmax_u，Yu }:
[01 30] Xmax u _ { Xmax_e , Xmax_e+1 , Xmax_e+2 ,......Xmax_E} j
[0131] Yu= {ye,ye+l,ye+2,......ye}；
[0132] 其中u = e，e+l，e+2,…E，ye为矩形区域I"AV中起始行的行坐标，yE为矩形区域I"AV 中终止行的行坐标；
[0133] 步骤5.4:对于矩形区域I"AV中的第U行，11 = e，e+1，e+2，…E，W该行灰度最大值对 应像素(max_u，u)为中屯、，计算中屯、两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝 对值的均值：
[0134]
[013引其中Nu同时满足条件
Xd为Γ AV中起始列的列坐标，XD为I" AV中 终止列的列坐标，Ψ为像素灰度差阔值，Φ为像素灰度最小阔值，Φ取值为背景灰度值Ibk; 本实施例中像素灰度差阔值Ψ取2;
[0136] 得到数据序列 Ay={S(ye)，S(ye+i)，S(ye+2)，......，S(yE)};
[0137] 步骤5.5:捜索数据序列Δ y的所有极小值：糾Λν, )，)，..........，S(打ν,);得 到所有极小值行坐标组成的新的数据序列，，化1%，·…，，.化Nf}，LNi~LNf是对应极 小值所在行在矩形区域I"Av中的行号;f是极小值个数;在数据序列Ω中任选两个相邻数据 LNp和LNp+i，得到倾斜角日为：
[013 引
[0139] 本实施例中，取p = int(f/2)，int〇为取整函数；
[0140] 步骤5.6:对矩形区域I"Av中，处于第LNp行和第LNp+l行之间的像素的灰度数据进行 合并，形成数据序列化S(q)}:
[0141] 当 α〉〇时，LS(q)为：
[014引其中
，int()为取整函数，modO为取余函数；q = l，2，···，RNX(D-d)，RN=LNp+广LNp;d为I"A冲起始列的像素编号，D为I"A冲终止列的像素 编号；
[0146] 步骤5.7:采用W下步骤对数据序列化S(q)}进行采样，得到数据序列化SF(t)}:
[0147] 步骤5.7.1:计算采样间隔INTI:
[014 引
[0149] 式中into为取整函数;S为采样率;S-般为大于1且小于5的整数，在本实施例中， 选取S为4进行计算；
[0150] 步骤5.7.2:根据采样间隔INT1计算MODq: MODq=mod (q/INT1)，mod ()为取余函数；
[0151] 步骤5.7.3:在数据序列化S(q)}捜索MODq为0的数据，组成数据序列化SF(t)}，数 据序列化SF( t)}中的数据个数为INT2;
[0152] 步骤5.8:根据W下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):
[0162] 式中，Fh为空间频率，h为整数，h = 0，l，2,3,……Κ，Κ= INT2; Δ v为采样间隔，d为 I"A冲起始列的像素编号，D为I"Av中终止列的像素编号，峽J系统放大倍率，Αχ为待测电视 成像系统的像素的水平尺寸，1为单狭缝祀标的狭缝宽度。
[0163] 本实施例中，在步骤5.8中，先采用高斯函数对数据序列化SF(t)}进行最小二乘法 数据拟合，得到拟合后的数据序列化SFfit(t)};再根据W下公式计算Fo的h倍空间频率Fh的 调制传递函数MTF(Fh)
[0169] 计算机处理系统(7)完成上述计算后，将MTF测量值及图像在屏幕上显示并将其存 储于内部的存储器中，至此完成了单狭缝宽度下电视成像系统MTF测量。
[0170] 本优选实施例实现了对0.4μπι~1. Ομπι波长范围内电视成像系统MTF准确测量， MTF 测量不确定度达到0.1。
【主权项】
1. 一种电视成像系统调制传递函数测试装置，其特征在于:包括积分球光源系统（I)、 靶标组(2)、滤光片组(3)、平面反射镜(4)、离轴抛物面反射镜(5)和计算机处理系统(7); 所述平面反射镜(4)和所述离轴抛物面反射镜(5)组成准直光学系统； 所述积分球光源系统（1)的输出光强可调； 所述靶标组(2)由单狭缝靶标、多狭缝靶标和全黑靶标组成;根据需要选择靶标组（2) 中的靶标固定在积分球光源系统（1)的出口位置;靶标处于离轴抛物面反射镜(5)的焦点位 置； 所述滤光片组(3)安装在积分球光源系统（1)的出口位置;将从积分球光源系统（1)输 出的光变为所需波长的光； 所述计算机处理系统(7)与积分球光源系统（1)以及待测电视成像系统（8)连接，待测 的电视成像系统(8)接收准直光学系统输出的平行光;计算机处理系统(7)控制积分球光源 系统（1)调整输出光强，计算机处理系统(7)控制待测电视成像系统(8)进行图像采集，计算 机处理系统(7)对采集的图像进行处理。2. 根据权利要求1所述一种电视成像系统调制传递函数测试装置，其特征在于:所述积 分球光源系统（1)由积分球（I -1 )、卤钨灯（1 -2 )、精密直流程控电源（1 -3)和电动可变光阑 (1-4)组成;计算机处理系统(7)控制精密直流程控电源（1-3)和/或电动可变光阑（1-4)调 整积分球光源系统（1)输出光强。3. 根据权利要求2所述一种电视成像系统调制传递函数测试装置，其特征在于:积分球 (1-1)出口端面安装有定位销;靶标组(2)的靶标边缘有缺口，靶标边缘缺口与定位销配合； 带有狭缝的靶标上的狭缝方向垂直于待测电视成像系统(8)的CCD靶面扫描方向。4. 根据权利要求1或2或3所述一种电视成像系统调制传递函数测试装置，其特征在于： 靶标组(2)中，单狭缝靶标的狭缝宽度选取为50μπι、100μπι或200μπι，单狭缝靶标和多狭缝靶 标的狭缝开口直线性优于〇. 5%。5. -种利用权利要求1所述装置进行电视成像系统调制传递函数测试的方法，其特征 在于:包括以下步骤： 步骤1:将靶标组(2)中的多狭缝靶标安装在积分球光源系统（1)出口位置，调节待测电 视成像系统(8)，对多狭缝靶标进行图像采集，计算机处理系统(7)计算成像所得图像数据 中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值I full_t的比值，若比值不满足设 定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;所述多狭缝靶标 上的狭缝个数为T; 步骤2:计算机处理系统(7)提取待测电视成像系统(8)对多狭缝靶标采集的图像，利用 该图像并按照以下步骤获得系统放大率β: 步骤2.1:采用二值化阈值Tbw对图像进行二值化处理，得到二值化图像Ibw; 步骤2.2:在二值化图像Ibw中选取一个矩形区域I·，所述矩形区域〖隱高度方向上的两 条边均处在狭缝上，宽度方向上的两条边分别处于T个狭缝整体的两侧； 步骤2.3:计算矩形区域Irq1内每一行中灰度值为1的像素个数，形成数据序列{Nw} :并依次计算系统放大率β:其中，v，w均为整数，v = a，a+l，a+2, '"Aj = I^b+!，b+l，…，B; (XvJw)为矩形区域Iroi中 像素(v，w)的坐标，（Xa，yb)为矩形区域Iro冲的像素起点(a，b)的坐标，（ XA，yB)为矩形区域 Iroi中的像素终点(A，B)的坐标;Irqi (Xv，yw)为矩形区域中（Xv，yw)处的灰度值；Δ X为待测电 视成像系统(8)的像素尺寸，L为多狭缝靶标的单个狭缝宽度;f '为待测电视成像系统(8)的 焦距，fo为准直光学系统焦距； 步骤3:采用全黑靶标替换多狭缝靶标安装在积分球光源系统（1)出口位置，调节待测 电视成像系统(8)，对全黑靶标进行图像采集;计算机处理系统(7)提取待测电视成像系统 (8)对全黑靶标采集的图像，得到背景灰度Ibk; 步骤4:采用单狭缝靶标替换全黑靶标安装在积分球光源系统（1)出口位置，调节待测 电视成像系统(8)，对单狭缝靶标进行图像采集，计算机处理系统（7)计算成像所得图像数 据中最大灰度值与待测电视成像系统(8)中探测器饱和灰度值I full_t的比值，若比值不满足 设定要求，则调节积分球光源系统（1)的输出光强，直至比值满足设定要求;而后待测电视 成像系统(8)对单狭缝靶标再次进行图像采集，计算机处理系统(7)对采集的若干帧图像进 行平均，得到平均灰度图像Lw; 步骤5:通过以下步骤测量传递函数： 步骤5.1:在平均灰度图像Iav中选取一个矩形区域V AV，所述矩形区域I^v高度方向上 的两条边均处在狭缝上，宽度方向上的两条边分别处于狭缝两侧； 步骤5.2:对V AV中的每一行，均采用高斯函数对该行的灰度数据序列{L^}进行最小二 乘法数据拟合，得到拟合后的灰度数据序列ILFIT1,丄其中表示I'AV中第j行的第i点的 灰度值，LFIT i,j表示F AV中第j行的第i点的拟合灰度值；由拟合后的灰度数据序列{LFITi,j} 得到矩形区域I〃av; 步骤5.3:搜索矩形区域I" AV中每一行的灰度最大值，得到矩形区域I〃AV中每一行灰度最 大值对应像素的坐标，其中矩形区域I〃m中第u行灰度最大值对应像素(maX_u，u)的坐标为 (Xmax_u，yu)，得到数据序列{Xmax_u，Yu}:其中u = e，e+l，e+2,…E，ye为矩形区域I〃AV中起始行的行坐标，yE为矩形区域I〃AV中终 止行的行坐标； 步骤5.4:对于矩形区域I" Av中的第u行，u = e，e+1，e+2，…E，以该行灰度最大值对应像 素(maX_u，u)为中心，计算中心两侧Nu个对称像素点的灰度值之差的绝对值，并得到绝对值 的均值：其中Nu同时满足条件Xd为1〃 AV中起始列的列坐标，XD为1〃 AV中 终止列的列坐标，Ψ为像素灰度差阈值，Φ为像素灰度最小阈值，Φ取值为背景灰度值Ibk; 得到数据序列 Ay={5(ye)，5(ye+1)，5(ye+2)，......,δ(γ Ε)}； 步骤5.5:搜索数据序列&7的所有极小值：6〇^1)，0\^ :)，.一.，01>^：>;得到所 有极小值行坐标组成的新的数据序列〇=_!.V1 ，Λ.Ν2，…，Ymf I^LN1-LNf是对应极小值 所在行在矩形区域I" Av中的行号;f是极小值个数;在数据序列ω中任选两个相邻数据υαρ LNP+i，得到倾斜角α为：步骤5.6:对矩形区域I〃AV中，处于第LNp行和第LNP+1行之间的像素的灰度数据进行合 并，形成数据序列{LS(q)}:其中> int()为取整函数，mod()为取余函数;q=l， 2，-"，1^\(0-(1)，1^ = 1^+1-1^;(1为1〃;《中起始列的像素编号，0为1〃;《中终止列的像素编 号； 步骤5.7:采用以下步骤对数据序列{1^((1)}进行采样，得到数据序列{1^?(〇}: 步骤5.7.1:计算采样间隔INTl:式中int()为取整函数;S为采样率； 步骤5.7.2:根据采样间隔INTl计算MODq: M0Dq=mod (q/INTl )，mod 〇为取余函数； 步骤5.7.3:在数据序列{LS(q)}搜索MODq为0的数据，组成数据序列{LSF(t)}，数据序列 {LSF(t)}中的数据个数为INT2; 步骤5.8:根据以下公式计算Fo的h倍空间频率Fh的调制传递函数MTF(Fh):其中：K= INT2; 1为单狭缝靶标的狭缝宽度。6. 根据权利要求5所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤1以及步骤4中比值的设定要求为30 %~80 %。7. 根据权利要求5或6所述一种进行电视成像系统调制传递函数测试的方法，其特征在 于：步骤4中，计算机处理系统（7)对采集的不少于10帧图像进行平均，得到平均灰度图像 Iav 〇8. 根据权利要求7所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤5.2中，拟合后的灰度数据序列{LFIT i, j}为：其中M为每行灰度数据序列{Li, j}中的数据个数。9. 根据权利要求8所述一种进行电视成像系统调制传递函数测试的方法，其特征在于： 步骤5.5中，P = int(f/2)，int〇为取整函数。10. 根据权利要求9所述一种进行电视成像系统调制传递函数测试的方法，其特征在 于:步骤5.8中，先采用高斯函数对数据序列{LSF(t)}进行最小二乘法数据拟合，得到拟合 后的数据序列{LSF flt(t)};再根据以下公式计算Fo的h倍空间频率Fh的调制传递函数MTF (Fh)其中：
【专利摘要】本发明提出一种电视成像系统调制传递函数测试装置及方法，属于光学测量技术领域。其特点是将待测电视成像系统放置于由发光强度可调、空间均匀性高的积分球光源照明的准直光学系统中，准直光学系统的焦面放置狭缝靶标。狭缝靶标经过准直光学系统和待测电视成像光学系统后，成像到电视成像系统的靶面上，图像记录与处理系统对狭缝图像进行一系列计算和处理后，获得狭缝靶标的线扩散函数分布以及MTF值，由此完成电视成像系统MTF的测量。
【IPC分类】G01M11/02
【公开号】CN105486489
【申请号】CN201510882364
【发明人】王雷, 康登魁, 姜昌录, 杨斌, 马世帮, 杨红
【申请人】西安应用光学研究所
【公开日】2016年4月13日
【申请日】2015年12月4日
最新回复(0)