TWI703962B - 介質黏彈性定量方法及裝置 - Google Patents
介質黏彈性定量方法及裝置 Download PDFInfo
- Publication number
- TWI703962B TWI703962B TW107121598A TW107121598A TWI703962B TW I703962 B TWI703962 B TW I703962B TW 107121598 A TW107121598 A TW 107121598A TW 107121598 A TW107121598 A TW 107121598A TW I703962 B TWI703962 B TW I703962B
- Authority
- TW
- Taiwan
- Prior art keywords
- angle
- slope
- medium
- position time
- time map
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N11/10—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material
- G01N11/16—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material by measuring damping effect upon oscillatory body
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0092—Visco-elasticity, solidification, curing, cross-linking degree, vulcanisation or strength properties of semi-solid materials
- G01N2203/0094—Visco-elasticity
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Vascular Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本發明揭露了介質黏彈性定量方法及裝置。在介質受到振動激發後,得到振動傳播的位置時間圖。利用角度投影確定位置時間圖中訊號能量最大的角度,該訊號能量最大的角度對應位置時間圖的斜率,而位置時間圖的斜率就是振動在介質中的傳播速度。由於振動在介質中的傳播速度和介質的黏彈性相關,所以在得到位置時間圖的斜率之後,就能定量計算出介質的黏彈性參數。本發明無需從位置時間圖中選取特徵點來計算位置時間圖的斜率,不會受到雜訊干擾且計算量小,可以高效準確的進行介質黏彈性定量。
Description
本發明涉及測量技術領域,特別涉及介質黏彈性定量方法及裝置。
對介質進行振動激發時,振動在介質中的傳播特性與介質的黏彈性有關,通過測量振動的傳播特性,可以對介質的黏彈性進行定量。
上述原理目前已被應用至多個技術領域,以醫學檢測為例,檢測肝臟、甲狀腺和肌肉等器官或組織時,通過對介質的黏彈性進行定量,可以對病變進行定位。
因此,如何進行高效準確的介質黏彈性定量,是一個需要解決的問題。
本發明實施例提供了介質黏彈性定量方法及裝置。為了對披露的實施例的一些方面有一個基本的理解,下面給出了簡單的概括。該概括部分不是泛泛評述,也不是要確定關鍵/重要組成元素或描繪這些實施例的保護範圍。其唯一目的是用簡單的形式呈現一些概念,以此作為後面的詳細說明的序言。
根據本發明實施例的第一方面,提供了一種介質黏彈性定量方法,該方法包括: 在介質受到振動激發後,得到振動傳播的位置時間圖; 沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率; 根據該斜率,得到該介質的黏彈性參數。
在該方法的基礎上,作為可選的實施例一,該沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括: 對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 將積分值最大的角度,確定為該位置時間圖的斜率線的斜率角; 利用該斜率角確定該斜率線的斜率。
在該方法的基礎上,作為可選的實施例二,該沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括: 針對該位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣; 得到該各個角度的影像紋理特徵; 利用該影像紋理特徵,將訊號能量最大的角度,確定為該位置時間圖的斜率線的斜率角; 利用該斜率角確定該斜率線的斜率。
在該方法、該實施例一或該實施例二的基礎上,作為可選的實施例三,該方法還包括: 在該角度投影之前,濾除該位置時間圖中的反射波。
在該實施例三的基礎上,作為可選的實施例四,該濾除該位置時間圖中的反射波,包括:對該位置時間圖進行方向濾波。
在該方法、該實施例一或該實施例二的基礎上,作為可選的實施例五,該得到振動傳播的位置時間圖,包括: 沿設定的振動傳播方向,得到該振動傳播的位置時間圖。
根據本發明實施例的第二方面,提供了一種介質黏彈性定量裝置,該裝置包括: 影像模組,用於在介質受到振動激發後,得到振動傳播的位置時間圖; 確定模組,用於沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率; 定量模組,用於根據該斜率,得到該介質的黏彈性參數。
在該裝置的基礎上,作為可選的實施例一,該確定模組,包括: 計算子模組,用於對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 確定子模組,用於將該積分子模組計算出的積分最大值的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的斜率。
在該裝置的基礎上,作為可選的實施例二,該確定模組,包括: 計算子模組,用於針對該位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣; 確定子模組,用於得到該各個角度的影像紋理特徵;利用該影像紋理特徵,將訊號能量最大的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的斜率。
在該裝置、該實施例一或該實施例二的基礎上,作為可選的實施例三,該裝置還包括: 濾波模組,用於在該角度投影之前,濾除該位置時間圖中的反射波。
在該裝置、該實施例一或該實施例二的基礎上,作為可選的實施例四,該影像模組沿設定的振動傳播方向,得到該振動傳播的位置時間圖。
根據本發明實施例的第三方面,提供了一種介質黏彈性定量裝置,該裝置包括: 記憶體,儲存執行指令; 處理器,被配置於讀取該執行指令,完成如下操作: 在介質受到振動激發後,得到振動傳播的位置時間圖; 沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率; 根據該斜率,得到該介質的黏彈性參數。
本發明實施例提供的技術方案可以包括以下有益效果: 利用角度投影確定位置時間圖中訊號能量最大的角度,該訊號能量最大的角度對應位置時間圖的斜率,而位置時間圖的斜率就是振動在介質中的傳播速度。由於振動在介質中的傳播速度和介質的黏彈性相關,所以在得到位置時間圖的斜率之後,就能定量計算出介質的黏彈性參數。本發明實施例無需從位置時間圖中選取特徵點來計算位置時間圖的斜率,不會受到雜訊干擾且計算量小,可以高效準確的進行介質黏彈性定量。
應當理解的是,以上的一般描述和後文的細節描述僅是示例性和解釋性的,並不能限制本發明。
以下描述和附圖充分地示出本發明的具體實施方案,以使本領域的技術人員能夠實踐它們。實施例僅代表可能的變化。除非明確要求,否則單獨的部件和功能是可選的,並且操作的順序可以變化。一些實施方案的部分和特徵可以被包括在或替換其他實施方案的部分和特徵。本發明的實施方案的範圍包括申請專利範圍的整個範圍,以及申請專利範圍的所有可獲得的等同物。在本文中,各實施方案可以被單獨地或總地用術語“發明”來表示,這僅僅是為了方便,並且如果事實上揭露了超過一個的發明,不是要自動地限制該應用的範圍為任何單個發明或發明構思。本文中,諸如第一和第二等之類的關係術語僅僅用於將一個實體或者操作與另一個實體或操作區分開來,而不要求或者暗示這些實體或操作之間存在任何實際的關係或者順序。而且,術語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含,從而使得包括一系列要素的程序、方法或者裝置不僅包括那些要素,而且還包括沒有明確列出的其他要素。本文中各個實施例採用遞進的方式描述,每個實施例重點說明的都是與其他實施例的不同之處,各個實施例之間相同相似部分互相參見即可。對於實施例揭露的結構、產品等而言,由於其與實施例揭露的部分相對應,所以描述的比較簡單,相關之處參見方法部分說明即可。
第1圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖。如第1圖所示,該方法包括如下步驟。
在步驟11中,在介質受到振動激發後,得到振動傳播的位置時間圖。
在步驟12中,沿位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度,該訊號能量最大的角度對應位置時間圖的斜率。
預設角度範圍,指根據實際情況,選取的進行角度投影的角度範圍。作為一種可選的實施方式,預設角度範圍可以是360度,此時需進行全角度的角度投影。作為另一種可選的實施方式,根據得到的位置時間圖的特性,來選擇進行角度投影的角度範圍。步驟11中得到的位置時間圖的橫軸指示時間、縱軸指示位置,如果振動只由振動激發起點處向遠處傳播,則振動傳播的速度無限大時,在位置時間圖上接近於與縱軸平行的直線,而振動傳播的速度無限小時,在位置時間圖上接近于與橫軸平行的直線,此時預設角度範圍是90度即可滿足需求,而無需做全角度投影,進而提高介質黏彈性定量的效率。如果振動除了由振動激發起點處向遠處傳播,還可沿相反方向繼續傳播,則預設角度範圍可以是180度。至於預設角度範圍的實際起始點和終止點,在直角坐標系保持不變的情況下,與0度的起始點及逆時針或順時針的旋轉方向相關,可以按需要選取,只要保證預設角度範圍即可。
各個角度,指在預設角度範圍內進行角度投影的各個角度。具體角度的選取根據時間精度需求和計算速度需求來確定,時間精度要求越高時,角度選取的精度要求越高,計算速度需求越高時,角度選取的精度要求越低。例如,可以從0.0.1度到1度的範圍內選擇。
角度投影,指標對設定角度進行影像特徵識別或提取,以確定訊號能量最大的角度。
在步驟13中,根據斜率,得到介質的黏彈性參數。
黏彈性參數包括黏性參數和彈性參數中的至少一項。
位置時間圖的斜率由單位時間的振動傳播的距離決定,即介質中振動傳播的速度。在均勻的介質中,振動傳播的速度和介質的黏彈性有關,在得到位置時間圖的斜率之後,就可以定量計算出介質的黏彈性參數,因此如何高效準確的獲得上述斜率成為介質黏彈性定量的關鍵。本示例性實施例利用角度投影確定訊號能量最大的角度,由於該訊號能量最大的角度對應位置時間圖的斜率,即相當於獲得了位置時間圖的斜率。這種方法無需從位置時間圖中選取波峰、波谷或振動的某一相位作為特徵點來計算位置時間圖的斜率,不會受到雜訊干擾且計算量小,是一種高效準確的進行介質黏彈性定量的方法。
在一示例性的實施例中,通過機械振動、聲輻射力或其他可以產生振動的方式,對介質進行振動激發後,介質產生振動,振動在介質中傳播。由於上述振動在介質中的傳播速度有限,因此可以利用檢測波對介質進行動態成像。上述檢測波可以是光波、超音波等。上述動態成像可以是一維成像、二維成像或三維成像等。
由於上述振動在介質中傳播時,在不同的時刻,波前會沿著傳播方向到達不同的位置。檢測波對介質成像產生的回波訊號會發生相位去相關,利用這種相位去相關的特性,可以通過互相關、自相關、光流等演算法得到介質的運動資訊,沿設定的振動傳播方向,得到位置時間圖。上述演算法可以是各種基於塊匹配的方法或非塊匹配的方法。上述設定的振動傳播方向,在振動只在一個傳播方向上傳播時,為振動的實際傳播方向,在振動在多個傳播方向上傳播時,為選擇出的某一個傳播方向。例如,當介質為一均勻薄片時,對介質進行振動激發後,振動將沿該薄片的延展方向進行傳播,此時設定的振動傳播方向就是該振動的實際傳播方向。又例如,當介質為立體不規則形狀時,振動傳播的波前為立體形狀,例如振動傳播的波前為橢球,則沿不同的振動傳播方向得到的位置時間圖是不同的,此時設定的振動傳播方向就是選擇出的某一個感興趣的傳播方向。上述感興趣的傳播方向根據實際要測量的方向來確定,例如可以是振動傳播最快的方向、振動傳播最慢的方向和振動傳播速度為某一區間的方向中的至少一個。
在一示例性的實施例中,由於振動在介質中傳播時,遇介質的邊緣或異物時會產生反射波,為提高後續處理的精度,如第2圖所示,在進行角度投影之前,還可以包括步驟11’,即對位置時間圖中的反射波進行濾除。濾除的方式可以有很多種,方向濾波是其中一種實現方式。
在一示例性的實施例中,通過角度投影來確定訊號能量最大的角度,進而得到位置時間圖的斜率,可以通過積分計算來實現。例如,對位置時間圖進行沿預設角度範圍內的各個角度的積分計算,當積分角度和振動傳播方向一致時能量聚集,此時得到的積分值最大,因此將積分值最大的角度為位置時間圖的斜率線的斜率角。根據得到的斜率角,再結合位置和時間資訊,即可以得到位置時間圖的斜率線的斜率。上述積分計算又稱為Radon變換。
在一示例性的實施例中,由於通過計算灰階共生矩陣可以獲得影像紋理特徵,影像紋理特徵又可以體現訊號能量的大小,因此可以用灰階共生矩陣來獲取訊號能量最大的角度的資訊。基於上述原理,通過角度投影來確定訊號能量最大的角度,進而得到位置時間圖的斜率,可以通過灰階共生矩陣的計算來實現。例如,針對位置時間圖,先沿預設角度範圍內的各個角度計算灰階共生矩陣。然後,利用灰階共生矩陣,得到各個角度的影像紋理特徵。接著,利用影像紋理特徵,將訊號能量最大的角度,確定為位置時間圖的斜率線的斜率角。最後,利用斜率角確定斜率線的斜率。
根據力學原理,介質的黏彈性決定了振動在其中的傳播速度,因此通過得到位置時間圖的斜率,可以獲知振動在介質中傳播的速度,進而根據力學原理,可以定量得出介質的黏彈性參數。這裡的黏彈性參數可以包括剪切模量、楊氏模量、剪切黏彈性、剪切黏度、機械阻抗、機械鬆弛時間、各向異性等。
下面以一種具體的應用場景給出本發明實施例中介質黏彈性定量方法的應用。
在對人體肝臟等黏彈性介質進行無損黏彈性檢測時,需要對介質黏彈性進行定量。檢測裝置中包括激發裝置和成像裝置,其中激發裝置對待檢測介質進行振動激發,成像裝置利用超音波對振動激發後的介質進行成像。振動在介質中傳播時,在不同的時刻,波前會沿著傳播方向到達不同的位置,形成位置時間圖。上述波前可以是波峰、波谷、或振動的同一相位中的一種。
如第3圖所示,在這種具體應用場景的介質黏彈性定量方法可以包括如下步驟。
在步驟31中,對介質進行振動激發。
在步驟32中,利用超音波對介質進行動態成像。
在步驟33中,從介質的成像中得到振動傳播的位置時間圖。
在步驟34中,對位置時間圖進行方向濾波。
在步驟35中,進行Radon變換,確定出位置時間圖的斜率。
在步驟36中,根據確定出的斜率和力學原理,計算得到介質的黏彈性參數。
在以上介質黏彈性定量方法的各個示例性實施例中,當設定的振動傳播方向為至少兩個時,每個設定的振動傳播方向會對應得到一個位置時間圖,進而會得到該位置時間圖對應的介質的黏彈性參數。綜合得到的至少兩套黏彈性參數,可以更全面的評價介質的黏彈性。
上述給出的介質黏彈性定量方法的各個示例性實施例,可以根據情況進行組合,這裡並不限定各個示例性實施例之間的組合關係。
第4圖是一示例性實施例示出的介質黏彈性定量裝置的框圖,該裝置可以位於介質黏彈性檢測裝置的控制主機中,例如在醫療檢測領域,可以位於肝臟無損檢測裝置的控制主機中。該裝置還可以位於雲端,介質黏彈性檢測裝置的檢測資料需要在雲端進行處理。
第4圖所示的裝置包括:影像模組41、確定模組42和定量模組43。
影像模組41,用於在介質受到振動激發後,得到振動傳播的位置時間圖。
確定模組42,用於沿位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度。上述訊號能量最大的角度對應位置時間圖的斜率。
定量模組43,用於根據位置時間圖的斜率,得到介質的黏彈性參數。
在一示例性實施例中,如第5圖所示,確定模組42採用Radon變換來進行角度投影並確定訊號能量最大的角度。此時,確定模組42包括:計算子模組421和確定子模組422。
計算子模組421,用於對位置時間圖進行沿預設角度範圍內的各個角度的積分計算。
確定子模組422,用於將積分子模組421計算出的積分最大值的角度,確定為位置時間圖的斜率線的斜率角;通過該斜率角確定位置時間圖的斜率線的斜率。
作為可選的實施方式,在使用灰階共生矩陣來確定斜率角時,計算子模組421可以用於針對位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣。確定子模組422可以用於得到該各個角度的影像紋理特徵;利用該影像紋理特徵,將訊號能量最大的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的斜率。
在一示例性實施例中,如第6圖所示,介質黏彈性定量裝置中還包括:濾波模組44,用於在確定模組42進行角度投影之前,濾除位置時間圖中的反射波。
在一示例性實施例中,影像模組41沿設定的振動傳播方向,得到振動傳播的位置時間圖。
下面以一種具體的應用場景給出本發明實施例中介質黏彈性定量裝置的應用。
在對人體肝臟等黏彈性介質進行無損黏彈性檢測時,需要對介質黏彈性進行定量。檢測裝置中包括激發裝置和成像裝置,其中激發裝置對待檢測介質進行振動激發,成像裝置利用超音波對振動激發後的介質進行成像。
振動在介質中傳播時,在不同的時刻,波前會沿著傳播方向到達不同的位置,形成位置時間圖。上述波前可以是波峰、波谷、或振動的同一相位中的一種。介質黏彈性定量裝置沿設定傳播方向,從介質的成像中得到振動傳播的位置時間圖。然後介質黏彈性定量裝置對位置時間圖沿各個角度進行積分計算,將積分值最大的角度確定為位置時間圖的斜率線的斜率角,進而確定出位置時間圖的斜率。最後介質黏彈性定量裝置根據確定出的斜率和力學原理,計算得到介質的黏彈性參數。
第7圖是根據一示例性實施例示出的介質黏彈性定量裝置的框圖,該裝置可以位於介質黏彈性檢測裝置的控制主機中,例如在醫療檢測領域,可以位於肝臟無損檢測裝置的控制主機中。該裝置還可以位於雲端,介質黏彈性檢測裝置的檢測資料需要在雲端進行處理。
上述給出的介質黏彈性定量裝置的各個示例性實施例,可以根據情況進行組合,這裡並不限定各個示例性實施例之間的組合關係。
第7圖所示的裝置包括:記憶體71和處理器72。
記憶體71中儲存有執行指令。
處理器72,被配置於讀取記憶體71中的執行指令,執行前文該介質黏彈性定量方法的各示例性實施例中的部分或全部步驟。處理器72可以由晶片實現。
如果第7圖所示的介質黏彈性定量裝置位於介質黏彈性檢測裝置的控制主機中,可以通過匯流排、無線等方式與介質黏彈性定量裝置中的激發裝置、成像裝置進行連接,此時該裝置上具備與實現上述連接的介面及相應通訊機制。
如果第7圖所示的介質黏彈性定量裝置位於雲端,可以通過網路與介質黏彈性檢測裝置進行通訊。
應當理解的是,本發明並不侷限於上面已經描述並在附圖中示出的流程及結構,並且可以在不脫離其範圍進行各種修改和改變。本發明的範圍僅由所附的申請專利範圍來限制。
11、11’、12、13、31、32、33、34、35、36‧‧‧步驟41‧‧‧影像模組42‧‧‧確定模組43‧‧‧定量模組44‧‧‧濾波模組71‧‧‧記憶體72‧‧‧處理器421‧‧‧計算子模組422‧‧‧確定子模組
此處的附圖被併入說明書中並構成本說明書的一部分,示出了符合本發明的實施例,並與說明書一起用於解釋本發明的原理。 第1圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第2圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第3圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第4圖是根據一示例性實施例示出的介質黏彈性定量裝置的框圖; 第5圖是第4圖中所示的確定模組的框圖; 第6圖是根據一示例性實施例示出的介質黏彈性定量裝置的框圖; 第7圖是根據一示例性實施例示出的介質黏彈性定量裝置的框圖。
11、12、13‧‧‧步驟
Claims (12)
- 一種介質黏彈性定量方法,其特徵在於,該方法包括: 在一介質受到振動激發後,得到一振動傳播的一位置時間圖; 沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的一斜率; 根據該斜率,得到該介質的黏彈性參數。
- 如申請專利範圍第1項所述的方法,其中,該沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括: 對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 將積分值最大的角度,確定為該位置時間圖的一斜率線的一斜率角; 利用該斜率角確定該斜率線的斜率。
- 如申請專利範圍第1項所述的方法,其中,該沿該位置時間圖上的預設角度範圍內的各個角度進行一角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括: 針對該位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣; 得到該各個角度的一影像紋理特徵; 利用該影像紋理特徵,將訊號能量最大角度,確定為該位置時間圖的斜率線的斜率角; 利用該斜率角確定該斜率線的斜率。
- 如申請專利範圍第1項至第3項中任一項所述的方法,其中,該方法還包括: 在該進行角度投影之前,濾除該位置時間圖中的反射波。
- 如申請專利範圍第4項所述的方法,其中,該濾除該位置時間圖中的反射波,包括: 對該位置時間圖進行方向濾波。
- 如申請專利範圍第1項至第3項中任一項所述的方法,其中,該得到振動傳播的位置時間圖,包括: 沿設定的振動傳播方向,得到該振動傳播的位置時間圖。
- 一種介質黏彈性定量裝置,其特徵在於,該裝置包括: 一影像模組,用於在一介質受到振動激發後,得到一振動傳播的一位置時間圖; 一確定模組,用於沿該位置時間圖上的預設角度範圍內的各個角度進行一角度投影,確定訊號能量最大的角度所對應的該位置時間圖的一斜率; 一定量模組,用於根據該斜率,得到介質的黏彈性參數。
- 如申請專利範圍第7項所述的裝置,其中,該確定模組,包括: 一計算子模組,用於對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 一確定子模組,用於將該積分子模組計算出的積分最大值的角度,確定為該位置時間圖的一斜率線的一斜率角;利用該斜率角確定該斜率線的斜率。
- 如申請專利範圍第7項所述的裝置,其中,該確定模組,包括: 一計算子模組,用於針對該位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣; 一確定子模組,用於得到該各個角度的一影像紋理特徵;利用該影像紋理特徵,將訊號能量最大的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的斜率。
- 如申請專利範圍第7項至第9項中任一項所述的裝置,其中,該裝置還包括: 一濾波模組,用於在進行該角度投影之前,濾除該位置時間圖中的反射波。
- 如申請專利範圍第7項至第9項中任一項所述的裝置,其中,該影像模組沿設定的振動傳播方向,得到該振動傳播的位置時間圖。
- 一種介質黏彈性定量裝置,其特徵在於,該裝置包括: 一記憶體,儲存一執行指令; 一處理器,被配置於讀取該執行指令,完成如下操作: 在一介質受到振動激發後,得到振動傳播的一位置時間圖; 沿該位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的一斜率; 根據該斜率,得到該介質的黏彈性參數。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710649552.9A CN107440740B (zh) | 2017-07-21 | 2017-07-21 | 介质粘弹性定量方法及装置 |
CN201710649552.9 | 2017-07-21 | ||
??201710649552.9 | 2017-07-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201907870A TW201907870A (zh) | 2019-03-01 |
TWI703962B true TWI703962B (zh) | 2020-09-11 |
Family
ID=60489354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW107121598A TWI703962B (zh) | 2017-07-21 | 2018-06-22 | 介質黏彈性定量方法及裝置 |
Country Status (13)
Country | Link |
---|---|
US (2) | US11719613B2 (zh) |
EP (1) | EP3656310B1 (zh) |
JP (1) | JP7006888B2 (zh) |
KR (1) | KR102317869B1 (zh) |
CN (1) | CN107440740B (zh) |
AU (1) | AU2018301988B2 (zh) |
BR (1) | BR112020001317A2 (zh) |
CA (1) | CA3070622C (zh) |
ES (1) | ES2941278T3 (zh) |
PL (1) | PL3656310T3 (zh) |
RU (1) | RU2744929C1 (zh) |
TW (1) | TWI703962B (zh) |
WO (1) | WO2019015397A1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107440740B (zh) * | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107505232B (zh) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
US11576051B2 (en) * | 2019-10-30 | 2023-02-07 | Samsung Electronics Co., Ltd. | Method and apparatus for indexing of SS/PBCH block in unlicensed spectrum |
KR20210115177A (ko) | 2020-03-12 | 2021-09-27 | 주식회사 엘지에너지솔루션 | 조립이 용이한 전지 모듈 및 이를 포함하는 전지 팩 |
Family Cites Families (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100081929A1 (en) * | 2005-09-14 | 2010-04-01 | Takao Suzuki | Position tracking method, position tracking device, and ultrasonograph |
FR2899336B1 (fr) * | 2006-03-29 | 2008-07-04 | Super Sonic Imagine | Procede et dispositif pour l'imagerie d'un milieu viscoelastique |
JP4797194B2 (ja) * | 2006-05-09 | 2011-10-19 | 独立行政法人産業技術総合研究所 | 超音波断層画像による生体組織評価システム |
FR2902879B1 (fr) | 2006-06-22 | 2008-10-10 | Michelin Soc Tech | Rheometre orthogonal |
CN100470215C (zh) * | 2007-01-15 | 2009-03-18 | 国家海洋技术中心 | 视频波浪测量方法和测量系统 |
US8118744B2 (en) * | 2007-02-09 | 2012-02-21 | Duke University | Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction |
JP5559788B2 (ja) * | 2009-07-07 | 2014-07-23 | 株式会社日立メディコ | 超音波診断装置 |
CN102283679B (zh) * | 2011-08-04 | 2014-05-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
MX344792B (es) * | 2011-09-25 | 2017-01-06 | Theranos Inc | Sistemas y métodos para múltiples análisis. |
WO2013183432A1 (ja) | 2012-06-07 | 2013-12-12 | 日立アロカメディカル株式会社 | 関心領域設定方法及び超音波診断装置 |
US9546981B2 (en) * | 2012-08-22 | 2017-01-17 | Neotek Biosciences Co. Ltd. | Device and method for performing blood thromboelastographic assays by magnetic sensing |
JP5753633B2 (ja) * | 2013-03-05 | 2015-07-22 | 日立アロカメディカル株式会社 | 超音波診断装置、及び送受信方法 |
US9345448B2 (en) * | 2013-04-01 | 2016-05-24 | Mayo Foundation For Medical Education And Research | System and method for non-invasive determination of tissue wall viscoelasticity using ultrasound vibrometry |
JP6002845B2 (ja) * | 2013-07-19 | 2016-10-05 | メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ | 多方向波動場からのせん断波速度を測定するシステム及び医療機器の作動方法 |
CN104640506B (zh) * | 2013-09-18 | 2017-06-30 | 东芝医疗系统株式会社 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
FR3017041B1 (fr) * | 2014-01-31 | 2016-03-04 | Centre Nat Rech Scient | Procede et dispositif ultrasonore de caracterisation des milieux mous anisotropes, et ensemble de sonde ultrasonore pour un tel dispositif de caracterisation |
BR112017008908B1 (pt) * | 2014-10-28 | 2022-10-18 | Supersonic Imagine | Métodos e aparelhos de formação de imagem para realizar formação de imagem por elastografia com onda de cisalhamento |
US9907539B2 (en) * | 2015-01-12 | 2018-03-06 | Siemens Medical Solutions Usa, Inc. | Sparse tracking in acoustic radiation force impulse imaging |
US11166698B2 (en) * | 2015-01-30 | 2021-11-09 | Canon Medical Systems Corporation | Ultrasonic diagnostic apparatus |
US9726647B2 (en) * | 2015-03-17 | 2017-08-08 | Hemosonics, Llc | Determining mechanical properties via ultrasound-induced resonance |
US10582911B2 (en) * | 2015-08-11 | 2020-03-10 | Siemens Medical Solutions Usa, Inc. | Adaptive motion estimation in acoustic radiation force imaging |
CN105232085B (zh) * | 2015-11-18 | 2017-10-31 | 中国人民解放军第三军医大学第三附属医院 | 基于动态孔径控制的超声剪切波弹性成像方法 |
US11071524B2 (en) * | 2015-12-04 | 2021-07-27 | Canon Medical Systems Corporation | Analyzing apparatus |
CN105455851B (zh) * | 2015-12-24 | 2018-03-13 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
CN105844645B (zh) * | 2016-03-31 | 2019-11-12 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
CN106037816A (zh) * | 2016-05-19 | 2016-10-26 | 深圳大学 | 一种超声粘弹性测量方法及系统 |
CN106419960A (zh) * | 2016-12-04 | 2017-02-22 | 无锡圣诺亚科技有限公司 | 具备超声穿刺导航的超声仪 |
CN106805997B (zh) * | 2016-12-26 | 2020-08-07 | 乐普(北京)医疗器械股份有限公司 | 一种弹性成像方法和装置 |
CN106618635B (zh) * | 2017-01-12 | 2019-11-08 | 清华大学 | 剪切波弹性成像方法和装置 |
CN107440740B (zh) | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107505232B (zh) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
-
2017
- 2017-07-21 CN CN201710649552.9A patent/CN107440740B/zh active Active
-
2018
- 2018-05-25 PL PL18835631.5T patent/PL3656310T3/pl unknown
- 2018-05-25 WO PCT/CN2018/088405 patent/WO2019015397A1/zh unknown
- 2018-05-25 ES ES18835631T patent/ES2941278T3/es active Active
- 2018-05-25 AU AU2018301988A patent/AU2018301988B2/en active Active
- 2018-05-25 EP EP18835631.5A patent/EP3656310B1/en active Active
- 2018-05-25 KR KR1020207004555A patent/KR102317869B1/ko active IP Right Grant
- 2018-05-25 JP JP2020503047A patent/JP7006888B2/ja active Active
- 2018-05-25 CA CA3070622A patent/CA3070622C/en active Active
- 2018-05-25 RU RU2020106885A patent/RU2744929C1/ru active
- 2018-05-25 BR BR112020001317-0A patent/BR112020001317A2/pt unknown
- 2018-06-22 TW TW107121598A patent/TWI703962B/zh active
-
2020
- 2020-01-17 US US16/746,698 patent/US11719613B2/en active Active
-
2023
- 2023-06-14 US US18/209,961 patent/US20230324273A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
EP3656310A4 (en) | 2021-04-07 |
WO2019015397A1 (zh) | 2019-01-24 |
EP3656310B1 (en) | 2023-03-15 |
BR112020001317A2 (pt) | 2020-08-11 |
RU2744929C1 (ru) | 2021-03-17 |
PL3656310T3 (pl) | 2023-07-10 |
JP7006888B2 (ja) | 2022-01-24 |
KR20200030580A (ko) | 2020-03-20 |
AU2018301988B2 (en) | 2021-03-18 |
CN107440740A (zh) | 2017-12-08 |
AU2018301988A1 (en) | 2020-02-20 |
EP3656310A1 (en) | 2020-05-27 |
CN107440740B (zh) | 2021-06-25 |
CA3070622A1 (en) | 2019-01-24 |
US11719613B2 (en) | 2023-08-08 |
US20200150013A1 (en) | 2020-05-14 |
KR102317869B1 (ko) | 2021-10-25 |
CA3070622C (en) | 2024-01-23 |
JP2020527723A (ja) | 2020-09-10 |
US20230324273A1 (en) | 2023-10-12 |
ES2941278T3 (es) | 2023-05-19 |
TW201907870A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI703962B (zh) | 介質黏彈性定量方法及裝置 | |
KR101914021B1 (ko) | 점탄성 매질의 점탄성 파라미터 검출 방법 및 기기 | |
CN107049360A (zh) | 剪切波弹性成像方法和装置 | |
US11517289B2 (en) | Method and apparatus for acquiring motion information | |
Jarvis et al. | Application of the distributed point source method to rough surface scattering and ultrasonic wall thickness measurement | |
CN106618635A (zh) | 剪切波弹性成像方法和装置 | |
CN105212968A (zh) | 弹性检测方法和设备 | |
TWI704906B (zh) | 介質黏彈性的測量方法和裝置 | |
JP5456367B2 (ja) | フェーズドアレイ開口合成処理方法 | |
TWI704905B (zh) | 介質黏彈性的測量方法和裝置 | |
CN110897656B (zh) | 超声剪切波弹性测量方法及装置 | |
Xu et al. | Research on ultrasonic damage imaging based on inverse scattering model in frequency domain | |
Schmerr, Jr et al. | Ultrasonic measurement models for imaging with phased arrays |