CN113295286A - 一种夏克-哈特曼波前传感器测量系统及相位恢复方法 - Google Patents
一种夏克-哈特曼波前传感器测量系统及相位恢复方法 Download PDFInfo
- Publication number
- CN113295286A CN113295286A CN202110552221.XA CN202110552221A CN113295286A CN 113295286 A CN113295286 A CN 113295286A CN 202110552221 A CN202110552221 A CN 202110552221A CN 113295286 A CN113295286 A CN 113295286A
- Authority
- CN
- China
- Prior art keywords
- image
- shack
- wavefront
- phase
- unit
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000011084 recovery Methods 0.000 title claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 238000010587 phase diagram Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000000926 separation method Methods 0.000 abstract description 2
- 230000010354 integration Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J2009/002—Wavefront phase distribution
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明公开了一种夏克‑哈特曼波前传感器测量系统及相位恢复方法,解决测量、处理、显示过程分离,测量效率低的问题,夏克‑哈特曼波前传感器测量系统包括夏克‑哈特曼波前传感器、GPU核心板、显示单元,夏克‑哈特曼波前传感器包括微透镜阵列和CMOS单元,GPU核心板包括DDR存储单元、GPU单元和CPU,通过GPU核心板经PCIE接口接收从CMOS单元传输的图像,对图像进行处理、计算,经显示单元显示;夏克‑哈特曼波前传感器相位恢复方法通过夏克‑哈特曼波前传感器实时采集图像,GPU核心板对采集的斑点图像经确定计算轮廓、归一化图像、阈值分割、质心计算、波前相位重构等步骤,实现高速并行计算,最终将带有波前起伏信息的原始波前相位实时显示在显示器上。
Description
技术领域
本发明涉及精密测量仪器领域,特别涉及一种夏克-哈特曼波前传感器测量系统及相位恢复方法。
背景技术
夏克-哈特曼波前传感器是一种能够检测波面形状的仪器,在大气测量、望远测量、远距离通信、等离子体测量等领域得到应用。
传统的夏克-哈特曼波前传感器的测量与处理往往是分离的;传统的夏克-哈特曼波前传感器在测量后,需要通过网线连接到上位机上进行处理并显示,并且上位机在处理中的计算速度慢影响测量效率,这通常会限制夏克-哈特曼波前传感器的应用范围,是一个难以克服的问题。
发明内容
本发明为解决现有技术存在夏克-哈特曼波前传感器的测量、处理、显示分离,测量效率低的问题,提出了一种高速集成化的夏克-哈特曼波前测量系统及方法。
本发明的技术解决方案是:
一种夏克-哈特曼波前传感器测量系统,其特殊性在于:
包括夏克-哈特曼波前传感器、GPU核心板、显示单元;
所述夏克-哈特曼波前传感器包括微透镜阵列和CMOS单元,微透镜阵列将被测光场图像调制成斑点/点阵图像,CMOS单元用于采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像;
所述GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构后,得到带有波前起伏信息的相位图,其包括相互连接的DDR存储单元和GPU单元、CPU,所述DDR存储单元用于图像数据的存储,所述GPU单元通过HDMI接口与显示单元连接,所述CPU用于波前相位重构;
所述CMOS单元通过PCIE接口与GPU单元连接;
所述显示单元用于显示波前相位重构得到的带有波前起伏信息的相位图。
进一步地,所述微透镜阵列为56*44的方形阵列。
进一步地,所述COMS单元的型号为安森美5000。
一种夏克-哈特曼波前传感器相位恢复方法,其特殊性在于:包括以下步骤:
1)微透镜阵列将实时采集的被测光场图像调制成斑点/点阵图像;
2)CMOS单元采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像I(x,y);
3)GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构:
3.1)遍历斑点图像I(x,y),得到像素最大值Imax(x,y)、像素最小值Imin(x,y),并以像素最大值Imax(x,y)坐标为斑点图像I(x,y)的中心;
3.2)以斑点图像I(x,y)的中心为圆心,以500个像素为半径,得到斑点图像的圆形计算轮廓;
3.3)对圆形计算轮廓内的斑点图像I(x,y)的每个像素乘以Imin(x,y)/Imax(x,y)系数,得到归一化图像;
3.4)取0.15*Imax(x,y)作为阈值,将归一化图像中小于阈值的像素值全部置零,完成阈值分割;
3.5)对阈值分割后的归一化图像按照透镜阵列的阵列数进行方块划分,然后同时对划分的所有方块进行质心计算,划分后的方块中每个像素记为当前像素,得到每个方块中的质心坐标(xfc,yfc),所述质心计算公式为:
3.6)计算每个方块中的质心坐标(xfc,yfc)相对于该方块的中心,计算质心坐标的偏移(Δxi,Δyi);
3.8)由波前斜率Sx,Sy与相位之间的关系式S=CΦ,求解出带有波前起伏信息的原始波前相位Φ=C+S;其中C为系数,C+为采用最小二乘解得到的系数C值;
4)将带有波前起伏信息的原始波前相位实时显示在显示器上。
本发明与现有技术相比,其有益效果如下:
(1)本发明的夏克-哈特曼波前传感器测量系统,相比传统的夏克-哈特曼波前传感器测量系统需要额外配置计算机进行上位机数据计算和结果显示,通过GPU核心板经PCIE接口接收从CMOS单元传输的图像,对图像进行处理、计算,提升了数据吞吐效率,使速度可以达到5GT/s,该夏克-哈特曼波前传感器测量系统实现了测量、计算与显示一体化,集成化程度高,提高测量效率。
(2)本发明的夏克-哈特曼波前传感器相位恢复方法通过夏克-哈特曼波前传感器实时采集图像,GPU核心板对采集的斑点图像经确定计算轮廓、归一化图像、阈值分割、质心计算、波前相位重构等步骤,实现高速并行计算,最终,将带有波前起伏信息的原始波前相位实时显示在显示器,该夏克-哈特曼波前传感器相位恢复方法实现高速的相位恢复。
附图说明
图1为本发明夏克-哈特曼波前传感器测量系统的结构示意图;
图2为本发明夏克-哈特曼波前传感器相位恢复方法的流程图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解为本发明的保护范围并不受具体实施方式的限制。
如图1所示,该夏克-哈特曼波前传感器测量系统,包括夏克-哈特曼波前传感器、GPU核心板、显示单元;
夏克-哈特曼波前传感器包括微透镜阵列和CMOS单元,微透镜阵列将被测光场图像调制成斑点/点阵图像,微透镜阵列为56*44的矩形阵列,CMOS单元用于采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像,COMS单元的型号为安森美5000;
GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构后,得到带有波前起伏信息的相位图,其包括相互连接的DDR存储单元和GPU单元、CPU,DDR存储单元用于图像数据的存储,GPU单元通过HDMI接口与显示单元连接,CPU用于波前相位重构;
CMOS单元通过PCIE接口与GPU单元连接;
显示单元用于显示波前相位重构得到的带有波前起伏信息的相位图。
工作原理:夏克-哈特曼波前传感器中的探测器实时采集被测光场,经过夏克-哈特曼波前传感器中微透镜阵列调制成为点阵图像,通过PCIE接口传输到GPU核心板中,在GPU单元中进行确定计算轮廓、归一化图像、阈值分割、质心计算、波前相位重构等步骤,达到高速并行计算,最后将波前相位结果通过HDMI接口传输到显示单元,并实时显示。
该夏克-哈特曼波前传感器相位恢复方法,包括以下步骤:
1)微透镜阵列将实时采集的被测光场图像调制成斑点/点阵图像;
2)CMOS单元采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像I(x,y);
3)GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构:
3.1)遍历斑点图像I(x,y),得到像素最大值Imax(x,y)、像素最小值Imin(x,y),并以像素最大值Imax(x,y)坐标为斑点图像I(x,y)的中心;
3.2)以斑点图像I(x,y)的中心为圆心,以500个像素为半径,得到斑点图像的圆形计算轮廓;
3.3)对圆形计算轮廓内的斑点图像I(x,y)的每个像素乘以Imin(x,y)/Imax(x,y)系数,得到归一化图像;
3.4)取0.15*Imax(x,y)作为阈值,将归一化图像中小于阈值的像素值全部置零,完成阈值分割;
3.5)对阈值分割后的归一化图像按照透镜阵列的阵列数进行方块划分,然后同时对划分的所有方块进行质心计算,划分后的方块中每个像素记为当前像素,得到每个方块中的质心坐标(xfc,yfc),质心计算公式为:
式中:xfc为质心的x坐标,yfc为质心的y坐标,xf为当前像素的x坐标,yf为当前像素的y坐标,Ifl为当前像素的像素值;
3.6)计算每个方块中的质心坐标(xfc,yfc)相对于该方块的中心,计算质心坐标的偏移(Δxi,Δyi);
3.8)由波前斜率Sx,Sy与相位之间的关系式S=CΦ,求解出带有波前起伏信息的原始波前相位Φ=C+S;其中C为系数,C+为采用最小二乘解得到的系数C值;
4)将带有波前起伏信息的原始波前相位实时显示在显示器上。
以上公开的仅为本发明的具体实施例,但是,本发明实施例并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。
Claims (4)
1.一种夏克-哈特曼波前传感器测量系统,其特征在于:
包括夏克-哈特曼波前传感器、GPU核心板、显示单元;
所述夏克-哈特曼波前传感器包括微透镜阵列和CMOS单元,微透镜阵列将被测光场图像调制成斑点/点阵图像,CMOS单元用于采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像;
所述GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构后,得到带有波前起伏信息的相位图,其包括相互依次连接的DDR存储单元、GPU单元、CPU,所述DDR存储单元用于图像数据的存储,所述GPU单元用于进行质心计算,所述CPU用于波前相位重构;
所述GPU单元通过HDMI接口与显示单元连接,所述CMOS单元通过PCIE接口与GPU单元连接;
所述显示单元用于显示波前相位重构得到的带有波前起伏信息的相位图。
2.如权利要求1所述的夏克-哈特曼波前传感器测量系统,其特征在于:所述微透镜阵列为56*44的矩形阵列。
3.如权利要求2所述的夏克-哈特曼波前传感器测量系统,其特征在于:所述COMS单元的型号为安森美5000。
4.一种夏克-哈特曼波前传感器相位恢复方法,其特征在于:包括以下步骤:
1)微透镜阵列将实时采集的被测光场图像调制成斑点/点阵图像;
2)CMOS单元采集微透镜阵列调制后的带有波前相位信息的被测光场的斑点图像I(x,y);
3)GPU核心板对带有波前相位信息的被测光场的斑点图像进行质心计算和波前相位重构:
3.1)遍历斑点图像I(x,y),得到像素最大值Imax(x,y)、像素最小值Imin(x,y),并以像素最大值Imax(x,y)坐标为斑点图像I(x,y)的中心;
3.2)以斑点图像I(x,y)的中心为圆心,以500个像素为半径,得到斑点图像的圆形计算轮廓;
3.3)对圆形计算轮廓内的斑点图像I(x,y)的每个像素乘以Imin(x,y)/Imax(x,y)系数,得到归一化图像;
3.4)取0.15*Imax(x,y)作为阈值,将归一化图像中小于阈值的像素值全部置零,完成阈值分割;
3.5)对阈值分割后的归一化图像按照透镜阵列的阵列数进行方块划分,然后同时对划分的所有方块进行质心计算,得到每个方块中的质心坐标(xfc,yfc),所述质心计算公式为:
式中:xfc为质心的x坐标,yfc为质心的y坐标,xf为每个方块中任意像素的x坐标,yf为每个方块中任意像素的的y坐标,Ifl为每个方块中任意像素的像素值;
3.6)计算每个方块中的质心坐标(xfc,yfc)相对于该方块的中心的质心坐标偏移量(Δxi,Δyi);
3.8)由波前斜率Sx,Sy与相位之间的关系式S=CΦ,求解出带有波前起伏信息的原始波前相位Φ=C+S;其中C为系数,C+为采用最小二乘解得到的系数C值;
4)将带有波前起伏信息的原始波前相位实时显示在显示器上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110552221.XA CN113295286A (zh) | 2021-05-20 | 2021-05-20 | 一种夏克-哈特曼波前传感器测量系统及相位恢复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110552221.XA CN113295286A (zh) | 2021-05-20 | 2021-05-20 | 一种夏克-哈特曼波前传感器测量系统及相位恢复方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113295286A true CN113295286A (zh) | 2021-08-24 |
Family
ID=77323164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110552221.XA Pending CN113295286A (zh) | 2021-05-20 | 2021-05-20 | 一种夏克-哈特曼波前传感器测量系统及相位恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113295286A (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020027640A1 (en) * | 2000-04-25 | 2002-03-07 | Campin John Alfred | Spatial filter for enhancing Hartmann-Shack images and associated methods |
CN101806957A (zh) * | 2010-03-10 | 2010-08-18 | 中国科学院光电技术研究所 | 自适应光学系统中波前传感器与校正器对准装置 |
CN102221842A (zh) * | 2011-05-18 | 2011-10-19 | 中国科学院长春光学精密机械与物理研究所 | 千单元可扩展自适应光学系统波前处理器 |
CN102749143A (zh) * | 2012-07-10 | 2012-10-24 | 中国科学院光电技术研究所 | 一种提高夏克-哈特曼波前传感器测量精度的波前重构方法 |
CN105116705A (zh) * | 2015-09-29 | 2015-12-02 | 北京工业大学 | 多波长自适应数字全息成像系统及方法 |
CN105264428A (zh) * | 2013-06-06 | 2016-01-20 | 浜松光子学株式会社 | 自适应光学系统的调整方法、自适应光学系统和存储自适应光学系统用程序的记录介质 |
CN107144419A (zh) * | 2017-04-17 | 2017-09-08 | 中国科学院光电研究院 | 一种基于夏克‐哈特曼波前传感器的光学系统波像差测量装置与方法 |
CN109633885A (zh) * | 2019-01-29 | 2019-04-16 | 清华大学 | 基于波前聚焦的梯度折射率透镜显微成像系统及方法 |
CN111031244A (zh) * | 2019-12-17 | 2020-04-17 | 中国科学院光电技术研究所 | 一种基于fpga+gpu异构计算平台的波前处理控制方法 |
-
2021
- 2021-05-20 CN CN202110552221.XA patent/CN113295286A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020027640A1 (en) * | 2000-04-25 | 2002-03-07 | Campin John Alfred | Spatial filter for enhancing Hartmann-Shack images and associated methods |
CN101806957A (zh) * | 2010-03-10 | 2010-08-18 | 中国科学院光电技术研究所 | 自适应光学系统中波前传感器与校正器对准装置 |
CN102221842A (zh) * | 2011-05-18 | 2011-10-19 | 中国科学院长春光学精密机械与物理研究所 | 千单元可扩展自适应光学系统波前处理器 |
CN102749143A (zh) * | 2012-07-10 | 2012-10-24 | 中国科学院光电技术研究所 | 一种提高夏克-哈特曼波前传感器测量精度的波前重构方法 |
CN105264428A (zh) * | 2013-06-06 | 2016-01-20 | 浜松光子学株式会社 | 自适应光学系统的调整方法、自适应光学系统和存储自适应光学系统用程序的记录介质 |
CN105116705A (zh) * | 2015-09-29 | 2015-12-02 | 北京工业大学 | 多波长自适应数字全息成像系统及方法 |
CN107144419A (zh) * | 2017-04-17 | 2017-09-08 | 中国科学院光电研究院 | 一种基于夏克‐哈特曼波前传感器的光学系统波像差测量装置与方法 |
CN109633885A (zh) * | 2019-01-29 | 2019-04-16 | 清华大学 | 基于波前聚焦的梯度折射率透镜显微成像系统及方法 |
CN111031244A (zh) * | 2019-12-17 | 2020-04-17 | 中国科学院光电技术研究所 | 一种基于fpga+gpu异构计算平台的波前处理控制方法 |
Non-Patent Citations (2)
Title |
---|
卫培霞等: "基于GPU的快速夏克-哈特曼波前重构方法研究", 《长春理工大学学报(自然科学版)》 * |
钱思羽等: "优化夏克哈特曼波前传感器的光斑质心探测方法研究", 《长春理工大学学报(自然科学版)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110634140B (zh) | 一种基于机器视觉的大直径管状物定位及内壁缺陷检测方法 | |
CN112818988B (zh) | 一种指针式仪表自动识别读数方法及系统 | |
CN109489620B (zh) | 一种单目视觉测距方法 | |
CN109883533B (zh) | 基于机器视觉的低频振动测量方法 | |
WO2019105044A1 (zh) | 一种镜头畸变矫正和特征提取的方法及系统 | |
CN110246124B (zh) | 基于深度学习的目标尺寸测量方法及系统 | |
CN102589516B (zh) | 一种基于双目线扫描摄像机的动态距离测量系统 | |
CN110672020A (zh) | 一种基于单目视觉的立木高度测量方法 | |
CN103630091B (zh) | 一种基于激光与图像处理技术的叶面积测量方法 | |
CN103617611A (zh) | 一种自动阈值分割光斑中心及尺寸检测方法 | |
CN111307046B (zh) | 基于半球图像的树高测量方法 | |
CN112233076B (zh) | 基于红色圆标靶图像处理的结构振动位移测量方法及装置 | |
CN111179335A (zh) | 一种基于双目视觉的立木测定方法 | |
CN109308462B (zh) | 一种指静脉和指节纹感兴趣区域定位方法 | |
CN110617772A (zh) | 一种非接触式线径测量装置及方法 | |
CN111243003A (zh) | 车载双目摄像机及其检测道路限高杆的方法、装置 | |
CN114580559A (zh) | 一种基于单目视觉系统的测速方法 | |
CN115717867A (zh) | 一种基于机载双相机和目标追踪的桥梁变形测量方法 | |
CN117372498A (zh) | 一种基于三维点云的多姿态螺栓尺寸测量方法 | |
CN112330667A (zh) | 一种基于形态学的激光条纹中心线提取方法 | |
CN113295286A (zh) | 一种夏克-哈特曼波前传感器测量系统及相位恢复方法 | |
WO2023134251A1 (zh) | 一种基于聚类的光条提取方法及装置 | |
CN111243006A (zh) | 一种基于图像处理的液滴接触角及尺寸的测量方法 | |
CN114882095B (zh) | 一种基于轮廓匹配的物体高度在线测量方法 | |
CN115307865A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210824 |
|
RJ01 | Rejection of invention patent application after publication |