CN110517263B - 确定内中膜厚度的方法、装置及存储介质 - Google Patents
确定内中膜厚度的方法、装置及存储介质 Download PDFInfo
- Publication number
- CN110517263B CN110517263B CN201910824455.8A CN201910824455A CN110517263B CN 110517263 B CN110517263 B CN 110517263B CN 201910824455 A CN201910824455 A CN 201910824455A CN 110517263 B CN110517263 B CN 110517263B
- Authority
- CN
- China
- Prior art keywords
- region
- energy
- interest
- determining
- 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
Images
Classifications
-
- 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
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- 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/10132—Ultrasound image
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请公开了一种确定内中膜厚度的方法及装置,属于超声诊断成像领域。在本申请中,终端在待检测图像中获取到感兴趣区域,可以确定该感兴趣区域内每行像素点的累积能量值,根据累积能量值判断感兴趣区域内是否包含血管壁,如果感兴趣区域内包含血管壁则进一步确定感兴趣区域内包含血管壁的部位,根据包含血管壁的部位,进而确定内中膜厚度。由此可见,在本申请实施例中,终端可以自动判断感兴趣区域内是否包含有血管壁,并自动确定感兴趣区域内包含的血管壁的部位,从而确定内中膜厚度,相较于相关技术,无需用户自身去识别待检测图像中的哪部分图像区域内包含的是血管壁的哪个部位,降低了操作的复杂性,提高了诊断效率。
Description
技术领域
本申请涉及超声诊断成像领域,特别涉及一种确定内中膜厚度的方法、装置及存储介质。
背景技术
随着人年龄的增长,脂质可以沉积在颈动脉血管壁上,形成颈动脉粥样硬化。颈动脉粥样硬化可以引发心脏病、脑卒中等心脑血管疾病,因此对颈动脉粥样硬化进行预防、诊断和监控具有非常重要的意义。颈动脉血管壁分三层,内膜、中膜和外膜,通常硬化发生于内膜和中膜之间。基于此,通常可以通过确定颈动脉血管壁的内中膜的厚度来确定是否存在颈动脉粥样硬化。其中,内中膜厚度是指内膜的内边界与中膜的外边界之间的厚度。
相关技术中,颈动脉血管壁可以分为前壁部位和后壁部位。在测量颈动脉血管壁的内中膜厚度时,用户可以首先在终端上选择测量选项,该测量选项用于指示通过前壁部位还是后壁部位来确定内中膜厚度。之后,终端可以显示采集到的包含有前壁部位和后壁部位的待检测图像。此时,用户需要在待检测图像中选择包含有测量选项所指示的部位的感兴趣区域。之后,终端可以在该感兴趣区域内,采用灰度极值搜索的方式,确定中膜的外边界和内膜的内边界,最后,根据中膜的外边界和内膜的内边界确定内中膜厚度。
由此可见,相关技术需要用户根据自身选择的测量选项所指示的部位从待检测图像中选择对应的感兴趣区域,这就要求用户自身能够识别待检测图像中的哪部分图像区域内包含的是血管壁的哪个部位,提高了操作的复杂性,且影响诊断效率。
发明内容
本申请实施例提供了一种确定内中膜厚度的方法、装置及存储介质,可以用于解决相关技术中测量内中膜厚度时操作复杂且诊断效率低的问题。所述技术方案如下:
一方面,提供了一种确定内中膜厚度的方法,所述方法包括:
获取待检测图像中的感兴趣区域;
确定所述感兴趣区域内每行像素点的累积能量值;
如果根据所述感兴趣区域内每行像素点的累积能量值确定所述感兴趣区域内包含血管壁,则确定所述感兴趣区域内包含的血管壁的部位,所述血管壁的部位包括前壁部位和/或后壁部位;
根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度。
可选地,所述确定所述感兴趣区域内每行像素点的累积能量值,包括:
根据所述感兴趣区域内每行像素点的灰度值,确定所述感兴趣区域内每行像素点的能量和;
根据所述感兴趣区域内的每行像素点的能量和,确定所述感兴趣区域内每行像素点的累积能量值。
可选地,所述根据所述感兴趣区域内的每行像素点的能量和,确定所述感兴趣区域内每行像素点的累积能量值,包括:
将所述感兴趣区域内的第一行像素点的能量和作为所述第一行像素点的累积能量值;
将第i行像素点的能量和与所述第i行像素点之前的每行像素点的能量和的总和作为所述第i行像素点的累积能量值,所述i为大于1的整数。
可选地,所述确定所述感兴趣区域内每行像素点的累积能量值之后,还包括:
从所述感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值;
根据所述最小累积能量值和所述最大累积能量值,确定能量斜率;
根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值;
根据确定的多个能量偏差值,判断所述感兴趣区域内是否包含血管壁。
可选地,所述根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值,包括:
根据所述能量斜率和每行像素点的累积能量值,通过下述公式计算每行像素点的能量偏差值:
sub[j]=maxf×j-allpower[j]
其中,sub[j]为第j行像素点的能量偏差值,maxf为能量斜率,allpower[j]为第j行像素点的累积能量值。
可选地,所述根据确定的多个能量偏差值,判断所述感兴趣区域内是否包含血管壁,包括:
当所述多个能量偏差值中包括大于0的多个第一能量偏差值和小于0的多个第二能量偏差值时,如果所述多个第一能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值和/或所述多个第二能量偏差值中存在绝对值大于所述第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均大于0时,如果所述多个能量偏差值中存在大于所述第一参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均小于0时,如果所述多个能量偏差值中存在绝对值大于所述第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁。
可选地,所述确定所述感兴趣区域内包含的血管壁的部位,包括:
如果所述一个或多个目标能量偏差值均大于0,则确定所述感兴趣区域内包含的血管壁的部位为前壁部位;
如果所述一个或多个目标能量偏差值均小于0,则确定所述感兴趣区域内包含的血管壁的部位为后壁部位;
如果所述一个或多个目标能量偏差值包括大于0的目标能量偏差值和小于0的目标能量偏差值,则确定所述感兴趣区域内包含的血管壁的部位包括前壁部位和后壁部位。
可选地,所述根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度,包括:
根据所述感兴趣区域内包含的血管壁的部位,在所述感兴趣区域内确定初始能量演化曲线,所述初始能量演化曲线为闭合曲线;
根据第一图像区域内像素点的灰度值和第二图像区域内像素点的灰度值,确定所述血管壁的内膜的内边界,所述第一图像区域是指所述初始能量演化曲线内的图像区域,所述第二图像区域是指所述初始能量演化曲线外的图像区域;
根据所述内膜的内边界和所述感兴趣区域内包含的血管壁的部位,确定所述血管壁的中膜的外边界;
根据所述内膜的内边界和所述中膜的外边界,确定所述内中膜厚度。
可选地,所述根据所述内膜的内边界和所述感兴趣区域内包含的血管壁的部位,确定所述血管壁的中膜的外边界,包括:
确定所述感兴趣区域内每个像素点的纵向灰度梯度值;
根据所述感兴趣区域内包含的血管壁的部位,确定查找方向;
按照所述查找方向,查找所述内膜的内边界上每个像素点对应的多个像素点中纵向灰度梯度值最大的像素点;
根据查找到的多个像素点,确定所述中膜的外边界。
另一方面,提供了一种确定内中膜厚度的装置,所述装置包括:
第一获取模块,用于获取待检测图像中的感兴趣区域;
第一确定模块,用于确定所述感兴趣区域内每行像素点的累积能量值;
第二确定模块,用于如果根据所述感兴趣区域内每行像素点的累积能量值确定所述感兴趣区域内包含血管壁,则确定所述感兴趣区域内包含的血管壁的部位,所述血管壁的部位包括前壁部位和/或后壁部位;
第三确定模块,用于根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度。
可选地,所述第一确定模块包括:
第一确定子模块,用于根据所述感兴趣区域内每行像素点的灰度值,确定所述感兴趣区域内每行像素点的能量和;
第二确定子模块,用于根据所述感兴趣区域内的每行像素点的能量和,确定所述感兴趣区域内每行像素点的累积能量值。
可选地,所述第二确定子模块具体用于:
将所述感兴趣区域内的第一行像素点的能量和作为所述第一行像素点的累积能量值;
将第i行像素点的能量和与所述第i行像素点之前的每行像素点的能量和的总和作为所述第i行像素点的累积能量值,所述i为大于1的整数。
可选地,所述装置还包括:
第二获取模块,用于从所述感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值;
第四确定模块,用于根据所述最小累积能量值和所述最大累积能量值,确定能量斜率;
第五确定模块,用于根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值;
第六确定模块,用于根据确定的多个能量偏差值,判断所述感兴趣区域内是否包含血管壁。
可选地,所述第五确定模块具体用于:
根据所述能量斜率和每行像素点的累积能量值,通过下述公式计算每行像素点的能量偏差值:
sub[j]=maxf×j-allpower[j]
其中,sub[j]为第j行像素点的能量偏差值,maxf为能量斜率,allpower[j]为第j行像素点的累积能量值。
可选地,所述第六确定模块具体用于:
当所述多个能量偏差值中包括大于0的多个第一能量偏差值和小于0的多个第二能量偏差值时,如果所述多个第一能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值和/或所述多个第二能量偏差值中存在绝对值大于所述第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均大于0时,如果所述多个能量偏差值中存在大于所述第一参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均小于0时,如果所述多个能量偏差值中存在绝对值大于所述第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁。
可选地,所述第六确定模块还用于:
如果所述一个或多个目标能量偏差值均大于0,则确定所述感兴趣区域内包含的血管壁的部位为前壁部位;
如果所述一个或多个目标能量偏差值均小于0,则确定所述感兴趣区域内包含的血管壁的部位为后壁部位;
如果所述一个或多个目标能量偏差值包括大于0的目标能量偏差值和小于0的目标能量偏差值,则确定所述感兴趣区域内包含的血管壁的部位包括前壁部位和后壁部位。
可选地,所述第三确定模块包括:
初始子模块,用于根据所述感兴趣区域内包含的血管壁的部位,在所述感兴趣区域内确定初始能量演化曲线,所述初始能量演化曲线为闭合曲线;
第三确定子模块,用于根据第一图像区域内像素点的灰度值和第二图像区域内像素点的灰度值,确定所述血管壁的内膜的内边界,所述第一图像区域是指所述初始能量演化曲线内的图像区域,所述第二图像区域是指所述初始能量演化曲线外的图像区域;
第四确定子模块,用于根据所述内膜的内边界和所述感兴趣区域内包含的血管壁的部位,确定所述血管壁的中膜的外边界;
第五确定子模块,用于根据所述内膜的内边界和所述中膜的外边界,确定所述内中膜厚度。
可选地,所述第四确定子模块具体用于:
确定所述感兴趣区域内每个像素点的纵向灰度梯度值;
根据所述感兴趣区域内包含的血管壁的部位,确定查找方向;
按照所述查找方向,查找所述内膜的内边界上每个像素点对应的多个像素点中纵向灰度梯度值最大的像素点;
根据查找到的多个像素点,确定所述中膜的外边界。
另一方面,提供了一种确定内中膜厚度的装置,所述装置包括处理器、通信接口、存储器和通信总线;
其中,所述处理器、所述通信接口和所述存储器通过所述通信总线完成相互间的通信;
所述存储器用于存放计算机程序;
所述处理器用于执行所述存储器上所存放的程序,以实现前述提供确定内中膜厚度的方法。
另一方面,提供了一种计算机可读存储介质,所述存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现前述提供的确定内中膜厚度的方法的步骤。
本申请实施例提供的技术方案带来的有益效果至少包括:
在本申请实施例中,终端在待检测图像中获取到感兴趣区域,可以确定该感兴趣区域内每行像素点的累积能量值,根据累积能量值判断感兴趣区域内是否包含血管壁,如果感兴趣区域内包含血管壁则进一步确定感兴趣区域内包含血管壁的部位,感兴趣区域内可以包含前壁部位和/或后壁部位,根据包含血管壁的部位,进而确定内中膜厚度。由此可见,在本申请实施例中,可以由终端自动判断感兴趣区域内是否包含有血管壁,并在确定包含血管壁的情况下,自动识别感兴趣区域内包含的血管壁的部位,从而根据包含的血管壁的部位来确定内中膜厚度,相较于相关技术,无需用户自身去识别待检测图像中的哪部分图像区域内包含的是血管壁的哪个部位,降低了操作的复杂性,提高了诊断效率。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种确定内中膜厚度的***示意图;
图2是本申请实施例提供的一种确定内中膜厚度的方法流程图;
图3是本申请实施例提供的一种血管壁的部位的划分示意图;
图4是本申请实施例提供的一种初始能量演化曲线的位置的示意图;
图5是本申请实施例提供的一种确定内中膜厚度的装置的结构示意图;
图6是本申请实施例提供的一种确定内中膜厚度的终端的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请实施方式作进一步地详细描述。
在对本申请实施例进行详细的解释说明之前,先对本申请实施例涉及的应用场景予以介绍。
颈动脉粥样硬化是一种可以引发心脏病、脑卒中等心脑血管疾病的疾病,颈动脉血管壁分为三层,内膜、中膜和外膜,通常硬化发生于内膜和中膜之间。因此,可以通过确定颈动脉血管壁的内中膜的厚度来确定是否存在颈动脉粥样硬化。颈动脉的内中膜厚度通常可以通过采集的超声图像来测量,本申请实施例提供的确定内中膜厚度的方法,即可以通过对采集到的超声图像进行处理来得到内中膜厚度。
接下来对本申请实施例提供的确定内中膜厚度的方法所涉及的***架构进行介绍。
图1是本申请实施例提供的一种确定内中膜厚度的***示意图。如图1所示,该***100包括探头101、超声成像及显示***102。探头101和超声成像及显示***102之间通过无线或有线方式连接以进行通信。
其中,该探头101为可以采集超声图像的探头。用户把探头101放在患者需要采集超声图像的部位,探头101可以对该部位进行扫描,从而得到超声信号。之后,探头101可以将获取到的超声信号发送给超声成像及显示***102,超声成像及显示***102可以根据该超声信号得到超声图像。其中,该超声图像即为待检测图像。
超声成像及显示***102在得到待检测图像之后,可以获取用户在待检测图像中选择的感兴趣区域,并自动判断获取的感兴趣区域内是否包含血管壁,如果该感兴趣区域内包含血管壁,则该终端可以识别所包含的血管壁的部位,并根据包含的血管壁的部位,进一步确定内中膜厚度。
在本申请实施例中,探头101可以为能够采集超声图像的超声探头,超声成像及显示***102可以为台式机、便携机等超声设备,本申请实施例对此不做限定。
可选的,在一种可能的情况中,超声成像及显示***102可以将得到的超声图像,即待检测图像发送给其他终端,该终端在接收到该待检测图像之后,可以获取用户在待检测图像中选择的感兴趣区域,并自动判断获取的感兴趣区域内是否包含血管壁,如果该感兴趣区域内包含血管壁,则该终端可以识别所包含的血管壁的部位,并根据包含的血管壁的部位,进一步确定内中膜厚度。
也即,在本申请实施例中,用于对待检测图像进行处理以得到内中膜厚度的设备可以是用于成像和显示的超声设备,也可以是除超声设备之外的其他具有数据处理功能的设备,本申请实施例对此不做限定。
接下来对本申请实施例提供的确定内中膜厚度的方法进行介绍。
图2是本申请实施例提供的一种确定内中膜厚度的方法的流程图,该方法可以应用于终端中,该终端可以是指前述***结构中介绍的超声成像及显示***102,也可以是指其他具有数据处理功能的终端。如图2所示,该方法包括以下步骤:
步骤201:获取待检测图像中的感兴趣区域。
在本申请实施例中,终端可以接收探头扫描得到的待检测图像,并将接收到的待检测图像进行显示。用户可以在终端显示的待检测图像中框选一个矩形区域。终端在检测到用户在待检测图像中框选的矩形区域之后,将该矩形区域之外的待检测图像进行裁剪,保留该矩形区域,此时,该矩形区域即为感兴趣区域。
需要说明的是,由于待检测图像为探头采集的超声图像,因此通常会含有大量的散斑噪声,并且会存在伪影,其中,散斑噪声为激光投射在粗糙表面上的光向各个方向反射时,在空间出现随机分布的亮斑和暗斑;伪影为被扫描的物体并不存在而在图像上却出现的各种形态的影像。因此,终端在获取感兴趣区域之后,还可以对感兴趣区域进行去噪处理,以此来提升感兴趣区域的图像质量,进而提高测量准确率。
其中,终端可以通过图像滤波对感兴趣区域进行去噪。图像滤波是一个加权平均的运算过程,滤波后图像中的每个像素点值都是原图像中该像素点邻域内多个像素点值的加权平均,不同的滤波器,最根本的差异就是权重不同。示例性地,在本申请实施例中,终端可以采用双边滤波对感兴趣区域进行去噪,双边滤波在滤波时同时考虑了像素间的几何距离和色彩距离。双边滤波通过空间域权重和像素值域权重两种权重对图像进行加权处理,既能够去除噪声,也能进行边缘保护。终端可以使用下述公式对感兴趣区域进行双边滤波:
其中,EF[1]p为双边滤波函数,Wp为一个标准量,Gσ(x)为二维高斯内核,p和q为感兴趣区域内不同的二维像素点,Ip为像素点p的像素灰度值,Iq为像素点q的像素灰度值,σd为空间域高斯函数的标准差,σr为像素值域高斯函数的标准差,为空间域权重,为像素值域权重。
其中,空间域权重可以通过邻域像素点和中心像素点之间的几何距离确定得到,像素值域权重可以通过邻域像素点和中心像素点之间的色彩距离确定得到。随着邻域像素点和中心像素点之间的几何距离和色彩距离的变化,空间域权重和像素值域权重也会随之变化。当(Ip-Iq)的值变化很小,则说明在图像的平坦区域,对应的像素值域权重接近于1,此时空间域权重起主要作用,对此区域进行高斯模糊,当(Ip-Iq)的值变化较大,则说明在图像的边缘区域,此时像素值域权重会减小,从而保持了边缘的细节信息。
步骤202:确定感兴趣区域内每行像素点的累积能量值。
在获取到感兴趣区域后,终端可以根据感兴趣区域内每行像素点的灰度值,确定感兴趣区域内每行像素点的能量和,根据感兴趣区域内的每行像素点的能量和,确定感兴趣区域内每行像素点的累积能量值。
由前述描述可知,感兴趣区域是一个矩形区域,其中,该矩形区域的长度可以用W来表示,高度可以用H来表示。对于该感兴趣区域内的每行像素点,终端可以将该行所有像素点的灰度值求和,从而得到这一行像素点的能量和。
在确定感兴趣区域内每行像素点的能量和之后,终端可以将感兴趣区域内的第一行像素点的能量和作为第一行像素点的累积能量值,将第i行像素点的能量和与第i行像素点之前的每行像素点的能量和的总和作为第i行像素点的累积能量值,i为大于1的整数。
需要说明的是,终端可以沿感兴趣区域的高度方向,将感兴趣区域的下边界的像素点作为第一行像素点,上边界为最后一行像素点。这样,终端可以将感兴趣区域内第一行像素点的能量和作为第一行像素点的累积能量值,将第二行像素点的能量和加上第一行像素点的能量和,得到第二行像素点的累积能量值,将第三行像素点的能量和加上第二行像素点的能量和,再加上第一行像素点的能量和,得到第三行像素点的累积能量值,以此类推。
步骤203:根据感兴趣区域内每行像素点的累积能量值,判断该感兴趣区域内是否包含有血管壁。
在确定该感兴趣区域内每行像素点的累积能量值之后,该终端可以根据每行像素点的累积能量值,来判断该感兴趣区域内是否包含有血管壁。
其中,终端可以从感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值,根据最小累积能量值和最大累积能量值,确定能量斜率,其中,该能量斜率为最大累积能量值与最小累积能量值之间的差值与该感兴趣区域内包括的像素点的行数的比值。根据能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值,根据确定的多个能量偏差值,判断该感兴趣区域内是否包含血管壁。
需要说明的是,将感兴趣区域包括的多行像素点的累积能量值两两依次进行比较,每次比较后保留较大的累积能量值,全部行的像素点的累积能量值比较完成后,当前保留的累积能量值即为最大累积能量值。同样的,将感兴趣区域包括的多行像素点的累积能量值两两依次进行比较,每次比较后保留较小的累积能量值,全部行的像素点的累积能量值比较完成后,当前保留的累积能量值即为最小累积能量值。在获取到最小累积能量值和最大累积能量值后,终端可以根据该最小累积能量值和最大累积能量值,通过下述公式计算能量斜率。
其中,maxf为能量斜率,allpower[H]为最大累积能量值,allpower[1]为最小累积能量值。H为该感兴趣区域内包含的像素点的行数。由前述确定累积能量值的过程可知,累积能量值为当前行像素点的能量和与该行之前的每行像素点的能量和的总和,因此,最大累积能量值即为该感兴趣区域内多行像素点中最后一行像素点的累积能量值,最小累积能量值即为该感兴趣区域内第一行像素点的能量累积值。
需要说明的是,如果该感兴趣区域内的各行像素点的能量和是相同的,也即,各行像素点的灰度值是均匀的,则从第一行像素点到最后一行像素点,累积能量值的变化率将是相等。在本申请实施例中,能量斜率即为假设各行像素点的能量和相同的情况下,确定得到的累积能量值的参考变化率。后续通过该能量斜率,即可以确定出各行像素点的参考累积能量值,进而将各行像素点的参考累积能量值与实际累积能量值进行比较,得到各行像素点的能量偏差值。
示例性地,终端可以根据该能量斜率和每行像素点的累积能量值通过下述公式计算每行像素点的能量偏差值:
sub[j]=maxf×j-allpower[j]
其中,sub[j]为能量偏差值,maxf为能量斜率,j用于指示该感兴趣区域内的第j行像素点,allpower[j]为第j行像素点的累积能量值。
由此可见,能量偏差值实际上表征了每行像素点的实际累积能量值与该行像素点的参考累积能量值之间的偏差,通过各行像素点的能量偏差值,可以反映出该感兴趣区域内的像素点的累积能量值的变化率。
需要说明的是,通常,在图像中,包含血管腔的图像区域的像素点的灰度值较包含血管壁的图像区域的像素点的灰度值小,而感兴趣区域内每行像素点的能量和是通过对感兴趣区域内每行像素点的灰度值求和得到,在这种情况下,如果感兴趣区域内某行像素点的能量和较小,则说明该行像素点对应的可能是血管腔,如果感兴趣区域内某行像素点的能量和较大,则说明该行像素点对应的可能是血管壁。由于感兴趣区域内每行像素点的累积能量值为当前行像素点的能量和与该行之前的每行像素点的能量和的总和,因此,当感兴趣区域包含血管壁时,在血管壁和血管腔的交界处,像素点的累积能量值的变化率较大,得到的能量偏差值也会较大。
值得注意的是,如果感兴趣区域内仅包含有血管腔不包含血管壁,则该感兴趣区域内的各行像素点的累积能量值的变化将较为均匀,在这种情况下,得到的能量偏差值也将较小。基于此,在得到多个能量偏差值之后,终端可以通过设定的参考阈值来判断该感兴趣区域内是否包含血管壁。另外,由于包含的血管壁部位的不同,因此,多个能量偏差值中可能有大于0的值,也可能有小于0的值。
其中,当多个能量偏差值中既存在大于0也存在小于0的值时,终端可以判断大于0的多个第一能量偏差值中是否存在大于第一参考阈值的目标能量偏差值,同时判断小于0的多个第二能量偏差值中是否存在绝对值大于第二参考阈值的能量偏差值。如果多个第一能量偏差值中存在大于第一参考阈值的能量偏差值,或者是多个第二能量偏差值中存在绝对值大于第二参考阈值的能量偏差值,或者是,多个第一能量偏差值中存在绝对值大于第一参考阈值的能量偏差值且多个第二能量偏差值中也存在绝对值大于第二参考阈值的能量偏差值,则说明该感兴趣区域内包含有血管壁。需要说明的是,在这种情况下,多个第一能量偏差值中大于第一参考阈值的能量偏差值和多个第二能量偏差值中绝对值大于第二参考阈值的能量偏差值均可以称为目标能量偏差值。
当确定的多个能量偏差值均大于0时,则终端可以直接判断该多个能量偏差值中是否存在大于第一参考阈值的能量偏差值,如果存在,则说明该感兴趣区域内存在血管壁。在这种情况下,可以将多个能量偏差值中大于第一参考阈值的能量偏差值称为目标能量偏差值。
当确定的多个能量偏差值均小于0时,则终端可以直接判断该多个能量偏差值中是否存在绝对值大于第二参考阈值的目标能量偏差值,如果存在,则可以确定该感兴趣区域内包含血管壁。在这种情况下,可以将多个能量偏差值中绝对值大于第二参考阈值的能量偏差值称为目标能量偏差值。
其中,第一参考阈值和第二参考阈值可以是根据先验经验设置的一个能量偏差值,该参考阈值可以等于或略大于图像中仅包含有血管腔时各行像素点的最大能量偏差值或平均能量偏差值。可选地,在另一种可能的实现方式中,第一参考阈值可以是确定出的多个能量偏差值中大于0的能量偏差值的平均值,第二参考阈值则可以为多个能量偏差值中小于0的能量偏差值的平均值。
可选地,如果确定的多个能量偏差值中大于0的能量偏差值中不存在大于第一参考阈值目标能量偏差值,且小于0的能量偏差值中不存在绝对值大于第二参考阈值的目标能量偏差值,则说明感兴趣区域中每行像素点的累积能量值间的变化率较小,也即,感兴趣区域内每行像素点间的灰度值的大小比较平均,此时,可以确定该感兴趣区域内包含的为血管腔部分,并没有包含血管壁。由于感兴趣区域内并未包含有血管壁,因此,在这种情况下,终端可以直接结束操作,或者是显示提示信息,以提示用户重新选择感兴趣区域,并重新通过本申请实施例的方法来进行内中膜厚度的测量。
步骤204:如果根据感兴趣区域内每行像素点的累积能量值确定感兴趣区域内包含血管壁,则确定感兴趣区域内包含的血管壁的部位。
其中,血管壁包括内膜、中膜和外膜。内膜围成了血管腔,中膜位于内膜外部,外膜位于中膜外部且与外层组织相邻。并且,内膜与血管腔的边界为内膜的内边界,中膜与外膜的边界为中膜的外边界。由于血管近似圆柱形,因此,对于包含有血管腔和血管壁的血管纵截面,以该纵截面的水平中线为分界,可以将血管壁划分为前壁部位和/或后壁部位。图3示出了一种血管壁的部位的划分示意图。如图3所示,水平中线以上部分的血管壁为前壁部位301,水平中线以下部分的血管壁为后壁部位303。其中,血管壁包含外膜304,中膜305,内膜306,内膜的内边界为3061,中膜的外边界3051。内膜306围成的空腔为血管腔302。内中膜厚度为中膜305与内膜306的厚度之和,测量内中膜厚度即为测量内膜的内边界3061与中膜的外边界3051之间的距离。另外,如图3中所示,对于血管壁中的前壁部位301,中膜305位于内膜306之上,外膜304位于中膜305之上。对于血管壁中的后壁部位303,中膜305位于内膜306之下,外膜304位于中膜305之下。
如果通过步骤203确定当前的感兴趣区域内包含血管壁,则终端可以通过前述的多个能量偏差值中存在的目标能量偏差值的正负来进一步确定该感兴趣区域内包含的血管壁的部位为前壁部位还是后壁部位或者是同时包含前壁部位和后壁部位。
其中,当获得的一个或多个目标能量偏差值均大于0时,则确定感兴趣区域内包含的血管壁的部位为前壁部位,当获得的一个或多个目标能量偏差值均小于0,则确定感兴趣区域内包含的血管壁的部位为后壁部位。当获得的一个或多个目标能量偏差值包括大于0的目标能量偏差值和小于0的目标能量偏差值,则确定感兴趣区域内包含的血管壁的部位包括前壁部位和后壁部位。
需要说明的是,当一个或多个目标能量偏差值均大于0时,说明从该感兴趣区域的下边界到上边界,存在累积能量值的增量突增的情况,也即,某几行的像素点的灰度值突然增大了,由于血管腔的灰度值较小,血管壁的灰度值较大,因此可知,沿着该感兴趣区域的下边界到上边界的方向,从血管腔过渡到了血管壁。也即,血管腔位于血管壁的下方,此时,基于前述图3中的介绍可知,当前该感兴趣区域中包含的是血管壁的前壁部位。
当一个或多个目标能量偏差值均小于0时,说明从该感兴趣区域的下边界到上边界,存在累积能量值的增量突减的情况,也即,某几行的像素点的灰度值突然减小了,由于血管腔的灰度值较小,血管壁的灰度值较大,因此可知,沿着该感兴趣区域的下边界到上边界的方向,从血管臂过渡到了血管腔。也即,血管腔位于血管壁的上方,此时,基于前述图3中的介绍可知,当前该感兴趣区域中包含的是血管壁的后壁部位。
当多个目标能量偏差值既包括大于0的能量偏差值也包括小于0的能量偏差值时,则基于前述介绍可知,当前该感兴趣区域同时包含血管壁的前壁部位和后壁部位。
步骤205:根据感兴趣区域内包含的血管壁的部位,确定内中膜厚度。
在确定感兴趣区域内包含的血管壁的部位之后,终端可以根据该感兴趣区域内包含的血管壁的部位,在感兴趣区域内确定初始能量演化曲线,初始能量演化曲线为闭合曲线;根据第一图像区域内像素点的灰度值和第二图像区域内像素点的灰度值,确定血管壁的内膜的内边界,第一图像区域是指初始能量演化曲线内的图像区域,第二图像区域是指所述初始能量演化曲线外的图像区域;根据内膜的内边界和感兴趣区域内包含的血管壁的部位,确定血管壁的中膜的外边界;根据内膜的内边界和所述中膜的外边界,确定内中膜厚度。
其中,初始能量演化曲线为一个闭合曲线,且该闭合曲线可以位于血管腔内。根据感兴趣区域内包含的血管壁的部位的不同,血管腔在感兴趣区域内所处的位置也不同,相应地,确定的初始能量演化曲线的位置不同。
示例性地,当感兴趣区域内包含血管壁的前壁部位时,如图3所示,前壁部位的中膜位于内膜上方,外膜位于中膜上方,也即血管腔位于前壁部位的内膜的下方,因此,感兴趣区域的下部将包含有血管腔,此时,初始能量演化曲线可以位于靠近感兴趣区域下边界的位置;当感兴趣区域内包含血管壁的后壁部位时,由于后壁部位的中膜位于内膜下方,外膜位于中膜下方,因此可知,血管腔位于后壁部位的内膜的上方,因此,感兴趣区域的上部将包含有血管腔,此时,初始能量演化曲线可以位于靠近感兴趣区域上边界的位置;当感兴趣区域内同时包含血管壁的前壁部分和后壁部分时,血管腔将位于前壁部位和后壁部位之间,此时,该血管腔将位于感兴趣区域的中部,初始能量演化曲线可以位于感兴趣区域中间的位置。
示例性地,假设感兴趣区域的左上点坐标为(0,0),右上点坐标为(W,0),左下点坐标为(0,H),右下点坐标为(W,H)。当感兴趣区域只包含血管壁的前壁部位时,参见图4,初始能量演化曲线401的左上点坐标可以为(m,H-m),右上点坐标可以为(W-m,H-m),左下点坐标可以为(m,H-n),右下点坐标可以为(W-m,H-n),其中,m大于n;当感兴趣区域只包含血管壁的后壁部位时,初始能量演化曲线的左上点坐标为(m,n),右上点坐标为(W-m,n),左下点坐标为(m,m),右下点坐标为(W-m,m);当感兴趣区域同时包含血管壁的后壁部位和后壁部位时,初始能量演化曲线的左上点坐标为(m,(H/2)-n),右上点坐标为(W-m,(H/2)-n),左下点坐标为(m,(H/2)+n),右下点坐标为(W-m,(H/2)+n),其中,m、n∈[2,5],即m和n可以为2至5这个闭区间的任意一个值。
通过该初始能量演化曲线,可以将感兴趣区域分为内侧、外侧和边缘三个区域。基于此,在确定初始能量演化曲线后,终端可以通过下述水平集函数初始化能量泛函:
其中,x表示感兴趣区域内的像素点,表示初始的水平集函数,水平集函数通常可以初始化为距离函数,p为任意整数,该整数表示x到初始能量演化曲线所围平面的距离,当x在该平面外时,为负号,当x在该平面内时,为正号,这样通过水平集函数随时间运动的情况,就可以用函数值为零的点的集合来描述能量演化曲线所围的平面,也即可以描述能量演化曲线。
确定能量演化曲线内侧和外侧的灰度均值后,终端可以确定有符号压力函数,示例性地,终端可以根据感兴趣区域中某一点的像素灰度值和能量演化曲线内侧的灰度均值和外侧的灰度均值,通过下述公式确定有符号压力函数。
其中,SPF为有符号压力函数,x为感兴趣区域内的某一点,I(x)为感兴趣区域中某一点的像素灰度值,c1为能量演化曲线内侧的灰度均值、c2为能量演化曲线外侧的灰度均值。
需要说明的是,当能量演化曲线外侧的灰度均值大于内侧的灰度均值时,有符号压力函数大于0,当能量演化曲线外侧的灰度均值小于内侧的灰度均值,则有符号压力函数小于0。有符号压力函数的值的符号不同,则能量演化曲线演化的方向不同。当有符号压力函数大于0时,能量演化曲线演化的方向为向外演化,当有符号压力函数小于0时,能量演化曲线演化的方向为向内演化。
在确定能量演化曲线演化的方向后,终端可以根据有符号压力函数和能量演化曲线上的每一点的像素梯度确定能量演化曲线演化的增量,也即演化的速度。其中,像素梯度为在感兴趣区域中像素点灰度值的变化率。当能量演化曲线上的像素梯度较小时,说明该能量演化曲线离内膜的内边界较远,此时,可以采用较大的演化速度对能量演化曲线进行演化,当像素梯度较大时,说明该能量演化曲线距离内膜的内边界较近,此时可以采用较小的演化速度对能量演化曲线进行演化。
示例性地,可以根据有符号压力函数和每一点的像素梯度通过下述公式计算能量演化曲线演化的速度:
其中,x表示感兴趣区域内的某一点,SPF表示有符号压力函数,α可以用于控制能演化曲线收缩或扩张的能力,α越大则能量演化曲线收缩或扩张的能力就越强,示例性地,α可以取1.5至3之间的任意实数,表示梯度即的变化量。
示例性地,经过演化的能量演化曲线可以通过下述公式表示:
需要说明的是,初始能量演化曲线经过一次演化后,得到一个新的能量演化曲线,将该新的能量演化曲线作为初始能量演化曲线,继续进行下一次的演化,以此类推,终端可以进行指定次数的演化。将最后一个演化得到的演化曲线作为内膜的内边界。其中。每一次演化能量演化曲线后,水平集函数经过演化后就不再是距离函数,而能量演化曲线演化的过程是一个关于的函数,因此不正确的就会导致不正确的能量演化曲线的演化。因此,为了避免这一影响,需要在每一次演化结束后对水平集函数进行重新初始化,将水平集函数重新变为距离函数。示例性地,可以在每一次演化结束后重新给赋值,对于令其他情况下令其中,指定次数可以是用户设置的次数,通常,该指定次数可以设置为100至150次。
可选地,在本申请实施例中,在每一次演化结束之后,终端还可以对每一次演化后得到的新的能量演化曲线进行优化,示例性地,终端可以通过高斯滤波函数将得到的能量演化曲线进行处理,从而使得能量演化曲线平滑。
在确定内膜的内边界后,终端可以确定感兴趣区域内每个像素点的纵向灰度梯度值,根据感兴趣区域内包含的血管壁的部位,确定查找方向,按照查找方向,查找内膜的内边界上每个像素点对应的多个像素点中纵向灰度梯度值最大的像素点,根据查找到的多个像素点,确定中膜的外边界。
需要说明的是,终端将得到的内膜的内边界离散为多个像素点,通过多个像素点分别向外查找灰度梯度最大值的像素点。其中,由图3可知,当感兴趣区域包含血管壁的前壁部位时,血管腔位于血管壁的前壁部位的下方,此时,血管壁中膜的外边界位于内膜的内边界的上方,因此,终端可以确定此时的查找方向为向上查找,也即,终端可以从内膜的内边界向上查找。当感兴趣区域包含的为血管壁的后壁部位时,血管壁中膜的外边界位于内膜的内边界的下方,此时,终端可以确定查找方向为向下查找,也即,终端可以从内膜的内边界开始向下查找。
示例性地,当感兴趣区域包含血管壁的前壁部位,终端可以从内膜的内边界上的每个像素点向上的第三个像素点作为起点,向上查找每个像素点对应的灰度梯度最大值的像素点,当感兴趣区域包含血管壁的后壁部位,终端可以从内膜的内边界上的每个像素点向下的第三个像素点作为起点,向下查找每个像素点对应的灰度梯度最大值的像素点,得到每个像素点对应的目标像素点。在获取到多个目标像素点后,终端可以根据多个目标像素点确定中膜的外边界。
其中,在一种可能的实现方式中,终端可以直接将该多个目标像素点连接起来,从而得到中膜的外边界。
可选地,在另一种可能的实现方式中,终端可以计算内膜的内边界上的每个像素点与该像素点对应的目标像素点之间的距离,从而得到多个距离。确定多个距离的距离平均值,并计算多个距离中的每个距离与该距离平均值之间的差值,得到多个差值。从该多个差值中获取大于预设阈值的差值。对于获取到的每个差值,终端可以按照各个差值对应的目标像素点所在列从左到右的顺序,依次对各个目标像素点进行优化。具体地,对于任一目标像素点,终端可以将该目标像素点用前一个像素点代替,也即用左邻侧已经经过优化的像素点代替。依次进行优化完成后,终端可以将优化后的多个像素点连接起来,从而得到中膜的外边界。
在确定内膜的内边界和中膜的外边界后,终端可以将内膜的内边界和中膜的外边界离散为像素点,之后,终端可以沿血管壁的法线方向确定内膜的内边界上的每个像素点和中膜的外边界上对应的像素点之间的像素距离,并按照图像坐标系和实际坐标系之间的转换关系将该多个像素距离转换为多个物理距离。之后,终端可以根据该多个物理距离确定内中膜厚度。示例性地,终端可以将该多个物理距离中的最大值作为内中膜厚度进行显示,或者,终端可以确定多个物理距离的平均值,将该平均值作为最终的内中膜厚度进行显示。
在本申请实施例中,终端在待检测图像中获取到感兴趣区域,可以确定该感兴趣区域内每行像素点的累积能量值,根据累积能量值判断感兴趣区域内是否包含血管壁,如果感兴趣区域内包含血管壁则进一步确定感兴趣区域内包含血管壁的部位,感兴趣区域内可以包含前壁部位和/或后壁部位,根据包含血管壁的部位,进而确定内中膜厚度。由此可见,在本申请实施例中,可以由终端自动判断感兴趣区域内是否包含有血管壁,并在确定包含血管壁的情况下,自动识别感兴趣区域内包含的血管壁的部位,从而根据包含的血管壁的部位来确定内中膜厚度,相较于相关技术,无需用户自身去识别待检测图像中的哪部分图像区域内包含的是血管壁的哪个部位,降低了操作的复杂性,提高了诊断效率。
接下来对本申请实施例提供的确定内中膜厚度的装置进行介绍。
参见图5,本申请实施例提供了一种确定内中膜厚度的装置500,该确定内中膜厚度的装置可以应用于客户端中,该装置500包括:
第一获取模块501,用于获取待检测图像中的感兴趣区域;
第一确定模块502,用于确定感兴趣区域内每行像素点的累积能量值;
第二确定模块503,用于如果根据感兴趣区域内每行像素点的累积能量值确定感兴趣区域内包含血管壁,则确定感兴趣区域内包含的血管壁的部位,血管壁的部位包括前壁部位和/或后壁部位;
第三确定模块504,用于根据感兴趣区域内包含的血管壁的部位,确定内中膜厚度。
可选地,该第一确定模块502包括:
第一确定子模块,用于根据感兴趣区域内每行像素点的灰度值,确定感兴趣区域内每行像素点的能量和;
第二确定子模块,用于根据感兴趣区域内的每行像素点的能量和,确定感兴趣区域内每行像素点的累积能量值。
可选地,该第二确定子模块具体用于:
将感兴趣区域内的第一行像素点的能量和作为第一行像素点的累积能量值;
将第i行像素点的能量和与第i行像素点之前的每行像素点的能量和的总和作为第i行像素点的累积能量值,i为大于1的整数。
可选地,该装置还包括:
第二获取模块,用于从感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值;
第四确定模块,用于根据最小累积能量值和最大累积能量值,确定能量斜率;
第五确定模块,用于根据能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值;
第六确定模块,用于根据确定的多个能量偏差值,判断感兴趣区域内是否包含血管壁。
可选地,该第五确定模块具体用于:
根据能量斜率和每行像素点的累积能量值,通过下述公式计算每行像素点的能量偏差值:
sub[j]=maxf×j-allpower[j]
其中,sub[j]为第j行像素点的能量偏差值,maxf为能量斜率,allpower[j]为第j行像素点的累积能量值。
可选地,该第六确定模块具体用于:
当多个能量偏差值中包括大于0的多个第一能量偏差值和小于0的多个第二能量偏差值时,如果多个第一能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值和/或多个第二能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定感兴趣区域内包含血管壁;
当多个能量偏差值均大于0时,如果多个能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值,则确定感兴趣区域内包含血管壁;
当多个能量偏差值均小于0时,如果多个能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定感兴趣区域内包含血管壁。
可选地,该第六确定模块还用于:
如果一个或多个目标能量偏差值均大于0,则确定感兴趣区域内包含的血管壁的部位为前壁部位;
如果一个或多个目标能量偏差值均小于0,则确定感兴趣区域内包含的血管壁的部位为后壁部位;
如果一个或多个目标能量偏差值包括大于0的目标能量偏差值和小于0的目标能量偏差值,则确定感兴趣区域内包含的血管壁的部位包括前壁部位和后壁部位。
可选地,该第三确定模块504包括:
初始子模块,用于根据感兴趣区域内包含的血管壁的部位,在感兴趣区域内确定初始能量演化曲线,初始能量演化曲线为闭合曲线;
第三确定子模块,用于根据第一图像区域内像素点的灰度值和第二图像区域内像素点的灰度值,确定血管壁的内膜的内边界,第一图像区域是指初始能量演化曲线内的图像区域,第二图像区域是指初始能量演化曲线外的图像区域;
第四确定子模块,用于根据内膜的内边界和感兴趣区域内包含的血管壁的部位,确定血管壁的中膜的外边界;
第五确定子模块,用于根据内膜的内边界和中膜的外边界,确定内中膜厚度。
可选地,该第四确定子模块具体用于:
确定感兴趣区域内每个像素点的纵向灰度梯度值;
根据感兴趣区域内包含的血管壁的部位,确定查找方向;
按照查找方向,查找内膜的内边界上每个像素点对应的多个像素点中纵向灰度梯度值最大的像素点;
根据查找到的多个像素点,确定中膜的外边界。
在本申请实施例中,终端在待检测图像中获取到感兴趣区域,可以确定该感兴趣区域内每行像素点的累积能量值,根据累积能量值判断感兴趣区域内是否包含血管壁,如果感兴趣区域内包含血管壁则进一步确定感兴趣区域内包含血管壁的部位,感兴趣区域内可以包含前壁部位和/或后壁部位,根据包含血管壁的部位,进而确定内中膜厚度。由此可见,在本申请实施例中,可以由终端自动判断感兴趣区域内是否包含有血管壁,并在确定包含血管壁的情况下,自动识别感兴趣区域内包含的血管壁的部位,从而根据包含的血管壁的部位来确定内中膜厚度,相较于相关技术,无需用户自身去识别待检测图像中的哪部分图像区域内包含的是血管壁的哪个部位,降低了操作的复杂性,提高了诊断效率。
需要说明的是:上述实施例提供的确定内中膜厚度的装置在确定内中膜厚度时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的确定内中膜厚度的方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
图6是根据一示例性实施例示出的一种确定内中膜厚度终端600的结构框图。其中,该终端600可以是笔记本电脑、台式电脑、智能手机或平板电脑等。
通常,终端600包括有:处理器601和存储器602。
处理器601可以包括一个或多个处理核心,比如4核心处理器、8核心处理器等。处理器601可以采用DSP(Digital Signal Processing,数字信号处理)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)、PLA(Programmable Logic Array,可编程逻辑阵列)中的至少一种硬件形式来实现。处理器601也可以包括主处理器和协处理器,主处理器是用于对在唤醒状态下的数据进行处理的处理器,也称CPU(Central ProcessingUnit,中央处理器);协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中,处理器601可以在集成有GPU(Graphics Processing Unit,图像处理器),GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中,处理器601还可以包括AI(Artificial Intelligence,人工智能)处理器,该AI处理器用于处理有关机器学习的计算操作。
存储器602可以包括一个或多个计算机可读存储介质,该计算机可读存储介质可以是非暂态的。存储器602还可包括高速随机存取存储器,以及非易失性存储器,比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中,存储器602中的非暂态的计算机可读存储介质用于存储至少一个指令,其中,,该至少一个指令用于被处理器601所执行以实现本申请中方法实施例提供的确定内中膜厚度的方法。
在一些实施例中,终端600还可选包括有:***设备接口603和至少一个***设备。处理器601、存储器602和***设备接口603之间可以通过总线或信号线相连。各个***设备可以通过总线、信号线或电路板与***设备接口603相连。具体地,***设备包括:射频电路604、触摸显示屏605、摄像头606、音频电路607、定位组件608和电源609中的至少一种。
***设备接口603可被用于将I/O(Input/Output,输入/输出)相关的至少一个***设备连接到处理器601和存储器602。在一些实施例中,处理器601、存储器602和***设备接口603被集成在同一芯片或电路板上;在一些其他实施例中,处理器601、存储器602和***设备接口603中的任意一个或两个可以在单独的芯片或电路板上实现,本实施例对此不加以限定。
射频电路604用于接收和发射RF(Radio Frequency,射频)信号,也称电磁信号。射频电路604通过电磁信号与通信网络以及其他通信设备进行通信。射频电路604将电信号转换为电磁信号进行发送,或者,将接收到的电磁信号转换为电信号。可选地,射频电路604包括:天线***、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路604可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于:万维网、城域网、内联网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity,无线保真)网络。在一些实施例中,射频电路604还可以包括NFC(Near Field Communication,近距离无线通信)有关的电路,本申请对此不加以限定。
显示屏605用于显示UI(User Interface,用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏605是触摸显示屏时,显示屏605还具有采集在显示屏605的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器601进行处理。此时,显示屏605还可以用于提供虚拟按钮和/或虚拟键盘,也称软按钮和/或软键盘。在一些实施例中,显示屏605可以为一个,设置终端600的前面板;在另一些实施例中,显示屏605可以为至少两个,分别设置在终端600的不同表面或呈折叠设计;在再一些实施例中,显示屏605可以是柔性显示屏,设置在终端600的弯曲表面上或折叠面上。甚至,显示屏605还可以设置成非矩形的不规则图形,也即异形屏。显示屏605可以采用LCD(LiquidCrystal Display,液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。需要说明的是,在本申请实施例中,当该终端600为横屏终端时,该终端600的显示屏的宽高比大于1,例如,该终端600的显示屏的宽高比可以为16:9或4:3。当该终端600为竖屏终端时,则该终端600的显示屏的宽高比小于1,例如,该终端600的显示屏的宽高比可以为9:18或3:4等。
摄像头组件606用于采集图像或视频。可选地,摄像头组件606包括前置摄像头和后置摄像头。通常,前置摄像头设置在终端的前面板,后置摄像头设置在终端的背面。在一些实施例中,后置摄像头为至少两个,分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种,以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality,虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中,摄像头组件606还可以包括闪光灯。闪光灯可以是单色温闪光灯,也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合,可以用于不同色温下的光线补偿。
音频电路607可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波,并将声波转换为电信号输入至处理器601进行处理,或者输入至射频电路604以实现语音通信。出于立体声采集或降噪的目的,麦克风可以为多个,分别设置在终端600的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器601或射频电路604的电信号转换为声波。扬声器可以是传统的薄膜扬声器,也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时,不仅可以将电信号转换为人类可听见的声波,也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中,音频电路607还可以包括耳机插孔。
定位组件608用于定位终端600的当前地理位置,以实现导航或LBS(LocationBased Service,基于位置的服务)。定位组件608可以是基于美国的GPS(GlobalPositioning System,全球定位***)、中国的北斗***或欧盟的伽利略***的定位组件。
电源609用于为终端600中的各个组件进行供电。电源609可以是交流电、直流电、一次性电池或可充电电池。当电源609包括可充电电池时,该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池,无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
在一些实施例中,终端600还包括有一个或多个传感器610。该一个或多个传感器610包括但不限于:加速度传感器611、陀螺仪传感器612、压力传感器613、指纹传感器614、光学传感器615以及接近传感器616。
加速度传感器611可以检测以终端600建立的坐标系的三个坐标轴上的加速度大小。比如,加速度传感器611可以用于检测重力加速度在三个坐标轴上的分量。处理器601可以根据加速度传感器611采集的重力加速度信号,控制触摸显示屏605以横向视图或纵向视图进行用户界面的显示。加速度传感器611还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器612可以检测终端600的机体方向及转动角度,陀螺仪传感器612可以与加速度传感器611协同采集用户对终端600的3D动作。处理器601根据陀螺仪传感器612采集的数据,可以实现如下功能:动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器613可以设置在终端600的侧边框和/或触摸显示屏605的下层。当压力传感器613设置在终端600的侧边框时,可以检测用户对终端600的握持信号,由处理器601根据压力传感器613采集的握持信号进行左右手识别或快捷操作。当压力传感器613设置在触摸显示屏605的下层时,由处理器601根据用户对触摸显示屏605的压力操作,实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器614用于采集用户的指纹,由处理器601根据指纹传感器614采集到的指纹识别用户的身份,或者,由指纹传感器614根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时,由处理器601授权该用户执行相关的敏感操作,该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器614可以被设置终端600的正面、背面或侧面。当终端600上设置有物理按键或厂商Logo时,指纹传感器614可以与物理按键或厂商Logo集成在一起。
光学传感器615用于采集环境光强度。在一个实施例中,处理器601可以根据光学传感器615采集的环境光强度,控制触摸显示屏605的显示亮度。具体地,当环境光强度较高时,调高触摸显示屏605的显示亮度;当环境光强度较低时,调低触摸显示屏605的显示亮度。在另一个实施例中,处理器601还可以根据光学传感器615采集的环境光强度,动态调整摄像头组件606的拍摄参数。
接近传感器616,也称距离传感器,通常设置在终端600的前面板。接近传感器616用于采集用户与终端600的正面之间的距离。在一个实施例中,当接近传感器616检测到用户与终端600的正面之间的距离逐渐变小时,由处理器601控制触摸显示屏605从亮屏状态切换为息屏状态;当接近传感器616检测到用户与终端600的正面之间的距离逐渐变大时,由处理器601控制触摸显示屏605从息屏状态切换为亮屏状态。
也即是,本申请实施例不仅提供了一种终端,包括处理器和用于存储处理器可执行指令的存储器,其中,处理器被配置为执行图2所示的确定内中膜厚度的方法,而且,本申请实施例还提供了一种计算机可读存储介质,该存储介质内存储有计算机程序,该计算机程序被处理器执行时可以实现图2所示的确定内中膜厚度的方法。
本申请实施例还提供了一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行上述图2所示实施例提供的确定内中膜厚度的方法。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本申请的可选实施例,并不用以限制本申请,凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (9)
1.一种确定内中膜厚度的方法,其特征在于,所述方法包括:
获取待检测图像中的感兴趣区域;
确定所述感兴趣区域内每行像素点的累积能量值;
从所述感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值;
根据所述最小累积能量值和所述最大累积能量值,确定能量斜率;
根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值;
根据确定的多个能量偏差值,判断所述感兴趣区域内是否包含血管壁;
如果根据所述感兴趣区域内每行像素点的累积能量值确定所述感兴趣区域内包含血管壁,则确定所述感兴趣区域内包含的血管壁的部位,所述血管壁的部位包括前壁部位和/或后壁部位;
根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度;
所述根据确定的多个能量偏差值,判断所述感兴趣区域内是否包含血管壁,包括:
当所述多个能量偏差值中包括大于0的多个第一能量偏差值和小于0的多个第二能量偏差值时,如果所述多个第一能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值和/或所述多个第二能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均大于0时,如果所述多个能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当所述多个能量偏差值均小于0时,如果所述多个能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁。
2.根据权利要求1所述的方法,其特征在于,所述确定所述感兴趣区域内每行像素点的累积能量值,包括:
根据所述感兴趣区域内每行像素点的灰度值,确定所述感兴趣区域内每行像素点的能量和;
根据所述感兴趣区域内的每行像素点的能量和,确定所述感兴趣区域内每行像素点的累积能量值。
3.根据权利要求2所述的方法,其特征在于,所述根据所述感兴趣区域内的每行像素点的能量和,确定所述感兴趣区域内每行像素点的累积能量值,包括:
将所述感兴趣区域内的第一行像素点的能量和作为所述第一行像素点的累积能量值;
将第i行像素点的能量和与所述第i行像素点之前的每行像素点的能量和的总和作为所述第i行像素点的累积能量值,所述i为大于1的整数。
4.根据权利要求1所述的方法,其特征在于,所述根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值,包括:
根据所述能量斜率和每行像素点的累积能量值,通过下述公式计算每行像素点的能量偏差值:
sub[j]=maxf×j-allpower[j]
其中,sub[j]为第j行像素点的能量偏差值,maxf为能量斜率,allpower[j]为第j行像素点的累积能量值。
5.根据权利要求1所述的方法,其特征在于,所述确定所述感兴趣区域内包含的血管壁的部位,包括:
如果所述一个或多个目标能量偏差值均大于0,则确定所述感兴趣区域内包含的血管壁的部位为前壁部位;
如果所述一个或多个目标能量偏差值均小于0,则确定所述感兴趣区域内包含的血管壁的部位为后壁部位;
如果所述一个或多个目标能量偏差值包括大于0的目标能量偏差值和小于0的目标能量偏差值,则确定所述感兴趣区域内包含的血管壁的部位包括前壁部位和后壁部位。
6.根据权利要求1-3任一所述的方法,其特征在于,所述根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度,包括:
根据所述感兴趣区域内包含的血管壁的部位,在所述感兴趣区域内确定初始能量演化曲线,所述初始能量演化曲线为闭合曲线;
根据第一图像区域内像素点的灰度值和第二图像区域内像素点的灰度值,确定所述血管壁的内膜的内边界,所述第一图像区域是指所述初始能量演化曲线内的图像区域,所述第二图像区域是指所述初始能量演化曲线外的图像区域;
根据所述内膜的内边界和所述感兴趣区域内包含的血管壁的部位,确定所述血管壁的中膜的外边界;
根据所述内膜的内边界和所述中膜的外边界,确定所述内中膜厚度。
7.根据权利要求6所述的方法,其特征在于,所述根据所述内膜的内边界和所述感兴趣区域内包含的血管壁的部位,确定所述血管壁的中膜的外边界,包括:
确定所述感兴趣区域内每个像素点的纵向灰度梯度值;
根据所述感兴趣区域内包含的血管壁的部位,确定查找方向;
按照所述查找方向,查找所述内膜的内边界上每个像素点对应的多个像素点中纵向灰度梯度值最大的像素点;
根据查找到的多个像素点,确定所述中膜的外边界。
8.一种确定内中膜厚度的装置,其特征在于,所述装置包括:
第一获取模块,用于获取待检测图像中的感兴趣区域;
第一确定模块,用于确定所述感兴趣区域内每行像素点的累积能量值;
第二确定模块,用于如果根据所述感兴趣区域内每行像素点的累积能量值确定所述感兴趣区域内包含血管壁,则确定所述感兴趣区域内包含的血管壁的部位,所述血管壁的部位包括前壁部位和/或后壁部位;
第三确定模块,用于根据所述感兴趣区域内包含的血管壁的部位,确定内中膜厚度;
所述装置还包括用于执行以下操作的模块:
从所述感兴趣区域包括的多行像素点的累积能量值中获取最小累积能量值和最大累积能量值;
根据所述最小累积能量值和所述最大累积能量值,确定能量斜率;
根据所述能量斜率和每行像素点的累积能量值,确定每行像素点的能量偏差值;
当确定的多个能量偏差值中包括大于0的多个第一能量偏差值和小于0的多个第二能量偏差值时,如果所述多个第一能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值和/或所述多个第二能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当确定的多个能量偏差值均大于0时,如果所述多个能量偏差值中存在大于第一参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁;
当确定的多个能量偏差值均小于0时,如果所述多个能量偏差值中存在绝对值大于第二参考阈值的一个或多个目标能量偏差值,则确定所述感兴趣区域内包含血管壁。
9.一种计算机可读存储介质,其特征在于,所述存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1-7任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910824455.8A CN110517263B (zh) | 2019-09-02 | 2019-09-02 | 确定内中膜厚度的方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910824455.8A CN110517263B (zh) | 2019-09-02 | 2019-09-02 | 确定内中膜厚度的方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110517263A CN110517263A (zh) | 2019-11-29 |
CN110517263B true CN110517263B (zh) | 2022-05-20 |
Family
ID=68630411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910824455.8A Active CN110517263B (zh) | 2019-09-02 | 2019-09-02 | 确定内中膜厚度的方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110517263B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833757A (zh) * | 2009-03-11 | 2010-09-15 | 深圳迈瑞生物医疗电子股份有限公司 | 血管图像组织结构的边缘检测及血管内膜检测方法和*** |
US20100312110A1 (en) * | 2006-01-11 | 2010-12-09 | Panasonic Corporation | Ultrasonograph |
CN104517277A (zh) * | 2013-09-30 | 2015-04-15 | 中国人民解放军第二军医大学 | 一种颈总动脉超声波纵切图像后处理装置及方法 |
CN104665872A (zh) * | 2014-12-29 | 2015-06-03 | 深圳开立生物医疗科技股份有限公司 | 一种基于超声图像的颈动脉内中膜厚度测量方法和装置 |
CN105380682A (zh) * | 2015-12-24 | 2016-03-09 | 飞依诺科技(苏州)有限公司 | 基于超声图像的双边内中膜的测量方法及*** |
CN106388867A (zh) * | 2016-09-28 | 2017-02-15 | 深圳华声医疗技术有限公司 | 血管的内中膜自动识别测量方法及超声仪 |
CN109919953A (zh) * | 2019-01-21 | 2019-06-21 | 深圳蓝韵医学影像有限公司 | 颈动脉内中膜厚度测量的方法、***和设备 |
CN110047086A (zh) * | 2019-04-24 | 2019-07-23 | 飞依诺科技(苏州)有限公司 | 颈动脉内膜厚度自动测量方法及*** |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102800087B (zh) * | 2012-06-28 | 2014-12-31 | 华中科技大学 | 超声颈动脉血管膜的自动分割方法 |
CN102982314B (zh) * | 2012-11-05 | 2016-05-25 | 深圳市恩普电子技术有限公司 | 一种血管内外膜识别、描记和测量的方法 |
CN106570856A (zh) * | 2016-08-31 | 2017-04-19 | 天津大学 | 水平集动态规划颈总动脉内中膜厚度测量装置和方法 |
CN106570871B (zh) * | 2016-11-02 | 2019-10-18 | 天津大学 | 模糊c均值颈动脉超声图像内中膜厚度测量方法及*** |
-
2019
- 2019-09-02 CN CN201910824455.8A patent/CN110517263B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100312110A1 (en) * | 2006-01-11 | 2010-12-09 | Panasonic Corporation | Ultrasonograph |
CN101833757A (zh) * | 2009-03-11 | 2010-09-15 | 深圳迈瑞生物医疗电子股份有限公司 | 血管图像组织结构的边缘检测及血管内膜检测方法和*** |
CN104517277A (zh) * | 2013-09-30 | 2015-04-15 | 中国人民解放军第二军医大学 | 一种颈总动脉超声波纵切图像后处理装置及方法 |
CN104665872A (zh) * | 2014-12-29 | 2015-06-03 | 深圳开立生物医疗科技股份有限公司 | 一种基于超声图像的颈动脉内中膜厚度测量方法和装置 |
CN105380682A (zh) * | 2015-12-24 | 2016-03-09 | 飞依诺科技(苏州)有限公司 | 基于超声图像的双边内中膜的测量方法及*** |
CN106388867A (zh) * | 2016-09-28 | 2017-02-15 | 深圳华声医疗技术有限公司 | 血管的内中膜自动识别测量方法及超声仪 |
CN109919953A (zh) * | 2019-01-21 | 2019-06-21 | 深圳蓝韵医学影像有限公司 | 颈动脉内中膜厚度测量的方法、***和设备 |
CN110047086A (zh) * | 2019-04-24 | 2019-07-23 | 飞依诺科技(苏州)有限公司 | 颈动脉内膜厚度自动测量方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN110517263A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110210571B (zh) | 图像识别方法、装置、计算机设备及计算机可读存储介质 | |
US11205282B2 (en) | Relocalization method and apparatus in camera pose tracking process and storage medium | |
WO2020221012A1 (zh) | 图像特征点的运动信息确定方法、任务执行方法和设备 | |
CN109815150B (zh) | 应用测试方法、装置、电子设备及存储介质 | |
CN110148178B (zh) | 相机定位方法、装置、终端及存储介质 | |
CN111091576A (zh) | 图像分割方法、装置、设备及存储介质 | |
KR20210111833A (ko) | 타겟의 위치들을 취득하기 위한 방법 및 장치와, 컴퓨터 디바이스 및 저장 매체 | |
CN109302632B (zh) | 获取直播视频画面的方法、装置、终端及存储介质 | |
CN109886208B (zh) | 物体检测的方法、装置、计算机设备及存储介质 | |
CN110570460A (zh) | 目标跟踪方法、装置、计算机设备及计算机可读存储介质 | |
CN113763228A (zh) | 图像处理方法、装置、电子设备及存储介质 | |
CN111104893A (zh) | 目标检测方法、装置、计算机设备及存储介质 | |
CN112749590B (zh) | 目标检测方法、装置、计算机设备和计算机可读存储介质 | |
CN109754439B (zh) | 标定方法、装置、电子设备及介质 | |
CN111784841A (zh) | 重建三维图像的方法、装置、电子设备及介质 | |
CN111598896A (zh) | 图像检测方法、装置、设备及存储介质 | |
CN108304841B (zh) | ***定位方法、装置及存储介质 | |
CN113191976A (zh) | 图像拍摄方法、装置、终端及存储介质 | |
CN110728744B (zh) | 一种体绘制方法、装置及智能设备 | |
CN110517263B (zh) | 确定内中膜厚度的方法、装置及存储介质 | |
CN111127539B (zh) | 视差确定方法、装置、计算机设备及存储介质 | |
CN114241055A (zh) | 一种改进型鱼眼镜头内参标定方法、***、终端及存储介质 | |
CN110443841B (zh) | 地面深度的测量方法、装置及*** | |
CN109472855B (zh) | 一种体绘制方法、装置及智能设备 | |
CN109685881B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |