CN107016210B - 一种菲涅尔衍射光学系统的成像质量仿真方法 - Google Patents
一种菲涅尔衍射光学系统的成像质量仿真方法 Download PDFInfo
- Publication number
- CN107016210B CN107016210B CN201710255306.5A CN201710255306A CN107016210B CN 107016210 B CN107016210 B CN 107016210B CN 201710255306 A CN201710255306 A CN 201710255306A CN 107016210 B CN107016210 B CN 107016210B
- Authority
- CN
- China
- Prior art keywords
- image
- diffraction
- image block
- formula
- mirror
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Eyeglasses (AREA)
Abstract
一种菲涅尔衍射光学系统的成像质量仿真方法,本发明属于光学成像仿真技术领域,所述方法步骤如下:步骤一:基于等晕区分块思想,将原始图像分割为若干个具有近似相同点扩散函数的图像块;步骤二:计算各图像块中心的子午面位置及其对应的视场角;步骤三:构造离散化的菲涅尔衍射公式,计算各图像块的点扩散函数;步骤四:利用图像不同区域的MTF对各图像块进行低通滤波,并添加随机噪声。本发明适用于以衍射镜为主镜的光学系统的成像质量仿真,可仿真生成具有低MTF、低对比度、低SNR、大尺寸PSF空变等衍射光学成像特性的图像,为菲涅尔衍射成像系统的优化设计及其图像处理算法研究与验证提供支持。
Description
技术领域
本发明属于光学成像仿真技术领域,具体涉及一种菲涅尔衍射光学系统的成像质量仿真方法。
背景技术
光学遥感成像质量仿真预估在侦察、测绘、监视等遥感任务预测、成像系统优化设计与性能评估以及图像处理算法验证等方面具有重要的指导意义和应用价值。随着高轨监视对空间分辨力需求的不断提高,迫切需要空间光学载荷同时具有超大口径、轻量化、加工周期短、成本低等特点,传统的反射式或折射式成像系统已无法全面满足这些要求,而衍射光学成像技术能够克服这一局限,因此急需针对光学衍射成像系统开展成像质量的仿真模型与方法研究。但从现有文献看,仿真方法方面的研究多是针对于反射式或折射式系统,关于衍射成像系统像质退化机理、仿真模型与方法的研究少见报道,而与传统载荷相比,衍射光学成像系统的图像质量退化更加严重,其图像的传递函数(MTF)和信噪比(SNR)均较低,特别是不同视场处的点扩散函数(PSF)空间变化性大,显然无法采用传统的反射式或折射式系统模型进行仿真。因此,从衍射光学系统像质退化机理出发,开展成像质量仿真模型与方法研究,对衍射光学成像系统的实际空间应用具有重要意义。
发明内容
本发明的目的是针对目前衍射光学成像系统成像质量仿真方法理论研究的空白,提出一种菲涅尔衍射光学系统的成像质量仿真方法。本发明适用于以衍射镜为主镜的光学系统的成像仿真,可仿真生成具有衍射成像特性的图像,为衍射光学成像系统的优化设计以及图像处理算法的研究与验证提供支持。
为实现上述目的,本发明采取的技术方案如下:
一种菲涅尔衍射光学系统的成像质量仿真方法,所述方法步骤如下:
步骤一:基于等晕区分块思想,将原始图像分割为若干个具有近似相同点扩散函数的图像块;
步骤二:计算各图像块中心的子午面位置及其对应的视场角;
步骤三:构造离散化的菲涅尔衍射公式,计算各图像块的点扩散函数;
步骤四:利用图像不同区域的MTF对各图像块进行低通滤波,并添加随机噪声。
本发明相对于现有技术的有益效果是:
(1)本发明提出的菲涅尔衍射光学系统的成像质量仿真方法,适用于以衍射镜为主镜的光学系统成像质量的仿真预估,能够获取具有低MTF、低对比度、低SNR、PSF空间移变大等衍射光学特性的仿真图像,可为衍射光学成像系统设计以及图像处理算法研究与验证提供必要的支持,填补了目前衍射光学成像仿真理论方法的空白。
(2)本发明将像元尺寸、镜面采样间隔、视场角等成像系统参量引入菲涅尔衍射公式,并采用等效衍射镜透过率函数同时结合观察平面坐标系变换,构造了衍射光学成像系统的PSF模型,该模型能够反映不同成像视场、方向以及衍射镜透过率、焦距、像元尺寸等参数下的PSF变化规律,从而更准确地揭示衍射光学成像系统的像质退化规律。
(3)本发明基于等晕区分块思想,将原始图像分割为若干个具有近似相同点扩散函数的图像块,并设计适当的重叠区域,针对相邻图像块拼接过程中可能造成重叠区域边缘附近产生振铃等伪像的问题,构建基于重叠区域图像灰度加权组合的伪像抑制策略,并利用由离散化菲涅尔衍射公式计算得出的点扩散函数对其进行低通滤波,从而能够获得具有大尺寸PSF空变、低对比度等衍射光学成像特性的仿真图像。
综上,本发明提供了一种菲涅尔衍射光学系统的成像质量仿真方法,适用于以衍射镜为主镜的光学系统的成像仿真,获取具有大尺寸PSF空变、低对比度等衍射光学成像特性的图像,可支持衍射光学成像系统的优化设计,同时还能够为衍射光学成像系统图像处理算法的研究与验证提供输入图像。
附图说明
图1为一种菲涅尔衍射光学系统的成像质量仿真方法流程图;
图2为子午面与视场角的几何示意图;
图3为衍射传播模型示意图;
图4为倾斜入射时光的传播示意图;
图5为观察平面位置转换示意图;
图6为坐标系变换示意图;
图7为中心视场PSF示意图;
图8为视场角2°×2°时PSF示意图;
图9为视场角-1°×4°时PSF示意图;
图10为图像不同区域PSF示意图;
图11为遥感图像(一);
图12为遥感图像(一)仿真结果图;
图13为遥感图像(二);
图14为遥感图像(二)仿真结果图;
具体实施方式
下面结合实施例和附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
具体实施方式一:如图1所示,本实施方式记载的是一种菲涅尔衍射光学系统的成像质量仿真方法,所述方法步骤如下:
步骤一:基于等晕区分块思想,将原始图像分割为若干个具有近似相同点扩散函数的图像块;
步骤二:计算各图像块中心的子午面位置及其对应的视场角;
步骤三:构造离散化的菲涅尔衍射公式,计算各图像块的点扩散函数;
步骤四:利用图像不同区域的MTF对各图像块进行低通滤波,并添加随机噪声。
具体实施方式二:如图1所示,具体实施方式一所述的一种菲涅尔衍射光学系统的成像质量仿真方法,所述步骤一的具体步骤如下:
(1)基于等晕区分块思想,将原始图像分割为若干个具有近似相同点扩散函数的图像块:
菲涅尔衍射成像系统具有明显的PSF空间移变特性,其退化模型可表示为:
Y=H(ωx,ωy)X+Noi
式中,Y为退化图像频谱;H(ωx,ωy)表示视场角为(ωx,ωy)时,光学系统传递函数(MTF);Noi为噪声频谱;X表示原始图像频谱;
(2)先将若干个具有近似相同点扩散函数的图像块的相邻图像块间设计重叠区域;以图像中心为原点,建立正交直角坐标系,标记每个图像块的中心坐标为(xa,yb)。
具体实施方式三:如图1所示,具体实施方式二所述的一种菲涅尔衍射光学系统的成像质量仿真方法,所述步骤二的具体步骤如下:
(1)计算各图像块在x轴、y轴方向视场角:
设衍射光学成像系统焦距为f,像元尺寸为P1,则中心坐标为(xa,yb)的图像块在x、y方向视场角ωx和ωy分别为:
(2)计算各图像块在子午面上视场角:
如图2所示,由空间几何关系可知子午面上视场角α和子午面与y轴夹角β分别为:
式中,sign(·)为符号函数。
具体实施方式四:如图1所示,具体实施方式三所述的一种菲涅尔衍射光学系统的成像质量仿真方法,所述步骤三的具体步骤如下:
(1)将菲涅尔衍射公式离散化,对镜面透过率函数进行采样,取采样间隔为P2mm,则平行光垂直入射时距离衍射镜zmm处的衍射场复振幅分布为:
式中,坐标系x0-y0和x-y分别位于光学系统的衍射平面和观察平面,如图3所示,且x,y,x0,y0均为整数;U表示观察平面光场复振幅分布;t为衍射镜复振幅透过率采样矩阵,其大小为M×N;M表示矩阵的行数,N表示矩阵的列数;镜面和像面坐标原点均位于矩阵中心;λ为入射光波长;j表示虚数;波数k=2π/λ,π为圆周率;
(2)斜入射时成像模型如图4所示,此时衍射镜的等效复振幅透过率t′相当于原透过率t在与入射方向垂直平面上的投影,设衍射镜透过率函数为:
(3)利用离散化的菲涅尔衍射公式可得到,图5中I′平面上有:
式中,坐标系原点为主光线与坐标平面交点;U′表示I′平面光场复振幅分布。但实际像面位于I平面,而不是I′。考虑到空间内某点的复振幅与观察屏位置无关,那么,在I平面和I′平面,O点具有相同的复振幅,即:
U(0,0)=U′(0,0)
由此推理,图5中I平面上任意一点Q(x,y)的复振幅,可以在I1平面求解,具体方法如下:如图6所示,将坐标轴旋转,得到Q点在弧矢-子午坐标系中坐标(lx,ly),即像面上任意一点Q(x,y)在弧矢和子午方向上投影lx和ly分别为:
lx=xcosβ-ysinβ
ly=xsinβ+ycosβ
则,I1平面作用距离z1表示为:
z1=z′-lysinα
此外,Q点在I1平面弧矢-子午坐标系中的坐标(lx1,ly1)为:
l1x=lx,l1y=lycosα
故,Q点在I1平面x-y坐标下坐标(x1,y1)为:
x1=lxcosβ+lycosαsinβ
y1=-lxsinβ+lycosαcosβ
则Q点与沿着光线入射方向作用距离为z1的观察平面上的Q1(x1,y1)点的光强相等;
(4)将等效条件代入离散化的菲涅尔衍射公式得:
式中,U(ωx,ωy;x,y)表示视场角为ωx×ωy时,衍射光学成像系统的振幅扩散函数,取其模的平方可得:
h(ωx,ωy;x,y)=|U(ωx,ωy;x,y)|2
式中,h(ωx,ωy;x,y)为点扩散函数,作傅里叶变换可得该视场下衍射主镜成像传递函数Hdiff,即:
Hdiff(ωx,ωy)=F{h(ωx,ωy;x,y)}
式中,F{·}表示傅里叶变换符;衍射主镜后端光学系统可视为线性空不变系统,设其传递函数为Haft,则衍射光学成像系统传不同视场递函数H(ωx,ωy)可表示为:
H(ωx,ωy)=Hdiff(ωx,ωy)·Haft
图7至图9是不同视场角下的点扩散函数仿真结果示意图。
具体实施方式五:如图1所示,具体实施方式四所述的一种菲涅尔衍射光学系统的成像质量仿真方法,所述步骤四的具体步骤如下:
(1)利用图像不同区域的MTF对各图像块进行低通滤波,并添加随机噪声:利用具体实施方式三的步骤二和具体实施方式四的步骤三可得到图像各个位置上的传递函数,如图10;将各图像块分别使用相应的MTF进行低通滤波:
Yab=H(ωx,ωy)Xab
式中,Xab为原始图像块频谱;Yab表示退化频谱,对其作傅里叶逆变换可得各图像块Imgab的低通滤波结果;
(2)将所有图像块滤波结果拼合成完整图像,记为Img′,重叠区域采用相邻图像块灰度的加权组合,距离某图像块中心距离越近,则该图像块灰度所占权重越大;最后对图像整体添加噪声,即可得到最终菲涅尔衍射成像系统图像仿真结果:
Imgdiff=Img′+n
式中,n表示随机噪声;Imgdiff为最终仿真结果。仿真结果展示如图11~图14所示。
Claims (2)
1.一种菲涅尔衍射光学系统的成像质量仿真方法,其特征在于:所述方法步骤如下:
步骤一:基于等晕区分块思想,将原始图像分割为若干个具有相同点扩散函数的图像块;所述步骤一的具体步骤如下:先将若干个具有相同点扩散函数的图像块的相邻图像块间设计重叠区域;以图像中心为原点,建立正交直角坐标系,标记每个图像块的中心坐标为(xa,yb);
步骤二:计算各图像块中心的子午面位置及其对应的视场角;所述步骤二的具体步骤如下:
(1)计算各图像块在x轴、y轴方向视场角:
设衍射光学成像系统焦距为f,像元尺寸为P1,则中心坐标为(xa,yb)的图像块在x、y方向视场角ωx和ωy分别为:
(2)计算各图像块在子午面上视场角:
由空间几何关系可知子午面上视场角α和子午面与y轴夹角β分别为:
式中,sign(·)为符号函数;
步骤三:构造离散化的菲涅尔衍射公式,计算各图像块的点扩散函数;所述步骤三的具体步骤如下:
(1)将菲涅尔衍射公式离散化,对镜面透过率函数进行采样,取采样间隔为P2mm,则平行光垂直入射时距离衍射镜zmm处的衍射场复振幅分布为:
式中,坐标系x0-y0和x-y分别位于光学系统的衍射平面和观察平面,且x,y,x0,y0均为整数;U表示观察平面光场复振幅分布;t为衍射镜复振幅透过率采样矩阵,其大小为M×N;M表示矩阵的行数,N表示矩阵的列数;镜面和像面坐标原点均位于矩阵中心;λ为入射光波长;j表示虚数;波数k=2π/λ,π为圆周率;
(2)设衍射镜透过率函数为:
(3)像面上任意一点Q(x,y)在弧矢和子午方向上投影lx和ly分别为:
lx=xcosβ-ysinβ
ly=xsinβ+ycosβ
在I1平面x-y坐标下有:
x1=lxcosβ+lycosαsinβ
y1=-lxsinβ+lycosαcosβ
z1=z'-lysinα
式中,(x1,y1)为Q点在I1平面x-y坐标下的坐标;z1为沿着光线入射方向作用距离;z'为I'平面内点的z轴坐标;则Q点光强可利用Q1(x1,y1)点求得;
(4)将等效条件代入离散化的菲涅尔衍射公式得:
式中,U(ωx,ωy;x,y)表示视场角为ωx×ωy时,衍射光学成像系统的振幅扩散函数,取其模的平方可得:
h(ωx,ωy;x,y)=|U(ωx,ωy;x,y)|2
式中,h(ωx,ωy;x,y)为点扩散函数,作傅里叶变换可得该视场下衍射主镜成像传递函数Hdiff,即:
Hdiff(ωx,ωy)=F{h(ωx,ωy;x,y)}
式中,F{·}表示傅里叶变换符;衍射主镜后端光学系统可视为线性空不变系统,设其传递函数为Haft,则衍射光学成像系统传不同视场递函数H(ωx,ωy)可表示为:
H(ωx,ωy)=Hdiff(ωx,ωy)·Haft;
步骤四:利用图像不同区域的MTF对各图像块进行低通滤波,并添加随机噪声。
2.根据权利要求1所述的一种菲涅尔衍射光学系统的成像质量仿真方法,其特征在于:所述步骤四的具体步骤如下:
(1)将各图像块分别使用相应的MTF进行低通滤波:
Yab=H(ωx,ωy)Xab
式中,Xab为原始图像块频谱;Yab表示退化频谱,对其作傅里叶逆变换可得各图像块的低通滤波结果;
(2)将各个图像块滤波结果拼合成完整图像,记为Img',重叠区域采用相邻图像块灰度的加权组合,距离某图像块中心距离越近,则该图像块灰度所占权重越大;最后对图像整体添加噪声,即可得到最终菲涅尔衍射成像系统图像仿真结果:
Imgdiff=Img'+n
式中,n表示随机噪声;Imgdiff为最终仿真结果。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2016111963079 | 2016-12-21 | ||
CN201611196307 | 2016-12-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016210A CN107016210A (zh) | 2017-08-04 |
CN107016210B true CN107016210B (zh) | 2020-08-07 |
Family
ID=59448363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710255306.5A Active CN107016210B (zh) | 2016-12-21 | 2017-04-18 | 一种菲涅尔衍射光学系统的成像质量仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016210B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108364250A (zh) * | 2018-01-30 | 2018-08-03 | 哈尔滨工业大学 | 一种基于矢量衍射理论的菲涅尔衍射光学系统成像仿真方法 |
CN109459589B (zh) * | 2018-09-26 | 2020-04-17 | 东莞市凯融光学科技有限公司 | 一种测量等高菲涅尔结构产品面型的方法 |
CN111428095B (zh) | 2020-06-11 | 2020-08-28 | 上海冰鉴信息科技有限公司 | 图数据质量验证方法及图数据质量验证装置 |
CN113219650B (zh) * | 2021-03-25 | 2022-08-12 | 北京空间机电研究所 | 一种高分辨率大视场空间光学遥感器的设计方法 |
CN114859550B (zh) * | 2022-05-13 | 2022-11-22 | 哈尔滨工业大学 | 一种针对菲涅尔单透镜计算成像系统的端到端设计方法 |
CN117045279B (zh) * | 2023-09-22 | 2024-03-29 | 哈尔滨工业大学 | 一种基于点扩散函数参数优化的掌上超声成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102499648A (zh) * | 2011-11-16 | 2012-06-20 | 清华大学 | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 |
CN102798472A (zh) * | 2012-05-23 | 2012-11-28 | 北京理工大学 | 基于fpa非制冷红外热成像的菲涅耳透镜光学读出系统 |
CN204666945U (zh) * | 2015-04-30 | 2015-09-23 | 北京空间机电研究所 | 一种采用直角反射棱镜的双通道成像光学系统 |
CN105572895A (zh) * | 2015-12-18 | 2016-05-11 | 哈尔滨工业大学 | 大口径谐衍射主镜成像系统及其应用 |
-
2017
- 2017-04-18 CN CN201710255306.5A patent/CN107016210B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102499648A (zh) * | 2011-11-16 | 2012-06-20 | 清华大学 | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 |
CN102798472A (zh) * | 2012-05-23 | 2012-11-28 | 北京理工大学 | 基于fpa非制冷红外热成像的菲涅耳透镜光学读出系统 |
CN204666945U (zh) * | 2015-04-30 | 2015-09-23 | 北京空间机电研究所 | 一种采用直角反射棱镜的双通道成像光学系统 |
CN105572895A (zh) * | 2015-12-18 | 2016-05-11 | 哈尔滨工业大学 | 大口径谐衍射主镜成像系统及其应用 |
Non-Patent Citations (3)
Title |
---|
Optical testing of a membrane diffractive optic for space-based solar imaging;Olha Asmolova等;《PROCEEDINGS OF SPIE》;20141231;第1-6页 * |
偏振对高NA光学系统成像影响的研究;庞武斌;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20131015(第10期);第3.4节 * |
成像光谱仪光学系统设计与像质评价研究;方煜;《中国博士学位论文全文数据库基础科学辑》;20140615(第06期);第4.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN107016210A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016210B (zh) | 一种菲涅尔衍射光学系统的成像质量仿真方法 | |
US10127687B2 (en) | Calibration device, calibration method, optical device, image-capturing device, projection device, measuring system, and measuring method | |
CN103201648B (zh) | 重叠关联成像术中探测的校准 | |
CN102509261B (zh) | 一种鱼眼镜头的畸变校正方法 | |
CN104155765B (zh) | 在拼接式集成成像显示器中校正三维图像的方法和设备 | |
CN103852060B (zh) | 一种基于单目视觉的可见光图像测距方法 | |
JP4873485B2 (ja) | 多数の基準面を用いた形状計測方法および形状計測装置 | |
CN101755190A (zh) | 校准方法、校准装置及具备该校准装置的校准系统 | |
CN114174892A (zh) | 光瞳扩展器校准 | |
CN105528500A (zh) | 一种分米级星载tdi ccd立体测绘相机成像仿真方法和系统 | |
CN109799073A (zh) | 一种光学畸变测量装置及方法、图像处理系统、电子设备和显示设备 | |
CN109544643A (zh) | 一种摄像机图像校正方法及装置 | |
CN109615664A (zh) | 一种用于光学透视增强现实显示器的标定方法与设备 | |
CN105306922A (zh) | 一种深度相机参考图的获取方法和装置 | |
Hayasaka et al. | Mobile visualization of density fields using smartphone background-oriented schlieren | |
CN106169076A (zh) | 一种基于透视变换的角度车牌图像库搭建方法 | |
Verbiest et al. | Modeling the effects of windshield refraction for camera calibration | |
CN115375775A (zh) | 基于前向投影折射模型的水下相机标定方法及装置 | |
CN110274752A (zh) | 中继镜成像质量的多功能测试卡及其测试方法 | |
CN105547286B (zh) | 一种复合三视场星敏感器星图模拟方法 | |
JP6899915B2 (ja) | 標高データグリッドのためのシャドウキャスティング | |
CN103208120B (zh) | 切向和径向双近似圆综合的全景环带图像校正展开方法 | |
CN102034234A (zh) | 一种多视图结构光系统自标定方法 | |
CN108196437B (zh) | 一种全息成像显示方法、装置、设备及存储介质 | |
CN110231021B (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 |