CN104414685A - 超声波诊断装置、超声波图像处理方法 - Google Patents

超声波诊断装置、超声波图像处理方法 Download PDF

Info

Publication number
CN104414685A
CN104414685A CN201410415059.7A CN201410415059A CN104414685A CN 104414685 A CN104414685 A CN 104414685A CN 201410415059 A CN201410415059 A CN 201410415059A CN 104414685 A CN104414685 A CN 104414685A
Authority
CN
China
Prior art keywords
frame
disease score
ultrasound
disease
image signals
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.)
Pending
Application number
CN201410415059.7A
Other languages
English (en)
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.)
Konica Minolta Inc
Konica Minolta Opto Inc
Original Assignee
Konica Minolta Opto Inc
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 Konica Minolta Opto Inc filed Critical Konica Minolta Opto Inc
Publication of CN104414685A publication Critical patent/CN104414685A/zh
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0875Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4209Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/462Displaying means of special interest characterised by constructional features of the display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4405Device being mounted on a trolley
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/486Diagnostic techniques involving arbitrary m-mode
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • G06T2207/101363D ultrasound image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Rheumatology (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Quality & Reliability (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

提供一种超声波诊断装置,在使用了超声波诊断装置的风湿诊断中,能够降低评价结果的检测者依赖性。超声波诊断装置(1100)具备:超声波图像取得部(2001),取得多个帧的超声波图像信号;评价对象判定部(2002),分析多个帧的超声波图像信号,判定在各帧中是否包含表示关节部位的对象图像部分;疾病分数计算部(2003),根据评价对象帧的超声波图像信号中的对象图像部分的超声波图像信号,计算将疾病的程度量化后的疾病分数;选择部(2004),对所计算出的疾病分数,基于规定的数值处理而选择1个以上的疾病分数;以及显示控制部(1006),使显示器(1008)显示所选择的疾病分数以及与该疾病分数对应的帧的超声波图像。

Description

超声波诊断装置、超声波图像处理方法
技术领域
本公开属于超声波图像处理方法的技术领域,尤其涉及对风湿疾病的程度进行定量化的超声波图像处理方法以及使用了该方法的超声波诊断装置。
背景技术
超声波诊断装置很早起就开始用于心脏、腹部、胎儿等的诊断、观察。近年来,超声波诊断装置的分辨率提升,也开始能够进行靠近体表面的部位的观察。与此相伴,即使在整形领域也开始使用超声波诊断装置,使用超声波诊断装置,进行骨、肌腱、肌肉等的观察、或风湿诊断中的关节部分的炎症的评价等。
风湿的炎症的评价使用B模式图像和多普勒模式图像而进行(例如,参照非专利文献1)。具体而言,使用B模式图像来观察关节腔的厚度、有无骨侵蚀和关节积液,并使用多普勒模式图像来观察有无新生血管。
在临床研究中,对炎症的程度进行定量化的方法受到关注,例如,提出了基于主观分类为4个阶段的半定量评价法、对矩形内的多普勒信号进行计数的BOX法、对将骨作为2个边的三角形内的多普勒信号进行计数的跟踪(trace)法等。这样的疾病的程度的定量化对于判定治疗效果是有用的。
现有技术文献
非专利文献
[非专利文献1]:小池隆夫,“超音波検査法を用いた関節リウマチの新しい診療”,p40-43,メディカルレビュー社,2010年3月10日
可是,根据从被检测体中的哪个断面取得超声波图像而进行评价,得到不同的评价结果。因此,存在评价结果依赖于检查者的技巧或主观那样的课题。
发明内容
本公开用于解决上述课题,其目的在于提供一种降低评价结果的检测者依赖性的超声波诊断装置、超声波图像处理方法以及计算机可读取的非暂时性记录介质。
为了达成上述目的,本公开的一个方式所涉及的超声波诊断装置是从经由超声波探头而从被检测体取得的多个帧的超声波图像信号中选择1个以上的超声波图像信号的帧,并在显示器上进行显示的超声波诊断装置,其特征在于,具备控制电路,所述控制电路具备:超声波图像取得部,取得所述多个帧的超声波图像信号;评价对象判定部,分析所述多个帧的超声波图像信号,在各帧中包含表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧;疾病分数计算部,根据所述评价对象帧的超声波图像信号中的所述对象图像部分的超声波图像信号,按每个所述评价对象帧计算将疾病的程度定量化后的疾病分数;选择部,对所述计算的疾病分数,基于规定的数值处理而选择1个以上的疾病分数;以及显示控制部,使显示器显示所述选择的疾病分数以及与该疾病分数对应的帧的超声波图像。
根据本公开,通过上述结构,在从多方向得到的超声波图像之中,能够通过基于规定的数值处理而选择出的断面的超声波图像,评价疾病的程度。因此,能够对在适当的技巧下得到的超声波图像进行选择从而客观地评价疾病的程度,并且能够降低疾病评价中的评价结果的检测者依赖性。
附图说明
图1是包含超声波诊断装置1100的超声波诊断系统1000的外观图。
图2是表示超声波诊断装置1100的内部结构的方框图。
图3是超声波诊断装置1100中的疾病定量化部1007的方框图。
图4是超声波诊断装置1100中的评价对象判定部2002的方框图。
图5是表示正在进行手指关节的风湿检查的状态的示意图。
图6是超声波诊断装置1100中的疾病分数计算处理的流程图。
图7是表示将手指关节作为测定对象而拍摄的超声波图像的示意图。
图8是表示超声波诊断装置1100中的关节检测处理的动作的流程图。
图9(a)、图9(b)是表示关节搜索处理所使用的关节检测窗的示意图。
图10是实施方式所涉及的超声波诊断装置1100所涉及的关节搜索处理中的动作说明图。
图11是表示超声波诊断装置1100中的基于模板匹配(テンプレートマッチング)法的关节检测处理的细节的流程图。
图12是表示使用机器学习的关节搜索器的一例的示意图。
图13是超声波诊断装置1100中的运动噪声(モーションノイズ)检测的流程图。
图14(a)、图14(b)是表示帧的B模式图像信号中产生了运动噪声的状态的一例的示意图。
图15是说明体表压迫度检测的处理的图。
图16是超声波诊断装置1100中的体表压迫检测处理的流程图。
图17是超声波诊断装置1100中的疾病分数计算处理的流程图。
图18是说明了检测用于表示骨的图像部分1、2的处理的示意图。
图19是说明了检测用于表示关节囊的边界的图像部分5的处理的示意图。
图20是说明了检测骨侵蚀7的程度的处理的示意图。
图21是示出了正常的骨表面的边界上的拟合(フィッティング)的一例的示意图。
图22是示出了骨侵蚀的骨表面的边界上的拟合的一例的示意图。
图23是用于注册检查顺序的操作输入画面。
图24是测量开始前的显示画面。
图25是测量中的显示画面。
图26是表示在图像中不包含关节部位的警告消息的显示画面。
图27是表示在多普勒模式图像中是由运动噪声引起的图像的警告消息的显示画面。
图28是表示超声波探头正在压迫被检测体体表的警告消息的显示画面。
图29是表示测量后的疾病分数的显示画面。
图30是表示测量完成的显示画面。
图31是校正操作中的显示画面。
图32是校正后的显示画面。
图33是示出了各手指关节中的检查结果的显示画面。
图34是表示检查结果的显示画面,该检查结果表示根据手指关节中的多个帧的超声波图像计算出的疾病分数。
标号说明
1000 超声波诊断系统
1001 超声波探头
1002 超声波发送接收部
1003 B模式处理部
1004 多普勒模式处理部
1005 存储器(存储装置)
1006 显示控制部
1007 疾病定量化部
1008 显示器
1009 操作单元
1010 控制电路
1100 超声波诊断装置
2001 超声波图像取得部
2002 评价对象判定部
2003 疾病分数计算部
2003A 形态定量化部
2003B 炎症定量化部
2004 选择部
3001 关节检测部
3002 压迫检测部
3003 运动噪声检测部
3004 技巧判定部
S1000 读出超声波图像信号的步骤
S1001、S1002、S1003 判定是否是评价对象帧的判定步骤
S1004 计算疾病分数的步骤
S1005 基于规定的数值处理选择1个以上的疾病分数的选择步骤
S1006 显示所选择出的疾病分数、与疾病分数对应的帧的超声波图像的步骤
具体实施方式
《实施方式》
<关于结构>
以下,关于实施方式所涉及的超声波诊断装置1100、超声波图像处理方法、包含超声波诊断装置1100的超声波诊断系统1000,参照附图进行说明。
1.超声波诊断系统1000的整体结构
说明实施方式所涉及的超声波诊断装置1100的概略结构。图1是包含超声波诊断装置1100的超声波诊断系统1000的外观图。超声波诊断系统1000由超声波诊断装置1100、超声波探头1001、显示器1008、操作单元1009构成。以下,说明构成超声波诊断系统1000的各装置。
(1)超声波探头1001
超声波探头1001将由超声波发送接收部1002生成的发送波向被检测体体内进行发送,并接收由被检测体内反射的反射超声波。由超声波探头1001接收到的反射超声波作为回声信号而输出至超声波发送接收部1002。
超声波探头1001例如具有由在一维方向(下面,设为“振子排列方向”)上排列的多个压电元件(不图示)构成的振动元件列,将从后述的超声波发送接收部1002供给的脉冲状的电信号(下面,设为“发送超声波信号”)变换为脉冲状的超声波。此时,超声波探头1001在使振动元件侧外表面与被检测体的皮肤表面接触了的状态下,将由从多个振动元件发出的超声波构成的超声波束向测定对象进行发送。然后,超声波探头1001接收来自被检测体的多个反射超声波,通过多个振动元件将这些反射超声波分别变换为电信号(下面,设为“接收超声波信号”),将接收超声波信号输出至超声波发送接收部。
经由这样的超声波探头1001向与振动元件列相对的被检测体的规定区域发送超声波束,进行用于接收从被检测体被反射的反射超声波的超声波扫描,从而能够取得断面中的超声波图像,该断面由振动元件列方向和垂直于该振动元件列方向的深度方向构成。将通过1次超声波扫描而取得的超声波图像的单位设为帧,表示为用于构筑1张断层图像所需的一组信号的单位。通过重复该超声波扫描,能够取得多个帧的超声波图像。
(2)显示器1008
显示器1008是所谓的图像显示用的显示装置,在显示画面上显示来自后述的显示控制部的图像输出。在显示器1008中,能够使用液晶显示器、CRT、有机EL显示器等。
(3)操作单元1009
操作单元1009是输入设备,接受操作者对于超声波诊断装置1100的各种设定/操作等的各种操作输入,输出至控制电路1010。由操作者所输入的信息是患者名、检查日期和时间、画面的再现/停止、保存、画质调整等。这些输入信息被存储在存储器1005中。
操作单元1009例如具体而言相当于键盘、轨迹球、触摸面板。在是触摸面板的情况下,也可以与显示器1008作为一体而构成。该情况下,对显示器1008所显示的操作键进行触摸操作或拖拽操作,从而能够进行超声波诊断装置1100的各种设定/操作,超声波诊断装置1100构成为能够通过该触摸面板来进行操作。
此外,操作单元1009也可以是具有各种操作用的键的键盘、或具有各种操作用的按钮、杆等的操作面板。此外,操作单元1009也可以是用于移动显示器1008所显示的光标显示的轨迹球、鼠标等。或者,也可以使用多个这些装置,也可以是组合了多个这些装置的结构的单元。
(4)超声波诊断装置1100
超声波诊断装置1100经由超声波探头1001发送接收用于超声波诊断的超声波信号,并对接收到的反射超声波信号形成图像。说明超声波诊断装置1100的内部结构。图2是表示超声波诊断装置1100的内部结构的功能框图。示出了对超声波诊断装置1100连接了超声波探头1001、显示器1008、操作单元1009的状态。
如图2所示,超声波诊断装置1100由控制电路1010构成,该控制电路1010具备了超声波发送接收部1002、B模式处理部1003、多普勒模式处理部1004、存储器1005、显示控制部1006、以及疾病定量化部1007。超声波发送接收部1002、B模式处理部1003、多普勒模式处理部1004、显示控制部1006、疾病定量化部1007分别由例如FPGA(Field Programmable Gate Array,现场可编程门阵列)、ASIC(Application Specific Integrated Circuit,专用集成电路)等的硬件电路来实现。或者,也可以是由CPU(Central Processing Unit,中央处理单元)或GPU(Graphics Processing Unit,图形处理单元)或处理器等的可编程设备和软件来实现的结构。这些结构要素能够设为一个电路部件,也能够设为多个电路部件的集合体。此外,能够对多个结构要素进行组合而设为一个电路部件,也能够设为多个电路部件的集合体。
构成为超声波探头1001能够连接到超声波发送接收部1002,显示器1008能够连接到显示控制部1006。此外,构成为接受来自操作者的输入的操作单元1009能够连接到控制电路1010。
上面是关于构成超声波诊断系统1000的各装置的说明。
2.超声波诊断装置1100的各部分结构
接着,说明超声波诊断装置1100中包含的各要素的结构。
(1)超声波发送接收部1002
超声波发送接收部1002与超声波探头1001连接。超声波发送接收部1002进行对超声波探头1001供给脉冲状的发送超声波信号的发送处理,用于使超声波探头1001发送超声波束。具体而言,超声波发送接收部1002具备例如时钟产生电路、脉冲产生电路、延迟电路。时钟产生电路是产生用于决定超声波束的发送定时的时钟信号的电路。脉冲产生电路是用于产生驱动各振动元件的脉冲信号的电路。延迟电路是用于将超声波束的发送定时对各振动元件的每个设定延迟时间,使超声波束的发送延迟延迟时间的量,并进行超声波束的聚焦或转向(ステアリング)的电路。
此外,超声波发送接收部1002进行基于由超声波探头1001输入的回声信号按超声波扫描的顺序而进行波束成形,对B模式处理部1003以及多普勒模式处理部1004输出声音线信号(音響線信号)的接收处理。
另外,在波束成形中,对将从超声波探头1001取得的接收超声波信号放大后进行AD变换后的RF信号进行整相加法运算,生成在深度方向上连续的声音线信号。在此,RF信号例如由振动元件排列方向和超声波的发送方向即与振动元件排列垂直的方向构成的多个信号而构成,各信号是对从反射超声波的振幅被变换的电信号进行了A/D变换后的数字信号。声音线信号是在构成整相加法运算处理后的RF信号的深度方向上连续的数据。
超声波发送接收部1002重复地连续进行发送处理以及接收处理。
(2)B模式处理部1003
在B模式处理部1003中,基于由超声波发送接收部1002输入的声音线信号,生成B模式图像信号。所生成的B模式图像信号在存储器1005中暂时地存储。
具体而言,B模式处理部1003对声音线信号实施包络线检波、对数压缩等的处理而变换为与其强度对应的亮度信号,并将该亮度信号对直角坐标系进行坐标变换,从而生成B模式图像信号。B模式处理部1003按每个帧依次进行该处理,所生成的帧B模式图像信号在每次进行超声波扫描时被输出至存储器1005。此外,如后述那样,B模式图像信号被供给至显示控制部1006,并在显示器1008上被显示。
(3)多普勒模式处理部1004
在多普勒模式处理部1004中,对从超声波发送接收部1002输入的声音线信号实施自相关运算,从而提取多普勒分量,该多普勒分量是显示活体内的血流的时间变化的多普勒信号的源数据。然后,生成通过色彩来表现平均速度或方差、动力(パワー)等的血流信息的强度的多普勒模式图像信号。
具体而言,多普勒模式处理部1004进行对于声音线信号的自相关运算,在变换为流速后,通过滤波而提取表示血流分量的速度分量。进而,对滤波后的血流分量的平均速度、方差以及动力进行计算,生成多普勒模式图像信号。所生成的多普勒模式图像信号在每次进行超声波扫描时被输出至存储器1005。此外,如后述那样,多普勒模式图像信号被供给至显示控制部1006,并在显示器1008上被显示。
(4)存储器1005
存储器1005是将所生成的帧B模式图像信号以及多普勒模式图像信号作为输入,在每次进行超声波扫描时进行暂时性存储的存储装置。此外,存储由后述的疾病定量化部1007计算的疾病分数。存储器1005由包含使用了半导体存储器的DRAM(Dynamic Random Access Memory,动态随机存取存储器)、SRAM(Static Random Access Memory,静态随机存取存储器)等的RAM而构成。此外,也可以是通过硬盘驱动器、光盘驱动器、磁存储装置等而实现的结构。另外,在本说明书中,设为B模式图像信号以及多普勒模式图像信号构成超声波图像信号。
(5)疾病定量化部1007
疾病定量化部1007分析存储器1005所存储的B模式图像信号和多普勒模式图像信号,对风湿疾病的程度进行定量化。关于定量化的方法,在后面叙述。进行了定量化的结果按超声波图像信号的每个帧而被存储到存储器1005。
(6)显示控制部1006
显示控制部1006基于存储器1005所存储的B模式图像信号或多普勒模式图像信号而生成显示画面,将检查者名、患者名、时间信息、超声波诊断装置的设定信息、由疾病定量化部1007计算出的分数等重叠在显示画面上,在外部连接的显示器1008上进行显示。另外,超声波诊断装置1100根据需要,也可以设为包含超声波探头1001、操作单元1009、显示器1008的一部分或者全部的结构。
3.疾病定量化部1007的详细结构
接着,参照附图,说明疾病定量化部1007的详细结构。图3是疾病定量化部1007的方框图。
疾病定量化部1007由超声波图像取得部2001、评价对象判定部2002、疾病分数计算部2003、选择部2004构成。疾病分数计算部2003由形态定量化部2003A以及炎症定量化部2003B构成。疾病定量化部1007将存储器1005所存储的B模式图像信号4001以及多普勒模式图像信号4002作为输入,将表示疾病的程度的疾病分数4003、疾病分数4003的例如为最大值的最大疾病分数4004输出到存储器1005。
(1)超声波图像取得部2001
超声波图像取得部2001从存储器1005所存储的多个帧的B模式图像信号4001以及多普勒模式图像信号4002按每个帧而读出B模式图像信号4001以及多普勒模式图像信号4002。
(2)评价对象判定部2002
评价对象判定部2002将从超声波图像取得部2001输出的B模式图像信号4001和多普勒模式图像信号4002作为输入而对两个图像进行分析,判定是否是在操作者的适当的技巧下取得的超声波图像,将判定结果输出至形态定量化部2003A以及炎症定量化部2003B。
接着,参照附图说明评价对象判定部2002的详细结构。图4是评价对象判定部2002的方框图。如图4所示,评价对象判定部2002由关节检测部3001、压迫检测部3002、运动噪声检测部3003、技巧判定部3004构成。
i)关节检测部3001
关节检测部3001将从超声波图像取得部2001输出的B模式图像信号4001作为输入,按超声波图像信号的每个帧判定在超声波图像信号中是否包含表示关节部位的对象图像部分。判定结果是〔有关节〕〔没有关节〕的其中一个,且输出至技巧判定部3004。关于判定方法在后面叙述。
ii)压迫检测部3002
压迫检测部3002将从超声波图像取得部2001输出的B模式图像信号4001作为输入,根据图像来判定超声波图像信号的帧是否是在没有超声波探头对被检测体的体表的压迫的状态下被取得的。判定结果是〔有压迫〕〔没有压迫〕的其中一个,且输出至技巧判定部3004。关于判定方法在后面叙述。
iii)运动噪声检测部3003
运动噪声检测部3003将从超声波图像取得部2001输出的B模式图像信号4001以及多普勒模式图像信号4002作为输入,判定超声波图像信号的帧中的多普勒模式图像信号4002是否是由运动噪声引起的图像信号。判定结果是〔有运动噪声〕〔没有运动噪声〕的其中一个,且输出至技巧判定部3004。关于判定方法在后面叙述。
iv)技巧判定部3004
技巧判定部3004将关节检测部3001、压迫检测部3002、运动噪声检测部3003各自的判定结果作为输入。然后,在判定结果为〔有关节囊〕时,判定为在适当的技巧下得到的图像。或者,也可以在判定结果为〔有关节囊〕,且进而满足〔没有压迫〕或者〔没有运动噪声〕的至少一方或者双方时,判定为在适当的技巧下得到的图像。此时,技巧判定部3004对形态定量化部2003A和炎症定量化部2003B输出〔执行定量化〕。在除此之外的情况下输出〔中止定量化〕。
(3)形态定量化部2003A
返回图3,说明疾病分数计算部2003的结构。如上述那样的疾病分数计算部2003由形态定量化部2003A以及炎症定量化部2003B构成。
首先,说明形态定量化部2003A的结构。形态定量化部2003A将从超声波图像取得部2001输出的B模式图像信号4001、和评价对象判定部2002的判定结果作为输入。在判定结果为〔中止定量化〕的情况下,形态定量化部2003A不实施定量化而将疾病分数设为无效。另一方面,在判定结果为〔执行定量化〕的情况下,基于B模式图像信号4001中表示关节腔的图像部分的大小或表示骨的图像部分的亮度图案而计算疾病分数。由形态定量化部2003A计算出的疾病分数称为肿胀分数(GS)。肿胀分数被存储到存储器1005。关于肿胀分数的计算方法在后面叙述。
(4)炎症定量化部2003B
接着,说明炎症定量化部2003B的结构。炎症定量化部2003B将从超声波图像取得部2001输出的多普勒模式图像信号4002、和评价对象判定部2002的判定结果、以及来自形态定量化部2003A的表示关节腔的信息作为输入。在判定结果为〔中止定量化〕的情况下,炎症定量化部2003B不执行定量化而将疾病分数设为无效。另一方面,在判定结果为〔执行定量化〕的情况下,根据多普勒模式图像信号4002中的表示关节囊的图像部分5和表示位于与骨表面之间的关节腔的图像部分6中的、被检测出多普勒信号的图像区域的大小,计算疾病分数。由炎症定量化部2003B计算出的疾病分数称为炎症分数(PD)。炎症分数被存储到存储器1005。关于炎症分数的计算方法在后面叙述。
(5)选择部2004
选择部2004将由形态定量化部2003A和炎症定量化部2003B计算出且存储在存储器1005中的疾病分数(肿胀分数、炎症分数)作为输入。然后,基于规定的数值处理而选择最佳的1个以上的疾病分数。例如,也可以将多个帧中的疾病分数的最大值选择作为最佳的疾病分数。此外,也可以设为选择多个帧中的疾病分数的中值、平均值的结构。所选择的1以上的最佳的疾病分数被输出到存储器1005并存储。
此外,也可以设为以下结构:在将多个帧中的疾病分数的最大值选择为最佳的疾病分数的情况下,选择部2004输入在形态定量化部2003A和炎症定量化部2003B中计算的炎症分数,对最大值进行保持,在全部帧中的疾病分数被输入后将所保持的最大值作为最佳的疾病分数而输出至存储器1005。
<关于动作>
1.超声波诊断装置1100的动作
接着,使用手指关节的例子,说明进行风湿检查的情况下的超声波诊断装置1100的动作。在手指关节的风湿检查中,对从全部的手指关节中选择的关节进行检查。在对于1个关节的检查中,在关节内的不同的多个断面中进行超声波检查。
图5是表示使用超声波诊断装置1100正在进行手指关节的风湿检查的状态的示意图。如图5所示,以振动元件列与指的纵向方向平行的状态,将超声波探头1001沿着手指关节进行配置。在该状态下,将指作为旋转轴而使超声波探头1001如rtl、rt2那样旋转,从而以使对于关节部位的倾斜角度不同的状态而进行多次超声波扫描,对1个手指关节拍摄多张超声波图像。然后,基于所取得的多张超声波图像进行对于1个手指关节的风湿疾病的检查。
(1)疾病分数计算处理的概要
图6是实施方式所涉及的超声波诊断装置1100中的疾病分数计算处理的流程图。
i)评价对象图像的判定
在步骤S1000中,超声波图像取得部2001从存储器1005所存储的多个帧的B模式图像信号4001以及多普勒模式图像信号4002中按每个帧读出B模式图像信号4001以及多普勒模式图像信号4002。
接着,在步骤S1001中,关节检测部3001将B模式图像信号4001作为输入,按超声波图像信号的每个帧判定在超声波图像信号中是否包含表示关节部位的对象图像部分。在包含表示关节部位的对象图像部分的情况下,前进至下一个步骤,在除此之外的情况下前进至步骤S1006。
接着,在步骤S1002中,运动噪声检测部3003判定1帧的超声波图像中的多普勒模式图像信号4002是否是由运动噪声引起的图像信号。在多普勒模式图像信号4002不是由其引起的图像信号的情况下,前进至下一个步骤,在除此之外的情况下前进至步骤S1006。
接着,在步骤S1003中,压迫检测部3002判定超声波探头1001是否没有压迫被检测体的体表。在超声波探头1001没有压迫被检测体的体表的情况下,前进至下一个步骤,在除此之外的情况下前进至步骤S1006。
即,技巧判定部3004在图6的步骤S1001的判定结果为〔有关节〕、步骤S1002的判定结果为〔没有运动噪声〕、步骤S1003的判定结果为〔没有压迫〕的情况下,判定为对象帧的图像是在适当的技巧下得到的图像,将〔执行定量化〕输出至形态定量化部2003A和炎症定量化部2003B。
ii)疾病分数的计算
在步骤S1004中,疾病分数计算部2003计算疾病分数。形态定量化部2003A将肿胀分数输出至存储器1005,存储器1005对其进行存储。炎症定量化部2003B将炎症分数输出至存储器1005,存储器1005对其进行存储。
iii)疾病分数的选择
在步骤S1005中,选择部2004将存储器1005所存储的疾病分数(肿胀分数、炎症分数)作为输入,基于规定的数值处理而选择1个以上的最佳的疾病分数。
在第2次以后实施关于作为检查对象的手指关节的步骤S1000~S1005的处理的情况下,选择部2004从多个疾病分数中基于规定的数值处理而选择最佳的疾病分数,该多个疾病分数是根据已进行检查且存储在存储器1005中的多个帧的超声波图像信号而取得的多个疾病分数。
由此,能够不依赖于进行检查时的检查者的技巧或主观,根据基于规定的数值处理而选择的最佳的疾病分数和作为其计算的基础的、与最佳的疾病分数对应的超声波图像信号的帧,对疾病的程度进行评价。
另外,在初次进行关于对象的手指关节的步骤S1000~S1005的处理的情况下,仅在存储器1005中存储根据1个帧的超声波图像信号而取得的疾病分数。在该情况下,选择部2004选择存储器1005所存储的1个疾病分数。
iv)疾病分数等的显示
在步骤S1006中,显示控制部1006使用存储器1005所存储的疾病分数、和在S1005中基于规定的数值处理而选择的最佳的疾病分数,生成显示画面,在外部连接的显示器1008上进行显示。也可以基于在S1005中所选择的疾病分数、和作为其计算的基础的超声波图像信号,生成显示画面,在外部连接的显示器1008上进行显示。
v)检查的结束
在步骤S1007中,判定在对于1个手指关节的检查中是否对全部的超声波图像的帧完成了评价。在对全部的超声波图像的帧完成了评价的情况下,结束处理。进而,在存在进行评价的超声波图像的帧的情况下,返回至S1001。
下面,说明各步骤中的超声波诊断装置1100的动作。
(2)关节检测的处理
使用附图说明关节检测部3001中的处理。
图7是表示使用超声波诊断装置1100将手指关节作为测定对象而拍摄的超声波图像的示意图,并且是将振动元件列与指的纵向方向平行地配置,对手指关节进行超声波扫描而取得的B模式图像。如图7所示,在B模式图像中,绘制了指的关节部位,包含表示骨的图像部分1、2、表示皮肤的图像部分3、表示肌腱的图像部分4、表示关节囊的图像部分5。骨或皮肤或肌腱是相对硬的组织,因此,表示骨的图像部分1、2和表示皮肤的图像部分3和表示肌腱的图像部分4在超声波图像上被以高亮度的方式绘制。超声波的大部分在骨表面被反射,因此,骨的内部没有被绘制而仅以高亮度绘制相当于骨表面的骨皮质的部分。与表示骨的图像部分1、2、或表示皮肤的图像部分3相比,表示关节腔的图像部分6被以低亮度的方式绘制。此外,滑膜和软骨部分几乎不具有亮度值而不被显示。
因此,在表示关节的超声波图像中被以相对高的亮度绘制的组织是皮肤、肌腱、骨表面。超声波诊断装置的关节检测部3001从B模式图像信号中检测表示关节囊的图像部分5、表示骨的图像部分1、2,从而确定表示两者所包围的关节腔的图像部分6。然后,确定表示关节部位的图像部分,该表示关节部位的图像部分由表示关节囊的图像部分5、表示骨的图像部分1、2、表示关节腔的图像部分6构成。
图8是表示关节检测处理的动作的流程图。
在此,关节检测部3001经由超声波图像取得部2001而取得存储器1005所存储的1个帧的B模式图像信号,进行表示关节部位的图像的搜索。
首先,在步骤S2001中,关节检测部3001按每个关节检测窗302计算表示关节图案可能性的评价值,搜索评价值成为最大的位置。评价值成为最大的关节检测窗302的位置的搜索如下面那样被进行。
图9(a)、图9(b)是表示用于关节搜索处理的关节检测窗的示意图。在B模式图像301中,设定关节检测窗302,关节检测窗302向扫描线303的方向移动,通过对各位置中的关节检测窗302内的B模式图像的部分进行搜索处理,从而进行B模式图像301整体的搜索处理。此外,为了能够检测各种大小的关节,也可以代替B模式图像301,而生成对B模式图像301进行了放大或者缩小的B模式图像304,并在B模式图像304内通过关节检测窗302进行扫描。
图10是超声波诊断装置1100中的关节搜索处理中的动作说明图。如图10所示,在实施方式所涉及的关节搜索处理中,使用的模板匹配法,即:在B模式图像301上设定关节检测窗302,与表示关节的典型的图像图案的模板305比较,将相似度(误差值、相关值等)作为表示关节检测窗302内的图像部分的关节图案可能性的评价值进行计算。对于模板305,能够使用例如表示关节的平均的图像图案的B模式图像。使关节检测窗302向扫描线303的方向移动,关于各位置中的关节检测窗302内的B模式图像的部分进行搜索处理。以关节检测窗302扫描B模式图像301整体的方式使关节检测窗302移动,并且重复该处理,确定在评价值成为最大的位置上的关节检测窗302内的图像部分。
图11是表示基于模板匹配法的关节检测处理的细节的流程图,且是表示图8的步骤S2001中的处理的一例的图。如图10所示,在B模式图像301上的左上角部上将关节检测窗302设定为初始位置。计算表示关节的典型的图像图案的模板305和关节检测窗302内的B模式图像部分的相似度(S20011)。按每个像素计算模板305和关节检测窗302内的B模式图像部分的亮度信息的误差值,在关节检测窗302中包含的全部像素中进行总计,从而计算相似度。或者,计算模板305和关节检测窗302内的B模式图像部分的亮度信息的相关值,从而求出相似度。
然后,在关节检测窗302位于初始位置的情况下将计算出的相似度设为最大值,在除此以外的情况下,计算出的相似度与相似度的最大值进行比较(S20012)。
在计算出的相似度比最大值更大的情况下,将新计算出的相似度设为最大值,对关节检测窗302的位置进行保存(S20013)。在除此之外的情况下前进至S20014。
接着,在B模式图像301上使关节检测窗302的位置如扫描线303所示那样向X方向或者Y方向移动微小量,计算各位置中的相似度。搜索范围是B模式图像301全体,当在B模式图像301的全部区域中完成了相似度的计算的情况下(S20015),将在S20013中保存的关节检测窗302的位置确定为相似度为最大的关节检测窗302的位置(S20016)。
除了上述的模板匹配法之外,评价值也可以通过例如由机器学习法得到的使用了Haar-Like滤波器的关节图案/非关节图案的识别器来取得。图12是表示使用机器学习法的关节检测器的一例的示意图。在使用了机器学习法的关节搜索处理中,代替表示关节部位的典型的图像图案的模板305,使用由多个Haar-Like滤波器306a~g构成的关节检测器306,进行与图11所示的流程相同的处理。各Haar-Like滤波器306a~g是检测亮度变化的滤波器,对于在检测窗内的其位置/大小/亮度变化的灵敏度/权重能够通过Adaboost那样的机器学习法而得到。该各滤波器的输出值在S20016中确定关节位置时被计算,将各滤波器的输出值和通过学习而得到的权重的线性组合的和作为相似度而进行计算。
返回至图8,在步骤S2002中,关节检测部3001对评价值的最大值和阈值的大小进行比较。关节检测部3001在评价值的最大值为阈值以上的情况下,判定为〔有关节〕(步骤S2003),在下面的情况下,判定为〔没有关节〕(步骤S2004),将判定结果输出至技巧判定部3004。此时,关节检测部3001将表示评价值为最大的关节检测窗302的位置的信息输出至技巧判定部3004。
(3)运动噪声检测的处理
在图6的步骤S1002中,运动噪声检测部3003判定1帧的超声波图像中的多普勒模式图像信号是否是由运动噪声引起的信号。在本说明书中,运动噪声是指操作者较大地移动超声波探头1001,从而因不同于血流的原因而产生的多普勒信号。
当在1帧的多普勒模式图像信号中多普勒信号本身没有出现的情况下,判定为〔没有运动噪声〕。
当在1帧的超声波图像中多普勒信号出现的情况下,通过
1)帧的B模式图像信号和前一帧的B模式图像信号之间的每个像素的亮度变动值(以后,设为“帧间亮度差分”)、
2)1帧的B模式图像信号中的出现多普勒信号的区域(以后,设为“多普勒信号检测区域”)占据的面积比、
3)1帧的B模式图像信号内的高亮度区域中的多普勒信号检测区域占据的面积比,从而判定存在运动噪声的可能性。
图13是运动噪声检测的流程图。
首先,在步骤S3001中,运动噪声检测部3003对帧的B模式图像信号和前帧的B模式图像信号之间的每个像素的亮度变动值进行计算。亮度变动值是帧间的相关值,若操作者较大地移动超声波探头则帧间的相关值变低。也可以代替相关值,使用像素值的误差和。
接着,在S3002中,运动噪声检测部3003对1帧的B模式图像信号中的多普勒信号检测区域占据的面积比进行计算。具体而言,计算1帧的多普勒模式图像信号中的多普勒信号出现的像素数相对于全部像素数的比率。
图14(a)、图14(b)是表示帧的B模式图像信号中的产生了运动噪声的状态的一例的示意图。在图14(a)中,B模式图像301中的较大的范围内存在多普勒信号检测区域301A。多普勒信号虽是表示血流流速的信号,但原理上是动态地进行捕捉的信号,若操作者较大地移动超声波探头1001,则如图14(a)所示那样,多普勒信号检测区域301A占据B模式图像信号301整体的较大的范围。因此,在B模式图像301中多普勒信号检测区域301A的面积比比较大的情况下,设为产生了运动噪声。
接着,在S3003中,运动噪声检测部3003对1帧的B模式图像信号内的高亮度区域中的多普勒信号检测区域占据的面积比进行计算。具体而言,计算在1帧的B模式图像信号中的示出规定的亮度的像素占据的区域中出现多普勒信号的像素的比率。在B模式图像中,高亮度的区域是被检测体中的组织较硬的区域,骨等符合于此。在这样的区域中,不存在血流。因此,如图14(b)所示,当在B模式图像301中的表示高亮度的区域的骨的图像部分1中存在多普勒信号检测区域301A的一部分的情况下,设为产生了运动噪声。
接着,在S3004中,运动噪声检测部3003计算运动噪声判定的评价值,从而与阈值比较大小。若1)将B模式图像信号中的帧间亮度差分设为Mx,2)将1帧的B模式图像信号中的多普勒信号检测区域占据的面积比设为My,3)将1帧的B模式图像信号内的高亮度区域中的多普勒信号检测区域占据的面积比设为Mz,将ma、mb、mc设为常数,则运动噪声判定的评价值p被通过[数1]而赋予。
【数1】
p=ma·Mx+mb·My+mc·Mz
Mx、My、Mz使用各自的最大值和最小值,也可以归一化为0~1。评价值p越大运动噪声的可能性越高。运动噪声检测部3003在运动噪声判定评价值p为阈值以上的情况下,判定为〔有运动噪声〕,在为阈值以下的情况下,判定为〔没有运动噪声〕,并将判定结果输出至技巧判定部3004。
(4)体表压迫度检测的处理
在图6的步骤S1003中,压迫检测部3002判定超声波探头1001是否正在压迫被检测体的体表。若在取得超声波图像时被检测体的体表正在被超声波探头1001压迫,则存在被检测体中的新生血管被压迫,在超声波图像中不出现炎症反应的可能性。因此,在被检测体的体表正在被超声波探头1001压迫的情况下所取得的超声波图像不适合作为风湿疾病的评价对象图像,从评价对象中除去。
步骤S1003中的具体的判定基于以下情况进行:在B模式图像信号中,在被检测体的体表和超声波探头1001的表面之间,是否存在表示满足规定的基准的大小的超声波凝胶介质层的图像部分。
图15是说明体表压迫度检测的处理的图。判定在位于B模式图像301的上端边缘部分的超声波探头1001和表示被检测体体表的皮肤的图像部分3之间,是否存在表示超声波凝胶介质层的图像部分7(斜线部分)。表示超声波凝胶介质层的图像部分7由于成为表示纯黑的低亮度且低方差,因此根据在B模式图像301的上端边缘部分是否存在低亮度且低方差区域,能够判定是否存在表示超声波凝胶介质层的图像部分7。因此,如图15所示,在B模式图像301的上端边缘部分邻近的预定的范围中设定矩形形状的凝胶判定区域307,对凝胶判定区域307内的像素的亮度的平均和方差进行计算,与阈值进行比较。作为位于上端边缘的矩形区域的凝胶判定区域307的范围Y7优选为设定为例如距离上端部3~5毫米的范围。有时在B模式图像301的上端部由于超声波探头1001内的超声波的多重反射而具有亮度。能够通过除去距离上端部不到3毫米的范围从而避免多重反射的影响,并且能够更准确地判定表示超声波凝胶介质层的图像部分7的存在与否。
图16是体表压迫检测处理的流程图。
在步骤S4001中,压迫检测部3002计算B模式图像中的凝胶判定区域307的亮度的平均值和方差。根据方差,提取超声波凝胶介质层的亮度分布均匀那样的特征。
接着,在步骤S4002中,计算凝胶判定的评价值,与阈值进行比较。若将凝胶判定区域307的亮度的平均值设为Gx,将方差设为Gy,将gd、ge设为常数,则凝胶判定评价值q被通过[数2]而赋予。
【数2】
q=ga·Gx+gb·Gy
评价值q越小,存在超声波凝胶介质层的可能性越高。若凝胶判定评价值q为阈值以下,则由于存在一定的超声波凝胶介质层,因此压迫检测部3002判定为〔没有压迫〕(步骤S4003)。另一方面,若为阈值以上,则由于不存在一定的超声波凝胶介质层而判定为〔有压迫〕(步骤S4004),将判定结果输出至技巧判定部3004。
(5)关于评价对象的判定
技巧判定部3004在图6的步骤S1001的判定结果为〔有关节〕、步骤S1002的判定结果为〔没有运动噪声〕、步骤S1003的判定结果为〔没有压迫〕的情况下,判定为对象帧的图像是在适当的技巧下被得到的图像,并将〔执行定量化〕输出至形态定量化部2003A和炎症定量化部2003B。一并将表示评价值为最大的关节检测窗302的位置的信息输出至形态定量化部2003A和炎症定量化部2003B。另一方面,在步骤S1001的判定结果为〔没有关节〕、步骤S1002的判定结果为〔有运动噪声〕、步骤S1003的判定结果为〔有压迫〕的其中一个的情况下,技巧判定部3004判定为对象帧的图像不是在适当的技巧下被得到的图像,将〔中止定量化〕输出至形态定量化部2003A和炎症定量化部2003B,结束处理。
由此,技巧判定部3004能够更准确地判定对象帧的超声波图像信号是否是在适当的技巧下被得到的图像,能够更加准确地选择用于计算疾病分数的评价对象帧。
(6)肿胀分数(GS)的计算
在图6的步骤S1004中,疾病分数计算部2003计算疾病分数。图17是疾病分数计算处理的流程图。
形态定量化部2003A根据B模式图像的关节腔的大小/亮度/骨侵蚀的程度而计算疾病分数(肿胀分数)。
i)关节腔的大小/亮度的计算
首先,形态定量化部2003A基于表示评价值为最大的关节检测窗302的位置的信息,确定表示关节部位的B模式图像301的部分(S5001),从关节检测窗302内的B模式图像的部分中检测表示骨的图像部分(S5002)。
图18是说明了检测表示骨的图像部分1、2的处理的示意图。如上所述,骨是相对硬的组织,因此,在超声波图像上也被以高亮度绘制。超声波的大部分在骨表面被反射,因此,骨的内部没有被描绘而仅以高亮度描绘相当于骨表面的骨皮质的部分。
表示骨的图像部分1、2存在于比关节检测窗302的垂直中心位置308更下方的关节检测窗302内。因此,将比关节检测窗302的垂直中心位置308更下方的区域302A内的B模式图像的部分设为搜索范围。将以关节检测窗302的垂直中心位置308为起点朝向下方的方向设为搜索方向309(图中箭头方向),检测搜索方向309的亮度变化(S5002)。通过在检测亮度从低亮度成为高亮度的方向上进行搜索处理,从而能够降低边界的误检测。B模式图像301中,表示骨的图像部分1、2为高亮度,因此使用动态轮廓模型(SNAKES)等,能够对表示骨的图像部分1、2和周围的图像部分的边界1A、2A进行检测。此时,如图18所示,所检测到的由亮度变化表示的边界lB、2B在搜索方向309上急剧地变化,且搜索方向309上的边界位置的变化在与搜索方向309垂直的方向上平稳地变化的情况下,能够确定为该边界1B、2B是表示骨的图像部分1、2的边界lA、2A。
接着,形态定量化部2003A检测表示关节囊的边界的图像部分。
图19是说明了检测表示关节囊的边界的图像部分5的处理的示意图。在关节检测窗302内的表示骨的图像部分1、2的上方,存在表示关节囊的图像部分5。因此,关节检测容302中,将位于比表示骨的图像部分1、2更上方的区域内的B模式图像的部分设为搜索范围。在此,也能够通过在检测亮度从低亮度成为高亮度的方向上进行搜索处理,从而降低边界的误检测。将以关节检测窗302的表示骨的图像部分1、2和周围的图像部分的边界1A、2A为起点朝向上方的方向设为搜索方向309(图中箭头方向),检测搜索方向309上的亮度变化(S5003)。B模式图像301中,位于表示关节囊的图像部分5上部的表示脂肪层或肌肉层的图像部分的亮度比表示关节腔的图像部分6亮度更高,使用动态轮廓模型(SNAKES)等,能够对表示关节囊的图像部分5和表示关节腔的图像部分6的边界5A进行检测。此时,如图19所示,所检测到的边界5A在搜索方向309上急剧地变化,且搜索方向309上的边界位置的变化在与搜索方向309垂直的方向上平稳地变化的情况下,能够确定为该边界5A是表示关节囊的图像部分5。此外,也可以通过基于亮度值的相似度的区域扩张法等的图像处理,确定图像部分5。
接着,形态定量化部2003A计算表示关节腔的图像部分6的大小(S5004)。在表示关节腔的图像部分6较大的情况下,能够评价为疾病的程度较高。表示关节腔的图像部分6的大小能够作为骨的边界lA、2A和表示关节囊的图像部分5(5A)的垂直方向的距离或面积而进行计算。计算沿着与搜索方向309垂直的方向的平均值以及最大值。
接着,形态定量化部2003A对表示关节腔的图像部分6的平均亮度进行计算(S5005)。在表示关节腔的图像部分6的平均亮度较低的情况下,能够评价为疾病的程度较高。表示关节腔的图像部分6的平均亮度能够通过对表示关节腔的图像部分6中包含的各像素示出的亮度进行平均而计算出,所述图像部分6由骨的边界lA、2A和表示关节囊的图像部分5(5A)包围。表示关节腔的图像部分6是指在图19中由表示关节腔的图像部分6的外周边缘6A包围的范围。此外,也可以代替计算平均亮度,而计算表示关节腔的图像部分6中的亮度的中值。
ii)骨侵蚀的程度的计算
接着,形态定量化部2003A计算骨侵蚀7的程度(S5006)。在骨侵蚀7的程度较高的情况下,能够评价为疾病的程度较高。
图20是说明了检测骨侵蚀7的程度的处理的示意图。骨侵蚀7的程度能够根据与表示关节腔的图像部分6之间的骨1、2的边界lA、2A的平滑度而进行评价。具体而言,如图20所示,从边界lA、2A的最深部起向左右分开在沿着边界lA、2A的搜索方向309上进行搜索处理,从而进行评价。计算边界lA、2A的拟合函数而根据拟合误差来评价平滑度。
图21是示出了正常的骨表面的边界上的拟合的一例的示意图,图22是示出了骨侵蚀的骨表面的边界上的拟合的一例的示意图。分别示出了表示骨边界的图像1A和其拟合曲线310的一例。在图22所示的骨侵蚀的骨表面的边界,与图21所示的正常的骨表面的边界相比,拟合误差310A变大,能够通过沿着拟合曲线310累积拟合误差310A从而评价骨侵蚀的程度。另外,用于拟合的函数是任意的,但是优选为以不与骨侵蚀的部分拟合的方式来进行阶数的选择。
此外,也可以设为以下结构:代替骨边界,对骨边界进行一阶微分或者二阶微分,对在沿着该骨表面的方向上的变化进行检测,从而检测骨表面的平滑度。能够排除骨表面的个体间的形状变化的影响,关于骨表面的平滑度进行评价。
iii)肿胀分数的计算
接着,形态定量化部2003A计算作为疾病分数之一的肿胀分数(S5007)。若将表示关节腔的图像部分6的大小设为GSx,将表示关节腔的图像部分6的亮度设为GSy,将骨侵蚀的程度设为GSz,将GSa、GSb、GSc设为常数,则肿胀分数被通过[数3]而赋予。
【数3】
GS=gsa·GSx+gsb·GSy+gsc·GSz
根据发明者们的研讨,使用与疾病的进展有较高的关联性的3个参数GSx、GSy、GSz进行了线性组合的肿胀分数(GS),从而能够更准确地评价疾病的程度。GSx、GSy、GSz可以使用各自的最大值和最小值,归一化为0~1。然后,如上所述,形态定量化部2003A将肿胀分数输出至存储器1005,存储器1005对其进行存储(步骤S1004)。形态定量化部2003A将示出用于表示关节腔的图像部分6的位置信息输出至炎症定量化部2003B。
通过上述结构,根据评价对象帧的B模式图像信号中的对象图像部分的信号,能够按每个评价对象帧计算将疾病的程度定量化后的疾病分数,能够提高评价的客观性。
(7)炎症分数(PD)的计算
在图17中,炎症定量化部2003B将来自形态定量化部2003A的示出用于表示关节腔的图像部分6的位置信息作为输入,基于多普勒模式图像信号而计算疾病分数之一的炎症分数(S5008)。
在S5005中,基于形态定量化部2003A确定了的表示关节腔的图像部分6的范围,炎症定量化部2003B将表示关节腔的图像部分6设定为关心区域(ROI:Region of Interest,兴趣区域)。然后,对关心区域占据的面积(PDy)进行计算。进而,炎症定量化部2003B计算在作为关心区域的表示关节腔的图像部分6中包含的各像素中被检测到多普勒信号的像素占据的面积(PDx)。在像素中是否被检测到多普勒信号的判定能够基于是否被检测到超过规定的阈值的多普勒信号而适宜地进行设定。炎症分数作为对面积(PDx)除以面积(PDy)的值而被通过[数4]而赋予。
【数4】
PD = PDx PDy
此时,也可以是如下结构:将包含表示关节腔的图像部分6的至少一部分的例如矩形形状的固定尺寸的图像部分设为关心区域,根据骨的检测结果而进行设定,对该关心区域的面积(PDy)和关心区域中的多普勒信号占据的面积(PDx)进行计算。该情况下,面积(PDy)成为固定值。
根据发明者们的研讨,使用基于是否检测出像素中的多普勒信号的炎症分数(PD),能够更有效且更准确地评价疾病的程度。如上所述,炎症定量化部2003B将炎症分数输出至存储器1005,存储器1005对其进行存储(步骤S1004)。
通过上述结构,根据评价对象帧的B模式图像信号和多普勒模式图像信号中的对象图像部分的信号,能够按每个评价对象帧而对定量化了疾病的程度的疾病分数进行计算,能够提高评价的客观性。
(8)关于疾病分数的选择、所选择的疾病分数的显示以及保存
如上所述,在图6的步骤S1005中,选择部2004将存储器1005所存储的疾病分数(肿胀分数、炎症分数)作为输入,基于规定的数值处理,对肿胀分数和炎症分数分别选择1个以上的最佳的疾病分数。在本实施方式中,选择部2004设为选择存储器1005所存储的疾病分数的最大值作为最佳的疾病分数的结构。
可是,也可以设为将表示疾病的程度平均的状态的疾病分数的平均值、或者表示疾病的程度居中的状态的疾病分数的中值选择作为最佳的疾病分数的结构。此外,所选择的疾病分数也可以是与从多个疾病分数的最大值、平均值、以及中值中选择的至少1个相同的疾病分数、或者与该数值最近的疾病分数的其中一个。由此,能够基于医生或病院的检查方针、疾病的状态、被检测者的特性等各种条件,适宜地设定基于用于选择最佳的疾病分数的规定的数值处理的选择基准。
此外,也可以设为具有多个基于规定的数值处理的选择基准,并选择来进行检查的结构。能够基于疾病的状态、被检验者的特性等各种条件来选择选择基准从而进行检查。该情况下,优选设为在检查中用于疾病分数的选择的选择基准被输出至存储器1005并进行存储的结构。
此外,如上所述,基于规定的数值处理而选择了的最佳的疾病分数被输出至存储器1005并进行存储。此时,进而,也可以设为用于疾病分数的选择的数值处理或者基于该数值处理的选择基准也被输出至存储器1005并进行存储的结构。通过所涉及的结构,能够使用与在过去计算出的疾病分数相同的选择基准,进行疾病分数的选择,并且在定期检查等的进展观察时与过去的评价结果的比较变得容易。
另外,在对疾病分数进行校正的情况下,通过操作单元1009取得操作者的校正信息,存储到存储器1005。在关节囊或多普勒信号被进行了校正的情况下,在形态定量化部2003A和炎症定量化部2003B中,基于校正信息而重新计算疾病分数。
2.超声波诊断装置1100的显示画面
(1)操作画面
接着,说明超声波诊断装置1100的操作画面。
下面,假设对超声波诊断装置1100的输入是对在显示器1008的显示画面上配置的触摸面板即操作单元1009的触摸操作,与显示画面所显示的图像一并进行说明。
在检查之前,操作者对手指关节的检查顺序进行注册。图23是用于注册检查顺序的超声波诊断装置1100的操作输入画面。对表示关节位置的图标201进行了触摸的顺序作为对图标201表示的关节进行检查的检查顺序而被注册。在图标201内显示用于表示检查顺序的数字。此外,在想要改变初始设定的顺序的情况下,在图标201内显示了用于表示检查顺序的数字的状态下,进一步对图标201进行了触摸的顺序被更新注册为检查顺序。
接着,说明对于1个手指关节的检查中的超声波诊断装置1100的操作画面。
图24是测量开始前的显示画面。如图24所示,在显示画面内显示B模式图像显示框101、多普勒模式图像显示框102、肿胀分数显示框103、炎症分数显示框104、帧号码106。此外,显示用于表示帧位置的滑块107、检查中的被选择的肿胀分数的例如最大值显示框109、检查中的被选择的炎症分数的例如最大值显示框110、测量开始以及停止按钮111、确定疾病分数并且前进至下一个关节的检查的下一个按钮112。
在该状态下,若操作者按压测量开始以及停止按钮111,则开始图6的S100l中的测量,跳转至图25所示的测量中的显示画面,并显示消息113。然后,在B模式图像显示框101中实时地显示所取得的B模式图像,在多普勒模式图像显示框102中实时地显示所取得的多普勒模式图像。此外,在各图像中,诊断装置检测出的关节囊边界105被重叠显示。
在S1001中判定为图像中不包含关节部位的情况下,如图26所示,显示在图像中不包含关节部位的警告消息113。此外,在S1002中判定为多普勒模式图像是由运动噪声引起的图像的情况下、或在S1003中判定为超声波探头1001正在压迫被检测体体表的情况下,都不计算疾病分数,而是如图27、图28所示,显示警告消息113。
然后,若在操作者对超声波探头1001进行了移动的结果捕捉到疾病进展的关节断面,则如图29所示,显示用于表示测量出的疾病分数的显示画面。具体而言,在多普勒模式图像上显示表示存在血流的区域的多普勒信号存在区域114,作为最佳的疾病分数的一例,更新肿胀分数的最大值和炎症分数的最大值。此外,分别显示用于表示各疾病分数(最大肿胀分数和最大炎症分数)成为最大的帧的图标108。进而,在S1004中计算疾病分数,在肿胀分数显示框103中显示肿胀分数,在炎症分数显示框104中显示炎症分数。进而,在S1005中,选择最大疾病分数作为最佳的疾病分数的一例,在肿胀分数的最大值显示框109中显示肿胀分数的最大值,在检查中的炎症分数的最大值显示框110中显示炎症分数的最大值。
在该状态下,若操作者按压测量开始以及停止按钮111则停止测量,如图30所示,显示用于表示测量完成的显示画面。进而,通过S1006,显示关于肿胀分数和炎症分数的最大疾病分数、各疾病分数为最大的帧的超声波图像、表示进行了检查的最终帧的图标图像116。
若触摸表示最终帧的图标116,则显示最终帧。由此,检查中,在通过诊断装置而新取得违反操作者的意思的不适当的帧,并且选择为评价对象帧的情况下,操作者能够变更用于评价的帧。此外,显示疾病分数的校正按钮115,操作者在没有得到适当的疾病分数的情况下,可以通过按下疾病分数的校正按钮115从而对测量结果进行校正。此外,在多普勒信号存在区域114不是对因炎症产生的血流,而是对血管捕捉的区域的情况下,可以删除多普勒信号存在区域114。在没有校正的需要的情况下,若按压下一个按钮112则跳转至下一个关节的测量。
图31是校正操作中的显示画面。若有想要删除的多普勒信号存在区域114,则操作者用指117对该区域进行触摸。由此,如图32所示,从多普勒模式图像中删除该多普勒信号存在区域114,显示用于表示删除的图标118,伴随删除而更新疾病分数。此外,虽然未图示,但是关节囊边界105也能够通过相同的方法进行校正,伴随校正而更新疾病分数。
在S1007中,在关于1个手指关节中的全部的超声波图像完成了评价的情况下,通过操作者的操作而完成检查,如图30所示,显示用于表示测量完成的显示画面。或者,继续显示如图30所示的显示画面。关于相同的手指关节进行新的检查的情况下,返回至S1001,显示如图25所示的显示画面。
(2)检查结果的显示画面
在关于全部的手指关节中的全部的超声波图像完成了评价的情况下,进行用于表示检查结果的显示。图33是示出了各手指关节中的检查结果的显示画面。如图33所示的显示画面可以在全部的手指关节的检查完成了后进行显示,也可以在每次按压了下一个按钮112后进行显示。
在图33中,示出用于表示各手指关节的疾病分数的图标202,即左侧的数字表示肿胀分数,右侧的数字表示炎症分数。也可以是表示疾病分数较高的关节的疾病分数的图标202被彩色显示等、被增强显示的结构。
若对表示手指关节的疾病分数的图标202的1个进行触摸,则显示被计算出该手指关节中的疾病分数的最佳值的帧的超声波图像和示出了疾病分数的表示测量完成的显示画面(图30)。
若长按用于表示各手指关节的疾病分数的图标202之一,则如图34所示,显示对疾病分数进行了图形化显示的检查结果,所述疾病分数根据手指关节中的多个帧的超声波图像而计算出。图34是显示了超声波探头1001的每个倾斜角度的疾患分数的显示方法的一例。显示与疾病分数的角度相关的积分值和最大值203、每个角度的肿胀分数204、每个角度的炎症分数205。这样,通过按每个角度而保留记录,从而,在定期检查等的进展观察中,容易对根据在相同倾斜角度上取得的超声波图像计算出的疾病分数之间进行比较。
<效果>
上面,如所说明那样,实施方式所涉及的超声波诊断装置1100是从经由超声波探头1001从被检测体所取得的多个帧的超声波图像信号中,选择1个以上的超声波图像信号的帧,并在显示器1008上进行显示的超声波诊断装置1100,所述超声波检测装置的特征在于,具备控制电路1010,该控制电路1010具备:超声波图像取得部2001,取得多个帧的超声波图像信号;评价对象判定部2002,分析多个帧的超声波图像信号,在各帧中包含用于表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧;疾病分数计算部2003,根据评价对象帧的超声波图像信号中的对象图像部分的超声波图像信号,对定量化了疾病的程度的疾病分数进行计算;选择部2004,对所计算出的疾病分数,基于规定的数值处理而选择1个以上的疾病分数;以及显示控制部1006,使显示器1008对所选择出的疾病分数以及与该疾病分数对应的帧的超声波图像进行显示。
通过上述结构,在从多方向得到的超声波图像中,能够根据基于规定的数值处理而选择了的断面的超声波图像,评价疾病的程度。
以往,存在如下的课题:根据从被检测体的哪个断面取得超声波图像而进行评价,评价结果依赖于检查者的技巧或主观。对此,通过上述结构,能够选择在适当的技巧下得到的超声波图像而客观地评价疾病的程度,因此,能够降低疾病评价中的评价结果的检测者依赖性。
<其他的变形例>
上面,说明了各实施方式所涉及的超声波诊断装置。另外,本发明并不限定于各实施方式,能够基于实施方案而适宜地进行变更。
在上述实施方式中,设为了如下的结构:技巧判定部3004在图6的步骤S1001的判定结果为〔有关节〕、步骤S1002的判定结果为〔没有运动噪声〕、步骤S1003的判定结果为〔没有压迫的〕的情况下,判定为对象帧的图像是在适当的技巧下被得到的图像,将〔执行定量化〕输出至形态定量化部2003A和炎症定量化部2003B。然后,在除此之外的情况下,技巧判定部3004设为了如下结构:判定为对象帧的图像不是在适当的技巧下被得到的图像,输出〔中止定量化〕,结束处理。
可是,也可以设为如下的结构:技巧判定部3004在步骤S1001的判定结果为〔有关节囊〕时,判定为对象帧的图像是在适当的技巧下被得到的图像。
或者,技巧判定部3004也可以在步骤S1001的判定结果为〔有关节囊〕,且进一步满足步骤S1002的判定结果为〔没有压迫〕或者步骤S1003的判定结果为〔没有运动噪声的〕的至少一个时,判定为是在适当的技巧下被得到的图像。
由此,技巧判定部3004能够更简便地判定对象帧的超声波图像信号是否是在适当的技巧下被得到的图像,从而更加迅速地进行用于计算疾病分数的评价对象帧的选择。
此外,在上述实施方式中,疾病分数设为了[数3][数4],但是只要是与风湿疾病关联的分数,则并不限定于此。
此外,在上述实施方式中,作为定量化了疾病的程度的疾病分数的一例,以风湿为例而说明了选择疾病分数并进行疾病的评价的方法。可是,能够使用本公开的方法的疾病的种类不限定于风湿,只要是根据超声波图像能够对疾病的程度进行定量化的疾病,则能够对其他的疾病进行应用。例如,可以将本公开的方法用于根据超声波图像对器官内的肿瘤的大小或肿瘤内的新生血管的面积比进行定量化,从而对癌症疾病的程度进行定量化。
在上述实施方式所涉及的超声波诊断装置1100中,设为了在超声波诊断装置1100内包含作为存储装置的存储器1005的结构,但是,存储装置不限定于此,也可以是半导体存储器、硬盘驱动器、光盘驱动器、磁存储装置等从外部连接到超声波诊断装置1100的结构。
此外,设为了超声波探头1001以及显示器1008被从外部连接到超声波诊断装置1100的结构,但是,也可以设为在超声波诊断装置1100内一体地具备这些的结构。
此外,也可以设为在超声波探头1001中设置例如角度传感器等的倾斜角度测定部件,测定超声波探头1001的倾斜角度,并与疾病分数一并记录到检查结果的结构。
此外,在上述实施方式中,超声波探头示出了多个压电元件被排列在一维方向上的超声波探头结构。可是,超声波探头的结构不限定于此,也能够使用例如多个压电变换元件以2维的方式被排列的超声波探头。在使用了以2维的方式被排列的超声波探头的情况下,通过分别使对压电变换元件提供的电压的定时或电压的值变化,能够控制要发送的超声波束的照射位置或照射方向。
此外,超声波探头也可以包含发送接收处理部的一部分的功能。例如,能够采用如下的结构:基于用于生成从发送接收处理部被输出的发送电信号的控制信号,在超声波探头内生成发送电信号,并将该发送电信号变换为超声波。一并将所接收到的反射超声波变换为接收电信号,在超声波探头内基于接收电信号而生成接收信号。
此外,各实施方式所涉及的超声波诊断装置中包含的各处理部,典型的是作为集成电路即LSI而被实现。这些处理部也可以单独地作为1个芯片,也可以包含一部分或者全部而作为1个芯片。
此外,在实施方式中,说明了各块由独立的硬件构成。可是,构成超声波诊断装置的各块不一定必须由独立的硬件而构成,也可以是例如,通过将各块根据需要设为一体的CPU以及软件,实现其功能的结构。
此外,超声波诊断装置的各功能块能够将各个的功能块的一部分或者全部的功能典型地以集成电路即LSI实现。这些功能块也可以独自地作为1个芯片,也可以包含一部分或者全部作为1个芯片。另外,LSI根据集成度的不同,有时也称为IC、系统LSI、Super LSI、Ultra LSI。
另外,集成电路化的方法不限定于LSI,也可以通过专用电路或者通用处理器来实现。也可以利用在LSI制造后能够进行编程的FPGA(FieldProgrammable Gate Array,现场可编程门阵列)、或能够对LSI内部的电路单元的连接或设定进行再次构成的可重构处理器(ReConfigurab1e Processor)。
进而,如果由于半导体技术的进步或者衍生的其他技术而出现置换LSI的集成电路化的技术,当然也可以使用该技术进行功能块的集成化。
此外,也可以通过CPU等的处理器执行程序,从而实现各实施方式所涉及的、超声波诊断装置的功能的一部分或者全部。
进而,本发明也可以是上述程序,也可以是记录了上述程序的非暂时性计算机可读取的记录介质。此外,当然,上述程序能够经由互联网等的传送介质而流通。
此外,方框图中的功能块的分割是一例,也可以将多个功能块作为1个功能块而实现,或者将1个功能块分割为多个,或者将一部分的功能转移到其他的功能块。此外,也可以单个硬件或者软件以并行或者时分的方式对具有类似的功能的多个功能块的功能进行处理。
此外,上述的步骤被执行的顺序是为了具体地说明本发明而示例的顺序,也可以是上述以外的顺序。此外,上述步骤的一部分也可以与其他的步骤同时(并行)地执行。
此外,也可以对各实施方式所涉及的超声波诊断装置、以及其变形例的功能之中至少一部分进行组合。
进而,对本实施方式实施了在本领域技术人员能够想到的范围内的变更的各种变形例,也被包含于本发明。
《总结》
上面,如所说明那样,本实施方式所涉及的超声波诊断装置是从经由超声波探头从被检测体所取得的多个帧的超声波图像信号中选择1个以上的超声波图像信号的帧,并在显示器上进行显示的超声波诊断装置,所述超声波诊断装置的特征在于,所述超声波检测装置具备控制电路,该控制电路具备:超声波图像取得部,取得所述多个帧的超声波图像信号;评价对象判定部,分析所述多个帧的超声波图像信号,在各帧中包含用于表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧;疾病分数计算部,根据所述评价对象帧的超声波图像信号中的所述对象图像部分的超声波图像信号,按每个所述评价对象帧对定量化了疾病的程度的疾病分数进行计算;选择部,对所述计算出的疾病分数,基于规定的数值处理而选择1个以上的疾病分数;以及显示控制部,使显示器对所述选择出的疾病分数以及与该疾病分数对应的帧的超声波图像进行显示。
通过所涉及的结构,在从多方向得到的超声波图像中,能够根据基于规定的数值处理而选择出断面的超声波图像,对疾病的程度进行评价,能够客观地对疾病的程度进行评价。因此,能够降低疾病评价中的评价结果的检测者依赖性。
此外,在其他的方式中,也可以是如下的结构:所述评价对象判定部进一步,在根据所述多个帧的超声波图像信号,检测出在取得该信号时没有所述超声波探头对所述被检测体体表的压迫,且检测出所述超声波图像信号的帧中的多普勒信号不是由运动噪声引起的信号时,判定为所述超声波图像信号的帧是评价对象帧。
通过所涉及的结构,能够更加正确地判定对象帧的超声波图像信号是否是在适当的技巧下被得到的图像,能够更加正确地进行用于计算疾病分数的评价对象帧的选择。
此外,在其他的方式中,也可以是如下的结构:所述规定的数值处理是如下的处理:从所计算出的所述疾病分数中,选择与从(a)表示疾病进展最严重的状态的疾病分数的最大值、(b)表示疾病的程度平均的状态的疾病分数的平均值、或者(c)表示疾病的程度居中的状态的疾病分数的中值中被选择出的至少1个数值相同的疾病分数,或者选择与该数值最近的疾病分数。
通过所涉及的结构,基于医生或病院的检查方针、疾病的状态、被检测者的特性等各种条件,能够对基于用于选择最佳的疾病分数的规定的数值处理的选择基准适宜地进行设定。
此外,在其他的方式中,也可以是如下的结构:所述超声波图像信号包含B模式图像信号以及多普勒模式图像信号。
此外,在其他的方式中,也可以是如下的结构:所述疾病是风湿。
此外,在其他的方式中,也可以是如下的结构:还具备:所述超声波探头;以及角度检测部,计算所述超声波探头的倾斜角度,所述显示控制部使所述显示器对所述疾病分数以及与该疾病分数对应的超声波图像信号的帧被取得时的所述超声波探头的倾斜角度进行显示。
通过所涉及的结构,通过按疾病分数的每个角度而保留记录,从而,容易对根据在相同倾斜角度上取得的超声波图像计算出的疾病分数之间进行比较。
此外,在其他的方式中,也可以是如下的结构:还具备:存储装置,保存所述多个帧的超声波图像信号,所述选择部进一步将所述选择出的疾病分数保存到所述存储装置。
通过所涉及的结构,在进展观察时,容易与根据在过去取得的超声波图像计算出的疾病分数进行比较。
此外,在其他的方式中,也可以是如下的结构:所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号,所述疾病分数计算部在将表示关节腔的图像部分的尺寸设为GSx,将表示关节腔的图像部分的平均亮度设为GSy,将骨侵蚀的程度设为GSz,将GSa、GSb、GSc设为常数时,对通过[数3]表示的GS进行计算。
【数3】
GS=gsa·GSx+gsb·GSy+gsc·GSz
通过所涉及的结构,根据评价对象帧的B模式图像信号中的对象图像部分的信号,能够按每个评价对象帧对定量化了疾病的程度的疾病分数进行计算,能够提高评价的客观性。
此外,在其他的方式中,也可以是如下的结构:所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号以及多普勒模式图像信号,所述疾病分数计算部在将占据表示关节腔的图像部分的一部分或者全部的关心区域的面积设为PDy,将所述关心区域中包含的像素中被检测出多普勒信号的像素占据的面积设为PDx时,对通过[数4]表示的PD进行计算。
【数4】
PD = PDx PDy
通过所涉及的结构,根据评价对象帧的B模式图像信号和多普勒模式图像信号中的对象图像部分的信号,能够按每个评价对象帧对定量化了疾病的程度的疾病分数进行计算,能够提高评价的客观性。
此外,在其他的方式中,也可以是如下的结构:所述存储装置中保存有在过去计算出的疾病分数以及基于所述规定的数值处理的选择基准,所述选择部参照所述过去的疾病分数,决定所述选择基准。
通过所涉及的结构,能够使用与在过去计算出的疾病分数相同的选择基准,进行疾病分数的选择,在定期检查等的进展观察时与过去的评价结果的比较变得容易。
此外,实施方式所涉及的超声波图像处理方法是从经由超声波探头从被检测体所取得的多个帧的超声波图像信号中选择1个以上的超声波图像信号的帧,并使显示器进行显示的超声波图像处理方法,所述超声波图像处理方法的特征在于,所述超声波图像处理方法具有:从保存了所述多个帧的超声波图像信号的存储装置读出所述多个帧的超声波图像信号的步骤;分析所述多个帧的超声波图像信号,在各帧中包含用于表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧的判定步骤;根据所述评价对象帧的超声波图像信号中的所述对象图像部分的超声波图像信号,对定量化了风湿的程度的疾病分数进行计算的步骤;对所述计算出的疾病分数,基于规定的数值处理而选择1个以上的疾病分数的选择步骤;以及使显示器对所述选择出的疾病分数以及与该疾病分数对应的帧的超声波图像进行显示的步骤。此外,在其他的方式中,也可以是一种计算机可读取的非暂时性的记录介质,所述记录介质记录了使计算机执行上述超声波图像处理方法的程序。
通过所涉及的结构,能够提供一种超声波图像处理方法,所述超声波图像处理方法能够选择在适当的技巧下得到的超声波图像,客观地评价疾病的程度,降低疾病评价中的评价结果的检测者依赖性。
《补充》
上面说明了的实施方式都是表示本发明的优选的具体例的实施方式。实施方式所示的数值、形状、材料、结构要素、结构要素的配置位置以及连接方式、步骤、步骤的顺序等是一例,并不限定本发明。此外,实施方式中的结构要素中,关于在表示本发明的最上位概念的独立权利要求中未记载的步骤,以构成更优选的方式的任意的结构要素而进行说明。
此外,为了容易理解发明,上述各实施方式中列举出的各图的结构要素的比例尺存在与实际的比例尺不同的情况。此外,本发明不根据上述各实施方式的记载而被限定,在不脱离本发明的主旨的范围内能够适宜地进行变更。
进而,在超声波诊断装置中,在基板上存在电路部件、引线等的构件,但是关于电气布线、电路,能够基于该技术领域中的通常的知识而对各种的方式进行实施,作为本发明的说明没有直接的关系,因此省略了说明。另外,上述示出了的各图是示意图,并不一定是严格地进行了图示的图。
产业上的利用可能性
根据本公开所涉及的超声波诊断装置以及超声波图像处理方法,通过设为上述结构,从而能够降低疾病评价中的评价结果的检测者依赖性。因此,在超声波诊断装置以及超声波图像处理方法中能够广泛地进行运用。

Claims (18)

1.一种超声波诊断装置,从经由超声波探头从被检测体所取得的多个帧的超声波图像信号中选择1个以上的超声波图像信号的帧,并在显示器上进行显示,其中,
所述超声波诊断装置具备控制电路,
该控制电路具备:
超声波图像取得部,取得所述多个帧的超声波图像信号;
评价对象判定部,分析所述多个帧的超声波图像信号,在各帧中包含用于表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧;
疾病分数计算部,根据所述评价对象帧的超声波图像信号中的所述对象图像部分的信号,按每个所述评价对象帧,计算将疾病的程度定量化后的疾病分数;
选择部,对所述计算出的疾病分数,基于规定的数值处理而选择1个以上的疾病分数;以及
显示控制部,使显示器显示所述选择的疾病分数以及与该疾病分数对应的帧的超声波图像。
2.如权利要求1所述的超声波诊断装置,其中,
所述评价对象判定部进一步在根据所述多个帧的超声波图像信号检测出在取得该信号时所述超声波探头没有对所述被检测体体表进行压迫,并且所述超声波图像信号的帧中的多普勒信号不是由运动噪声引起的信号时,判定为所述超声波图像信号的帧是评价对象帧。
3.如权利要求1所述的超声波诊断装置,其中,
所述选择部中的规定的数值处理,是从计算出的所述疾病分数中选择与从下列数值中选择的至少1个数值相同的疾病分数,或者与该数值最近的疾病分数的处理:(a)表示疾病进展最严重的状态的疾病分数的最大值、(b)表示疾病的程度平均的状态的疾病分数的平均值、或者(c)表示疾病的程度居中的状态的疾病分数的中值。
4.如权利要求1所述的超声波诊断装置,其中,
所述超声波图像信号包含B模式图像信号以及多普勒模式图像信号。
5.如权利要求1所述的超声波诊断装置,其中,
所述疾病是风湿。
6.如权利要求1所述的超声波诊断装置,其中,
所述超声波诊断装置还具备:所述超声波探头;以及角度检测部,计算所述超声波探头的倾斜角度,
所述显示控制部使所述显示器显示所述疾病分数以及与该疾病分数对应的超声波图像信号的帧被取得时的所述超声波探头的倾斜角度。
7.如权利要求1所述的超声波诊断装置,其中,
所述超声波诊断装置还具备:存储装置,保存所述多个帧的超声波图像信号,
所述选择部将所述选择的疾病分数保存到所述存储装置中。
8.如权利要求4所述的超声波诊断装置,其中,
所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号,
所述疾病分数计算部在将表示关节腔的图像部分的大小设为GSx,将表示关节腔的图像部分的平均亮度设为GSy,将骨侵蚀的程度设为GSz,将GSa、GSb、GSc设为常数时,计算由【数3】表示的GS,
【数3】
GS=gsa·GSx+gsb·GSy+gsc·GSz。
9.如权利要求4所述的超声波诊断装置,其中,
所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号以及多普勒模式图像信号,
所述疾病分数计算部在将占据表示关节腔的图像部分的一部分或者全部的关心区域的面积设为PDy,将所述关心区域中包含的像素中检测出多普勒信号的像素所占据的面积设为PDx时,计算由【数4】表示的PD,
【数4】
PD = PDx PDy .
10.如权利要求9所述的超声波诊断装置,其中,
所述PDy是固定值。
11.如权利要求7所述的超声波诊断装置,其中,
所述存储装置中保存有在过去计算出的疾病分数以及基于所述规定的数值处理的选择基准,
所述选择部参照所述过去的疾病分数,决定所述选择基准。
12.一种超声波图像处理方法,从经由超声波探头从被检测体取得了的多个帧的超声波图像信号中选择1个以上的超声波图像信号的帧,并使显示器进行显示,其中,
所述超声波图像处理方法具有:
从保存了所述多个帧的超声波图像信号的存储装置中读出所述多个帧的超声波图像信号的步骤;
分析所述多个帧的超声波图像信号,在各帧中包含用于表示关节部位的对象图像部分时,判定为该超声波图像信号的帧是评价对象帧的判定步骤;
根据所述评价对象帧的超声波图像信号中的所述对象图像部分的信号,计算将风湿的程度定量化后的疾病分数的步骤;
对所述计算的疾病分数,基于规定的数值处理而选择1个以上的疾病分数的选择步骤;以及
使显示器显示所述选择的疾病分数以及与该疾病分数对应的帧的超声波图像的步骤。
13.如权利要求12所述的超声波图像处理方法,其中,
在所述判定步骤中,进一步根据所述多个帧的超声波图像信号,检测出在取得了该信号时所述超声波探头没有对所述被检测体体表进行压迫,并且所述超声波图像信号的帧中的多普勒信号不是由运动噪声引起的信号时,判定为所述超声波图像信号的帧是评价对象帧。
14.如权利要求12所述的超声波图像处理方法,其中,
所述选择判定步骤中的规定的数值处理,是从计算出的所述疾病分数中选择与从下列数值中选择的至少1个数值相同的疾病分数,或者与该数值最近的疾病分数的处理:(a)表示疾病进展最严重的状态的疾病分数的最大值、(b)表示疾病的程度平均的状态的疾病分数的平均值、或者(c)表示疾病的程度居中的状态的疾病分数的中值。
15.如权利要求12所述的超声波图像处理方法,其中,
所述超声波图像信号的帧包含B模式图像信号以及多普勒模式图像信号。
16.如权利要求15所述的超声波图像处理方法,其中,
所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号,
所述疾病分数计算部在将表示关节腔的图像部分的大小设为GSx,将表示关节腔的图像部分的平均亮度设为GSy,将骨侵蚀的程度设为GSz,将GSa、GSb、GSc设为常数时,计算由【数3】表示的GS,
【数3】
GS=gsa·GSx+gsb·GSy+gsc·GSz。
17.如权利要求15所述的超声波图像处理方法,其中,
所述对象图像部分的超声波图像信号是表示关节腔的图像部分的B模式图像信号以及多普勒模式图像信号,
所述疾病分数计算部在将占据表示关节腔的图像部分的一部分或者全部的关心区域的面积设为PDy,将所述关心区域中包含的像素中检测出多普勒信号的像素所占据的面积设为PDx时,计算由【数4】表示的PD,
【数4】
PD = PDx PDy .
18.一种计算机可读取的非暂时性的记录介质,所述记录介质记录了使计算机执行权利要求12所述的超声波图像处理方法的程序。
CN201410415059.7A 2013-08-21 2014-08-21 超声波诊断装置、超声波图像处理方法 Pending CN104414685A (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2013-171036 2013-08-21
JP2013171036 2013-08-21
JP2014164845A JP2015061592A (ja) 2013-08-21 2014-08-13 超音波診断装置、超音波画像処理方法およびコンピュータ読み取り可能な非一時的な記録媒体
JP2014-164845 2014-08-13

Publications (1)

Publication Number Publication Date
CN104414685A true CN104414685A (zh) 2015-03-18

Family

ID=52480972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410415059.7A Pending CN104414685A (zh) 2013-08-21 2014-08-21 超声波诊断装置、超声波图像处理方法

Country Status (3)

Country Link
US (1) US9801614B2 (zh)
JP (1) JP2015061592A (zh)
CN (1) CN104414685A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106913357A (zh) * 2015-12-25 2017-07-04 通用电气公司 关节超声成像系统及其方法
CN109069119A (zh) * 2016-04-26 2018-12-21 皇家飞利浦有限公司 用于超声胎儿成像的3d图像合成
CN110772280A (zh) * 2018-07-31 2020-02-11 佳能医疗系统株式会社 超声波诊断装置和方法以及图像处理装置和方法
WO2021087687A1 (zh) * 2019-11-04 2021-05-14 深圳迈瑞生物医疗电子股份有限公司 超声图像分析方法、超声成像系统和计算机存储介质
CN114206224A (zh) * 2019-07-23 2022-03-18 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015131100A (ja) * 2013-12-12 2015-07-23 コニカミノルタ株式会社 超音波診断装置、超音波画像記録方法、及びコンピュータ読み取り可能な非一時的な記録媒体
CN106462967B (zh) * 2014-05-14 2020-06-30 皇家飞利浦有限公司 用于超声图像的基于模型的分割的采集取向相关特征
KR101579740B1 (ko) * 2014-09-01 2015-12-23 삼성메디슨 주식회사 초음파 진단장치, 그에 따른 초음파 진단 방법 및 그에 따른 컴퓨터 판독 가능한 저장매체
US20160135782A1 (en) * 2014-11-14 2016-05-19 General Electric Company Finger joint ultrasound imaging
JP7377739B2 (ja) * 2015-07-31 2023-11-10 株式会社ニデック 端末装置、端末制御プログラム、端末装置に操作される眼科装置、および検眼システム
JP6648587B2 (ja) * 2016-03-23 2020-02-14 コニカミノルタ株式会社 超音波診断装置
CN110622034B (zh) * 2017-05-11 2023-10-20 皇家飞利浦有限公司 超声诊断图像中的混响伪影消除
CN109223031A (zh) * 2017-07-11 2019-01-18 中慧医学成像有限公司 一种获得人体骨架的成像方法
WO2019154667A1 (en) * 2018-02-12 2019-08-15 Koninklijke Philips N.V. Workflow assistance for medical doppler ultrasound evaluation
JP6993907B2 (ja) * 2018-03-09 2022-01-14 富士フイルムヘルスケア株式会社 超音波撮像装置
US20210251602A1 (en) * 2018-08-22 2021-08-19 Koninklijke Philips N.V. System, device and method for constraining sensor tracking estimates in interventional acoustic imaging
CN113040823B (zh) * 2019-12-27 2023-03-21 深圳迈瑞生物医疗电子股份有限公司 一种超声成像设备及超声图像的分析方法
JP7326151B2 (ja) * 2019-12-27 2023-08-15 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用情報処理装置、及び医用情報処理プログラム
JP7457571B2 (ja) 2020-05-19 2024-03-28 富士フイルムヘルスケア株式会社 超音波診断装置及び診断支援方法
CN114259257A (zh) * 2020-09-16 2022-04-01 深圳迈瑞生物医疗电子股份有限公司 确定面积的方法、超声装置及计算机存储介质
WO2024030649A2 (en) * 2022-08-05 2024-02-08 Ochalski Pawel G Leaky joint syndrome diagnostics and treatments
KR20240065446A (ko) * 2022-10-28 2024-05-14 주식회사 엠티이지 초음파 영상 처리 방법 및 디바이스

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060239523A1 (en) * 2005-04-05 2006-10-26 Bradley University Radiographic imaging display apparatus and method
US20080187544A1 (en) * 2005-05-10 2008-08-07 Burkly Linda C Treating and evaluating inflammatory disorders
US20100113930A1 (en) * 2008-11-04 2010-05-06 Fujifilm Corporation Ultrasonic diagnostic device
US20100256504A1 (en) * 2007-09-25 2010-10-07 Perception Raisonnement Action En Medecine Methods and apparatus for assisting cartilage diagnostic and therapeutic procedures
CN102596052A (zh) * 2009-10-30 2012-07-18 株式会社日立医疗器械 超声波诊断装置、被检体的诊断对象部位的疾病评价用图像生成方法、及被检体的诊断对象部位的疾病评价用图像生成程序
CN102973293A (zh) * 2011-09-06 2013-03-20 通用电气公司 基于超声的病理自动化检测、量化和追踪的方法和系统
CN103202714A (zh) * 2012-01-16 2013-07-17 株式会社东芝 超声波诊断装置、医用图像处理装置以及医用图像处理方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5714328A (en) * 1980-06-30 1982-01-25 Tokyo Shibaura Electric Co Ultrasonic diagnosing device
JP2004057804A (ja) * 2002-06-05 2004-02-26 Fuji Photo Film Co Ltd 骨関節評価方法、装置およびそのためのプログラム
JP2004113644A (ja) * 2002-09-27 2004-04-15 Konica Minolta Holdings Inc 診断支援装置、診断支援方法、プログラム及び記憶媒体
US9498530B2 (en) * 2002-12-24 2016-11-22 Rinat Neuroscience Corp. Methods for treating osteoarthritis pain by administering a nerve growth factor antagonist and compositions containing the same
EP1852431A4 (en) * 2005-02-07 2009-11-04 Mitsubishi Tanabe Pharma Corp OPTICALLY ACTIVE TETRAHYDRONAPHTALENE DERIVATIVE
JP2007007434A (ja) * 2006-08-21 2007-01-18 Toshiba Corp 超音波診断装置
JP2010506147A (ja) * 2006-09-29 2010-02-25 エフ.ホフマン−ラ ロシュ アーゲー 関節リウマチを有する患者のための疾患進行のリスク評価
DK2308027T3 (da) * 2008-07-24 2013-06-24 Univ Gent Scoringssystem til monitorering af naturlig eller lægemiddelmodificeret sygdomsprogression i erosiv osteoarthritis i fingres interphalangealled
JP2011072522A (ja) * 2009-09-30 2011-04-14 Fujifilm Corp 超音波診断装置及び方法
JP2013051998A (ja) * 2011-08-31 2013-03-21 Toshiba Corp 超音波診断装置及び超音波診断装置の制御プログラム
US20140066469A1 (en) * 2012-08-21 2014-03-06 Department Of Veterans Affairs Treatment of diseases associated with inflammation

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060239523A1 (en) * 2005-04-05 2006-10-26 Bradley University Radiographic imaging display apparatus and method
US20080187544A1 (en) * 2005-05-10 2008-08-07 Burkly Linda C Treating and evaluating inflammatory disorders
US20100256504A1 (en) * 2007-09-25 2010-10-07 Perception Raisonnement Action En Medecine Methods and apparatus for assisting cartilage diagnostic and therapeutic procedures
US20100113930A1 (en) * 2008-11-04 2010-05-06 Fujifilm Corporation Ultrasonic diagnostic device
CN102596052A (zh) * 2009-10-30 2012-07-18 株式会社日立医疗器械 超声波诊断装置、被检体的诊断对象部位的疾病评价用图像生成方法、及被检体的诊断对象部位的疾病评价用图像生成程序
CN102973293A (zh) * 2011-09-06 2013-03-20 通用电气公司 基于超声的病理自动化检测、量化和追踪的方法和系统
CN103202714A (zh) * 2012-01-16 2013-07-17 株式会社东芝 超声波诊断装置、医用图像处理装置以及医用图像处理方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106913357A (zh) * 2015-12-25 2017-07-04 通用电气公司 关节超声成像系统及其方法
CN109069119A (zh) * 2016-04-26 2018-12-21 皇家飞利浦有限公司 用于超声胎儿成像的3d图像合成
CN110772280A (zh) * 2018-07-31 2020-02-11 佳能医疗系统株式会社 超声波诊断装置和方法以及图像处理装置和方法
CN110772280B (zh) * 2018-07-31 2023-05-23 佳能医疗系统株式会社 超声波诊断装置和方法以及图像处理装置和方法
CN114206224A (zh) * 2019-07-23 2022-03-18 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
CN114206224B (zh) * 2019-07-23 2023-09-08 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
WO2021087687A1 (zh) * 2019-11-04 2021-05-14 深圳迈瑞生物医疗电子股份有限公司 超声图像分析方法、超声成像系统和计算机存储介质

Also Published As

Publication number Publication date
JP2015061592A (ja) 2015-04-02
US9801614B2 (en) 2017-10-31
US20150057544A1 (en) 2015-02-26

Similar Documents

Publication Publication Date Title
CN104414685A (zh) 超声波诊断装置、超声波图像处理方法
Araki et al. PCA-based polling strategy in machine learning framework for coronary artery disease risk assessment in intravascular ultrasound: A link between carotid and coronary grayscale plaque morphology
US8840555B2 (en) System and method of ultrasound image processing
TWI750583B (zh) 醫學圖像分割方法、裝置、系統及圖像分割方法
KR102202617B1 (ko) 의료 영상 기반 복부 질환 분석 방법 및 장치
JP2017108769A (ja) 画像処理装置、画像処理方法、および、画像処理装置を備える超音波診断装置
US11826201B2 (en) Ultrasound lesion assessment and associated devices, systems, and methods
Rizi et al. Carotid wall longitudinal motion in ultrasound imaging: an expert consensus review
JP6648587B2 (ja) 超音波診断装置
Soleimani et al. Carotid artery wall motion estimation from consecutive ultrasonic images: Comparison between block-matching and maximum-gradient algorithms
US20190392579A1 (en) Method and system of performing medical treatment outcome assessment or medical condition diagnostic
US20190298304A1 (en) Medical diagnosis apparatus, medical image processing apparatus, and image processing method
US20240050062A1 (en) Analyzing apparatus and analyzing method
US11250564B2 (en) Methods and systems for automatic measurement of strains and strain-ratio calculation for sonoelastography
Zahnd et al. Dynamic block matching to assess the longitudinal component of the dense motion field of the carotid artery wall in B‐mode ultrasound sequences—Association with coronary artery disease
JP2016112285A (ja) 超音波診断装置
KR101059824B1 (ko) 초음파 영상을 이용한 경동맥 혈관의 내막두께와 중막두께의 비율 측정방법
US11672503B2 (en) Systems and methods for detecting tissue and shear waves within the tissue
WO2022103659A1 (en) System and method for detecting medical conditions
CN112515705A (zh) 用于投影轮廓启用的计算机辅助检测(cad)的方法和系统
JP2020081742A (ja) 診断装置及び診断方法
JP6286926B2 (ja) 超音波診断装置、画像処理方法、およびプログラム
US20230157661A1 (en) Ultrasound image analysis apparatus, ultrasound diagnostic apparatus, and control method for ultrasound image analysis apparatus
JP6931888B2 (ja) 解析装置及び解析プログラム
Cheung A machine learning approach for the objective sonographic assessment of patellar tendinopathy in collegiate basketball athletes

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
AD01 Patent right deemed abandoned

Effective date of abandoning: 20181102

AD01 Patent right deemed abandoned