CN104111080B - 一种凝视卫星面阵ccd相机的mtf在轨检测方法 - Google Patents
一种凝视卫星面阵ccd相机的mtf在轨检测方法 Download PDFInfo
- Publication number
- CN104111080B CN104111080B CN201410286149.0A CN201410286149A CN104111080B CN 104111080 B CN104111080 B CN 104111080B CN 201410286149 A CN201410286149 A CN 201410286149A CN 104111080 B CN104111080 B CN 104111080B
- Authority
- CN
- China
- Prior art keywords
- ccd camera
- array ccd
- area array
- point source
- point
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
一种凝视卫星面阵CCD相机的MTF在轨检测方法。首先建立点源物分布与像分布之间端对端的点扩散函数模型,并对点源进行设计和布设;其次求取点源区每个像元的点扩散函数值,并利用最小范数求解得到面阵CCD相机的点扩散函数参数,进而得到光学传递函数;最后光学传递函数取模,即为调制传递函数,并以第一个调制传递函数值为基准进行归一化处理,通过坐标变换得到空间任意方向θ上的调制传递函数;将空间任意方向θ上Nyquist频率处的调制传递函数值作为凝视卫星面阵CCD相机的MTF指标评价值。本发明解决了背景技术中的凝视卫星面阵CCD相机MTF无在轨检测的技术问题。
Description
技术领域
本发明涉及凝视卫星面阵CCD相机的MTF在轨检测方法,特别涉及一种地球同步静止轨道卫星面阵相机的MTF检测方法。
背景技术
凝视卫星可实时、定点的对视场内发生的现象进行连续的观测,也可根据用户需求直接定制连续的观测,其灵活连续的观测方式获得比传统对地观测卫星更多的关于“感兴趣目标地区”的动态信息。随着凝视卫星面阵CCD相机的灵敏度和分辨率提高,其获得的连续高精度的地物空间信息、辐射信息和光谱信息将具有极广的应用领域,而这些信息的应用价值更大地取决于凝视卫星面阵CCD数据的成像质量。但由于受到卫星发射震动、空间环境变化以及成像过程中光学、探测器和电子学子系统衰减等因素的影响,凝视卫星面阵CCD相机的成像质量会使逐渐退化和衰变,从而降低了凝视卫星面阵CCD相机数据的应用价值。因而,对凝视卫星面阵CCD相机的成像质量评价是十分必要的。
调制传递函数(MTF)是评价凝视卫星面阵CCD相机数据质量的一个重要物理参数,它既能客观地反映一个成像系统的空间频率特性,又能真实地反映成像系统对不同空间目标的调制传递特性,是国际上通行的评价成像质量的方法。目前,针对线阵CCD相机的MTF已经有对比度法、刃边法和脉冲法等多种在轨检测方法,但是由于面阵CCD光电接收阵列是受空间和探元耦合影响很大的离散采集器件,其光学传递函数(OTF)不满足线性移不变的前提,使得线阵CCD相机较为成熟的MTF在轨测试方法无法直接应用于凝视卫星面阵CCD相机的MTF在轨检测。并且由于凝视卫星面阵CCD相机的遥感应用起步较晚,到目前为止,尚未见可查阅的文献资料中有关这方面的检测方法。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种凝视卫星面阵CCD相机的MTF在轨检测方法,解决了背景技术中的凝视卫星面阵CCD相机MTF无在轨检测的技术问题。
本发明的技术方案是:一种凝视卫星面阵CCD相机的MTF在轨检测方法,步骤如下:
1)建立点源物分布与像分布之间端对端的点扩散函数模型:
式中,x为点源像元在凝视卫星垂轨方向的位置;y为点源像元在凝视卫星沿轨方向的位置;μx为点源像元在x方向的均值;所述的x方向为凝视卫星垂轨方向;μy为点源像元在y方向的均值;所述的y方向为凝视卫星沿轨方向;σx为点源像元在x方向的标准差;σy为点源像元在y方向的标准差;ρxy为点源在x方向和y方向的相关系数;
2)根据凝视卫星CCD相机地面分辨率的大小,基于Nyquist采样定律,在点源布设区进行点源的设计和布设;
21)点源设计:选取黑白方块靶标作为点源材料;所述的白色方块靶标和黑色方块靶标的反射率分别为60%和5%,对比度为12:1;每块黑色方块靶标或白色方块靶标的成像区对应凝视卫星CCD相机的一个像元;
22)点源布设:每个点源由中心为1个白色方块靶标,四周为多个黑色方块靶标的矩形或长方形组成,且每个矩形或长方形点源的靶标数目M×N≥5×5;其中M为矩形或长方形的行靶标数,N为矩形或长方形的列靶标数;点源间的距离大于等于5倍的凝视卫星CCD相机地面分辨率;
3)获取凝视卫星面阵CCD相机过点源区的星地同步成像数据,将该数据转换为表观辐亮度数据,并对点源区每个像元的表观辐亮度数据进行归一化,得到点源区每个像元的点扩散函数值;
4)对步骤3)获得的点扩散函数值进行变换,进一步将点源区每个像元的点扩散函数变换整理为矩阵方程的向量表达式,并通过最小范数求解得到面阵CCD相机的点扩散函数参数,即步骤1)点扩散函数模型中的μx、μy、σx、σy、ρxy;
5)将步骤4)得到的面阵CCD相机点扩散函数参数代入步骤1)中的点扩散函数模型,并对点扩散函数进行傅里叶变换,得到光学传递函数;
6)对步骤5)得到的光学传递函数取模,即为调制传递函数,并以第一个调制传递函数值为基准进行归一化处理,通过坐标变换得到空间任意方向θ上的调制传递函数;将空间任意方向θ上Nyquist频率处的调制传递函数值作为凝视卫星面阵CCD相机的MTF指标评价值。
本发明与现有技术相比的优点在于:
(1)本发明针对凝视卫星面阵CCD相机接收阵列的空间移变特点,建立点源物分布与像分布之间端对端的点扩散函数模型,通过布设星地同步的点源,实现了凝视卫星面阵CCD相机目前无在轨MTF检测的技术问题。
(2)本发明依据点扩散函数的原理,基于凝视卫星面阵CCD相机的成像特点和参数,对点源进行了合理可行性的设计及布设,有效地控制了像元间临近效应的影响,提高了在轨检测的点源精度,确保了该方法更具有很强普适性。
(3)本发明确立的凝视卫星面阵CCD相机二维空间调制传递函数,避免了仅通过凝视卫星面阵CCD相机沿轨方向或垂轨方向建立的一维调制传递函数作为评价指标的单一性和片面性,从而全面地评价了凝视卫星面阵CCD相机的成像质量。
(4)本发明为凝视卫星面阵CCD相机MTF在轨检测提出的方法,兼顾了凝视卫星面阵CCD成像特点与点源靶标布设可实行性,对凝视卫星面阵CCD相机成像质量评价具有很强的工程应用意义。
附图说明
图1为本发明方法流程图。
具体实施方式
(1)点扩散函数确立:依据几何光学在凝视卫星成像中的原理,结合面阵CCD接收阵列的离散特性,考虑光在传播中受大气的影响,建立点源物分布与像分布之间端对端的点扩散函数(PSF)模型;
式中:x为点源像元在凝视卫星垂轨方向的位置;
y为点源像元在凝视卫星沿轨方向的位置;
μx为点源像元在x方向的均值;所述的x方向为凝视卫星垂轨方向;
μy为点源像元在y方向的均值;所述的y方向为凝视卫星沿轨方向;
σx为点源像元在x方向的标准差;
σy为点源像元在y方向的标准差;
ρxy为点源在x方向和y方向的相关系数。
其中,点扩散函数模型中5个待求解参数(μx、μy、σx、σy、ρxy)的大小确定了点源物分布与像分布间的扩散程度。
(2)点源设计及布设:依据凝视卫星CCD相机的地面分辨率(GSD)大小,基于Nyquist采样定律,考虑临近像元影响,进行点源的设计和布设;
(2.1)点源材料选取
考虑到外场星地同步工作的特点,点源的材料选择适应野外环境变化和野外收展与运输方便的常用黑白方块靶标(白色方块靶标和黑色方块靶标的反射率为60%和5%,对比度为12:1)。
(2.2)黑白方块靶标大小
每块黑色方块靶标或白色方块靶标的大小与凝视卫星CCD相机的地面分辨率(GSD)相同,即每块黑色方块靶标或白色方块靶标的成像区对应凝视卫星CCD相机的一个像元。
(2.3)点源设计
每个点源由中心为1个白色方块靶标和四周为多个黑色方块靶标组成的矩形或长方形,考虑到Nyquist采用定律和临近像元影响,要求每个点源的靶标数目M×N≥5×5(中间方块为白色其余为黑色)。
(2.4)点源布设
为了避免单个点源的测量误差,测试中可依据情况布设多个点源;考虑到Nyquist采用定律和临近像元影响,每个点源的间距D≥5×GSD。
(3)点源区每个像元的点扩散函数值提取:读取凝视卫星面阵CCD相机过点源区的星地同步成像数据,对该数据进行处理将其转换为表观辐亮度(即相机入瞳处辐射能量值)数据,并对点源区每个像元的表观辐亮度数据归一化得到点源区每个像元的点扩散函数值;
(3.1)点源区每个像元的表观辐亮度
读取面阵CCD相机的数据,依据数据附带文件中的绝对定标系数,按照公式②计算点源区每个像元的表观辐亮度Lk(x,y):
Lk(x,y)=Gain×DNR,k(x,y)+Bias; (2)
式中:k为第k点源,k∈[1,K](K为点源区总数);
DNR,k(x,y)为第k个点源区内,像元在x,y处的相对辐射校正后的数字DN值;
Gain和Bias为面阵CCD相机的绝对辐射定标系数。
(3.2)点源区每个像元的表观辐亮度归一化
寻找K个点源区中所有像元的最大值记为Lmax,按照公式(3)将其归一化记为
(3.3)点源区每个像元的点扩散函数值
点源区每个像元的表观辐亮度归一化即为点源区每个像元的点扩散函数值,代入方程(1)为:
式中:μk,x为第k个点源像元在x方向的均值;μk,y为第k个点源像元在y方向的均值;σk,x为第k个点源像元在x方向的标准差;σk,y为第k个点源像元在y方向的标准差;ρk,xy为第k个点源像元在x方向和y方向的相关系数。
(4)点扩散函数中参数的确定:对点扩散函数进行变换,进一步将点源区每个像元的点扩散函数变换形式整理为矩阵方程的向量表达式,通过最小范数解得到面阵CCD相机的点扩散函数参数;
(4.1)点扩散函数的变换形式
对方程④两边同时取对数,展开平方项进一步整理为方程(5)式:
式中:
(4.2)点扩散函数的向量表达式
将点源区中的M×N个像元点,利用矩阵方程形式将方程(5)式整理为向量表达式(6):
(4.3)点扩散函数中参数的求解
第k个点源区中参与运算的像元数M×N>5个,即方程(6)为超定方程,为此通过最小范数解方程,解为:
Tk=(Ak TAk)-1Ak TZk;
解出向量Tk后,从中可以得到第k个点源的点扩散函数参数:
对K个点源的点扩散函数参数(μk,x、μk,y、σk,x、σk,y、ρk,xy,k∈[1,K])进行平均,最终得到整个凝视卫星面阵CCD相机的点扩散函数参数:
(5)光学传递函数计算:将得到的面阵CCD相机点扩散函数参数代入点扩散函数(PSF)模型,对点扩散函数进行傅里叶变换,得到光学传递函数(OTF);
式中:fx为点源像元在卫星垂轨方向(x方向)的空间频率;
fy为点源像元在卫星沿轨方向(y方向)的空间频率。
(6)调制传递函数计算:对光学传递函数取模即为调制传递函数,并以第一个调制传递函数值为基准作归一化处理,通过坐标变换得到空间任意方向(θ)上的调制传递函数;将空间任意方向(θ)上Nyquist频率处的调制传递函数值作为凝视卫星面阵CCD相机的MTF指标评价值。
式中:为空间位置(x,y)与y方向的夹角;
为空间位置(x,y)在θ方向上的空间频率。
本发明说明书中未详细说明部分属于本领域技术人员公知常识。
Claims (1)
1.一种凝视卫星面阵CCD相机的MTF在轨检测方法,其特征在于步骤如下:
1)建立点源物分布与像分布之间端对端的点扩散函数模型:
式中,x为点源像元在凝视卫星垂轨方向的位置;y为点源像元在凝视卫星沿轨方向的位置;μx为点源像元在x方向的均值;所述的x方向为凝视卫星垂轨方向;μy为点源像元在y方向的均值;所述的y方向为凝视卫星沿轨方向;σx为点源像元在x方向的标准差;σy为点源像元在y方向的标准差;ρxy为点源在x方向和y方向的相关系数;
2)根据凝视卫星面阵CCD相机地面分辨率的大小,基于Nyquist采样定律,在点源布设区进行点源的设计和布设;
21)点源设计:选取黑白方块靶标作为点源材料;所述的白色方块靶标和黑色方块靶标的反射率分别为60%和5%,对比度为12:1;每块黑色方块靶标或白色方块靶标的成像区对应凝视卫星面阵CCD相机的一个像元;
22)点源布设:每个点源由中心为1个白色方块靶标,四周为多个黑色方块靶标的矩形组成,且每个矩形点源的靶标数目M×N≥5×5;其中M为矩形的行靶标数,N为矩形的列靶标数;点源间的距离大于等于5倍的凝视卫星面阵CCD相机地面分辨率;
3)获取凝视卫星面阵CCD相机过点源区的星地同步成像数据,将该数据转换为表观辐亮度数据,并对点源区每个像元的表观辐亮度数据进行归一化,得到点源区每个像元的点扩散函数值;
4)对步骤3)获得的点扩散函数值进行变换,进一步将点源区每个像元的点扩散函数变换整理为矩阵方程的向量表达式,并通过最小范数求解得到面阵CCD相机的点扩散函数参数,即步骤1)点扩散函数模型中的μx、μy、σx、σy、ρxy;
5)将步骤4)得到的面阵CCD相机点扩散函数参数代入步骤1)中的点扩散函数模型,并对点扩散函数进行傅里叶变换,得到光学传递函数;
6)对步骤5)得到的光学传递函数取模,即为调制传递函数,并以第一个调制传递函数值为基准进行归一化处理,通过坐标变换得到空间任意方向θ上的调制传递函数;将空间任意方向θ上Nyquist频率处的调制传递函数值作为凝视卫星面阵CCD相机的MTF指标评价值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410286149.0A CN104111080B (zh) | 2014-06-24 | 2014-06-24 | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410286149.0A CN104111080B (zh) | 2014-06-24 | 2014-06-24 | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104111080A CN104111080A (zh) | 2014-10-22 |
CN104111080B true CN104111080B (zh) | 2017-05-24 |
Family
ID=51707967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410286149.0A Active CN104111080B (zh) | 2014-06-24 | 2014-06-24 | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104111080B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181150A (zh) * | 2015-09-14 | 2015-12-23 | 中国资源卫星应用中心 | 一种中波红外焦平面阵列探测器的绝对辐射定标方法 |
CN108174186B (zh) * | 2017-12-12 | 2019-08-09 | 北京空间机电研究所 | 一种空间面阵红外相机的点扩散函数测试方法 |
CN111220971B (zh) * | 2019-10-22 | 2021-10-19 | 北京理工大学 | 不受倾斜角影响的高精度测量绝对距离的方法 |
CN114112325B (zh) * | 2021-11-22 | 2023-07-14 | 北京市遥感信息研究所 | 一种自动化高分辨星载相机在轨mtf评价系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281250A (zh) * | 2007-04-04 | 2008-10-08 | 南京理工大学 | 基于图像要素的在轨卫星遥感器调制传递函数的监测方法 |
CN102692273A (zh) * | 2012-05-31 | 2012-09-26 | 中国资源卫星应用中心 | 一种干涉型高光谱成像仪的mtf在轨检测方法 |
-
2014
- 2014-06-24 CN CN201410286149.0A patent/CN104111080B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281250A (zh) * | 2007-04-04 | 2008-10-08 | 南京理工大学 | 基于图像要素的在轨卫星遥感器调制传递函数的监测方法 |
CN102692273A (zh) * | 2012-05-31 | 2012-09-26 | 中国资源卫星应用中心 | 一种干涉型高光谱成像仪的mtf在轨检测方法 |
Non-Patent Citations (4)
Title |
---|
Effects of point-spread function on calibration and radiometric accuracy of CCD camera;Hong Du 等;《APPLIED OPTICS》;20040120;第43卷(第3期);第665-670页 * |
In-flight refocusing and MTF assessment of SPOT5 HRG and HRS cameras;Dominique Leger 等;《Proc.SPIE 4881,Sensors,Systems,and Next-Generation Satellites VI》;20030430;第224-231页 * |
基于线性特征的中分辨率航天遥感器MTF在轨评价及其在图像恢复中的应用;孙吉娟 等;《航天返回与遥感》;20060630;第27卷(第2期);第28-33页 * |
高分辨光学卫星传感器在轨MTF检测;徐伟伟 等;《大气与环境光学学报》;20140331;第9卷(第2期);第97-111页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104111080A (zh) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Duggin et al. | Assumptions implicit in remote sensing data acquisition and analysis | |
CN106896069B (zh) | 一种基于彩色数码相机单幅rgb图像的光谱重建方法 | |
US10274420B2 (en) | Compact multifunctional system for imaging spectroscopy | |
Honkavaara et al. | Processing and assessment of spectrometric, stereoscopic imagery collected using a lightweight UAV spectral camera for precision agriculture | |
Brauers et al. | Multispectral filter-wheel cameras: Geometric distortion model and compensation algorithms | |
Kelcey et al. | Sensor correction and radiometric calibration of a 6-band multispectral imaging sensor for UAV remote sensing | |
Slocum et al. | Simulated imagery rendering workflow for UAS-based photogrammetric 3D reconstruction accuracy assessments | |
CN104111080B (zh) | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 | |
CN108446596A (zh) | 基于可见光相机矩阵的虹膜3d四维数据采集系统及方法 | |
Zheng et al. | Joint image and depth estimation with mask-based lensless cameras | |
CN103905746A (zh) | 亚像素级图像偏移定位及叠加方法和装置以及视频设备 | |
Draper et al. | Gemini planet imager observational calibrations IX: least-squares inversion flux extraction | |
CN105841819A (zh) | 一种有云条件下的地表温度的估算方法和装置 | |
Xu et al. | Hyperspectral image reconstruction based on the fusion of diffracted rotation blurred and clear images | |
Wang et al. | The impact of variable illumination on vegetation indices and evaluation of illumination correction methods on chlorophyll content estimation using UAV imagery | |
CN109087341A (zh) | 一种近距离高光谱相机与测距传感器的融合方法 | |
CN104833425B (zh) | 一种偏振成像系统中偏振片角度的组合方法 | |
Gašparović et al. | Testing of image quality parameters of digital cameras for photogrammetric surveying with unmanned aircrafts | |
Marton et al. | The Herschel/PACS Point Source Catalogue Explanatory Supplement | |
Eskelinen | Computational methods for hyperspectral imaging using Fabry–Perot interferometers and colour cameras | |
Soszyńska et al. | Feasibility study of hyperspectral line-scanning camera imagery for remote sensing purposes | |
De Oliveira et al. | Metrological analysis of the three-dimensional reconstruction based on close-range photogrammetry and the fusion of long-wave infrared and visible-light images | |
Klein | Multispectral imaging and image processing | |
Mengual et al. | The PANIC software system | |
Kleinschmidt et al. | Spatial fusion of different imaging technologies using a virtual multimodal camera |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |