CN104156563B - 基于多重欠采样实现光学遥感器像能量超频采样的方法 - Google Patents
基于多重欠采样实现光学遥感器像能量超频采样的方法 Download PDFInfo
- Publication number
- CN104156563B CN104156563B CN201410339156.2A CN201410339156A CN104156563B CN 104156563 B CN104156563 B CN 104156563B CN 201410339156 A CN201410339156 A CN 201410339156A CN 104156563 B CN104156563 B CN 104156563B
- Authority
- CN
- China
- Prior art keywords
- sampling
- sample frequency
- energy
- under
- sampled
- 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
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]
其中
P11 P12...P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2...xkn...]为像能量在采样频率f2下的采样值;
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22...P2n]
其中
P21 P22...P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1 Pk2...Pkn]
其中
Pk1 Pk2...Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2...SDk组成采样数据矩阵SD:
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2...xkn...],具体计算方法如下:
令得到
将代入下式得到在采样频率f2下的采样值X=[x1,x2...xkn...]:
其中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]
其中
P11 P12...P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2...xkn...]为像能量在采样频率f2下的采样值;
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21 P22...P2n]
其中
P21 P22...P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1 Pk2...Pkn]
其中
Pk1 Pk2...Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2...SDk组成采样数据矩阵SD:
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2...xkn...];具体计算方法如下:
首先将采样数据矩阵SD相邻行相减,得到矩阵TM1:
将TM1的相邻列两两相加,得到矩阵TM2:
一般数据采样区域边缘成像能量已非常微弱,可认为成像能量在边缘像元处分布均匀,所以可以令得到
将代入下式得到在采样频率f2下的采样值X=[x1,x2...xkn...]:
其中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:
(4)对SD进行数据处理,得到采样频率f2下的采样值X=[x1 x2 x3 x4 x5 x6]:
令得到
由x1计算x2,的余数,
因此可得:
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=[P11P12…P1n]
其中
P11=x1+x2+…+xk
P12=xk+1+xk+2+…+x2k
.
.
.
P1n=xkn-k+1+xkn-k+2+…+xkn
P11,P12,…,P1n为像能量在采样频率f1下的采样值,n为采样像元数,取正整数,X=[x1,x2…xkn…]为像能量在采样频率f2下的采样值;
(3b)将成像信号相对于光学遥感器的焦面移动dx,再利用采样频率f1对像能量进行采样得到采样数据SD2,记为:
SD2=[P21P22…P2n]
其中
P21=x2+x3+…+xk+1
P22=xk+2+xk+3+…+x2k+1
.
.
.
P2n=xkn-k+2+xkn-k+3+…+xkn+1
P21,P22,…,P2n为像能量在采样频率f1下的采样值;
(3c)以此类推,将成像信号相对于光学遥感器的焦面第k-1次移动dx,再利用采样频率f1对像能量进行采样得到采样数据SDk,记为:
SDk=[Pk1Pk2…Pkn]
其中
Pk1=xk+xk+1+…+x2k-1
Pk2=x2k+x2k+1+…+x3k-1
.
.
.
Pkn=xkn+xkn+1+…+xkn+k-1
Pk1,Pk2,…,Pkn为像能量在采样频率f1下的采样值;
(3d)将步骤(3a)-(3c)得到采样数据SD1,SD2…SDk组成采样数据矩阵SD:
其中,r为小于k的正整数,c为小于n的正整数;
(4)将步骤(3)中得到的采样数据矩阵SD进行数据处理,计算得到采样频率f2下的采样值X=[x1,x2…xkn…],具体计算方法如下:
令得到
将代入下式得到在采样频率f2下的采样值X=[x1,x2…xkn…]:
Xs为采样频率f2下序号从k到kn的采样值x;
其中Pij为采样矩阵SD中的值,k<s≤kn,s为正整数,的余数,
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 CN104156563A (zh) | 2014-11-19 |
CN104156563B true 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 (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 江南大学 | 基于小波变换多分辨率预测的无参考图像质量评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02145913A (ja) * | 1988-11-26 | 1990-06-05 | Mitsubishi Electric Corp | 光学センサ装置 |
-
2014
- 2014-07-16 CN CN201410339156.2A patent/CN104156563B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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)
Title |
---|
提高星载成像传感器空间分辨率的方法研究;徐杰;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20090715(第07期);第C031-66页 * |
随机振动对TDICCD相机像质影响分析;李勃;《中国优秀硕士学位论文全文数据库(信息科技辑)》;20140315(第03期);第I140-705页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104156563A (zh) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10048226B2 (en) | Imaging method and apparatus based on magnetic flux leakage testing | |
CN102221341B (zh) | 基于随机并行梯度下降优化技术的快速数字图像相关测量方法 | |
CN102607526B (zh) | 双介质下基于双目视觉的目标姿态测量方法 | |
US10032273B2 (en) | Machine vision system calibration using inaccurate calibration targets | |
CN103913293B (zh) | 编码结构光系统光学传递模型测量系统 | |
CN103971353A (zh) | 采用激光辅助大型锻件测量图像数据的拼接方法 | |
JP2015072713A (ja) | マルチスペクトルデータの幾何学的リファレンシング | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
CN107367329A (zh) | 一种图像、光谱、偏振态一体化获取装置及探测方法 | |
CN103234454B (zh) | 一种影像测量仪的自标定方法 | |
CN103559708A (zh) | 基于方靶模型的工业定焦相机参数标定装置 | |
CN102509304A (zh) | 基于智能优化的摄像机标定方法 | |
CN103743485A (zh) | 用于同步探测地物光与天空光的全偏振光谱成像系统 | |
Zhao et al. | Displacement monitoring for slope stability evaluation based on binocular vision systems | |
CN103115631A (zh) | 遥感相机成像参数误差校正系统及方法 | |
CN103900504A (zh) | 纳米尺度下的实时三维视觉信息反馈方法 | |
CN103424087A (zh) | 一种大尺度钢板三维测量拼接系统及方法 | |
CN104123726A (zh) | 基于消隐点的大锻件测量系统标定方法 | |
CN102306382A (zh) | 基于道路环境固有特征的事故现场图像校正方法 | |
CN104156563B (zh) | 基于多重欠采样实现光学遥感器像能量超频采样的方法 | |
CN107329190B (zh) | 静止气象卫星恒星敏感的成像试验方法 | |
CN106407987B (zh) | 基于图像互相关匹配的电光晶体光轴出露点提取方法 | |
CN103134664B (zh) | 一种基于凸面反射镜的在轨光学卫星相机mtf测量方法 | |
CN104931140A (zh) | 四次测量型全Stokes参数偏振仪的参数优化方法 | |
Ma et al. | Variable subset DIC algorithm for measuring discontinuous displacement based on pixel-level ZNCC value distribution map |
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 |