TWI686180B - 運動資訊獲取方法及裝置 - Google Patents
運動資訊獲取方法及裝置 Download PDFInfo
- Publication number
- TWI686180B TWI686180B TW107121597A TW107121597A TWI686180B TW I686180 B TWI686180 B TW I686180B TW 107121597 A TW107121597 A TW 107121597A TW 107121597 A TW107121597 A TW 107121597A TW I686180 B TWI686180 B TW I686180B
- Authority
- TW
- Taiwan
- Prior art keywords
- signal
- slope
- vibration
- frequency domain
- angle
- 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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices 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
-
- 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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0048—Detecting, measuring or recording by applying mechanical forces or stimuli
- A61B5/0051—Detecting, measuring or recording by applying mechanical forces or stimuli by applying vibrations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0097—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying acoustic waves and detecting light, i.e. acoustooptic measurements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
- A61B5/4222—Evaluating particular parts, e.g. particular organs
- A61B5/4244—Evaluating particular parts, e.g. particular organs liver
-
- 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/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
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/008—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means by using ultrasonic waves
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/262—Analysis of motion using transform domain methods, e.g. Fourier domain methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Biochemistry (AREA)
- Physiology (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Gastroenterology & Hepatology (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Endocrinology (AREA)
Abstract
本發明揭露了運動資訊獲取方法及裝置。對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號,接著去除該頻域訊號中位於設定振動速度範圍之外的訊號,即只保留振動的訊號,再沿設定的振動傳播方向得到位置時間圖。無需藉由複雜的計算對振動的傳播進行運動估計,只需藉由在頻域進行處理來判斷振動的有無,即可得到位置時間圖,是一種高效的運動資訊獲取方式。
Description
本發明涉及測量技術領域,特別涉及運動資訊獲取方法及裝置。
對介質進行振動激發時,振動在介質中的傳播特性與介質的黏彈性有關,藉由測量振動的傳播特性,可以對介質的黏彈性進行定量。獲取振動的傳播特性,又需要先利用針對該振動的檢測訊號,來獲取該振動的運動資訊。
上述原理目前已被應用至多個技術領域,以醫學檢測為例,檢測肝臟、甲狀腺和肌肉等器官或組織時,藉由對介質的黏彈性進行定量,可以對病變進行定位。
因此,如何高效的獲取在介質中傳播的振動的運動資訊,是一個需要解決的問題。
本發明實施例提供了運動資訊獲取方法及裝置。為了對披露的實施例的一些方面有一個基本的理解,下面給出了簡單的概括。該概括部分不是泛泛評述,也不是要確定關鍵/重要組成元素或描繪這些實施例的保護範圍。其唯一目的是用簡單的形式呈現一些概念,以此作為後面的詳細說明的序言。
根據本發明實施例的第一方面,提供了一種運動資訊獲取方法,該方法包括: 對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號; 去除該頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號; 利用該處理訊號得到該振動的位置時間圖。
根據該方法,作為可選的實施例一,該去除該頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號,包括: 對該頻域訊號進行濾波或特徵值選取,得到處理訊號; 該濾波的參數與該設定振動速度範圍相關; 該特徵值選取與該設定振動速度範圍相關。
根據該方法,作為可選的實施例二,該利用該處理訊號得到該振動的位置時間圖,包括: 根據設定的振動傳播方向,利用該處理訊號得到該振動的位置時間圖。
根據該方法,作為可選的實施例三,該方法還包括: 對該位置時間圖進行影像分割; 提取影像特徵; 利用該影像特徵進行線性擬合,得到該位置時間圖的斜率線的斜率; 根據該斜率,計算該介質的黏彈性參數。
根據該方法,作為可選的實施例四,該方法還包括: 沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率; 根據該斜率,得到該介質的黏彈性參數。
根據該實施例四,作為可選的實施例五,該沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括: 對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 將積分值最大的角度,確定為該位置時間圖的斜率線的斜率角; 利用該斜率角確定該斜率線的斜率。
根據本發明實施例的第二方面,提供了一種運動資訊獲取裝置,該裝置包括: 第一處理模組,用於對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號; 第二處理模組,用於去除該頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號; 獲取模組,用於利用該處理訊號得到該振動的位置時間圖。
根據該方法,作為可選的實施例一,該第二處理模組對該頻域訊號進行濾波或特徵值選取,得到處理訊號; 該濾波的參數與該設定振動速度範圍相關; 該特徵值選取與該設定振動速度範圍相關。
根據該裝置,作為可選的實施例二,該獲取模組根據設定的振動傳播方向,利用該處理訊號得到該振動的位置時間圖。
根據該裝置,作為可選的實施例三,該裝置還包括: 黏彈性定量模組,用於對該位置時間圖進行影像分割;提取影像特徵;利用該影像特徵進行線性擬合,得到該位置時間圖的斜率線的斜率; 根據該斜率,計算該介質的黏彈性參數。
根據該裝置,作為可選的實施例四,該裝置還包括: 黏彈性定量模組,用於沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率;根據該斜率,得到該介質的黏彈性參數。
根據該實施例四,作為可選的實施例五,該黏彈性定量模組,包括: 計算子模組,用於對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算; 確定子模組,用於將該計算子模組計算出的積分最大值的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的斜率; 定量子模組,用於根據該斜率,得到該介質的黏彈性參數。
根據本發明實施例的第三方面,提供了一種運動資訊獲取裝置,該裝置包括: 記憶體,儲存執行指令; 處理器,被配置於讀取該執行指令,完成如下操作: 用於對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號; 去除該頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號; 利用該處理訊號得到該振動的位置時間圖。
本發明實施例提供的技術方案可以包括以下有益效果: 對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號,接著去除該頻域訊號中位於設定振動速度範圍之外的訊號,即只保留振動的訊號,再得到振動的位置時間圖。無需藉由複雜的計算對振動的傳播進行運動估計,只需藉由在頻域進行處理來判斷振動的有無,即可得到位置時間圖,是一種高效的運動資訊獲取方式。
應當理解的是,以上的一般描述和後文的細節描述僅是示例性和解釋性的,並不能限制本發明。
以下描述和附圖充分地示出本發明的具體實施方案,以使本領域的技術人員能夠實踐它們。實施例僅代表可能的變化。除非明確要求,否則單獨的部件和功能是可選的,並且操作的順序可以變化。一些實施方案的部分和特徵可以被包括在或替換其他實施方案的部分和特徵。本發明的實施方案的範圍包括申請專利範圍的整個範圍,以及申請專利範圍的所有可獲得的等同物。在本文中,各實施方案可以被單獨地或總地用術語“發明”來表示,這僅僅是為了方便,並且如果事實上揭露了超過一個的發明,不是要自動地限制該應用的範圍為任何單個發明或發明構思。本文中,諸如第一和第二等之類的關係術語僅僅用於將一個實體或者操作與另一個實體或操作區分開來,而不要求或者暗示這些實體或操作之間存在任何實際的關係或者順序。而且,術語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含,從而使得包括一系列要素的程序、方法或者裝置不僅包括那些要素,而且還包括沒有明確列出的其他要素。本文中各個實施例採用遞進的方式描述,每個實施例重點說明的都是與其他實施例的不同之處,各個實施例之間相同相似部分互相參見即可。對於實施例揭露的結構、產品等而言,由於其與實施例揭露的部分相對應,所以描述的比較簡單,相關之處參見方法部分說明即可。
第1圖是根據一示例性實施例示出的運動資訊獲取方法的流程圖。如第1圖所示,該方法包括如下步驟。
在步驟11中,對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號。
在步驟12中,去除頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號。
在步驟13中,利用處理訊號得到振動的位置時間圖。
介質受到振動激發後,振動在介質中傳播,在不同的時刻,波前會沿著傳播方向到達不同的位置,這種位置和時刻的對應關係就是振動的運動資訊。目前常用的運動資訊獲取方法,利用的是在介質振動時,針對該振動的檢測訊號會發生相位去相關,根據這種相位去相關的特性,可以藉由互相關、自相關、光流等演算法,得到該振動的位置時間圖。上述演算法可以選擇任意基於塊匹配的方法。這種常用的運動資訊獲取方法,相當於需要藉由介質的位移和應變等資訊先對振動的傳播進行運動估計,然後才能得到該振動的位置時間圖。
在本示例性實施例中,利用檢測訊號中包含振動資訊產生多普勒效應的特性,對檢測訊號在成像時間維度進行頻域變換,得到頻域訊號,接著去除該頻域訊號中位於設定振動速度範圍之外的訊號,即去除相對靜止或振動速度低的訊號,再得到振動的位置時間圖。可見,本示例性實施例中的運動資訊獲取方法,不需要複雜的計算,藉由在頻域進行處理,就可以得到不以位移或應變為特徵的位置時間圖。這種方法相當於無需先對振動的傳播進行運動估計,只需判斷振動的有無即可得到位置時間圖,是一種高效的運動資訊的獲取方法。
在一示例性實施例中,上述步驟11中,可以採用多種方式進行頻域變換,例如採用傅裡葉變換或奇異值分解來進行頻域變換。
在一示例性實施例中,上述步驟12中,去除頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號,可以藉由對頻域訊號濾波或進行特徵值選取來實現。以採用濾波實現上述去除操作為例,為了只保留振動速度在設定範圍之內的訊號,考慮訊號在空間、時間上的取樣速率,結合設定的振動速度範圍,例如0.1米每秒至30米每秒,來設置濾波器的濾波參數,就可以基於該濾波去除頻域訊號中位於設定振動速度範圍之外的訊號。在採用特徵值選取實現上述去除操作時,同樣可以設置特徵值的選取與設定的振動速度範圍相關,從而實現去除頻域訊號中位於設定振動速度範圍之外的訊號。藉由步驟12,可以只將振動的訊號保留下來,提高後續形成位置時間圖時的準確度。
在一示例性實施例中,藉由機械振動、聲輻射力或其他可以產生振動的方式,對介質進行振動激發後,介質產生振動,振動在介質中傳播。由於上述振動在介質中的傳播速度有限,因此可以利用檢測訊號對介質進行動態成像。上述檢測波可以是光波、超音波等。上述動態成像可以是一維成像、二維成像或三維成像等。無論成像方式是哪種,都可以根據設定的振動傳播方向,利用藉由去除操作之後的處理訊號得到振動的位置時間圖。上述設定的振動傳播方向,在振動只在一個傳播方向上傳播時,為振動的實際傳播方向,在振動在多個傳播方向上傳播時,為選擇出的某一個傳播方向。例如,當介質為一均勻薄片時,對介質進行振動激發後,振動將沿該薄片的延展方向進行傳播,此時設定的振動傳播方向就是該振動的實際傳播方向。又例如,當介質為立體不規則形狀時,振動傳播的波前為立體形狀,例如振動傳播的波前為橢球,則沿不同的振動傳播方向得到的位置時間圖是不同的,此時設定的振動傳播方向就是選擇出的某一個感興趣的傳播方向。上述感興趣的傳播方向根據實際要測量的方向來確定,例如可以是振動傳播最快的方向、振動傳播最慢的方向和振動傳播速度為某一區間的方向中的至少一個。
在一示例性實施例中,第1圖所示的方法中還可以包括介質黏彈性定量的步驟。第2圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖,該流程基於第1圖所示流程實現,包括如下步驟。
在步驟21中,沿位置時間圖上的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的位置時間圖的斜率。
預設角度範圍,指根據實際情況,選取的進行角度投影的角度範圍。作為一種可選的實施方式,預設角度範圍可以是360度,此時需進行全角度的角度投影。作為另一種可選的實施方式,根據得到的位置時間圖的特性,來選擇進行角度投影的角度範圍。步驟11中得到的位置時間圖的橫軸指示時間、縱軸指示位置,如果振動只由振動激發起點處向遠處傳播,則振動傳播的速度無限大時,在位置時間圖上接近於與縱軸平行的直線,而振動傳播的速度無限小時,在位置時間圖上接近於與橫軸平行的直線,此時預設角度範圍是90度即可滿足需求,而無需做全角度投影,進而提高介質黏彈性定量的效率。如果振動除了由振動激發起點處向遠處傳播,還可沿相反方向繼續傳播,則預設角度範圍可以是180度。至於預設角度範圍的實際起始點和終止點,在直角坐標系保持不變的情況下,與0度的起始點及逆時針或順時針的旋轉方向相關,可以按需要選取,只要保證預設角度範圍即可。
各個角度,指在預設角度範圍內進行角度投影的各個角度。具體角度的選取根據時間精度需求和計算速度需求來確定,時間精度要求越高時,角度選取的精度要求越高,計算速度需求越高時,角度選取的精度要求越低。例如,可以從0.0.1度到1度的範圍內選擇。
角度投影,指標對設定角度進行影像特徵識別或提取,以確定訊號能量最大的角度。
在步驟22中,根據斜率,得到介質的黏彈性參數。
黏彈性參數包括黏性參數和彈性參數中的至少一項。
位置時間圖的斜率由單位時間的振動傳播的距離決定,即介質中振動傳播的速度。在均勻的介質中,振動傳播的速度和介質的黏彈性有關,在得到位置時間圖的斜率之後,就可以定量計算出介質的黏彈性參數,因此如何高效準確的獲得上述斜率成為介質黏彈性定量的關鍵。本示例性實施例利用角度投影確定訊號能量最大的角度,由於該訊號能量最大的角度對應位置時間圖的斜率,即相當於獲得了位置時間圖的斜率。這種方法無需從位置時間圖中選取波峰、波谷或振動的某一相位作為特徵點來計算位置時間圖的斜率,不會受到雜訊干擾且計算量小,是一種高效準確的進行介質黏彈性定量的方法。
由於振動在介質中傳播時,遇介質的邊緣或異物時會產生反射波,為提高後續處理的精度,如第3圖所示,在進行角度投影之前,還可以包括步驟21’,即對位置時間圖中的反射波進行濾除。濾除的方式可以有很多種,方向濾波是其中一種實現方式。
作為一種可選的實施方式,藉由角度投影來確定訊號能量最大的角度,進而得到位置時間圖的斜率,可以藉由積分計算來實現。例如,對位置時間圖進行沿預設角度範圍內的各個角度的積分計算,當積分角度和振動傳播方向一致時能量聚集,此時得到的積分值最大,因此將積分值最大的角度為位置時間圖的斜率線的斜率角。根據得到的斜率角,再結合位置和時間資訊,即可以得到位置時間圖的斜率線的斜率。上述積分計算又稱為Radon變換。
作為另一種可選的實施方式,藉由角度投影來確定訊號能量最大的角度,進而得到位置時間圖的斜率,還可以藉由計算灰階共生矩陣來實現。由於藉由計算灰階共生矩陣可以獲得影像紋理特徵,影像紋理特徵又可以體現訊號能量的大小,因此可以用灰階共生矩陣來獲取訊號能量最大的角度的資訊。基於上述原理,藉由角度投影來確定訊號能量最大的角度,進而得到位置時間圖的斜率,可以藉由灰階共生矩陣的計算來實現。例如,針對位置時間圖,先沿預設角度範圍內的各個角度計算灰階共生矩陣。然後,利用灰階共生矩陣,得到各個角度的影像紋理特徵。接著,利用影像紋理特徵,將訊號能量最大的角度,確定為位置時間圖的斜率線的斜率角。最後,利用斜率角確定斜率線的斜率。
在一示例性實施例中,第1圖所示的方法中還可以包括介質黏彈性定量的步驟。第4圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖,該流程基於第1圖所示流程實現,包括如下步驟。
在步驟41中,對位置時間圖進行影像分割。
在步驟42中,提取影像特徵。
影像特徵可以是中軸線、峰值、谷值和零交火點中的至少一個。上述中軸線指的是位置時間圖上圖紋的骨架,上述零交火點指的是斜率最大值點或二階導的最大值點。提取出的影像特徵中均包含有振動的資訊。
在步驟43中,利用影像特徵進行線性擬合,得到位置時間圖的斜率線的斜率。
在步驟44中,根據斜率,計算介質的黏彈性參數。
在前述兩個介質黏彈性定量的示例性實施例中,根據力學原理,介質的黏彈性決定了振動在該介質中的傳播速度,因此藉由得到位置時間圖的斜率,可以獲知振動在介質中傳播的速度,進而根據力學原理,可以定量得出介質的黏彈性參數。這裡的黏彈性參數可以包括剪切模量、楊氏模量、黏性模量、剪切黏彈性、剪切黏度、機械阻抗、機械鬆弛時間、各向異性等。
可選的,在使用線性擬合來實現介質黏彈性定量時,也可以先濾除位置時間圖中的反射波,以達到更精確的定量效果。
下面以一種具體的應用場景給出本發明實施例中運動資訊獲取方法的應用。
在對人體肝臟等黏黏彈性介質進行無損黏彈性檢測時,需要對介質黏彈性進行定量,定量之前就需要先獲取運動資訊。檢測裝置中包括激發裝置和成像裝置,其中激發裝置對待檢測介質進行振動激發,成像裝置利用超音波對振動激發後的介質進行成像。振動在介質中傳播時,在不同的時刻,波前會沿著傳播方向到達不同的位置,形成位置時間圖。上述波前可以是波峰、波谷、或振動的同一相位中的一種。
如第5圖所示,在這種具體應用場景的運動資訊獲取方法可以包括如下步驟。
在步驟51中,對介質進行振動激發。
在步驟52中,對在介質中傳播的振動的超聲檢測訊號進行頻域變換,得到頻域訊號。
在步驟53中,對頻域訊號進行濾波或特徵值選取,得到處理訊號。
在步驟54中,根據設定的振動傳播方向,利用處理訊號得到振動的位置時間圖。
在步驟55中,利用線性擬合或Radon變換得到位置時間圖的斜率線的斜率。
在步驟56中,根據斜率,計算介質的黏彈性參數。
在以上運動資訊獲取方法的各個示例性實施例中,針對介質黏彈性定量的步驟,當設定的振動傳播方向為至少兩個時,每個設定的振動傳播方向會對應得到一個位置時間圖,進而會得到該位置時間圖對應的介質的黏彈性參數。綜合得到的至少兩套黏彈性參數,可以更全面的評價介質的黏彈性。
上述給出的運動資訊獲取方法的各個示例性實施例,可以根據情況進行組合,這裡並不限定各個示例性實施例之間的組合關係。
第6圖是一示例性實施例示出的運動資訊獲取裝置的框圖,該裝置可以位於介質黏彈性檢測裝置的控制主機中,例如在醫療檢測領域,可以位於肝臟無損檢測裝置的控制主機中。該裝置還可以位於雲端,介質黏彈性檢測裝置的檢測資料需要在雲端進行處理。
第6圖所示的裝置包括:第一處理模組61、第二處理模組62和獲取模組63。
第一處理模組61,用於對在介質中傳播的振動的檢測訊號進行頻域變換,得到頻域訊號。第一處理模組61在進行頻域變換時,可以採用傅裡葉變換、奇異值分解等多種方法。
第二處理模組62,用於去除頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號。
獲取模組63,用於利用處理訊號得到振動的位置時間圖。
在一示例性實施例中,第二處理模組62對頻域訊號進行濾波或特徵值選取,得到處理訊號。其中,濾波的參數與設定振動速度範圍相關,特徵值選取與設定振動速度範圍相關。
在一示例性實施例中,獲取模組63根據設定的振動傳播方向,利用處理訊號得到振動的位置時間圖。
在一示例性實施例中,如第7圖所示,第6圖所示的裝置中還可以包括:黏彈性定量模組64,用於對位置時間圖進行影像分割;提取影像特徵;利用影像特徵進行線性擬合,得到位置時間圖的斜率線的斜率;根據斜率,計算介質的黏彈性參數。此處影像特徵的含義與前文方法中所述的相同。
作為另一種可選的實施方式,黏彈性定量模組64還可以使用角度投影來實現相同功能。黏彈性定量模組64,用於沿位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的位置時間圖的斜率;根據斜率,得到介質的黏彈性參數。
進一步可選的,如第8圖所示,黏彈性定量模組64可以包括:計算子模組641和確定子模組642。
計算子模組641,用於對位置時間圖進行沿預設角度範圍內的各個角度的積分計算。
確定子模組642,用於將計算子模組641計算出的積分最大值的角度,確定為位置時間圖的斜率線的斜率角。利用斜率角確定該斜率線的斜率。
定量子模組643,用於根據斜率,得到介質的黏彈性參數。
作為另一種可選的實施方式,除了積分計算的方式之外,黏彈性定量模組64還可以藉由計算灰階共生矩陣來確定斜率。此時,計算子模組641用於針對該位置時間圖,沿預設角度範圍內的各個角度計算灰階共生矩陣。確定子模組642用於得到各個角度的影像紋理特徵;利用影像紋理特徵,將訊號能量最大的角度,確定為位置時間圖的斜率線的斜率角;利用斜率角確定斜率線的斜率。
進一步可選的,如第9圖所示,運動資訊獲取裝置中還包括:濾波模組65,用於在黏彈性定量模組64進行角度投影之前,濾除位置時間圖中的反射波。當然,在採用線性擬合來實現介質黏彈性定量時,也可以先由濾波模組65來濾除位置時間圖中的反射波。
第10圖是根據一示例性實施例示出的運動資訊獲取裝置的框圖,該裝置可以位於介質黏彈性檢測裝置的控制主機中,例如在醫療檢測領域,可以位於肝臟無損檢測裝置的控制主機中。該裝置還可以位於雲端,介質黏彈性檢測裝置的檢測資料需要在雲端進行處理。
第10圖所示的裝置包括:記憶體101和處理器102。
記憶體101中儲存有執行指令。
處理器102,被配置於讀取記憶體101中的執行指令,執行前文該運動資訊獲取方法的各示例性實施例中的部分或全部步驟。處理器102可以由晶片實現。
如果第10圖所示的運動資訊獲取裝置位於介質黏彈性檢測裝置的控制主機中,可以藉由匯流排、無線等方式與介質黏彈性定量裝置中的激發裝置、成像裝置進行連接,此時該裝置上具備與實現上述連接的介面及相應通訊機制。
如果第10圖所示的運動資訊獲取裝置位於雲端,可以藉由網路與介質黏彈性檢測裝置進行通訊。
應當理解的是,本發明並不侷限於上面已經描述並在附圖中示出的流程及結構,並且可以在不脫離其範圍進行各種修改和改變。本發明的範圍僅由所附的申請專利範圍來限制。
11、12、13、21、22、21’、41、42、43、44、51、52、53、54、55、56‧‧‧步驟61、62‧‧‧處理模組63‧‧‧獲取模組64‧‧‧黏彈性定量模組65‧‧‧濾波模組101‧‧‧記憶體102‧‧‧處理器641‧‧‧計算子模組
642:確定子模組
643:定量子模組
此處的附圖被併入說明書中並構成本說明書的一部分,示出了符合本發明的實施例,並與說明書一起用於解釋本發明的原理。 第1圖是根據一示例性實施例示出的運動資訊獲取方法的流程圖; 第2圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第3圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第4圖是根據一示例性實施例示出的介質黏彈性定量方法的流程圖; 第5圖是根據一示例性實施例示出的運動資訊獲取方法的流程圖; 第6圖是根據一示例性實施例示出的運動資訊獲取裝置的框圖; 第7圖是根據一示例性實施例示出的運動資訊獲取裝置的框圖; 第8圖是第7圖中所示黏彈性定量模組的框圖; 第9圖是根據一示例性實施例示出的運動資訊獲取裝置的框圖; 第10圖是根據一示例性實施例示出的運動資訊獲取裝置的框圖。
11、12、13‧‧‧步驟
Claims (13)
- 一種運動資訊獲取方法,其特徵在於,該方法包括:對在一介質中傳播的一振動的檢測訊號進行頻域變換,得到一頻域訊號;去除該頻域訊號中位於一設定振動速度範圍之外的訊號,得到一處理訊號,該處理訊號為一頻域處理訊號;以及利用該處理訊號,藉由在一頻域進行處理,得到該振動的一位置時間圖。
- 如申請專利範圍第1項所述的方法,其中,該去除該頻域訊號中位於設定振動速度範圍之外的訊號,得到處理訊號,包括:對該頻域訊號進行一濾波或一特徵值選取,得到處理訊號;該濾波的參數與該設定振動速度範圍相關;該特徵值選取與該設定振動速度範圍相關。
- 如申請專利範圍第1項所述的方法,其中,該利用該處理訊號得到該振動的位置時間圖,包括:根據設定的振動傳播方向,利用該處理訊號得到該振動的位置時間圖。
- 如申請專利範圍第1項所述的方法,其中,該方法還包括:對該位置時間圖進行影像分割;提取一影像特徵;利用該影像特徵進行線性擬合,得到該位置時間圖的斜率線的一斜率;根據該斜率,計算該介質的黏彈性參數。
- 如申請專利範圍第1項所述的方法,其中,該方法還包括:沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率;根據該斜率,得到該介質的黏彈性參數。
- 如申請專利範圍第5項所述的方法,其中,該沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的斜率,包括:對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算;將積分值最大的角度,確定為該位置時間圖的一斜率線的一斜率角;利用該斜率角確定該斜率線的斜率。
- 一種運動資訊獲取裝置,其中,該裝置包括:一第一處理模組,用於對在介質中傳播的一振動的檢測訊號進行頻域變換,得到一頻域訊號;一第二處理模組,用於去除該頻域訊號中位於一設定振動速度範圍之外的訊號,得到一處理訊號,該處理訊號為一頻域處理訊號;以及一獲取模組,用於利用該處理訊號,藉由在一頻域進行處理,得到該振動的一位置時間圖。
- 如申請專利範圍第7項所述的裝置,其中,該第二處理模組對該頻域訊號進行一濾波或一特徵值選取,得到處理訊號;該濾波的參數與該設定振動速度範圍相關;該特徵值選取與該設定振動速度範圍相關。
- 如申請專利範圍第7項所述的裝置,其中,該獲取模組根據設定的振動傳播方向,利用該處理訊號得到該振動的位置時間圖。
- 如申請專利範圍第7所述的裝置,其中,該裝置還包括:一黏彈性定量模組,用於對該位置時間圖進行影像分割;提取一影像特徵;利用該影像特徵進行線性擬合,得到該位置時間圖的斜率線的一斜率;根據該斜率,計算該介質的黏彈性參數。
- 如申請專利範圍第7項所述的裝置,其中,該裝置還包括: 一黏彈性定量模組,用於沿該位置時間圖的預設角度範圍內的各個角度進行角度投影,確定訊號能量最大的角度所對應的該位置時間圖的一斜率;根據該斜率,得到該介質的黏彈性參數。
- 如申請專利範圍第11項所述的裝置,其中,該黏彈性定量模組,包括:一計算子模組,用於對該位置時間圖進行沿預設角度範圍內的各個角度的積分計算;一確定子模組,用於將該計算子模組計算出的積分最大值的角度,確定為該位置時間圖的斜率線的斜率角;利用該斜率角確定該斜率線的一斜率;一定量子模組,用於根據該斜率,得到該介質的黏彈性參數。
- 一種運動資訊獲取裝置,其中,該裝置包括:一記憶體,儲存一執行指令;一處理器,被配置於讀取該執行指令,完成如下操作:用於對在介質中傳播的振動的檢測訊號進行頻域變換,得到一頻域訊號;去除該頻域訊號中位於一設定振動速度範圍之外的訊號,得到一處理訊號,該處理訊號為一處理頻域訊號;利用該處理訊號,藉由在一頻域進行處理,得到該振動的位置時間圖。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
??201710649554.8 | 2017-07-21 | ||
CN201710649554.8 | 2017-07-21 | ||
CN201710649554.8A CN107505232B (zh) | 2017-07-21 | 2017-07-21 | 运动信息获取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201907869A TW201907869A (zh) | 2019-03-01 |
TWI686180B true TWI686180B (zh) | 2020-03-01 |
Family
ID=60690286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW107121597A TWI686180B (zh) | 2017-07-21 | 2018-06-22 | 運動資訊獲取方法及裝置 |
Country Status (11)
Country | Link |
---|---|
US (1) | US11517289B2 (zh) |
EP (1) | EP3656311A4 (zh) |
JP (1) | JP7403163B2 (zh) |
KR (1) | KR102235737B1 (zh) |
CN (1) | CN107505232B (zh) |
AU (2) | AU2018301989A1 (zh) |
BR (1) | BR112019025264B1 (zh) |
CA (1) | CA3065575C (zh) |
RU (1) | RU2020106886A (zh) |
TW (1) | TWI686180B (zh) |
WO (1) | WO2019015398A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505232B (zh) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
CN107440740B (zh) | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN113075120B (zh) * | 2021-02-03 | 2022-08-30 | 湖南大学 | 一种实时土体类别识别方法、系统及一种盾构机 |
Family Cites Families (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS55165587A (en) * | 1979-06-11 | 1980-12-24 | Hitachi Cable | Method of forming rubber and plastic insulated cable connector |
US5099848A (en) * | 1990-11-02 | 1992-03-31 | University Of Rochester | Method and apparatus for breast imaging and tumor detection using modal vibration analysis |
US5810731A (en) | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
AU5247198A (en) | 1997-11-10 | 1999-05-31 | Medacoustics, Inc. | Non-invasive turbulent blood flow imaging system |
US20040034304A1 (en) * | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
US8211019B2 (en) | 2005-01-21 | 2012-07-03 | Chikayoshi Sumi | Clinical apparatuses |
WO2006124603A2 (en) | 2005-05-12 | 2006-11-23 | The Trustees Of Columbia University In The City Of New York | System and method for electromechanical wave imaging of body structures |
US7464583B1 (en) | 2005-06-10 | 2008-12-16 | Carnegie Mellon University | Methods and apparatuses using proximal probes |
CN100470215C (zh) * | 2007-01-15 | 2009-03-18 | 国家海洋技术中心 | 视频波浪测量方法和测量系统 |
EP2310829B1 (en) * | 2008-07-30 | 2016-04-06 | Centre Hospitalier De L'Universite de Montreal | A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance |
DE202008014792U1 (de) * | 2008-11-07 | 2010-03-25 | Qass Gmbh | Vorrichtung zum Bewerten von Zerspanungsprozessen |
CN101699280B (zh) | 2009-10-15 | 2011-08-17 | 北京索瑞特医学技术有限公司 | 超声无损检测粘弹性介质弹性的方法及其装置 |
US8469891B2 (en) * | 2011-02-17 | 2013-06-25 | Siemens Medical Solutions Usa, Inc. | Viscoelasticity measurement using amplitude-phase modulated ultrasound wave |
US9237878B2 (en) * | 2011-04-22 | 2016-01-19 | Mayo Foundation For Medical Education And Research | Generation and assessment of shear waves in elasticity imaging |
WO2014055410A1 (en) * | 2012-10-01 | 2014-04-10 | Mayo Foundation For Medical Education And Research | Shear wave attenuation from k-space analysis system |
KR101307526B1 (ko) * | 2012-11-19 | 2013-09-12 | (주)나다에스앤브이 | 실시간 진동 모니터링 장치 |
CN103054552B (zh) * | 2012-12-24 | 2014-12-10 | 深圳先进技术研究院 | 生物组织粘弹性测量方法和系统 |
US9895137B2 (en) | 2013-03-05 | 2018-02-20 | Hitachi, Ltd. | Ultrasonic diagnosis device and transmission/reception method |
JP6645961B2 (ja) | 2013-06-19 | 2020-02-14 | ザ ジェネラル ホスピタル コーポレイション | カテーテルによって全方位表示を得るための装置、デバイスおよび方法 |
WO2015040710A1 (ja) * | 2013-09-18 | 2015-03-26 | 株式会社 東芝 | 超音波診断装置、医用画像処理装置および医用画像処理方法 |
US10278671B2 (en) | 2013-09-30 | 2019-05-07 | Siemens Medical Solutions Usa, Inc. | Shear wave detection in medical ultrasound imaging |
JP6223169B2 (ja) | 2013-12-19 | 2017-11-01 | キヤノン株式会社 | 情報処理装置、情報処理方法およびプログラム |
CN103995053A (zh) * | 2014-03-13 | 2014-08-20 | 天津大学 | 声表面波技术探测各向同性材料表层杨氏模量的检测方法 |
US10722174B2 (en) | 2014-07-11 | 2020-07-28 | The Board Of Trustees Of The Leland Stanford Junior University | Skin-conformal sensors |
US10376242B2 (en) | 2015-04-16 | 2019-08-13 | Siemens Medical Solutions Usa, Inc. | Quantitative viscoelastic ultrasound imaging |
US10582911B2 (en) * | 2015-08-11 | 2020-03-10 | Siemens Medical Solutions Usa, Inc. | Adaptive motion estimation in acoustic radiation force imaging |
KR101620507B1 (ko) * | 2015-10-07 | 2016-05-13 | 한국지질자원연구원 | 볼 밀 시스템의 진동 모니터링 장치 및 방법 |
CN105455851B (zh) * | 2015-12-24 | 2018-03-13 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
CN107505233A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107440740B (zh) * | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107505232B (zh) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
-
2017
- 2017-07-21 CN CN201710649554.8A patent/CN107505232B/zh active Active
-
2018
- 2018-05-25 AU AU2018301989A patent/AU2018301989A1/en not_active Abandoned
- 2018-05-25 BR BR112019025264-9A patent/BR112019025264B1/pt active IP Right Grant
- 2018-05-25 WO PCT/CN2018/088406 patent/WO2019015398A1/zh unknown
- 2018-05-25 KR KR1020207004623A patent/KR102235737B1/ko active IP Right Grant
- 2018-05-25 EP EP18835507.7A patent/EP3656311A4/en active Pending
- 2018-05-25 RU RU2020106886A patent/RU2020106886A/ru unknown
- 2018-05-25 CA CA3065575A patent/CA3065575C/en active Active
- 2018-05-25 JP JP2020503053A patent/JP7403163B2/ja active Active
- 2018-06-22 TW TW107121597A patent/TWI686180B/zh active
-
2019
- 2019-11-19 US US16/689,022 patent/US11517289B2/en active Active
-
2022
- 2022-01-25 AU AU2022200463A patent/AU2022200463A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
WO2019015398A1 (zh) | 2019-01-24 |
US11517289B2 (en) | 2022-12-06 |
KR20200026302A (ko) | 2020-03-10 |
KR102235737B1 (ko) | 2021-04-01 |
CN107505232A (zh) | 2017-12-22 |
US20200100769A1 (en) | 2020-04-02 |
RU2020106886A (ru) | 2021-08-23 |
BR112019025264A2 (pt) | 2020-06-16 |
CA3065575A1 (en) | 2019-01-24 |
AU2022200463A1 (en) | 2022-02-17 |
BR112019025264B1 (pt) | 2021-05-11 |
JP2020527429A (ja) | 2020-09-10 |
TW201907869A (zh) | 2019-03-01 |
EP3656311A1 (en) | 2020-05-27 |
JP7403163B2 (ja) | 2023-12-22 |
AU2018301989A1 (en) | 2019-12-19 |
CA3065575C (en) | 2023-08-01 |
CN107505232B (zh) | 2019-09-03 |
EP3656311A4 (en) | 2021-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zahiri-Azar et al. | Motion estimation in ultrasound images using time domain cross correlation with prior estimates | |
US11719613B2 (en) | Method and device for quantifying viscoelasticity of a medium | |
TWI686180B (zh) | 運動資訊獲取方法及裝置 | |
CN107049360A (zh) | 剪切波弹性成像方法和装置 | |
Kijanka et al. | Robust phase velocity dispersion estimation of viscoelastic materials used for medical applications based on the multiple signal classification method | |
Kijanka et al. | Two point method for robust shear wave phase velocity dispersion estimation of viscoelastic materials | |
Kijanka et al. | Two-point frequency shift method for shear wave attenuation measurement | |
US20200163649A1 (en) | Estimating Phase Velocity Dispersion in Ultrasound Elastography Using a Multiple Signal Classification | |
Ghosh et al. | Modified error in constitutive equations (MECE) approach for ultrasound elastography | |
JP2020168232A (ja) | 超音波診断装置、信号処理装置、および、プログラム | |
Quang-Huy et al. | Two-dimensional viscoelastic imaging using an enhanced FDTD-AHI approach | |
de Almeida et al. | A Comparison of Time-Delay Estimators for Speckle Tracking Echocardiography | |
Touil et al. | Motion decorrelation in echocardiography: Analysis from a realistic simulation | |
Jiao et al. | A shear wave propagation tracking method based on Modal Assurance Criterion in acoustic radiation force impulse imaging | |
CN116784876A (zh) | 可穿戴实时监测的新型便携式剪切波超声检测系统及方法 |