CN108174186B - 一种空间面阵红外相机的点扩散函数测试方法 - Google Patents
一种空间面阵红外相机的点扩散函数测试方法 Download PDFInfo
- Publication number
- CN108174186B CN108174186B CN201711320928.8A CN201711320928A CN108174186B CN 108174186 B CN108174186 B CN 108174186B CN 201711320928 A CN201711320928 A CN 201711320928A CN 108174186 B CN108174186 B CN 108174186B
- Authority
- CN
- China
- Prior art keywords
- asterism
- battle array
- infrared camera
- face battle
- 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.)
- Active
Links
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
-
- 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
- G01M11/02—Testing optical properties
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radiation Pyrometers (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Image Processing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
一种空间面阵红外相机的点扩散函数测试方法,涉及航天光学遥感技术领域;包括如下步骤:步骤(一)、建立靶标;步骤(二)、获取第M帧靶标图像中的第i星点的星点灰度值图像;步骤(三)、定义5行×5列的矩阵E;步骤(四)、选取所需星点灰度值图像;步骤(五)、将所需星点灰度值图像的灰度值转换为电子数;步骤(六)、拟合得到面阵红外相机的二维热扩散分布函数;步骤(七)、图像复原,提高遥感图像质量;本发明有效提高空间面阵红外PSF测量精度,填补了空间面阵红外相机的PSF测试方法空白,具有很好的使用价值。
Description
技术领域
本发明涉及一种航天光学遥感技术领域,特别是一种空间面阵红外相机的点扩散函数测试方法。
背景技术
光学***在理想状态下,物空间一点发出的光能量在像空间也集中在一点上,但是实际的光学***成像时,由于衍射和像差以及其它因素的影响,物空间一点发出的光在像空间是分布在一定的区域内,其分布曲线称为点扩散函数PSF。高慧婷等人2014年在国家发明专利《一种星载推扫式相机点扩散函数的实验室测量装置及方法》一文中给出适用于推扫式可见光相机的PSF测试方法;刘薇等人2014年在国家发明专利《一种星载面阵CCD相机点扩散函数的测量装置及方法》给出凝视型面阵可见光相机的PSF测试方法。上述两篇专利都是将靶标设与数据处理一体化考虑,同时考虑了相位精度问题,但是这两种方法一个是针对推扫式可见光相机,一个是针对凝视型面阵可见光相机的PSF测试方法,红外相机PSF测试面临相位更难找准以及扩散和时间相关等问题,上述方法并没有考虑。
发明内容
本发明的目的在于克服现有技术的上述不足,提供一种空间面阵红外相机的点扩散函数测试方法,通过靶标设计和数据处理方法解决高精度相位对准和时间扩散问题,保证空间面阵红外相机PSF的测试准确性,最大程度上满足了空间面阵红外相机PSF测试的需求。
本发明的上述目的是通过如下技术方案予以实现的:
一种空间面阵红外相机的点扩散函数测试方法,包括如下步骤:
步骤(一)、建立靶标,在靶标上设置100个星点,100个星点呈10行×10列均匀对称分布;并采用面阵红外相机对靶标进行图像采集,采集帧数为N,N>100,且N为正整数;
步骤(二)、对第M帧靶标图像中的第i星点及该星点相邻像素进行采样;得到星点灰度值子图mi;其中,1≤M≤N,且N为正整数;1≤i≤100,且i为正整数;
步骤(三)、定义一个5行×5列的矩阵E,
步骤(四)、将步骤(三)中的矩阵E与步骤(二)中的星点灰度值子图mi进行计算,求每一帧中每个星的星点灰度值J(i);
得到N帧中每一帧每个星的星点灰度值J(i);取每一帧中的最大星点灰度值J(i),即为星点灰度值图像;
步骤(五)、将步骤(四)得到的星点灰度值图像的灰度值转换为电子数xe;
步骤(六)、针对步骤(四)中得到的最大星点灰度值J(i)所对应的第i星点,建立平面坐标系oxy;x与y方向垂直,交点为o;将第i星点设置在坐标系oxy的原点;将第i星点的相邻像素对应设置在平面坐标系oxy中,相邻像素点坐标值间隔为1;对步骤(五)中得到的电子数进行拟合;
将最大星点灰度值J(i)所对应的第i星点的灰度值图像中的所有像素进行拟合,计算得到所需第一参数σ1和第二参数σ2;即得到面阵红外相机的二维热扩散分布函数u(t,σ1,σ2);
步骤(七)、通过步骤(六)得到的面阵红外相机的二维热扩散分布函数N(t,σh1,σh2)对红外面阵相机的成像数据进行图像复原,实现提高遥感图像质量。
在上述的一种空间面阵红外相机的点扩散函数测试方法,所述步骤(一)中,对靶标进行图像采集时,每一帧图像为100个星点的整体图像。
在上述的一种空间面阵红外相机的点扩散函数测试方法,所述步骤(二)中,对第i星点相邻像素进行采样的方法为;以第i星点为中心,采集5行×5列的方形相邻像素。
在上述的一种空间面阵红外相机的点扩散函数测试方法,所述步骤(四)中,每个星的星点灰度值J(i)的计算方法为:
式中,J(i)为第i星点灰度值;
j为矩阵E的行数;
k为矩阵E的列数。
在上述的一种空间面阵红外相机的点扩散函数测试方法,所述步骤(五)中,星点灰度值图像的灰度值转换为电子数xe的转换方法为:
式中,x为图像灰度值序列中采样点的灰度值;
bs为面阵红外相机自身量化位数;
L为面阵红外相机的饱和电子数。
在上述的一种空间面阵红外相机的点扩散函数测试方法,所述步骤(六)中,电子数的拟合方法为:
式中,t为面阵红外相机的积分时间;
σ1为二维热扩散分布函数第一参数;
σ2为二维热扩散分布函数第一参数;
u(t,σ1,σ2)为红外面阵相机不同积分时间的点扩散函数;
xq为第i星相邻第q像素对应在平面坐标系oxy横坐标;1≤q≤25,q为正整数;
yq为第i星相邻第q像素对应在平面坐标系oxy纵坐标;1≤q≤25,q为正整数。
本发明与现有技术相比具有如下优点:
(1)本发明通过高精度相位配准方法,克服了红外相机噪声严重导致PSF测试不准的缺陷;
(2)本发明采用时间相关的拟合函数,克服红外相机热扩散导致的PSF测试不准的问题;
(3)本发明通过星点位置和形状的设计,使得靶标中至少有一个星点完全落在一个像元的中心位置,从而找到用于计算PSF函数的目标点序列,增加了测试装置和方法的准确性;
(4)本发明通过统计方法对多次测量结果进行了平差处理,降低了噪声影响,减小了测量误差。
附图说明
图1为本发明面阵红外相机的点扩散函数测试流程图;
图2为本发明靶标建立示意图;
图3为本发明第i星点相邻像素采样示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示为面阵红外相机的点扩散函数测试流程图,由图可知,一种空间面阵红外相机的点扩散函数测试方法,包括如下步骤:
步骤(一)、建立靶标,如图2所示为靶标建立示意图,由图可知,在靶标上设置100个星点,100个星点呈10行×10列均匀对称分布;并采用面阵红外相机对靶标进行图像采集,采集帧数为N,N>100,且N为正整数;对靶标进行图像采集时,每一帧图像为100个星点的整体图像。
步骤(二)、如图3所示为第i星点相邻像素采样示意图,由图可知,对第M帧靶标图像中的第i星点及该星点相邻像素进行采样;得到星点灰度值子图mi;其中,1≤M≤N,且N为正整数;1≤i≤100,且i为正整数;对第i星点相邻像素进行采样的方法为;以第i星点为中心,采集5行×5列的方形相邻像素。
步骤(三)、定义一个5行×5列的矩阵E,
步骤(四)、将步骤(三)中的矩阵E与步骤(二)中的星点灰度值子图mi进行计算,求每一帧中每个星的星点灰度值J(i):
式中,J(i)为第i星点灰度值;
j为矩阵E的行数;
k为矩阵E的列数;
得到N帧中每一帧每个星的星点灰度值J(i);取每一帧中的最大星点灰度值J(i),即为星点灰度值图像;
步骤(五)、将步骤(四)得到的星点灰度值图像的灰度值转换为电子数xe,转换公式为:
式中,x为图像灰度值序列中采样点的灰度值;
bs为面阵红外相机自身量化位数;
L为面阵红外相机的饱和电子数;
步骤(六)、针对步骤(四)中得到的最大星点灰度值J(i)所对应的第i星点,建立平面坐标系oxy;x与y方向垂直,交点为o;将第i星点设置在坐标系oxy的原点;将第i星点的相邻像素对应设置在平面坐标系oxy中,相邻像素点坐标值间隔为1;对步骤(五)中得到的电子数进行拟合,拟合函数为:
式中,t为面阵红外相机的积分时间;
σ1为二维热扩散分布函数第一参数;
σ2为二维热扩散分布函数第一参数;
u(t,σ1,σ2)为红外面阵相机不同积分时间的点扩散函数;
xq为第i星相邻第q像素对应在平面坐标系oxy横坐标;1≤q≤25,q为正整数;
yq为第i星相邻第q像素对应在平面坐标系oxy纵坐标;1≤q≤25,q为正整数;
将最大星点灰度值J(i)所对应的第i星点的灰度值图像中的所有像素代入公式(三)拟合,计算得到所需第一参数σ1和第二参数σ2;即得到面阵红外相机的二维热扩散分布函数u(t,σ1,σ2);
步骤(七)、通过步骤(六)得到的面阵红外相机的二维热扩散分布函数N(t,σh1,σh2)对红外面阵相机的成像数据进行图像复原,实现提高遥感图像质量。
本发明设计了新的相位匹配和PSF拟合方法,实现了面阵红外相机PSF函数的测试,可以直观地描述相机性能,可应用于空间面阵红外相机的点扩散函数测试。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (4)
1.一种空间面阵红外相机的点扩散函数测试方法,其特征在于:包括如下步骤:
步骤(一)、建立靶标,在靶标上设置100个星点,100个星点呈10行×10列均匀对称分布;并采用面阵红外相机对靶标进行图像采集,采集帧数为N,N>100,且N为正整数;
步骤(二)、对第M帧靶标图像中的第i星点及该星点相邻像素进行采样;得到星点灰度值子图mi;其中,1≤M≤N,且N为正整数;1≤i≤100,且i为正整数;
步骤(三)、定义一个5行×5列的矩阵E,
步骤(四)、将步骤(三)中的矩阵E与步骤(二)中的星点灰度值子图mi进行计算,求每一帧中每个星的星点灰度值J(i);每个星的星点灰度值J(i)的计算方法为:
式中,J(i)为第i星点灰度值;
j为矩阵E的行数;
k为矩阵E的列数;
得到N帧中每一帧每个星的星点灰度值J(i);取每一帧中的最大星点灰度值J(i),即为星点灰度值图像;
步骤(五)、将步骤(四)得到的星点灰度值图像的灰度值转换为电子数xe;
步骤(六)、针对步骤(四)中得到的最大星点灰度值J(i)所对应的第i星点,建立平面坐标系oxy;x与y方向垂直,交点为o;将第i星点设置在坐标系oxy的原点;将第i星点的相邻像素对应设置在平面坐标系oxy中,相邻像素点坐标值间隔为1;对步骤(五)中得到的电子数进行拟合;
电子数的拟合方法为:
式中,t为面阵红外相机的积分时间;
σ1为二维热扩散分布函数第一参数;
σ2为二维热扩散分布函数第一参数;
u(t,σ1,σ2)为红外面阵相机不同积分时间的点扩散函数;
xq为第i星相邻第q像素对应在平面坐标系oxy横坐标;1≤q≤25,q为正整数;
yq为第i星相邻第q像素对应在平面坐标系oxy纵坐标;1≤q≤25,q为正整数;
将最大星点灰度值J(i)所对应的第i星点的灰度值图像中的所有像素进行拟合,计算得到所需第一参数σ1和第二参数σ2;即得到面阵红外相机的二维热扩散分布函数u(t,σ1,σ2);
步骤(七)、通过步骤(六)得到的面阵红外相机的二维热扩散分布函数N(t,σh1,σh2)对红外面阵相机的成像数据进行图像复原,实现提高遥感图像质量。
2.根据权利要求1所述的一种空间面阵红外相机的点扩散函数测试方法,其特征在于:所述步骤(一)中,对靶标进行图像采集时,每一帧图像为100个星点的整体图像。
3.根据权利要求1所述的一种空间面阵红外相机的点扩散函数测试方法,其特征在于:所述步骤(二)中,对第i星点相邻像素进行采样的方法为;以第i星点为中心,采集5行×5列的方形相邻像素。
4.根据权利要求1所述的一种空间面阵红外相机的点扩散函数测试方法,其特征在于:所述步骤(五)中,星点灰度值图像的灰度值转换为电子数xe的转换方法为:
式中,x为图像灰度值序列中采样点的灰度值;
bs为面阵红外相机自身量化位数;
L为面阵红外相机的饱和电子数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711320928.8A CN108174186B (zh) | 2017-12-12 | 2017-12-12 | 一种空间面阵红外相机的点扩散函数测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711320928.8A CN108174186B (zh) | 2017-12-12 | 2017-12-12 | 一种空间面阵红外相机的点扩散函数测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108174186A CN108174186A (zh) | 2018-06-15 |
CN108174186B true CN108174186B (zh) | 2019-08-09 |
Family
ID=62525731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711320928.8A Active CN108174186B (zh) | 2017-12-12 | 2017-12-12 | 一种空间面阵红外相机的点扩散函数测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108174186B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115586638A (zh) * | 2022-10-10 | 2023-01-10 | 长春理工大学 | 可见光宽波段含单层衍射元件***点扩散函数构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103957404A (zh) * | 2014-04-22 | 2014-07-30 | 北京空间机电研究所 | 一种星载推扫式相机点扩散函数的实验室测量装置及方法 |
CN104111080A (zh) * | 2014-06-24 | 2014-10-22 | 中国资源卫星应用中心 | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 |
CN104236859A (zh) * | 2014-09-19 | 2014-12-24 | 北京空间机电研究所 | 一种星载面阵ccd相机点扩散函数测量装置及方法 |
CN106651928A (zh) * | 2016-12-07 | 2017-05-10 | 北京空间机电研究所 | 一种空间面阵相机亚像元点扩散函数测量装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100691268B1 (ko) * | 2005-08-02 | 2007-03-12 | 삼성전기주식회사 | Psf를 이용하여 이미지를 처리하는 광학 시스템 및이미지 처리방법 |
-
2017
- 2017-12-12 CN CN201711320928.8A patent/CN108174186B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103957404A (zh) * | 2014-04-22 | 2014-07-30 | 北京空间机电研究所 | 一种星载推扫式相机点扩散函数的实验室测量装置及方法 |
CN104111080A (zh) * | 2014-06-24 | 2014-10-22 | 中国资源卫星应用中心 | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 |
CN104236859A (zh) * | 2014-09-19 | 2014-12-24 | 北京空间机电研究所 | 一种星载面阵ccd相机点扩散函数测量装置及方法 |
CN106651928A (zh) * | 2016-12-07 | 2017-05-10 | 北京空间机电研究所 | 一种空间面阵相机亚像元点扩散函数测量装置及方法 |
Non-Patent Citations (2)
Title |
---|
基于点源阵列的星载面阵CMOS相机静态PSF测量;高慧婷等;《航天返回与遥感》;20170228;第38卷(第1期);53-60 * |
点扩散函数高斯拟合估计与遥感图像恢复;杨利红等;《中国光学》;20120430;第5卷(第2期);181-188 * |
Also Published As
Publication number | Publication date |
---|---|
CN108174186A (zh) | 2018-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106556461B (zh) | 一种基于自适应光学的光谱成像装置 | |
CN106228532B (zh) | 一种lcd屏幕亚像素级缺陷检测方法 | |
CN104236859B (zh) | 一种星载面阵ccd相机点扩散函数测量装置及方法 | |
Abreu et al. | Prediction of diffuse horizontal irradiance using a new climate zone model | |
CN106971408B (zh) | 一种基于时空转换思想的摄像机标定方法 | |
CN113177512B (zh) | 一种星星间交叉辐射定标的匹配阈值分析方法 | |
CN106352806A (zh) | 一种立体视觉三维数字图像相关测量的高精度标定方法 | |
CN103984979A (zh) | 无透镜衍射成像的藻类自动检测计数装置及方法 | |
CN106289532A (zh) | 一种红外热图像的温度提取方法及装置 | |
CN105547342A (zh) | 基于液晶面板的工业镜头测试装置及方法 | |
CN107562791A (zh) | 一种基于大数据统计的遥感卫星相对辐射定标处理方法 | |
CN106709944A (zh) | 卫星遥感图像配准方法 | |
CN102752504B (zh) | 一种宽视场线阵ccd相机的相对辐射校正方法 | |
CN103389310A (zh) | 一种基于辐射标定的亚像素光学元件损伤在线检测方法 | |
CN111076698A (zh) | 用于标定光电产品安装偏差的观察评估装置及解算方法 | |
CN107796595A (zh) | 低温光学***能量集中度测试***及方法 | |
CN107560541A (zh) | 图像中心偏差的测量方法及装置 | |
CN103957404B (zh) | 一种星载推扫式相机点扩散函数的实验室测量装置及方法 | |
CN108174186B (zh) | 一种空间面阵红外相机的点扩散函数测试方法 | |
CN113850868B (zh) | 一种波浪爬高图像识别方法 | |
CN106651928B (zh) | 一种空间面阵相机亚像元点扩散函数测量装置及方法 | |
CN101726316A (zh) | 内方位元素及畸变测试仪 | |
CN109064515B (zh) | 一种对月观测数据提取方法及提取装置 | |
CN103778637A (zh) | 一种基于直方图平衡的side-slither辐射定标处理方法 | |
CN103913126B (zh) | 用于角秒级三维变形的光学准直测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |