CN104833480A - 一种测试刃边倾角对调制传递函数影响的装置及方法 - Google Patents
一种测试刃边倾角对调制传递函数影响的装置及方法 Download PDFInfo
- Publication number
- CN104833480A CN104833480A CN201510205022.6A CN201510205022A CN104833480A CN 104833480 A CN104833480 A CN 104833480A CN 201510205022 A CN201510205022 A CN 201510205022A CN 104833480 A CN104833480 A CN 104833480A
- Authority
- CN
- China
- Prior art keywords
- target
- camera
- sword
- image
- transfer function
- 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
Abstract
本发明一种测试刃边倾角对调制传递函数影响的装置及方法,装置包括光源系统、靶标、平行光管、相机、图像采集系统;靶标放置在平行光管的焦面上,平行光管与相机的光轴共轴放置,靶标被光源系统产生的均匀光照亮,经过平行光管和相机,在相机的焦面上成像,形成靶标的图像,图像采集系统对图像进行采集和处理。由于在设计中考虑了使不同倾斜角度的刃边真值具有较好的一致性,从而使各角度MTF测量值差异能够反映刃边角度的变化。本发明可用于相机MTF测量和靶标角度优选,填补MTF测试方法相关空白。
Description
技术领域
本发明属于航天光学遥感技术领域,涉及一种测试刃边倾角对调制传递函数影响的装置及方法,特别涉及靶标设计、实验方法和实验数据处理方法。
背景技术
光学系统在理想状态下,物空间一点发出的光能量在像空间也集中在一点上,但是实际的光学系统成像时,由于衍射和像差以及其它工艺的影响,物空间一点发出的光在像空间是分布在一定的区域内,其分布曲线称为PSF。为了简化处理,假设PSF各向同性,因此可以不考虑相位信息,只考虑PSF的幅值,称作调制传递函数(MTF),MTF是评价光学系统成像质量的重要客观指标,为了定量化评价相机成像质量,发射之前都进行实验室的MTF测试,发射之后还要进行在轨MTF测试。Becker等人2005年在ISPRS文章中总结提出了四种MTF测试方法,第一种是点脉冲法,第二种基于线脉冲法,第三种是刃边法,第四种是辐射状靶标法,黄巧林等人2006年在《航天返回与遥感》上发表的《航天光学遥感器MTF测试技术研究》一文中还给出了第五种方法,高对比度矩形靶标测试方法。卫星发射后,基于可以使用的靶标——实际地面特征物的可及性,刃边法测试MTF应用最为广泛,但是目前尚无直观描述刃边角度影响变化规律的实验方法。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供一种测试刃边倾角对调制传递函数影响的装置及方法,通过靶标设计和数据处理方法可同时研究多个刃边角度的变化,刻画相机静态调制传递函数的变化规律。
本发明的技术方案是:一种测试刃边倾角对调制传递函数影响的装置,包括光源系统、靶标、平行光管、相机、图像采集系统;靶标放置在平行光管的焦面上,平行光管与相机的光轴共轴放置,靶标被光源系统产生的均匀光照亮,经过平行光管和相机,成像在相机的焦面上,形成靶标的图像,图像采集系统对图像进行采集和处理;
所述的靶标整体为矩形,长宽比为2:1;内部分解为2行共16个同样大小的矩形单元,每行8个矩形单元,每个矩形单元长宽比为2:1,均为黑白刃边单元;各矩形单元内部,刃边与短边方向夹角为α;从第一行第一个矩形单元开始至第一行最后一个矩形单元,α以1°为步长依次取值为74°-81°;从第二行第一个矩形单元开始至第二行最后一个矩形单元,α以1°为步长依次取值为82°-89°;所述的靶标上所有矩形单元内部刃边一侧透光另一侧不透光,两部分对比度不低于100:1。
一种利用上述装置进行测试的方法,步骤如下:
1)将靶标放置在相机视场的中心位置,并固定靶标的位置;
2)利用图像采集系统采集得到某帧靶标图像;
3)在靶标图像上,按照上述靶标的内部构成分割出全部16个刃边单元;
4)从第一个刃边单元开始,根据图像灰度提取该刃边单元内刃边点的亚像元位置;
5)依据步骤4)获得的刃边点的亚像元位置提取采样数据,拟合边缘扩展函数曲线;
6)对步骤5)求得的边缘扩展函数曲线一次求导,得出线扩展函数;
7)对步骤6)得到的线扩展函数曲线做傅立叶变换得到该刃边单元的调制传递函数;
8)重复步骤4)-步骤7),得到该帧靶标图像的16个刃边单元的调制传递函数;
9)根据步骤8)得出的调制传递函数,得出调制传递函数随角度α的变化曲线。
本发明与现有技术相比的优点在于:
(1)本发明克服星载相机地面测试阶段刃边靶标角度单一导致测量样本不足的缺陷,通过具有连续角度的刃边靶标进行MTF测试,可以同时获得多个角度的MTF;
(2)本发明通过靶标设计和数据处理方法保证测量精度,靶标中包含16个角度靶标,基本覆盖了常用的刃边角度范围;
(3)本发明可以利用现有MTF测试系统设备,易于工程实现。
附图说明
图1为本发明MTF测试靶标设计图;
图2为本发明MTF实验装置示意图;
图3为本发明数据处理流程图。
具体实施方式
本发明涉及的实验装置如图2所示,包括光源系统1、靶标2、平行光管3、相机4、图像采集系统5;靶标2放置在平行光管3的焦面上,平行光管3与相机4的光轴共轴放置,靶标2被光源系统1产生的均匀光照亮,经过平行光管3和相机4,成像在相机4的焦面上,形成靶标2的图像,图像采集系统5对成像结果进行采集和处理;
如图1所示,靶标2整体为正方形。内部分解为16个同样大小的黑白刃边单元,均为矩形,每个矩形长宽比为2:1。各矩形内部,刃边与短边方向夹角为α,α以1°为步长依次在74°-89°之间取值。
靶标2上所有矩形单元内部刃边一侧透光另一侧不透光,两部分对比度不低于100:1。
本发明一种实验室中利用上述测量装置测量相机系统MTF的方法,分为测量步骤和对靶标图像进行处理步骤,具体如下:
一、测量步骤:
(一)将靶标2安装在平行光管3的焦平面上,调整好光源系统1的位置;
(二)将相机4摆放在平行光管3出口前,通过相机调整装置调整好相机的位置,保证相机入光口与平行光管3的出光口对准,平行光管3的有效口径充满相机入瞳;
(三)连接好相机4与图像采集系统5之间的电缆,并通电检查,确保工作正常;
(四)打开光源系统1;
(五)固定靶标2的位置;
(六)调整光源系统1亮度,使靶标图像亮区输出灰度值不饱和且不低于饱和值的0.6倍;
(七)调整相机焦面至最佳焦面位置;
利用图像采集系统5采集得到靶标2的图像采集并存储图像并处理,如图3步骤如下:
(一)在靶标图像上,按照上述靶标2的内部构成分割出全部16个刃边单元;
(二)从第一个刃边单元开始,根据图像灰度提取该刃边单元内刃边点的亚像元位置;
(三)依据步骤(二)获得的刃边点的亚像元位置提取采样数据,拟合边缘扩展函数曲线;
(四)对步骤(三)求得的边缘扩展函数曲线一次求导,得出线扩展函数;
(五)对步骤(四)得到的线扩展函数曲线做傅立叶变换得到该刃边单元的调制传递函数;
(六)重复步骤(一)-步骤(四),得到该帧靶标图像的16个刃边单元的调制传递函数;
(七)根据步骤(六)得出的调制传递函数,得出调制传递函数随角度α的变化曲线。
本发明说明书中未作详细描述的内容属本领于技术人员的公知技术。
Claims (2)
1.一种测试刃边倾角对调制传递函数影响的装置,其特征在于:包括光源系统(1)、靶标(2)、平行光管(3)、相机(4)、图像采集系统(5);靶标(2)放置在平行光管(3)的焦面上,平行光管(3)与相机(4)的光轴共轴放置,靶标(2)被光源系统(1)产生的均匀光照亮,经过平行光管(3)和相机(4),成像在相机(4)的焦面上,形成靶标(2)的图像,图像采集系统(5)对图像进行采集和处理;
所述的靶标(2)整体为矩形,长宽比为2:1;内部分解为2行共16个同样大小的矩形单元,每行8个矩形单元,每个矩形单元长宽比为2:1,均为黑白刃边单元;各矩形单元内部,刃边与短边方向夹角为α;从第一行第一个矩形单元开始至第一行最后一个矩形单元,α以1°为步长依次取值为74°-81°;从第二行第一个矩形单元开始至第二行最后一个矩形单元,α以1°为步长依次取值为82°-89°;所述的靶标(2)上所有矩形单元内部刃边一侧透光另一侧不透光,两部分对比度不低于100:1。
2.一种利用权利要求1所述装置进行测试的方法,其特征在于步骤如下:
1)将靶标(2)放置在相机(4)视场的中心位置,并固定靶标(2)的位置;
2)利用图像采集系统(5)采集得到某帧靶标图像;
3)在靶标图像上,按照上述靶标(2)的内部构成分割出全部16个刃边单元;
4)从第一个刃边单元开始,根据图像灰度提取该刃边单元内刃边点的亚像元位置;
5)依据步骤4)获得的刃边点的亚像元位置提取采样数据,拟合边缘扩展函数曲线;
6)对步骤5)求得的边缘扩展函数曲线一次求导,得出线扩展函数;
7)对步骤6)得到的线扩展函数曲线做傅立叶变换得到该刃边单元的调制传递函数;
8)重复步骤4)-步骤7),得到该帧靶标图像的16个刃边单元的调制传递函数;
9)根据步骤8)得出的调制传递函数,得出调制传递函数随角度α的变化曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510205022.6A CN104833480A (zh) | 2015-04-27 | 2015-04-27 | 一种测试刃边倾角对调制传递函数影响的装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510205022.6A CN104833480A (zh) | 2015-04-27 | 2015-04-27 | 一种测试刃边倾角对调制传递函数影响的装置及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104833480A true CN104833480A (zh) | 2015-08-12 |
Family
ID=53811504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510205022.6A Pending CN104833480A (zh) | 2015-04-27 | 2015-04-27 | 一种测试刃边倾角对调制传递函数影响的装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104833480A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427267A (zh) * | 2015-10-28 | 2016-03-23 | 中国科学院光电研究院 | 获取刃边法测量光学遥感载荷在轨mtf精度的方法 |
CN106289319A (zh) * | 2016-07-29 | 2017-01-04 | 北京空间机电研究所 | 一种基于在轨序列图像的卫星系统传函测算方法 |
CN109990985A (zh) * | 2019-04-25 | 2019-07-09 | 绍兴文理学院 | 一种品字型线列红外探测器调制传递函数测试方法 |
WO2021026515A1 (en) * | 2019-08-07 | 2021-02-11 | Agilent Technologies, Inc. | Optical imaging performance test system and method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050254041A1 (en) * | 2004-05-12 | 2005-11-17 | Edmund Optics, Inc. | Tilted edge for optical-transfer-function measurement |
JP2007309764A (ja) * | 2006-05-18 | 2007-11-29 | Acutelogic Corp | Mtf測定装置、mtf測定方法およびmtf測定プログラム |
CN101458156A (zh) * | 2007-12-12 | 2009-06-17 | 鸿富锦精密工业(深圳)有限公司 | 调制传递函数值的测量装置及其测量方法 |
CN104236859A (zh) * | 2014-09-19 | 2014-12-24 | 北京空间机电研究所 | 一种星载面阵ccd相机点扩散函数测量装置及方法 |
-
2015
- 2015-04-27 CN CN201510205022.6A patent/CN104833480A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050254041A1 (en) * | 2004-05-12 | 2005-11-17 | Edmund Optics, Inc. | Tilted edge for optical-transfer-function measurement |
JP2007309764A (ja) * | 2006-05-18 | 2007-11-29 | Acutelogic Corp | Mtf測定装置、mtf測定方法およびmtf測定プログラム |
CN101458156A (zh) * | 2007-12-12 | 2009-06-17 | 鸿富锦精密工业(深圳)有限公司 | 调制传递函数值的测量装置及其测量方法 |
CN104236859A (zh) * | 2014-09-19 | 2014-12-24 | 北京空间机电研究所 | 一种星载面阵ccd相机点扩散函数测量装置及方法 |
Non-Patent Citations (3)
Title |
---|
戴奇燕 等: "刀刃法在轨MTF测量性能分析", 《航天返回与遥感》 * |
相春昌: "CCD相机调制传递函数的测量", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
赵占平 等: "基于刃边法的航天光学遥感器在轨MTF测试研究", 《航天返回与遥感》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427267A (zh) * | 2015-10-28 | 2016-03-23 | 中国科学院光电研究院 | 获取刃边法测量光学遥感载荷在轨mtf精度的方法 |
CN105427267B (zh) * | 2015-10-28 | 2018-03-23 | 中国科学院光电研究院 | 获取刃边法测量光学遥感载荷在轨mtf精度的方法 |
CN106289319A (zh) * | 2016-07-29 | 2017-01-04 | 北京空间机电研究所 | 一种基于在轨序列图像的卫星系统传函测算方法 |
CN106289319B (zh) * | 2016-07-29 | 2019-06-18 | 北京空间机电研究所 | 一种基于在轨序列图像的卫星系统传函测算方法 |
CN109990985A (zh) * | 2019-04-25 | 2019-07-09 | 绍兴文理学院 | 一种品字型线列红外探测器调制传递函数测试方法 |
CN109990985B (zh) * | 2019-04-25 | 2021-01-05 | 绍兴文理学院 | 一种品字型线列红外探测器调制传递函数测试方法 |
WO2021026515A1 (en) * | 2019-08-07 | 2021-02-11 | Agilent Technologies, Inc. | Optical imaging performance test system and method |
CN114174791A (zh) * | 2019-08-07 | 2022-03-11 | 安捷伦科技有限公司 | 光学成像性能测试系统和方法 |
AU2020325162B2 (en) * | 2019-08-07 | 2023-08-31 | Agilent Technologies, Inc. | Optical imaging performance test system and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108764257B (zh) | 一种多视角的指针式仪表识别方法 | |
CN105841810B (zh) | 全光成像系统的物空间校准 | |
CN106596073A (zh) | 一种检测光学系统像质的方法和系统及一种测试标板 | |
CN104833480A (zh) | 一种测试刃边倾角对调制传递函数影响的装置及方法 | |
CN106885622B (zh) | 一种大视场多点三维振动测量方法 | |
JP6594294B2 (ja) | 顕微鏡画像の画像品質評価 | |
CN110428762B (zh) | 基于像素点亮度的oled面板老化测试发光特性检测方法 | |
CN104236859A (zh) | 一种星载面阵ccd相机点扩散函数测量装置及方法 | |
CN109186759B (zh) | 一种光栅光谱仪像质测量方法和装置 | |
CN105842062A (zh) | 一种裂纹扩展实时监测装置及方法 | |
TW202004149A (zh) | 檢測浸於溶液中眼用鏡片屈光度及厚度之系統及方法 | |
CN104111165A (zh) | 一种近眼显示光学镜头像质的评价装置 | |
CN110763600A (zh) | 悬浮颗粒实时在线检测装置 | |
CN109712089A (zh) | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 | |
CN105510364A (zh) | 基于x射线的工业部件缺陷的无损检测系统及其检测方法 | |
CN108254072A (zh) | 一种新型高光谱视频成像仪 | |
CN103234475A (zh) | 一种基于激光三角测量法的亚像素表面形态检测方法 | |
CN110108715A (zh) | 一种透明平行平板的缺陷检测方法 | |
CN110958450A (zh) | 一种成像系统空间测试装置、对比度及频率测试方法 | |
CN103957404A (zh) | 一种星载推扫式相机点扩散函数的实验室测量装置及方法 | |
CN106773158A (zh) | 一种移动式自寻位液晶屏像素质量分析装置和方法 | |
CN108205210A (zh) | 基于傅里叶梅林及特征匹配的lcd缺陷检测系统及方法 | |
CN106651928B (zh) | 一种空间面阵相机亚像元点扩散函数测量装置及方法 | |
Boden et al. | Accuracy of measurement of dynamic surface deformations by the image pattern correlation technique | |
CN109632269A (zh) | 基于图像灰度信息检测光学衍射元件性能的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150812 |
|
RJ01 | Rejection of invention patent application after publication |