CN103778637B - 一种基于直方图平衡的在轨90°辐射定标处理方法 - Google Patents
一种基于直方图平衡的在轨90°辐射定标处理方法 Download PDFInfo
- Publication number
- CN103778637B CN103778637B CN201410036433.2A CN201410036433A CN103778637B CN 103778637 B CN103778637 B CN 103778637B CN 201410036433 A CN201410036433 A CN 201410036433A CN 103778637 B CN103778637 B CN 103778637B
- Authority
- CN
- China
- Prior art keywords
- calibration data
- histogram
- calibration
- radiation correction
- relative
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 48
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 238000001514 detection method Methods 0.000 claims abstract description 15
- 238000013519 translation Methods 0.000 claims abstract description 8
- 238000010276 construction Methods 0.000 claims abstract description 7
- 238000012937 correction Methods 0.000 claims description 41
- 238000000034 method Methods 0.000 claims description 21
- 238000003384 imaging method Methods 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 230000004907 flux Effects 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 5
- 238000013139 quantization Methods 0.000 claims description 4
- 238000009533 lab test Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 abstract 1
- 239000000284 extract Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明一种基于直方图平衡的在轨90°辐射定标处理方法,(1)选择原始定标数据,对原始定标数据做平移处理,提取有效定标数据;(2)根据步骤(1)中提取的有效定标数据,计算所有探测元的标准累计直方图进而得到每个探测元的累计直方图;(3)由上步得到的累计直方图结合实验室相对辐射校正模型构建相对辐射校正查找表;(4)利用上述查找表对原始图像进行相对辐射校正。采用本发明提高了side‑slither定标的适用性,实现了基于直方图平衡的side‑slither相对辐射校正。
Description
技术领域
本发明涉及一种基于直方图平衡的side-slither辐射定标处理方法。
背景技术
为获得高质量的遥感图像,航天光学遥感相机需要进行相对辐射定标工作。对于有较强敏捷特性的推扫式成像光学遥感卫星,进行在轨相对辐射定标时,可以使用将相机绕卫星偏航轴旋转一定角度(通常为90°)的方法来获取相对均匀图像的方法来进行不同像元间的非均匀性校正。这种在轨相对辐射定标方法被称为“side-slither”。
side-slither定标方法原理图如附图1所示。推扫式光学遥感相机正常成像模式下,焦面探测器线列与成像方向垂直;进行side-slither定标时,将焦面旋转90°,探测器线列方向与成像方向平行。这样,在不考虑其他影响因素的情况下,理论上探测器线列上每个像元都依次同样的地面区域成像。而side-slither定标时焦面探测线列输出得到的图像中,每个像元获得的辐射能量信息是一样的,这样就相当于满足了相对辐射定标所需要的给焦面探测线列提供均匀辐照度场的条件,在此基础上进行焦面像元非均匀性校正工作。
side-slither定标方法典型的辐射定标方法为:通过获取一幅均匀地区例如大部分由冰面覆盖的区域,绿地或南极大陆,这样就能得到相同光照条件下的探测单元的响应,然后可以将这种校正方法用于每一次测试。这种方法对于目前大多数光电线性***来说都是非常有效的,这种方法的一个主要不足是这种极地区域的云层覆盖,这种现象在成像过程中很难避免。每一个探测单元的响应必须测量有效,这需要相当长的均衡过程,通常要2至3周。
Philippe KUBIK使用基于直方图均衡的方法对Pleiades-HR的side-slither定标数据进行处理,从而进行相对辐射校正,但是这个方法有两点不足:1)实际成像***的线性响应模型与Philippe KUBIK使用的模型并不一致(Philippe KUBIK的线性模型:在低端DN<DNs时为一个线性响应模型,在中高端DN>DNs为一个线性响应模型);2)在轨side-slither定标时成像区域很难覆盖响应低端DN<DNs,此时不能保证低端的校正效果。
发明内容
本发明的技术解决问题是:针对推扫型成像遥感器,提供了一种基于直方图平衡的定标数据处理方法,实现了在轨基于直方图的side-slither相对辐射定标。
本发明的技术解决方案是:一种基于直方图平衡的side-slither辐射定标处理方法,步骤如下:
(1)选择原始定标数据,对原始定标数据做平移处理,提取有效定标数据;
(2)根据步骤(1)中提取的有效定标数据,计算所有探测元的标准累计直方图进而得到每个探测元的累计直方图;
(3)由上步得到的累计直方图结合实验室相对辐射校正模型构建相对辐射校正查找表;
(4)利用上述查找表对原始图像进行相对辐射校正。
所述步骤(1)中原始定标数据的选择原则为:图像成像行数M必须足够多,假设探测器总元数为N,图像数据的量化位数为n,那么M>N+100*2n。
所述步骤(3)中构建相对辐射校正查找表依据每个探测元的DN值确定,具体过程如下:
(3.1)对于第j个探测元,判断该探测元DN值所在的区间,当DN∈[DNminj,DNmaxj]时,转步骤(3.2);当DN∈[0,DNminj]或者DN∈[DNmaxj,2n-1]时,转步骤(3.3);[DNminj,DNmaxj]为第j个探测元有效定标数据的DN灰度范围;
(3.2)由的对应关系建立此范围内的辐射校正查找表;的对应关系如下:
其中,k=1,2,3,...,K,K=DNmax-DNmin+1,[DNmin,DNmax]为Mi的DN值范围;为Mi中DN值小于的像素个数,h[j] -1为第j个探测元的累积直方图h[j][X]的逆函数,j=1,2,3,...,N;Mi为步骤(1)提取的有效定标数据第i行的均值;
(3.3)利用实验室测试的***线性响应模型构建实验室相对辐射校正模型:
其中,为第j个探测元在辐射通量为Φr时的输出,Yr为第j个探测元在辐射通量为Φr相对校正后的输出;
(3.4)利用步骤(3.2)中的对应关系和步骤(3.3)中构建的实验室相对辐射校正模型,统计入轨后的校正模型变化量Zk:
Yk′为令时,利用步骤(3.3)构建的实验室相对辐射校正模型得到的第j个探测元在辐射通量为Φr相对校正后的输出;
(3.5)利用步骤(3.3)中构建的实验室相对辐射校正模型以及步骤(3.4)中的校正模型变化量Zk,确定Yr′之间的关系,利用该关系构建[0,DNminj]和[DNmaxj,2n-1]范围内的查找表。
本发明与现有技术相比有益效果为:
(1)本发明不但适用于线性响应***,也可以有效的校正非线性响应***,并可以保证低端和高端的校正效果。;
(2)本发明的定标不再需要针对均匀地区成像,此方法定标时只要保证定标时所有探测单元严格通过景区内的同样点即可,定标周期短;
(3)本发明可以有效校正响应低端及高端的辐射非均匀性;
(4)本发明实现方便,更新校正系数周期短,可以不定时地对遥感相机进行在轨相对辐射定标;
附图说明
图1为side-slither定标方法原理图;
图2为本发明的处理流程图;
图3为步骤一中对角线平移过程示意图,其中,图3a为原始定标图像,图3b为平移后的定标图像;
图4为视场拼接探测器线阵位置示意图;
图5为某个探测元的累计直方图曲线;
图6为***线性响应模型。
具体实施方式
下面就结合图2所示流程图对本发明做进一步介绍。
(1)选择原始定标数据,对原始定标数据做平移处理,提取有效定标数据;
首先,选择原始side-slither定标数据F,选择的原则为:图像的灰度范围覆盖足够宽,最好覆盖动态范围;图像中尽量没有饱和的景物;图像成像行数M必须足够多,假设探测器总元数为N,图像数据的量化位数为n,那么M>N+100*2n;
然后,针对原始side-slither定标数据F,根据焦面探测器的像元数及排列方式、拼接方式和搭接的像元数,得到像元的成像时间顺序,把不同像元对同一景物的成像数据重新排列,最终在同一行上的数据为不同像元对同一地物成像的数据,平移前后对比图如图3所示。
像元排列方法如下(假设每片探测器的像元数为N1,共有N2片探测器通过光学拼接排列成一长线阵,N3为片间搭接像元数,N4为采用视场拼接的偶数线阵相对于奇数线阵在沿飞行方向上延时成像的像元数):
a)针对原始定标数据进行如下的平移处理:
G(i,N1×(l-1)+j)=F(i+j-1,N1×(l-1)+j)
l=1,2,3,...,N2,i=1,2,3,...,M-N1,j=1,2,3,...,N1
b)根据片间探测器的拼接方式及搭接的像元数,对不同片间的定标数据做如下平移处理:
当采用光学拼接时:
G′(i,N1×(l-1)+j)=G(N3×(l-1)+i,N1×(l-1)+j)
l=1,2,3,...,N2,i=1,2,3,...,N5,N5=M-N1-(N2-1)×N3,j=1,2,3,...,N1
当采用视场拼接时(如图4所示):
l为奇数时:G′(i,N1×(l-1)+j)=G(N3×(l-1)+i,N1×(l-1)+j)
i=1,2,3,...,N5,N5=M-N1-(N2-1)×N3,j=1,2,3,...,N1
l为偶数时:G′(i,N1×(l-1)+j)=G(N3×(l-1)+N4+i,N1×(l-1)+j)
i=1,2,3,...,N5,N5=M-N1-(N2-1)×N3-N4,j=1,2,3,...,N1
G′为所有像元排列后的图像;
c)剔除平移后的图像数据中的饱和数据,并最终得到有效的定标数据G″:
计算每一行的均值:
当Mi>2n-1时此行数据为饱和数据,假设总共有N6行数据为饱和数据,此时剔除后的有效定标数据G″行数为Nb=N5-N6,列数N=N1·N2。
(2)根据步骤(1)中提取的有效定标数据,计算所有探测元的标准累计直方图进而得到每个探测元的累计直方图;
(2.1)针对上一步得到的图像G″,计算每行数据的均值Mi,i=1,2,3,...,Nb;
(2.2)构建标准累计直方图q1,q2,...,qK∈[0,1],构建方法如下:
其中,k=1,2,3,...,K,K=DNmax-DNmin+1,[DNmin,DNmax]为Mi的DN值范围;为Mi中DN值小于的像素个数。DNminj
(2.3)对于第j个探测元计算归一化的累计直方图如图5所示:
其中,Nx为G″中第j个探测元所有DN值小于X的像素个数,[DNminj,DNmaxj]为G″中第j个探测元的DN值范围,j=1,2,3,...,N;X为变量,取值范围[DNminj,DNmaxj]。
(2.4)构建序列使得即第j元中小于的像元总数占所有像元总素的比例为qk:
其中,h[j] -1为h[j][X]的逆函数,j=1,2,3,...,N。
(3)由上步得到的累计直方图结合实验室辐射校正模型构建相对辐射校正查找表;
对于第j个探测元:
a.当DN值X∈[DNminj,DNmaxj]时,由的对应关系建立此探测元的辐射校正查找表(即利用公式1、2、4可以得到与之间的一一对应关系):
其中k=1,2,3,...,K,j=1,2,3,...,N。
b.当DN值X∈[0,DNminj]或当DN值X∈[DNmaxj,2n-1]时,构建[0,DNminj]和[DNmaxj,2n-1]范围内的辐射校正查找表:
i.由实验室的***线性响应模型(如图6所示),构建实验室的相对辐射校正模型构建方法如下:
其中,为第j个探测元在辐射通量为Φr时的输出,Φr为辐射通量,Yr第j个探测元在辐射通量为Φr相对校正后的输出。
ii.当DN值X∈[DNminj,DNmaxj],由的对应关系和实验室相对辐射校正模型统计入轨后的校正模型变化量Zk:
令:由(即公式5),得到此时第j个探测元在辐射通量为Φr相对校正后的输出记为Yk′,进而得到(Yk′二者之间的一一对应关系);
其中,k=1,2,3,...,K,j=1,2,3,...,N。
iii.构建[0,DNminj]和[DNmaxj,2n-1]范围内的查找表:
其中,r=0,1,2,...,DNminj,或者r=DNmaxj,DNmaxj+1,DNmaxj+2,...,2n-1,
Yr′=Yr·Z,
(4)利用上述查找表对原始图像进行相对辐射校正。
对于第j个探测元:
当原始图像的DN值X∈[DNminj,DNmaxj],使用辐射校正查找表找到X对应的DN值Y,把Y值赋值给校正后的对应位置的DN值进行相对辐射校正;
当原始图像的DN值X∈[0,DNminj]或X∈[DNmaxj,2n-1]的时候,使用辐射校正查找表找到X对应的DN值Y′,把Y′赋值给校正后的对应位置的DN值进行相对辐射校正。
本发明未详细说明部分属于本领域技术人员公知常识。
Claims (2)
1.一种基于直方图平衡的在轨90°辐射定标处理方法,其特征在于步骤如下:
(1)选择原始定标数据,对原始定标数据做平移处理,提取有效定标数据;
(2)根据步骤(1)中提取的有效定标数据,计算所有探测元的标准累计直方图进而得到每个探测元的累计直方图;
(3)由上步得到的累计直方图结合实验室相对辐射校正模型构建相对辐射校正查找表;所述构建相对辐射校正查找表依据每个探测元的DN值确定,具体过程如下:
(3.1)对于第j个探测元,判断该探测元DN值所在的区间,当DN∈[DNminj,DNmaxj]时,转步骤(3.2);当DN∈[0,DNminj]或者DN∈[DNmaxj,2n-1]时,转步骤(3.3);[DNminj,DNmaxj]为第j个探测元有效定标数据的DN灰度范围;所述的n为图像数据的量化位数;
(3.2)由序列的对应关系建立此范围内的辐射校正查找表;的对应关系如下:
其中,k=1,2,3,...,K,K=DNmax-DNmin+1,[DNmin,DNmax]为Mi的DN值范围;为Mi中DN值小于的像素个数,h[j] -1为第j个探测元的累积直方图h[j][X]的逆函数,j=1,2,3,...,N;Mi为步骤(1)提取的有效定标数据第i行的均值,qk为第j元中小于的像元总数占所有像元总素的比例;Nb为有效定标数据G″的行数;
(3.3)利用实验室测试的***线性响应模型构建实验室相对辐射校正模型:
其中,为第j个探测元在辐射通量为Φr时的输出,Yr为第j个探测元在辐射通量为Φr相对校正后的输出;
(3.4)利用步骤(3.2)中的对应关系和步骤(3.3)中构建的实验室相对辐射校正模型,统计入轨后的校正模型变化量Zk:
Yk′为令时,利用步骤(3.3)构建的实验室相对辐射校正模型得到的第j个探测元在辐射通量为Φr相对校正后的输出;
(3.5)利用步骤(3.3)中构建的实验室相对辐射校正模型以及步骤(3.4)中的校正模型变化量Zk,确定Yr′之间的关系,利用该关系构建[0,DNminj]和[DNmaxj,2n-1]范围内的查找表;所述的Yr′=Yr·Z,
(4)利用上述查找表对原始图像进行相对辐射校正。
2.如权利要求1所述的一种基于直方图平衡的在轨90°辐射定标处理方法,其特征在于:所述步骤(1)中原始定标数据的选择原则为:图像成像行数M必须足够多,假设探测器总元数为N,图像数据的量化位数为n,那么M>N+100*2n。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036433.2A CN103778637B (zh) | 2014-01-24 | 2014-01-24 | 一种基于直方图平衡的在轨90°辐射定标处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036433.2A CN103778637B (zh) | 2014-01-24 | 2014-01-24 | 一种基于直方图平衡的在轨90°辐射定标处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103778637A CN103778637A (zh) | 2014-05-07 |
CN103778637B true CN103778637B (zh) | 2017-05-31 |
Family
ID=50570833
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410036433.2A Active CN103778637B (zh) | 2014-01-24 | 2014-01-24 | 一种基于直方图平衡的在轨90°辐射定标处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103778637B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106851141B (zh) * | 2016-12-14 | 2019-06-28 | 中国资源卫星应用中心 | 一种遥感图像的非均匀性校正方法 |
CN108776955B (zh) * | 2018-04-16 | 2020-08-18 | 国家卫星气象中心 | 一种遥感图像实时校正方法及校正装置 |
CN109712089A (zh) * | 2018-12-14 | 2019-05-03 | 航天恒星科技有限公司 | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 |
CN110009688A (zh) * | 2019-03-19 | 2019-07-12 | 北京市遥感信息研究所 | 一种红外遥感影像相对辐射定标方法、***及遥感平台 |
CN113469899B (zh) * | 2021-06-04 | 2023-12-29 | 中国资源卫星应用中心 | 一种基于辐射能量重构的光学遥感卫星相对辐射校正方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799871A (zh) * | 2010-01-18 | 2010-08-11 | 浙江林学院 | 基于回归参数变换的遥感图像专题信息提取方法 |
CN103226819A (zh) * | 2013-03-29 | 2013-07-31 | 武汉大学 | 一种基于分段统计的相对辐射校正方法 |
-
2014
- 2014-01-24 CN CN201410036433.2A patent/CN103778637B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799871A (zh) * | 2010-01-18 | 2010-08-11 | 浙江林学院 | 基于回归参数变换的遥感图像专题信息提取方法 |
CN103226819A (zh) * | 2013-03-29 | 2013-07-31 | 武汉大学 | 一种基于分段统计的相对辐射校正方法 |
Non-Patent Citations (4)
Title |
---|
CBERS-01/02卫星CCD图像相对辐射校正研究;郭建宁 等;《中国科学E辑:信息科学》;20051230(第2005年卷第S1期);第15-16页、第21页 * |
CBERS-02B卫星TDICCD相机的相对辐射定标方法及结果;何红艳 等;《航天返回与遥感》;20100831;第31卷(第4期);第38-44页 * |
Destriping multiple sensor imagery by improved histogram matching;M Wegener et al.;《International Journal of Remote Sensing》;19901231;第11卷(第5期);第859-875页 * |
基于探元直方图匹配的CBERS-01星CCD数据相对辐射校正方法;潘志强 等;《武汉大学学报(信息科学版)》;20051030;第30卷(第10期);第925-927页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103778637A (zh) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103778637B (zh) | 一种基于直方图平衡的在轨90°辐射定标处理方法 | |
Murthy et al. | SkySat-1: very high-resolution imagery from a small satellite | |
CN103528694B (zh) | 一种用红外热像仪测量目标物体的温度的方法 | |
CN109870239B (zh) | 非制冷红外焦平面探测器自适应定标方法 | |
CN103487058B (zh) | 一种提高aps星敏感器动态性能的方法 | |
Chary et al. | The nature of faint 24 micron sources seen in spitzer space telescope observations of elais-n1 | |
CN107255521B (zh) | 一种红外图像非均匀性校正方法及*** | |
CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及*** | |
CN109712089A (zh) | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 | |
CN104298887B (zh) | 一种多片线阵ccd相机的相对辐射定标方法 | |
CN108645401B (zh) | 基于姿态关联图像叠加的全天时星敏感器星点提取方法 | |
CN104977024B (zh) | 一种日盲紫外遥感相机绝对辐射定标系数在轨修正方法 | |
KR102071879B1 (ko) | 하늘 영상을 이용한 미세먼지 측정 시스템 | |
CN110411585B (zh) | 一种高精度红外辐射测量方法 | |
CN102938137A (zh) | 基于图像序列分析的线扫描图像动态非均匀性校正方法 | |
CN106878700B (zh) | 广角镜头校正***及其方法 | |
CN102752504B (zh) | 一种宽视场线阵ccd相机的相对辐射校正方法 | |
Plazas et al. | Transverse electric fields' effects in the Dark Energy Camera CCDs | |
CN107271046A (zh) | 一种基于局部统计相关性的非均匀性校正方法及装置 | |
Mammoli et al. | An experimental method to merge far-field images from multiple longwave infrared sensors for short-term solar forecasting | |
CN105092043B (zh) | 一种基于场景的变积分时间的非均匀性校正方法 | |
CN103312983B (zh) | 一种太赫兹成像仪镜头补偿方法 | |
Savoy et al. | Geo-referencing and stereo calibration of ground-based whole sky imagers using the sun trajectory | |
CN102410844B (zh) | 一种高动态星敏感器图像非均匀校正方法及装置 | |
CN112710397B (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 |