TWI704906B - 介質黏彈性的測量方法和裝置 - Google Patents
介質黏彈性的測量方法和裝置 Download PDFInfo
- Publication number
- TWI704906B TWI704906B TW107121599A TW107121599A TWI704906B TW I704906 B TWI704906 B TW I704906B TW 107121599 A TW107121599 A TW 107121599A TW 107121599 A TW107121599 A TW 107121599A TW I704906 B TWI704906 B TW I704906B
- Authority
- TW
- Taiwan
- Prior art keywords
- echo signal
- time
- medium
- ultrasonic echo
- sample
- 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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本發明揭露了一種介質黏彈性的測量方法和裝置,涉及測量技術領域。包括,獲取介質在振動激發下的檢測波回波訊號,將檢測波回波訊號與資料集中儲存的檢測波回波訊號樣本匹配,確定介質的黏彈性資訊,其中,資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本。本發明通過在資料集中查找與待測量介質在振動激發下的檢測回波訊號相似程度最高的已知黏彈性的介質樣本在振動激發下的檢測波回波訊號樣本,來確定待測量介質的黏彈性,可以在不進行運動估計和特徵點選取的情況下,獲得介質的黏彈性,運算量大幅降低,測量速度顯著加快。
Description
本發明涉及測量技術領域,特別涉及介質黏彈性的測量方法和裝置。
機械振動或剪切波振動在介質中傳播時,在不同的時刻,波前會沿著傳播方向到達不同的位置,這樣波前在不同時間到達不同位置這個資訊可以組成位置時間圖,在均勻的組織中,位置時間圖通常為直的斜紋狀圖。針對這樣的斜紋狀圖,其斜率由單位時間的振動傳播的距離決定,即振動傳播的速度。在均勻的介質中,振動傳播的速度和介質黏彈性有關,因而對介質進行振動激發時,通過測量振動的傳播特性,可以對介質黏彈性進行測量。
先前技術中,在對介質黏彈性進行測量時,需要進行特徵點選取,對不同時刻、不同位置的特徵點進行線性擬合,得出動態超聲黏彈性成像的定量參數。通常情況下,特徵點由振動在介質中傳播時介質的位移或應變等決定,這時需要先進行運動估計,得到介質的位移或應變等資訊,然後再進行後續的運算。但這種方法存在運動估計所需的運算量較大,特徵點選取複雜等問題。
本發明實施例提供了一種介質黏彈性的測量方法和裝置。旨在解決先前技術中介質黏彈性測量中運動估計所需的運算量較大,特徵點選取複雜的問題。為了對披露的實施例的一些方面有一基本的理解,下面給出了簡單的概括。該概括部分不是泛泛評述,也不是要確定關鍵/重要組成元素或描繪這些實施例的保護範圍。其唯一目的是用簡單的形式呈現一些概念,以此作為後面的詳細說明的序言。
第一方面,本發明實施例提供了一種介質黏彈性的測量方法,包括: 獲取介質在振動激發下的檢測波回波訊號; 將檢測波回波訊號與資料集中儲存的檢測波回波訊號樣本匹配,確定介質的黏彈性資訊; 其中,資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本。
可選的,振動激發為剪切波激發,檢測波回波訊號為超聲回波訊號,資料集為超聲訊號資料集。
可選的,還包括: 在第一時刻t1獲取超聲回波資料; 將超聲回波資料載入傳播模式資料集中儲存的傳播模式資料,產生超聲回波訊號樣本; 產生超聲訊號資料集; 其中,傳播模式資料包括已知黏彈性資訊的介質樣本在給定的振動激發下的各個時刻的傳播狀態資料。
可選的,獲取介質在振動激發下的超聲回波訊號,包括:從第二時刻t2獲取介質在振動激發下的超聲回波訊號; 將檢測波回波訊號與資料集中儲存的檢測波回波訊號樣本匹配,包括:將從第二時刻t2獲取的超聲回波訊號中的一或複數與超聲回波訊號樣本進行匹配; 第二時刻t2在振動激發之後。
可選的,還包括: 在匹配之前,將超聲回波訊號與超聲回波訊號樣本進行時間對準及空間對準。
可選的,時間對準包括:確定時間差,時間差為第二時刻t2與振動激發起始時刻的差、及第二時刻t2與第一時刻t1的差中的較小值;根據時間差進行時間對準; 空間對準包括:確定介質的邊界和深度,進行空間對準;或,基於預設的深度進行空間對準。
可選的,傳播模式資料集包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
可選的,振動激發包括:連續或脈動方式的機械振動和聲輻射力脈動中的一項; 介質在振動激發下的運動範圍為0.01微米至10 毫米; 介質在振動激發下的運動頻率為20赫茲至2500 赫茲,持續時間為50 微秒至1秒。
可選的,檢測波回波訊號為超聲回波訊號;振動激發的脈動重複發生頻率為10至20000 赫茲;檢測波回波訊號包括至少連續採集的100訊框檢測波回波資料。
可選的,黏彈性包括黏性和彈性中的至少一項;當黏彈性為黏性時,黏彈性資訊包括黏性資訊、或黏性資訊和彈性資訊;當黏彈性為彈性時,黏彈性資訊包括彈性資訊、或黏性資訊和彈性資訊;當黏彈性為黏性和彈性時,黏彈性資訊包括黏性資訊和彈性資訊。
第二方面,本發明實施例提供了一種介質黏彈性的測量裝置,包括: 探頭,用於獲取介質在振動激發下的檢測波回波訊號; 記憶體,用於儲存資料集,資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本; 處理器,用於將檢測波回波訊號與資料集中的檢測波回波訊號樣本匹配,確定介質的黏彈性資訊。
可選的,振動激發為剪切波激發,檢測波回波訊號為超聲回波訊號,資料集為超聲訊號資料集。
可選的,探頭還用於,在第一時刻t1獲取超聲回波資料; 處理器還用於,將超聲回波資料載入傳播模式資料集中儲存的傳播模式資料,產生超聲回波訊號樣本;產生超聲訊號資料集; 其中,傳播模式資料包括已知黏彈性資訊的介質樣本在給定的振動激發下的各個時刻的傳播狀態資料。
可選的,探頭還用於,從第二時刻t2獲取介質在振動激發下的超聲回波訊號; 處理器還用於,將從第二時刻t2獲取的超聲回波訊號中的一或複數與超聲回波訊號樣本進行匹配; 其中,第二時刻t2在振動激發之後。
可選的,處理器還用於:在匹配之前,將超聲回波訊號與超聲回波訊號樣本進行時間對準及空間對準。
可選的,時間對準包括:確定時間差,時間差為第二時刻t2與振動激發起始時刻的差、及第二時刻t2與第一時刻t1的差中的較小值;根據時間差進行時間對準; 空間對準包括:確定介質的邊界和深度,進行空間對準;或,基於預設的深度進行空間對準。
可選的,傳播模式資料集包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
本發明實施例揭露的介質黏彈性的測量方法,通過在資料集中查找與待測量介質在振動激發下的檢測波回波訊號相似程度最高的已知黏彈性的介質樣本在振動激發下的檢測波回波訊號樣本,來確定待測量介質的黏彈性,可以在不進行運動估計和特徵點選取的情況下,獲得介質的黏彈性,運算量大幅降低,測量速度顯著加快。
應當理解的是,以上的一般描述和後文的細節描述僅是示例性和解釋性的,並不能限制本發明。
以下描述和附圖充分地示出本發明的具體實施方案,以使本領域的技術人員能夠實踐它們。實施例僅代表可能的變化。除非明確要求,否則單獨的部件和功能是可選的,並且操作的順序可以變化。一些實施方案的部分和特徵可以被包括在或替換其他實施方案的部分和特徵。本發明的實施方案的範圍包括申請專利範圍的整個範圍,以及申請專利範圍的所有可獲得的等同物。在本文中,各實施方案可以被單獨地或總地用術語“發明”來表示,這僅僅是為了方便,並且如果事實上揭露了超過一的發明,不是要自動地限制該應用的範圍為任何單個發明或發明構思。本文中,諸如第一和第二等之類的關係術語僅僅用於將一實體或者操作與另一實體或操作區分開來,而不要求或者暗示這些實體或操作之間存在任何實際的關係或者順序。而且,術語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含,從而使得包括一系列要素的程序、方法或者裝置不僅包括那些要素,而且還包括沒有明確列出的其他要素。本文中各個實施例採用遞進的方式描述,每個實施例重點說明的都是與其他實施例的不同之處,各個實施例之間相同相似部分互相參見即可。對於實施例揭露的結構、產品等而言,由於其與實施例揭露的部分相對應,所以描述的比較簡單,相關之處參見方法部分說明即可。
本發明實施例揭露了一種介質黏彈性的測量方法,如第1圖所示,包括如下步驟。
S101、獲取介質在振動激發下的檢測波回波訊號。
S102、將檢測波回波訊號與資料集中儲存的檢測波回波訊號樣本匹配,確定介質的黏彈性資訊,其中,資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本。
在S101中,振動激發包括但不限定為機械振動激發或聲輻射力激發。進一步的,機械振動激發和聲輻射力激發可以為連續或脈動等方式。示例性的,介質在振動激發下的運動範圍可以為0.01微米至10 毫米,介質在振動激發下的運動頻率為20赫茲至2500 赫茲,持續時間為50 微秒至1秒。
檢測波包括但不限定為超音波或光波,可選的,通過機械振動、聲輻射力或其他可以產生振動的方式,對介質進行振動激發後,介質產生振動,振動在介質中傳播。由於振動在介質中的傳播速度有限,因此可以利用超音波、光波或其他傳播速度較快的波對振動的傳播模式進行追蹤。通過獲取介質在振動激發下的檢測波回波訊號,反映振動在介質中的傳播模式。進一步的,檢測波回波訊號可以為超聲回波訊號,振動激發的脈動重複發生頻率為10至20000 赫茲。
一段時間內的檢測波回波訊號,由這段時間內的一系列時刻下的檢測波回波資料組成。可選的,檢測波回波訊號可以包括至少連續採集的100訊框檢測波回波資料。
介質黏彈性是決定振動在介質中傳播模式的重要影響因素。介質在振動激發下,在不同的時刻會產生相應的傳播狀態,傳播狀態沿時間的變化則構成傳播模式。傳播狀態是指某一時刻的介質運動狀態,傳播模式則是介質沿時間的傳播狀態的變化。
當介質具有不同的黏彈性時,則產生的傳播狀態不同,並進一步導致檢測波回波訊號的不同。
在S102中,資料集儲存有介質樣本的黏彈性、及介質樣本在振動激發下的檢測波回波訊號樣本,根據振動激發和檢測波回波訊號,在資料集中進行查找匹配即可確定相應的介質黏彈性。
可選的,介質黏彈性可以包括剪切模量、楊氏模量、剪切黏彈性、剪切黏度、機械阻抗、機械鬆弛時間和各向異性等中的至少一種或多種的組合。
可選的,黏彈性包括黏性和彈性中的至少一項。示例性的,當黏彈性為黏性時,黏彈性資訊可以包括黏性資訊、或黏性資訊和彈性資訊。示例性的,當黏彈性為彈性時,黏彈性資訊可以包括彈性資訊、或黏性資訊和彈性資訊。示例性的,當黏彈性為黏性和彈性時,黏彈性資訊可以包括黏性資訊和彈性資訊。
可見,本發明實施例揭露的介質黏彈性的測量方法,通過在資料集中查找與待測量介質在振動激發下的檢測波回波訊號相似程度最高的已知黏彈性的介質樣本在振動激發下的檢測波回波訊號樣本,來確定待測量介質的黏彈性,可以在不進行運動估計和特徵點選取的情況下,獲得介質的黏彈性,運算量大幅降低,測量速度顯著加快。
示例性的,當振動激發為剪切波激發,檢測波為超音波,檢測波回波訊號為超聲回波訊號,資料集為超聲訊號資料集時,本發明實施例還揭露了另一種介質黏彈性的測量方法,如第2圖所示,包括如下步驟。
S201、獲取介質在剪切波激發下的超聲回波訊號。
S202、將超聲回波訊號與資料集中儲存的超聲回波訊號樣本匹配,確定介質的黏彈性,其中,資料集儲存有介質樣本的黏彈性資訊、及介質樣本在剪切波激發下的超聲回波訊號樣本。
資料集可以儲存多組資料,每組資料分別記錄了介質在一種剪切波激發下產生的超聲回波訊號。進一步的,超聲訊號資料集可以儲存至少一種介質樣本分別在至少一種剪切波激發下產生的超聲回波訊號樣本,且資料集中,介質樣本的包括介質樣本黏彈性在內的一種或多種屬性參數已預先儲存。
獲得的超聲回波訊號與超聲回波訊號樣本進行匹配,以確定目標超聲回波訊號樣本的方法可以為任一種模式匹配方法,本發明實施例對此並不限定。
在S201中,可選的,本領域技術人員可以參照本發明實施揭露的內容,在預設時間範圍內,直接監測並獲得超聲回波訊號。示例性的,超聲回波訊號的一實例如第3圖所示,為了直觀說明,本申請人將超聲射頻資料用其包絡幅度的對數壓縮顯示(顯示範圍60 dB)。
可選的,在執行S202之前,還可以先產生超聲訊號資料集,包括: S203、在第一時刻t1獲取超聲回波資料; S204、將超聲回波資料載入傳播模式資料集中儲存的傳播模式資料,產生超聲回波訊號樣本; S205、產生超聲訊號資料集,其中,傳播模式資料包括介質樣本在振動激發下的各個時刻的傳播狀態資料。
具體的,根據t1時刻的超聲回波資料和傳播模式資料集,確定至少一超聲回波訊號樣本,根據確定的至少一超聲回波訊號樣本,確定超聲訊號資料集。
可選的,在第一時刻t1獲取超聲回波資料,可以為介質在未加振動激發時產生的超聲回波資料,也可以為介質在振動激發下產生的超聲回波資料。示例性的,未加振動激發時產生的超聲回波資料的一實例如第4圖所示。
可選的,傳播模式可以包括複數傳播狀態,即,傳播模式可以包括複數時刻對應的傳播狀態。由此,進一步可選的,傳播模式資料集可以儲存有在剪切波激發下介質的傳播模式,例如傳播模式資料集包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
具體的,是通過預先對仿體或者在體進行實際測試或者模擬,對不同給定黏彈性特性的組織施加特定的激發,將得到的後續不同時刻的系列超聲訊號呈現的傳播狀態所組合成的集合儲存為傳播模式資料集。
可選的,超聲回波訊號樣本可以根據傳播模式確定。示例性的,超聲回波訊號樣本可以與傳播模式具有一一對應的關係,或超聲回波訊號樣本可以對應於複數傳播模式,或超聲回波訊號樣本與傳播模式具有其他形式的對應關係。示例性的,超聲回波資料集中可以至少包括如第5圖至第8圖所示的超聲回波訊號樣本,其中,第5圖至第8圖所示的超聲回波訊號樣本分別由第4圖示意的超聲回波資料載入系統中儲存的傳播模式資料集產生。同樣,為了直觀說明,本申請人也將產生的超聲射頻資料用其包絡幅度的對數壓縮顯示(顯示範圍60 dB)
進一步的,S202中,將超聲回波訊號與超聲回波訊號樣本匹配,具體的,將第3圖所示的超聲回波訊號可以與第5圖至第8圖所示的超聲回波訊號樣本進行查找或比對。示例性的,在查找或比對時求超聲回波訊號和超聲回波訊號樣本的差的平方和,將計算結果最小值對應的超聲回波樣本,確定為與超聲回波訊號匹配的超聲回波訊號樣本,即確定第3圖所示的超聲回波訊號與第7圖所示的超聲回波訊號樣本匹配。但本發明並不限制進行比對的方法。
可選的,S201可以包括,在第二時刻t2獲取介質在振動激發下的超聲回波訊號,即將獲取超聲回波訊號的時刻記為t2。
可選的,S202可以包括,將第三時刻t3獲取的超聲回波訊號與超聲回波訊號樣本進行匹配,其中,第三時刻t3大於等於第二時刻t2。
其中,t3時刻可以與t2時刻重合,也可以晚於t2時刻。
可選的,S202之前,還可以包括: S206、將超聲回波訊號與超聲回波訊號樣本進行時間對準。
進一步可選的,S206可以包括: S2061、確定時間差,其中,時間差為第二時刻t2與振動激發起始時刻的差、及第二時刻t2與第一時刻t1的差中的較小值; S2062、根據時間差進行時間對準。
具體的,S2061中,當t1時刻無振動激發,即振動激發起始時刻晚於t1時刻時,第二時刻t2與振動激發起始時刻的差較小,此時時間差為第二時刻t2與振動激發起始時刻的差;當t1時刻有振動激發,即振動激發起始時刻早於t1時刻時,第二時刻t2與第一時刻t1的差較小,此時時間差為第二時刻t2與第一時刻t1的差。
S2062中,示例性的,記時間差為∆t,當t1時刻早於t2時刻時,可以將t3時刻獲取的超聲回波訊號與∆t時刻的超聲回波訊號樣本進行匹配,進一步的,此時t3時刻可以與t2時刻重合;當t1時刻晚於t2時刻時,可以將t3時刻獲取的超聲回波訊號與超聲回波訊號樣本進行匹配,進一步的,此時t3時刻晚於t2時刻的時間差為∆t。
可選的,S202之前,還可以包括,將超聲回波訊號與超聲回波訊號樣本進行空間對準,其中,空間對準可以包括:確定介質的邊界和深度,進行空間對準;或者基於預設的感興趣深度進行空間對準。
可選的,在S201中,可以從t2時刻起獲取一或複數超聲回波訊號,在S202中,可以將一或複數超聲回波訊號與超聲回波訊號樣本進行匹配。其中,將複數超聲回波訊號與超聲回波訊號樣本進行匹配後,選取匹配度最高的超聲回波訊號樣本作為用於確定介質黏彈性的超聲回波訊號樣本。相比於僅獲得一超聲回波訊號用於進行匹配,獲得複數超聲回波訊號並確定匹配度最高的超聲回波訊號樣本的方案能夠增強強健性。
本發明實施例於系統中僅預先儲存傳播模式資料集,包括一或者複數介質樣本的傳播模式資料。該些傳播模式資料是通過對一或者複數已知黏彈性資訊的仿體或者在體施加特定的激發,通過實際測量或者模擬後,得到的後續不同時刻的系列超聲回波訊號呈現的傳播狀態的集合。當需要測量時,以本次測量中某一時刻獲取的超聲回波資料載入傳播模式資料集中預先儲存的已知黏彈性資訊的介質樣本的傳播模式資料產生超聲訊號資料集。該超聲訊號資料集中包括一或複數介質樣本在當前激發下的對應的超聲回波樣本序列。將後續獲取的超聲回波訊號與超聲訊號資料集中的超聲訊號樣本進行匹配,通過簡單比對相似度最高的樣本,即可獲得介質黏彈性的測量結果,避免較為複雜的運動估計和特徵點選取,且需要儲存以及獲取的資訊較少,實現簡單,成本較低。
本發明實施例還揭露了一種介質黏彈性的測量裝置30,如第13圖所示,包括: 探頭301,用於獲取介質在振動激發下的檢測波回波訊號; 記憶體302,用於儲存資料集,資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本; 處理器303,用於將檢測波回波訊號與資料集中的檢測波回波訊號樣本匹配,確定介質的黏彈性資訊。
其中,探頭301可以用於執行S101或S201的步驟,處理器303可以用於執行S102或S202的步驟,記憶體302儲存的資料集,相關描述可以參考前述實施例中揭露的內容,此處不再贅述。
可選的,探頭301獲取的介質在振動激發下的檢測波回波訊號,可以為介質在剪切波激發下的超聲回波訊號,記憶體302儲存的資料集可以為超聲訊號資料集。
可選的,探頭301還可以用於,在第一時刻t1獲取超聲回波資料; 處理器303還可以用於,將超聲回波資料載入傳播模式資料集中儲存的傳播模式資料,產生超聲回波訊號樣本;產生超聲訊號資料集。其中,傳播模式資料包括已知黏彈性資訊的介質樣本在給定的振動激發下的各個時刻的傳播狀態資料。
可選的,探頭301還可以用於,從第二時刻t2獲取介質在振動激發下的超聲回波訊號; 處理器還可以用於,將從第二時刻t2獲取的超聲回波訊號中的一或複數與超聲回波訊號樣本進行匹配,其中,第二時刻t2在振動激發之後。
可選的,處理器303還可以用於:在匹配之前,將超聲回波訊號與超聲回波訊號樣本進行時間對準及空間對準,其中,時間對準用於確定第三時刻t3。
可選的,時間對準包括:確定時間差,時間差為第二時刻t2與振動激發起始時刻的差、及第二時刻t2與第一時刻t1的差中的較小值,根據時間差進行時間對準。
可選的,空間對準包括:確定介質的邊界和深度,進行空間對準;或,基於預設的深度進行空間對準。
可選的,傳播模式資料集可以包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
需要注意的是,第13圖所示的測量裝置,可以用於執行第1圖和第2圖所示的測量方法中的全部步驟,相關內容已在前述實施例中予以描述,此處不再贅述。
應當理解的是,本發明並不侷限於上面已經描述並在附圖中示出的流程及結構,並且可以在不脫離其範圍進行各種修改和改變。本發明的範圍僅由所附的申請專利範圍來限制。
30‧‧‧測量裝置301‧‧‧探頭302‧‧‧記憶體303‧‧‧處理器
此處的附圖被併入說明書中並構成本說明書的一部分,示出了符合本發明的實施例,並與說明書一起用於解釋本發明的原理。 第1圖是本發明實施例揭露的一種介質黏彈性的測量方法的流程圖; 第2圖是本發明實施例揭露的另一種介質黏彈性的測量方法的流程圖; 第3圖是本發明實施例揭露的一種超聲回波訊號的示意圖; 第4圖是本發明實施例揭露的一種超聲回波資料的示意圖; 第5圖是本發明實施例揭露的一種超聲回波訊號樣本的示意圖; 第6圖是本發明實施例揭露的一種超聲回波訊號樣本的示意圖; 第7圖是本發明實施例揭露的一種超聲回波訊號樣本的示意圖; 第8圖是本發明實施例揭露的一種超聲回波訊號樣本的示意圖; 第9圖是本發明實施例揭露的一種介質黏彈性的測量裝置的示意圖。
Claims (17)
- 一種介質黏彈性的測量方法,其特徵在於,包括:獲取一介質在一振動激發下的一檢測波回波訊號;將該檢測波回波訊號與一資料集中儲存的一檢測波回波訊號樣本匹配,確定該介質的黏彈性資訊;其中,該資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的檢測波回波訊號樣本。
- 如申請專利範圍第1項所述之方法,其中,該振動激發為一剪切波激發,該檢測波回波訊號為一超聲回波訊號,該資料集為一超聲訊號資料集。
- 如申請專利範圍第2項所述之方法,其中,該方法還包括:在一第一時刻t1獲取一超聲回波資料;將該超聲回波資料載入傳播模式資料集中儲存的一傳播模式資料,產生一超聲回波訊號樣本;產生該超聲訊號資料集;其中,該傳播模式資料包括已知黏彈性資訊的介質樣本在給定的振動激發下的各個時刻的傳播狀態資料。
- 如申請專利範圍第3項所述之方法,其中,該獲取該介質在振動激發下的該超聲回波訊號,包括:從一第二時刻t2獲取該介質在振動激發下的該超聲回波訊號;該將該檢測波回波訊號與資料集中儲存的檢測波回波訊號樣本匹配,包括:將從第二時刻t2獲取的超聲回波訊號中的一或者複數與該超聲回波訊號樣本進行一匹配;該第二時刻t2在振動激發之後。
- 如申請專利範圍第4項所述之方法,其中,該方法還包括:在該匹配之前,將該超聲回波訊號與該超聲回波訊號樣本進行一時間對準及一空間對準。
- 如申請專利範圍第5項所述之方法,其中,該時間對準包括:確定一時間差,該時間差為該第二時刻t2與振動激發起始時刻的差、及該第二時刻t2與該第一時刻t1的差中的較小值;根據該時間差進行時間對準;該空間對準包括:確定該介質的邊界和深度,進行空間對準;或,基於預設的深度進行空間對準。
- 如申請專利範圍第3項所述之方法,其中,該傳播模式資料集包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
- 如申請專利範圍第1項所述之方法,其中:該振動激發包括:連續或脈動方式的機械振動和聲輻射力脈動中的一項;該介質在該振動激發下的運動範圍為0.01微米至10毫米;該介質在該振動激發下的運動頻率為20赫茲至2500赫茲,持續時間為50微秒至1秒。
- 如申請專利範圍第1項所述之方法,其中:該檢測波回波訊號為一超聲回波訊號;該振動激發的脈動重複發生頻率為10至20000赫茲;該檢測波回波訊號包括至少連續採集的100訊框檢測波回波資料。
- 如申請專利範圍第1項所述之方法,其中,該黏彈性包括黏性和彈性中的至少一項; 當該黏彈性為黏性時,該黏彈性資訊包括一黏性資訊、或一黏性資訊和一彈性資訊;當該黏彈性為彈性時,該黏彈性資訊包括一彈性資訊、或一黏性資訊和一彈性資訊;當該黏彈性為黏性和彈性時,該黏彈性資訊包括一黏性資訊和一彈性資訊。
- 一種介質黏彈性的測量裝置,其特徵在於,包括:一探頭,用於獲取一介質在一振動激發下的一檢測波回波訊號;一記憶體,用於儲存一資料集,該資料集儲存有已知黏彈性資訊的介質樣本在振動激發下的一檢測波回波訊號樣本;一處理器,用於將該檢測波回波訊號與該資料集中的該檢測波回波訊號樣本匹配,確定該介質的黏彈性資訊。
- 如申請專利範圍第11項所述之測量裝置,其中,該振動激發為一剪切波激發,該檢測波回波訊號為一超聲回波訊號,該資料集為一超聲訊號資料集。
- 如申請專利範圍第12項所述之測量裝置,其中,該探頭還用於,在一第一時刻t1獲取一超聲回波資料;該處理器還用於,將該超聲回波資料載入一傳播模式資料集中儲存的一傳播模式資料,產生一超聲回波訊號樣本;產生該超聲訊號資料集;其中,該傳播模式資料包括已知黏彈性資訊的介質樣本在給定的振動激發下的各個時刻的傳播狀態資料。
- 如申請專利範圍第13項所述之測量裝置,其中,該探頭還用於,從一第二時刻t2獲取該介質在振動激發下的該超聲回波訊號; 該處理器還用於,將從第二時刻t2獲取的超聲回波訊號中的一或複數與該超聲回波訊號樣本進行匹配;其中,該第二時刻t2在振動激發之後。
- 如申請專利範圍第14項所述之測量裝置,其中,該處理器還用於:在該匹配之前,將該超聲回波訊號與該超聲回波訊號樣本進行一時間對準及一空間對準。
- 如申請專利範圍第15項所述之測量裝置,其中,該時間對準包括:確定一時間差,該時間差為該第二時刻t2與振動激發起始時刻的差、及該第二時刻t2與該第一時刻t1的差中的較小值;根據該時間差進行時間對準;該空間對準包括:確定該介質的邊界和深度,進行空間對準;或,基於預設的深度進行空間對準。
- 如申請專利範圍第13項所述之測量裝置,其中,該傳播模式資料集包括不同黏彈性的不同介質樣本分別在至少一種剪切波激發下,在一或者複數時刻產生的傳播狀態資料。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710650202.4A CN107495986A (zh) | 2017-07-21 | 2017-07-21 | 介质粘弹性的测量方法和装置 |
??201710650202.4 | 2017-07-21 | ||
CN201710650202.4 | 2017-07-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201907871A TW201907871A (zh) | 2019-03-01 |
TWI704906B true TWI704906B (zh) | 2020-09-21 |
Family
ID=60690305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW107121599A TWI704906B (zh) | 2017-07-21 | 2018-06-22 | 介質黏彈性的測量方法和裝置 |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN107495986A (zh) |
TW (1) | TWI704906B (zh) |
WO (1) | WO2019015396A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107495986A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性的测量方法和装置 |
CN110045371A (zh) * | 2019-04-28 | 2019-07-23 | 软通智慧科技有限公司 | 一种鉴定方法、装置、设备及存储介质 |
CN115856083A (zh) * | 2023-02-27 | 2023-03-28 | 中国汽车技术研究中心有限公司 | 汽车碰撞假人皮肤的性能测试方法、装置、设备和介质 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2877097A4 (en) * | 2012-07-27 | 2016-04-13 | Univ Leland Stanford Junior | MANIPULATION OF IMAGING SENSOR IN MEDICAL PROCEDURES |
CN103006274A (zh) * | 2012-12-21 | 2013-04-03 | 深圳大学 | 一种超声检测角膜粘弹性的方法和系统 |
US9726647B2 (en) * | 2015-03-17 | 2017-08-08 | Hemosonics, Llc | Determining mechanical properties via ultrasound-induced resonance |
CN104825195A (zh) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | 剪切波粘弹性成像方法及系统 |
CN107708575A (zh) * | 2015-06-01 | 2018-02-16 | 杜克大学 | 用于单个跟踪位置剪切波弹性成像的方法、系统和计算机程序产品 |
CN105266851B (zh) * | 2015-10-29 | 2018-01-05 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘度参数检测方法和设备 |
CN105740763A (zh) * | 2016-01-21 | 2016-07-06 | 珠海格力电器股份有限公司 | 身份的识别方法及装置 |
CN107505393A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性的测量方法和装置 |
CN107478723A (zh) * | 2017-07-21 | 2017-12-15 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性的测量方法和装置 |
CN107495986A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性的测量方法和装置 |
-
2017
- 2017-07-21 CN CN201710650202.4A patent/CN107495986A/zh active Pending
-
2018
- 2018-05-25 WO PCT/CN2018/088404 patent/WO2019015396A1/zh active Application Filing
- 2018-06-22 TW TW107121599A patent/TWI704906B/zh active
Also Published As
Publication number | Publication date |
---|---|
WO2019015396A1 (zh) | 2019-01-24 |
CN107495986A (zh) | 2017-12-22 |
TW201907871A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102309904B1 (ko) | 비등방성 매질을 이미지화하기 위한 전단파 엘라스토그라피 방법 및 기기 | |
AU2016347143B2 (en) | Elasticity detection method and device | |
McLaskey et al. | Hertzian impact: Experimental study of the force pulse and resulting stress waves | |
CN102641137B (zh) | 使用幅度-相位调制超声波的粘弹性测量 | |
TWI704906B (zh) | 介質黏彈性的測量方法和裝置 | |
Mariani et al. | Compensation for temperature-dependent phase and velocity of guided wave signals in baseline subtraction for structural health monitoring | |
Viola et al. | Sonorheometry: a noncontact method for the dynamic assessment of thrombosis | |
JP2010526626A (ja) | 関心領域の粘弾性の平均値を測定するための方法および装置 | |
JP2018531138A6 (ja) | 弾性の検出方法及び機器 | |
CN105266851A (zh) | 粘弹性介质的粘度参数检测方法和设备 | |
JP2013529098A (ja) | せん断波を使用する撮像方法および装置 | |
CN105662473A (zh) | 组织参数检测方法和系统 | |
US11719613B2 (en) | Method and device for quantifying viscoelasticity of a medium | |
Urban et al. | Generalized response of a sphere embedded in a viscoelastic medium excited by an ultrasonic radiation force | |
TWI704905B (zh) | 介質黏彈性的測量方法和裝置 | |
Long et al. | Clinical acceptance testing and scanner comparison of ultrasound shear wave elastography | |
Zhou et al. | Acoustic emission source localization using coupled piezoelectric film strain sensors | |
JP4899049B2 (ja) | 青果物の粘性測定方法、およびその装置 | |
JP7295126B2 (ja) | ハーモニックエラストグラフィにより誘導される超音波減衰パラメータを測定するための方法、本方法を実施するためのプローブおよびデバイス | |
Jensen et al. | Rapid measurements of intensities for safety assessment of advanced imaging sequences | |
Ostiguy et al. | In situ characterization technique to increase robustness of imaging approaches in structural health monitoring using guided waves | |
EP3505927A1 (en) | Method and system for measuring a characteristic loop sensitivity for an acoustic transducer | |
CN107505393A (zh) | 介质粘弹性的测量方法和装置 | |
RU2786332C2 (ru) | Способ для измерения параметра затухания ультразвука методом гармонической эластографии; зонд и устройство для осуществления способа | |
Mahat et al. | Natural frequency informed finite element modal analysis method for estimating elastic properties of solid materials |