CN103927559B - 超声图像胎儿颜面部标准切面自动识别方法及系统 - Google Patents
超声图像胎儿颜面部标准切面自动识别方法及系统 Download PDFInfo
- Publication number
- CN103927559B CN103927559B CN201410155814.2A CN201410155814A CN103927559B CN 103927559 B CN103927559 B CN 103927559B CN 201410155814 A CN201410155814 A CN 201410155814A CN 103927559 B CN103927559 B CN 103927559B
- Authority
- CN
- China
- Prior art keywords
- feature
- vector
- gaussian
- rootsift
- module
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000001605 fetal effect Effects 0.000 title abstract description 9
- 208000035126 Facies Diseases 0.000 title abstract 3
- 239000013598 vector Substances 0.000 claims abstract description 64
- 238000000605 extraction Methods 0.000 claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000009826 distribution Methods 0.000 claims description 30
- 238000010606 normalization Methods 0.000 claims description 26
- 238000007781 pre-processing Methods 0.000 claims description 17
- 239000000203 mixture Substances 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 10
- 210000003754 fetus Anatomy 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 5
- 238000002604 ultrasonography Methods 0.000 abstract description 17
- 238000003745 diagnosis Methods 0.000 abstract description 9
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000000513 principal component analysis Methods 0.000 description 16
- 238000013507 mapping Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000012706 support-vector machine Methods 0.000 description 6
- 230000009977 dual effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011524 similarity measure Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000002216 heart Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 210000002307 prostate Anatomy 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种超声图像胎儿颜面部标准切面自动识别方法,包括如下步骤:对原始超声图像进行预处理;对预处理后的超声图像采用RootSIFT进行特征提取;对上述提取的特征利用混合高斯模型进行变换,对变换后的混合高斯模型进行Fisher向量编码,并获得向量编码的直方图的值;对直方图进行归一化以得到特征向量;根据所述特征向量采用基于随机双协调机制的SVM分类器对标准切面进行学习和识别。本发明还涉及一种超声图像胎儿颜面部标准切面自动识别系统。本发明能够提高医生诊断结果的准确性,减少工作量及缩短诊断时间。
Description
技术领域
本发明涉及一种超声图像胎儿颜面部标准切面自动识别方法及系统。
背景技术
超声成像由于低消耗、实时成像及无辐射广泛应用于妊娠诊断中。在超声成像检查中,标准切面的获取是生物测量和诊断的先决条件。在超声诊断的过程中,临床医生最先标注标准切面,通过标准切面检查主要解剖结构,然后再进行进一步的诊断和检查,根据胎儿生长情况做出分析和诊断。
在目前的超声诊断过程中,标准切面的获取过度地依赖于医生的临床经验以及解剖结构知识,标准切面的获取不仅对新手极具挑战性,对有经验的医生也很耗费时间。而且在欠发达国家,有经验和专业知识的医生非常稀缺。因此,精确自动的识别标准切面不仅非常有用,而且适用于有经验和无经验的操作者,这项技术对贫困国家非常有利。
近年来,从超声图像中自动识别标准切面获得了广泛的关注和发展。由于标准切面与非标准切面之间的差异很小,标准切面很高的类内变化主要由于各种妊娠年龄、不同的胎儿姿势和各种扫描方向以及超声图像中斑纹和噪声的存在。目前,标准切面的自动识别仍然非常具有挑战性,也是一个难题。
发明内容
有鉴于此,有必要提供一种超声图像胎儿颜面部标准切面自动识别方法及系统。
本发明提供一种超声图像胎儿颜面部标准切面自动识别方法,该方法包括如下步骤:a.对原始超声图像进行预处理;b.对预处理后的超声图像采用RootSIFT进行特征提取;c.对上述提取的特征利用混合高斯模型进行变换,对变换后的混合高斯模型进行Fisher向量编码,并获得所述向量编码的直方图的值;d.对直方图进行归一化以得到特征向量;e.根据所述特征向量采用基于随机双协调机制的SVM分类器对标准切面进行学习和识别。
其中,所述的预处理包括图像降噪和图像增强。
所述的步骤b包括:对于预处理后的超声图像,手动选择感兴趣区域;对所述感兴趣区域密集采样;用RootSIFT对密集采样后的感兴趣区域进行特征提取。
所述的混合高斯模型基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用Fisher向量编码成一个独立的矢量。
所述的标准切面包括矢状面、冠状面和柱状面。
本发明还提供一种超声图像胎儿颜面部标准切面自动识别系统,包括相互电性连接的预处理模块、特征提取模块、编码模块、归一化模块及识别模块。其中,所述预处理模块用于对原始超声图像进行预处理;所述特征提取模块用于对预处理后的超声图像采用RootSIFT进行特征提取;所述编码模块用于对上述提取的特征利用混合高斯模型进行变换,对变换后的混合高斯模型进行Fisher向量编码,并获得所述向量编码的直方图的值;所述归一化模块用于对直方图进行归一化以得到特征向量;所述识别模块用于根据所述特征向量采用基于随机双协调机制的SVM分类器对标准切面进行学习和识别。
其中,所述的预处理包括图像降噪和图像增强。
所述的特征提取模块具体用于:对于预处理后的超声图像,手动选择感兴趣区域;对所述感兴趣区域密集采样;用RootSIFT对密集采样后的感兴趣区域进行特征提取。
所述的混合高斯模型基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用Fisher向量编码成一个独立的矢量。
所述的标准切面包括矢状面、冠状面和柱状面。
本发明超声图像胎儿颜面部标准切面自动识别方法及系统,能够实现超声图像中胎儿颜面部标准切面的成功识别。在超声图像中,标准切面的自动识别是提高医生诊断效果,减少繁琐工作量以及诊断过程时间的有效方式。此外,本发明也可应用到其他领域的分类和检测中,也可以用于其他器官的标准切面(如:腹部,心脏,前列腺,肺和肝脏)的分类,以及预测和识别癌细胞。
附图说明
图1为本发明超声图像胎儿颜面部标准切面自动识别方法的流程图;
图2为本发明超声图像胎儿颜面部标准切面自动识别系统的硬件架构图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细的说明。
参阅图1所示,是本发明超声图像胎儿颜面部标准切面自动识别方法较佳实施例的作业流程图。
步骤S401,对原始超声图像进行预处理。具体而言,在识别标准切面之前对原始超声图像进行预处理,所述预处理包括图像降噪和图像增强。
步骤S402,对预处理后的超声图像采用RootSIFT进行特征提取。对于预处理后的超声图像,可以手动选择感兴趣区域(Regionof Interest,ROI),以减少搜索范围。所述感兴趣区域包括矢状面区域、冠状面区域和柱状面区域。对所述感兴趣区域密集采样,而后用RootSIFT进行特征提取。将所述超声图像划分成小块,用RootSIFT块述符来代表每个小块。具体步骤如下:
本实施例中,所述特征提取是基于RootSIFT并利用SIFT(Scale-invariantfeature transform,利用尺度不变特征转换)的非线性变换,表示如下:
RootSIFT=sqrt(SIFT/sum(SIFT)).
利用RootSIFT的特征提取类似于在原始的SIFT特征中运用Hellinger(即Hel.)核。RootSIFT优于SIFT,因为欧式距离在RootSIFT向量相当于最初SIFT中用到的Hel.核,或者Chi2距离证明在识别性能上优于欧式距离。
由于欧式距离对距离远的特别敏感,而Hel.距离经常受小距离控制。Hel.距离表示为l1的规一化x和y(n-矢量数)计算如下:
假设||x||2=||y||2=1,欧式距离定义为:
将x用x'替换(x'是元素的平方根),x'是l2的规一化,在特征图谱空间的欧式距离就相当于原始空间的Hel.距离。
x'Ty'=H(x,y).
使用RootSIFT描述符,基于SIFT的欧式距离每一步都可以很容易的进行转换,没有额外的成本、存储和处理空间。
步骤S403,对上述提取的特征利用混合高斯模型(Gaussian Mixture Model,GMM)进行变换,对变换后的混合高斯模型进行Fisher向量(FV)编码,并获得所述向量编码的直方图的值。所述GMM基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用FV编码成一个独立的矢量。具体而言:
本实施例中采用GMM模型,以提高识别性能。用K-means学习得到编码{μk,k=1,...,K},一组局部描述符:{xm,m=1,...,N}。提取特征向量的步骤如下:
邻近分配:
计算vk:
关联vk和规一化所有的特征向量。
在图形表示中,一个固定长度向量vk的维数依赖于参数的数量。为了选择最合适数据的参数,需要关联到更多的数据统计。总之,特征向量对应混合高斯(GMM)模型。
用FV对混合高斯模型的衍生对数概度进行编码。高斯均值和方差的一次和二次求导,密集特征与混合高斯模型(GMM)中心之间的计算如下:
其中,{wk,μk,σk}是GMM的混合权重、均值和对角斜方差,γm(k)是第K个高斯分布的第m个特征的软分配权重。通过将不向的向量的关联到一起:得到FVφ。编码主要的目的,区分一个特定测试图像与所训练图像之间的分布差异。FV是软分布VLAD的高阶统计和BoVW的重要扩展。对于D维特征向量,所述BoVW和FV的主要区别可以表示为:
φBoVW(xm)=[0,...,0,1,0,...,0],
FV的维数比传统的BoVW方法要高,因此用PCA(主成分分析)来降低特征向量的维数以及处理时间。由于不相关特征和GMM协方差对角矩阵假设是一致的,PCA降维的之后也是满足协方差对角矩阵假设的。在本实施例中,每个图像像素共生成128个特征向量。通过PCA,特征维数由128降到64。
步骤S404,对直方图进行归一化以得到特征向量,也即,用直方图表示FV聚类的K-means算法。具体步骤如下:
考虑到训练数据集的多样性,对特征值进行归一化以提高识别是非常有效的。选取合适的lp归一化特征值对提高分类结果非常有帮助,这样可以有效地去除背景信息。每个矩阵除以相应的规一化值,lp归一化后特征矩阵分布在半径为1的球面上。lp幅值用实数表示,当p≥1时定义如下:
其中,p=1为l1归一化,p=2是l2也就是欧式距离归一化。
假设一个数据集,ui~U(p,n)遵循p分布,在n维空间的单位球面上使用i.i.d.模型和特征矩阵的分布,封闭形式来解规范化坐标:表示为:
其中,Γ(·)是伽马函数,上式中当p=2时符合高斯分布。
如果数据是一个广义的高斯分布,用lp归一化是有效方法。lp幅值可以定义为寻求一个最大熵的pdf.
其中,(σp)p=constant常量和(σp)p=1.为了减少交叉空间维数的差异,特征缩放比例非高斯级用类似高斯的特征来表示。l2可以消除与内容无关的背景信息。
一般来说,一些突发的能量值会影响由SVM分类器的相似性得分,所述突发的能量值可能导致次优相似性得分。因此,在特征矩阵中的较大值相应应该得到抑制。一个简单而有效的方法,如功率归一化定义为:
z←sign(z)|z|ρ,0≤ρ≤1.
在本实施例中,ρ设置为:ρ=0.5,它表示平方根或者更简单的平方根,进行功率规一化可以减少稀疏化数据的突发性特征的影响。除以之外,功率规一化被认为是Hel.核的明显数据显示。值得注意的是,由平方根变换的功率归一化是不特定并对标准切面的识别有利。
步骤S405,采用随机双协调(SDCA)机制对标准切面进行学习和识别,以提升SVM分类器进行分类。一对多的SVM分类器基于SDCA,根据余弦相似性度量对超声图像进行分类。具体步骤如下:
支持向量机已应用广泛的识别算法,使用统计学习理论和结构风险最小化原则找到全局最优解。SVM最主要的优势是可以处理大尺寸数据,用于解决识别问题。一对多的评分方案来识别超声图像中的不同切面。评分函数在SVM分类超平面H定义为:
H:wTxi+b=0,i=1,2,...,n,
其中,x1,x2,...,xn是输入的向量RD维内,b∈R是偏置参数,w1,w2,...,wn是向量的权重,T表示转置运算符。支持向量机的主要目的是获取最优w1,w2,...,wn值,SDCA是探索,因为其能够通过高精确度双重目标获得这些最佳值。同时,SVM中的目标函数通过不用损失函数的SDCA进行最小化。标签y1,...,yn属于{±1},支持向量机问题的线性内核和无偏差项定义为:
Ψi(α)=max{0,1-yiα}.
这个问题可以转化为解决最小其中:
上述Ψi(α)=max{0,1-yiα}.的对偶问题的由SDCA解决:
其中
令初始化w(0)=w(α(0)),目标函数每次迭代t获得:
更新规则如下:
α(t)←α(t-1)+Δαiei.
w(t)←w(t-1)+(λn)-1Δαixi.
最后,输出结果通过SDCA方法输出的α和w的平均值获得。
本实施例中,PCA首先应用于学习线性投影:W∈RD'×D,D'<<D,从高维到低维:PCA映射后,如果i和j来自相同的图像,则两幅图像之间的欧式距离的平方值或者学习后的阈值:会变小;反之,则会变大。SVM问题可通过下面的约束解决:
公式1
其中,当且仅当图像来自同一类时,yi,j=1。这种投影可以被视为一个在原始空间中的低秩度量:
公式2
其中,WTW是马氏矩阵。由于分解,这个矩阵的秩等于D',小于满秩D。因此,学习映射矩阵W就和学习低秩矩阵WTW一样。将公式1中的约束方程结合,化简后的方程可表示为:
公式3
初始化后,映射矩阵将满足如下条件:
其中,是不同向量的外积,μ是学习比例参数。公式3中的目标函数在W中不是一个凸矩阵,初始化W是非常关键的。一个实际的案例就是提取D'的最大PCA维数。此外,如果不使用PCA标准,主特征值的明显度将被降低,因为多样性越差的模型越容易被忽略。
参阅图2所示,是本发明超声图像胎儿颜面部标准切面自动识别系统的硬件架构图。该系统包括相互电性连接的预处理模块、特征提取模块、编码模块、归一化模块及识别模块。
所述预处理模块用于对原始超声图像进行预处理。具体而言,在识别标准切面之前对原始超声图像进行预处理,所述预处理包括图像降噪和图像增强。
所述特征提取模块用于对预处理后的超声图像采用RootSIFT进行特征提取。对于预处理后的超声图像,可以手动选择感兴趣区域(Regionof Interest,ROI),以减少搜索范围。所述感兴趣区域包括矢状面区域、冠状面区域和柱状面区域。对所述感兴趣区域密集采样,而后用RootSIFT进行特征提取。将所述超声图像划分成小块,用RootSIFT块述符来代表每个小块。具体如下:
本实施例中,所述特征提取是基于RootSIFT并利用SIFT(Scale-invariantfeature transform,利用尺度不变特征转换)的非线性变换,表示如下:
RootSIFT=sqrt(SIFT/sum(SIFT)).
利用RootSIFT的特征提取类似于在原始的SIFT特征中运用Hellinger(即Hel.)核。RootSIFT优于SIFT,因为欧式距离在RootSIFT向量相当于最初SIFT中用到的Hel.核,或者Chi2距离证明在识别性能上优于欧式距离。
由于欧式距离对距离远的特别敏感,而Hel.距离经常受小距离控制。Hel.距离表示为l1的规一化x和y(n-矢量数)计算如下:
假设||x||2=||y||2=1,欧式距离定义为:
将x用x'替换(x'是元素的平方根),x'是l2的规一化,在特征图谱空间的欧式距离就相当于原始空间的Hel.距离。
x'Ty'=H(x,y).
使用RootSIFT描述符,基于SIFT的欧式距离每一步都可以很容易的进行转换,没有额外的成本、存储和处理空间。
所述编码模块用于对上述提取的特征利用混合高斯模型(Gaussian MixtureModel,GMM)进行变换,对变换后的混合高斯模型进行Fisher向量(FV)编码,并获得所述向量编码的直方图的值。所述GMM基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用FV编码成一个独立的矢量。具体而言:
本实施例中采用GMM模型,以提高识别性能。用K-means学习得到编码{μk,k=1,...,K},一组局部描述符:{xm,m=1,...,N}。提取特征向量的步骤如下:
邻近分配:
计算vk:
关联vk和规一化所有的特征向量。
在图形表示中,一个固定长度向量vk的维数依赖于参数的数量。为了选择最合适数据的参数,需要关联到更多的数据统计。总之,特征向量对应混合高斯(GMM)模型。
用FV对混合高斯模型的衍生对数概度进行编码。高斯均值和方差的一次和二次求导,密集特征与混合高斯模型(GMM)中心之间的计算如下:
其中,{wk,μk,σk}是GMM的混合权重、均值和对角斜方差,γm(k)是第K个高斯分布的第m个特征的软分配权重。通过将不向的向量的关联到一起:得到FVφ。编码主要的目的,区分一个特定测试图像与所训练图像之间的分布差异。FV是软分布VLAD的高阶统计和BoVW的重要扩展。对于D维特征向量,所述BoVW和FV的主要区别可以表示为:
φBoVW(xm)=[0,...,0,1,0,...,0],
FV的维数比传统的BoVW方法要高,因此用PCA(主成分分析)来降低特征向量的维数以及处理时间。由于不相关特征和GMM协方差对角矩阵假设是一致的,PCA降维的之后也是满足协方差对角矩阵假设的。在本实施例中,每个图像像素共生成128个特征向量。通过PCA,特征维数由128降到64。
所述归一化模块用于对直方图进行归一化以得到特征向量,也即,用直方图表示FV聚类的K-means算法。具体步骤如下:
考虑到训练数据集的多样性,对特征值进行归一化以提高识别是非常有效的。选取合适的lp归一化特征值对提高分类结果非常有帮助,这样可以有效地去除背景信息。每个矩阵除以相应的规一化值,lp归一化后特征矩阵分布在半径为1的球面上。lp幅值用实数表示,当p≥1时定义如下:
其中,p=1为l1归一化,p=2是l2也就是欧式距离归一化。
假设一个数据集,ui~U(p,n)遵循p分布,在n维空间的单位球面上使用i.i.d.模型和特征矩阵的分布,封闭形式来解规范化坐标:表示为:
其中,Γ(·)是伽马函数,上式中当p=2时符合高斯分布。
如果数据是一个广义的高斯分布,用lp归一化是有效方法。lp幅值可以定义为寻求一个最大熵的pdf.
其中,(σp)p=constant常量和(σp)p=1.为了减少交叉空间维数的差异,特征缩放比例非高斯级用类似高斯的特征来表示。l2可以消除与内容无关的背景信息。
一般来说,一些突发的能量值会影响由SVM分类器的相似性得分,所述突发的能量值可能导致次优相似性得分。因此,在特征矩阵中的较大值相应应该得到抑制。一个简单而有效的方法,如功率归一化定义为:
z←sign(z)|z|ρ,0≤ρ≤1.
在本实施例中,ρ设置为:ρ=0.5,它表示平方根或者更简单的平方根,进行功率规一化可以减少稀疏化数据的突发性特征的影响。除以之外,功率规一化被认为是Hel.核的明显数据显示。值得注意的是,由平方根变换的功率归一化是不特定并对标准切面的识别有利。
所述识别模块用于采用随机双协调(SDCA)机制对标准切面进行学习和识别,以提升SVM分类器进行分类。一对多的SVM分类器基于SDCA,根据余弦相似性度量对超声图像进行分类。具体步骤如下:
支持向量机已应用广泛的识别算法,使用统计学习理论和结构风险最小化原则找到全局最优解。SVM最主要的优势是可以处理大尺寸数据,用于解决识别问题。一对多的评分方案来识别超声图像中的不同切面。评分函数在SVM分类超平面H定义为:
H:wTxi+b=0,i=1,2,...,n,
其中,x1,x2,...,xn是输入的向量RD维内,b∈R是偏置参数,w1,w2,...,wn是向量的权重,T表示转置运算符。支持向量机的主要目的是获取最优w1,w2,...,wn值,SDCA是探索,因为其能够通过高精确度双重目标获得这些最佳值。同时,SVM中的目标函数通过不用损失函数的SDCA进行最小化。标签y1,...,yn属于{±1},支持向量机问题的线性内核和无偏差项定义为:
Ψi(α)=max{0,1-yiα}.
这个问题可以转化为解决最小其中:
上述Ψi(α)=max{0,1-yiα}.的对偶问题的由SDCA解决:
其中
令初始化w(0)=w(α(0)),目标函数每次迭代t获得:
更新规则如下:
α(t)←α(t-1)+Δαiei.
w(t)←w(t-1)+(λn)-1Δαixi.
最后,输出结果通过SDCA方法输出的α和w的平均值获得。
本实施例中,PCA首先应用于学习线性投影:W∈RD'×D,D'<<D,从高维到低维:PCA映射后,如果i和j来自相同的图像,则两幅图像之间的欧式距离的平方值或者学习后的阈值:会变小;反之,则会变大。SVM问题可通过下面的约束解决:
公式1
其中,当且仅当图像来自同一类时,yi,j=1。这种投影可以被视为一个在原始空间中的低秩度量:
公式2
其中,WTW是马氏矩阵。由于分解,这个矩阵的秩等于D',小于满秩D。因此,学习映射矩阵W就和学习低秩矩阵WTW一样。将公式1中的约束方程结合,化简后的方程可表示为:
公式3
初始化后,映射矩阵将满足如下条件:
其中,是不同向量的外积,μ是学习比例参数。公式3中的目标函数在W中不是一个凸矩阵,初始化W是非常关键的。一个实际的案例就是提取D'的最大PCA维数。此外,如果不使用PCA标准,主特征值的明显度将被降低,因为多样性越差的模型越容易被忽略。
虽然本发明参照当前的较佳实施方式进行了描述,但本领域的技术人员应能理解,上述较佳实施方式仅用来说明本发明,并非用来限定本发明的保护范围,任何在本发明的精神和原则范围之内,所做的任何修饰、等效替换、改进等,均应包含在本发明的权利保护范围之内。
Claims (8)
1.一种超声图像胎儿颜面部标准切面自动识别方法,其特征在于,该方法包括如下步骤:
a.对原始超声图像进行预处理;
b.对预处理后的超声图像采用RootSIFT进行特征提取;
c.对上述提取的特征利用混合高斯模型进行变换,对变换后的混合高斯模型进行Fisher向量编码,并获得所述向量编码的直方图的值;
d.对直方图进行归一化以得到特征向量;
e.根据所述特征向量采用基于随机双协调机制的SVM分类器对标准切面进行学习和识别;
所述的步骤b包括:
对于预处理后的超声图像,手动选择感兴趣区域;
对所述感兴趣区域密集采样;
用RootSIFT对密集采样后的感兴趣区域进行特征提取;
所述特征提取是基于RootSIFT并利用SIFT(Scale-invariant feature transform,利用尺度不变特征转换)的非线性变换,表示如下:
RootSIFT=sqrt(SIFT/sum(SIFT)),
每个特征矩阵除以相应的规一化值,lp归一化后特征矩阵分布在半径为1的球面上,lp幅值用实数表示,当p≥1时定义如下:
其中,p=1为l1归一化,p=2是l2也就是欧式距离归一化,x为n-矢量数;
假设一个数据集,ui~U(p,n)遵循p分布,在n维空间的单位球面上 使用i.i.d.模型和特征矩阵的分布,封闭形式来解规范化坐标:表示为:
其中,Γ(·)是伽马函数,上式中当p=2时符合高斯分布;
如果数据是一个广义的高斯分布,lp幅值定义为寻求一个最大熵的pdf,
其中,(σp)p=constant常量和(σp)p=1,为了减少交叉空间维数的差异,特征缩放比例非高斯级用高斯的特征来表示;
功率归一化定义为:
z←sign(z)|z|ρ,0≤ρ≤1,
ρ设置为:ρ=0.5,它表示平方根或者更简单的平方根。
2.如权利要求1所述的方法,其特征在于,所述的预处理包括图像降噪和图像增强。
3.如权利要求1所述的方法,其特征在于,所述的混合高斯模型基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用Fisher向量编码成一个独立的矢量。
4.如权利要求1所述的方法,其特征在于,所述的标准切面包括矢状面、冠状面和柱状面。
5.一种超声图像胎儿颜面部标准切面自动识别系统,其特征在于,该系统包括相互电性连接的预处理模块、特征提取模块、编码模块、归一化模块及识别模块,其中:
所述预处理模块用于对原始超声图像进行预处理;
所述特征提取模块用于对预处理后的超声图像采用RootSIFT进行特征提取;
所述编码模块用于对上述提取的特征利用混合高斯模型进行变换,对变换后的混合高斯模型进行Fisher向量编码,并获得所述向量编码的直方图的值;
所述归一化模块用于对直方图进行归一化以得到特征向量;
所述识别模块用于根据所述特征向量采用基于随机双协调机制的SVM分类器对标准切面进行学习和识别;
所述特征提取模块具体用于:
对于预处理后的超声图像,手动选择感兴趣区域;
对所述感兴趣区域密集采样;
用RootSIFT对密集采样后的感兴趣区域进行特征提取;
所述特征提取是基于RootSIFT并利用SIFT(Scale-invariant feature transform,利用尺度不变特征转换)的非线性变换,表示如下:
RootSIFT=sqrt(SIFT/sum(SIFT)),
所述归一化模块在进行归一化时,每个特征矩阵除以相应的规一化值,lp归一化后特征矩阵分布在半径为1的球面上,lp幅值用实数表示,当p≥1时定义如下:
其中,p=1为l1归一化,p=2是l2也就是欧式距离归一化,x为n-矢量数;
假设一个数据集,ui~U(p,n)遵循p分布,在n维空间的单位球面上使用i.i.d.模型和特征矩阵的分布,封闭形式来解规范化坐标:表示为:
其中,Γ(·)是伽马函数,上式中当p=2时符合高斯分布;
如果数据是一个广义的高斯分布,lp幅值定义为寻求一个最大熵的pdf,
其中,(σp)p=constant常量和(σp)p=1,为了减少交叉空间维数的差异,特征缩放比例非高斯级用高斯的特征来表示;
功率归一化定义为:
z←sign(z)|z|ρ,0≤ρ≤1,
ρ设置为:ρ=0.5,它表示平方根或者更简单的平方根。
6.如权利要求5所述的系统,其特征在于,所述的预处理包括图像降噪和图像增强。
7.如权利要求5所述的系统,其特征在于,所述的混合高斯模型基于对角协方差矩阵假设产生K个高斯分布,对一组局部描述符用Fisher向量编码成一个独立的矢量。
8.如权利要求5所述的系统,其特征在于,所述的标准切面包括矢状面、冠状面和柱状面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410155814.2A CN103927559B (zh) | 2014-04-17 | 2014-04-17 | 超声图像胎儿颜面部标准切面自动识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410155814.2A CN103927559B (zh) | 2014-04-17 | 2014-04-17 | 超声图像胎儿颜面部标准切面自动识别方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103927559A CN103927559A (zh) | 2014-07-16 |
CN103927559B true CN103927559B (zh) | 2017-06-16 |
Family
ID=51145775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410155814.2A Active CN103927559B (zh) | 2014-04-17 | 2014-04-17 | 超声图像胎儿颜面部标准切面自动识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103927559B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989094B (zh) * | 2015-02-12 | 2020-09-01 | 中国科学院西安光学精密机械研究所 | 基于隐层语义中层表达的图像检索方法 |
CN104881651B (zh) * | 2015-05-29 | 2017-12-15 | 南京信息工程大学 | 一种基于随机投影和Fisher向量的人物行为识别方法 |
CN106056159B (zh) * | 2016-06-03 | 2019-03-26 | 西安电子科技大学 | 基于Fisher Vector的图像精细分类方法 |
CN106504226B (zh) * | 2016-09-26 | 2019-07-19 | 深圳大学 | 超声图像膀胱脱垂自动分级系统 |
US11534133B2 (en) | 2017-04-27 | 2022-12-27 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Ultrasonic detection method and ultrasonic imaging system for fetal heart |
CN111383207B (zh) * | 2018-12-11 | 2024-03-15 | 深圳开立生物医疗科技股份有限公司 | 一种肌骨超声图像处理方法、系统、装置及可读存储介质 |
CN110335624A (zh) * | 2019-07-29 | 2019-10-15 | 吉林大学 | 基于功率归一化倒谱系数特征的帕金森病语音检测方法 |
CN110570409B (zh) * | 2019-09-05 | 2022-02-18 | 南方医科大学 | 基于混合卷积网络获取胎儿四腔心切面心动周期视频方法 |
CN110464380B (zh) * | 2019-09-12 | 2021-10-29 | 李肯立 | 一种对中晚孕期胎儿的超声切面图像进行质量控制的方法 |
CN110652317B (zh) * | 2019-09-24 | 2020-12-29 | 深圳度影医疗科技有限公司 | 一种产前胎儿超声容积图像中标准切面的自动定位方法 |
CN111210909A (zh) * | 2020-01-13 | 2020-05-29 | 青岛大学附属医院 | 一种基于深度神经网络的直肠癌t分期自动诊断系统及其构建方法 |
CN111310851B (zh) * | 2020-03-03 | 2023-04-28 | 四川大学华西第二医院 | 一种人工智能超声辅助系统及其应用 |
CN111428713B (zh) * | 2020-03-20 | 2023-04-07 | 华侨大学 | 一种基于特征融合的超声图像自动分类法 |
CN111428778B (zh) * | 2020-03-20 | 2023-04-07 | 华侨大学 | 一种胎儿颜面部标准切面自动分类法 |
CN111949003B (zh) * | 2020-07-17 | 2021-09-03 | 浙江浙能技术研究院有限公司 | 一种基于SFA与Hellinger距离的闭环控制回路性能评价方法 |
CN112155603B (zh) * | 2020-09-24 | 2023-06-09 | 广州爱孕记信息科技有限公司 | 胎儿结构特征的权重值确定方法及装置 |
CN112215829B (zh) * | 2020-10-21 | 2021-12-14 | 深圳度影医疗科技有限公司 | 一种髋关节标准切面的定位方法和计算机设备 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1966761A2 (en) * | 2005-12-22 | 2008-09-10 | Philips Intellectual Property & Standards GmbH | Adaptive point-based elastic image registration |
US20110257505A1 (en) * | 2010-04-20 | 2011-10-20 | Suri Jasjit S | Atheromatic?: imaging based symptomatic classification and cardiovascular stroke index estimation |
US20120220875A1 (en) * | 2010-04-20 | 2012-08-30 | Suri Jasjit S | Mobile Architecture Using Cloud for Hashimoto's Thyroiditis Disease Classification |
-
2014
- 2014-04-17 CN CN201410155814.2A patent/CN103927559B/zh active Active
Non-Patent Citations (1)
Title |
---|
智能超声扫查与细胞学筛查;张灵;《中国博士学位论文全文数据库 医药卫生科技辑》;20130815(第08(2013)期);E060-4 * |
Also Published As
Publication number | Publication date |
---|---|
CN103927559A (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103927559B (zh) | 超声图像胎儿颜面部标准切面自动识别方法及系统 | |
Maraci et al. | A framework for analysis of linear ultrasound videos to detect fetal presentation and heartbeat | |
WO2016091017A1 (zh) | 一种高光谱图像分类中光谱向量互相关特征的抽取方法 | |
Sridar et al. | Decision fusion-based fetal ultrasound image plane classification using convolutional neural networks | |
CN110338844B (zh) | 三维成像数据的显示处理方法和三维超声成像方法及系统 | |
Si et al. | Learning hybrid image templates (hit) by information projection | |
US9076197B2 (en) | Probability density function estimation | |
US7995820B2 (en) | System and method for detection of fetal anatomies from ultrasound images using a constrained probabilistic boosting tree | |
US10762630B2 (en) | System and method for structures detection and multi-class image categorization in medical imaging | |
CN107527023B (zh) | 基于超像素和主题模型的极化sar图像分类方法 | |
CN101551855B (zh) | 自适应核匹配追踪辅助诊断系统及其辅助诊断方法 | |
CN103413146B (zh) | 基于Freeman熵和自学习的极化SAR图像精细分类方法 | |
Lei et al. | Automatic recognition of fetal facial standard plane in ultrasound image via fisher vector | |
CN111583320B (zh) | 融合深度卷积网络和影像组学特征的乳腺癌超声图分型方法、系统及存储介质 | |
CN106096517A (zh) | 一种基于低秩矩阵与特征脸的人脸识别方法 | |
CN112767355A (zh) | 一种甲状腺结节Tirads分级自动识别模型构建方法及装置 | |
Lei et al. | Automatic recognition of fetal standard plane in ultrasound image | |
Li et al. | Human sperm health diagnosis with principal component analysis and K-nearest neighbor algorithm | |
CN104123563B (zh) | 一种基于Cloude分解的极化SAR图像分类方法 | |
Lanfredi et al. | Comparing radiologists' gaze and saliency maps generated by interpretability methods for chest x-rays | |
CN117275080B (zh) | 基于计算机视觉的眼部状态识别方法及系统 | |
CN114373079A (zh) | 一种快速准确的探地雷达目标检测方法 | |
Yu et al. | Establishment of computer-aided diagnosis system for liver tumor CT based on SVM | |
Li et al. | Multi-scale feature based medical image classification | |
Zhang et al. | Detection and retrieval of cysts in joint ultrasound b-mode and elasticity breast images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210203 Address after: 518000 5/6 Pingshan Private Enterprise Science Park, Taoyuan Street, Nanshan District, Shenzhen City, Guangdong Province Patentee after: SHENZHEN WISONIC MEDICAL TECHNOLOGY Co.,Ltd. Address before: 518000 No. 3688 Nanhai Road, Shenzhen, Guangdong, Nanshan District Patentee before: SHENZHEN University |
|
TR01 | Transfer of patent right |