CN104156563A - 基于多重欠采样实现光学遥感器像能量超频采样的方法 - Google Patents

基于多重欠采样实现光学遥感器像能量超频采样的方法 Download PDF

Info

Publication number
CN104156563A
CN104156563A CN201410339156.2A CN201410339156A CN104156563A CN 104156563 A CN104156563 A CN 104156563A CN 201410339156 A CN201410339156 A CN 201410339156A CN 104156563 A CN104156563 A CN 104156563A
Authority
CN
China
Prior art keywords
centerdot
sampling
sample frequency
energy
under
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
CN201410339156.2A
Other languages
English (en)
Other versions
CN104156563B (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 CN201410339156.2A priority Critical patent/CN104156563B/zh
Publication of CN104156563A publication Critical patent/CN104156563A/zh
Application granted granted Critical
Publication of CN104156563B publication Critical patent/CN104156563B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Transforming Light Signals Into Electric Signals (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明基于多重欠采样实现光学遥感器像能量超频采样的方法,包括下列步骤:确定光学遥感器焦面像元尺寸以及在该像元尺寸下的采样频率f1和所需采样频率f2;计算成像信号的移动步数k和移动步长dx;利用采样频率f1对光学遥感器像能量进行采样得到采样数据矩阵SD;将得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值。本发明基于欠采样焦面对真实像能量分布的进行不同位置的采样,从而实现欠采样焦面的超分辨率采样。在提高采样频率的同时,采样数据都是真实像能量分布上的值,本发明有较高的实用价值,在实现光学遥感器光学传递函数的计算有着重要的应用。

Description

基于多重欠采样实现光学遥感器像能量超频采样的方法
技术领域
本发明涉及一种基于多重欠采样实现光学遥感器像能量超频采样的方法,属于光学测试技术领域。
背景技术
高分辨率、高光学传递函数值是光学遥感器长久发展以来一直追求的目标。由于其光学镜头的衍射极限的限制,焦面器件的采样间隔都不会小于光学镜头的极限分辨率,采样频率受到限制,在整机状态下焦面器件的采样频率都会低于重建实际像能量分布所需的频率,也即低于像能量分布的乃奎斯特频率,因此无法全面的评价遥感器整机的光学性能。
现有技术一是利用高倍率的显微物镜将像放大,再用焦面器件接收放大后的像能量,从而实现高频率的采样,但是这种方法在光学遥感器整机测试中因其整机为封闭结构,无法加入显微物镜。二是使用图像处理方法实现高频采样,主要有曲面拟合和模板运算,但这种方法得到的高频采样值并非真实像能量分布值,精度较差。
发明内容
本发明的技术解决问题是:针对现有技术的不足,提供了一种基于多重欠采样实现光学遥感器像能量超频采样的方法,本发明采用移相采样的方式实现光学遥感器的像能量分布的超频采样,且采样值为实际像能量分布值,为光学遥感器整机光学性能评价和实现超分辨率成像奠定了基础。
本发明的技术解决方案是:
基于多重欠采样实现光学遥感器像能量超频采样的方法,包括下列步骤:
(1)基于欠采样焦面器件的响应正比于落于其上的像能量在空间上的积分值的特性,再采用移相采样的方式实现像能量分布的超频采样。确定光学遥感器焦面像元尺寸a以及在该像元尺寸下的欠采样频率f1和所需采样频率f2,f2>f1(f1为欠采样低频率,f2为超频采样频率,通过利用f1进行多次采样,求的在f2采样下的能量分布);
(2)计算成像信号的移动步数k和移动步长dx:
k=f2/f1,且向上取整
dx=a/k
相当于将每个像元等分成k个更小的宽度为dx的单元,将像相对于焦面以dx为步长移动,共移动k-1次。
(3)利用采样频率f1对光学遥感器像能量进行采样得到采样数据矩阵SD;具体实现方法如下:
(3a)利用采样频率f1对像能量进行采样,得到采样数据SD1,记为:
SD1=[P11 P12...P1n]
其中
P 11 = x 1 + x 2 + . . + x k P 12 = x k + 1 + x k + 2 + . . + x 2 k · · P 1 n = x kn - k + 1 + x kn - k + 2 + . . + x kn
P11 P12...P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2...xkn...]为像能量在采样频率f2下的采样值;
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22...P2n]
其中
P 21 = x 2 + x 3 + . . + x k + 1 P 22 = x k + 2 + x k + 3 + . . + x 2 k + 1 · · P 2 n = x kn - k + 2 + x kn - k + 3 + . . + x kn + 1
P21 P22...P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1 Pk2...Pkn]
其中
P k 1 = x k + x k + 1 + . . + x 2 k - 1 P k 2 = x 2 k + x 2 k + 1 + . . + x 3 k - 1 · · P kn = x kn + x kn + 1 + . . + x kn + k - 1
Pk1 Pk2...Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2...SDk组成采样数据矩阵SD:
SD = SD 1 · · SD r · · SD k = P 11 . . P 1 c . . P 1 n · · · · · · P r 1 . . P rc . . P rn · · · · · · P k 1 . . P kc . . P kn
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2...xkn...],具体计算方法如下:
x kn + 1 x kn + 2 · · x kn + k - 1 = P 1 n / k P 1 n / k · · P 1 n / k , 得到
x 1 x 2 · · x k - 1 = P 1 n / k - ( Σ i = 1 n P 2 i - Σ i = 1 n P 1 i ) P 1 n / k - ( Σ i = 1 n P 3 i - Σ i = 1 n P 2 i ) · · P 1 n / k - ( Σ i = 1 n P ki - Σ i = 1 n P ( k - 1 ) i )
x 1 x 2 · · x k - 1 代入下式得到在采样频率f2下的采样值X=[x1,x2...xkn...]:
X s = P ij - Σ h = s - k s - 1 x h
其中Pij为采样矩阵SD中的值,s=k:1:kn,的余数,
本发明与现有技术相比的有益效果是:
(1)对于光学遥感器的整机测试,焦面器件对测试信号的采样会形成欠采样,也就是采样频率达不到测试计算的要求;对于整机光学传递函数测试来说,需要使用焦面器件对星点像进行采样,但只能形成欠采样,本发明能够实现焦面器件的超频采样,从而实现完整的光学传递函数曲线的测试,进而对整机状态下的光学遥感器的成像性能进行评价,本发明适用性强,测试成本低,效果高。
(2)本发明能够实现相机的超分辨率成像,通过移动成像位置,利用多幅欠采样图像,计算得到更高分辨率的图像;同时本发明在光学测试和相机超分辨率成像上有着重要的作用,并且易于实现,能够大大降低光学遥感器的制造成本。
附图说明
图1为本发明流程图;
图2为本发明的原理示意图。
具体实施方式
下面结合附图对本发明的工作原理和工作过程作进一步解释和说明。
本发明使用欠采样焦面器件,实现了超频采样,并且超频采样值为实际像能量分布值。在实现光学遥感器光学传递函数的计算有着重要的应用。光学遥感器在整机状态下起焦面器件的采样频率都会低于重建实际像能量分布所需的频率,也即低于像能量分布的乃奎斯特频率,因此无法全面的评价遥感器整机的光学性能,本发明采用移相采样的方式实现光学遥感器的像能量分布的超频采样,且采样值为实际像能量分布值,为光学遥感器整机光学性能评价和实现超分辨率成像奠定了基础。
如图1所示,基于多重欠采样实现光学遥感器像能量超频采样的方法,包括下列步骤:
(1)基于欠采样焦面器件的响应正比于落于其上的像能量在空间上的积分值的特性,再采用移相采样的方式实现像能量分布的超频采样;确定光学遥感器焦面像元尺寸a以及在该像元尺寸下的欠采样频率f1和所需采样频率f2,f2>f1(f1为欠采样低频率,f2为超频采样频率,通过利用f1进行多次采样,求的在f2采样下的能量分布);
(2)计算成像信号的移动步数k和移动步长dx:
k=f2/f1,且向上取整
dx=a/k
相当于将每个像元等分成k个更小的宽度为dx的单元,将像相对于焦面以dx为步长移动,共移动k-1次。
(3)利用采样频率f1对光学遥感器像能量进行采样得到采样数据矩阵SD;具体实现方法如下:
(3a)利用采样频率f1对像能量进行采样,得到采样数据SD1,记为:
SD1=[P11 P12...P1n]
其中
P 11 = x 1 + x 2 + . . + x k P 12 = x k + 1 + x k + 2 + . . + x 2 k · · P 1 n = x kn - k + 1 + x kn - k + 2 + . . + x kn
P11 P12...P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2...xkn...]为像能量在采样频率f2下的采样值;
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22...P2n]
其中
P 21 = x 2 + x 3 + . . + x k + 1 P 22 = x k + 2 + x k + 3 + . . + x 2 k + 1 · · P 2 n = x kn - k + 2 + x kn - k + 3 + . . + x kn + 1
P21 P22...P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1 Pk2...Pkn]
其中
P k 1 = x k + x k + 1 + . . + x 2 k - 1 P k 2 = x 2 k + x 2 k + 1 + . . + x 3 k - 1 · · P kn = x kn + x kn + 1 + . . + x kn + k - 1
Pk1 Pk2...Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2...SDk组成采样数据矩阵SD:
SD = SD 1 · · SD r · · SD k = P 11 . . P 1 c . . P 1 n · · · · · · P r 1 . . P rc . . P rn · · · · · · P k 1 . . P kc . . P kn
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2...xkn...];具体计算方法如下:
首先将采样数据矩阵SD相邻行相减,得到矩阵TM1
TM 1 = x k + 1 - x 1 x 2 k + 1 - x k + 1 . . x kn + 1 - x kn - k + 1 x k + 2 - x 2 x 2 k + 2 - x k + 2 . . x kn + 2 - x kn - k + 2 · · · · · · · · x 2 k - 1 - x k - 1 x 3 k - 1 - x 2 k - 1 . . x kn + k - 1 - x kn - 1 = P 21 - P 11 P 22 - P 12 . . P 2 n - P 1 n P 31 - P 21 P 32 - P 22 . . P 3 n - P 2 n · · · · · · · · P k 1 - P ( k - 1 ) 1 P k 2 - P ( k - 1 ) 2 . . P kn - P ( k - 1 n )
将TM1的相邻列两两相加,得到矩阵TM2
TM 2 = x kn + 1 - x 1 x kn + 2 - x 2 · · x kn + k - 1 - x k - 1 = Σ i = 1 n P 2 i - Σ i = 1 n P 1 i Σ i = 1 n P 3 i - Σ i = 1 n P 2 i · · Σ i = 1 n P ki - Σ i = 1 n P ( k - 1 ) i
一般数据采样区域边缘成像能量已非常微弱,可认为成像能量在边缘像元处分布均匀,所以可以令 x kn + 1 x kn + 2 · · x kn + k - 1 = P 1 n / k P 1 n / k · · P 1 n / k , 得到
x 1 x 2 · · x k - 1 = P 1 n / k - ( Σ i = 1 n P 2 i - Σ i = 1 n P 1 i ) P 1 n / k - ( Σ i = 1 n P 3 i - Σ i = 1 n P 2 i ) · · P 1 n / k - ( Σ i = 1 n P ki - Σ i = 1 n P ( k - 1 ) i )
x 1 x 2 · · x k - 1 代入下式得到在采样频率f2下的采样值X=[x1,x2...xkn...]:
X s = P ij - Σ h = s - k s - 1 x h
其中Pij为采样矩阵SD中的值,s=k:1:kn,的余数,例如根据P11=x1+x2+..+xk可以求出xk,再利用P21=x2+x3+..+xk+1可以求出xk+1,依次类推能够全部求出。
下面以一个具体实施例对本发明的工作过程和工作原理作进一步解释:
如图2所示,焦面像元尺寸为a,所需采样频率为f2
(1)确定在该焦面像元尺寸的采样频率使
(2)计算成像信号的移动步数k和移动步长dx:
k=f2/f1=2
dx=a/k=a/2
(3)利用采样频率f1对光学遥感器像能量进行采样得到采样数据矩阵SD:
(3a)利用采样频率f1对像能量进行采样,得到采样数据SD1,记为:
SD1=[P11 P12 P13]
其中:
P11=x1+x2
P12=x3+x4
P13=x5+x6
(3b)将成像信号相对于光学遥感器的焦面移动再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22 P23]
其中:
P21=x2+x3
P22=x4+x5
P23=x6+x7
得到欠采样数据矩阵SD:
SD = SD 1 SD 2 = P 11 P 12 P 13 P 21 P 22 P 23
(4)对SD进行数据处理,得到采样频率f2下的采样值X=[x1 x2 x3 x4 x5 x6]:
x 7 = P 13 2 , 得到
x 1 = P 13 2 - ( Σ 1 3 P 2 i - Σ 1 3 P 1 i )
由x1计算x2, s = 2 , i = s - ( k - 1 ) k = 2 - ( 2 - 1 ) 2 的余数, i = 1 , j = s - ( k - 1 ) - i k + 1 = 1 ,
因此可得:
x2=P11-x1
同理可得:
x3=P21-x2
x4=P12-x3
x5=P22-x4
x6=P13-x5
由此就得到了在采样频率f2下的采样值X。
本发明已经经过专家评定,应用在航天某探测卫星型号上,效果良好。
本发明未公开的部分都是公知常识。

Claims (1)

1.基于多重欠采样实现光学遥感器像能量超频采样的方法,其特征在于包括下列步骤:
(1)确定光学遥感器焦面像元尺寸a以及在该像元尺寸下的采样频率f1和所需采样频率f2,f2>f1
(2)计算成像信号的移动步数k和移动步长dx:
k=f2/f1,且向上取整
dx=a/k
(3)利用采样频率f1对光学遥感器像能量进行采样得到采样数据矩阵SD;具体实现方法如下:
(3a)利用采样频率f1对像能量进行采样,得到采样数据SD1,记为:
SD1=[P11 P12...P1n]
其中
P 11 = x 1 + x 2 + . . + x k P 12 = x k + 1 + x k + 2 + . . + x 2 k · · P 1 n = x kn - k + 1 + x kn - k + 2 + . . + x kn
P11 P12...P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2...xkn...]为像能量在采样频率f2下的采样值。
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22...P2n]
其中
P 21 = x 2 + x 3 + . . + x k + 1 P 22 = x k + 2 + x k + 3 + . . + x 2 k + 1 · · P 2 n = x kn - k + 2 + x kn - k + 3 + . . + x kn + 1
P21 P22...P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1 Pk2...Pkn]
其中
P k 1 = x k + x k + 1 + . . + x 2 k - 1 P k 2 = x 2 k + x 2 k + 1 + . . + x 3 k - 1 · · P kn = x kn + x kn + 1 + . . + x kn + k - 1
Pk1 Pk2...Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2...SDk组成采样数据矩阵SD:
SD = SD 1 · · SD r · · SD k = P 11 . . P 1 c . . P 1 n · · · · · · P r 1 . . P rc . . P rn · · · · · · P k 1 . . P kc . . P kn
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2...xkn...],具体计算方法如下:
x kn + 1 x kn + 2 · · x kn + k - 1 = P 1 n / k P 1 n / k · · P 1 n / k , 得到
x 1 x 2 · · x k - 1 = P 1 n / k - ( Σ i = 1 n P 2 i - Σ i = 1 n P 1 i ) P 1 n / k - ( Σ i = 1 n P 3 i - Σ i = 1 n P 2 i ) · · P 1 n / k - ( Σ i = 1 n P ki - Σ i = 1 n P ( k - 1 ) i )
x 1 x 2 · · x k - 1 代入下式得到在采样频率f2下的采样值X=[x1,x2...xkn...]:
X s = P ij - Σ h = s - k s - 1 x h
其中Pij为采样矩阵SD中的值,s=k:1:kn,的余数,
CN201410339156.2A 2014-07-16 2014-07-16 基于多重欠采样实现光学遥感器像能量超频采样的方法 Active CN104156563B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410339156.2A CN104156563B (zh) 2014-07-16 2014-07-16 基于多重欠采样实现光学遥感器像能量超频采样的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410339156.2A CN104156563B (zh) 2014-07-16 2014-07-16 基于多重欠采样实现光学遥感器像能量超频采样的方法

Publications (2)

Publication Number Publication Date
CN104156563A true CN104156563A (zh) 2014-11-19
CN104156563B CN104156563B (zh) 2017-04-19

Family

ID=51882061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410339156.2A Active CN104156563B (zh) 2014-07-16 2014-07-16 基于多重欠采样实现光学遥感器像能量超频采样的方法

Country Status (1)

Country Link
CN (1) CN104156563B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02145913A (ja) * 1988-11-26 1990-06-05 Mitsubishi Electric Corp 光学センサ装置
CN101140468A (zh) * 2007-10-11 2008-03-12 中国科学院长春光学精密机械与物理研究所 航天光学遥感器调焦控制电路的仿真测试方法
EP2060933A2 (en) * 2007-11-13 2009-05-20 Optoplan AS Non-uniform Sampling to Extend Dynamic Range of Interferometric Sensors
CN102034229A (zh) * 2010-11-03 2011-04-27 中国科学院长春光学精密机械与物理研究所 高分辨多光谱空间光学遥感器的实时图像融合方法
CN102572499A (zh) * 2010-12-07 2012-07-11 江南大学 基于小波变换多分辨率预测的无参考图像质量评价方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02145913A (ja) * 1988-11-26 1990-06-05 Mitsubishi Electric Corp 光学センサ装置
CN101140468A (zh) * 2007-10-11 2008-03-12 中国科学院长春光学精密机械与物理研究所 航天光学遥感器调焦控制电路的仿真测试方法
EP2060933A2 (en) * 2007-11-13 2009-05-20 Optoplan AS Non-uniform Sampling to Extend Dynamic Range of Interferometric Sensors
CN102034229A (zh) * 2010-11-03 2011-04-27 中国科学院长春光学精密机械与物理研究所 高分辨多光谱空间光学遥感器的实时图像融合方法
CN102572499A (zh) * 2010-12-07 2012-07-11 江南大学 基于小波变换多分辨率预测的无参考图像质量评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐杰: "提高星载成像传感器空间分辨率的方法研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 *
李勃: "随机振动对TDICCD相机像质影响分析", 《中国优秀硕士学位论文全文数据库(信息科技辑)》 *

Also Published As

Publication number Publication date
CN104156563B (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
CN102221341B (zh) 基于随机并行梯度下降优化技术的快速数字图像相关测量方法
CN104048744A (zh) 一种非接触式的基于影像的实时在线振动测量方法
CN104535998B (zh) 旋转镜像综合孔径辐射计及测量方法
CN104268876A (zh) 基于分块的摄像机标定方法
CN103063300B (zh) 一种实现全偏振成像的微偏振调制阵列
CN104089697A (zh) 一种基于线程池并发技术的实时在线的视觉振动测量方法
Xia et al. Global calibration of multi-cameras with non-overlapping fields of view based on photogrammetry and reconfigurable target
CN103428500A (zh) 一种超分辨率大视场红外成像方法
CN104391295A (zh) 一种图像熵最优的压缩传感sar稀疏自聚焦成像方法
CN103791899B (zh) 一种基于成像传感器视差的卫星姿态颤振探测方法
CN103743485A (zh) 用于同步探测地物光与天空光的全偏振光谱成像系统
CN105869141A (zh) 图像配准方法及装置
CN103900504A (zh) 纳米尺度下的实时三维视觉信息反馈方法
CN103439001A (zh) 一种非均匀矢量偏振光的测量与评价方法及装置
CN109752844A (zh) 一种基于随机光强涨落的成像方法和系统
CN104504691A (zh) 基于低秩纹理的摄像机位置和姿态测量方法
CN103424087A (zh) 一种大尺度钢板三维测量拼接系统及方法
CN103399030B (zh) 一种实现空管3/a模式下应答信号相位差检测的系统和方法
CN103776381A (zh) 一种mems微结构平面位移测量方法
CN104391286A (zh) 一种逆合成孔径雷达方位定标方法
CN103679673B (zh) 一种宽视场线阵ccd影像几何畸变模拟方法
CN104111080B (zh) 一种凝视卫星面阵ccd相机的mtf在轨检测方法
CN107152969B (zh) 一种离岸式傅里叶成像光谱仪数据处理方法
CN103955956B (zh) 一种面向压缩感知的图像联合重构方法
CN103913114A (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