CN113960778A - 一种基于中频滤波的动态步长聚焦方法 - Google Patents
一种基于中频滤波的动态步长聚焦方法 Download PDFInfo
- Publication number
- CN113960778A CN113960778A CN202111153799.4A CN202111153799A CN113960778A CN 113960778 A CN113960778 A CN 113960778A CN 202111153799 A CN202111153799 A CN 202111153799A CN 113960778 A CN113960778 A CN 113960778A
- Authority
- CN
- China
- Prior art keywords
- search
- image
- current
- focusing
- definition 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000001914 filtration Methods 0.000 title claims abstract description 18
- 238000011156 evaluation Methods 0.000 claims abstract description 15
- 238000010845 search algorithm Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 description 13
- 210000004027 cell Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000000349 chromosome Anatomy 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000031864 metaphase Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/24—Base structure
- G02B21/241—Devices for focusing
- G02B21/244—Devices for focusing using image analysis techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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/10—Image acquisition modality
- G06T2207/10056—Microscopic image
- G06T2207/10061—Microscopic image from scanning electron microscope
-
- 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/20048—Transform domain processing
- G06T2207/20052—Discrete cosine transform [DCT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Optics & Photonics (AREA)
- Quality & Reliability (AREA)
- Automatic Focus Adjustment (AREA)
Abstract
本发明提出了一种基于中频滤波的动态步长聚焦方法,主要是分为两部分,第一部分是聚焦评价函数,在进行聚焦评价前先对图像做中频滤波的预处理,滤除图像的高频和低频信息,再使用拉普拉斯算子计算图像的清晰度值;第二部分是基于第一部分计算清晰度值的方法,找到最佳聚焦平面的动态步长搜索算法,在每轮搜索的过程中记录搜索区域,搜索步长,并根据上一轮的搜索区域和搜索步长动态更新本轮搜索的初始搜索步长和搜索区域,有效降低了局部极值对搜索结果的影响,可以有效的找到清晰度值最大的聚焦平面。
Description
技术领域
本发明涉及岩石薄片显微图像的自动聚焦方法,尤其涉及一种基于中频滤波的动态步长聚焦方法。
背景技术
在石油地质勘探中,通过将钻井取出来的岩石制作成岩石薄片,再对岩石薄片进行分析就可以得出石油的含量情况。这一整个岩石薄片图像采集过程中,对于岩石薄片的聚焦清晰度的评判直接影响后续的图像拼接和图像分析的结果。而手动聚焦耗费大量的人力物力,因此,对于岩石薄片图像采集使用自动聚焦方法。自动聚焦又分为主动聚焦和被动聚焦两种方式,主动聚焦是通过光学仪器测量镜头和待测物体之间的距离,来调整二者的距离达到聚焦的目的;被动聚焦则是通过图像处理算法来评价当前焦平面的清晰度,移动载物台获取清晰度评价最好的焦平面,主动聚焦对光学仪器的精度要求较高,实现起来有一定难度,所以被动聚焦的应用范围更加广泛。
对于不同的应用场景,自动聚焦有着不同的技术难题。在被动式的岩石薄片自动聚焦问题主要是清晰度评价函数和聚焦搜索两个方面,对此,国内外学者都提出了各自的解决方案。梁隆恺针对采集的图片分辨率高,计算量大,根据图像灰度特点提出了灰度非零值统计函数和低灰度值统计函数的聚焦算法。He C对低倍率显微图像中的细胞定位,提出了基于八邻域算子和最小二乘平面拟合的自动聚焦算法。还有Qiu Y等人针对高分辨率下的中期染色体聚焦算法,Mo C等人针对视频监控和跟踪系统的自动对焦算法。另外,对于聚焦区域的选择,中外学者都针对相应的应用场景提出不同的聚焦区域选择算法,其中,由于细胞薄片的感兴趣区域分布不均匀,Zhang F等人与江旻珊等人分别就人工鱼群算法进行改进,以获取合适的聚焦窗口。Tsai T H等人分别就运动检测获取感兴趣区域,在对这个区域选择合适的焦窗,在连续图像中对该窗口进行跟踪,Jeon J等人就图像的熵和边缘信息选择最优的点扩散函数获取聚焦区域。为高效的获取聚焦平面,江旻珊提出了混合搜索法来提高搜索效率。
在岩石薄片图像的采集中,由于岩石薄片中不同矿物本身的岩性不同以及岩石薄片的制作工艺等问题,容易引进一些噪声。传统的聚焦评价函数若直接用于对岩石显微图像的评价,对评价结果的准确性往往不理想,因此需要进行改进。
因此需要研发出一种基于中频滤波的动态步长聚焦方法来解决上述问题。
发明内容
本发明的目的就在于为了解决上述问题设计了一种基于中频滤波的动态步长聚焦方法。
本发明通过以下技术方案来实现上述目的:
一种基于中频滤波的动态步长聚焦方法,包括以下步骤:
S1、将当前视场的图像做中频离散余弦变换,即将图像与中频离散余弦变换算子做卷积,滤除图像中的高频和低频信息,保留图像中的中频信息;
S2、对S1中的处理结果的图像使用Laplacian算子计算图像的清晰度值;
S3、根据S2中的清晰度值,进行动态步长搜索找到最佳聚焦平面。
具体地,步骤S2中对S1中的中频信息图像使用Laplacian算子计算图像梯度,再根据计算出的梯度之和作为图像清晰度值。
具体地,步骤S3中使用动态步长搜索算法,找到最佳聚焦平面,动态步长搜索算法如下:
S31、先选定一个初始步长L为起始搜索步长;
S32、每次跨L步长向后搜索,并且每次到位置后计算清晰度评价值F[i],i表示当前位置,同时判断F[i]是否是单调递增的,若是,则在本轮搜索完成后用M记录本次搜索的清晰度值最大的位置,若出现F[i]<F[i-1],则继续向后搜索K次,K为F[i]的位置与本次搜索起始位置之差的一半,在后的K步搜索中若得到的清晰度值为递减的,则终止此次搜索前往S33步,若不是则按规则更新K值后继续搜索;
S33、根据本次搜索的最大清晰度值的位置M与搜索起始位置、搜索终止位置计算搜索区域D,D为M与搜索起始位置之差和M与搜索结束位置之差较小者,同时用diff表示本次搜索得到的清晰度最大值与上一次搜索得到的清晰度最大值之差,用end表示结束搜索标志,end为本次搜索的清晰度最大值乘以终止搜索系数t,若diff小于end则终止搜索,本次得到的最大清晰度值的位置即为最佳聚焦平面,若diff大于等于end,则转S34;
S34、由S33中得到的D更新搜索区域,跟新为[M-D/2,M+D/2],M为上次搜索得到的最大聚焦评价值的位置,同时更新新一轮搜索的步长L,L为上次搜索步长*本次搜索区域D/上次搜索区域D,转S32;
算法结束即找到最佳聚焦平面。
本发明的有益效果在于:
本发明提出的基于中频滤波的动态步长聚焦方法,就岩石薄片的特性先将图像用中频滤波算子(中频离散余弦变换)预处理,滤除高频和低频信息,保留图像的基本结构,再根据Laplacian图像梯度算子计算图像聚焦评价值。在计算出聚焦评价值的基础上使用动态步长搜索算法记录每次的搜索结果,并通过这些搜索结果动态更新搜索区域和搜索步长,可以更准确的找到最佳聚焦平面。
附图说明
图1是基于中频滤波的动态步长聚焦方法的流程图;
图2是无局部极值的两阶段搜索法;
图3是无局部极值的动态步长搜索法;
图4是有局部极值的两阶段搜索法;
图5是有局部极值的动态步长搜索法;
图6是动态步长聚焦搜索算法流程图。
具体实施方式
下面结合附图对本发明作进一步说明:
如图1所示,一种基于中频滤波的动态步长聚焦方法,包括以下步骤:
S1、将当前视场的图像做中频离散余弦变换,即将图像与中频离散余弦变换算子做卷积,滤除图像中的高频和低频信息,保留图像中的中频信息;
S2、对S1中的处理结果的图像使用Laplacian算子计算图像的清晰度值;步骤S2中对S1中的中频信息图像使用Laplacian算子计算图像梯度,再根据计算出的梯度之和作为图像清晰度值。
S3、根据S2中的清晰度值,进行动态步长搜索找到最佳聚焦平面。
如图6所示,步骤S3中使用动态步长搜索算法,找到最佳聚焦平面,动态步长搜索算法如下:
S31、先选定一个初始步长L为起始搜索步长;
S32、每次跨L步长向后搜索,并且每次到位置后计算清晰度评价值F[i],i表示当前位置,同时判断F[i]是否是单调递增的,若是,则在本轮搜索完成后用M记录本次搜索的清晰度值最大的位置,若出现F[i]<F[i-1],则继续向后搜索K次,K为F[i]的位置与本次搜索起始位置之差的一半,在后的K步搜索中若得到的清晰度值为递减的,则终止此次搜索前往S33步,若不是则按规则更新K值后继续搜索;
S33、根据本次搜索的最大清晰度值的位置M与搜索起始位置、搜索终止位置计算搜索区域D,D为M与搜索起始位置之差和M与搜索结束位置之差较小者,同时用diff表示本次搜索得到的清晰度最大值与上一次搜索得到的清晰度最大值之差,用end表示结束搜索标志,end为本次搜索的清晰度最大值乘以终止搜索系数t,若diff小于end则终止搜索,本次得到的最大清晰度值的位置即为最佳聚焦平面,若diff大于等于end,则转S34;
S34、由S33中得到的D更新搜索区域,跟新为[M-D/2,M+D/2],M为上次搜索得到的最大聚焦评价值的位置,同时更新新一轮搜索的步长L,L为上次搜索步长*本次搜索区域D/上次搜索区域D,转S32;
算法结束即找到最佳聚焦平面。
为验证本发明的可行性,选取了两组图像样本,一组是无局部极值点的图像,另一组是有局部极值点的图像,从图2和图3的结果可以看出两阶段搜索法对大小步长的设置比较依赖,找到的结果也不是最佳聚焦平面,而本发明的动态步长可以动态调整搜索步长,最终找到最佳聚焦平面;从图4和图5来看,对于存在局部极值点的图像样本,两阶段搜索法容易受到局部极值点的影响,错把局部极值点当成最佳聚焦平面,而本发明的动态步长搜索由于可以动态的调整步长和调整搜索区域,可以有效避免局部极值的影响,找到最终的最佳聚焦平面。
本发明的技术方案不限于上述具体实施例的限制,凡是根据本发明的技术方案做出的技术变形,均落入本发明的保护范围之内。
Claims (3)
1.一种基于中频滤波的动态步长聚焦方法,其特征在于,包括以下步骤:
S1、将当前视场的图像做中频离散余弦变换,即将图像与中频离散余弦变换算子做卷积,滤除图像中的高频和低频信息,保留图像中的中频信息;
S2、对S1中的处理结果的图像使用Laplacian算子计算图像的清晰度值;
S3、根据S2中的清晰度值,进行动态步长搜索找到最佳聚焦平面。
2.根据权利要求1所述的基于中频滤波的动态步长聚焦方法,其特征在于,步骤S2中对S1中的中频信息图像使用Laplacian算子计算图像梯度,再根据计算出的梯度之和作为图像清晰度值。
3.根据权利要求1所述的基于中频滤波的动态步长聚焦方法,其特征在于,步骤S3中使用动态步长搜索算法,找到最佳聚焦平面,动态步长搜索算法如下:
S31、先选定一个初始步长L为起始搜索步长;
S32、每次跨L步长向后搜索,并且每次到位置后计算清晰度评价值F[i],i表示当前位置,同时判断F[i]是否是单调递增的,若是,则在本轮搜索完成后用M记录本次搜索的清晰度值最大的位置,若出现F[i]<F[i-1],则继续向后搜索K次,K为F[i]的位置与本次搜索起始位置之差的一半,在后的K步搜索中若得到的清晰度值为递减的,则终止此次搜索前往S33步,若不是则按规则更新K值后继续搜索;
S33、根据本次搜索的最大清晰度值的位置M与搜索起始位置、搜索终止位置计算搜索区域D,D为M与搜索起始位置之差和M与搜索结束位置之差较小者,同时用diff表示本次搜索得到的清晰度最大值与上一次搜索得到的清晰度最大值之差,用end表示结束搜索标志,end为本次搜索的清晰度最大值乘以终止搜索系数t,若diff小于end则终止搜索,本次得到的最大清晰度值的位置即为最佳聚焦平面,若diff大于等于end,则转S34;
S34、由S33中得到的D更新搜索区域,跟新为[M-D/2,M+D/2],M为上次搜索得到的最大聚焦评价值的位置,同时更新新一轮搜索的步长L,L为上次搜索步长*本次搜索区域D/上次搜索区域D,转S32;
算法结束即找到最佳聚焦平面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111153799.4A CN113960778A (zh) | 2021-09-29 | 2021-09-29 | 一种基于中频滤波的动态步长聚焦方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111153799.4A CN113960778A (zh) | 2021-09-29 | 2021-09-29 | 一种基于中频滤波的动态步长聚焦方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113960778A true CN113960778A (zh) | 2022-01-21 |
Family
ID=79463198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111153799.4A Pending CN113960778A (zh) | 2021-09-29 | 2021-09-29 | 一种基于中频滤波的动态步长聚焦方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113960778A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012058665A (ja) * | 2010-09-13 | 2012-03-22 | Sony Corp | 顕微鏡制御装置及び処理範囲決定方法 |
EP2448246A1 (en) * | 2010-10-28 | 2012-05-02 | Axis AB | Method for focusing |
CN104459940A (zh) * | 2013-09-25 | 2015-03-25 | 北京环境特性研究所 | 一种快速自适应自动聚焦方法 |
CN105354817A (zh) * | 2015-09-25 | 2016-02-24 | 济南中维世纪科技有限公司 | 一种噪声图像自动聚焦方法 |
CN109785323A (zh) * | 2019-01-25 | 2019-05-21 | 淮阴师范学院 | 一种基于中频滤波的图像聚焦测度实现方法 |
CN110166692A (zh) * | 2019-05-31 | 2019-08-23 | 济南中维世纪科技有限公司 | 一种提高摄像机自动聚焦准确率和速度的方法及装置 |
WO2021148465A1 (en) * | 2020-01-22 | 2021-07-29 | Intuity Media Lab GmbH | Method for outputting a focused image through a microscope |
-
2021
- 2021-09-29 CN CN202111153799.4A patent/CN113960778A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012058665A (ja) * | 2010-09-13 | 2012-03-22 | Sony Corp | 顕微鏡制御装置及び処理範囲決定方法 |
EP2448246A1 (en) * | 2010-10-28 | 2012-05-02 | Axis AB | Method for focusing |
CN104459940A (zh) * | 2013-09-25 | 2015-03-25 | 北京环境特性研究所 | 一种快速自适应自动聚焦方法 |
CN105354817A (zh) * | 2015-09-25 | 2016-02-24 | 济南中维世纪科技有限公司 | 一种噪声图像自动聚焦方法 |
CN109785323A (zh) * | 2019-01-25 | 2019-05-21 | 淮阴师范学院 | 一种基于中频滤波的图像聚焦测度实现方法 |
CN110166692A (zh) * | 2019-05-31 | 2019-08-23 | 济南中维世纪科技有限公司 | 一种提高摄像机自动聚焦准确率和速度的方法及装置 |
WO2021148465A1 (en) * | 2020-01-22 | 2021-07-29 | Intuity Media Lab GmbH | Method for outputting a focused image through a microscope |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102053339B (zh) | 自动聚焦方法 | |
CN105578048B (zh) | 一种快速对焦方法和装置、移动终端 | |
US8553137B2 (en) | Image data processing method and apparatus | |
AU2011357735B2 (en) | Fast auto-focus in microscopic imaging | |
CN107664899B (zh) | 自动对焦的方法、装置以及系统 | |
CN112099217B (zh) | 一种显微镜自动对焦方法 | |
CN110531484B (zh) | 一种对焦过程模型可设定的显微镜自动对焦方法 | |
CN110824689B (zh) | 全自动显微图像景深拓展系统及其方法 | |
CN102062929A (zh) | 一种显微镜系统的自动聚焦方法和装置 | |
CN109361849B (zh) | 一种自动对焦的方法 | |
CN112019751B (zh) | 基于标定信息的自动聚焦方法 | |
CN110324536B (zh) | 一种用于显微相机的图像变化自动感知对焦方法 | |
CN107911599B (zh) | 一种红外图像全局自动对焦方法及装置 | |
CN110363734B (zh) | 厚样本显微荧光图像重构方法及系统 | |
CN110830710A (zh) | 一种超长焦镜头的聚焦方法 | |
CN113705298A (zh) | 图像采集方法、装置、计算机设备及存储介质 | |
CN116594168A (zh) | 一种基于可变焦透镜的快速大景深显微成像系统及方法 | |
CN113538545B (zh) | 一种基于电液可调焦镜头的单目深度估计方法及相应的相机和存储介质 | |
CN109218624B (zh) | 一种光电跟踪系统的温度调焦补偿方法 | |
CN113960778A (zh) | 一种基于中频滤波的动态步长聚焦方法 | |
CN117671033A (zh) | 一种基于夜间灯光追踪的相机像主点快速定标方法和系统 | |
CN117170081A (zh) | 一种基于图像清晰度识别的自动对焦方法 | |
CN116405775A (zh) | 自动聚焦方法和装置、电子设备及拍摄系统 | |
CN112689853A (zh) | 图像处理方法、装置、拍摄设备、可移动平台和存储介质 | |
CN108431660B (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 |