CN106022369A - 一种基于超声背向散射信号成像的骨质评价方法 - Google Patents

一种基于超声背向散射信号成像的骨质评价方法 Download PDF

Info

Publication number
CN106022369A
CN106022369A CN201610325883.2A CN201610325883A CN106022369A CN 106022369 A CN106022369 A CN 106022369A CN 201610325883 A CN201610325883 A CN 201610325883A CN 106022369 A CN106022369 A CN 106022369A
Authority
CN
China
Prior art keywords
echo
bone
ultrasonic
signal
image
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
CN201610325883.2A
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.)
Nanjing University
Original Assignee
Nanjing University
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 Nanjing University filed Critical Nanjing University
Priority to CN201610325883.2A priority Critical patent/CN106022369A/zh
Publication of CN106022369A publication Critical patent/CN106022369A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/759Region-based matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images
    • G06V2201/033Recognition of patterns in medical or anatomical images of skeletal patterns

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Rheumatology (AREA)
  • Biophysics (AREA)
  • Evolutionary Biology (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明为一种基于超声背向散射信号成像的骨质评价方法,应用于医学超声领域。本方法使用单个超声探头对待测骨骼沿直线进行扫描,在每个扫描位置,对超声探头进行激励发射出超声波,超声波进入松质骨中受到其微结构的影响,形成散射波。其中,与入射方向相反的散射波被同一探头接收,即为背向散射信号。将每个扫描位置获得的回波信号映射为图像中的某一列,待扫描完成后即可获得回波图像。选取图像中背向散射部分的两个目标区域,并对其均值进行差分归一化运算,得到平均像素差分值来评价骨质情况。本方法相对于传统的基于透射法和轴向传输法的骨质评价方法,能更好地反映骨微结构的信息;相对于基于单点背向散射信号的骨质评价方法,能更好地反映区域骨质的整体情况。

Description

一种基于超声背向散射信号成像的骨质评价方法
技术领域
本发明涉及医学超声定量测量技术领域,具体涉及超声骨密度仪中一种基于背向散射信号成像的骨质评价方法。
背景技术
骨质疏松症是一种骨强度下降引起的疾病,会引起骨微结构的退化,导致骨折概率的增加。伴随着我国人口老龄化的加剧,骨质疏松症将会成为越来越严重的社会问题。目前,骨质疏松症最标准的诊断方法是使用双能X射线测量骨密度,但是由于其具有电离辐射而且仪器成本高,因此难以作为日常的检测手段。而超声诊断技术由于其无电离辐射、便捷、性价比高等特点,在骨密度检测方面得到了越来越多的应用。
目前常用的超声诊断方法是超声透射法和超声轴向传输法。相比而言,超声背向散射法具有以下优点:(1)超声透射法和超声轴向传输法需要多个探头或多个换能器,对于测量骨骼的位置有很大的限制,而超声背向散射法使用同一个探头对超声波进行收发,测量位置更加灵活。(2)骨质疏松不仅仅带来骨矿密度的下降,也和骨骼微结构的退化相关。超声透射法和超声轴向传输法只能够检测骨矿密度的变化,而超声背向散射信号中含有骨微结构的信息,因此能够更加准确地诊断骨质疏松症。
骨质分为皮质骨和松质骨,松质骨由许多骨小梁交织而成,当超声波进入松质骨内部时,由于波长接近骨小梁的厚度,会产生各个方向上的散射波,其中被发射探头接收的散射波即为背向散射信号。由于松质骨的各向异性,松质骨中各点的特性并不一致,因此通过单点的背向散射信号来分析骨质情况难以获得整体骨质情况,因此对区域内各点的背向散射信号成像,再对成像结果进行分析,就能够掌握区域内整体骨质的情况,这使得测量结果更有参考意义。
发明内容
本发明要解决的问题是:超声透射法和超声轴向传输法这两种传统的超声测量骨密度方法具有测量位置限制和无法获得骨微结构信息两个缺点,而基于单点超声背向散射信号的测量方法难以回避松质骨各向异性的特点导致的结果不精确。本发明提供的方法在弥补了传统超声方法缺点的基础之上,将超声背向散射信号成像来对区域内骨质情况进行综合评价,提高了骨质测量的精确性和稳定性。
本发明的技术方案由以下步骤实现:
(1)使用超声探头沿直线扫描待测骨,在每个扫描位置,对超声探头进行激励,向待测骨头发射超声波。超声波传播到松质骨中受到其微结构的影响形成散射波,其中的背向散射波被同一探头接收;
(2)在每个扫描位置上,将接收到的回波信号映射到图像中的一列,待扫描完成后获得扫描区域的回波图像;
(3)选取图像中背向散射部分的两个目标区域的均值进行差分归一化运算得到平均像素差分值,根据结果对区域内骨质情况进行评价。
所述的测量系统包括探头、前端采集电路和PC机处理软件。测量系统使用的探头采用单阵元非聚焦超声换能器,探头频率为2.25MHz,测量部位为根骨。前端采集电路的作用是产生探头激励脉冲、回波处理、模数转换和上传数据到PC机。其中回波处理部分包括可变增益放大器,确保收到的回波信号经过处理后进入模数转换时处于合适的量程范围内。模数转换器采用ANALOG公司的AD6645,采样位数为14bit,最大采样速率为105MSPS。前端采集电路使用USB总线和PC机通信,PC机软件读取USB总线上传的数据并对其进行解析,绘制成图像,并将从图像中提取的参数和数据库中的标准值进行对比,随后把测量结果反馈给测试者。
所述步骤(1)中,先将探头置于根骨部位,使用医用超声耦合剂进行耦合,在测量区域内移动探头,在移动过程中,探头不间断发射超声波并接收回波信号,前端采集电路完成回波信号的预处理之后进行模数转换,并将模数转换的结果和对应可变增益放大器的增益值通过USB总线上传到PC机进行下一步的数据处理。
所述步骤(2)中,在每个扫描位置上,将接收到的回波信号映射到图像中的一列,这样在扫描完成后所有列拼接在一起形成了完整的二维回波图像。图像的横轴代表了每个扫描位置,图像的纵轴代表回波信号的时间轴,图像的每一个像素值的大小代表了每个扫描位置的回波中信号点的幅度。图像采用256阶灰度图表示,每个像素的取值范围在0到255之间。在每个扫描位置上获得的回波信号先根据增益值换算到回波信号的绝对值,再将其整体线性映射到0到255之间,填充到图像中新的一列上,如此往复即可获得整个扫描区域内的二维回波图像。
所述步骤(3)中,由于受到表面反射波和皮肤等中间组织的影响,回波的初始部分并非所需要的背向散射信号,因此需要在二维回波图像中选取背向散射信号部分的两个目标区域,两个目标区域均为矩形,长宽相同,面积相同。目标区域的长度和回波图像的横向长度相同,目标区域的宽度为时间Tw,并且两个目标区域间隔时间为Ts。选取到合适的目标区域后,将两个目标区域内的像素均值相减,再除以归一化系数Tw+Ts,得到最终的平均像素差分值。由于骨密度的不同,两个目标区域内的回波信号的衰减程度也不一样,当骨密度越大时,两个目标区域内回波信号的衰减程度相差也越大,因此两个目标区域的平均像素差分值与骨密度相关。最后将结果与数据库中的标准值进行比对,获得区域内的骨质情况信息。
本发明所述的一种基于超声背向散射信号成像的骨质评价方法应用于骨密度的超声检测中,既避免了传统超声透射法和超声轴向传输法带来的测量位置限制和测量结果误差,又避免了松质骨各向异性带来的单点测量的误差,实现了对区域内骨质情况的整体评价。
附图说明
图1是本发明测量系统整体结构框图。
图2是本发明采集信号示意图。
图3是本发明成像结果示意图。
具体实施方式
本发明使用的测量系统整体结构框图如图1所示,测量系统采用同一个探头进行超声波的发射和回波的接收,测量系统由FPGA来进行控制。由FPGA产生的激励信号经过高压模块放大,连接到探头对其进行激励产生超声波。回波信号被探头接收,转换成相应的电信号传输到预处理模块。预处理模块完成与高压电路的隔离和回波信号的前置放大,放大后的回波信号进入可控放大模块,由FPGA根据其峰峰值控制增益的大小,防止信号幅度太大超过模数转换器的量程或者太小影响模数转换结果的精度,以保证信号处于合适的模数转换的量程范围内。模数转换器采用ANALOG公司的14位单芯片模数转换器AD6645,最大采样速率可以达到105MSPS。其后回波信号经过模数转换器由模拟信号转换成数字信号,数字信号和对应增益值由FPGA通过USB总线上传送给PC机。
测量开始时,将超声探头置于测试者根骨部位,沿着根骨进行直线扫描,在每个扫描位置,由测量系统不断激励超声探头,超声探头产生的超声波传播到松质骨中,受到其微结构的影响产生散射波,其中与发射超声波方向相反的背向散射波被超声探头所接收,与表面反射波等信号一起构成了超声探头接收到的回波信号。回波信号如图2所示,前4us为表面反射波,其后的波形为背向散射波。
测量系统将探头接收到的回波转换成数字信号经过USB总线传输到PC机,在PC机上读取USB总线上的数据进行处理。在PC机上每读取到的一帧信号代表探头在一个扫描位置接收到的回波,将其映射到图像的一列上,其中的像素值与信号幅度相对应。图像用256阶灰度图表示,每个像素取值范围在0到255之间。在每个扫描位置上获得的一帧回波信号先将其根据增益值换算到回波信号的绝对值上,再将回波信号绝对值线性映射到0到255的区间上。假设某帧回波信号绝对值中幅度最大值为Max,最小值为Min,对于某一个幅度为A的点映射后的像素值P为:
P = A - M i n M a x - M i n * 255
在探头沿着扫描区域移动的过程中,不断发射超声波以及接收回波信号,不断有新的列添加到图像中,图像的长度也在不断扩展。待扫描完成后,便可形成回波信号图像,图像的横轴代表了每个扫描位置,图像的纵轴代表了回波信号的时间轴,图像的像素值代表了回波信号的幅度。
回波图像如图3所示,最上层的部分为表面反射波的图像,需要排除在分析之外。其余的部分为分析需要的背向散射信号图像部分。选取背向散射信号图像的两个矩形目标区域进行分析,如图3中的两个白色边框矩形所示,两个目标区域的长度和宽度相同,长度取整个图像的横轴范围,宽度取时间长度Tw,两个目标区域间隔时间长度Ts,并且令Tw=1us,Ts=1us。选取到合适的目标区域后,分别计算出两个目标区域内像素值的均值为M1和M2,则最终测量结果Val为:
V a l = M 1 - M 2 T w + T s
Val即为获得评价骨质情况的平均像素差分值。当骨密度越大,两个目标区域内的信号衰减程度将会越明显,因此两个目标区域内像素均值差距将会越大,得到的最终结果也会越大。将得到平均像素差分值和数据库中的标准值比对即可获得骨质情况的评价结果。

Claims (6)

1.一种基于超声背向散射信号成像的骨质评价方法,其特征在于具体步骤如下:
(1)使用超声探头沿直线扫描待测骨,在每个扫描位置,对超声探头进行激励,向待测骨头发射超声波。超声波传播到松质骨中受到其微结构的影响形成散射波,其中的背向散射波被同一探头接收;
(2)在每个扫描位置上,将接收到的回波信号映射到图像中的一列,待扫描完成后获得扫描区域的回波图像;
(3)选取图像中背向散射部分的两个目标区域,对两个目标区域的均值进行差分归一化运算得到平均像素差分值,根据平均像素差分值对区域内骨质情况进行评价。
2.根据权利要求1所述的方法,其特征是超声探头采用单阵元非聚焦超声换能器,测量部位为根骨。
3.根据权利要求1所述的方法,其特征是采用的测量系统包括探头、前端采集电路和PC机处理程序。前端采集电路的作用是产生探头激励脉冲、回波处理、模数转换和上传数据到PC机。回波处理部分包括可变增益放大器,确保收到回波信号经过处理后进入模数转换时处于合适的量程范围内。
4.根据权利要求1所述的方法,其特征是在步骤(2)中,得到的回波信号根据其增益值换算成回波信号的绝对值,再通过绝对值信号的最大值和最小值线性映射到256阶灰度图中。
5.根据权利要求1所述的方法,其特征是在步骤(3)中,回波图像的长度代表扫描距离,宽度代表接收到回波的时间。在回波图像中除去表面反射波的部分,在剩下的背向散射信号图像中选取宽度为Tw并间隔Ts的两个目标区域进行运算。
6.根据权利要求1所述的方法,其特征是在步骤(3)中,对两个目标区域内像素值的均值相减,再除以归一化系数Tw+Ts,得到平均像素差分值。若骨密度越大,回波信号随时间的衰减将会越明显,两个目标区域内像素均值差异将会越大,得到的最终结果也会越大。将结果与数据库中的标准值比对,即可得到骨质情况的评价信息。
CN201610325883.2A 2016-05-12 2016-05-12 一种基于超声背向散射信号成像的骨质评价方法 Pending CN106022369A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610325883.2A CN106022369A (zh) 2016-05-12 2016-05-12 一种基于超声背向散射信号成像的骨质评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610325883.2A CN106022369A (zh) 2016-05-12 2016-05-12 一种基于超声背向散射信号成像的骨质评价方法

Publications (1)

Publication Number Publication Date
CN106022369A true CN106022369A (zh) 2016-10-12

Family

ID=57097557

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610325883.2A Pending CN106022369A (zh) 2016-05-12 2016-05-12 一种基于超声背向散射信号成像的骨质评价方法

Country Status (1)

Country Link
CN (1) CN106022369A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107854141A (zh) * 2017-11-24 2018-03-30 南京大学 一种基于超声背向散射信号的骨小梁间距测量系统
CN110769754A (zh) * 2017-06-21 2020-02-07 夏里特柏林大学医学院 用于测定皮质骨的系统、方法和计算机程序产品
CN114287967A (zh) * 2021-11-24 2022-04-08 中国科学院深圳先进技术研究院 一种基于单阵元超声脉冲回波的血压测量系统及方法
WO2023092329A1 (zh) * 2021-11-24 2023-06-01 中国科学院深圳先进技术研究院 一种基于单阵元超声脉冲回波的血压测量系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101658434A (zh) * 2009-09-10 2010-03-03 复旦大学 用于定征松质骨微结构的超声频谱偏移量参数成像方法
CN102198009A (zh) * 2011-06-14 2011-09-28 复旦大学 基于超声背散射信号参量的松质骨诊断系统
CN104490418A (zh) * 2014-09-25 2015-04-08 深圳市恩普电子技术有限公司 一种基于信号统计分析的超声图像自动优化方法
CN105433945A (zh) * 2015-06-30 2016-03-30 北京大学 一种骨密度检测设备和检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101658434A (zh) * 2009-09-10 2010-03-03 复旦大学 用于定征松质骨微结构的超声频谱偏移量参数成像方法
CN102198009A (zh) * 2011-06-14 2011-09-28 复旦大学 基于超声背散射信号参量的松质骨诊断系统
CN104490418A (zh) * 2014-09-25 2015-04-08 深圳市恩普电子技术有限公司 一种基于信号统计分析的超声图像自动优化方法
CN105433945A (zh) * 2015-06-30 2016-03-30 北京大学 一种骨密度检测设备和检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐峰 等: "超声背散射法评价松质骨的Barker码方法研究", 《声学技术》 *
黄凯 等: "松质骨超声背散射信号频移参数成像", 《声学技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110769754A (zh) * 2017-06-21 2020-02-07 夏里特柏林大学医学院 用于测定皮质骨的系统、方法和计算机程序产品
CN110769754B (zh) * 2017-06-21 2023-06-27 夏里特柏林大学医学院 用于测定皮质骨的系统、方法和计算机程序产品
CN107854141A (zh) * 2017-11-24 2018-03-30 南京大学 一种基于超声背向散射信号的骨小梁间距测量系统
CN114287967A (zh) * 2021-11-24 2022-04-08 中国科学院深圳先进技术研究院 一种基于单阵元超声脉冲回波的血压测量系统及方法
WO2023092329A1 (zh) * 2021-11-24 2023-06-01 中国科学院深圳先进技术研究院 一种基于单阵元超声脉冲回波的血压测量系统及方法
CN114287967B (zh) * 2021-11-24 2023-12-29 中国科学院深圳先进技术研究院 一种基于单阵元超声脉冲回波的血压测量系统及方法

Similar Documents

Publication Publication Date Title
CN104605891B (zh) 检测剪切波在生物组织中传播速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法
US7535797B2 (en) High-resolution ultrasound displacement measurement apparatus and method
EP0066343B1 (en) Method and apparatus for measuring ultrasonic attenuation characteristics
US9492139B2 (en) Non-imaging low frequency ultrasonic testing and diagnostic evaluation system
CN103269639B (zh) 利用质心估计剪切波速度
JP4928989B2 (ja) 超音波診断装置および超音波診断装置の制御プログラム
EP0123427B1 (en) Ultrasonic medium characterization
CN106022369A (zh) 一种基于超声背向散射信号成像的骨质评价方法
CN102123668A (zh) 使用未聚焦发送波束的高帧率定量多普勒流成像
DE102017202141B4 (de) Schallgeschwindigkeitsbildgebung unter Verwendung von Scherwellen
CN102640014A (zh) 图像生成装置、图像生成方法和程序
CN101923072A (zh) 线性超声相控阵成像方法
FR3072870A1 (fr) Estimation viscoelastique d'un tissu a partir d'une vitesse de cisaillement dans une imagerie medicale a ultrasons
JP5905080B2 (ja) オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成
CN109900805A (zh) 基于频域稀疏反演的tofd盲区内缺陷定量检测方法
WO2019083491A1 (en) METHOD AND APPARATUS FOR ULTRASONIC MEASUREMENT AND IMAGING OF THE ELASTICITY OF REAL-TIME BIOLOGICAL TISSUES
CN101658434A (zh) 用于定征松质骨微结构的超声频谱偏移量参数成像方法
WO2024213075A1 (zh) 一种超声波成像设备和超声波成像系统
CN108508093A (zh) 一种工件缺陷高度的检测方法及系统
CN111528921A (zh) 一种信息采集装置
JPH08511184A (ja) 超音波走査のための方法及び装置
Bigelow et al. Scatterer size estimation in pulse-echo ultrasound using focused sources: Theoretical approximations and simulation analysis
CN104644215A (zh) 一种多模三维超声探测与成像装置及方法
CN102551793A (zh) 超声高帧率组织多普勒成像方法及装置
EP3179269A1 (en) Analysis methods of ultrasound echo signals based on statistics of scatterer distributions

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20161012