CN113532320B - 一种基于图像的光斑衍射环分析方法、存储介质及芯片 - Google Patents

一种基于图像的光斑衍射环分析方法、存储介质及芯片 Download PDF

Info

Publication number
CN113532320B
CN113532320B CN202110820767.9A CN202110820767A CN113532320B CN 113532320 B CN113532320 B CN 113532320B CN 202110820767 A CN202110820767 A CN 202110820767A CN 113532320 B CN113532320 B CN 113532320B
Authority
CN
China
Prior art keywords
ring
image
diffraction
grid
spot
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
Application number
CN202110820767.9A
Other languages
English (en)
Other versions
CN113532320A (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.)
Hubei Optics Valley Laboratory
Wuhan Huagong Laser Engineering Co Ltd
Original Assignee
Hubei Optics Valley Laboratory
Wuhan Huagong Laser Engineering Co Ltd
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 Hubei Optics Valley Laboratory, Wuhan Huagong Laser Engineering Co Ltd filed Critical Hubei Optics Valley Laboratory
Priority to CN202110820767.9A priority Critical patent/CN113532320B/zh
Publication of CN113532320A publication Critical patent/CN113532320A/zh
Application granted granted Critical
Publication of CN113532320B publication Critical patent/CN113532320B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明涉及光学分析技术领域，为一种基于图像的光斑衍射环分析方法、存储介质及芯片，先获取具有衍射环的光斑图像，以最小步长为步长，根据中心光斑的中心点坐标、旋转角以及椭圆率，得到中心网格环，并通过每次增加一个步长重复上述步骤以提取得到光斑图像的所有的网格环；然后以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制得到环均值曲线；最后在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格以分离出各网格环以供分析。该方案仅基于图像分析即可实现光束衍射环的提取和测量，除一台相机外不需要额外的测量器件，即插即用，具有直观、简单、便携、高效的特点，并可以实现对具有不同圆度和角度光斑的测量。

Description

一种基于图像的光斑衍射环分析方法、存储介质及芯片
技术领域
本发明涉及光学分析技术领域，具体涉及一种基于图像的光斑衍射环分析方法、存储介质及芯片。
背景技术
可形成衍射环光斑的激光光束一般具有中心光斑并附带多级衍射环，其在激光加工工业的领域得以广泛应用。申请号为CN201510255809.3的发明专利公开了一种基于激光束衍射光斑特性的超精密车削加工表面三维微观形貌测量方法，其步骤如下：一、将激光器预热；二、将激光器、线性衰减片、小孔光阑、透镜、CCD相机调整在同一高度；三、将工件安装在回转工作台上，打开激光器，激光器输出的激光束依次经线性衰减片和小孔光阑后照射到工件表面，转动回转工作台以调整激光束的入射角度；四、激光束在工件表面发生衍射现象，产生衍射光斑，衍射光斑按照不同级次分布开，通过透镜后，被调整为平行光束；五、CCD相机采集衍射光斑图像，得到各级衍射光斑的强度和位置关系，再利用光栅方程计算出表面三维微观形貌的大小。
在对上述可形成衍射环光斑的光束进行应用时，需要对每个衍射环的环能量占比、尺寸大小等参数进行分析，从而了解具有不同衍射环的光束对加工材料的作用有何区别，使其更好的服务于生产。但目前尚无对单个衍射环的相关参数进行分析的方法和工具。
发明内容
本发明提供了一种基于图像的光斑衍射环分析方法、存储介质及芯片，解决了以上所述的目前尚无对单个衍射环的相关参数进行分析的方法和工具的技术问题。
本发明为解决上述技术问题提供了一种基于图像的光斑衍射环分析方法，包括以下步骤：
S1，获取具有衍射环的光斑图像，所述光斑图像包括中心光斑及由内而外依次环绕分布于所述中心光斑的***的多个衍射环，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，设定两轴方向的最小步长dx和dy；
S2，以最小步长为步长，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，得到中心网格环，并通过每次增加一个步长重复上述步骤以提取得到光斑图像的所有的网格环；
S3，以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制得到环均值曲线；
S4，在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格以分离出各网格环，提取出各网格环的数据及图像以供分析。
可选地，所述最小步长dx和dy的大小与椭圆率相等或成倍数关系。
可选地，多个衍射环包括第一衍射环P1、第二衍射环P2、第三衍射环P3...及第N衍射环PN；
根据S2方法分别获取各衍射环对应的网格环，其中所述第一衍射环P1、第二衍射环P2、第三衍射环P3...及第N衍射环PN一一对应的各椭圆图像为第一***椭圆S1、第二***椭圆S2、第三***椭圆S3...及第N***椭圆SN，N为大于或等于1的正整数。
可选地，所述S2具体包括：
S21，生成两个具有相同中心、旋转角及椭圆率的椭圆图像MASK_a和椭圆图像MASK_b，其中MASK_b的两轴分别比MASK_a大一个步长，椭圆图像的内部像素值填充1，外部像素值全部为0，将椭圆图像MASK_a和椭圆图像MASK_b进行异或得到中心网格环的mask图像；
S22，根据mask图像中每个像素灰度是否为0，提取输入图像相同坐标像素的灰度值，不为0则提取，为0则不提取，得到网格环上所有像素的数据，再对数据求平均，得到网格环的均值，同时，计算网格环的大小且满足
Figure GDA0004084542660000031
S23，以原有椭圆图像MASK_b作为新的椭圆图像MASK_a，并生成具有大一个步长的椭圆图像作为新的椭圆图像MASK_b，重复上述S21～S22，直到提取完所有的网格环图像。
可选地，所述S23之后还包括：
S24，直接将对应的mask图像与原图像相乘即可提取对应网格的图像信息。
可选地，各所述网格环的间距为固定定值。
可选地，所述阈值为给定的常量，为光斑图像的峰值乘以一个小于1的比例，所述比例取0.135.或者1/e^2，用来划定光斑图像的有效范围。
可选地，所述S4之后还包括S5：
在分析激光时，结合激光的加工效果，通过分析网格环的分布情况与加工质量的关系以帮助改进激光加工工艺；
在光学研究中，通过定量观察网格环变化以辅助调整光路。
本发明还提供了一种存储介质，所述存储介质用于存储计算机程序，所述计算机程序包括：用于执行基于图像的光斑衍射环分析方法的指令。
本发明还提供了一种芯片，包括：处理器，用于从存储器中调用并运行计算机程序，所述计算机程序包括：用于执行基于图像的光斑衍射环分析方法的指令。
有益效果：本发明提供了一种基于图像的光斑衍射环分析方法、存储介质及芯片，先获取具有衍射环的光斑图像，以最小步长为步长，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，得到中心网格环，并通过每次增加一个步长重复上述步骤以提取得到光斑图像的所有的网格环；然后以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制得到环均值曲线；最后在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格以分离出各网格环，提取出各网格环的数据及图像以供分析。通过采用具有设定步长的环形网格计算的方法，得到光斑图像的网格强度均值-环大小关系曲线。根据该曲线和设定的阈值来找出衍射环的边界，再从光斑图像中分离出各个衍射环，从而实现分析光斑各级衍射环的目的。
该方案仅基于图像分析即可实现光束衍射环的提取和测量，除一台相机外不需要额外的测量器件，即插即用，具有直观、简单、便携、高效的特点，并可以实现对具有不同圆度和角度光斑的测量。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。本发明的具体实施方式由以下实施例及其附图详细给出。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的流程示意图；
图2为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的光斑图像；
图3为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的衍射环示意图；
图4a为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的椭圆图像示意图；
图4b为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的网格环示意图；
图5为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的网格环均值原理图；
图6为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的衍射环分析结果图；
图7为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的环均值曲线图；
图8为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的网格环分离过程原理图；
图9为本发明基于图像的光斑衍射环分析方法、存储介质及芯片的分离后的衍射环图像。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。在下列段落中参照附图以举例方式更具体地描述本发明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
需要说明的是，当组件被称为“固定于”另一个组件，它可以直接在另一个组件上或者也可以存在居中的组件。当一个组件被认为是“连接”另一个组件，它可以是直接连接到另一个组件或者可能同时存在居中组件。当一个组件被认为是“设置于”另一个组件，它可以是直接设置在另一个组件上或者可能同时存在居中组件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
如图1至图9所示，本发明提供了一种基于图像的光斑衍射环分析方法，包括以下步骤：
S1，获取具有衍射环的光斑图像，所述光斑图像包括中心光斑及由内而外依次环绕分布于所述中心光斑的***的多个衍射环，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，设定两轴方向的最小步长dx和dy；这里的最小步长是通过自动计算时的参数获取方法获得，自动计算的参数包括中心点坐标、旋转角、椭圆率和步长。根据测试需要，也可人为设置上述的所有参数。
S2，以最小步长为步长，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，得到中心网格环，并通过每次增加一个步长重复上述步骤以提取得到光斑图像的所有的网格环；
S3，以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制得到环均值曲线；
S4，在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格以分离出各网格环，提取出各网格环的数据及图像以供分析。
具体地操作过程如下：
步骤一、输入光斑图像，如图2所示，通过手动输入或计算机自动进行光斑分析的方法，确定光斑的中心点坐标、旋转角和椭圆率。设定两轴方向的最小步长dx和dy，其比例应该与椭圆率相等。
步骤二、提取网格：
1、生成两个图像MASK_a和图像MASK_b，其中包含具有相同中心、旋转角的椭圆率，其中MASK_b的两轴分别比MASK_a大一个步长，椭圆的内部像素值填充1，外部像素值全部为0。将两个图像即图像MASK_a和图像MASK_b异或，得到网格环的mask图像。
2、根据mask图像中每个像素灰度是否为0，提取输入图像相同坐标像素的灰度值，不为0则提取，为0则不提取，得到网格环上所有像素的数据，再对数据求平均，得到网格环的均值。同时，计算网格环的大小且满足
Figure GDA0004084542660000081
此时计算的是第一个网格环，第N个网格环相应的乘以系数N即可，代表SN的大小，即均值曲线的自变量。
3、以原有MASK_b作为新的MASK_a，并生成具有大一个步长的椭圆的图像作为新的MASK_b，重复上述步骤，直到提取完所有的网格环图像。图3的示例中展示了所有的网格环图像。
4、如果要提取网格的图像，则直接将对应的mask图像与原图相乘即可。
步骤三、以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制成环均值曲线，如图7所示。
步骤四、在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格，如图6，然后以内外边界对应的各个网格的覆盖范围生成图像MASK_a和MASK_b，按照与步骤二中提取网格时相同的方法提取出各个衍射环的数据以及图像，其中分离的衍射环图像如图9所示。
步骤五、通过环上所有像素的数据进行分析计算。
在分析激光时，结合激光的加工效果，通过分析网格环的分布情况与加工质量的关系以帮助改进激光加工工艺；
在光学研究中，通过定量观察网格环变化以辅助调整光路。
在一个具体的实施场景中，包括如下步骤：
第一步，采用成像设备(如高清相机等)获取具有衍射环的光斑图像(如图2所示)，且所述光斑图像具有中心光斑O1以及环绕所述中心光斑O1由内而外依次分布的若干衍射环，其中，所述若干衍射环包括由内而外依次分布的第一衍射环P1、第二衍射环P2、第三衍射环P3...第N衍射环PN，N为大于或等于1的正整数；根据光束的衍射性质，所述中心光斑O的边界以及每一衍射环均为椭圆形，全部椭圆具有相同的中心。
第二步，如图3所示，获取所述光斑图像的中心光斑O1的中心、旋转角以及椭圆率，根据中心光斑O1的中心、旋转角以及椭圆率拟合出中心椭圆O2，且所述中心椭圆O2的长半轴为dx、短半轴为dy，且长半轴dx、短半轴dy要与椭圆率匹配；所述长半轴dx、短半轴dy可根据后续识别衍射环边界的精度需求而定，长半轴dx、短半轴dy越小，则后续识别衍射环边界的精度越高。
第三步，如图3所示，根据所述中心椭圆O2确定若干由内而外依次分布的***椭圆，即，第一***椭圆S1、第二***椭圆S2、第三***椭圆S3...第N***椭圆SN，所述中心椭圆O2以及每一***椭圆具有相同的中心，且由内而外，每一***椭圆的长半轴依次对应为M*dx、短半轴依次对应为M*dy，其中，M为大于或等于2的正整数，且M依次对应为2,3,4...，例如，由内而外，第一***椭圆S1的长半轴为2dx、短半轴为2dy，第二***椭圆S2的长半轴为3dx、短半轴为3dy...；其中，中心椭圆O2与第一***椭圆S1、相邻两个***椭圆的长半轴均相差dx、短半轴均相差dy。.
如上所述，因所述中心椭圆O2以及每一***椭圆为同心椭圆，因此，当确定中心椭圆O2的参数信息(包括中心、旋转角、椭圆率、长半轴dx、短半轴dy)后，通过增加长半轴dx、短半轴dy即可推算出每一***椭圆的位置。
第四步，如图4a-4b所示，根据所述中心椭圆O2生成中心椭圆图像A1，以及根据每一***椭圆生成对应的***椭圆图像，即根据第一***椭圆S1生成第一***椭圆图像B1、根据第二***椭圆S2生成第二***椭圆图像B2、根据第三***椭圆S3生成第三***椭圆图像B3...根据第N***椭圆SN生成第N***椭圆图像BN。
以及确定每相邻两个椭圆图像边界之间(即中心椭圆图像A1与第一***椭圆图像B1的边界之间、相邻两个***椭圆图像的边界之间)、且由内而外依次分布的若干网格环，所述若干网格环包括第一网格环L1、第二网格环L2...第N网格环LN，进一步的，根据公式(1)-(4)对应确定每一网格环的大小D1、D2...DN。
Figure GDA0004084542660000111
dσx(z)＝2*dx*N  (2)；
dσy(z)＝z*y*N  (3)；
DN＝dσ(Z)  (4)；
其中，dσ(Z)为与第N网格环外部边界对应的光束直径，其推导过程参见“ISO11145-2018光学和光子学-激光和激光相关设备-词汇和符号”，dσx(Z)为dσ(Z)在X方向的分量，dσy(Z)为dσ(Z)在y方向的分量；由此，获得第一网格环L1的大小
Figure GDA0004084542660000112
第二网格环L2的大小
Figure GDA0004084542660000113
第N网格环LN的大小
Figure GDA0004084542660000114
具体的，如图5所示，第四步具体包括：
S41、在所述中心椭圆O2内部填充第一像素值、外部填充第二像素值，以生成中心椭圆图像A1；以及在第一***椭圆S1内部填充第一像素值、外部填充第二像素值，以生成第一***椭圆图像B1；所述第一像素值为1，第二像素值为0。
S42、将中心椭圆图像A1、第一***椭圆图像B1异或，以获得两者之间的第一网格环L1。
S43、在第二***椭圆S2内部填充第一像素值、外部填充第二像素值，以生成第二***椭圆图像B2。
S44、将第一***椭圆图像B1、第二***椭圆图像B2异或，以获得两者之间的第二网格环L2。
S45、重复步骤S43-S44，以确定每相邻两个椭圆图像边界之间、且由内而外依次分布的若干网格环。
第五步，将当前网格环内的网格环像素点与光斑图像O进行对应，且从光斑图像O中提取出与网格环像素点对应的图像像素点，并计算所有图像像素点的像素平均值，以作为当前网格环所对应的光斑图像O的图像像素点环均值；依次重复，以获得每一网格环的所对应的光斑图像O的图像像素点环均值。
具体地，第五步具体包括：
S51、将一环网格图像通过函数转换为像素值组成的环网格图像二维数组，并通过两个for循环对环网格图像二维数组中的每个像素值进行索引，以得到包含该环网格图像内所有像素值的数列；
S52、将环形区域与原图点对应的实现过程，就是创建一个空数列，然后遍历原图所有的像素值。在遍历过程中索引某一个像素的同时，也索引二值图像L1中相同坐标的像素值。如果索引的L1的像素值等于1，则把原图中索引到的像素值添加到数列中去；如果等于0则跳过本像素，索引下一个，直到遍历完两个图的所有像素。(原图和L1的分辨率相同)。
例如，第一环形区域图像L1中像素值为1的时候，在原图中找到相同坐标的点(如图6所示)，将原图中所有这些点提取出来，然后求平均，作为环均值；
第六步，如图7所示，根据每一环形区域图像的大小dN以及每一环形区域图像的环均值建立曲线，其中，以每一环形区域图像的环均值为纵坐标、每一环形区域图像的大小dN为横坐标；且在所述曲线中确定衍射环边界判定阈值；
第七步，如图7至图9所示，根据所述判定阈值确定衍射环的边界，并据此提取出衍射环的图像；图7中的横向分布的阈值比例*峰值横向即为判定阈值曲线，将该图对应嵌入至图8进行对比分析。图8展示了从环网格到衍射环的对应关系，其分析过程是环均值曲线。竖向箭头表示网格环到衍射环边界的对应关系，即哪些网格环被最后判定为了衍射环的边界，它们正对应图7中的强度均值曲线与判定阈值曲线两条曲线相交的点。竖向箭头也刚好经过这些交点。根据衍射环边界提取环形区域即可得到衍射环的图像如图9所示。
第八步，对提取出的衍射环图像进行分析，以获取衍射环数据，所述衍射环数据包括环能量。
其中，通过mask图像像素值是否等于1来决定是否提取原图相同坐标的像素点，文中的L1至LN即各个mask图像。
具体工作过程原理为：通过计算机程序完成环均值的计算，计算内容包括图像中的通过mask图像标记的区域内的像素，计算所有这些像素的强度的均值。
标记方法：mask图像与原图分辨率相同，通过mask图像中像素强度是否为0，判断原图像素是否被标记。为0则认为原图相同坐标的像素为未标记的像素，不包括在计算平均值的像素点中；不为0则认为原图中同坐标的像素为标记的像素，包括在计算平均值的像素点中，即使强度为0也应当包括其中。。
本发明实施例还提供了一种存储介质，所述存储介质用于存储计算机程序，所述计算机程序包括：用于执行如前所述的基于图像的光斑衍射环分析方法的指令。
本发明还提供了一种芯片，包括：处理器，用于从存储器中调用并运行计算机程序，所述计算机程序包括：用于执行如前所述的基于图像的光斑衍射环分析方法的指令。
该方案仅基于图像分析即可实现光束衍射环的提取和测量，除一台相机外不需要额外的测量器件，即插即用，具有直观、简单、便携、高效的特点，并可以实现对具有不同圆度和角度光斑的测量。
本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (8)

1.一种基于图像的光斑衍射环分析方法，其特征在于，包括以下步骤：
S1，获取具有衍射环的光斑图像，所述光斑图像包括中心光斑及由内而外依次环绕分布于所述中心光斑的***的多个衍射环，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，设定两轴方向的最小步长dx和dy；
S2，以最小步长为步长，根据所述中心光斑的中心点坐标、旋转角以及椭圆率，得到中心网格环，并通过每次增加一个步长重复上述步骤以提取得到光斑图像的所有的网格环；具体包括：
S21，生成两个具有相同中心、旋转角及椭圆率的椭圆图像MASK_a和椭圆图像MASK_b，其中MASK_b的两轴分别比MASK_a大一个步长，椭圆图像的内部像素值填充1，外部像素值全部为0，将椭圆图像MASK_a和椭圆图像MASK_b进行异或得到中心网格环的mask图像；
S22，根据mask图像中每个像素灰度是否为0，提取输入图像相同坐标像素的灰度值，不为0则提取，为0则不提取，得到网格环上所有像素的数据，再对数据求平均，得到网格环的均值，同时，计算第N网格环LN的大小
Figure FDA0004084542650000011
N为大于或等于1的正整数；
S23，以原有椭圆图像MASK_b作为新的椭圆图像MASK_a，并生成具有大一个步长的椭圆图像作为新的椭圆图像MASK_b，重复上述S21～S22，直到提取完所有的网格环图像；
S3，以网格环均值为纵坐标、网格环大小为横坐标，将得到的网格环均值绘制得到环均值曲线；
S4，在环均值曲线中，根据阈值找出各个衍射环边界所对应的网格以分离出各网格环，提取出各网格环的数据及图像以供分析。
2.根据权利要求1所述的基于图像的光斑衍射环分析方法，其特征在于，所述最小步长dx和dy的大小与椭圆率相等或成倍数关系。
3.根据权利要求1所述的基于图像的光斑衍射环分析方法，其特征在于，多个衍射环包括第一衍射环P1、第二衍射环P2、第三衍射环P3...及第N衍射环PN；
根据S2方法分别获取各衍射环对应的网格环，其中所述第一衍射环P1、第二衍射环P2、第三衍射环P3...及第N衍射环PN一一对应的各椭圆图像为第一***椭圆S1、第二***椭圆S2、第三***椭圆S3...及第N***椭圆SN，N为大于或等于1的正整数。
4.根据权利要求1所述的基于图像的光斑衍射环分析方法，其特征在于，各所述网格环的间距为固定定值。
5.根据权利要求1所述的基于图像的光斑衍射环分析方法，其特征在于，所述阈值为给定的常量，为光斑图像的峰值乘以一个小于1的比例，所述比例取0.135或者1/e^2，用来划定光斑图像的有效范围。
6.根据权利要求1所述的基于图像的光斑衍射环分析方法，其特征在于，所述S4之后还包括S5：
在分析激光时，结合激光的加工效果，通过分析网格环的分布情况与加工质量的关系以帮助改进激光加工工艺；
在光学研究中，通过定量观察网格环变化以辅助调整光路。
7.一种存储介质，其特征在于：所述存储介质用于存储计算机程序，所述计算机程序包括：用于执行权利要求1至6任一项所述的基于图像的光斑衍射环分析方法的指令。
8.一种芯片，其特征在于，包括：处理器，用于从存储器中调用并运行计算机程序，所述计算机程序包括：用于执行权利要求1至6任一项所述的基于图像的光斑衍射环分析方法的指令。
CN202110820767.9A 2021-07-20 2021-07-20 一种基于图像的光斑衍射环分析方法、存储介质及芯片 Active CN113532320B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110820767.9A CN113532320B (zh) 2021-07-20 2021-07-20 一种基于图像的光斑衍射环分析方法、存储介质及芯片

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110820767.9A CN113532320B (zh) 2021-07-20 2021-07-20 一种基于图像的光斑衍射环分析方法、存储介质及芯片

Publications (2)

Publication Number Publication Date
CN113532320A CN113532320A (zh) 2021-10-22
CN113532320B true CN113532320B (zh) 2023-04-28

Family

ID=78100516

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110820767.9A Active CN113532320B (zh) 2021-07-20 2021-07-20 一种基于图像的光斑衍射环分析方法、存储介质及芯片

Country Status (1)

Country Link
CN (1) CN113532320B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114694147B (zh) * 2022-02-28 2024-04-19 华南理工大学 分割椭圆形图案中环绕文字的方法和装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006059452A (ja) * 2004-08-20 2006-03-02 Victor Co Of Japan Ltd 光ピックアップ装置
CN100442140C (zh) * 2005-12-01 2008-12-10 上海交通大学 投影仪生成平移面阵条纹实现累积三维成像的方法
CN103390277B (zh) * 2013-07-23 2016-07-20 济南凌空信息科技有限公司 图像分析方法
CN105954301B (zh) * 2016-06-22 2019-08-06 广东省智能制造研究所 一种基于机器视觉的瓶口质量检测方法
US10417495B1 (en) * 2016-08-08 2019-09-17 Google Llc Systems and methods for determining biometric information
CN107798330B (zh) * 2017-11-10 2021-07-20 上海电力学院 一种焊缝图像特征信息提取方法
CN108305320B (zh) * 2018-02-09 2021-06-04 重庆大学 用于提高大视野全息成像质量的自适应滑动窗重建方法
CN109523565B (zh) * 2018-11-15 2020-07-10 湖北工业大学 一种无衍射光莫尔条纹中心定位方法和***
CN111551567B (zh) * 2020-07-02 2021-06-25 厦门大学 一种基于条纹投影的物体表面缺陷检测方法及***
CN112116619B (zh) * 2020-09-16 2022-06-10 昆明理工大学 一种基于结构约束的多线结构光***条纹中心线提取方法

Also Published As

Publication number Publication date
CN113532320A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
Moru et al. A machine vision algorithm for quality control inspection of gears
WO2019119567A1 (en) Systems and methods for block based edgel detection with false edge elimination
CN111223053A (zh) 基于深度图像的数据增强方法
US20100066737A1 (en) Dynamic-state estimating apparatus, dynamic-state estimating method, and program
WO2022126870A1 (en) Three-dimensional imaging method and method based on light field camera and three-dimensional imaging measuring production line
CN113532320B (zh) 一种基于图像的光斑衍射环分析方法、存储介质及芯片
CN113012174B (zh) 一种图像融合方法、***及设备
CA3157444A1 (en) Method for depth estimation for a variable focus camera
CN109035170A (zh) 基于单网格图分段映射的自适应广角图像校正方法及装置
CN114651276A (zh) 制造方法
CN112329880A (zh) 一种基于相似性度量和几何特征的模板快速匹配方法
CN112465952A (zh) 一种基于光场相机微纳关联成像传感的重建方法
CN110223356A (zh) 一种基于能量生长的单目相机全自动标定方法
CN111325218B (zh) 基于光场图像的Hog特征检测与匹配方法
WO2013191132A1 (ja) 繰返パターン画像・三次元繰返構造データの解析システム及び方法
CN110675337A (zh) 一种图像降噪的扩散式遍历方法
CN115719414A (zh) 基于任意四边形回归的目标检测与精确定位方法
CN114283081A (zh) 基于金字塔加速的深度恢复方法、电子设备及存储介质
CN109801334B (zh) 工件定位方法、标准点确定方法、装置及设备
JP2018081378A (ja) 画像処理装置、撮像装置、画像処理方法および画像処理プログラム
CN113012132A (zh) 一种图像相似度确定方法、装置及计算设备、存储介质
CN111862098A (zh) 基于光场语义的个体匹配方法、装置、设备和介质
CN112287740B (zh) 基于YOLOv3-tiny实现的输电线路的目标检测方法、装置、无人机
CN114581448B (zh) 图像检测方法、装置、终端设备以及存储介质
Wu et al. Image stitching based on attention mechanism

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230327

Address after: 430000 Huagong technology precision micro nano intelligent manufacturing industrial park, No.66, future 2nd Road, Donghu New Technology Development Zone, Wuhan City, Hubei Province

Applicant after: WUHAN HUAGONG LASER ENGINEERING Co.,Ltd.

Applicant after: Hubei Optics Valley Laboratory

Address before: 430223 No. 66, future 2nd Road, Donghu New Technology Development Zone, Wuhan City, Hubei Province (Wuhan area of free trade zone)

Applicant before: WUHAN HUAGONG LASER ENGINEERING Co.,Ltd.

GR01 Patent grant
GR01 Patent grant