CN107274359B - 一种面阵用虚拟像元内插细分与信号平滑化的方法 - Google Patents
一种面阵用虚拟像元内插细分与信号平滑化的方法 Download PDFInfo
- Publication number
- CN107274359B CN107274359B CN201710374595.0A CN201710374595A CN107274359B CN 107274359 B CN107274359 B CN 107274359B CN 201710374595 A CN201710374595 A CN 201710374595A CN 107274359 B CN107274359 B CN 107274359B
- Authority
- CN
- China
- Prior art keywords
- pixel
- pixels
- dimensional
- sub
- adjacent
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000009499 grossing Methods 0.000 title claims abstract description 12
- 238000003491 array Methods 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 abstract description 10
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种面阵用虚拟像元内插细分与信号平滑化的方法。对于每个二维像元或者相邻邻接的四个二维像元分别内插构建子像元并计算子像元的光电信号值,内插细分后子像元的光电信号面密度依据内插前的相近像元的平均面密度计算,实现了对光电信号的平滑化处理;对于单个二维像元,内插构建小正方形的中心子像元,小正方形的中心子像元是由单个二维像元的四条边的边长中点相依次连接构成;相邻邻接的四个二维像元内插构建小正方形的连接子像元,小正方形的连接子像元是由相邻四个二维像元形成的四条共同边的边长中点相依次连接构成。本发明能使减少虚拟像元间距,能使信号显著地平滑化,减小了各像元原来转换率误差,进而减小峰位标准差。
Description
技术领域
本发明属于技术领域,具体涉及一种面阵用虚拟像元内插细分与信号平滑化的方法。
背景技术
由于制造技术限制及使用情况等原因,使面阵各个像元的几何尺寸相对于面阵像元平均尺寸产生具有一定随机性的定值误差,使各像元等效几何中心坐标相对于理想几何中心坐标产生具有一定随机性的定值误差,使各像元光电信号转换率相对于面阵信号转换率平均值也产生具有一定随机性的定值误差。
对于面阵采集的由干涉、衍射等机理形成的近似平行的直条纹或大曲率半径的同心圆环条纹,常常需要求出在条纹的准垂线方向某一线段上的光电信号分布规律,进而用回归或其它方法求出信号峰位坐标的细分后的值。对于近似平行的条纹或同心圆环条纹,在求垂直于条纹方向的光电信号峰值位置坐标(简称峰位)时,需利用这一方向上多个相邻像元的光电信号值,按光强分布规律用回归等方法求出条纹信号的峰位值及其标准差。
现有技术中由于像元间距和转换率误差等都受限于面阵生产工艺等因素的限制,到时求出条纹的信号的峰位不够准确,峰位标准差较大。
发明内容
为了解决背景技术中存在的问题,本发明目的在于提供一种面阵用虚拟像元内插细分与信号平滑化的方法。
为实现上述目的,本发明采用的技术方案如下:
原面阵由像元间隔为w的像元阵列构成时,对于每个二维像元或者相邻邻接在一起的四个二维像元分别内插构建子像元并计算子像元的光电信号值,内插细分后的子像元重构面阵中,子像元的光电信号面密度依据面阵上内插前的相近像元的平均面密度按一定规律计算,实现了对光电信号的平滑化处理。
本发明的内插细分后计算面阵的平均面密度时相当于对光电信号做了平滑化处理,减小了原始面阵信号中有随机性的未定系统误差分量的影响。
对于单个二维像元,内插构建小正方形的中心子像元,小正方形的中心子像元是由单个二维像元的四条边的边长中点相依次连接构成。中心子像元的光电信号值Zmid为:
其中,ZO表示为单个二维像元的光电信号值,ZW表示与信号值为ZO的单个二维像元周围四边分别相邻的四个二维像元的光电信号的总和。
对于相邻邻接在一起的四个二维像元,内插构建小正方形的连接子像元,小正方形的连接子像元是由相邻四个二维像元形成的四条共同边的边长中点相依次连接构成。对于沿条纹方向,连接子像元的光电信号值Zadj为:
其中,ZW1表示相邻四个二维像元中位于沿条纹方向两侧的两个二维像元的光电信号的总和,ZW2表示相邻四个二维像元中沿条纹方向的两个二维像元的光电信号的总和。
若要获得条纹上的光电信号分布规律,就要对条纹垂直方向某一直线上的一系列相邻的中心子像元和连接子像元为基本单位的光电信号值作统计分析,,就能获得条纹上的光电信号分布规律的特征值,例如用回归或其它方法能求出信号峰位坐标的细分后的值及其标准差。
本发明的有益效果:
本发明的虚拟像元内插细分能使信号显著地平滑化,减小了像元间距,又减小了各像元原来的转换率误差的影响,进而就能减小峰位标准差,减弱了面阵生产工艺不完善的影响。实施例中采用本发明后使峰位标准差平均减小到原来的45%。
附图说明
图1是在像元B2中内插小正方形虚拟像元IJKL的示意图。
图2是以像元B2与C3的接点为中心内插小正方形虚拟像元JNOK的示意图
图3是相邻虚拟像元IJKL和JNOK的示意图。
图4是实施例中面阵获取的干涉产生的同心圆环图。
具体实施方式
以下通过具体实施例对本发明的技术做详细说明,但不局限于此。
本发明的实施例如下:
如图1和图2所示,用短划线表示边界线的二维像元A1~A3、B1~B3和C1~C3组成的3×3面阵,像元间隔为w,面阵中各二维虚拟像元的光电信号值分别为ZA1~ZA3、ZB1~ZB3和ZC1~ZC3。
如图1和图2所示,需要求取光电信号的条纹经过二维像元对角线方向,为图示的y”轴方向,相对于原直角坐标系y轴逆时针转过45度角。
(1)对于单个二维虚拟像元
如图1所示,对于原像元B2构建内接小正方形IJKL(图1中的粗实线表示)作为中心子像元,
其中,ZB2表示像元B2的光电信号值。
(2)对于相邻邻接在一起的四个二维虚拟像元
如图2所示,在像元C3与像元B2之间内插一个小正方形JNOK像元作为连接子像元(图2中的粗实线表示)。小正方形JNOK的边长为小正方形JNOK的光电信号面密度取长方形FMPG光电信号面密度的平均值,小正方形JNOK光电信号值为:
综合以上两个式子,在斜转45度的一条对角线(y”轴向)上就能产生一系列相邻的边长为的小正方形虚拟像元,如正方形IJKL与正方形JNOK相邻,如图3所示。统计分析条纹垂直方向某一直线上的一系列相邻的中心子像元和连接子像元为基本单位的光电信号值,就能获得条纹上的光电信号分布规律的特征值,例如用回归或其它方法能求出信号峰位坐标的细分后的值及其标准差。
类似地,当条纹方向与x”方向垂直时,在x”方向上也作同样的虚拟像元内插细分。
如小正方形虚拟像元IJKL全部在的像元B2之内,像元B2的几何尺寸相对于面阵像元平均尺寸具有一定随机性的定值误差,B2的等效几何中心坐标相对于理想几何中心坐标具有一定随机性的定值误差,B2的光电信号转换率相对于面阵信号转换率平均值也具有一定随机性的定值误差。因此按照本发明实施例实施后面阵中虚拟像元的等效像元间隔为组成的新构面阵的列或行的像元数增加到原来的约倍,同一尺寸面阵的虚拟像元总数增加约1倍。
将综合以上两个式子(1)和(2)的内插子像元计算面阵的平均面密度,实现了对光电信号的平滑化处理,减小了原始面阵信号中有随机性的未定系统误差分量的影响。对于式(1)的计算,其中以像元B2的信号为主同时考虑四个相邻像元的信号值影响这样使式(1)的结果不仅仅依据B2的信号,减小了B2信号中有一定随机性的未定系统误差的影响,起到了平滑化的效果。对于式(2)的计算,其中虚拟像元小正方形JNOK的信号值,平滑化的效果比较直观。
由于像元间距愈小则峰位标准差愈小,本发明用虚拟像元内插细分能使虚拟像元间距减少到原来像元间距的由于光电信号转换率误差愈小则峰位标准差愈小,但转换率误差受限于面阵生产工艺的限制,虚拟像元内插细分时能使信号显著地平滑化,减小了各像元原来转换率误差的影响,减弱了面阵生产工艺的限制。
本实施例用波长λ≈546.1nm的光透过间距d≈2.0mm的F-P标准具,在焦距f≈75mm的物镜焦平面上形一系列同心圆环,如图4所示。面阵成像器件采用了奥林巴斯EN-F型相机的面阵,尺寸为17.4×13mm2,面阵像元数为10368×7776,面阵的平均像元间距为w≈1.675μm。对从内向外计数的第10个圆环,相对直径D10/w≈4753.15。
对平行于x,y轴方向各取一条近似直径(具体方法是:对于x轴方向或者y轴方向,找出第i个圆环上沿所在轴方向坐标值最大与最小的光电信号极值点,两个光电信号极值点的坐标值取平均后再取整数值作为所在轴向的近似圆心坐标值;过近似圆心坐标定出的近似圆心点作平行于x,y轴方向的直线段作为近似直径,直线段与圆环相交),在每条近似直径两侧约±35w的范围内取71条间隔为w相邻的平行线,两方向各71条平行线与圆环相交共形成284个小线段。对于284个线段的每个线段,采用申请日为申请号为201510217472.7的专利申请说明书或者PCT国际申请号为PCT/CN2016/078164的专利申请说明书中的对光强峰值坐标值的求取方法来求出每个线段上的峰位坐标。然后再根据所有线段的峰位坐标统计所有线段峰位坐标的标准差的平均值及方均根值计算结果为及
设将x,y轴逆时针旋转45度之后的新坐标轴为x",y"轴。对平行于x",y"轴方向,通过上文定出的近似圆心点各取一条近似直径。采用本发明方法,在每条近似直径两侧±35w的范围内取101条间隔为相邻的平行线,两方向各101条平行线与圆环相交共形成404个小线段。对于404个线段的每个线段,采用申请日为申请号为201510217472.7的专利申请说明书或者PCT国际申请号为PCT/CN2016/078164的专利申请说明书中的对光强峰值坐标值的求取方法来求出每个线段上的峰位坐标。然后再根据所有线段的峰位坐标统计所有线段峰位坐标的标准差的平均值及方均根值计算结果为及
Claims (2)
1.一种面阵用虚拟像元内插细分与信号平滑化的方法,其特征在于:原面阵由像元间隔为w像元阵列构成时,对于每个二维像元或者相邻邻接在一起的四个二维像元分别内插构建子像元并计算子像元的光电信号值,内插细分后子像元的光电信号面密度依据内插前的相近原像元的光电信号平均面密度计算,实现了对光电信号的平滑化处理;
对于单个二维像元,内插构建小正方形的中心子像元,小正方形的中心子像元是由单个二维像元的四条边的边长中点相依次连接构成,中心子像元的光电信号值Zmid为:
其中,ZO表示为单个二维像元的光电信号,ZW表示与信号值为ZO的单个二维像元周围四边相邻的四个二维像元的光电信号总和;
对于相邻邻接在一起的四个二维像元,内插构建小正方形的连接子像元,小正方形的连接子像元是由相邻四个二维像元形成的四条共同边的边长中点相依次连接构成;对于沿条纹方向,连接子像元的光电信号值Zadj为:
其中,ZW1表示相邻四个二维像元中位于沿条纹方向两侧的两个二维像元的光电信号总和,ZW2表示相邻四个二维像元中沿条纹方向的两个二维像元的光电信号总和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710374595.0A CN107274359B (zh) | 2017-05-24 | 2017-05-24 | 一种面阵用虚拟像元内插细分与信号平滑化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710374595.0A CN107274359B (zh) | 2017-05-24 | 2017-05-24 | 一种面阵用虚拟像元内插细分与信号平滑化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107274359A CN107274359A (zh) | 2017-10-20 |
CN107274359B true CN107274359B (zh) | 2020-06-30 |
Family
ID=60064191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710374595.0A Active CN107274359B (zh) | 2017-05-24 | 2017-05-24 | 一种面阵用虚拟像元内插细分与信号平滑化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107274359B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4422019A (en) * | 1982-07-12 | 1983-12-20 | Tektronix, Inc. | Apparatus for providing vertical as well as horizontal smoothing of convergence correction signals in a digital convergence system |
CN103279927A (zh) * | 2013-05-16 | 2013-09-04 | 中国电子科技集团公司第四十一研究所 | 一种多面阵图像生成方法及多面阵图像输出装置 |
CN106651784A (zh) * | 2016-10-11 | 2017-05-10 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于去除高斯固定图案噪声的非均匀校正方法 |
-
2017
- 2017-05-24 CN CN201710374595.0A patent/CN107274359B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4422019A (en) * | 1982-07-12 | 1983-12-20 | Tektronix, Inc. | Apparatus for providing vertical as well as horizontal smoothing of convergence correction signals in a digital convergence system |
CN103279927A (zh) * | 2013-05-16 | 2013-09-04 | 中国电子科技集团公司第四十一研究所 | 一种多面阵图像生成方法及多面阵图像输出装置 |
CN106651784A (zh) * | 2016-10-11 | 2017-05-10 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于去除高斯固定图案噪声的非均匀校正方法 |
Non-Patent Citations (2)
Title |
---|
Scene-based nonuniformity correction algorithm based on interframe registration;Chao Zuo, Qian Chen, Guohua Gu, and Xiubao Sui;《Journal of the Optical Society of America A》;20110520;第28卷(第6期);1164-1176 * |
自适应虚拟电子微扫描红外焦平面像素倍增技术;顾国华等;《光电子技术》;20041231;第24卷(第4期);243-253 * |
Also Published As
Publication number | Publication date |
---|---|
CN107274359A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107014312B (zh) | 一种振镜式线激光扫描三维测量系统的整体标定方法 | |
CN111476242B (zh) | 一种激光点云语义分割方法及装置 | |
CN112381847B (zh) | 管路端头空间位姿测量方法及系统 | |
CN109313807B (zh) | 使用光场相机的三维重建方法 | |
CN104380342A (zh) | 图像处理设备、成像设备以及图像处理方法 | |
CN104320565B (zh) | 一种多镜头探测器阵列曲面像面拼接方法 | |
CN108692655A (zh) | 基于光场成像的三维测量方法 | |
CN104240229B (zh) | 一种红外双目相机自适应极线校正方法 | |
CN111182179A (zh) | 奇偶透镜线阵交替分布的分段平面侦察成像系统和方法 | |
CN115205466B (zh) | 一种基于结构光的输电通道三维重建方法和系统 | |
JP2016218115A (ja) | 光学素子の設計方法、光学素子アレイ、センサアレイおよび撮像装置 | |
CN103440619B (zh) | 一种斜模式采样建模与超分辨率重建方法 | |
CN104318583A (zh) | 一种可见光宽带光谱图像配准方法 | |
KR20150063010A (ko) | 플렌옵틱 카메라로 취득된 장면의 뷰와 연관된 시차를 추정하기 위한 방법 및 디바이스 | |
CN107274359B (zh) | 一种面阵用虚拟像元内插细分与信号平滑化的方法 | |
CN108322736B (zh) | 一种多个线阵相机绕视轴旋转角度校准的校准板和校准方法 | |
JP7169363B2 (ja) | 撮像装置 | |
CN101929836B (zh) | 一种物体三维定位方法及摄像机 | |
CN110708532B (zh) | 一种普适光场单元图像生成方法和系统 | |
JP7308045B2 (ja) | 立体画像生成システム、並びに、立体画像生成装置及びそのプログラム | |
CN115471616A (zh) | 焊缝三维重建方法、装置、处理器以及系统 | |
CN111089563B (zh) | 平面斜视图像测距定位方法及测距定位系统 | |
CN107767341A (zh) | 基于阵列镜头的图像超分辨率重建方法 | |
CN107101583A (zh) | 基于纵横转换放大光栅尺的图像处理方法及其应用 | |
CN110473261B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |