CN113160040A - 用于医疗内窥镜成像的管腔图像展开方法 - Google Patents

用于医疗内窥镜成像的管腔图像展开方法 Download PDF

Info

Publication number
CN113160040A
CN113160040A CN202110493564.3A CN202110493564A CN113160040A CN 113160040 A CN113160040 A CN 113160040A CN 202110493564 A CN202110493564 A CN 202110493564A CN 113160040 A CN113160040 A CN 113160040A
Authority
CN
China
Prior art keywords
image
blank
lumen
pixel points
rectangular image
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.)
Granted
Application number
CN202110493564.3A
Other languages
English (en)
Other versions
CN113160040B (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.)
Alton Shanghai Medical Instruments Co ltd
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN202110493564.3A priority Critical patent/CN113160040B/zh
Publication of CN113160040A publication Critical patent/CN113160040A/zh
Application granted granted Critical
Publication of CN113160040B publication Critical patent/CN113160040B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G06T3/073

Abstract

本发明涉及一种用于医疗内窥镜成像的管腔图像展开方法,该方法利用内窥镜前端的摄像头采集生物体管腔图像,并由管腔图像展开算法对所采集的图像进行处理,将环状管腔图像展开为矩形平铺图像,管腔图像展开算法的步骤包括:(1)以内窥镜图像的中心为圆心,获取内窥镜图像的最大内接圆;(2)建立空白矩形图像,以最大内接圆周长为空白矩形图像的长边,以最大内接圆半径的1/2作为空白矩形图像的短边;(3)计算空白矩形图像上每一个像素点对应于原图像的亮度值,填入空白矩形图像,最终形成完整的平铺图像。该方法可将生物体管腔内壁图像转化成连续的管腔内壁平铺图像,方便医生观察管腔内壁情况,降低了手术难度,提高了手术成功率。

Description

用于医疗内窥镜成像的管腔图像展开方法
技术领域
本发明涉及医疗内窥镜图像展开领域,具体是指一种用于医疗内窥镜成像的管腔图像展开方法。
背景技术
随着科技的发展,医用电子内窥镜已成为当前应用非常广泛的一种医疗仪器,医生通过电子内窥镜不仅能直接观察到人体内脏器官的组织形态及病变情况,而且还能对内窥镜图像做进一步处理,以达到更好的视觉和诊断效果。但由于内窥镜安装部位受限,医生只能观察某一时刻内窥镜所在管腔的环状内壁图像,无法观察管腔内壁平铺图像,而且镜头拍摄的图像中可能存在病灶点变形等一系列影响医生诊断的问题。因此,需要一种用于医疗内窥镜成像的管腔图像展开方法,为医生提供良好的内窥镜诊断视觉效果,降低医生误诊的可能性。
发明内容
本发明的目的是克服了上述现有技术的缺点,提供了一种用于医疗内窥镜成像的管腔图像展开方法,该方法可将所述环状内壁图像展开为所述平铺图像,方便医生观察管腔内壁的情况。
为了实现上述目的,本发明的用于医疗内窥镜成像的管腔图像展开方法如下:
该用于医疗内窥镜成像的管腔图像展开方法,其主要特点是,所述的方法包括以下步骤:
(1)以内窥镜图像的中心为圆心,获取内窥镜图像的最大内接圆;
(2)建立空白矩形图像,以最大内接圆周长为空白矩形图像的长边,以最大内接圆半径的1/2作为空白矩形图像的短边;
(3)计算空白矩形图像上每一个像素点对应于原图像的亮度值,填入空白矩形图像,最终形成完整的平铺图像。
较佳的,所述内窥镜图像的中心是内窥镜镜头光轴的中心。
较佳的,所述的步骤(2)的建立空白矩形图像的步骤具体包括以下处理过程:
(2-1)从圆心到最大内接圆之间包含若干个同心圆,选取半径为最大内接圆半径1/2的同心圆到最大内接圆之间的区域进行展开,从剖开线将同心圆区域展开为平铺矩形图像;
(2-2)建立空白矩形图像,以最大内接圆周长为展开平铺矩形图像的长边,以最大内接圆半径的1/2作为展开平铺矩形图像的短边;
较佳的,所述的步骤(3)的计算空白矩形图像上每一个像素点亮度值的步骤具体包括以下处理过程:
(3-1)计算获得空白矩形图像上每一个像素点对应于原图像上的坐标;
(3-2)利用原图中上述坐标点周围四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值;
较佳的,所述的步骤(3-1)具体包括以下处理过程:
(3-1.1)空白矩形图像的长对应原图像最大内接圆的一周,可计算出空白矩形每一列像素点的区域宽度对应的弧度;
(3-1.2)空白矩形图像第一行像素点对应于原图上同心圆的半径为Rmax,由Rmax依次递减可得到每一行像素点对应于原图上同心圆的半径,最后一行对应的半径为Rmax/2;
(3-1.3)由空白矩形图像上每一个像素点所对应的半径和弧度,计算所述像素点对应于原图像上的位置坐标。
较佳的,所述的步骤(3-2)的利用坐标点周围的四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值的步骤具体包括以下处理过程:
(3-2.1)判断横纵坐标是否为整数,如果是,则对应原图像上的像素点,取原图像素点的亮度值填入空白图像对应的像素点;否则,继续步骤(3-2.2);
(3-2.2)通过该坐标点周围的四个像素点的亮度值进行权重计算。
较佳的,所述的步骤(3-2.2)的通过该坐标点周围的四个像素点的亮度值进行权重计算,具体为:
根据以下公式计算亮度值:
Gij=(1-α)×(1-β)×GA+α×(1-β)×GB+(1-α)×β×GC+α×β×GD
其中,Gij为空白矩形图像第i行第j列的像素点的亮度值,GA、GB、GC、GD分别为坐标点周围的四个像素点A、B、C、D的亮度值,α和β为根据坐标点离周围四个像素点的距离而确定的权重系数。
采用了本发明的用于医疗内窥镜成像的管腔图像展开方法,可以将所述环状管腔图像展开为平铺图像,呈现出所述连续的管腔内壁平铺图。该发明可方便医生观察管腔内壁的情况。本发明实时性好,精度高,稳定性好,有效降低手术难度,提高手术成功率。
附图说明
图1为本发明的用于医疗内窥镜成像的管腔图像展开方法的示意图。
图2~图5为本发明的用于医疗内窥镜成像的管腔图像展开方法的空白矩形图像中像素点坐标位于不同区域时的示意图。
图6为本发明的用于医疗内窥镜成像的管腔图像展开方法的空白矩形图像中像素点亮度值计算方法的示意图。
具体实施方式
为了能够更清楚地描述本发明的技术内容,下面结合具体实施例来进行进一步的描述。
该用于医疗内窥镜成像的管腔图像展开方法,其主要特点是,所述的方法包括以下步骤:
(1)以内窥镜图像的中心为圆心,获取内窥镜图像的最大内接圆;
(2)建立空白矩形图像,以最大内接圆周长为空白矩形图像的长边,以最大内接圆半径的1/2作为空白矩形图像的短边;
(3)计算空白矩形图像上每一个像素点对应于原图像的亮度值,填入空白矩形图像,最终形成完整的平铺图像。
作为本发明的优选实施方式,所述内窥镜图像的中心是内窥镜镜头光轴的中心。
作为本发明的优选实施方式,所述的步骤(2)的建立空白矩形图像的步骤具体包括以下处理过程:
(2-1)从圆心到最大内接圆之间包含若干个同心圆,选取半径为最大内接圆半径1/2的同心圆到最大内接圆之间的区域进行展开,从剖开线将同心圆区域展开为平铺矩形图像;
(2-2)建立空白矩形图像,以最大内接圆周长为展开平铺矩形图像的长边,以最大内接圆半径的1/2作为展开平铺矩形图像的短边;
作为本发明的优选实施方式,所述的步骤(3)的计算空白矩形图像上每一个像素点亮度值的步骤具体包括以下处理过程:
(3-1)计算获得空白矩形图像上每一个像素点对应于原图像上的坐标;
(3-2)利用原图中上述坐标点周围四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值;
作为本发明的优选实施方式,所述的步骤(3-1)具体包括以下处理过程:
(3-1.1)空白矩形图像的长对应原图像最大内接圆的一周,可计算出空白矩形每一列像素点的区域宽度对应的弧度;
(3-1.2)空白矩形图像第一行像素点对应于原图上同心圆的半径为Rmax,由Rmax依次递减可得到每一行像素点对应于原图上同心圆的半径,最后一行对应的半径为Rmax/2;
(3-1.3)由空白矩形图像上每一个像素点所对应的半径和弧度,可计算出所述像素点对应于原图像上的位置坐标。
作为本发明的优选实施方式,所述的步骤(3-2)的利用坐标点周围的四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值的步骤具体包括以下处理过程:
(3-2.1)判断横纵坐标是否为整数,如果是,则对应原图像上的像素点,取原图像素点的亮度值插入空白图像对应的像素点;否则,继续步骤(3-2.2);
(3-2.2)通过该坐标点周围的四个像素点的亮度值进行权重计算。
作为本发明的优选实施方式,所述的步骤(3-2.2)的通过该坐标点周围的四个像素点的亮度值进行权重计算,具体为:
根据以下公式计算亮度值:
Gij=(1-α)×(1-β)×GA+α×(1-β)×GB+(1-α)×β×GC+α×β×GD
其中,Gij为空白矩形图像第i行第j列的像素点的亮度值,GA、GB、GC、GD分别为坐标点周围的四个像素点A、B、C、D的亮度值,α和β为根据坐标点离周围四个像素点的距离而确定的权重系数。
在本发明的一具体实施方式中,本技术方案具体通过以下方式予以实现:
(S1)以内窥镜图像的中心为圆心,获取图像的最大内接圆;
(S2)从圆心到最大内接圆之间包含若干个同心圆,选取半径为最大内接圆半径1/2的同心圆到最大内接圆之间的区域进行展开,从剖开线将同心圆区域展开为平铺矩形图像;
(S3)建立空白矩形图像,以最大内接圆周长为展开平铺矩形图像的长边,以最大内接
圆半径的1/2作为展开平铺矩形图像的短边;
(S4)计算空白矩形图像上每一个像素点对应于原图像上的坐标,具体步骤如下:
(S4-1)以原图图像中心为原点建立坐标系,以y轴负方向为0度,每个像素点的坐标都是整数,原点处坐标为(0,0);
(S4-2)空白矩形图像的长对应原图像最大内接圆的一周,可计算出空白矩形每一列像素点的区域宽度对应的弧度;
进一步地,所述弧度计算公式如下:
Figure BDA0003053378260000041
Rmax为同心圆区域最大半径。
(S4-3)空白矩形图像每一行像素点对应于同心圆的一层,第一行像素点对应于原图上同心圆的半径为Rmax,由Rmax依次递减可得到每一行像素点对应于原图上同心圆的半径,最后一行对应的半径为Rmax/2;
进一步地,设空白图像第i行元素对应的同心圆半径为Ri,则R1=Rmax,R2=Rmax-1,Ri=Rmax-i+1,依此类推。
(S4-4)由空白矩形图像上每一个像素点所对应的半径和弧度,可计算出所述像素点对应于原图像上的位置坐标。
进一步地,将原图像从以y轴负方向起顺时针将图像分为4个区域,
Figure BDA0003053378260000051
为I区,
Figure BDA0003053378260000052
为II区,
Figure BDA0003053378260000053
为III区,
Figure BDA0003053378260000054
为IV区。
当计算的坐标点位于I区时,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000055
Figure BDA0003053378260000056
Figure BDA00030533782600000516
Figure BDA00030533782600000517
以y轴负方向为0度算起时,
Figure BDA0003053378260000057
为正,
Figure BDA0003053378260000058
为正,而坐标点对应的横坐标为负,纵坐标为负,因而要取负号。权重系数α、β计算为正。
当计算的坐标点位于II区时,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000059
Figure BDA00030533782600000510
Figure BDA00030533782600000518
Figure BDA00030533782600000519
以y轴负方向为0度算起时,
Figure BDA00030533782600000511
为正,
Figure BDA00030533782600000512
为负,而坐标点对应的横坐标为负,纵坐标为正,因而要取负号。权重系数α、β计算为正。
当计算的坐标点位于III区时,坐标点和权重系数的计算公式具体为:
Figure BDA00030533782600000513
Figure BDA00030533782600000514
Figure BDA00030533782600000515
Figure BDA00030533782600000616
以y轴负方向为0度算起时,
Figure BDA0003053378260000061
为负,
Figure BDA0003053378260000062
为负,而坐标点对应的横坐标为正,纵坐标为正,因而要取负号。权重系数α、β计算为正。
当计算的坐标点位于IV区时,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000063
Figure BDA0003053378260000064
为负;
Figure BDA0003053378260000065
Figure BDA0003053378260000066
为正;
Figure BDA0003053378260000067
Figure BDA0003053378260000068
以y轴负方向为0度算起时,
Figure BDA0003053378260000069
为负,
Figure BDA00030533782600000610
为正,而坐标点对应的横坐标为正,纵坐标为负,因而要取负号。权重系数α、β计算为正。
故无论在哪个区域,坐标点和权重系数的计算公式都为:
Figure BDA00030533782600000611
Figure BDA00030533782600000612
Figure BDA00030533782600000613
Figure BDA00030533782600000614
其中,Xij、Yij分别为空白矩形图像上第i行第j列的像素点对应于原图像上的位置坐标,
Figure BDA00030533782600000615
表示从原点指向坐标点的向量与y轴负方向的夹角,α、β为上述坐标点于像素点A的距离。
(S5)利用坐标点周围的四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值,具体包括以下步骤:
(S5-1)当计算出来的横纵坐标Xij、Yij为整数时,正好可对应原图像上的像素点,取原图像素点的亮度值插入空白图像对应的像素点即可;
(S5-2)当计算出来的横纵坐标Xij、Yij不为整数时,利用该坐标点周围的的四个像素点的亮度值进行权重计算,像素点离得越近权重越大。
(S5-3)不妨设上述坐标点位于坐标分别为(P,Q)、(P,Q+1)、(P+1,Q)、(P+1,Q+1)的四个像素点A、B、C、D之间(P,Q都为整数),A、B、C、D的亮度值分别为GA、GB、GC、GD
(S5-4)上述坐标点亮度值的权重系数计算方法具体如下:
Figure BDA0003053378260000072
Figure BDA0003053378260000073
其中,Xij,Yij为上述坐标点的横纵坐标,
Figure BDA0003053378260000074
分别为横纵坐标向下取整。α、β为上述坐标点于像素点A的距离,由于A、B、C、D四点构成了一个边长为单位l的正方形,四个像素点与坐标点之间的距离都可用α、β来表示。像素点距离越近影响越大,A、B、C、D四个像素点对坐标点亮度值的影响权重系数也可用α、β来表示。
(S5-5)通过上述权重系数,计算出空白矩形图像上第i行第j列的亮度值;
进一步地,上述算法具体如下:
Gij=(1-α)×(1-β)×GA+α×(1-β)×GB+(1-α)×β×GC+α×β×GD
其中,Gij为空白矩形图像上第i行第j列的亮度值,GA、GB、GC、GD为所计算坐标(Xij,Yij)周围四个像素点的亮度值。
在本发明的一具体实施方式中,请参阅图1所示,该用于医疗内窥镜成像的管腔图像展开方法,其中所述的内窥镜图像展开为平铺矩形图像的方法具体包括以下步骤:
(S-1)以内窥镜图像的中心为圆心,获取图像的最大内接圆;
(S-2)从圆心到最大内接圆之间包含若干个同心圆,选取半径为最大内接圆半径1/2的同心圆到最大内接圆之间的区域进行展开,从剖开线将同心圆区域展开为平铺矩形图像;
(S-3)建立空白矩形图像,以最大内接圆周长为展开平铺矩形图像的长边,以最大内接圆半径的1/2作为展开平铺矩形图像的短边;
请参阅图2~图5所示,该用于医疗内窥镜成像的管腔图像展开方法,其中所述的空白矩形图像上每一个像素点对应于原图像上的坐标算法具体包括以下步骤:
(S4-1)以原图图像中心为原点建立坐标系,以y轴负方向为0度,每个像素点的坐标都是整数,原点处坐标为(0,0);
(S4-2)空白矩形图像的长对应原图像最大内接圆的一周,可计算出空白矩形每一列像素点的区域宽度对应的弧度;
进一步地,所述弧度计算公式如下:
Figure BDA0003053378260000071
Rmax为同心圆区域最大半径。
(S4-3)空白矩形图像每一行像素点对应于同心圆的一层,第一行像素点对应于原图上同心圆的半径为Rmax,由Rmax依次递减可得到每一行像素点对应于原图上同心圆的半径,最后一行对应的半径为Rmax/2;
进一步地,设空白图像第i行元素对应的同心圆半径为Ri,则R1=Rmax,R2=Rmax-1,Ri=Rmax-i+1,依此类推。
(S4-4)由空白矩形图像上每一个像素点所对应的半径和弧度,可计算出所述像素点对应于原图像上的位置坐标。
进一步地,将原图像从以y轴负方向起顺时针将图像分为4个区域,
Figure BDA0003053378260000081
为I区,
Figure BDA0003053378260000082
为II区,
Figure BDA0003053378260000083
为III区,
Figure BDA0003053378260000084
为IV区。
当计算的坐标点位于I区时,如图2,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000085
Figure BDA0003053378260000086
Figure BDA00030533782600000821
Figure BDA00030533782600000822
以y轴负方向为0度算起时,
Figure BDA0003053378260000087
为正数,
Figure BDA0003053378260000088
为正数,而坐标点对应的横坐标为负数,纵坐标为负数,因而要取负号。权重系数α、β计算为正数。
当计算的坐标点位于II区时,如图3,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000089
Figure BDA00030533782600000810
Figure BDA00030533782600000811
Figure BDA00030533782600000812
以y轴负方向为0度算起时,
Figure BDA00030533782600000813
为正数,
Figure BDA00030533782600000814
为负数,而坐标点对应的横坐标为负数,纵坐标为正数,因而要取负号。权重系数α、β计算为正数。
当计算的坐标点位于III区时,如图4,坐标点和权重系数的计算公式具体为:
Figure BDA00030533782600000815
Figure BDA00030533782600000816
Figure BDA00030533782600000817
Figure BDA00030533782600000818
以y轴负方向为0度算起时,
Figure BDA00030533782600000819
为负数,
Figure BDA00030533782600000820
为负数,而坐标点对应的横坐标为正数,纵坐标为正数,因而要取负号。权重系数α、β计算为正数。
当计算的坐标点位于IV区时,如图5,坐标点和权重系数的计算公式具体为:
Figure BDA0003053378260000091
为负数;
Figure BDA0003053378260000092
为正数;
Figure BDA0003053378260000093
Figure BDA0003053378260000094
以y轴负方向为0度算起时,
Figure BDA0003053378260000095
为负数,
Figure BDA0003053378260000096
为正数,而坐标点对应的横坐标为正数,纵坐标为负数,因而要取负号。权重系数α、β计算为正数。
故无论在哪个区域,坐标点和权重系数的计算公式都为:
Figure BDA0003053378260000097
Figure BDA0003053378260000098
Figure BDA0003053378260000099
Figure BDA00030533782600000910
其中,Xij、Yij分别为空白矩形图像上第i行第j列的像素点对应于原图像上的位置坐标,
Figure BDA00030533782600000911
表示从原点指向坐标点的向量与y轴负方向的夹角,α、β为上述坐标点于像素点A的距离。
请参阅图6所示,所述的利用坐标点周围的四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值具体包括以下步骤:
(S5-1)当计算出来的横纵坐标Xij、Yij为整数时,正好可对应原图像上的像素点,取原图像素点的亮度值插入空白图像对应的像素点即可;
(S5-2)当计算出来的横纵坐标Xij、Yij不为整数时,利用该坐标点周围的的四个像素点的亮度值进行权重计算,像素点离得越近权重越大。
(S5-3)不妨设上述坐标点位于坐标分别为(P,Q)、(P,Q+1)、(P+1,Q)、(P+1,Q+1)的四个像素点A、B、C、D之间(P,Q都为整数),A、B、C、D的亮度值分别为GA、GB、GC、GD
(S5-4)上述坐标点亮度值的权重系数计算方法具体如下:
Figure BDA00030533782600000912
Figure BDA00030533782600000913
其中,Xij,Yij为上述坐标点的横纵坐标,
Figure BDA00030533782600000914
分别为横纵坐标向下取整。α、β为上述坐标点于像素点A的距离,由于A、B、C、D四点构成了一个边长为单位l的正方形,四个像素点与坐标点之间的距离都可用α、β来表示。像素点距离越近影响越大,A、B、C、D四个像素点对坐标点亮度值的影响权重系数也可用α、β来表示。
(S5-5)通过上述权重系数,计算出空白矩形图像上第i行第j列的亮度值;
进一步地,上述算法具体如下:
Gij=(1-α)×(1-β)×GA+α×(1-β)×GB+(1-α)×β×GC+α×β×GD
其中,Gij为空白矩形图像上第i行第j列的亮度值,GA、GB、GC、GD分别为所计算坐标(Xij,Yij)周围四个像素点A、B、C、D的亮度值。α、β是根据坐标点离周围四个像素点的距离而确定的权重系数。该图以坐标点在第一象限为例,但不管坐标点位于第几象限,权重系数总为正,该公式总是适用。
采用了本发明的用于医疗内窥镜连续成像的管腔图像展开方法,可以将所述环状管腔图像展开为平铺矩形图像,并进行图像展开,呈现出所述连续的管腔内壁平铺图。该发明可方便医生观察管腔内壁的情况,可使医生实时观察完整的病人管腔内壁图像。本发明实时性好,精度高,稳定性好,降低了手术复杂度,有效降低手术难度,提高手术成功率。
在此说明书中,本发明已参照其特定的实施例作了描述。但是,很显然仍可以作出各种修改和变换而不背离本发明的精神和范围。因此,说明书和附图应被认为是说明性的而非限制性的。

Claims (7)

1.一种用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的方法包括以下步骤:
(1)以内窥镜图像的中心为圆心,获取内窥镜图像的最大内接圆;
(2)建立空白矩形图像,以最大内接圆周长为空白矩形图像的长边,以最大内接圆半径的1/2作为空白矩形图像的短边;
(3)计算空白矩形图像上每一个像素点对应于原图像的亮度值,填入空白矩形图像,最终形成完整的平铺图像。
2.根据权利要求1所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述内窥镜图像的中心是内窥镜镜头光轴的中心。
3.根据权利要求1所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的步骤(2)的建立空白矩形图像的步骤具体包括以下处理过程:
(2-1)从圆心到最大内接圆之间包含若干个同心圆,选取半径为最大内接圆半径1/2的同心圆到最大内接圆之间的区域进行展开,从剖开线将同心圆区域展开为平铺矩形图像;
(2-2)建立空白矩形图像,以最大内接圆周长为展开平铺矩形图像的长边,以最大内接圆半径的1/2作为展开平铺矩形图像的短边。
4.根据权利要求1所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的步骤(3)的计算空白矩形图像上每一个像素点亮度值的步骤具体包括以下处理过程:
(3-1)计算获得空白矩形图像上每一个像素点对应于原图像上的坐标;
(3-2)利用原图中上述坐标点周围四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值。
5.根据权利要求4所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的步骤(3-1)具体包括以下处理过程:
(3-1.1)计算出空白矩形每一列像素点的区域宽度对应的弧度;
(3-1.2)计算出空白矩形每一行像素点的对应的同心圆半径;
(3-1.3)由空白矩形图像上每一个像素点所对应的半径和弧度,计算所述像素点对应于原图像上的位置坐标。
6.根据权利要求4所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的步骤(3-2)的利用坐标点周围的四个像素点的亮度值计算出对应的空白矩形图像上像素点的亮度值的步骤具体包括以下处理过程:
(3-2.1)判断横纵坐标是否为整数,如果是,则对应原图像上的像素点,取原图像素点的亮度值填入空白图像对应的像素点;否则,继续步骤(3-2.2);
(3-2.2)通过该坐标点周围的四个像素点的亮度值进行权重计算。
7.根据权利要求6所述的用于医疗内窥镜成像的管腔图像展开方法,其特征在于,所述的步骤(3-2.2)的通过该坐标点周围的四个像素点的亮度值进行权重计算,具体为:
根据以下公式计算亮度值:
Gij=(1-α)×(1-β)×GA+α×(1-β)×GB+(1-α)×β×GC+α×β×GD
其中,Gij为空白矩形图像第i行第j列的像素点的亮度值,GA、GB、GC、GD分别为坐标点周围的四个像素点A、B、C、D的亮度值,α和β为根据坐标点离周围四个像素点的距离而确定的权重系数。
CN202110493564.3A 2021-05-07 2021-05-07 用于医疗内窥镜成像的管腔图像展开方法 Active CN113160040B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110493564.3A CN113160040B (zh) 2021-05-07 2021-05-07 用于医疗内窥镜成像的管腔图像展开方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110493564.3A CN113160040B (zh) 2021-05-07 2021-05-07 用于医疗内窥镜成像的管腔图像展开方法

Publications (2)

Publication Number Publication Date
CN113160040A true CN113160040A (zh) 2021-07-23
CN113160040B CN113160040B (zh) 2022-08-02

Family

ID=76873773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110493564.3A Active CN113160040B (zh) 2021-05-07 2021-05-07 用于医疗内窥镜成像的管腔图像展开方法

Country Status (1)

Country Link
CN (1) CN113160040B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101065052A (zh) * 2004-12-27 2007-10-31 奥林巴斯株式会社 医疗图像处理装置和医疗图像处理方法
WO2008152938A1 (ja) * 2007-06-11 2008-12-18 Hitachi Medical Corporation 医用画像表示装置、医用画像表示方法及び医用画像表示プログラム
CN105261052A (zh) * 2015-11-03 2016-01-20 沈阳东软医疗系统有限公司 管腔图像展开绘制方法及装置
US20160278612A1 (en) * 2013-09-27 2016-09-29 Olympus Corporation Endoscope system
CN111192198A (zh) * 2019-12-26 2020-05-22 台州学院 一种基于管道机器人的管道全景扫描方法
CN111861928A (zh) * 2020-07-24 2020-10-30 深圳开立生物医疗科技股份有限公司 一种确定内窥镜图像亮度的方法、装置及相关设备
CN112089404A (zh) * 2020-10-10 2020-12-18 南京航空航天大学 一种内窥OCT-Raman双模成像装置及成像方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101065052A (zh) * 2004-12-27 2007-10-31 奥林巴斯株式会社 医疗图像处理装置和医疗图像处理方法
WO2008152938A1 (ja) * 2007-06-11 2008-12-18 Hitachi Medical Corporation 医用画像表示装置、医用画像表示方法及び医用画像表示プログラム
US20160278612A1 (en) * 2013-09-27 2016-09-29 Olympus Corporation Endoscope system
CN105261052A (zh) * 2015-11-03 2016-01-20 沈阳东软医疗系统有限公司 管腔图像展开绘制方法及装置
CN111192198A (zh) * 2019-12-26 2020-05-22 台州学院 一种基于管道机器人的管道全景扫描方法
CN111861928A (zh) * 2020-07-24 2020-10-30 深圳开立生物医疗科技股份有限公司 一种确定内窥镜图像亮度的方法、装置及相关设备
CN112089404A (zh) * 2020-10-10 2020-12-18 南京航空航天大学 一种内窥OCT-Raman双模成像装置及成像方法

Also Published As

Publication number Publication date
CN113160040B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
JP4585048B2 (ja) 内視鏡システム
US8939892B2 (en) Endoscopic image processing device, method and program
JP5094036B2 (ja) 内視鏡挿入方向検出装置
US7878971B2 (en) Endoscope inserting direction detecting method
US7830378B2 (en) Medical image processing apparatus and medical image processing method
JP7480477B2 (ja) 医療用観察システム、制御装置及び制御方法
US11954860B2 (en) Image matching method and device, and storage medium
JP7125479B2 (ja) 医療画像処理装置、医療画像処理装置の作動方法及び内視鏡システム
Liu et al. Global and local panoramic views for gastroscopy: an assisted method of gastroscopic lesion surveillance
US20210274089A1 (en) Method and Apparatus for Detecting Missed Areas during Endoscopy
WO2010087057A1 (ja) 内視鏡システム
WO2021171465A1 (ja) 内視鏡システム及び内視鏡システムによる管腔走査方法
CN113160040B (zh) 用于医疗内窥镜成像的管腔图像展开方法
JPH03165732A (ja) 内視鏡の暗部領域の検出方法
US20210169305A1 (en) Image processing apparatus, image processing method, and image processing system
US20240087113A1 (en) Recording Medium, Learning Model Generation Method, and Support Apparatus
JP2005258062A (ja) 内視鏡システム、及び内視鏡装置
CN113317874B (zh) 一种医学图像处理装置及介质
CN110969603A (zh) 病变位置的相对定位方法、装置及终端设备
JP5530406B2 (ja) 内視鏡挿入方向検出装置の作動方法および内視鏡挿入方向検出装置
EP4177664A1 (en) Program, information processing method, and endoscope system
CN115311317A (zh) 一种基于Scaleformer类算法的腹腔镜图像分割方法及系统
CN113808057A (zh) 基于无监督学习的内窥镜图像增强方法
JP5687583B2 (ja) 内視鏡挿入方向検出方法及び内視鏡挿入方向検出装置
JP7155670B2 (ja) 医用画像処理装置、医用画像処理方法、プログラム、及びデータ作成方法

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

Effective date of registration: 20210824

Address after: 200949 building 24, No. 1688, Jinshao Road, Baoshan District, Shanghai

Applicant after: ALTON (SHANGHAI) MEDICAL INSTRUMENTS Co.,Ltd.

Address before: No.99 SHANGDA Road, Baoshan District, Shanghai 200044

Applicant before: Shanghai University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant