CN116597016A - 一种光纤内窥镜图像标定方法 - Google Patents

一种光纤内窥镜图像标定方法 Download PDF

Info

Publication number
CN116597016A
CN116597016A CN202310516632.2A CN202310516632A CN116597016A CN 116597016 A CN116597016 A CN 116597016A CN 202310516632 A CN202310516632 A CN 202310516632A CN 116597016 A CN116597016 A CN 116597016A
Authority
CN
China
Prior art keywords
image
calibration plate
function
calibration
extracting
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
Application number
CN202310516632.2A
Other languages
English (en)
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.)
General Engineering Research Institute China Academy of Engineering Physics
Original Assignee
General Engineering Research Institute China Academy of Engineering Physics
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 General Engineering Research Institute China Academy of Engineering Physics filed Critical General Engineering Research Institute China Academy of Engineering Physics
Priority to CN202310516632.2A priority Critical patent/CN116597016A/zh
Publication of CN116597016A publication Critical patent/CN116597016A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10068Endoscopic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种光纤内窥镜图像标定方法,包括对标定板图像进行高斯滤波;通过二值化滤掉经上步处理后的图像中的黑色背景和蜂窝网格,消除图像中无用的信息,提取标定板轮廓和标记点的特征;对标定板图像中的多个标记点进行识别和提取,对标定板里的标记点进行筛选,提取标记点区域,然后利用重心计算方式获得标记点中心坐标;根据提取出来的标记点中心坐标数组,采用点阵样板标定的方法对光纤内窥镜图像进行矫正。本发明鲁棒性强,该标定方法对不同距离、不同角度和混合多种噪声干扰下标定板的标记点识别有较好的鲁棒性;精度高,运用该标定方法的图像位移非接触式测试系统,测量误差不大于10μm,运用亚像素、分割等处理方法误差可达数个μm。

Description

一种光纤内窥镜图像标定方法
技术领域
本发明涉及图像标定技术领域,尤其涉及一种光纤内窥镜图像标定方法。
背景技术
光纤内窥镜因其直径小,探头可弯曲等特点被广泛应用于复杂结构狭小空间里的裂纹,腐蚀等情况的探查。探查时,将光纤内窥镜伸入被测物体内部,并在内窥镜目镜端安装一个高分辨率的数字相机,对光纤内窥镜里的影像进行拍摄存储,利用软件对采集的图像进行实时分析。
然而相机经光纤内窥镜采集到的图像受外界环境和自身缺陷的影响会有噪声产生,造成图像的降质,与常规镜头拍摄的图像相比会产生对比度明显降低、画质模糊、背景和目标物难以区分等现象,给后续的分析带来困难;如图1所示为常规镜头标定板图像;如图2所示为光纤镜头标定板图像;同时由于光纤镜的固有特性,为使拍摄到的靶面尽可能大,光线内窥镜的物镜采用了小孔径、大视场和大景深的设计原则,因此相机在目镜处采集到的图像会存在30%左右的桶形畸变,为了精确提取标记点中心坐标,需要对图像进行畸变校正。
基于上述两种原因,常规的图像标定计算方法对光纤内窥镜图像难以适用,会给分析带来较大的误差。
所以研发出一种光纤内窥镜图像标定方法来解决上述问题。
发明内容
本发明的目的就在于为了解决上述问题设计了一种光纤内窥镜图像标定方法。
本发明通过以下技术方案来实现上述目的:
一种光纤内窥镜图像标定方法,包括以下步骤:
S1、对相机通过光纤镜拍摄的标定板图像进行高斯滤波;
S2、通过二值化方法滤掉经步骤S1处理后的图像中的黑色背景和蜂窝网格,消除图像中无用的信息,提取标定板轮廓和标记点的特征;
S3、对经步骤S2处理后的标定板图像中的多个标记点进行识别和提取,对标定板里的标记点进行筛选,提取标记点区域,然后利用重心计算方式获得标记点中心坐标;
S4、根据第三步提取出来的标记点中心坐标数组,采用点阵样板标定的方法对光纤内窥镜图像进行矫正。
具体地,在步骤S1中,根据标定板图像中的蜂窝网格的噪声类型用gauss_filter()函数对图像进行高斯滤波。
具体地,在步骤S2中,采用局部阈值的二值化方法来提取标定板轮廓和标记点的特征。
进一步地,在步骤S2中,运用local_threshold()函数,对高斯滤波后的图像进行局部阈值的二值化操作,具体包括:
S21、通过调节local_threshold()函数中mask_size的值来设置掩码大小;
S22、通过设置local_threshold()函数中range的最大假设值,来对标记点的细节进一步提取;
S23、利用fill_up()函数对区域进行填充,将标定板填充成方形,再通select_shape()函数,根据标定板的实际长宽参数将标定板矩形区域提取出来,然后通过difference()函数,对提取的标定板矩形区域和原有区域做两次减法,将标定板轮廓和标记点提取出来;
S24、在提取出标定板轮廓和标记点后,运用gen_image_const()函数生成一张和原图同样大小的空白图像,再用paint_region()函数将提取的标定板轮廓和标记点打印到新图像中,至此提取出了标定板轮廓和标记点,过滤掉了蜂窝网格。
作为一种优选,在步骤S21中,将mask_size设置为大于要分割的标记点的当量直径,且为奇数。
作为另一种优选,在步骤S21中,mask_size为偶数,则默认使用下一个较大的奇数值。
优选地,在步骤S22中,用min_max_gray()函数获取图像的灰度值范围,range取GrayRange的0.5倍以内。
具体地,在步骤S4中,先通过calibrate Camera()函数得到图像的畸变系数,然后运用ReadParam()函数,读取畸变系数,当函数参数alpha为1时视场大小不变,小于1时需要缩放视场,再运用initUndistortRectifyMap()函数时根据畸变系数输出得到map1和map2两个参数,最后运用remap()函数获得矫正后图像。
本发明的有益效果在于:
1、通用性强,该标定方法可以针对不同直径、不同像元的光纤内窥镜拍摄的图像进行矫正识别;
2、鲁棒性强,该标定方法对不同距离、不同角度和混合多种噪声干扰下标定板的标记点识别有较好的鲁棒性;
3、精度高,运用该标定方法的图像位移非接触式测试系统,测量误差不大于10μm,运用亚像素、分割等处理方法误差可达数个μm。
附图说明
图1为常规镜头标定板图像;
图2为光纤镜头标定板图像;
图3-1为经高斯滤波前的标定板图像;
图3-2为经高斯滤波后的标定板图像;
图4为经自动阈值二值化后的标定板图像;
图5为经局部阈值二值化的标定板图像;
图6为经local_threshold()函数中不同的mask_size值处理的标定板图像;其中,a为mask_size=15;b为mask_size=31;
图7为经local_threshold()函数中不同的range最大假设值处理的标定板图像;其中,a为range=92.5;b为range=20.5;
图8为经过fill_up()函数处理后的标定板图像;
图9为经过select_shape()函数处理后的标定板图像;
图10为经过difference()函数第一次处理后的标定板图像;
图11为经过difference()函数第二次处理后的标定板图像;
图12为经过paint_region()函数处理后的标定板图像;
图13为经findCirclesGrid()函数中不同Minarea面积参数值处理的效果对比图,其中,a为Minarea=300;b为Minarea=100;
图14为标记点识别成功与提取的标定板图像;
图15为矫正过后的标定板图像;
图16为标记点图像区域提取示意图;
图17为激光位移计和图像位移测试系统的位移数据比较示意图;
图18为标定板距离镜头1cm,与镜头平行的标定矫正结果;其中,a为原图;b为标定矫正结果示意图;
图19为标定板距离镜头3cm,与镜头平行的标定矫正结果;其中,a为原图;b为标定矫正结果示意图;
图20为标定板距离镜头0.5cm,与镜头呈15°角度且有较多噪声干扰的标定矫正结果;其中,a为原图;b为标定矫正结果示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中,需要理解的是,术语“上”、“下”、“内”、“外”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,或者是本领域技术人员惯常理解的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,“设置”、“连接”等术语应做广义理解,例如,“连接”可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接连接,也可以通过中间媒介间接连接,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
下面结合附图,对本发明的具体实施方式进行详细说明。
如图1所示,一种光纤内窥镜图像标定方法,包括以下步骤:
S1、对相机通过光纤镜拍摄的标定板图像进行高斯滤波;光纤镜头由无数根光纤束缠绕聚合而成,相机通过光纤镜拍摄的图像会呈现出蜂窝网格状,对标记点的识别产生的干扰,因此需要对图像中的蜂窝网格进行滤波。根据该蜂窝网格的噪声类型用gauss_filterlocal_threshold()函数对图像进行高斯滤波,滤波效果如图3-1和3-2所示,图3-1为经高斯滤波前的标定板图像;图3-2为经高斯滤波后的标定板图像;
S2、通过二值化方法滤掉经步骤S1处理后的图像中的黑色背景和蜂窝网格,消除图像中无用的信息,增强感兴趣信息的可检测性,为标定板上的标记点和背景的分离提供可靠的保证;采用局部阈值的二值化方法来提取标定板轮廓和标记点的特征;
常用的二值化方法有自动阈值和局部阈值两种方法,自动阈值二值化auto_threshold()函数可以显示标定板、标记点的细节,但是对标记点和标定板轮廓的有效特征却无法提取,如图4所示。因此本标定方法采用局部阈值的二值化方法来提取标定板轮廓和标记点的特征,即运用local_threshold()函数,对高斯滤波后的图像进行局部阈值的二值化操作,如图5所示。运用local_threshold()函数,对高斯滤波后的图像进行局部阈值的二值化操作具体包括:
S21、通过调节local_threshold()函数中的mask_size来设置掩码大小;即局部阈值二值化方法中邻域的大小。mask_size值越小,标记点分割的越细,故必须将mask_size设置为大于要分割的标记点的当量直径,一般为奇数,如果mask_size为偶数,则默认使用下一个较大的奇数值。如图6所示,为经不同local_threshold()函数中的mask_size值处理的标定板图像;其中,a为mask_size=15;b为mask_size=31;
S22、通过设置local_threshold()函数中range的最大假设值,来对标记点的细节进一步提取;用min_max_gray()函数获取图像的灰度值范围,range取GrayRange的0.5倍以内,range值越大,提取点的越细小,range值越小,点越粗。根据经验range值取range=GrayRange/RangeDev,其中RangeDev为系数,取值9。如图7所示,为不同range值情况下的标定板图像;其中,a为range=92.5;b为range=20.5;
S23、利用fill_up()函数对区域进行填充,将标定板填充成方形,如图8所示;再通select_shape()函数,根据标定板的实际长宽参数将标定板矩形区域提取出来,如图9所示;然后通过difference()函数,对提取的标定板矩形区域和原有区域做两次减法,将标定板轮廓和标记点提取出来,两次处理分别如图10、11所示;为经过fill_up()处理后的标定板图像;
S24、在提取出标定板轮廓和标记点后,运用gen_image_const()函数生成一张和原图同样大小的空白图像,再用paint_region()函数将提取的标定板轮廓和标记点打印到新图像中,至此成功提取出了标定板轮廓和标记点,过滤掉了蜂窝网格,如图12所示。
S3、运用findCirclesGrid()函数对经步骤S2处理后的标定板图像中的49个标记点进行识别和提取,提取标记点区域,然后利用重心计算方式获得标记点中心坐标;提取应用如图13所示。
当最小面积参数Minarea不合适时,标定板标记点中心坐标无法正确提取;需要根据标记点所占像素面积大小进行设置,只有设置合适的最小面积参数时,标定板标记点才可以正确的识别出来,并且按一定顺序排列。如图14为标记点识别成功与提取的标定板图像。
S4、根据第三步提取出来的标记点中心坐标数组,采用点阵样板标定的方法对光纤内窥镜图像进行矫正。
先通过calibrate Camera()函数得到图像的畸变系数,然后运用ReadParam()函数,读取畸变系数,当函数参数alpha为1时视场大小不变,小于1时需要缩放视场,再运用initUndistortRectifyMap()函数时根据畸变系数输出得到map1和map2两个参数,最后运用remap()函数获得矫正后图像;该标定方法可消除桶形畸变95%左右,误差1个像素,根据标定关系与采用的亚像素计算方法,实际约为5μm,如图15所示为矫正过后的图像。
一般标记点区域只占整张图片的一小部分,因此可以设置图像感兴趣区域(ROI)来框选需要的标记点,只对该ROI里的图像进行处理计算。使用draw_rectangle1()函数在图像中手动画出感兴趣的区域后,再使用gen_rectangle1()函数生成一个矩形区域,最后使用reduce_domain()函数获得感兴趣的图像。
经过区域提取获得感兴趣图像后,识别该区域标记点的横纵坐标,求取位移变化。利用精密位移平台制造微位移,分析计算位移变化前后标记点的形心坐标,根据已知的位移值可以得出像素和位移的对应关系。然后采用激光位移计与图像位移测试系统在精密位移平台上进行比对验证,验证方法:通过光纤内窥镜进行标记点图像采集,设计的标记点形状、布局如图16所示,标记点当量直径0.5mm,横向间距1mm,纵向间距2mm,后背胶粘贴,运用该标定方法对采集的标记点图像进行计算分析,将图像位移测试系统得出的位移值与激光位移计测出的值进行比对。
激光位移计的测量误差为1μm,数据如表1所示(对应地,如图17所示,为激光位移计和图像位移测试系统的位移数据比较示意图),其中T代表时间,利用位移平台制造产生8次微位移D1-D8,C代表像素坐标系下标记点形心坐标变化值,D代表平均像素对应位移值,JD代表激光位移计测量的数据,TD代表图像位移测试系统数据。根据标定关系与采用的亚像素计算方法,从表中数据可以分析得出1像素对应位移约10μm,且激光位移计测量数据和图像位移测试系统数据两者基本重合,最大差值2.4μm,说明该图像标定方法和图像位移测试系统可行,测试结果真实可靠。
表1激光位移计和图像位移测试系统的位移数据
运用设计的标定系统对光纤内窥镜拍摄的不同类型的标定板图像进行标定矫正,图18为标定板距离镜头1cm,与镜头平行的标定矫正结果;图19为标定板距离镜头3cm,与镜头平行的标定矫正结果;图20为标定板距离镜头0.5cm,与镜头呈一定角度且有较多噪声干扰的标定矫正结果。
由以上图例可以看出,该标定方法的通用性和鲁棒性较强,适合多种场景下、不同质量图像的矫正识别。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种光纤内窥镜图像标定方法,其特征在于,包括以下步骤:
S1、对相机通过光纤镜拍摄的标定板图像进行高斯滤波;
S2、通过二值化方法滤掉经步骤S1处理后的图像中的黑色背景和蜂窝网格,消除图像中无用的信息,提取标定板轮廓和标记点的特征;
S3、对经步骤S2处理后的标定板图像中的多个标记点进行识别和提取,对标定板里的标记点进行筛选,提取标记点区域,然后利用重心计算方式获得标记点中心坐标;
S4、根据第三步提取出来的标记点中心坐标数组,采用点阵样板标定的方法对光纤内窥镜图像进行矫正。
2.根据权利要求1所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S1中,根据标定板图像中的蜂窝网格的噪声类型用gauss_filter()函数对图像进行高斯滤波。
3.根据权利要求1所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S2中,采用局部阈值的二值化方法来提取标定板轮廓和标记点的特征。
4.根据权利要求3所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S2中,运用local_threshold()函数,对高斯滤波后的图像进行局部阈值的二值化操作,具体包括:
S21、通过调节local_threshold()函数中mask_size的值来设置掩码大小;
S22、通过设置local_threshold()函数中range的最大假设值,来对标记点的细节进一步提取;
S23、利用fill_up()函数对区域进行填充,将标定板填充成方形,再通select_shape()函数,根据标定板的实际长宽参数将标定板矩形区域提取出来,然后通过difference()函数,对提取的标定板矩形区域和原有区域做两次减法,将标定板轮廓和标记点提取出来;
S24、在提取出标定板轮廓和标记点后,运用gen_image_const()函数生成一张和原图同样大小的空白图像,再用paint_region()函数将提取的标定板轮廓和标记点打印到新图像中,至此提取出了标定板轮廓和标记点,过滤掉了蜂窝网格。
5.根据权利要求4所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S21中,将mask_size设置为大于要分割的标记点的当量直径,且为奇数。
6.根据权利要求4所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S21中,mask_size为偶数,则默认使用下一个较大的奇数值。
7.根据权利要求4所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S22中,用min_max_gray()函数获取图像的灰度值范围,range取GrayRange的0.5倍以内。
8.根据权利要求1所述的一种光纤内窥镜图像标定方法,其特征在于,在步骤S4中,先通过calibrate Camera()函数得到图像的畸变系数,然后运用ReadParam函数,读取畸变系数,当函数参数alpha为1时视场大小不变,小于1时需要缩放视场,再运用initUndistortRectifyMap()函数时根据畸变系数输出得到map1和map2两个参数,最后运用remap()函数获得矫正后图像。
CN202310516632.2A 2023-05-09 2023-05-09 一种光纤内窥镜图像标定方法 Pending CN116597016A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310516632.2A CN116597016A (zh) 2023-05-09 2023-05-09 一种光纤内窥镜图像标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310516632.2A CN116597016A (zh) 2023-05-09 2023-05-09 一种光纤内窥镜图像标定方法

Publications (1)

Publication Number Publication Date
CN116597016A true CN116597016A (zh) 2023-08-15

Family

ID=87592997

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310516632.2A Pending CN116597016A (zh) 2023-05-09 2023-05-09 一种光纤内窥镜图像标定方法

Country Status (1)

Country Link
CN (1) CN116597016A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117710250A (zh) * 2024-02-04 2024-03-15 江苏无右微创医疗科技有限公司 一种消除纤维镜成像蜂窝状结构的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117710250A (zh) * 2024-02-04 2024-03-15 江苏无右微创医疗科技有限公司 一种消除纤维镜成像蜂窝状结构的方法
CN117710250B (zh) * 2024-02-04 2024-04-30 江苏无右微创医疗科技有限公司 一种消除纤维镜成像蜂窝状结构的方法

Similar Documents

Publication Publication Date Title
CN110866924B (zh) 一种线结构光中心线提取方法及存储介质
CN112308832B (zh) 一种基于机器视觉的轴承质量检测方法
CN110706224B (zh) 基于暗场图像的光学元件弱划痕检测方法、系统和装置
CN111260615B (zh) 基于激光和机器视觉融合的无人机桥梁表观病害检测方法
CN115170669B (zh) 基于边缘特征点集配准的识别定位方法及系统、存储介质
CN112686920A (zh) 一种圆形零件几何尺寸参数视觉测量方法及系统
CN115880373B (zh) 基于新型编码特征的立体视觉系统的标定板及标定方法
CN113034474A (zh) 一种oled显示器晶圆图的测试方法
CN112200790B (zh) 布料缺陷检测方法、设备和介质
CN116597016A (zh) 一种光纤内窥镜图像标定方法
CN111354047A (zh) 一种基于计算机视觉的摄像模组定位方法及系统
CN110705553B (zh) 一种适用于车辆远景图像的划痕检测方法
CN112836726B (zh) 一种基于视频信息的指针类仪表示数读取方法及装置
CN113112396B (zh) 导电粒子的检测方法
CN116563298B (zh) 基于高斯拟合的十字线中心亚像素检测方法
CN116883987A (zh) 一种用于变电站无人巡检的指针仪表读数识别方法
CN116977316A (zh) 一种复杂外形构件损伤缺陷全场检测与量化评估方法
CN113554688B (zh) 一种基于单目视觉的o型密封圈尺寸测量方法
CN116125489A (zh) 一种室内物体三维检测方法、计算机设备及存储介质
CN113012121B (zh) 裸片扫描结果的处理方法、装置、电子设备与存储介质
CN114387232A (zh) 一种晶圆圆心定位、晶圆缺口定位及晶圆定位校准方法
CN114821035A (zh) 一种电力设备红外测温设备距离参数识别方法
CN110660073B (zh) 一种直线边缘识别设备
CN117115488B (zh) 一种基于图像处理的水表检测方法
CN117593302B (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