CN106530301B - 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 - Google Patents
卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 Download PDFInfo
- Publication number
- CN106530301B CN106530301B CN201611084674.XA CN201611084674A CN106530301B CN 106530301 B CN106530301 B CN 106530301B CN 201611084674 A CN201611084674 A CN 201611084674A CN 106530301 B CN106530301 B CN 106530301B
- Authority
- CN
- China
- Prior art keywords
- interference
- satellite
- remote sensing
- noise
- channel
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000003384 imaging method Methods 0.000 title claims abstract description 16
- 230000000007 visual effect Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 238000013139 quantization Methods 0.000 claims description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000005057 refrigeration Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 description 6
- 230000003321 amplification Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000011056 performance test Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
Abstract
本发明提供了一种卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，包括以下步骤：步骤一、在卫星有效载荷观测视场为暗背景条件，卫星正常工作的地面测试环境下，对仪器数据进行实时采集；步骤二、对有效载荷科学源报数据进行包格式解析和判读；步骤三、对图像数据进行频域分析，判读在频域内不存在固定频率的干扰，或者识别出存在固定频率干扰；步骤四、各通道遥感数据以DN值方式整幅图像显示，对明显的噪声干扰进行目视识别；步骤五、逐行逐探元计算探测图像的均值和标准差，以标准差的变化作为空域噪声评价的指标。本发明可以最大限度的保证对各种干扰情况进行快速识别无疏漏，并结合噪声定量评价，确保不受其他主客观因素影响干扰识别结果。
Description
技术领域
本发明涉及航天技术领域，具体地，涉及一种卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法。
背景技术
对地观测卫星一般由平台服务系统和有效载荷系统两大部分组成，有效载荷是卫星效能实现的关键，也是地面测试的重要环节。对于光学有效载荷测试，包括性能测试和接口测试。
性能测试主要考核仪器的定标性能、光学性能、几何性能与用户使用要求的符合性，需要专用的测试设备，一般在卫星测试中只进行功能和光电通路的检查；而接口测试是有效载荷卫星测试重点考核的内容，除常规的机械接口、电接口、热接口和软件通讯接口的匹配性外，更为重要的是电子学的兼容性，即卫星正常工作状态下，其他分系统与载荷可同时正常工作，相互之间不受影响。因此对有效载荷的干扰和噪声情况进行实时快速无疏漏的评价识别是非常重要的卫星测试项目。
图像干扰与噪声是相互关联的，干扰一般为外界引入的具有一定周期性的，其结果可能会引起噪声水平增高，也可能淹没在数据周期性采样中；噪声可以是外界引入，也包括系统内部探测器、电子学固有的噪声。本发明主要评价的是外界引入的干扰和噪声水平的放大。
传统目视评价方法仅能通过肉眼对图像中明显存在的周期性干扰条纹进行识别，其弊端主要在主观性强、受图像显示拉伸甚至显示器分辨率限制、对小干扰或非周期性随机干扰可能漏判等。
发明内容
本发明的目的是针对上述背景技术中的不足，提出一种卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法。
根据本发明的一个方面，提供一种卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，包括以下步骤：
步骤一、在卫星有效载荷观测视场为暗背景条件，卫星正常工作的地面测试环境下，对仪器数据进行实时采集；
步骤二、对有效载荷科学源报数据进行包格式解析和判读，分通道提取并按数据量化位数要求及探测区域顺序排列图像数据；
步骤三、对图像数据进行频域分析，判读在频域内不存在固定频率的干扰，或者识别出存在固定频率干扰，并可精确定位干扰频率位置和幅值；
步骤四、各通道遥感数据以DN值方式整幅图像显示，对明显的噪声干扰进行目视识别；
步骤五、逐行逐探元计算探测图像的均值和标准差，以标准差的变化作为空域噪声评价的指标。
优选地，所述步骤一关于卫星暗背景环境的建立，对于可见通道通过遮挡相机通光孔径方式实现。
优选地，所述步骤一关于卫星暗背景环境的建立，对于红外通道利用真实探测器、液氮制冷方式实现或通过模拟电阻盒替代真实探测器实现。
优选地，所述步骤二的有效载荷科学源报数据解析需分通道提取出光学观测遥感数据内容和辅助信息内容，并根据辅助信息进行按原始探测行列、探元和时间顺序排列。
优选地，所述步骤三的频域分析为将各通道各探测器单元数据按时间轴为序列进行FFT(傅里叶)变换。
优选地，所述步骤四对通道遥感数据进行整幅图全量化位数归一化拉伸后显示。
优选地，所述步骤五采用均值计算公式和标准差计算公式，均值计算公式为如下式：
所述标准差计算公式为如下式：
s为标准差，n为像元总数，xi为第i个像元的DN值，为所有像元的平均值。
优选地，所述步骤五包括以下评价准则：各通道多幅图像计算整幅均值基本稳定；各通道多幅图像计算整幅标准差基本稳定；同一通道各探元多幅图像计算探元内均值基本稳定；同一通道各探元多幅图像计算探元内标准差基本稳定；各通道各探元计算标准差值一般应在20以内。
与现有技术相比，本发明具有如下的有益效果：一，本方法可以最大限度的保证对各种干扰情况进行快速识别无疏漏，并结合噪声定量评价，确保不受其他主客观因素影响干扰识别结果。二，本方法可实现多通道多探元遥感数据并行判读和识别，提高效率的同时，也可对相互之间噪声和干扰来源进行快速定位。本发明方法综合进行频域、空域和目视判读，从而识别和评价遥感图像噪声水平以及是否存在干扰。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法的流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
一般认为，多通道成像类遥感仪器需考核其电子学自兼容性和卫星电子学兼容性，在系统级和验收测试中，需完成其自兼容性进行考核；而在装星后测试中，重点考察与卫星其他分系统电子学兼容性，以及卫星接地状态、电缆屏蔽状态差异可能引发的图像噪声水平放大。
本发明方法包含如下步骤：
步骤一、在卫星有效载荷观测视场为暗背景条件，卫星正常工作的地面测试环境下，对仪器数据进行实时采集；
可见通道通过遮挡相机通光孔径方式实现仪器暗背景成像条件的建立，外通道可使用真实红外探测器，通过液氮对探测器制冷使其正常探测；也可通过模拟电阻盒替代真实探测器实现。经验证，上述方法所建立的暗背景成像条件可真实模拟仪器底噪声环境。
步骤二、对有效载荷科学源报数据进行包格式解析和判读，分通道提取并按数据量化位数要求及探测区域顺序排列图像数据；
原始科学源包数据包括辅助信息部分和遥感数据部分，按照CCSDS协议要求进行有效编排，遥感数据部分排列以时间为基准所有通道混合下传。因此数据解析部分需根据辅助信息中行列信息、对各通道数据进行重新组包，并按照约定的量化位数和通道探元数解析、排列。
步骤三、对图像数据进行频域分析，可判读在频域内不存在固定频率的干扰，或者可识别出存在固定频率干扰，并可精确定位干扰频率位置和幅值；
具体的频域分析方法为，各通道各探测元按时间周序列进行FFT变换，具体变化公式为如下式：
ωN＝e(-2πi)/N
各通道各探测器单元频域数据实时更新显示，如出现异常尖峰，则对特定频率进行识别和分析。
步骤四、各通道遥感数据以DN值方式整幅图像显示，对明显的噪声干扰进行目视识别；
光学遥感数据需按照探测场景的行列顺序，按通道进行整幅图全量化位数归一化拉伸后显示。目视识别中明显的噪声和干扰类型包括横条带(条纹)、竖条纹(条纹)、斜条带(条纹)、多像素图像DN值跳变等。
步骤五、逐行逐探元计算探测图像的均值和标准差，以标准差的变化作为空域噪声评价的指标。
对全量化位数的遥感数据，分别按照整幅图像、逐像元图像、逐行图像计算图像DN值的均值和标准差。均值计算公式为：
标准差计算公式为：
s为标准差，n为像元总数，xi为第i个像元的DN值，为所有像元的平均值。
计算结果对噪声水平的定量化评价标准包括：各通道多幅图像计算整幅均值基本稳定；各通道多幅图像计算整幅标准差基本稳定；同一通道各探元多幅图像计算探元内均值基本稳定；同一通道各探元多幅图像计算探元内标准差基本稳定；各通道各探元计算标准差值一般应在20以内。根据上述评价标准，可对连续多幅图像空域噪声水平进行计算评估，并可与分系统测试结果进行对比，如对比结果存在较大差异，则可以考虑对干扰和噪声进行排查。
频域干扰监测结果、目视判读结果和空域噪声水平定量计算结果可联合作为多通道成像类遥感仪器是否在卫星测试中受到干扰的判据。如出现干扰情况，排查过程可以采用星上各分系统分别加电、物理屏蔽措施、更改接地方式等，分阶段采用本发明方法评估并定位上述干扰和噪声来源。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。
Claims (8)
1.一种卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，包括以下步骤：
步骤一、在卫星有效载荷观测视场为暗背景条件，卫星正常工作的地面测试环境下，对仪器数据进行实时采集；
步骤二、对有效载荷科学源报数据进行包格式解析和判读，分通道提取并按数据量化位数要求及探测区域顺序排列图像数据；
步骤三、对图像数据进行频域分析，判读在频域内不存在固定频率的干扰，或者识别出存在固定频率干扰，并精确定位干扰频率位置和幅值；
步骤四、各通道遥感数据以DN值方式整幅图像显示，对明显的噪声干扰进行目视识别；
步骤五、逐行逐探元计算探测图像的均值和标准差，以标准差的变化作为空域噪声评价的指标。
2.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤一关于卫星暗背景环境的建立，对于可见通道通过遮挡相机通光孔径方式实现。
3.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤一关于卫星暗背景环境的建立，对于红外通道利用真实探测器、液氮制冷方式实现或通过模拟电阻盒替代真实探测器实现。
4.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤二的有效载荷科学源报数据解析需分通道提取出光学观测遥感数据内容和辅助信息内容，并根据辅助信息进行按原始探测行列、探元和时间顺序排列。
5.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤三的频域分析为将各通道各探测器单元数据按时间轴为序列进行FFT变换。
6.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤四对通道遥感数据进行整幅图全量化位数归一化拉伸后显示。
7.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤五采用均值计算公式和标准差计算公式，均值计算公式为如下式：
所述标准差计算公式为如下式：
s为标准差，n为像元总数，xi为第i个像元的DN值，为所有像元的平均值。
8.根据权利要求1所述的卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法，其特征在于，所述步骤五包括以下评价准则：各通道多幅图像计算整幅均值基本稳定；各通道多幅图像计算整幅标准差基本稳定；同一通道各探元多幅图像计算探元内均值基本稳定；同一通道各探元多幅图像计算探元内标准差基本稳定；各通道各探元计算标准差值应在20以内。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611084674.XA CN106530301B (zh) | 2016-11-30 | 2016-11-30 | 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611084674.XA CN106530301B (zh) | 2016-11-30 | 2016-11-30 | 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106530301A CN106530301A (zh) | 2017-03-22 |
| CN106530301B true CN106530301B (zh) | 2019-06-07 |
Family
ID=58355374
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611084674.XA Active CN106530301B (zh) | 2016-11-30 | 2016-11-30 | 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106530301B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111536991B (zh) * | 2020-04-26 | 2022-09-27 | 哈尔滨工业大学 | 遥感卫星多模式成像载荷模拟装置 |
| CN116486277A (zh) * | 2023-03-28 | 2023-07-25 | 武汉卓目科技有限公司 | 遥感载荷图像判读方法及遥感载荷快视系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010148269A3 (en) * | 2009-06-18 | 2011-02-10 | Smith International, Inc. | Cyclic noise removal in borehole imaging |
| CN103279950A (zh) * | 2013-05-14 | 2013-09-04 | 中国科学院光电研究院 | 一种基于行列噪声标准差的遥感图像信噪比评估方法 |
| CN104539384A (zh) * | 2014-11-19 | 2015-04-22 | 华中科技大学 | 一种基于卫星被动微波遥感数据的射频干扰检测方法 |
| CN104680541A (zh) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | 基于相位一致性的遥感图像质量评价方法 |
| CN105096319A (zh) * | 2015-09-10 | 2015-11-25 | 北京空间机电研究所 | 一种基于凝视成像的卫星在轨信噪比测试方法 |
| CN105427349A (zh) * | 2015-11-18 | 2016-03-23 | 国家海洋局第一海洋研究所 | 遥感图像条带噪声定量评价方法 |
-
2016
- 2016-11-30 CN CN201611084674.XA patent/CN106530301B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010148269A3 (en) * | 2009-06-18 | 2011-02-10 | Smith International, Inc. | Cyclic noise removal in borehole imaging |
| CN103279950A (zh) * | 2013-05-14 | 2013-09-04 | 中国科学院光电研究院 | 一种基于行列噪声标准差的遥感图像信噪比评估方法 |
| CN104539384A (zh) * | 2014-11-19 | 2015-04-22 | 华中科技大学 | 一种基于卫星被动微波遥感数据的射频干扰检测方法 |
| CN104680541A (zh) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | 基于相位一致性的遥感图像质量评价方法 |
| CN105096319A (zh) * | 2015-09-10 | 2015-11-25 | 北京空间机电研究所 | 一种基于凝视成像的卫星在轨信噪比测试方法 |
| CN105427349A (zh) * | 2015-11-18 | 2016-03-23 | 国家海洋局第一海洋研究所 | 遥感图像条带噪声定量评价方法 |
Non-Patent Citations (2)
| Title |
|---|
| A New method for noise estimation in single-band remote sensing images;Peng Fu et al.;《2012 9th International Conference on Fuzzy Systems and Knowledge Discovery》;20121231;第1664-1668页 |
| 一种中波红外图像条纹处理方法;陈强 等;《遥感学报》;20121231;第16卷(第1期);第93-104页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106530301A (zh) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110617873B (zh) | 电缆的振动检测方法及相关产品 | |
| CN103499774B (zh) | 输电绝缘子闪络在线监视安全预警装置 | |
| CN102461192B (zh) | 三维显示器的光学特性评价装置和光学特性评价方法 | |
| US11856182B2 (en) | Icon-based home certification, in-home leakage testing, and antenna matching pad | |
| CN105100789A (zh) | 一种视频质量评价方法 | |
| CN106530301B (zh) | 卫星地面测试中成像类定量遥感仪器干扰和噪声评价方法 | |
| CN109389589A (zh) | 用于统计人数的方法和装置 | |
| CN106771924A (zh) | 一种利用光电场传感器检测绝缘子缺陷的检测系统及方法 | |
| CN104729718B (zh) | 一种用于红外成像系统netd的处理系统及方法 | |
| KR20140017733A (ko) | 지아이에스 매핑 정보를 이용한 불량 전력설비 검출 방법 및 시스템 | |
| CN111695512B (zh) | 一种无人值守文物监测方法及装置 | |
| CN103412246A (zh) | 基于紫外视频的电力设备异常放电自动检测方法 | |
| CN108881896A (zh) | 抓拍机的测试方法、系统以及存储介质 | |
| CN115453286A (zh) | Gis局部放电诊断方法、模型训练方法、装置及系统 | |
| CN103200349B (zh) | 一种扫描图像色偏自动检测方法 | |
| Montanya et al. | Potential use of space-based lightning detection in electric power systems | |
| CN110148176A (zh) | 红外高光谱多像元载荷的定位一致性的检验方法和装置 | |
| CN113762161A (zh) | 一种障碍物智能监测方法及系统 | |
| CN113391776A (zh) | 物品信息显示方法、装置、电子装置和存储介质 | |
| CN106204531B (zh) | 一种同步检测彩色图像中噪点与边缘点的方法 | |
| Leka et al. | On Identifying and Mitigating Bias in Inferred Measurements for Solar Vector Magnetic-Field Data | |
| CN112837286A (zh) | 一种用于支架变形的检测方法和系统 | |
| KR101743595B1 (ko) | 부분 방전 진단 방법 및 시스템, 이를 이용한 몰드 변압기 열화 모니터링 시스템. | |
| CN105044471A (zh) | 计算机数字视频信息电磁泄漏检测方法 | |
| Hernandez-Mejia | Time Domain Reflectometry (TDR) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| 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 |