CN117218088B - 一种前臂x光影像的处理方法 - Google Patents
一种前臂x光影像的处理方法 Download PDFInfo
- Publication number
- CN117218088B CN117218088B CN202311197811.0A CN202311197811A CN117218088B CN 117218088 B CN117218088 B CN 117218088B CN 202311197811 A CN202311197811 A CN 202311197811A CN 117218088 B CN117218088 B CN 117218088B
- Authority
- CN
- China
- Prior art keywords
- radius
- forearm
- center
- ulna
- gypsum
- 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
- 210000000245 forearm Anatomy 0.000 title claims abstract description 75
- 238000003672 processing method Methods 0.000 title description 5
- 229910052602 gypsum Inorganic materials 0.000 claims abstract description 40
- 239000010440 gypsum Substances 0.000 claims abstract description 40
- 210000000623 ulna Anatomy 0.000 claims abstract description 38
- 238000003062 neural network model Methods 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000002372 labelling Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 abstract 1
- 208000010392 Bone Fractures Diseases 0.000 description 53
- 206010017076 Fracture Diseases 0.000 description 53
- 239000011505 plaster Substances 0.000 description 27
- 230000000875 corresponding effect Effects 0.000 description 26
- 230000000694 effects Effects 0.000 description 7
- 206010037802 Radius fracture Diseases 0.000 description 6
- 208000017013 ulna fracture Diseases 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本申请涉及图像数据处理或产生技术领域,特别是涉及一种前臂X光影像的处理方法。所述方法包括以下步骤:获取目标前臂X光影像P;获取P中桡骨的骨折区域A1的最小矩形包围盒B1;获取B1的第一位置(x’1,y’1);获取P中桡骨近端中心的位置和获取P中桡骨远端中心的位置;获取桡骨对应的第一比例δ1;获取尺骨对应的第二比例δ2;获取预设的石膏固定范围列表H;遍历H,将(δ1,δ2)和hr输入经训练的目标神经网络模型,获取(δ1,δ2)和hr的匹配度ζr,并将ζr追加至预设的匹配度集合G;将max(G)对应的石膏固定范围确定为与P匹配的石膏固定范围。本发明可为用户推荐一个石膏固定范围。
Description
技术领域
本发明涉及图像数据处理或产生技术领域,特别是涉及一种前臂X光影像的处理方法。
背景技术
石膏固定具有能够根据肢体的形状进行塑形且塑形之后的固定作用较为可靠的优点,现有技术中当前臂出现骨折时,一般是给出现骨折的整个前臂都打上石膏,但是,石膏往往比较沉重,对骨折患者来说也是一种负担。如何为用户推荐一个既不影响对前臂的固定效果又较小的石膏固定范围,是亟待解决的技术问题。
发明内容
本发明目的在于,提供前臂X光影像的处理方法,以为用户推荐一个既不影响对前臂的固定效果又较小的石膏固定范围。
根据本发明的第一方面,提供了一种前臂X光影像的处理方法,包括以下步骤:
S100,获取目标前臂X光影像P;P={p1,p2,…,pi,…,pn},pi为P中第i个像素点,i的取值范围为1到n,n为P中像素点的数量。
S200,获取P中桡骨的骨折区域A1的最小矩形包围盒B1,A1={p1 1,p1 2,…,p1 a1,…,p1 N1},p1 a1为P中桡骨中第a1个被判定为骨折像素点的像素点,B1={b1 1,b1 2,…,b1 j1,…,b1 m1},b1 j1为B1中第j1个子区域,b1 j1={b1 j1,1,b1 j1,2,…,b1 j1,k1,…,b1 j1,u1},b1 j1,k1为b1 j1中第k1个像素点,a1的取值范围为1到N1,N1为P中桡骨被判定为骨折像素点的像素点的数量,j1的取值范围为1到m1,m1为B1中子区域的数量,k1的取值范围为1到u1,u1为b1 j1中像素点的数量。
S300,获取B1的第一位置(x’1,y’1),x’1为B1的第一位置的横坐标,y’1为B1的第一位置的纵坐标,x’1=∑m1 j1=1(qj1×xj1/u1)/m1,y’1=∑m1 j1=1(qj1×yj1/u1)/m1,qj1为A1∩b1 j1中像素点的数量,xj1为b1 j1的中心位置的横坐标,yj1为b1 j1的中心位置的纵坐标。
S400,获取P中桡骨近端中心的位置(x- 1,y- 1)和获取P中桡骨远端中心的位置(x+ 1,y+ 1),x- 1为P中桡骨近端中心的位置的横坐标,y- 1为P中桡骨近端中心的位置的纵坐标,x+ 1为P中桡骨远端中心的位置的横坐标,y+ 1为P中桡骨远端中心的位置的纵坐标。
S500,获取桡骨对应的第一比例δ1,δ1=(d2 1,2-d2 1,1)0.5/d1,3,d1,1为(x’1,y’1)到(x- 1,y- 1)与(x+ 1,y+ 1)的连线的距离,d1,2为(x’1,y’1)与(x- 1,y- 1)的距离,d1,3为(x- 1,y- 1)与(x+ 1,y+ 1)的距离。
S600,获取尺骨对应的第二比例δ2。
S700,获取预设的石膏固定范围列表H,H=(h1,h2,…,hr,…,hR),hr为H中第r个石膏固定范围,r的取值范围为1到R,R为H中石膏固定范围的数量。
S800,遍历H,将(δ1,δ2)和hr输入经训练的目标神经网络模型,获取(δ1,δ2)和hr的匹配度ζr,并将ζr追加至预设的匹配度集合G,G的初始化为空集。
S900,将max(G)对应的石膏固定范围确定为与P匹配的石膏固定范围,max()为取最大值。
本发明与现有技术相比,其至少具有以下有益效果:
本发明获取了目标前臂X光影像,并获取了该影像中桡骨的骨折区域的最小矩形包围盒,基于该最小矩形包围盒的每一子区域的中心位置的坐标和该子区域中包括的骨折像素点的数量,本发明获取了最小矩形包围盒的第一位置,该第一位置可以用于表示桡骨骨折的严重区域的位置;基于该第一位置、P中桡骨近端中心的位置和P中桡骨远端中心的位置,本发明获取了桡骨对应的第一比例,该第一比例可以用于表征桡骨骨折相对位置;本发明的尺骨对应的第二比例可以表征尺骨骨折相对位置;该第一比例结合第二比例可以唯一表征前臂的骨折情况,本发明将第一比例、第二比例和预设的石膏固定范围输入经训练的目标神经网络模型,可以得到该第一比例和第二比例是否与该预设的石膏固定范围匹配;本发明通过将第一比例、第二比例和预设的石膏固定范围列表中的每一石膏固定范围输入到经训练的目标神经网络模型,本发明得到了预设的石膏固定范围列表中与该第一比例和第二比例最匹配的石膏固定范围,并将该最匹配的石膏固定范围作为与目标前臂X光影像匹配的石膏固定范围。本发明利用经训练的神经网络模型得到了与目标前臂X光影像匹配的石膏固定范围,该石膏固定范围是根据目标前臂X光影像对应的桡骨骨折位置和尺骨骨折位置得到的,其是一个在不影响对前臂的固定效果的前提下的相较于将整个前臂都打上石膏的范围而言相对较小的范围,将该石膏固定范围输出在用户界面上可以作为用户在考量石膏固定范围时的一个参考。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的前臂X光影像的处理方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
根据本发明,提供了一种前臂X光影像的处理方法,如图1所示,包括以下步骤:
S100,获取目标前臂X光影像P;P={p1,p2,…,pi,…,pn},pi为P中第i个像素点,i的取值范围为1到n,n为P中像素点的数量。
S200,获取P中桡骨的骨折区域A1的最小矩形包围盒B1,A1={p1 1,p1 2,…,p1 a1,…,p1 N1},p1 a1为P中桡骨中第a1个被判定为骨折像素点的像素点,B1={b1 1,b1 2,…,b1 j1,…,b1 m1},b1 j1为B1中第j1个子区域,b1 j1={b1 j1,1,b1 j1,2,…,b1 j1,k1,…,b1 j1,u1},b1 j1,k1为b1 j1中第k1个像素点,a1的取值范围为1到N1,N1为P中桡骨被判定为骨折像素点的像素点的数量,j1的取值范围为1到m1,m1为B1中子区域的数量,k1的取值范围为1到u1,u1为b1 j1中像素点的数量。
可选的,利用经训练的第一神经网络模型识别P中桡骨中的骨折像素点。本实施例中骨折像素点即对应位置为骨折位置的像素点,非骨折像素点即对应位置不为骨折位置的像素点。可选的,对第一神经网络模型的训练过程包括:
S210,获取前臂X光影像样本集合。
S220,对前臂X光影像样本集合中每一前臂X光影像样本中桡骨中的骨折像素点进行标注。
可选的,利用人工对每一前臂X光影像样本中桡骨中的骨折像素点进行标注。
S230,将前臂X光影像样本集合和标注结果输入到第一神经网络模型进行训练。
本实施例中经训练的第一神经网络模型具有识别前臂X光影像中桡骨中的骨折像素点的功能。
本领域技术人员知悉,现有技术中任何的获取最小矩形包围盒的方法均落入本发明的保护范围。
S300,获取B1的第一位置(x’1,y’1),x’1为B1的第一位置的横坐标,y’1为B1的第一位置的纵坐标,x’1=∑m1 j1=1(qj1×xj1/u1)/m1,y’1=∑m1 j1=1(qj1×yj1/u1)/m1,qj1为A1∩b1 j1中像素点的数量,xj1为b1 j1的中心位置的横坐标,yj1为b1 j1的中心位置的纵坐标。
本实施例中m1=round(v1/v0)×round(w1/w0),v1为B1的长度,w1为B1的宽度,v0为预设的子区域的初始长度,w0为预设的子区域的初始宽度,round()为四舍五入取整。本实施例中w0和v0为经验值,B1中每一子区域的长度为v1/round(v1/v0),B1中每一子区域的宽度为w1/round(w1/w0),如果出现某子区域中包括不完整的像素点的情况,则舍弃该子区域中不完整的像素点,将舍弃后得到的完整的像素点的集合作为该子区域的像素点的集合,并将该舍弃后得到的完整的像素点的集合与A1求交集后得到的交集中像素点的数量作为该子区域对应的qj1;本实施例中不完整的像素点指的是只有部分面积落入子区域中的像素点,完整的像素点指的是全部面积落入子区域中的像素点。
本实施例B1中每一子区域的面积相较于单一像素点的面积较大,m1小于B1中像素点的数量,由此,本实施例获取B1的第一位置的过程中执行求和运算的时间较短,提高了获取B1的第一位置的效率。
S400,获取P中桡骨近端中心的位置(x- 1,y- 1)和获取P中桡骨远端中心的位置(x+ 1,y+ 1),x- 1为P中桡骨近端中心的位置的横坐标,y- 1为P中桡骨近端中心的位置的纵坐标,x+ 1为P中桡骨远端中心的位置的横坐标,y+ 1为P中桡骨远端中心的位置的纵坐标。
可选的,利用经训练的第二神经网络模型同时识别P中桡骨中的骨折像素点、P中桡骨近端中心的位置和P中桡骨远端中心的位置。可选的,对第二神经网络模型的训练过程包括:
S410,获取前臂X光影像样本集合。
S420,对前臂X光影像样本集合中每一前臂X光影像样本中桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
可选的,利用人工对每一前臂X光影像样本中桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
S430,将前臂X光影像样本集合和标注结果输入到第二神经网络模型进行训练。
本实施例中经训练的第二神经网络模型同时具有识别前臂X光影像中桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置的功能。
本领域技术人员知悉,现有技术中任何的可用于多目标识别的语义分割模型均落入本发明的保护范围。
S500,获取桡骨对应的第一比例δ1,δ1=(d2 1,2-d2 1,1)0.5/d1,3,d1,1为(x’1,y’1)到(x- 1,y- 1)与(x+ 1,y+ 1)的连线的距离,d1,2为(x’1,y’1)与(x- 1,y- 1)的距离,d1,3为(x- 1,y- 1)与(x+ 1,y+ 1)的距离。
S600,获取尺骨对应的第二比例δ2。
本实施例中获取尺骨对应的第二比例δ2的过程与上述获取桡骨对应的第一比例δ1的过程类似,包括以下步骤:
S610,获取P中尺骨的骨折区域A2的最小矩形包围盒B2,A2={p0 1,p0 2,…,p0 a2,…,p0 N2},p0 a2为P中尺骨中第a2个被判定为骨折像素点的像素点,B2={b0 1,b0 2,…,b0 j2,…,b0 m2},b0 j2为B2中第j2个子区域,b0 j2={b0 j2,1,b0 j2,2,…,b0 j2,k2,…,b0 j2,u2},b0 j2,k2为b0 j2中第k2个像素点,a2的取值范围为1到N2,N2为P中尺骨被判定为骨折像素点的像素点的数量,j2的取值范围为1到m2,m2为B2中子区域的数量,k2的取值范围为1到u2,u2为b0 j2中像素点的数量。
可选的,利用经训练的第三神经网络模型同时识别P中尺骨中的骨折像素点、P中桡骨中的骨折像素点、P中桡骨近端中心的位置和P中桡骨远端中心的位置。可选的,对第三神经网络模型的训练过程包括:
S611,获取前臂X光影像样本集合。
S612,对前臂X光影像样本集合中每一前臂X光影像样本中尺骨中的骨折像素点、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
可选的,利用人工对每一前臂X光影像样本中尺骨中的骨折像素点、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
S613,将前臂X光影像样本集合和标注结果输入到第三神经网络模型进行训练。
本实施例中经训练的第三神经网络模型同时具有识别前臂X光影像中尺骨中的骨折像素点、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置的功能。
本领域技术人员知悉,现有技术中任何的可用于多目标识别的语义分割模型均落入本发明的保护范围。
本领域技术人员知悉,现有技术中任何的获取最小矩形包围盒的方法均落入本发明的保护范围。
S620,获取B2的第一位置(x’2,y’2),x’2为B2的第一位置的横坐标,y’2为B2的第一位置的纵坐标,x’2=∑m2 j2=1(qj2×xj2/u2)/m2,y’2=∑m2 j2=1(qj2×yj2/u2)/m2,qj2为A2∩b0 j2中像素点的数量,xj2为b0 j2的中心位置的横坐标,yj2为b0 j2的中心位置的纵坐标。
本实施例中m2=round(v2/v0)×round(w2/w0),v2为B2的长度,w2为B2的宽度。本实施例中B2中每一子区域的长度为v2/round(v2/v0),B2中每一子区域的宽度为w2/round(w2/w0),如果出现某子区域中包括不完整的像素点的情况,则舍弃该子区域中不完整的像素点,将舍弃后得到的完整的像素点的集合作为该子区域的像素点的集合,并将该舍弃后得到的完整的像素点的集合与A2求交集后得到的交集中像素点的数量作为该子区域对应的qj2。
本实施例B2中每一子区域的面积相较于单一像素点的面积较大,m2小于B2中像素点的数量,由此,本实施例获取B2的第一位置的过程中执行求和运算的时间较短,提高了获取B2的第一位置的效率。
S630,获取P中尺骨近端中心的位置(x- 2,y- 2)和获取P中尺骨远端中心的位置(x+ 2,y+ 2),x- 2为P中尺骨近端中心的位置的横坐标,y- 2为P中尺骨近端中心的位置的纵坐标,x+ 2为P中尺骨远端中心的位置的横坐标,y+ 2为P中尺骨远端中心的位置的纵坐标。
可选的,利用经训练的第四神经网络模型同时识别P中尺骨中的骨折像素点、P中尺骨近端中心的位置、P中尺骨远端中心的位置、P中桡骨中的骨折像素点、P中桡骨近端中心的位置和P中桡骨远端中心的位置。可选的,对第四神经网络模型的训练过程包括:
S631,获取前臂X光影像样本集合。
S632,对前臂X光影像样本集合中每一前臂X光影像样本中尺骨中的骨折像素点、尺骨近端中心的位置、尺骨远端中心的位置、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
可选的,利用人工对每一前臂X光影像样本中尺骨中的骨折像素点、尺骨近端中心的位置、尺骨远端中心的位置、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置进行标注。
S633,将前臂X光影像样本集合和标注结果输入到第四神经网络模型进行训练。
本实施例中经训练的第四神经网络模型同时具有识别前臂X光影像中尺骨中的骨折像素点、尺骨近端中心的位置、尺骨远端中心的位置、桡骨中的骨折像素点、桡骨近端中心的位置和桡骨远端中心的位置的功能。
本领域技术人员知悉,现有技术中任何的可用于多目标识别的语义分割模型均落入本发明的保护范围。
S640,获取尺骨对应的第二比例δ2。
具体的,δ2=(d2 2,2-d2 2,1)0.5/d2,3,d2,1为(x’2,y’2)到(x- 2,y- 2)与(x+ 2,y+ 2)的连线的距离,d2,2为(x’2,y’2)与(x- 2,y- 2)的距离,d2,3为(x- 2,y- 2)与(x+ 2,y+ 2)的距离,x- 2为P中尺骨近端中心的位置的横坐标,y- 2为P中尺骨近端中心的位置的纵坐标,x+ 2为P中尺骨远端中心的位置的横坐标,y+ 2为P中尺骨远端中心的位置的纵坐标,x’2为B2的第一位置的横坐标,y’2为B2的第一位置的纵坐标,B2为P中尺骨的骨折区域的最小矩形包围盒。
S700,获取预设的石膏固定范围列表H,H=(h1,h2,…,hr,…,hR),hr为H中第r个石膏固定范围,r的取值范围为1到R,R为H中石膏固定范围的数量。
本实施例中H中任意两石膏固定范围不相同,每一hr为[0,1]的子范围。例如,[0,1]中的0表示从前臂的近端,1表示前臂的远端,当石膏固定范围为[0,1]时,表示给前臂都打上石膏;[0,1/2]中的1/2前臂的近端与前臂的远端的中间位置;[0,1/2]中的1/2表示前臂的近端与前臂的远端的中间位置,当石膏固定范围为[0,1/2]时,表示给前臂的近端到该中间位置之间的部分打上石膏;[1/4,3/4]中的1/4表示前臂上距离前臂的近端1/4倍的前臂总长的位置,3/4表示前臂上距离前臂的近端3/4倍的前臂总长的位置,当石膏固定范围为[1/4,3/4]时,表示给前臂上距离前臂的近端1/4倍的前臂总长的位置到3/4倍的前臂总长的位置之间的部分打上石膏;以此类推。
S800,遍历H,将(δ1,δ2)和hr输入经训练的目标神经网络模型,获取(δ1,δ2)和hr的匹配度ζr,并将ζr追加至预设的匹配度集合G,G的初始化为空集。
具体的,目标神经网络模型的训练过程包括:
S810,获取第一样本数据S,S=(s1,s2,…,sc,…,se),sc为第c个第一样本数据,sc=(δc,1,δc,2),c的取值范围为1到e,e为样本的数量,δc,1为第c个样本的桡骨对应的第一比例,δc,2为第c个样本的尺骨对应的第二比例。
S820,获取第一样本数据F,F=(f1,f2,…,fc,…,fe),fc为第c个样本对应的石膏固定范围。
优选的,fc=[min(max(0,δc,1-d0/d1,3),max(0,δc,2-d0/d2,3)),min(1,max(δc,1+d0/d1,3,δc,2+d0/d2,3))],d0为预设的石膏固定最小长度阈值,min()为取最小值。本实施例中d0对应的是石膏相较于骨折位置的延伸长度,可用于保证石膏对骨裂位置的固定作用,可选的,d0为经验值;本实施例中在桡骨和尺骨均存在骨折的情况下,结合桡骨的骨折位置和尺骨的骨折位置优选能够同时满足桡骨骨折位置固定和尺骨骨折位置固定的石膏固定范围,保证了石膏对骨裂位置的固定效果。
S830,获取标签E,E=(g1,g2,…,gc,…,ge),gc为sc的标签,0≤gc≤1,gc与sc和fc之间的匹配度正相关。
可选的,gc为人工的经验值。当gc=0时,表示sc和fc之间的匹配度为0,即sc和fc之间完全不匹配;当gc=1时,表示sc和fc之间的匹配度为1,即sc和fc之间完全匹配。本实施例中gc越大,表示sc和fc之间越匹配。
S840,利用S、F和E对目标神经网络模型进行训练。
本实施例中经训练的目标神经网络模型可以用于判定输入的桡骨对应的第一比例和尺骨对应的第二比例与输入的石膏固定范围的匹配度。
可选的,本实施例的目标神经网络模型包括第一编码器、第二编码器和判别器,所述第一编码器的用于对桡骨对应的第一比例和尺骨对应的第二比例进行编码,得到第一特征向量;所述第二编码器用于对石膏固定范围进行编码,得到第二特征向量;所述判别器用于获取第一特征向量和第二特征向量的匹配度。
在对目标神经网络模型进行训练的过程中,判别器的输入为第一特征向量、第二特征向量以及对应的标签;经训练的目标神经网络模型中判别器具有为第一特征向量和第二特征向量打标签的功能,即具有判别第一特征向量和第二特征向量之间匹配度的功能。
本领域技术人员知悉,可用于实现上述功能的编码器和判别器的结构为现有技术,此处不再赘述。
S900,将max(G)对应的石膏固定范围确定为与P匹配的石膏固定范围,max()为取最大值。
本发明获取了目标前臂X光影像,并获取了该影像中桡骨的骨折区域的最小矩形包围盒,基于该最小矩形包围盒的每一子区域的中心位置的坐标和该子区域中包括的骨折像素点的数量,本发明获取了最小矩形包围盒的第一位置,该第一位置可以用于表示桡骨骨折的严重区域的位置;基于该第一位置、P中桡骨近端中心的位置和P中桡骨远端中心的位置,本发明获取了桡骨对应的第一比例,该第一比例可以用于表征桡骨骨折相对位置;本发明的尺骨对应的第二比例可以表征尺骨骨折相对位置;该第一比例结合第二比例可以唯一表征前臂的骨折情况,本发明将第一比例、第二比例和预设的石膏固定范围输入经训练的目标神经网络模型,可以得到该第一比例和第二比例是否与该预设的石膏固定范围匹配;本发明通过将第一比例、第二比例和预设的石膏固定范围列表中的每一石膏固定范围输入到经训练的目标神经网络模型,本发明得到了预设的石膏固定范围列表中与该第一比例和第二比例最匹配的石膏固定范围,并将该最匹配的石膏固定范围作为与目标前臂X光影像匹配的石膏固定范围。本发明利用经训练的神经网络模型得到了与目标前臂X光影像匹配的石膏固定范围,该石膏固定范围是根据目标前臂X光影像对应的桡骨骨折位置和尺骨骨折位置得到的,其是一个在不影响对前臂的固定效果的前提下的相较于将整个前臂都打上石膏的范围而言相对较小的范围,将该石膏固定范围输出在用户界面上可以作为用户在考量石膏固定范围时的一个参考。
虽然已经通过示例对本发明的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本发明的范围。本发明的范围由所附权利要求来限定。
Claims (6)
1.一种前臂X光影像的处理方法,其特征在于,包括以下步骤:
S100,获取目标前臂X光影像P;P={p1,p2,…,pi,…,pn},pi为P中第i个像素点,i的取值范围为1到n,n为P中像素点的数量;
S200,获取P中桡骨的骨折区域A1的最小矩形包围盒B1,A1={p1 1,p1 2,…,p1 a1,…,p1 N1},p1 a1为P中桡骨中第a1个被判定为骨折像素点的像素点,B1={b1 1,b1 2,…,b1 j1,…,b1 m1},b1 j1为B1中第j1个子区域,b1 j1={b1 j1,1,b1 j1,2,…,b1 j1,k1,…,b1 j1,u1},b1 j1,k1为b1 j1中第k1个像素点,a1的取值范围为1到N1,N1为P中桡骨被判定为骨折像素点的像素点的数量,j1的取值范围为1到m1,m1为B1中子区域的数量,k1的取值范围为1到u1,u1为b1 j1中像素点的数量;
S300,获取B1的第一位置(x’1,y’1),x’1为B1的第一位置的横坐标,y’1为B1的第一位置的纵坐标,x’1=∑m1 j1=1(qj1×xj1/u1)/m1,y’1=∑m1 j1=1(qj1×yj1/u1)/m1,qj1为A1∩b1 j1中像素点的数量,xj1为b1 j1的中心位置的横坐标,yj1为b1 j1的中心位置的纵坐标;
S400,获取P中桡骨近端中心的位置(x- 1,y- 1)和获取P中桡骨远端中心的位置(x+ 1,y+ 1),x- 1为P中桡骨近端中心的位置的横坐标,y- 1为P中桡骨近端中心的位置的纵坐标,x+ 1为P中桡骨远端中心的位置的横坐标,y+ 1为P中桡骨远端中心的位置的纵坐标;
S500,获取桡骨对应的第一比例δ1,δ1=(d2 1,2-d2 1,1)0.5/d1,3,d1,1为(x’1,y’1)到(x- 1,y- 1)与(x+ 1,y+ 1)的连线的距离,d1,2为(x’1,y’1)与(x- 1,y- 1)的距离,d1,3为(x- 1,y- 1)与(x+ 1,y+ 1)的距离;
S600,获取尺骨对应的第二比例δ2;δ2=(d2 2,2-d2 2,1)0.5/d2,3,d2,1为(x’2,y’2)到(x- 2,y- 2)与(x+ 2,y+ 2)的连线的距离,d2,2为(x’2,y’2)与(x- 2,y- 2)的距离,d2,3为(x- 2,y- 2)与(x+ 2,y+ 2)的距离,x- 2为P中尺骨近端中心的位置的横坐标,y- 2为P中尺骨近端中心的位置的纵坐标,x+ 2为P中尺骨远端中心的位置的横坐标,y+ 2为P中尺骨远端中心的位置的纵坐标,x’2为B2的第一位置的横坐标,y’2为B2的第一位置的纵坐标,B2为P中尺骨的骨折区域的最小矩形包围盒;
S700,获取预设的石膏固定范围列表H,H=(h1,h2,…,hr,…,hR),hr为H中第r个石膏固定范围,r的取值范围为1到R,R为H中石膏固定范围的数量;
S800,遍历H,将(δ1,δ2)和hr输入经训练的目标神经网络模型,获取(δ1,δ2)和hr的匹配度ζr,并将ζr追加至预设的匹配度集合G,G的初始化为空集;
S900,将max(G)对应的石膏固定范围确定为与P匹配的石膏固定范围,max()为取最大值。
2.根据权利要求1所述的前臂X光影像的处理方法,其特征在于,目标神经网络模型的训练过程包括:
S810,获取第一样本数据S,S=(s1,s2,…,sc,…,se),sc为第c个第一样本数据,sc=(δc,1,δc,2),c的取值范围为1到e,e为样本的数量,δc,1为第c个样本的桡骨对应的第一比例,δc,2为第c个样本的尺骨对应的第二比例;
S820,获取第一样本数据F,F=(f1,f2,…,fc,…,fe),fc为第c个样本对应的石膏固定范围;
S830,获取标签E,E=(g1,g2,…,gc,…,ge),gc为sc的标签,0≤gc≤1,gc与sc和fc之间的匹配度正相关;
S840,利用S、F和E对目标神经网络模型进行训练。
3.根据权利要求2所述的前臂X光影像的处理方法,其特征在于,fc=[min(max(0,δc,1-d0/d1,3),max(0,δc,2-d0/d2,3)),min(1,max(δc,1+d0/d1,3,δc,2+d0/d2,3))],d0为预设的石膏固定最小长度阈值,min()为取最小值。
4.根据权利要求1所述的前臂X光影像的处理方法,其特征在于,所述目标神经网络模型包括第一编码器、第二编码器和判别器,所述第一编码器用于对桡骨对应的第一比例和尺骨对应的第二比例进行编码,得到第一特征向量;所述第二编码器用于对石膏固定范围进行编码,得到第二特征向量;所述判别器用于获取第一特征向量和第二特征向量的匹配度。
5.根据权利要求1所述的前臂X光影像的处理方法,其特征在于,利用经训练的第一神经网络模型获取p1 a1,所述第一神经网络模型的训练过程包括:
S210,获取前臂X光影像样本集合;
S220,对前臂X光影像样本集合中每一前臂X光影像样本中桡骨中的骨折像素点进行标注;
S230,将前臂X光影像样本集合和标注结果输入到第一神经网络模型进行训练。
6.根据权利要求1所述的前臂X光影像的处理方法,其特征在于,m1=round(v1/v0)×round(w1/w0),v1为B1的长度,w1为B1的宽度,v0为预设的子区域的初始长度,w0为预设的子区域的初始宽度,round()为四舍五入取整。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311197811.0A CN117218088B (zh) | 2023-09-15 | 2023-09-15 | 一种前臂x光影像的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311197811.0A CN117218088B (zh) | 2023-09-15 | 2023-09-15 | 一种前臂x光影像的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117218088A CN117218088A (zh) | 2023-12-12 |
CN117218088B true CN117218088B (zh) | 2024-03-29 |
Family
ID=89049117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311197811.0A Active CN117218088B (zh) | 2023-09-15 | 2023-09-15 | 一种前臂x光影像的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117218088B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009100943A (ja) * | 2007-10-23 | 2009-05-14 | Aloka Co Ltd | 骨測定装置および骨画像処理方法 |
JP2015065999A (ja) * | 2013-09-26 | 2015-04-13 | 日立アロカメディカル株式会社 | 骨密度測定装置 |
RU2645616C1 (ru) * | 2017-04-28 | 2018-02-26 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Нижегородская государственная медицинская академия" Министерства здравоохранения Российской Федерации (ФГБОУ ВО НижГМА Минздрава России) | Способ фиксации дистального сухожилия двуглавой мышцы плеча к бугристости лучевой кости |
CN113920088A (zh) * | 2021-10-11 | 2022-01-11 | 浙江康源医疗器械有限公司 | 一种基于深度学习的桡骨密度检测方法、系统、装置 |
CN114742760A (zh) * | 2022-03-09 | 2022-07-12 | 宁波工程学院 | 一种x光片桡骨远端骨折自动快速诊断方法 |
CN115170562A (zh) * | 2022-09-06 | 2022-10-11 | 天河超级计算淮海分中心 | 一种x光影像处理系统 |
CN116468652A (zh) * | 2023-02-06 | 2023-07-21 | 北京理工大学 | 一种前臂骨标志点智能识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110317871A1 (en) * | 2010-06-29 | 2011-12-29 | Microsoft Corporation | Skeletal joint recognition and tracking system |
-
2023
- 2023-09-15 CN CN202311197811.0A patent/CN117218088B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009100943A (ja) * | 2007-10-23 | 2009-05-14 | Aloka Co Ltd | 骨測定装置および骨画像処理方法 |
JP2015065999A (ja) * | 2013-09-26 | 2015-04-13 | 日立アロカメディカル株式会社 | 骨密度測定装置 |
RU2645616C1 (ru) * | 2017-04-28 | 2018-02-26 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Нижегородская государственная медицинская академия" Министерства здравоохранения Российской Федерации (ФГБОУ ВО НижГМА Минздрава России) | Способ фиксации дистального сухожилия двуглавой мышцы плеча к бугристости лучевой кости |
CN113920088A (zh) * | 2021-10-11 | 2022-01-11 | 浙江康源医疗器械有限公司 | 一种基于深度学习的桡骨密度检测方法、系统、装置 |
CN114742760A (zh) * | 2022-03-09 | 2022-07-12 | 宁波工程学院 | 一种x光片桡骨远端骨折自动快速诊断方法 |
CN115170562A (zh) * | 2022-09-06 | 2022-10-11 | 天河超级计算淮海分中心 | 一种x光影像处理系统 |
CN116468652A (zh) * | 2023-02-06 | 2023-07-21 | 北京理工大学 | 一种前臂骨标志点智能识别方法 |
Non-Patent Citations (2)
Title |
---|
Tomographic PIV analysis of physiological flow conditions in a patient‑specific arteriovenous fistula;Sanjiv Gunasekera;《Experiments in Fluids》;20201231;第1-14页 * |
基于尺桡骨有限元模型的蒙医整骨手法分析;李忠贤;《中国组织工程研究》;20190731;第23卷(第20期);第3153-3157页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117218088A (zh) | 2023-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108056786B (zh) | 一种基于深度学习的骨龄检测方法和装置 | |
JP5031641B2 (ja) | パターンの位置合わせ方法、照合方法及び照合装置 | |
US11017210B2 (en) | Image processing apparatus and method | |
CN108197592B (zh) | 信息获取方法和装置 | |
US20080187183A1 (en) | Identification of Image Characteristics | |
CN110263755B (zh) | 眼底图像识别模型训练方法、眼底图像识别方法和设备 | |
JP2021089730A (ja) | 畳み込みニューラルネットワークによって、入力画像内で可視である少なくとも1つの対象要素を検出する方法 | |
JP6221770B2 (ja) | 画像処理装置、およびプログラム | |
CN112614586B (zh) | 一种基于医学影像和区块链的远程疾病智能诊断系统 | |
WO2018198500A1 (ja) | 照合装置、照合方法および照合プログラム | |
CN115337044B (zh) | 核酸采样监测方法、装置、系统与计算机可读存储介质 | |
CN117218088B (zh) | 一种前臂x光影像的处理方法 | |
Arambula Cosio et al. | Automated prostate recognition: a key process for clinically effective robotic prostatectomy | |
WO2022088729A1 (zh) | 点定位方法及相关装置、设备、介质及计算机程序 | |
KR101471814B1 (ko) | 홍채 인식에 의한 식별 | |
Sadik et al. | Variability in reference levels for Deauville classifications applied to lymphoma patients examined with 18F-FDG-PET/CT | |
CN117689891A (zh) | 一个两阶段模型用于分割和分类目标腰椎骨结构的方法 | |
WO2008047934A1 (fr) | Dispositif d'enregistrement, dispositif de collation, procédé d'authentification et programme d'authentification | |
US8515971B2 (en) | Method for assisting in making a decision on biometric data | |
CN110801204B (zh) | 一种基于人体框架的平衡检测方法 | |
Morris et al. | Segmentation of the finger bones as a prerequisite for the determination of bone age | |
JP2005305048A (ja) | 患者位置決め装置および患者位置決め方法 | |
JP2005339389A (ja) | 画像処理方法および画像処理装置 | |
US20240320911A1 (en) | Method for forensic identification by automatically comparing a 3d model of a skull with one or more photos of a face | |
Margeta et al. | Nautilus: A Clinical Tool for the Segmentation of Intra-Cochlear Structures and Related Applications |
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 |