CN112365411A - 一种基于位置点匹配的手腕骨兴趣区域修复方法 - Google Patents

一种基于位置点匹配的手腕骨兴趣区域修复方法 Download PDF

Info

Publication number
CN112365411A
CN112365411A CN202011124718.3A CN202011124718A CN112365411A CN 112365411 A CN112365411 A CN 112365411A CN 202011124718 A CN202011124718 A CN 202011124718A CN 112365411 A CN112365411 A CN 112365411A
Authority
CN
China
Prior art keywords
bone
matching
roi
coordinates
bones
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
CN202011124718.3A
Other languages
English (en)
Other versions
CN112365411B (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202011124718.3A priority Critical patent/CN112365411B/zh
Publication of CN112365411A publication Critical patent/CN112365411A/zh
Application granted granted Critical
Publication of CN112365411B publication Critical patent/CN112365411B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/77Retouching; Inpainting; Scratch removal
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30008Bone

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)

Abstract

一种基于位置点匹配的手腕骨兴趣区域修复方法,包括:根据人手腕图像以及对应的参照骨ROI标注信息,准备匹配数据集;计算参照骨ROI中心点坐标,作为匹配位置点;根据位置点坐标,计算样本与匹配数据集中所有骨骼区域的匹配距离;将数据集中匹配距离最短的骨骼区域作为目标参照骨ROI,并填补参照骨ROI集合。本发明通过参照骨的位置点大数据样本,自动匹配并填补参照骨ROI集合中缺失的参照骨类型,准确率高、泛化能力强;本方法可以完整且准确地提取出人手腕中所有参照骨ROI,具备高度扩展性、可操作性。

Description

一种基于位置点匹配的手腕骨兴趣区域修复方法
技术领域
本发明涉及一种基于位置点匹配的手腕骨兴趣区域修复方法。
背景技术
骨龄也称骨骼年龄,用于表示骨骼的生长发育程度,是评价个体发育程度的重要指标之一。骨龄评估在众多领域具有广泛的应用,如鉴定嫌疑人或受害者的生理年龄、运动员选材和儿童成年后身高预测等。
骨龄评估的方法有很多种,现在国际上比较常用的有图谱法和计分法。图谱法将手骨图片与标准手骨图谱比较,选择最接近的标准手骨图谱的骨龄作为该手骨图片的骨龄;计分法则根据判定每块参照骨的成熟指征,通过成熟指征组合计算出骨龄。其中,完整地提取出手腕参照骨的兴趣区域(Region of Interest,ROI)是保证骨龄精确评估的关键。
将深度学习用于手腕骨兴趣区域的提取,已经越来越成为发展趋势。然而现有的方法大多数基于传统图像处理方法,如形态学、边缘处理、区域生长等,操作复杂度较高,复现起来比较困难。而且普通的卷积神经网络在参照骨的兴趣区域提取中,存在部分参照骨漏判的情况,导致参照骨ROI数量不完整,极大地影响了后续骨龄评估的准确率。
发明内容
本发明要克服现有技术的上述缺点,提供一种基于位置点匹配的手腕骨兴趣区域修复方法。
本发明包括:根据人手腕图像以及对应的参照骨ROI标注信息,准备匹配数据集;计算参照骨ROI中心点坐标,作为匹配位置点;根据位置点坐标,计算样本与匹配数据集中所有骨骼区域的匹配距离;将数据集中匹配距离最短的骨骼区域作为目标参照骨ROI,并填补参照骨ROI集合。本发明通过参照骨的位置点大数据样本,自动匹配并填补参照骨ROI集合中缺失的参照骨类型,准确率高、泛化能力强;本方法可以完整且准确地提取出人手腕中所有参照骨ROI,具备高度扩展性、可操作性。
本发明解决技术问题采用如下技术方案:
一种基于位置点匹配的手腕骨兴趣区域修复方法,包括以下步骤:
步骤1:准备匹配数据集;
步骤2:计算参照骨ROI位置点坐标;
步骤3:计算参照骨ROI匹配距离;
步骤4:填补参照骨ROI。
作为本发明的进一步优化,所述的步骤1具体包括:
1)确定骨龄评估计分法中要提取的所有参照骨的类型BON;
2)将所有参照骨按骨骼位置划分为4个骨骼区域,分别为指骨Ⅰ、指骨Ⅲ、指骨Ⅴ、腕骨区域,记为ZHO,ZHT,ZHF,WAN。每个骨骼区域内含若干不同参照骨类型;
BON=ZHO∪ZHT∪ZHF∪WAN (1-1)
REG={ZHO,ZHT,ZHF,WAN} (1-2)
3)定义参照骨顶点坐标集合为B,B中第n块参照骨坐标为bn,n∈BON,其左上角顶点坐标记为(xtop,ytop),右下角顶点坐标记为(xbot,ybot),则bn的顶点坐标表示为:
bn=(xtop,ytop,xbot,ybot) (1-3)
4)确定匹配数据集BONES,BONES={Bi|1≤i≤W},W为数据集的数量;
作为本发明的更进一步优化,所述的步骤2具体包括:
1)给定数据样本t,t中包含的参照骨集合为T,
Figure BDA0002733215740000021
其顶点坐标集合为Bt
Figure BDA0002733215740000031
2)确定t中缺失的参照骨集合U:
Figure BDA0002733215740000032
3)对于U中任意缺失参照骨u,确定其所属的骨骼区域Z,Z∈REG。对于Z区域内所有参照骨顶点坐标bn,计算其位置点坐标dn
Figure BDA0002733215740000033
4)将顶点坐标集合Bt,匹配数据集BONES,分别转换为位置点坐标集合Mt,位置点匹配坐标集合REGIONS;
作为本发明的更进一步优化,所述的步骤3具体包括:
1)以Mt中第一块参照骨的类型r为基准,r∈Z。分别计算Mt与REGIONS中,以r类型参照骨dr为参照原点的相对位置坐标fn,公式为:
fn=Relative(dn,dr) (3-1)
Relative(x,y)=(x[0]-y[0],x[1]-y[1]) (3-2)
2)使用Similar函数计算两个骨骼区域之间的匹配距离,公式为:
Figure BDA0002733215740000034
Figure BDA0002733215740000035
其中,
Figure BDA0002733215740000037
表示骨骼区域x内,参照骨的类型为n的相对位置坐标;
作为本发明的更进一步优化,所述的步骤4具体包括:
1)确定REGIONS中与Mt匹配距离最短的骨骼区域target,公式为:
Figure BDA0002733215740000036
2)根据target在REGIONS中的索引INDEX,填补Bt中缺失的参照骨ROI。
Bt=Bt∪BONES[INDEX][u] (4-2)
本发明具有如下有益效果:
(1)通过参照骨的位置点大数据样本,自动匹配并填补参照骨ROI集合中缺失的参照骨类型,准确率高、泛化能力强;
(2)本方法可以完整且准确地提取出人手腕中所有参照骨ROI,具备高度扩展性、可操作性。
附图说明
图1是本发明方法的骨骼区域划分图。
图2是本发明方法的位置点坐标转换流程图。
具体实施方式
下面结合附图进一步说明本发明的技术方案。
一种基于位置点匹配的手腕骨兴趣区域修复方法,其特征在于,包括以下步骤:
步骤1:准备匹配数据集;具体包括:
11)确定骨龄评估计分法中要提取的所有参照骨的类型BON;
12)将所有参照骨按骨骼位置划分为4个骨骼区域,如图1所示,分别为指骨Ⅰ、指骨Ⅲ、指骨Ⅴ、腕骨区域,记为ZHO,ZHT,ZHF,WAN。每个骨骼区域内含若干不同参照骨类型;
BON=ZHO∪ZHT∪ZHF∪WAN (1-1)
REG={ZHO,ZHT,ZHF,WAN} (1-2)
13)定义参照骨顶点坐标集合为B,B中第n块参照骨坐标为bn,n∈BON,其左上角顶点坐标记为(xtop,ytop),右下角顶点坐标记为(xbot,ybot),则bn的顶点坐标表示为:
bn=(xtop,ytop,xbot,ybot) (1-3)
14)确定匹配数据集BONES,BONES={Bi|1≤i≤W},W为数据集的数量;
步骤2:计算参照骨ROI位置点坐标;具体包括:
21)给定数据样本t,t中包含的参照骨集合为T,
Figure BDA0002733215740000051
其顶点坐标集合为Bt
Figure BDA0002733215740000052
22)确定t中缺失的参照骨集合U:
Figure BDA0002733215740000053
23)对于U中任意缺失参照骨u,确定其所属的骨骼区域Z,Z∈REG。对于Z区域内所有参照骨顶点坐标bn,计算其位置点坐标dn
Figure BDA0002733215740000054
24)将顶点坐标集合Bt,匹配数据集BONES,分别转换为位置点坐标集合Mt,位置点匹配坐标集合REGIONS,如图2所示;
步骤3:计算参照骨ROI匹配距离;具体包括:
31)以Mt中第一块参照骨的类型r为基准,r∈Z。分别计算Mt与REGIONS中,以r类型参照骨dr为参照原点的相对位置坐标fn,公式为:
fn=Relative(dn,dr) (3-1)
Relative(x,y)=(x[0]-y[0],x[1]-y[1]) (3-2)
32)使用Similar函数计算两个骨骼区域之间的匹配距离,公式为:
Figure BDA0002733215740000055
Figure BDA0002733215740000056
其中,
Figure BDA0002733215740000057
表示骨骼区域x内,参照骨的类型为n的相对位置坐标;
步骤4:填补参照骨ROI;具体包括:
41)确定REGIONS中与Mt匹配距离最短的骨骼区域target,公式为:
Figure BDA0002733215740000061
42)根据target在REGIONS中的索引INDEX,填补Bt中缺失的参照骨ROI。
Bt=Bt∪BONES[INDEX][u] (4-2)
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (5)

1.一种基于位置点匹配的手腕骨兴趣区域修复方法,包括以下步骤:
步骤1:准备匹配数据集;
步骤2:计算参照骨ROI位置点坐标;
步骤3:计算参照骨ROI匹配距离;
步骤4:填补参照骨ROI。
2.根据权利要求1所述的一种基于位置点匹配的手腕骨兴趣区域修复方法,其特征在于:步骤1具体包括:
11)确定骨龄评估计分法中要提取的所有参照骨的类型BON;
12)将所有参照骨按骨骼位置划分为4个骨骼区域,分别为指骨Ⅰ、指骨Ⅲ、指骨Ⅴ、腕骨区域,记为ZHO,ZHT,ZHF,WAN。每个骨骼区域内含若干不同参照骨类型;
BON=ZHO∪ZHT∪ZHF∪WAN (1-1)
REG={ZHO,ZHT,ZHF,WAN} (1-2)
13)定义参照骨顶点坐标集合为B,B中第n块参照骨坐标为bn,n∈BON,其左上角顶点坐标记为(xtop,ytop),右下角顶点坐标记为(xbot,ybot),则bn的顶点坐标表示为:
bn=(xtop,ytop,xbot,ybot) (1-3)
14)确定匹配数据集BONES,BONES={Bi|1≤i≤W},W为数据集的数量。
3.根据权利要求1所述的一种基于位置点匹配的手腕骨兴趣区域修复方法,其特征在于:步骤2具体包括:
21)给定数据样本t,t中包含的参照骨集合为T,
Figure FDA0002733215730000011
其顶点坐标集合为Bt
Figure FDA0002733215730000012
22)确定t中缺失的参照骨集合U:
Figure FDA0002733215730000021
23)对于U中任意缺失参照骨u,确定其所属的骨骼区域Z,Z∈REG。对于Z区域内所有参照骨顶点坐标bn,计算其位置点坐标dn
Figure FDA0002733215730000022
24)将顶点坐标集合Bt,匹配数据集BONES,分别转换为位置点坐标集合Mt,位置点匹配坐标集合REGIONS。
4.根据权利要求1所述的一种基于位置点匹配的手腕骨兴趣区域修复方法,其特征在于:步骤3具体包括:
31)以Mt中第一块参照骨的类型r为基准,r∈Z。分别计算Mt与REGIONS中,以r类型参照骨dr为参照原点的相对位置坐标fn,公式为:
fn=Relative(dn,dr) (3-1)
Relative(x,y)=(x[0]-y[0],x[1]-y[1]) (3-2)
32)使用Similar函数计算两个骨骼区域之间的匹配距离,公式为:
Figure FDA0002733215730000023
Figure FDA0002733215730000024
其中,
Figure FDA0002733215730000025
表示骨骼区域x内,参照骨的类型为n的相对位置坐标。
5.根据权利要求1所述的一种基于位置点匹配的手腕骨兴趣区域修复方法,其特征在于:所述的步骤1具体包括:
41)确定REGIONS中与Mt匹配距离最短的骨骼区域target,公式为:
Figure FDA0002733215730000026
42)根据target在REGIONS中的索引INDEX,填补Bt中缺失的参照骨ROI,
Bt=Bt∪BONES[INDEX][u] (4-2)。
CN202011124718.3A 2020-10-20 2020-10-20 一种基于位置点匹配的手腕骨兴趣区域修复方法 Active CN112365411B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011124718.3A CN112365411B (zh) 2020-10-20 2020-10-20 一种基于位置点匹配的手腕骨兴趣区域修复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011124718.3A CN112365411B (zh) 2020-10-20 2020-10-20 一种基于位置点匹配的手腕骨兴趣区域修复方法

Publications (2)

Publication Number Publication Date
CN112365411A true CN112365411A (zh) 2021-02-12
CN112365411B CN112365411B (zh) 2024-02-20

Family

ID=74508050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011124718.3A Active CN112365411B (zh) 2020-10-20 2020-10-20 一种基于位置点匹配的手腕骨兴趣区域修复方法

Country Status (1)

Country Link
CN (1) CN112365411B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109961044A (zh) * 2019-03-22 2019-07-02 浙江工业大学 基于形状信息和卷积神经网络的chn法兴趣区域提取方法
CN111598839A (zh) * 2020-04-22 2020-08-28 浙江工业大学 一种基于孪生网络的手腕骨等级分类方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109961044A (zh) * 2019-03-22 2019-07-02 浙江工业大学 基于形状信息和卷积神经网络的chn法兴趣区域提取方法
CN111598839A (zh) * 2020-04-22 2020-08-28 浙江工业大学 一种基于孪生网络的手腕骨等级分类方法

Also Published As

Publication number Publication date
CN112365411B (zh) 2024-02-20

Similar Documents

Publication Publication Date Title
CN108492343B (zh) 一种扩充目标识别的训练数据的图像合成方法
JP7489732B2 (ja) 深層学習に基づく脊椎mri映像キーポイント検出方法
AU2008345628B2 (en) Remote sensing and probabilistic sampling based method for determining the carbon dioxide volume of a forest
US20170076442A1 (en) Generating Image-Based Diagnostic Tests By Optimizing Image Analysis and Data Mining Of Co-Registered Images
CN109325398A (zh) 一种基于迁移学习的人脸属性分析方法
CN109785337A (zh) 一种基于实例分割算法的栏内哺乳动物清点方法
CN108846828A (zh) 一种基于深度学习的病理图像目标区域定位方法及系统
CN103310457B (zh) 一种基于抛物线修正凸包的肺实质分割方法
CN110648322A (zh) 一种子宫颈异常细胞检测方法及系统
CN110097974A (zh) 一种基于深度学习算法的鼻咽癌远端转移预测系统
CN106971397B (zh) 基于改进jseg算法的城市高分辨率遥感影像分割方法
CN106408566A (zh) 一种胎儿超声图像质量控制方法及系统
CN106846484A (zh) 一种基于双视图三维重建的食物体积估计方法
CN104573685A (zh) 一种基于线性结构提取的自然场景文本检测方法
CN109708658B (zh) 一种基于卷积神经网络的视觉里程计方法
CN109284760A (zh) 一种基于深度卷积神经网络的家具检测方法及装置
CN108416304B (zh) 一种利用上下文信息的三分类人脸检测方法
CN117636160A (zh) 一种基于半监督学习的高分遥感耕地地块自动更新方法
CN106324708B (zh) 雨量记录图纸的数字化方法、装置
CN112686872B (zh) 基于深度学习的木材计数方法
CN107729863B (zh) 人体指静脉识别方法
CN112365411A (zh) 一种基于位置点匹配的手腕骨兴趣区域修复方法
CN108230322A (zh) 一种基于弱样本标记的眼底特征检测装置
CN113222939B (zh) 一种基于拇指甲标定的食物图像体积计算方法
CN113989269B (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