CN106600646A - 一种红外摆扫相机在轨图像均匀性校正方法 - Google Patents

一种红外摆扫相机在轨图像均匀性校正方法 Download PDF

Info

Publication number
CN106600646A
CN106600646A CN201611066858.3A CN201611066858A CN106600646A CN 106600646 A CN106600646 A CN 106600646A CN 201611066858 A CN201611066858 A CN 201611066858A CN 106600646 A CN106600646 A CN 106600646A
Authority
CN
China
Prior art keywords
overbar
image
sigma
centerdot
value
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.)
Granted
Application number
CN201611066858.3A
Other languages
English (en)
Other versions
CN106600646B (zh
Inventor
曹世翔
李岩
晋利宾
聂云松
何红艳
江澄
周楠
张炳先
李方琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Space Research Mechanical and Electricity
Original Assignee
Beijing Institute of Space Research Mechanical and Electricity
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Space Research Mechanical and Electricity filed Critical Beijing Institute of Space Research Mechanical and Electricity
Priority to CN201611066858.3A priority Critical patent/CN106600646B/zh
Publication of CN106600646A publication Critical patent/CN106600646A/zh
Application granted granted Critical
Publication of CN106600646B publication Critical patent/CN106600646B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/80Calibration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077Imaging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radiation Pyrometers (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明公开了一种红外摆扫相机在轨图像均匀性校正方法,针对卫星入轨运行特点,获取星上黑体定标数据及多景地物图像后,以星上定标得到的相对辐射及绝对辐射校正系数为主,采用在轨统计的方法修正摆扫相机像元间存在的响应差异,代替传统均匀场图像修正。该校正方法不要求覆盖足够大的均匀地物进行高低温成像,同时对成像气候条件具有很好的鲁棒性,克服了在轨初期快速参数调整与传统校正方法效率低的难题,且取得的均匀性处理效果更好。本发明适用于各种红外摆扫相机,亦可推广到线阵推扫红外成像,满足遥感产品生产及应用需求。

Description

一种红外摆扫相机在轨图像均匀性校正方法
技术领域
本发明涉及一种红外摆扫相机在轨图像均匀性校正方法,属于图像处理技术领域。
背景技术
红外成像技术是通过探测器探测目标物体自身热辐射或获得目标图像,通过该技术能反映出可见光波之外的目标信息。另外,红外成像技术具有隐蔽性好、作用距离远、抗干扰性强、成像质量好、全天候使用等明显优势,在过去六十年间得到了长足发展。
红外探测器的结构形式分为光机扫描成像探测器,线阵扫描结构探测器和阵列探测器。应用在遥感对地成像时,考虑到成像幅宽的要求以及长线阵与大型面阵器件的工艺难度,一般采用限长线阵的红外探测器并设计摆扫机构实现大幅宽对地成像,如图1所示,称为红外摆扫相机。
由于制造工艺的限制,红外摆扫相机各像元响应度几乎都不一样,导致了叠加在图像上的非均匀性噪声,严重影响了系统成像质量,因此非均匀校正是最为关键的红外相机图像处理技术。
目前广泛应用于实践的非均匀性校正方法包括一点法和两点法。其中,两点法综合考虑了目标温度的高端及低端,处理效果明显。因此,为了量化不同像元间的差异,大多数遥感红外相机都配备了黑体定标装置,包括两个不同温度的均匀黑体。在对地成像前,先进行内定标,得到体现内定标非均匀状态的黑体成像DN(数字计数)值图像,计算相应的相对定标系数及绝对定标系数。利用地面真空环境下得到内定标光路到入瞳处能量的推算关系,反演入瞳前的辐亮度(温度)数据。
入轨后,由于器件本身的性能变化,且大气等传输媒介对能量的影响,采用高低温黑体直接计算得到的定标系数处理对地成像时,DN值产品和辐亮度产品残留非均匀性较大,直接影响对目标的判断。用户迫切期望能够高效地获取到DN值产品和辐亮度产品的稳定非均匀性状态,并且缩短调试周期,对成像计划的需求也尽可能少。而传统处理方式要求用户在入轨初期多次协调成像计划,得到不同均匀场的数据,从而根据高低温在轨均匀地物(子块)数据形成修正系数,然后处理其他景数据得到DN值产品和辐亮度产品。在实际操作过程中存在许多缺陷,具体地:
1)红外器件分辨率通常较低,因此要求均匀地物的范围巨大,难以寻找到合适的温度靶场;即使存在满足条件的温度靶场,安排成像时也有可能受天气影响(云雾遮挡等)而造成数据无法使用;
2)要求用户额外增加成像计划,且若天气(如多云)原因导致均匀场获取失败,仍需再安排成像,影响在轨调试进度;
3)均匀场数据中挑选的子块计算得到的非均匀性统计特性与实际情况存在差异,导致处理后全景均匀性精度不高;
4)在轨维护阶段,若由于器件特性漂移,做再次修正时,均匀场时相条件发生变化而可能无法使用,需再次寻找适合的场景,效率低下。
发明内容
本发明解决的技术问题是:克服现有技术的不足,本发明提供了一种红外摆扫相机在轨图像均匀性校正方法,针对红外摆扫相机采用内定标黑体进行校正后残留的非均匀性问题,结合在轨工作条件,克服现有依赖大面积均匀场统计的非均匀性校正方法,提出场景统计方法,仅约束景内统计量偏差,降低对均匀场的依赖,对天气条件等外界影响鲁棒,满足快速性和有效性的要求。
本发明的技术解决方案是:
一种红外摆扫相机在轨图像均匀性校正方法,包括如下步骤:
第一步,从相同成像条件下获取的图像中选择符合修正用的DN值低端原始图像IL和DN值高端原始图像IH
第二步,根据星上高端温度TH和低端温度TL以及对应温度下的黑体成像数据BH和BL,按两点法计算对应的星上相对辐射校正系数KRi、CRi和绝对辐射校正系数KAi、CAi、ΔDNAi,若DN值低端原始图像IL和DN值高端原始图像IH同轨,则共用一组系数,否则需计算各自对应的系数;
第三步,利用星上定标系数,对DN值低端原始图像IL和DN值高端原始图像IH进行一次校正,得到DN值图像ILDN和IHDN,按单扫的像元个数n及摆扫方向分别对DN值图像ILDN和IHDN进行叠加,形成景内统计数据块ILMDN和IHMDN;同理得到一次校正辐亮度图像ILLe和IHLe,并得到叠加数据块ILMLe和IHMLe
第四步,对叠加数据块ILMLe和IHMLe分别按两点法计算DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数
第五步,对典型温度下获取的包含均匀景物的原始图像G,计算对应轨的星上相对辐射校正系数K′Ri、C′Ri和绝对辐射校正系数K′Ai、C′Ai及ΔDN′Ai;利用DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数计算最终用于DN值产品DNG的校正系数KRFi、CRFi和辐亮度产品LeG的校正系数KAFi、CAFi和ΔDNAFi,并生成相应的产品;
第六步,对DN值产品DNG和辐亮度产品LeG进行均匀性检验,若不满足要求,则重新回到第一步并选择其他图像生成修正系数;
第七步,若摆扫相机正反扫均成像,则第三步至第六步需按扫的方向独立计算。
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第一步中,符合修正用的DN值低端图像IL和DN值高端图像IH挑选原则如下:
星上定标黑体的高端温度TH和低端温度TL为红外相机设计的温度线性范围,用于在控制均匀性校正精度的同时降低数据寻找难度,DN值低端图像IL的直方图峰值应在BL的灰度均值±20%之间;同理,DN值高端图像IH的直方图峰值应在BH灰度均值±20%之间。
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第二步中,计算星上辐射校正系数方法如下:
假定红外线阵探测器由n个探测像元组成,星上定标图像摆扫宽为W,对黑体成像数据BH和BL
i=1,2,…,n
其中,BL(i,j)表示BL第i个探元第j列的DN值,绝对定标系数的两点法:
其中,LeT为黑体温度,即绝对温度为T时的入瞳等效光谱辐亮度;L(λ,T)为黑体温度为T时的光谱辐亮度;R(λ)为红外相机的相对光谱响应;ε为黑体发射率;r(λ)为光学系统效率;λ1、λ2是起始和截止波长;ΔDNi为变温修正量,卫星发射前可测定ΔDNi与温度T的关系。
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第三步中,一次校正及数据叠加方法如下:
设DN值低端原始图像IL和DN值高端原始图像IH包含m扫数据,对每一扫数据:
s=1,2,…,m
按扫进行叠加,形成统计上均匀的数据块:
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第四步中,DN值数据和辐亮度数据的均匀性修正系数计算方法如下:
DN值数据均匀性修正系数满足:
i=1,2,…,n
同理,辐亮度数据的均匀性修正系数满足:
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第五步中,均匀景物的原始图像G生成DN值产品DNG和辐亮度LeG的方法如下:
按第二步所述方法计算图像G对应轨的星上相对辐射校正系数K′Ri、C′Ri和绝对辐射校正系数K′Ai、C′Ai及ΔDN′Ai
DN值产品DNG的校正系数KRFi、CRFi满足:
辐亮度产品LeG的校正系数KAFi、CAFi和ΔDNAi满足:
最后,根据第三步所述方法得到DN值产品DNG、辐亮度产品LeG
在上述的一种红外摆扫相机在轨图像均匀性校正方法中,所述第六步中,检验DN值产品DNG、辐亮度LeG均匀性精度的方法如下:
截取DNG中任意大小均匀的图像子块SUBDNG,计算图像子块SUBDNG的最大值DNmax、最小值DNmin以及均值DNmean表示DNG的均匀性水平;同理表示辐亮度产品LeG的均匀性水平。
本发明与现有技术相比的有益效果是:
【1】本发明针对卫星入轨运行特点,获取星上黑体定标数据及多景地物图像后,以星上定标得到的相对辐射及绝对辐射校正系数为主,采用在轨统计的方法修正摆扫相机像元间存在的响应差异,代替传统均匀场图像修正,不要求图像内有足够大的均匀地物。
【2】本发明在轨修正时,不严格要求成像气候条件,具有很好的鲁棒性,满足了有效性的设计需求。
【3】本发明解决了在轨初期快速参数调整与传统校正方法效率低的矛盾,且其取得的均匀性效果更好,维护性强。
【4】本发明步骤精简、方式多样、易于实施,兼顾了红外摆扫相机正反扫均成像等多种情况,具有广阔的应用前景。
【5】本发明逻辑通顺、思路清晰、设计合理,本领域技术人员按照本发明的步骤进行实验,能够快速实现红外摆扫相机在轨图像的均匀性校正,显著提升了工作效率。
附图说明
图1为摆扫相机工作示意图
图2为本发明流程图
图3为本发明同扫数据叠加示意图
具体实施方式
下面结合附图说明和具体实施例对本发明作进一步描述:
根据图1所示的摆扫相机成像原理,由于单景内同类型地物温度整体是均匀的,因此可以将景内的各扫数据进行叠加形成统计意义上的均匀景物。
由于扫间数据叠加得到的是探元在不同温度下的成像统计,因此它更能反映成像的整体效应,从而脱离了严格地物依赖。业务流程中,要求采用的均匀性校正图像覆盖高低端温度,因此形成的均匀性校正系数在一段时间内是稳定可靠的,能够适应于不同阶段的生产实际。
本发明采用星上辐射校正和入轨均匀性修正的处理模式,满足在轨后的DN值产品和辐亮度产品的均匀性需求。
以下结合说明书附图详细说明本发明的技术方案:
如图2所示,一种红外摆扫相机在轨图像均匀性校正方法,包括如下步骤:
第一步,从入轨后对高端温度TH和低端温度TL下的星上定标图像BH和BL,计算均值作为温度对应的DN值水平。为控制均匀性校正精度同时尽可能降低数据寻找的难度,挑选满足以下条件的对地成像数据IH和IL作为均匀性校正图像,所有参数的下角标中,含有H的均对应高端温度TH成像,含有L的均对应低端温度TL成像。
IL图像的直方图峰值应在BL的灰度均值±20%之间;同理,IH图像的直方图峰值应在BH灰度均值±20%之间;
第二步,根据星上高端温度TH和低端温度TL以及对应温度下的黑体成像数据BH和BL,按两点法计算对应的星上相对辐射校正系数KRi、CRi和绝对辐射校正系数KAi、CAi及ΔDNAi,若DN值低端原始图像IL和DN值高端原始图像IH同轨,则共用一组系数。否则需计算各自对应的系数。
假定红外线阵探测器由n个探测像元组成,星上定标图像摆扫宽为W,对BH和BL
i=1,2,…,n
其中,BL(i,j)表示BL第i个探元第j列的DN值。绝对定标系数的两点法
其中,LeT为黑体温度(绝对温度)为T时的入瞳等效光谱辐亮度;L(λ,T)为黑体温度为T时的光谱辐亮度;R(λ)为红外相机的相对光谱响应;ε黑体发射率;r(λ)为光学系统效率;λ1、λ2是起始和截止波长;ΔDNi为变温修正量,卫星发射前可测定ΔDNi与温度T的关系;
第三步,利用星上定标系数,对DN值低端原始图像IL和DN值高端原始图像IH进行一次校正,得到DN值图像ILDN和IHDN,按单扫的像元个数n及摆扫方向分别对DN值图像ILDN和IHDN进行叠加,形成景内统计数据块ILMDN和IHMDN;同理得到一次校正辐亮度图像ILLe和IHLe,并得到叠加数据块ILMLe和IHMLe
设DN值图像IL和IH包含m扫数据,对每一扫数据:
s=1,2,…,m
按扫进行叠加,形成统计上均匀的数据块,如附图3所示:
第四步,对叠加数据块,分别按两点法计算DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数
DN值数据的均匀性修正系数满足:
i=1,2,…,n
同理,辐亮度数据的均匀性修正系数满足:
第五步,对典型温度下获取的包含均匀景物(如附图3中的曲线所围的均匀地物)的原始图像G,计算其对应轨的星上相对辐射校正系数K′Ri、C′Ri和绝对辐射校正系数K′Ai、C′Ai及ΔDN′Ai。利用DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数计算最终用于DN值产品DNG的校正系数KRFi、CRFi和辐亮度产品LeG的校正系数KAFi、CAFi和ΔDNAFi,并生成相应的产品。
DN值产品DNG的校正系数KRFi、CRFi满足:
辐亮度产品LeG的校正系数KAFi、CAFi和ΔDNAi满足:
最后,根据第三步所述方法得到DN值产品DNG、辐亮度产品LeG
第六步,对DN值产品DNG和辐亮度产品LeG进行均匀性检验,如不满足要求,则重新回到第一步并选择其他图像生成修正系数。
具体为:截取DNG中任意大小均匀的图像子块SUBDNG,计算SUBDNG的最大值DNmax、最小值DNmin以及均值DNmean表示DNG的均匀性水平;
同理表示辐亮度产品LeG的均匀性水平;
第七步,若摆扫相机正反扫均成像,则第三步到第六步需按扫的方向独立计算。如图3所示,图中S0和S1表示的正反扫均成像,则需要单独计算。
对本实施例的分析:均匀性校正后的DN值产品DNG和辐亮度产品LeG目标可分辨能力、信噪比等各项质量指标均有大幅度提高,均匀性精度达到5%以内的技术水平。
该发明技术内容适用于各种红外摆扫相机,亦可推广到线阵推扫红外成像,满足遥感产品生产及应用需求。
本发明说明书中未详细描述的内容为本领域技术人员公知技术。

Claims (7)

1.一种红外摆扫相机在轨图像均匀性校正方法,其特征在于,包括如下步骤:
第一步,从相同成像条件下获取的图像中选择符合修正用的DN值低端原始图像IL和DN值高端原始图像IH
第二步,根据星上高端温度TH和低端温度TL以及对应温度下的黑体成像数据BH和BL,按两点法计算对应的星上相对辐射校正系数KRi、CRi和绝对辐射校正系数KAi、CAi、ΔDNAi,若DN值低端原始图像IL和DN值高端原始图像IH同轨,则共用一组系数,否则需计算各自对应的系数;
第三步,利用星上定标系数,对DN值低端原始图像IL和DN值高端原始图像IH进行一次校正,得到DN值图像ILDN和IHDN,按单扫的像元个数n及摆扫方向分别对DN值图像ILDN和IHDN进行叠加,形成景内统计数据块ILMDN和IHMDN;同理得到一次校正辐亮度图像ILLe和IHLe,并得到叠加数据块ILMLe和IHMLe
第四步,对叠加数据块ILMLe和IHMLe分别按两点法计算DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数
第五步,对典型温度下获取的包含均匀景物的原始图像G,计算对应轨的星上相对辐射校正系数K′Ri、C′Ri和绝对辐射校正系数K′Ai、C′Ai及ΔDN′Ai;利用DN值数据的均匀性修正系数及辐亮度数据的均匀性修正系数计算最终用于DN值产品DNG的校正系数KRFi、CRFi和辐亮度产品LeG的校正系数KAFi、CAFi和ΔDHAFi,并生成相应的产品;
第六步,对DN值产品DNG和辐亮度产品LeG进行均匀性检验,若不满足要求,则重新回到第一步并选择其他图像生成修正系数;
第七步,若摆扫相机正反扫均成像,则第三步至第六步需按扫的方向独立计算。
2.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第一步中,符合修正用的DN值低端图像IL和DN值高端图像IH挑选原则如下:
星上定标黑体的高端温度TH和低端温度TL为红外相机设计的温度线性范围,用于在控制均匀性校正精度的同时降低数据寻找难度,DN值低端图像IL的直方图峰值应在BL的灰度均值±20%之间;同理,DN值高端图像IH的直方图峰值应在BH灰度均值±20%之间。
3.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第二步中,计算星上辐射校正系数方法如下:
假定红外线阵探测器由n个探测像元组成,星上定标图像摆扫宽为W,对黑体成像数据BH和BL
D N ‾ l = K R i · B L ( i ) ‾ + C R i , D N ‾ h = K R i · B H ( i ) ‾ + C R i
D N ‾ l = 1 n · W Σ i = 1 n Σ j = 1 W B L ( i , j ) , B L ( i ) ‾ = 1 W Σ j = 1 W B L ( i , j )
D N ‾ h = 1 n · W Σ i = 1 n Σ j = 1 W B H ( i , j ) , B H ( i ) ‾ = 1 W Σ j = 1 W B H ( i , j )
i=1,2,…,n
其中,BL(i,j)表示BL第i个探元第j列的DN值,绝对定标系数的两点法:
B L ( i ) ‾ = K A i · Le T L + C A i + ΔDN i
B H ( i ) ‾ = K A i · Le T H + C A i + ΔDN i
Le T = ϵ ∫ λ 1 λ 2 r ( λ ) L ( λ , T ) R ( λ ) d λ ∫ λ 1 λ 2 R ( λ ) d λ
其中,LeT为黑体温度,即绝对温度为T时的入瞳等效光谱辐亮度;L(λ,T)为黑体温度为T时的光谱辐亮度;R(λ)为红外相机的相对光谱响应;ε为黑体发射率;r(λ)为光学系统效率;λ1、λ2是起始和截止波长;ΔDNi为变温修正量,卫星发射前可测定ΔDNi与温度T的关系。
4.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第三步中,一次校正及数据叠加方法如下:
设DN值低端原始图像IL和DN值高端原始图像IH包含m扫数据,对每一扫数据:
I LDN s = I L s · K R i + C R i , I HDN s = I H s · K R i + C R i
I LLe s = I L s - ΔDN i - C A i K A i , I HLe s = I H s - ΔDN i - C A i K A i
s=1,2,…,m
按扫进行叠加,形成统计上均匀的数据块:
I L M D N = Σ s = 1 m I LDN s , I H M D N = Σ s = 1 m I HDN s
I L M L e = Σ s = 1 m I LLe s , I H M L e = Σ s = 1 m I HLe s .
5.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第四步中,DN值数据和辐亮度数据的均匀性修正系数计算方法如下:
DN值数据均匀性修正系数满足:
D N ‾ L = K R i C · I L M D N ( i ) ‾ + C R i C , D N ‾ H = K R i C · I H M D N ( i ) ‾ + C R i C
D N ‾ L = 1 n · W Σ i = 1 n Σ j = 1 W I L M D N ( i , j ) , I L M D N ( i ) ‾ = 1 W Σ j = 1 W I L M D N ( i , j )
D N ‾ H = 1 n · W Σ i = 1 n Σ j = 1 W I H M D N ( i , j ) , I H M D N ( i ) ‾ = 1 W Σ j = 1 W I H M D N ( i , j )
i=1,2,…,n
同理,辐亮度数据的均匀性修正系数满足:
L e ‾ L = K A i C · I L M L e ( i ) ‾ + C A i C , L e ‾ H = K A i C · I H M L e ( i ) ‾ + C A i C
L e ‾ L = 1 n · W Σ i = 1 n Σ j = 1 W I L M L e ( i , j ) , I L M L e ( i ) ‾ = 1 W Σ j = 1 W I L M L e ( i , j )
L e ‾ H = 1 n · W Σ i = 1 n Σ j = 1 W I H M L e ( i , j ) , I H M L e ( i ) ‾ = 1 W Σ j = 1 W I H M L e ( i , j ) .
6.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第五步中,均匀景物的原始图像G生成DN值产品DNG和辐亮度产品LeG的方法如下:
按第二步所述方法计算图像G对应轨的星上相对辐射校正系数K′Ri、C′Ri和绝对辐射校正系数K′Ai、C′Ai及ΔDN′Ai
DN值产品DNG的校正系数KRFi、CRFi满足:
K R F i = K R i C · K R i ′ , C R F i = K R i C · C R i ′ + C R i C
辐亮度产品LeG的校正系数KAFi、CAFi和ΔDNAi满足:
K A F i = K A i ′ K A i C , C A F i = C A i ′ - K A i ′ K A i C · C A i C , ΔDN A F i = ΔDN A i ′
最后,根据第三步所述方法得到DN值产品DNG、辐亮度产品LeG
7.根据权利要求1所述的一种红外摆扫相机在轨图像均匀性校正方法,其特征在于:所述第六步中,检验产品DNG、LeG均匀性精度的方法如下:
截取DNG中任意大小均匀的图像子块SUBDNG,计算图像子块SUBDNG的最大值DNmax、最小值DNmin以及均值DNmean表示DNG的均匀性水平;同理表示辐亮度产品LeG的均匀性水平。
CN201611066858.3A 2016-11-25 2016-11-25 一种红外摆扫相机在轨图像均匀性校正方法 Active CN106600646B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611066858.3A CN106600646B (zh) 2016-11-25 2016-11-25 一种红外摆扫相机在轨图像均匀性校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611066858.3A CN106600646B (zh) 2016-11-25 2016-11-25 一种红外摆扫相机在轨图像均匀性校正方法

Publications (2)

Publication Number Publication Date
CN106600646A true CN106600646A (zh) 2017-04-26
CN106600646B CN106600646B (zh) 2019-05-24

Family

ID=58595350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611066858.3A Active CN106600646B (zh) 2016-11-25 2016-11-25 一种红外摆扫相机在轨图像均匀性校正方法

Country Status (1)

Country Link
CN (1) CN106600646B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107796518A (zh) * 2017-09-22 2018-03-13 上海卫星工程研究所 线阵扫描相机对复杂地球背景成像试验方法
CN108873924A (zh) * 2018-07-26 2018-11-23 中国人民解放军陆军工程大学 一种机载视频对地摆扫面积计算方法
CN110906906A (zh) * 2019-12-02 2020-03-24 王睿 一种增加遥感图像测绘精度方法
CN111257238A (zh) * 2020-02-05 2020-06-09 国家卫星气象中心(国家空间天气监测预警中心) 一种基于星载太阳漫射板的探测元件间相对定标方法
CN112526506A (zh) * 2020-11-17 2021-03-19 中国科学院长春光学精密机械与物理研究所 目标搜索跟踪方法及目标追踪装置
CN114252161A (zh) * 2021-11-29 2022-03-29 湖北久之洋信息科技有限公司 一种基于线列红外探测器的成像方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101852623A (zh) * 2010-06-02 2010-10-06 中国资源卫星应用中心 一种卫星光学遥感相机内方元素在轨检校方法
CN102622739A (zh) * 2012-03-30 2012-08-01 中国科学院光电技术研究所 一种Bayer滤波阵列彩色相机图像非均匀性校正方法
CN103033270A (zh) * 2012-12-17 2013-04-10 无锡艾立德智能科技有限公司 一种红外热像仪非均匀校正系数生成及控制方法
CN104580944A (zh) * 2013-10-10 2015-04-29 中国科学院光电研究院 对ccd图像进行相对辐射校正的方法
CN105160631A (zh) * 2015-07-02 2015-12-16 山东大学 一种求辐射校正系数的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101852623A (zh) * 2010-06-02 2010-10-06 中国资源卫星应用中心 一种卫星光学遥感相机内方元素在轨检校方法
CN102622739A (zh) * 2012-03-30 2012-08-01 中国科学院光电技术研究所 一种Bayer滤波阵列彩色相机图像非均匀性校正方法
CN103033270A (zh) * 2012-12-17 2013-04-10 无锡艾立德智能科技有限公司 一种红外热像仪非均匀校正系数生成及控制方法
CN104580944A (zh) * 2013-10-10 2015-04-29 中国科学院光电研究院 对ccd图像进行相对辐射校正的方法
CN105160631A (zh) * 2015-07-02 2015-12-16 山东大学 一种求辐射校正系数的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KANG CHANG-QING: "Non uniformity correction algorithm for IR images based on midway equalization histogram", 《LASER & INFRARED》 *
W. JIN: "Infrared nonuniformity correction and radiometric calibration technology using U-shaped blackbody", 《SOCIETY OF PHOTO-OPTICAL INSTRUMENTATION ENGINEERS CONFERENCE》 *
黄本雄等: "红外图像自适应的非均匀校正算法", 《华中科技大学学报(自然科学版)》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107796518A (zh) * 2017-09-22 2018-03-13 上海卫星工程研究所 线阵扫描相机对复杂地球背景成像试验方法
CN107796518B (zh) * 2017-09-22 2019-08-09 上海卫星工程研究所 线阵扫描相机对复杂地球背景成像试验方法
CN108873924A (zh) * 2018-07-26 2018-11-23 中国人民解放军陆军工程大学 一种机载视频对地摆扫面积计算方法
CN108873924B (zh) * 2018-07-26 2021-06-29 中国人民解放军陆军工程大学 一种机载视频对地摆扫面积计算方法
CN110906906A (zh) * 2019-12-02 2020-03-24 王睿 一种增加遥感图像测绘精度方法
CN111257238A (zh) * 2020-02-05 2020-06-09 国家卫星气象中心(国家空间天气监测预警中心) 一种基于星载太阳漫射板的探测元件间相对定标方法
CN111257238B (zh) * 2020-02-05 2023-03-14 国家卫星气象中心(国家空间天气监测预警中心) 一种基于星载太阳漫射板的探测元件间相对定标方法
CN112526506A (zh) * 2020-11-17 2021-03-19 中国科学院长春光学精密机械与物理研究所 目标搜索跟踪方法及目标追踪装置
CN112526506B (zh) * 2020-11-17 2024-03-01 中国科学院长春光学精密机械与物理研究所 目标搜索跟踪方法及目标追踪装置
CN114252161A (zh) * 2021-11-29 2022-03-29 湖北久之洋信息科技有限公司 一种基于线列红外探测器的成像方法
CN114252161B (zh) * 2021-11-29 2023-08-08 湖北久之洋信息科技有限公司 一种基于线列红外探测器的成像方法

Also Published As

Publication number Publication date
CN106600646B (zh) 2019-05-24

Similar Documents

Publication Publication Date Title
CN106600646A (zh) 一种红外摆扫相机在轨图像均匀性校正方法
CN104217404B (zh) 雾霾天视频图像清晰化处理方法及其装置
CN102324098B (zh) 一种结合实验室定标和均匀景统计的相对辐射定标方法
CN103983343B (zh) 一种基于多光谱影像的卫星平台震颤检测方法及系统
CN103673995B (zh) 一种线阵推扫式相机在轨光学畸变参数标定方法
CN106871925B (zh) 一种在轨综合动态调整的遥感卫星相对辐射定标处理方法
CN104299228B (zh) 一种基于精确点位预测模型的遥感影像密集匹配方法
CN105139350A (zh) 一种无人机侦察图像地面实时重建处理系统
CN109389630B (zh) 可见光图像与红外图像特征点集确定、配准方法及装置
CN106372268B (zh) 基于热模型的实时红外仿真方法
CN104008543A (zh) 一种图像融合质量评价方法
CN102341827A (zh) 利用大气和太阳照射建模来融合图像对的系统和方法
CN104318583B (zh) 一种可见光宽带光谱图像配准方法
CN108038893A (zh) 一种基于混合相机阵列的十亿像素视频的生成方法
CN105095608A (zh) 一种卫星杂散光的测试方法
CN108830811A (zh) 一种飞行参数与相机内参相结合的航空影像实时校正方法
Zhou et al. Scene-based nonuniformity correction for airborne point target detection systems
CN106530226A (zh) 获取高分辨率高清晰工业图像的实现方法
CN105139432B (zh) 基于高斯模型的红外弱小目标图像仿真方法
Jo et al. Dense thermal 3d point cloud generation of building envelope by drone-based photogrammetry
CN103870847A (zh) 一种低照度环境下对地监控的运动目标检测方法
CN113762161A (zh) 一种障碍物智能监测方法及系统
CN108076341A (zh) 一种视频卫星成像在轨实时稳像方法及系统
CN116109540B (zh) 基于粒子群优化灰度曲线匹配的图像配准融合方法和系统
CN111595781A (zh) 一种曲面拟合的地面高光谱影像反射率校正方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Cao Shixiang

Inventor after: Li Yan

Inventor after: Jin Libing

Inventor after: Nie Yunsong

Inventor after: He Hongyan

Inventor after: Jiang Cheng

Inventor after: Zhou Nan

Inventor after: Zhang Bingxian

Inventor after: Li Fangqi

Inventor before: Cao Shixiang

Inventor before: Li Yan

Inventor before: Jin Libin

Inventor before: Nie Yunsong

Inventor before: He Hongyan

Inventor before: Jiang Cheng

Inventor before: Zhou Nan

Inventor before: Zhang Bingxian

Inventor before: Li Fangqi

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant