CN113096205B - 一种应用于光斑整形的二值图像的生成方法 - Google Patents
一种应用于光斑整形的二值图像的生成方法 Download PDFInfo
- Publication number
- CN113096205B CN113096205B CN202110254737.6A CN202110254737A CN113096205B CN 113096205 B CN113096205 B CN 113096205B CN 202110254737 A CN202110254737 A CN 202110254737A CN 113096205 B CN113096205 B CN 113096205B
- Authority
- CN
- China
- Prior art keywords
- shaping
- image
- light
- spot
- gray value
- 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
- 238000007493 shaping process Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 29
- 230000002238 attenuated effect Effects 0.000 claims 1
- 239000004973 liquid crystal related substance Substances 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 6
- 239000000463 material Substances 0.000 abstract description 5
- 230000000737 periodic effect Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/09—Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
- G02B27/0927—Systems for changing the beam intensity distribution, e.g. Gaussian to top-hat
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Facsimile Image Signal Circuits (AREA)
- Control Of Indicators Other Than Cathode Ray Tubes (AREA)
Abstract
本发明提供一种应用于光斑整形的二值图像的生成方法,生成的二值调制图像亮暗像素点随机分布且满足调制要求的空间分布比例,不带有过多细节,没有周期性明暗分布,不能对输入光斑产生明显的衍射作用,整形后输出光斑中没有明显衍射场强度分布,提升整形质量;同时,亮暗像素的随机分布对输入光斑具有局部匀光作用,可使输出光斑的光强分布更加均匀;最后,本发明所提供的技术同样适用于传统的利用液晶阵列、照相胶片、干版等可呈现灰度图的显示设备及材料作为显像介质对光斑进行振幅调制整形的***。
Description
技术领域
本发明属于光斑整形领域,尤其涉及一种应用于光斑整形的二值图像的生成方法。
背景技术
光斑整形是一种光学领域常用的技术,用于将光斑的原始振幅分布转换为需要的振幅分布。目前,实现光斑整形的技术有很多种,例如利用数字微镜阵列(DMD)对光斑进行振幅调制从而实现整形,然而利用数字微镜阵列、液晶阵列等显示设备及照相胶片、干版等显像材料呈现灰度图来调制原始光斑空间振幅分布的方式具有很强的自主性,理论上可以实现任意振幅分布光斑的整形。而数字微镜阵列、液晶阵列等显示设备可以动态显示调制图像,所以利用数字微镜阵列、液晶阵列并添加反馈还可实现光斑的实时动态整形。
上述调制技术的典型***如下图所示,上图为透射式***,主要应用于透射式显示器件(例如透射式液晶空间光调制器)或者显像材料(例如照相胶片)作为调制图像显示介质的情况。下图为反射式***,主要应用于反射式空间光振幅调制器件(例如硅上液晶器件LCoS、数字微镜阵列DMD)作为调制图像显示介质的情况,输入光斑经过准直***、扩束***后变为尺寸与调制图像尺寸相匹配的准平行光,之后投射到显示介质呈现的调制图像上,调制图像对光斑进行空间振幅调制,最终实现对输入光斑的整形。
上述调制图像分为两种:灰度图像和二值图像,以灰度级0-255为例,灰度图像每个像素点可呈现0-255共256个灰度级,而二值图像只有0和255两个灰度级,当利用图像对输入光斑进行调制时,不同的灰度级对光的衰减不同,灰度值越高说明此像素点对投射其上的光衰减越小,则经此像素点调制后输出的光强则越强,灰度值越低说明此像素点对投射其上的光衰减越大,则经此像素点调制后输出的光强则越弱,而DMD是通过时序脉冲式地转动微镜,使投射到显像方向地光通量形成时间积分上的差异从而呈现出灰度的效果,即当灰度值为0时,微镜的方向处于完全关闭状态,光无法投射到显像方向,当灰度值为255时,微镜的方向处于完全开启状态,光全部投射到显像方向,而当灰度介于0和255之间时,微镜的方向则不断地在关闭和开启之间转换,并通过控制开启时间占比来控制投射到显像方向的光通量从而呈现出需要的灰度。如上所述,灰度图像在应用于使用液晶阵列、照相底片等无机械振动的设备及材料作为显像介质进行光斑整形时更具优势,而对于使用DMD设备作为显像介质时,因为其呈现灰度值需要时序转动微镜方向,而这一过程存在机械转动,会使整形过程产生杂散光,同时会影响输入光斑的光束质量,而二值图像只有0和255两个灰度级,分别对应于DMD微镜的完全关闭和完全开启状态,所以用DMD呈现二值图像时不存在微镜的机械转动,不会引入杂散光,不影响光斑的光束质量,整形质量更好。
如上所述,原始光斑经准直扩束后投射到调制图像显示介质,调制图像显示介质上呈现的调制图像对光斑进行调制,从而改变光斑的光强空间分布进而实现对输入光斑的整形。显示介质呈现的调制图像为所述技术的核心部分,直接影响光斑整形的质量。因为二值图像只有0和255两个灰度级,只能通过0和255两个灰度值的空间分布比例来调制入射光斑的空间强度分布,常见的二值调制图像通常为使用Floyd-Stenberg方法、Hauck方法、Barnard方法等误差扩散算法将灰度图像转化而成的二值图像,因为上述误差算法最初是应用于图像处理的,所以转化而成的二值图像保留了很多原灰度图像的细节,同时会使生成的二值图像带有周期性明暗分布,当这类二值图像用于光斑调制时这些细节和明暗分布往往会对光形成衍射,使调制后的输出光斑带有衍射场分布,降低整形质量。
图3所示为使用Floyd-Stenberg方法将图2所示灰度图像转化而成的二值图像,明显可见其保留了原始灰度图像的太多细节且明暗分布周期性太强,这会使输出光斑带有明显的衍射场强度分布。
发明内容
为解决上述问题,本发明提供一种应用于光斑整形的二值图像的生成方法,生成的二值图像不包含过多细节,不会使整形后光斑带有衍射斑,同时可使整形后的光斑强度分布更加均匀,提升整形质量。
一种应用于光斑整形的二值图像的生成方法,包括以下步骤:
S1:假设调制图像包含M×N个像素,按照设定整形要求确定调制图像上各像素点的灰度值,得到大小为M×N的灰度值矩阵F1;
S2:随机生成一个大小为M×N的随机矩阵Frand,其中,随机矩阵Frand的每一个元素均为0~255之间的常数;
S3:分别比较灰度值矩阵F1与随机矩阵Frand相同位置上的元素值大小,并根据比较结果获取二值调制图像F2,其中,若F1(i,j)≥Frand(i,j),则F2(i,j)=255;若F1(i,j)<Frand(i,j),则F2(i,j)=0,其中,F1(i,j)为灰度值矩阵F1第i行第j列的灰度值,Frand(i,j)为随机矩阵Frand第i行第j列的元素值,F2(i,j)为二值调制图像F2第i行第j列的灰度值,且i=1,2,…,M,j=1,2,…,N。
进一步地,所述M=400,N=400。
进一步地,所述灰度值矩阵F1通过MATLAB按照设定整形要求计算得到。
进一步地,待整形的光斑为高斯光斑,且所述整形要求为将高斯光斑的光强自中心向边缘从强到弱进行衰减。
进一步地,所述灰度值矩阵F1的灰度值分布F1=0.5/A,其中,A表示高斯分布。
有益效果:
本发明提供一种应用于光斑整形的二值图像的生成方法,生成的二值调制图像亮暗像素点随机分布且满足调制要求的空间分布比例,不带有过多细节,没有周期性明暗分布,不能对输入光斑产生明显的衍射作用,整形后输出光斑中没有明显衍射场强度分布,提升整形质量;同时,亮暗像素的随机分布对输入光斑具有局部匀光作用,可使输出光斑的光强分布更加均匀;最后,本发明所提供的技术同样适用于传统的利用液晶阵列、照相胶片、干版等可呈现灰度图的显示设备及材料作为显像介质对光斑进行振幅调制整形的***。
附图说明
图1为光斑振幅调制整形***的结构示意图;
图2为灰度调制图像;
图3为使用Floyd-Stenberg方法将图2所示的灰度图像转化而成的二值图像的示意图;
图4为本发明提供的一种应用于光斑整形的二值图像的生成方法的流程图;
图5为采用本发明方法将图2所示的灰度图像转化而成的二值图像的示意图;
图6为高斯光斑光强分布;
图7为灰度值矩阵F1=0.5/A时对应的灰度图像示意图;
图8为本发明随机矩阵对应的随机分布灰度图像的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
实施例一
如图4所示,一种应用于光斑整形的二值图像的生成方法,包括以下步骤:
S1:假设调制图像包含M×N个像素,按照设定整形要求确定调制图像上各像素点的灰度值,得到大小为M×N的灰度值矩阵F1;
S2:随机生成一个大小为M×N的随机矩阵Frand,其中,随机矩阵Frand的每一个元素均为0~255之间的常数;
S3:分别比较灰度值矩阵F1与随机矩阵Frand相同位置上的元素值大小,并根据比较结果获取二值调制图像F2,其中,若F1(i,j)≥Frand(i,j),则F2(i,j)=255;若F1(i,j)<Frand(i,j),则F2(i,j)=0,其中,F1(i,j)为灰度值矩阵F1第i行第j列的灰度值,Frand(i,j)为随机矩阵Frand第i行第j列的元素值,F2(i,j)为二值调制图像F2第i行第j列的灰度值,且i=1,2,…,M,j=1,2,…,N。
也就是说,本发明提出的二值调制图像生成方法是在保证亮暗像素点空间分布比例符合调制要求的前提下实现亮暗像素的随机分布,其整体生成步骤为:第一步按照整形要求确定每个像素点的灰度值生成一幅灰度图像,假设调制图像包含M×N个像素,此处以一个M×N矩阵F1表示此步骤生成图像所有像素点的灰度值分布,同时生成一幅M×N大小且每个像素点灰度在0~255之间随机取值的图像。第二步对第一步生成的两幅图像的像素点灰度值逐个进行比较,同时生成一幅新的图像M×N大小的图像F2,此处分别用i和j表示像素点在整幅图像中的所在行位置和所在列位置,如果F1(i,j)≥Frand(i,j)则令F2(i,j)=255,否则令F2(i,j)=0,F2即为最终生成的二值调制图像,利用本发明提出的方法将图2所示灰度图像转化而成的二值图像如图5所示,可见其不包过多细节且没有周期性明暗分布。
由此可见,本发明生成的二值图像亮暗像素点随机分布且满足调制要求的空间分布比例,不带有过多细节,没有周期性明暗分布,不能对输入光斑产生明显的衍射作用,整形后输出光斑中没有明显衍射场强度分布,提升整形质量。除此之外,亮暗像素的随机分布对输入光斑具有局部匀光作用,可使输出光斑的光强分布更加均匀。
实施例二
下面以削减高斯光斑中心强度使其光斑更为平整的光斑整形实例来说明本发明提供的算法:
高斯光斑的强度分布如图6所示自中心向边缘逐渐减弱,为使其中心光强更为平整,需调制其振幅空间分布将其光强进行自中心向边缘从强到弱的衰减,从而实现整形的目的。
假设调制图像像素点个数为400×400,则用一个400×400的矩阵F1表示输入高斯光斑光强空间分布,第一步取调制图像的灰度值分布为F1=0.5/A,其中A表示高斯分布,则得到的图像如图7所示,同时生成一幅400×400大小且每个像素灰度值在0-255之间随机取值的灰度图像Frand,如图8所示。
第二步,对第一步生成的两幅图像的像素点灰度值逐个进行比较,同时生成一幅新的图像400×400大小的图像F2,如果F1(i,j)≥Frand(i,j)则令F2(i,j)=255,否则令F2(i,j)=0,F2即为最终生成的二值调制图像,如图5所示,至此利用本发明提出的算法生成调制图像完成。
最后利用图1所示***由最终调制图像对输入高斯光斑进行调制即可得到整形后中心光强平整的光斑。
由此可见,本发明所提算法生成的二值图像不包含过多细节,不会使整形后光斑带有衍射斑,同时可使整形后的光斑强度分布更加均匀,提升整形质量。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当然可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (3)
1.一种应用于光斑整形的二值图像的生成方法,其特征在于,包括以下步骤:
S1:假设调制图像包含M×N个像素,按照设定整形要求确定调制图像上各像素点的灰度值,得到大小为M×N的灰度值矩阵F1;其中,待整形的光斑为高斯光斑,且所述整形要求为将高斯光斑的光强自中心向边缘从强到弱进行衰减;所述灰度值矩阵F1的灰度值分布F1=0.5/A,其中,A表示高斯分布;
S2:随机生成一个大小为M×N的随机矩阵Frand,其中,随机矩阵Frand的每一个元素均为0~255之间的常数;
S3:分别比较灰度值矩阵F1与随机矩阵Frand相同位置上的元素值大小,并根据比较结果获取二值调制图像F2,其中,若F1(i,j)≥Frand(i,j),则F2(i,j)=255;若F1(i,j)<Frand(i,j),则F2(i,j)=0,其中,F1(i,j)为灰度值矩阵F1第i行第j列的灰度值,Frand(i,j)为随机矩阵Frand第i行第j列的元素值,F2(i,j)为二值调制图像F2第i行第j列的灰度值,且i=1,2,…,M,j=1,2,…,N。
2.如权利要求1所述的一种应用于光斑整形的二值图像的生成方法,其特征在于,所述M=400,N=400。
3.如权利要求1所述的一种应用于光斑整形的二值图像的生成方法,其特征在于,所述灰度值矩阵F1通过MATLAB按照设定整形要求计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110254737.6A CN113096205B (zh) | 2021-03-09 | 2021-03-09 | 一种应用于光斑整形的二值图像的生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110254737.6A CN113096205B (zh) | 2021-03-09 | 2021-03-09 | 一种应用于光斑整形的二值图像的生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113096205A CN113096205A (zh) | 2021-07-09 |
CN113096205B true CN113096205B (zh) | 2024-06-04 |
Family
ID=76666579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110254737.6A Active CN113096205B (zh) | 2021-03-09 | 2021-03-09 | 一种应用于光斑整形的二值图像的生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113096205B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4651287A (en) * | 1984-06-14 | 1987-03-17 | Tsao Sherman H | Digital image processing algorithm for output devices with discrete halftone gray scale capability |
DE19840769A1 (de) * | 1998-09-07 | 2000-03-09 | Holger Guenther Frenzel | Vorrichtung zur Verbesserung des Strahlproduktes und Kaskaden von mehreren dieser Vorrichtungen |
CN102141682A (zh) * | 2010-01-29 | 2011-08-03 | 中国科学院力学研究所 | 一种激光整形方法及整形后激光硬化处理设备及方法 |
CN102692714A (zh) * | 2012-05-31 | 2012-09-26 | 北京国科世纪激光技术有限公司 | 一种可调式光斑截取装置和方法 |
CN106127842A (zh) * | 2016-06-15 | 2016-11-16 | 北京工业大学 | 一种结合光源分布与反射特性的面曝光3d打印的方法及*** |
CN106228598A (zh) * | 2016-07-25 | 2016-12-14 | 北京工业大学 | 一种面向面曝光3d打印的模型自适应光照均匀化方法 |
CN107003785A (zh) * | 2014-12-09 | 2017-08-01 | 巴斯夫欧洲公司 | 光学检测器 |
CN107343200A (zh) * | 2017-07-06 | 2017-11-10 | 中南大学 | 一种基于光学方法的多彩色图像压缩及解压缩方法及装置 |
CN109785245A (zh) * | 2018-12-06 | 2019-05-21 | 江苏大学 | 一种光斑图像修整方法 |
CN111103744A (zh) * | 2018-10-26 | 2020-05-05 | 深圳光峰科技股份有限公司 | 显示设备与显示*** |
-
2021
- 2021-03-09 CN CN202110254737.6A patent/CN113096205B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4651287A (en) * | 1984-06-14 | 1987-03-17 | Tsao Sherman H | Digital image processing algorithm for output devices with discrete halftone gray scale capability |
DE19840769A1 (de) * | 1998-09-07 | 2000-03-09 | Holger Guenther Frenzel | Vorrichtung zur Verbesserung des Strahlproduktes und Kaskaden von mehreren dieser Vorrichtungen |
CN102141682A (zh) * | 2010-01-29 | 2011-08-03 | 中国科学院力学研究所 | 一种激光整形方法及整形后激光硬化处理设备及方法 |
CN102692714A (zh) * | 2012-05-31 | 2012-09-26 | 北京国科世纪激光技术有限公司 | 一种可调式光斑截取装置和方法 |
CN107003785A (zh) * | 2014-12-09 | 2017-08-01 | 巴斯夫欧洲公司 | 光学检测器 |
CN106127842A (zh) * | 2016-06-15 | 2016-11-16 | 北京工业大学 | 一种结合光源分布与反射特性的面曝光3d打印的方法及*** |
CN106228598A (zh) * | 2016-07-25 | 2016-12-14 | 北京工业大学 | 一种面向面曝光3d打印的模型自适应光照均匀化方法 |
CN107343200A (zh) * | 2017-07-06 | 2017-11-10 | 中南大学 | 一种基于光学方法的多彩色图像压缩及解压缩方法及装置 |
CN111103744A (zh) * | 2018-10-26 | 2020-05-05 | 深圳光峰科技股份有限公司 | 显示设备与显示*** |
CN109785245A (zh) * | 2018-12-06 | 2019-05-21 | 江苏大学 | 一种光斑图像修整方法 |
Non-Patent Citations (1)
Title |
---|
《激光光斑整形算法》;王锋 等;《红外与激光工程》;第43卷(第4期);3937-3940 * |
Also Published As
Publication number | Publication date |
---|---|
CN113096205A (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Amako et al. | Speckle-noise reduction on kinoform reconstruction using a phase-only spatial light modulator | |
Seldowitz et al. | Synthesis of digital holograms by direct binary search | |
Takaki et al. | Improvement of grayscale representation of the horizontally scanning holographic display | |
US8867113B2 (en) | Laser processing device and laser processing method | |
TW200923839A (en) | Apparatus and method for projection image | |
Wang et al. | Snapshot ptychography on array cameras | |
Qian et al. | Adaptive compressed 3D ghost imaging based on the variation of surface normals | |
Ye et al. | Robust phase retrieval with green noise binary masks | |
CN113096205B (zh) | 一种应用于光斑整形的二值图像的生成方法 | |
CN113325657B (zh) | 投影***、投影显示方法及终端设备 | |
Hu et al. | Quantitative phase imaging based on wavefront correction of a digital micromirror device | |
EP3817377B1 (en) | Device, system and method for modulating light using a phase light modulator and a spatial light modulator | |
Tudela et al. | Wavefront reconstruction by adding modulation capabilities of two liquid crystal devices | |
Luo et al. | Efficient glare suppression with Hadamard-encoding-algorithm-based wavefront shaping | |
Gene et al. | Defect corrections for coherent optical information processing of grayscale images in a DMD-based 4f-system using a collimated light source | |
Dashiell et al. | Nonlinear optical processing: analysis and synthesis | |
Melnychuck et al. | Fourier spectra of digital halftone images containing dot-position errors | |
CN112987320B (zh) | 一种应用于光斑振幅调制整形的调制图像生成方法 | |
Velez-Zea et al. | Alternative constraints for improved multiplane hologram generation | |
JP2022161038A (ja) | 液晶表示装置を使用するホログラム表示 | |
Shimobaba et al. | Comparison of double-phase hologram and binary amplitude encoding: holographic projection and vortex beam generation | |
CN107796594A (zh) | 一种基于空间光调制器的大气湍流模拟器 | |
Carlsen | Holographic page synthesis for sequential input of data | |
Ohana et al. | Linear optical random projections without holography | |
Cohn et al. | Ternary pseudorandom encoding of Fourier transform holograms |
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 |