CN106289319B - 一种基于在轨序列图像的卫星***传函测算方法 - Google Patents
一种基于在轨序列图像的卫星***传函测算方法 Download PDFInfo
- Publication number
- CN106289319B CN106289319B CN201610617345.0A CN201610617345A CN106289319B CN 106289319 B CN106289319 B CN 106289319B CN 201610617345 A CN201610617345 A CN 201610617345A CN 106289319 B CN106289319 B CN 106289319B
- Authority
- CN
- China
- Prior art keywords
- target
- satellite
- orbit
- letter
- measuring method
- 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)
- Image Analysis (AREA)
Abstract
本发明一种基于在轨序列图像的卫星***传函测算方法,采集含有脉冲靶标的序列图像,在每帧图像中定位亚像元边缘位置,根据脉冲靶标的灰度分布提取采样点,拟合线扩展函数,对序列统计求平均,再做傅里叶变换得到传函曲线,最后结合靶标像面宽度对传函曲线进行修正。该方法通过连续时间采样弥补空间采样的不足,能够弥补靶标规模不足对卫星传函测算的限制,显著降低成本,基于以上优点,针对高轨遥感卫星分辨率相对低轨卫星较低但是可以实现凝视成像的特点,在高轨遥感卫星相机在轨性能评定中具有广泛的应用前景和重大应用价值。
Description
技术领域
本发明属于卫星性能测试评价领域,涉及一种基于在轨序列图像的卫星***传函测算方法。
背景技术
调制传递函数(MTF)是评价光学卫星成像***性能的重要客观指标,为了定量化评价相机成像质量,发射之前都进行实验室的MTF测试,发射之后还要进行在轨MTF测试。Becker等人2005年在ISPRS文章中总结提出了四种MTF测试方法,第一种是点脉冲法,第二种基于线脉冲法,第三种是刃边法,第四种是辐射状靶标法,黄巧林等人2006年在《航天返回与遥感》上发表的《航天光学遥感器MTF测试技术研究》一文中给出了第五种方法,高对比度矩形靶标测试方法。地球静止轨道遥感卫星可以实现凝视成像,显著提高对同一地球观测的时间分辨率,高分辨对地观测***重大工程专项中的高分四号卫星就是一颗地球静止轨道遥感卫星,与既有太阳同步轨道卫星的观测能力具有很强的互补性。但是同时,由于成像距离远,地球静止轨道遥感卫星分辨率相对太阳同步轨道卫星低,以高分四号卫星为例,如果运行在太阳同步轨道,分辨率可以达到0.7m,但是在地球静止轨道分辨率仅为50m。分辨率降低会给卫星在轨测试带来不利影响,对于传函测试,如果采用常用的刃边法或脉冲法,特征地物面积需要提高一两个数量级,布设靶标的工程量和成本剧增,因此迫切需要研究提出有效的解决方案。综合考虑上述几种方法,脉冲法可以直接从图像中提取线扩展函数,跳过边缘扩展函数拟合的问题,因此比刃边法计算结果精度要高。但是与之相应的是对采样图像的质量要求较高,简单从图像上选取样区困难,有必要通过人工靶标试验获取高质量的样本图像。
本发明解决的技术问题是:克服现有技术的不足,提出一种基于在轨序列图像的卫星***传函测算方法,本发明方法能够仅以约10个像元大小的特征地物测量出分辨率在几十米量级的卫星在轨传函,计算结果能够客观反映卫星性能。
本发明的技术解决方案是:一种基于在轨序列图像的卫星***传函测算方法,其步骤如下:
1)确定用于在轨传函测试的地面均匀背景区,所述的地面均匀背景区反射率为ρ1,空间范围水平方向不小于M个像元GSDx,垂直方向应不小于N个像元GSDy;其中GSDx为卫星水平方向像元分辨率,GSDy为卫星中心视场垂直方向像元分辨率;地面上均匀背景区布设高亮度靶标,靶标反射率为ρ2,靶标形状为矩形,其中短边长度不少于m倍实际地面采样距离,长边长度不少于n倍实际地面采样距离;
2)调整卫星姿态,使成像中心视场对准地面均匀背景区中心位置,连续成像K次形成K帧图像,其中K为正整数;
3)在第一帧图像中提取靶标中垂线处的线扩展函数;
4)重复步骤3)并叠加各帧结果获得整个序列的线扩展函数;
5)将线扩展函数的傅里叶变换与宽度为靶标宽度的方波的傅里叶变换做比值处理,消除实际靶标宽度的影响,得到最终的传函曲线。
所述的K、M、N、m、n均为正整数,且K≥20>M>10>N≥n≥5>m≥2。
所述的ρ1、ρ2均为正实数,且1>ρ2>3*ρ1>0。
本发明与现有技术相比的有益效果是:
本发明方法从理论上解决了分辨率在几十米量级的卫星在轨传函中由于靶标面积要求过大而不可行的问题,通过时间连续采样弥补靶标图像空间样本不足,可以反映卫星图像真实的传函,从而客观评价地球同步轨道遥感卫星成像能力,因此在卫星地面***工程开展在轨测试和在轨运行阶段图像质量评价中具有实际应用价值。
附图说明
图1为本发明方法流程图。
具体实施方式
下面参照附图对本方法进行说明。
如图1所示,一种基于在轨序列图像的卫星***传函测算方法,步骤如下:
1)确定用于在轨传函测试的地面均匀背景区,若卫星中心视场地面分辨率水平方向为50m,垂直方向为50m,则矩形靶标短边长度120m,长边长度300m,靶标短边方向背景长度600m,靶标长边方向背景长度300m;若靶标材料反射率超过0.6,背景反射率应不超过0.2;
2)调整卫星姿态,使成像中心视场对准地面均匀背景区中心位置,连续成像30次,得到30帧图像;
3)在第一帧图像中提取靶标中垂线处的线扩展函数;
4)重复步骤3)并叠加30桢结果获得整个序列的线扩展函数;
5)将线扩展函数的傅里叶变换与宽度为靶标宽度的方波的傅里叶变换做比值处理,消除实际靶标宽度的影响,得到最终的传函曲线。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (3)
1.一种基于在轨序列图像的卫星***传函测算方法,其特征在于步骤如下:
1)确定用于在轨传函测试的地面均匀背景区,所述的地面均匀背景区反射率为ρ1,空间范围水平方向不小于M个像元GSDx,垂直方向应不小于N个像元GSDy;其中GSDx为卫星中心视场水平方向像元分辨率,GSDy为卫星中心视场垂直方向像元分辨率;地面上均匀背景区布设高亮度靶标,靶标反射率为ρ2,靶标形状为矩形,其中短边长度不少于m倍实际地面采样距离,长边长度不少于n倍实际地面采样距离;
2)调整卫星姿态,使成像中心视场对准地面均匀背景区中心位置,连续成像K次形成K帧图像,其中K为正整数;
3)在第一帧图像中提取靶标中垂线处的线扩展函数;
4)按照步骤3)中的方法在除第一帧外的每一帧图像中提取靶标中垂线处的线扩展函数,并叠加各帧结果获得整个序列的线扩展函数;
5)将步骤4)获得的整个序列的线扩展函数的傅里叶变换与宽度为靶标宽度的方波的傅里叶变换做比值处理,消除实际靶标宽度的影响,得到最终的传函曲线。
2.根据权利要求1所述的一种基于在轨序列图像的卫星***传函测算方法,其特征在于:所述的K、M、N、m、n均为正整数,且K≥20>M>10>N≥n≥5>m≥2。
3.根据权利要求1所述的一种基于在轨序列图像的卫星***传函测算方法,其特征在于:所述的ρ1、ρ2均为正实数,且1>ρ2>3*ρ1>0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610617345.0A CN106289319B (zh) | 2016-07-29 | 2016-07-29 | 一种基于在轨序列图像的卫星***传函测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610617345.0A CN106289319B (zh) | 2016-07-29 | 2016-07-29 | 一种基于在轨序列图像的卫星***传函测算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106289319A CN106289319A (zh) | 2017-01-04 |
CN106289319B true CN106289319B (zh) | 2019-06-18 |
Family
ID=57663588
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610617345.0A Active CN106289319B (zh) | 2016-07-29 | 2016-07-29 | 一种基于在轨序列图像的卫星***传函测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106289319B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108470099B (zh) * | 2018-03-15 | 2021-11-02 | 长沙天玖卫星科技有限公司 | 光学成像类小卫星成像能力分析与姿态控制指标分析方法 |
CN114112325B (zh) * | 2021-11-22 | 2023-07-14 | 北京市遥感信息研究所 | 一种自动化高分辨星载相机在轨mtf评价*** |
CN116399820B (zh) * | 2023-06-07 | 2023-08-04 | 中国科学院空天信息创新研究院 | 植被遥感产品真实性校验方法、装置、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833480A (zh) * | 2015-04-27 | 2015-08-12 | 北京空间机电研究所 | 一种测试刃边倾角对调制传递函数影响的装置及方法 |
CN105096319A (zh) * | 2015-09-10 | 2015-11-25 | 北京空间机电研究所 | 一种基于凝视成像的卫星在轨信噪比测试方法 |
-
2016
- 2016-07-29 CN CN201610617345.0A patent/CN106289319B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833480A (zh) * | 2015-04-27 | 2015-08-12 | 北京空间机电研究所 | 一种测试刃边倾角对调制传递函数影响的装置及方法 |
CN105096319A (zh) * | 2015-09-10 | 2015-11-25 | 北京空间机电研究所 | 一种基于凝视成像的卫星在轨信噪比测试方法 |
Non-Patent Citations (1)
Title |
---|
"基于卫星遥感图像的MTF计算和分析";蔡新明;《中国优秀硕士学位论文全文数据库 信息科技辑》;20080115;正文第17-19页 |
Also Published As
Publication number | Publication date |
---|---|
CN106289319A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106767714B (zh) | 提高卫星图像定位精度的等效失配模型多级标校方法 | |
US10578735B2 (en) | Multilook coherent change detection | |
Wang et al. | Cloud base height from sky imager and cloud speed sensor | |
CN105096319B (zh) | 一种基于凝视成像的卫星在轨信噪比测试方法 | |
CN106289319B (zh) | 一种基于在轨序列图像的卫星***传函测算方法 | |
CN103983343B (zh) | 一种基于多光谱影像的卫星平台震颤检测方法及*** | |
GB2554796A (en) | Testing 3D imaging systems | |
CN109712089A (zh) | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 | |
CN105824022A (zh) | 一种电网不良地质体三维形变监测方法 | |
CN101706951B (zh) | 一种基于特征融合的气动光学图像质量客观评价方法、装置及*** | |
CN102855628A (zh) | 多源多时相高分辨率卫星遥感影像自动匹配方法 | |
Kuhn et al. | Determination of the optimal camera distance for cloud height measurements with two all-sky imagers | |
US20130093880A1 (en) | Height Measurement Apparatus And Method | |
CN108416091A (zh) | 一种简易的相机地面分辨率与无人机飞行高度关系的测量方法 | |
CN109064515B (zh) | 一种对月观测数据提取方法及提取装置 | |
CN104111080B (zh) | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 | |
CN103957404A (zh) | 一种星载推扫式相机点扩散函数的实验室测量装置及方法 | |
CN106651928A (zh) | 一种空间面阵相机亚像元点扩散函数测量装置及方法 | |
Drolshagen et al. | Distance determination method of dust particles using Rosetta OSIRIS NAC and WAC data | |
CN106568437B (zh) | 静止卫星成像导航与配准恒星南北参数提取方法及*** | |
Baillarin et al. | Automatic registration of optical images, a stake for future missions: Application to ortho-rectification, time series and mosaic products | |
Zhang et al. | Fusion of ascending and descending polarimetric SAR data for color orthophoto generation | |
Türkmenoğlu et al. | On-orbit modulation transfer function estimation For RASAT | |
CN106528498B (zh) | 恒星遥感影像的点扩散函数提取方法 | |
CN109556715A (zh) | 一种航空高光谱多航带影像辐射畸变校正方法 |
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 |