ES2986863T3 - Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica; sonda y dispositivo para la puesta en práctica del procedimiento - Google Patents
Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica; sonda y dispositivo para la puesta en práctica del procedimiento Download PDFInfo
- Publication number
- ES2986863T3 ES2986863T3 ES19706655T ES19706655T ES2986863T3 ES 2986863 T3 ES2986863 T3 ES 2986863T3 ES 19706655 T ES19706655 T ES 19706655T ES 19706655 T ES19706655 T ES 19706655T ES 2986863 T3 ES2986863 T3 ES 2986863T3
- Authority
- ES
- Spain
- Prior art keywords
- ultrasonic
- probe
- acquisitions
- ultrasonic attenuation
- viscoelastic medium
- 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.)
- Active
Links
- 239000000523 sample Substances 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000005259 measurement Methods 0.000 title claims description 50
- 238000002091 elastography Methods 0.000 claims abstract description 57
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 description 20
- 230000006870 function Effects 0.000 description 13
- 238000002604 ultrasonography Methods 0.000 description 12
- 210000004185 liver Anatomy 0.000 description 8
- 230000033001 locomotion Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 4
- 238000002592 echocardiography Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000003187 abdominal effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000007863 steatosis Effects 0.000 description 2
- 231100000240 steatosis hepatitis Toxicity 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 208000004930 Fatty Liver Diseases 0.000 description 1
- 206010019668 Hepatic fibrosis Diseases 0.000 description 1
- 206010019708 Hepatic steatosis Diseases 0.000 description 1
- 206010020843 Hyperthermia Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000005520 electrodynamics Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 208000019423 liver disease Diseases 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
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—Clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
- A61B8/0833—Clinical applications involving detecting or locating foreign bodies or organic structures
- A61B8/085—Clinical applications 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/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- 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/54—Control of the diagnostic device
-
- 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
- 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52073—Production of cursor lines, markers or indicia by electronic means
-
- 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/52079—Constructional features
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Vascular Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
La invención se refiere a un procedimiento de medida de un parámetro de atenuación ultrasónica, guiado por elastografía armónica, que comprende una etapa de aplicación, mediante un vibrador contenido en una sonda en contacto con un medio viscoelástico, de una vibración continua de baja frecuencia, generando dicha vibración continua de baja frecuencia una onda elástica en el interior del medio viscoelástico, y generación (CW), durante la propagación de la onda elástica y mediante un transductor ultrasónico en contacto con el medio viscoelástico, de una serie de adquisiciones ultrasónicas, comprendiendo dicha serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, generándose dichos grupos de adquisiciones ultrasónicas a una velocidad de repetición, comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición; el parámetro de atenuación ultrasónica se mide a partir de las adquisiciones ultrasónicas obtenidas durante la aplicación de la vibración continua de baja frecuencia. (Traducción automática con Google Translate, sin valor legal)
Description
DESCRIPCIÓN
Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica; sonda y dispositivo para la puesta en práctica del procedimiento
Ámbito técnico
La invención pertenece al ámbito de la medición de la atenuación ultrasónica y más particularmente de la utilización de la elastografía armónica para guiar la determinación de la atenuación ultrasónica de un medio que presenta una señal ultrasónica tras la iluminación ultrasónica. Un primer objeto de la invención es un procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica. Un segundo objeto de la invención es una sonda de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica. Un tercer objeto de la invención es un dispositivo para la medición de la atenuación ultrasónica guiada por elastografía armónica. El procedimiento de medición de la atenuación ultrasónica guiado por elastografía armónica es particularmente adecuado para determinar las propiedades de un medio viscoelástico tal como el hígado humano o animal.
Estado de la técnica
Se sabe que la atenuación ultrasónica está correlacionada con la cantidad de grasa contenida en el hígado. Ésta por tanto puede ser utilizada para medir la cantidad de esteatosis en el hígado.
La solicitante ha desarrollado y comercializado un dispositivo que cuantifica de manera no invasiva la atenuación ultrasónica denominada CAP. La técnica utilizada es la medición de la atenuación ultrasónica guiada por elastografía transitoria controlada por vibración (o Vibration Controlled Transient Elastography, “VCTE”). La técnica VCTE se describe en el documento « T ransient Elastography: a new noninvasive method for assessment of hepatic fibrosis» de L. Sandrin y otros., publicado en Ultrasound in Medicine and Biology, vol. 29, páginas 1705-1713, 2003. El estudio de las señales ultrasónicas adquiridas para la medición de la elasticidad permite llegar a la atenuación ultrasónica del medio como se explica en los documentos:
• "The controlled attenuation parameter (CAP): A novel tool for the non-invasive évaluation of steatosis using Fibroscan®" de M. Sasso y otros, publicada en Clinical Research in Hepatology and Gastroenterology, 2011;
• "Controlled Atténuation Parameter (CAP): A Novel VCTE™ Guided Ultrasonic Atténuation Measurement for the Evaluation of Hepatic Steatosis: Preliminary Study and Validation in a Cohort of Patients with Chronic Liver Disease from Various Causes" de M. Sasso y otros, publicada en Ultrasound in Medecine and Biology, 2010;
• "Liver steatosis assessed by controlled attenuation parameter (cap) measured with the xl probe of the fibroscan: a pilot study assessing diagnostic accuracy" de M. Sasso y otros publicado en Ultrasound in Medecine and Biology, 2015
El dispositivo que pone en práctica esta técnica, denominado Fibroscan®, es capaz de medir la elasticidad y la atenuación ultrasónica del hígado humano de modo rápido, no invasivo y reproducible. En el FibroScan, la medición de la atenuación ultrasónica (CAP) es validada por la medición de la elasticidad: cuando la medición de la elasticidad es válida, el valor de CAP asociado se considera válido. Por tanto, no se trata de un guiado a priori sino de una validación a posteriori. Para medir con precisión la atenuación ultrasónica, el FibroScan valida la medición de CAP a partir del resultado de la elastografía por impulsos.
En tal dispositivo para elastografía transitoria, la onda de cizalladura por impulsos es generada por un vibrador colocado en contacto con el medio que haya que caracterizar. La propagación de la onda de cizalladura es seguida después con la ayuda de una serie de adquisiciones ultrasónicas realizadas por un transductor ultrasónico con una tasa de repetición elevada. Cada adquisición ultrasónica corresponde al menos a una emisión ultrasónica. Cada emisión ultrasónica puede ser asociada a la detección y registro sobre la marcha de los ecos generados por las partículas reflectantes presentes en el entorno estudiado para una gama de profundidades definida. Las señales ultrasónicas reflejadas son tratadas por correlación para llegar a los movimientos del tejido generados por la propagación de la onda de cizalladura en función del tiempo y de la posición en el medio. El estudio de estos movimientos permite llegar a la velocidad de propagación de la onda de cizalladura en el interior del medio viscoelástico y después a la elasticidad.
La medición del CAP (atenuación ultrasónica) por la técnica VCTE presenta varias limitaciones.
El primer inconveniente de la medición del CAP (atenuación ultrasónica) por la técnica VCTE es la dificultad de predecir que la sonda está efectivamente colocada enfrente del hígado. En efecto, es necesario realizar una medición válida por elastografía por impulsos para obtener una medición de la atenuación ultrasónica.
El segundo inconveniente de la medición del CAP por la técnica VCTE es el coste del dispositivo que requiere la utilización de vibración por impulsos, cuyo coste es elevado.
El tercer inconveniente de la medición del CAP por la técnica VCTE es la necesidad de realizar aproximadamente una decena de mediciones de elasticidad por elastografía por impulsos, cuya duración individual es de al menos 1 segundo, incluido el tiempo de cálculo, lo que conduce a una duración de examen de aproximadamente un minuto.
Hoy en día es posible utilizar los ultrasonidos para guiar el posicionamiento del vibrador para la elastografía transitoria. Es posible utilizar por ejemplo la formación de imágenes por ultrasonidos o una herramienta de orientación como la que se describe en la solicitud de patente. EP2739211 A1. Sin embargo, estas soluciones no son satisfactorias porque no permiten predecir directamente que la sonda esté bien situada con respecto al órgano estudiado. Además, estas técnicas dependen en gran medida del operador. Finalmente, las señales ultrasónicas pueden no ser muy específicas del hígado.
Existen además técnicas de elastografía denominadas armónicas. Estas técnicas se basan en la aplicación de una vibración continua que tenga una frecuencia comprendida entre 30 Hz y 100 Hz. Las ondas elásticas creadas en el interior del medio son ondas casi estacionarias, superposiciones de ondas de cizalladura y de ondas de compresión que no son tan eficientes como las técnicas de elastografía por impulsos para la medición precisa de las propiedades viscoelásticas pero que, sin embargo, utilizan ondas de cizalladura y pueden medir su propagación.
Finalmente, se sabe que las ondas de cizalladura de baja frecuencia se propagan bien en el hígado. Esta observación es cierta para las ondas de cizalladura por impulsos y armónicas.
Además, la técnica de elastografía armónica puede ser utilizada para guiar los procedimientos de tratamiento. Se trata, por ejemplo, de tratar tumores localizados por la técnica de elastografía armónica por medo de procedimientos de tipo hipertermia.
Problema técnico
Ninguna técnica de medición de la atenuación ultrasónica asegura un guiado óptimo y en tiempo real de la sonda ultrasónica utilizada con respecto al tejido que haya que caracterizar. En consecuencia, la medición puede resultar difícil de realizar habida cuenta de la imposibilidad de afirmar que el órgano estudiado está bien situado enfrente de la sonda. Esto no es favorable para una puesta en práctica con dispositivos pequeños y de utilización sencilla.
Resumen de la invención
Para resolver al menos parcialmente estos problemas, la presente invención describe una nueva técnica de medición de un parámetro de atenuación ultrasónica que es guiada por elastografía armónica.
Un objeto de la presente invención es un procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica que comprende una etapa de aplicación, con la ayuda de un vibrador comprendido en una sonda en contacto con un medio viscoelástico, de una vibración continua de baja frecuencia, generando la vibración continua de baja frecuencia una onda elástica en el interior del medio viscoelástico y de generación, durante la propagación de la onda elástica, con la ayuda de un transductor ultrasónico en contacto con el medio viscoelástico, de una serie de adquisiciones ultrasónicas, comprendiendo la citada serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, siendo generados los grupos de adquisiciones ultrasónicas con una tasa de repetición, comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición;
Siendo medido el parámetro de atenuación ultrasónica a partir de las adquisiciones ultrasónicas.
Según un modo de realización, las adquisiciones ultrasónicas son realizadas durante la aplicación de la vibración continua de baja frecuencia.
Se entiende por medición de atenuación ultrasónica guiada por elastografía armónica un procedimiento que comprende al menos una etapa de aplicación de una vibración continua y la medición de un parámetro que refleja la atenuación ultrasónica. En otras palabras, el procedimiento según la invención comprende a la vez una generación de una vibración continua, lo que es característico de una técnica de elastografía armónica, y la medición de un parámetro que refleja la atenuación ultrasónica.
Se entiende por atenuación ultrasónica cualquier parámetro que refleje la atenuación ultrasónica: Atenuación Ultrasónica en Banda Ancha (BUA, en dB/cm/MHz), la atenuación medida a una frecuencia concreta (en dB/cm), el Parámetro de Atenuación Controlada (CAP), etc.
Se entiende por vibración continua de baja frecuencia la reproducción continua de un motivo de forma de onda. Este motivo puede ser, por ejemplo, una sinusoide perfecta; se habla entonces de vibración monocromática. La vibración puede estar constituida igualmente por la reproducción de un motivo arbitrario. La vibración continua es interrumpida para detener el proceso de medición o cuando las condiciones de medición ya no sean satisfactorias. Las condiciones de medición pueden ser, por ejemplo, una condición sobre la fuerza de contacto con el medio estudiado. La frecuencia central de la vibración continua de baja frecuencia está comprendida entre 5 Hz y 500 Hz.
Se entiende por onda elástica la superposición de ondas de compresión y de ondas de cizalladura.
Se entiende por adquisición ultrasónica la emisión de un disparo ultrasónico. La citada emisión ultrasónica puede estar asociada a la detección y el registro sobre la marcha de los ecos generados por las partículas reflectantes presentes en una gama de profundidades definida del medio estudiado.
La serie de adquisiciones ultrasónicas está formada por tanto por una repetición de grupos de adquisiciones. Un grupo de adquisiciones comprende al menos una adquisición ultrasónica. Los grupos de adquisiciones son emitidos o generados con una primera tasa de repetición. La primera tasa de repetición es denominada igualmente tasa de repetición entre grupos. La primera tasa de repetición está comprendida típicamente entre 5 Hz y 500 Hz.
Cuando cada grupo de adquisiciones está formado por al menos dos adquisiciones ultrasónicas, las adquisiciones ultrasónicas que forman un mismo grupo son transmitidas o generadas con una tasa de repetición dentro del grupo comprendida típicamente entre 500 Hz y 100 kHz.
Ventajosamente, la utilización de una tasa de repetición baja durante la aplicación de la vibración continua permite medir los desplazamientos del tejido viscoelástico al tiempo que se limita la energía acústica enviada al propio tejido de manera que no se superen los límites de potencia acústica pico y media.
El término desplazamiento es considerado en el sentido amplio en este documento. El mismo engloba cualquier parámetro de movimiento como el desplazamiento, la velocidad, la deformación, la tasa de deformación, la velocidad de deformación y cualquier transformación matemática aplicada a estos parámetros.
Durante la aplicación de una vibración continua, se genera una onda elástica en el interior del medio viscoelástico.
La serie de adquisiciones ultrasónicas es utilizada para estudiar la propagación de la onda elástica en el interior del medio viscoelástico. Es posible detectar los ecos o señales ultrasónicas reflejadas por el medio viscoelástico y calcular, a partir de estas señales ultrasónicas reflejadas, los desplazamientos del medio viscoelástico provocados por la propagación de la onda elástica en el interior del medio viscoelástico generada por la vibración continua.
Por ejemplo, es posible calcular los desplazamientos del medio viscoelástico aplicando una técnica de correlación a las adquisiciones ultrasónicas que componen un mismo grupo de adquisiciones de la serie de adquisiciones ultrasónicas.
Es importante señalar que la propagación de la onda elástica evoluciona en función de la posición de la sonda en contacto con el medio estudiado.
Es entonces posible medir una propiedad de la propagación de la onda elástica en el interior del medio y calcular en tiempo real un indicador de posicionamiento a partir de las propiedades medidas. Idealmente, este indicador de posicionamiento en tiempo real es visualizado en tiempo real para guiar al operador con el fin de que el mismo encuentre la mejor posición para la sonda. Ejemplos de propiedades medidas para calcular el indicador de posicionamiento son la amplitud y la fase de la onda elástica, medidas en función de la profundidad en el tejido que haya que caracterizar. Es igualmente posible calcular la velocidad de fase de la onda elástica.
En lo que sigue de este documento, «indicador de posicionamiento en tiempo real» e «indicador de posicionamiento» se refieren al mismo indicador de posicionamiento en tiempo real.
Se entiende por tiempo real un indicador cuya visualización es actualizada regularmente durante el examen. En general, la tasa de actualización es de aproximadamente 20 Hz, pero también puede ser del orden de 1 Hz.
En otras palabras, el indicador de posicionamiento da una probabilidad de la presencia del medio viscoelástico que haya que caracterizar con respecto al transductor ultrasónico.
Es importante señalar que la vibración continua es utilizada ante todo para verificar la posición de la sonda utilizada para la medición de la atenuación ultrasónica. Por ejemplo, la vibración continua puede ser utilizada para comprobar la presencia del parénquima hepático delante de la sonda. En otras palabras, durante la etapa de aplicación de una vibración continua, la medición indirecta de las propiedades viscoelásticas del medio es posible pero no indispensable. Sin embargo, a partir de la velocidad de fase de la onda elástica se puede deducir efectivamente un valor de elasticidad. En efecto, la medición de los desplazamientos generados en el medio viscoelástico por la citada propagación permite llegar a continuación a la velocidad de propagación de la onda elástica. Si se hace la hipótesis de que la onda elástica es principalmente una onda de cizalladura, se puede calcular entonces la elasticidad del medio utilizando la fórmula E=3pVs2 donde E es la elasticidad o módulo de Young, p la densidad y Vs la velocidad de la onda elástica.
El procedimiento de guiado de una medición de atenuación ultrasónica por elastografía armónica según la invención permite por tanto predecir si el posicionamiento de la sonda es favorable con la ayuda de una técnica de elastografía armónica. Se trata de un guiado. Las mediciones de atenuación ultrasónica realizadas en el transcurso del tiempo pueden ser almacenadas en una memoria. En particular, cada medición de atenuación ultrasónica puede ser asociada a un coeficiente de calidad. En general, el coeficiente de calidad será un coeficiente comprendido entre 0 y 1. El valor 0 corresponde a una calidad mala y el valor 1 a una calidad de medición muy buena.
El cálculo del coeficiente de calidad puede ser efectuado a partir de las propiedades de propagación de la onda elástica. El mismo puede igualmente incluir características de las señales ultrasónicas y criterios de calidad que emanan del cálculo de la atenuación ultrasónica.
La medición de atenuación ultrasónica final retenida puede ser calculada a partir de las mediciones de atenuación ultrasónica y de los coeficientes de calidad almacenados en la memoria.
En otras palabras, el guiado por elastografía armónica permite asociar un coeficiente de calidad a cada medición de atenuación ultrasónica y guiar el posicionamiento de la sonda con respecto al tejido que haya que caracterizar al proporcionar al operador un indicador de posicionamiento predictivo de la presencia del órgano estudiado.
Por ejemplo, la medición de atenuación ultrasónica retenida es la media de las mediciones de atenuación ultrasónica almacenadas ponderadas por sus coeficientes de calidad. El número total de mediciones de atenuación ultrasónica puede estar comprendido típicamente entre 1 y varios miles de mediciones. Por ejemplo, si las mediciones son acumuladas durante 60 segundos a una velocidad de 20 mediciones por segundo, N será igual a 1200.
Ventajosamente, el procedimiento de guiado de una medición de la atenuación ultrasónica por elastografía armónica según la invención permite realizar una medición de atenuación ultrasónica del tejido que haya que caracterizar de modo fiable y reproducible al asegurar un posicionamiento óptimo de la sonda de modo simple y preciso gracias a una técnica de elastografía armónica.
Ventajosamente, el procedimiento de guiado de una medición de atenuación ultrasónica por elastografía armónica según la invención permite reducir los costes de fabricación por la utilización de un sistema vibratorio más sencillo que en elastografía por impulsos.
Ventajosamente, el procedimiento de guiado de una medición de atenuación ultrasónica por elastografía armónica según la invención permite proporcionar al operador un indicador de posicionamiento en tiempo real.
El procedimiento de guiado de una medición de atenuación ultrasónica por elastografía armónica según la invención puede igualmente presentar una o varias de las características siguientes, consideradas individualmente o según todas las combinaciones técnicamente posibles:
- el procedimiento según la invención comprende además una etapa de determinación, a partir de la serie de adquisiciones ultrasónicas, de al menos una propiedad de la propagación de la onda elástica en el interior del medio viscoelástico;
- la propiedad de propagación de la onda elástica en el interior del medio viscoelástico es utilizada para calcular un indicador de posicionamiento en tiempo real de la sonda con respecto al medio viscoelástico que haya que estudiar;
- el procedimiento según la invención comprende además una etapa de visualización en tiempo real del indicador de posicionamiento en tiempo real;
- la etapa de aplicación de una vibración continua de baja frecuencia es activada únicamente si la fuerza de contacto entre el vibrador y el medio viscoelástico es superior a un umbral inferior predeterminado;
- la etapa de aplicación de una vibración continua de baja frecuencia es activada únicamente si la fuerza de contacto entre el vibrador y el medio viscoelástico es superior a un umbral inferior e inferior a un umbral superior predeterminado;
- la serie de adquisiciones ultrasónicas está formada por una repetición de grupos que comprenden al menos dos adquisiciones ultrasónicas que tienen una tasa de repetición dentro del grupo comprendida entre 500 Hz y 10 kHz y una tasa de repetición comprendida entre 10 Hz y 10 kHz;
- la tasa de repetición es menor que la frecuencia de vibración continua;
- se calcula y se visualiza un parámetro de la onda ultrasónica que refleja la atenuación ultrasónica;
- el parámetro de atenuación ultrasónica es un parámetro instantáneo y se calcula un coeficiente de calidad asociado con la medición del parámetro de atenuación ultrasónica instantáneo a partir de una propiedad de la propagación de la onda elástica en el interior del medio viscoelástico y/o de las propiedades de la onda ultrasónica; Se entiende por parámetro instantáneo una medición obtenida en cada adquisición ultrasónica o a partir de una sola adquisición ultrasónica;
- se calcula un parámetro de atenuación ultrasónica medio a partir de una pluralidad de parámetros de atenuación ultrasónica instantáneos y de los coeficientes de calidad asociados con cada parámetro de atenuación ultrasónica instantáneo;
- los parámetros de atenuación ultrasónica, asociados o no a coeficientes de calidad calculados sobre la base de las propiedades de la propagación de la onda elástica en el interior del medio viscoelástico, son almacenados en una memoria;
- el umbral inferior predeterminado de fuerza de contacto para la aplicación de la vibración continua es elegido típicamente igual a 1 N;
- la frecuencia de la vibración de baja frecuencia, cSWF, aplicada por el vibrador está comprendida entre 5 Hz y 500 Hz;
- la amplitud de la vibración de baja frecuencia aplicada por el vibrador está comprendida entre 10 pm y 5 mm; - la serie de adquisiciones ultrasónicas está formada por una repetición de grupos que comprenden al menos dos adquisiciones ultrasónicas que tienen una tasa de repetición dentro del grupo comprendida entre 500 Hz y 10 kHz y una primera tasa de repetición comprendida entre 10 Hz y 10 kHz;
- la primera tasa de repetición es menor que la frecuencia de vibración continua;
- el parámetro de atenuación ultrasónica es elegido entre
o la BUA
o la atenuación de una frecuencia ultrasónica particular;
o parámetro de atenuación ultrasónica controlada (CAP);
- la propiedad de la propagación de la onda elástica es derivada al operador;
- la propiedad de propagación de la onda elástica es elegida entre:
o la amplitud de la onda elástica;
o la fase de la onda elástica;
o el módulo de Young del medio viscoelástico;
o el módulo de cizalladura del medio viscoelástico;
o la velocidad de cizalladura del medio viscoelástico.
Otro objeto de la presente invención es una sonda para la puesta en práctica del procedimiento de elastografía híbrida según la invención. La sonda según la invención comprende:
- un vibrador configurado para aplicar una vibración continua de baja frecuencia al medio viscoelástico, generando la vibración continua de baja frecuencia una onda elástica en el interior del medio viscoelástico; - un transductor ultrasónico configurado para emitir una serie de adquisiciones ultrasónicas, comprendiendo la citada serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, siendo generados los grupos de adquisiciones ultrasónicas con una tasa de repetición, comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición;
- medios de cálculo y de visualización en tiempo real de un indicador de posicionamiento de la sonda, siendo calculado el citado indicador de posicionamiento a partir de una propiedad de la propagación de la onda elástica, siendo determinada la citada propiedad de propagación de la onda elástica a partir de la serie de adquisiciones ultrasónicas;
siendo medido el parámetro de atenuación ultrasónica a partir de las adquisiciones ultrasónicas realizadas durante la aplicación de la vibración continua de baja frecuencia.
Según un modo de realización, la citada sonda está configurada además para aplicar la vibración continua cuando la fuerza de contacto entre la sonda y el medio viscoelástico es superior a un valor predeterminado.
La sonda según la invención permite la puesta en práctica del procedimiento según la invención.
El transductor ultrasónico es utilizado para enviar la serie de adquisiciones ultrasónicas al interior del medio viscoelástico. El mismo transductor ultrasónico detecta las señales ultrasónicas reflejadas en cada adquisición. Las señales ultrasónicas son tratadas a continuación para detectar los desplazamientos del medio viscoelástico generados por las ondas elásticas.
Se entiende por medios de cálculo al menos un microprocesador y una memoria destinada a almacenar las adquisiciones ultrasónicas y los resultados de los cálculos, tal como un indicador de posicionamiento de la sonda o una propiedad de propagación de la onda elástica.
Se entiende por medio de visualización una pantalla o un indicador configurado para visualizar el indicador de posicionamiento. El indicador puede ser, por ejemplo, un indicador luminoso tal como un diodo o un indicador sonoro. La sonda de elastografía híbrida según la invención puede presentar igualmente una o varias de las características siguientes, consideradas individualmente o según todas las combinaciones técnicamente posibles:
- el vibrador es un motor eléctrico o una bobina de audio o un actuador electrodinámico;
- el transductor ultrasónico está montado en el eje del vibrador;
- la sonda según la invención comprende además medios para activar la acumulación de mediciones;
- el transductor ultrasónico es circular con un diámetro comprendido entre 2 mm y 15 mm;
- el transductor ultrasónico tiene una frecuencia de funcionamiento comprendida entre 1 MHz y 15 MHz;
- la sonda está configurada para calcular un parámetro de atenuación ultrasónica medio, siendo calculado el citado parámetro de atenuación ultrasónica medio a partir de una multiplicidad de parámetros de atenuación ultrasónica instantáneos y de coeficientes de calidad asociados a los parámetros de atenuación ultrasónica instantáneos;
- el transductor ultrasónico es una sonda abdominal convexa.
Otro objeto de la presente invención es un dispositivo para elastografía híbrida que ponga en práctica el procedimiento de elastografía híbrida según la invención.
Tal dispositivo híbrido según la invención comprende:
- una sonda según la invención;
- una unidad central conectada a la sonda y que comprende al menos medios de cálculo para el tratamiento de las señales ultrasónicas reflejadas, medios de visualización y medios de control y/o de entrada.
En una puesta en práctica particular, la unidad central está colocada en el interior de la sonda.
Lista de las figuras
Otras características y ventajas de la invención se desprenderán claramente de la descripción que de la misma se hace a continuación, a título indicativo y en modo alguno limitativo, con referencia a las figuras entre las cuales: - la figura 1 ilustra las etapas del procedimiento de elastografía híbrida según la invención;
- la figura 2 muestra esquemáticamente las vibraciones aplicadas por el vibrador y las adquisiciones ultrasónicas durante la puesta en práctica del procedimiento según la invención ilustrado en la figura 1;
- la figura 3 muestra esquemáticamente un modo de realización particular del procedimiento de elastografía ilustrado en la figura 1, denominado modo estroboscópico;
- la figura 4 muestra los resultados obtenidos al poner en práctica la parte del procedimiento según la invención relativa al posicionamiento del vibrador;
- la figura 5 representa una sonda para elastografía híbrida según la invención;
- la figura 6a representa un modo de realización particular de la sonda para elastografía híbrida según la invención;
- la figura 6b representa un dispositivo para elastografía híbrida según la invención.
Descripción detallada
La figura 1 ilustra las etapas del procedimiento de elastografía híbrida P según la invención.
El procedimiento P según la invención comprende una etapa CW de aplicación de una vibración continua de baja frecuencia con la ayuda de un vibrador comprendido en una sonda en contacto con el medio viscoelástico.
La frecuencia central de la vibración continua está comprendida entre 5 Hz y 500 Hz.
La etapa CW del procedimiento P comprende además la generación, por el transductor ultrasónico, de una serie de adquisiciones ultrasónicas. La serie de adquisiciones ultrasónicas comprende grupos de adquisiciones ultrasónicas. Los grupos de adquisiciones ultrasónicas son emitidos con una tasa de repetición comprendida entre 5 Hz y 500 Hz, comprendiendo cada grupo al menos una adquisición ultrasónica.
La tasa de repetición de los grupos ultrasónicos es denominada igualmente tasa de repetición entre grupos.
Una adquisición ultrasónica comprende la emisión de un disparo ultrasónico seguida de la detección y el registro de las señales ultrasónicas reflejadas o ecos.
La aplicación de una vibración continua al medio viscoelástico genera una onda elástica en el interior del citado medio. La onda elástica comprende una superposición de ondas de cizalladura y de ondas de compresión. El estudio de las propiedades de esta onda elástica permite obtener informaciones concernientes al correcto posicionamiento de la sonda con respecto al medio viscoelástico.
El medio viscoelástico que haya que caracterizar difunde al menos parcialmente los disparos ultrasónicos. Por tanto, es posible detectar las señales ultrasónicas reflejadas durante la emisión de la primera serie de adquisiciones ultrasónicas.
La detección de señales ultrasónicas reflejadas puede ser realizada con la ayuda del mismo transductor ultrasónico utilizado para la transmisión.
A partir de las señales ultrasónicas reflejadas se puede determinar un parámetro de atenuación ultrasónica. Por ejemplo, se puede determinar el valor CAP_I de la atenuación ultrasónica correspondiente a una adquisición ultrasónica determinada. El valor CAP_I es denominado igualmente valor individual o instantáneo de atenuación ultrasónica o parámetro instantáneo de atenuación ultrasónica.
Las señales ultrasónicas reflejadas detectadas durante la etapa CW son tratadas sucesivamente durante una etapa de determinación de al menos una propiedad de la propagación de la onda elástica en el interior del medio viscoelástico CW_P.
Típicamente, durante esta etapa, las señales ultrasónicas reflejadas se correlacionan entre sí de modo que miden los desplazamientos del medio viscoelástico generados por la onda elástica generada por la aplicación de la vibración continua, según una técnica conocida en el ámbito de la elastografía y más generalmente de los ultrasonidos.
A partir de los desplazamientos medidos en el interior del medio viscoelástico, es posible calcular propiedades de la onda elástica tales como la amplitud y la fase en función de la posición en el interior del medio viscoelástico. La posición de un punto en el interior del medio viscoelástico se mide como la distancia entre el transductor ultrasónico y el citado punto calculada a lo largo de la dirección de propagación de los ultrasonidos emitidos por el transductor. Por esta razón, la posición de un punto en el interior del medio viscoelástico es denominada generalmente profundidad.
Es igualmente posible determinar otros parámetros de la onda elástica en el interior del medio viscoelástico, tales como la velocidad de fase o la atenuación de la onda elástica.
Se pueden calcular las variaciones de la amplitud y de la fase de la onda elástica en función de la profundidad en el interior del tejido. Al realizar un ajuste entre el modelo teórico y las propiedades medidas, es posible extraer un parámetro de calidad del ajuste. A partir de este parámetro de calidad del ajuste y/o de las otras propiedades de la onda elástica, es posible calcular un indicador de posicionamiento RT_IP de la sonda con respecto al tejido que haya que caracterizar.
Por ejemplo, uno de los modelos teóricos utilizados prevé una variación lineal del retardo de fase a la frecuencia central de la onda elástica con la profundidad en el medio que haya que caracterizar. En este caso el ajuste es lineal y el parámetro de calidad del ajuste traduce la linealidad de la fase en función de la profundidad en el medio. Un posible indicador es el coeficiente de determinación R2 que da la calidad de la predicción de la regresión lineal de la curva del retardo de fase en función de la profundidad en la gama de profundidad estudiada.
Según un modo de realización, la etapa CW_P de determinación de al menos una propiedad de la onda elástica en el interior del tejido es efectuada al mismo tiempo que la etapa de aplicación de la vibración continua CW y de detección de las primeras señales ultrasónicas reflejadas.
Gracias al procedimiento P según la invención es posible por tanto medir en tiempo real las propiedades de la onda elástica en el interior del tejido y obtener en tiempo real el indicador de posicionamiento de la sonda RT_IP.
Ventajosamente, una tasa de repetición baja permite reducir el tamaño de los datos registrados durante la etapa de generación de la serie de adquisiciones ultrasónicas en CW y tratar estos datos en tiempo real para obtener el indicador de posicionamiento RT_IP. El valor de este indicador está comprendido típicamente entre 0 y 1. El valor 0 corresponde a un indicador malo y el valor 1 a un indicador bueno.
Ventajosamente, el indicador de posicionamiento es facilitado al operador para ayudarle a encontrar un punto de medición satisfactorio. Éste puede ser facilitado, por ejemplo (lista no exhaustiva), en forma de una visualización de un indicador coloreado, en forma de una barra más o menos larga, etc.
A partir de la señal ultrasónica se calcula igualmente un coeficiente de calidad de la medición de la atenuación ultrasónica, CAP_C. El valor de este está comprendido típicamente entre 0 y 1. El valor 0 corresponde a una calidad baja y el valor 1 a una calidad elevada.
El coeficiente CAP_C está asociado al valor individual de la atenuación ultrasónica CAP_I obtenido a partir de los datos ultrasónicos en curso de adquisición.
El coeficiente de calidad CAP_C puede calcularse, por ejemplo, únicamente a partir de las propiedades de la señal ultrasónica. Éste puede ser igualmente una combinación de la calidad de la señal ultrasónica y de propiedades de la onda elástica.
Según un modo de realización, el valor individual de la atenuación ultrasónica sólo se conserva si el indicador de posicionamiento RT_IP es bueno. En el caso en que el indicador de posicionamiento RT_IP sea malo, el coeficiente CAP_C correspondiente es por ejemplo puesto a cero.
Según un modo de realización, la vibración continua sólo se activa si la fuerza de contacto F entre el vibrador y el tejido viscoelástico es superior a un umbral inferior predeterminado. El valor del umbral es típicamente de 1 N. Ventajosamente, este umbral inferior asegura un acoplamiento suficiente entre la sonda y el medio viscoelástico. Según un modo de realización, la vibración continua sólo se activa si la fuerza de contacto F entre el vibrador y el tejido viscoelástico es inferior a un umbral superior predeterminado. El valor del umbral es típicamente de 10 N. Ventajosamente, este umbral superior asegura que la vibración no se deforme y que el medio estudiado no se dañe. A causa del movimiento vibratorio continuo del vibrador, la determinación de la fuerza F de contacto entre el vibrador y el medio es más compleja que en el caso de un procedimiento de elastografía transitoria estándar. En presencia de la vibración continua de baja frecuencia, la fuerza de contacto entre el vibrador y el medio viscoelástico viene dada por la fórmula siguiente:
F = k(x A x cos(2 nfbaJat))
En esta fórmula x es el desplazamiento del vibrador, k la constante elástica del muelle colocado en la sonda, A la amplitud de la vibración continua y fbaja la frecuencia de la vibración continúa.
La fuerza F puede ser medida con la ayuda de un sensor de fuerza colocado en la sonda para elastografía híbrida. Sucesivamente, aplicando un filtro de paso bajo a la señal así medida, es posible eliminar la parte de baja frecuencia y deducir la fuerza media de contacto:
F media fc(x)
Ventajosamente, el valor de la fuerza media aplicada es facilitado al operador para que lo adapte de modo que continúen la vibración de baja frecuencia y la adquisición de datos.
Ventajosamente, los valores individuales CAP_I son acumulados en una memoria y utilizados para calcular un valor medio CAP_M. CAP_M puede ser calculado de varias formas. Por ejemplo
£ f=1(CAP_I(i) x CAP_C(i))
CAP_M =
Z f=1(CAP_C(i))
Los valores CAP_C son utilizados entonces para ponderar los valores CAP_I individuales medidos. El valor CAP_M se mantiene al final del examen como el valor de atenuación ultrasónica medido. La unidad del valor CAP_M es, por ejemplo, el dB/m.
La figura 2 muestra esquemáticamente:
- la vibración continua de baja frecuencia cSW aplicada por el vibrador durante la etapa CW ilustrada en la figura 1;
- la serie de adquisiciones ultrasónicas PA formada por grupos G de adquisiciones y generadas por el transductor ultrasónico durante la etapa CW ilustrada en la figura 1.
Durante la etapa de aplicación de una vibración CW continua, el vibrador oscila a una frecuencia comprendida entre 5 Hz y 500 Hz, con una amplitud comprendida entre 10 pm y 5 mm.
Ventajosamente, gracias a la baja amplitud y a la baja frecuencia de la vibración continua, un operador puede mantener fácilmente la sonda en contacto con el medio viscoelástico.
Al mismo tiempo que la aplicación de una vibración continua de baja frecuencia, el transductor ultrasónico emite adquisiciones ultrasónicas PA formadas por grupos G de adquisiciones ultrasónicas. En el ejemplo ilustrado en la figura 2, cada grupo G comprende dos adquisiciones ultrasónicas.
Los grupos G de adquisiciones ultrasónicas son emitidos con una tasa de repetición LPRF comprendida entre 10 Hz y 500 Hz o tasa de repetición entre grupos o simplemente tasa de repetición. Las adquisiciones ultrasónicas pertenecientes al mismo grupo G son emitidas con una tasa de repetición dentro del grupo HPRF comprendida entre 500 Hz y 10 kHz.
El transductor ultrasónico detecta igualmente las señales ultrasónicas reflejadas durante la generación de las adquisiciones ultrasónicas PA, como se explica con referencia a la etapa CW ilustrada en la figura 1. A partir de la primera serie de adquisiciones ultrasónicas PA, es posible calcular, por una correlación CORR entre señales ultrasónicas pertenecientes a un mismo grupo G, los desplazamientos del medio viscoelástico. Los citados desplazamientos del medio viscoelástico son generados por la propagación de la onda elástica generada por la vibración continua aplicada por el vibrador.
Es importante señalar que existe una gran cantidad de secuencias ultrasónicas posibles para la puesta en práctica este procedimiento y que los elementos indicados no constituyen en modo alguno una lista exhaustiva del ámbito de las posibles.
Ventajosamente, aplicando una técnica de correlación a las adquisiciones ultrasónicas que pertenecen al mismo grupo G - y por tanto próximas temporalmente - es posible detectar pequeños desplazamientos del orden de 1 pm a 10 pm.
Como se explica con referencia a la etapa CW_P ilustrada en la figura 1, los desplazamientos del medio viscoelástico son utilizados después para calcular propiedades de la onda elástica tales como su amplitud y su fase en función de la profundidad en el medio. Comparando las propiedades medidas con un modelo teórico es posible deducir en tiempo real un indicador de posicionamiento RT_IP.
Por ejemplo, el indicador de posicionamiento puede estar vinculado a la linealidad de la fase de la onda elástica en función de la profundidad en el medio que haya que caracterizar. El indicador depende entonces de la calidad del ajuste de la evolución de la fase en función de la profundidad por una línea recta.
Por ejemplo, el indicador de posicionamiento puede estar relacionado con la disminución de la amplitud de la onda elástica en función de la profundidad en el medio que haya que caracterizar. El indicador depende entonces de la calidad del ajuste en 1/Zn donde Z es la profundidad y n un coeficiente entero comprendido entre 1 y 3.
Por ejemplo, el valor del indicador de posicionamiento en tiempo real RT_IP está comprendido entre 0 y 1, con valores próximos a 1 si la sonda está bien posicionada con respecto al medio viscoelástico de interés.
La figura 3 muestra un modo de realización particular de las etapas CW y CW_P del procedimiento P según la invención, denominado modo estroboscópico.
La línea sinusoidal en trazo continuo representa esquemáticamente la vibración continua cSW aplicada por el primer vibrador. La vibración continua cSW tiene por ejemplo una frecuencia central cSWF de 50 Hz correspondiente a un periodo de 20 ms.
Las líneas verticales continuas representan los grupos G de adquisiciones ultrasónicas que forman la primera serie de adquisiciones ultrasónicas PA. Los grupos G son emitidos con una primera tasa de repetición LPRF. Según el modo de adquisición estroboscópica, la primera tasa de repetición LPRF es menor que la frecuencia central de la vibración continua cSWF.
La tasa de repetición dentro del grupo está comprendida entre 500 Hz y 100 kHz, lo que permite medir pequeños desplazamientos del orden de 1 pm.
Los círculos blancos y las flechas a lo largo de la vibración continua cSW corresponden a los muestreos realizados por cada grupo G de adquisiciones ultrasónicas.
Gracias al hecho de que la tasa de repetición LPRF de los grupos G es inferior a la frecuencia central de la vibración continua cSW, es posible muestrear de modo completo la vibración continua cSW después de algunos períodos de oscilaciones, como lo muestran los círculos blancos.
Ventajosamente, el modo estroboscópico permite muestrear de modo completo la vibración continua cSW utilizando una primera tasa de repetición LPRF baja. La utilización de una tasa de repetición baja permite tratar las señales reflejadas en tiempo real y por tanto obtener el indicador de posicionamiento R<t>_IP en tiempo real.
La figura 4 muestra esquemáticamente los resultados obtenidos al poner en práctica la parte del procedimiento P según la invención relativa al posicionamiento del vibrador.
El gráfico CW_Disp muestra el desplazamiento (o cualquier otro parámetro de movimiento como la velocidad, la deformación, la tasa de deformación) del medio viscoelástico en una región de interés ROI en función de la profundidad Z en el medio y del tiempo T. Los desplazamientos se representan utilizando una escala de colores falsos, representando los colores menos oscuros un desplazamiento según la dirección positiva del eje D. Los desplazamientos son causados por la vibración continua de baja frecuencia aplicada por el vibrador y son medidos por el transductor ultrasónico UT colocado en contacto con la superficie del medio, en Z=0.
A partir de los desplazamientos medidos CW_Disp en la región de interés ROI en el interior del medio viscoelástico, es posible extraer en tiempo real informaciones RT _Info concernientes a la onda elástica que se propaga en el interior del medio y generada por la vibración continua. Ejemplos de tales propiedades son la amplitud A y la fase Ph de la onda elástica en función de la profundidad en el interior del medio.
Comparando los valores de A y Ph medidos con umbrales predeterminados es posible determinar un indicador de posicionamiento del vibrador con respecto al medio viscoelástico.
Alternativamente, es posible obtener un parámetro de calidad del ajuste AJ entre las cantidades medidas A y Ph y un modelo teórico que describe la amplitud y fase de una onda elástica que se propaga en el interior del medio. En este caso, el indicador de posicionamiento es obtenido a partir del parámetro de calidad del ajuste AJ. Por ejemplo, un parámetro de calidad del ajuste es el coeficiente de determinación R2 que da la calidad de la predicción de la regresión lineal de la curva del retardo de fase en función de la profundidad en la gama de profundidad estudiada.
Según un modo de realización, el parámetro de calidad del ajuste AJ está comprendido entre 0 y 1.
Una vez calculado, el indicador de posicionamiento se puede mostrar en forma de un número o de una letra o utilizando una escala de colores. Alternativamente, el indicador de posicionamiento puede ser una simple indicación visual de tipo disco coloreado. Alternativamente, el indicador de posicionamiento puede ser una simple indicación visual de tipo «Posicionamiento correcto» que indica que el operador puede activar la etapa de elastografía transitoria.
Según un modo de realización, la velocidad de propagación de la onda elástica se toma como una medición de la elasticidad del medio.
Durante la puesta en práctica del procedimiento P según la invención, se calculan y muestran simultáneamente los gráficos CW_Disp, RT_Info y el indicador de posicionamiento del vibrador.
Ventajosamente, gracias a la estructura de la serie de adquisiciones ultrasónicas, el indicador de posicionamiento RT_IP así como el gráfico RT_lnfo pueden ser calculados y visualizados en tiempo real.
La figura 5 representa esquemáticamente una sonda PR para la medición de la atenuación ultrasónica guiada por elastografía armónica.
La sonda PR comprende:
- un vibrador VIB configurado para aplicar una vibración continua de baja frecuencia al medio viscoelástico, generando la vibración continua de baja frecuencia una onda elástica en el interior del medio viscoelástico;
- un transductor ultrasónico TUS configurado para emitir una serie de adquisiciones ultrasónicas, comprendiendo la citada serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, siendo generados los grupos de adquisiciones ultrasónicas con una tasa de repetición, comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición;
Según el modo de realización ilustrado en la figura 6, el transductor ultrasónico TUS está montado en el eje del vibrador VIB.
Según un modo de realización, la sonda PR comprende medios de cálculo para calcular en tiempo real el indicador de posicionamiento RT_IP a partir de las adquisiciones ultrasónicas.
Según un modo de realización, la sonda PR comprende medios de cálculo y de visualización I del indicador de posicionamiento en tiempo real RT_IP.
Según un modo de realización, la tasa de la actualización de la visualización del indicador de posicionamiento es superior a 5 Hz.
Ventajosamente, la visualización del indicador de posicionamiento en tiempo real a nivel de la sonda permite al operador optimizar el posicionamiento de la sonda sin apartar la vista de la sonda y del cuerpo del paciente. Esto simplifica la operación de posicionamiento de la sonda.
Según un modo de realización, el transductor ultrasónico TUS puede ser fijado al cuerpo de la sonda con la ayuda de una punta PT.
El vibrador VIB hace oscilar la sonda PR. Durante esta oscilación, el transductor ultrasónico TUS es empujado contra el medio viscoelástico aplicando la vibración continua de baja frecuencia y generando la onda elástica en el interior del medio.
Según un modo de realización, el vibrador VIB para la aplicación de la vibración de baja frecuencia comprende un anillo vibratorio colocado alrededor del transductor ultrasónico TUS o alrededor de la punta de la sonda PT.
Según un modo de realización, la punta de la sonda PT es móvil y puede ser accionada por el vibrador VIB. El transductor ultrasónico TUS es empujado entonces contra el medio viscoelástico para aplicar la vibración, según la dirección de la flecha en la figura 6.
Según un segundo modo de realización, ilustrado en la figura 7a, la sonda PR es una sonda inercial sin partes móviles. En este caso, el movimiento del vibrador VIB dentro de la sonda PR provoca el movimiento de la sonda y la vibración continua o por impulsos es aplicada nuevamente empujando el transductor TUS contra el medio viscoelástico.
El eje de movimiento del vibrador A es preferentemente un eje de simetría del transductor ultrasónico TUS. Por ejemplo, el transductor ultrasónico TUS puede tener una sección circular, pasando el eje A por el centro del transductor ultrasónico TUS.
Según un modo de realización, la sonda PR comprende medios de control TOG para activar la medición.
La sonda PR según la figura 6a comprende, por tanto, un vibrador destinado a aplicar la vibración continua de baja frecuencia.
Según un modo de realización, el diámetro del transductor ultrasónico está comprendido entre 2 mm y 15 mm. Según un modo de realización, la frecuencia central del transductor ultrasónico está comprendida entre 1 MHz y 15 MHz.
Según un modo de realización, el transductor ultrasónico TUS es una sonda abdominal convexa.
Según un modo de realización, la sonda comprende un indicador de posicionamiento que se activa cuando la sonda está colocada correctamente. Este indicador puede ser un indicador visual, por ejemplo un cambio de color de los diodos. Alternativamente, el indicador puede ser un indicador sonoro o táctil tal como un cambio de tipo o de amplitud de una vibración.
La figura 6b ilustra un dispositivo de elastografía híbrida según la invención DEV.
El dispositivo DEV según la invención comprende:
- una sonda PR según la invención;
- una unidad central UC conectada a la sonda PR.
La unidad central puede comprender:
- medios de cálculo para el tratamiento de las señales ultrasónicas reflejadas;
- una pantalla SC para visualizar los resultados obtenidos en las diferentes etapas del procedimiento P según la invención;
- medios de control o de entrada ENT para el control del dispositivo por parte del operador.
La unidad central UC puede estar conectada a la sonda PR por conexión cableada o por medios de comunicación inalámbricos.
Según un modo de realización, la pantalla SC está adaptada para la visualización de los resultados ilustrados en la figura 5. La pantalla SC puede visualizar igualmente en tiempo real el indicador de posición RT_IP calculado durante la etapa CW_P del procedimiento P según la invención.
Según un modo de realización, la unidad central comprende medios configurados para activar automáticamente la aplicación de un impulso de baja frecuencia sobre la base del valor del indicador de posicionamiento RT_IP calculado y visualizado en tiempo real.
Según un modo de realización, la unidad central está comprendida en la sonda PR.
Claims (14)
1. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica que comprende una etapa de aplicación, con la ayuda de un vibrador comprendido en una sonda en contacto con un medio viscoelástico, de una vibración continua de baja frecuencia que genera una onda elástica en el interior del medio viscoelástico y de generación (CW), durante la propagación de la onda elástica, con la ayuda de un transductor ultrasónico en contacto con el medio viscoelástico, de una serie de adquisiciones ultrasónicas, comprendiendo la citada serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, siendo generados los grupos de adquisiciones ultrasónicas con una tasa de repetición (LPRF), comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición; siendo medido el parámetro de atenuación ultrasónica a partir de las adquisiciones ultrasónicas;
el citado procedimiento comprende además una etapa de determinación (CW_P), a partir de la serie de adquisiciones ultrasónicas, de al menos una propiedad de la propagación de la onda elástica en el interior del medio viscoelástico, y de utilización de la propiedad de la propagación de la onda elástica en el interior del medio viscoelástico para calcular un indicador de posicionamiento en tiempo real (RT_IP) de la sonda con respecto al medio viscoelástico que haya que estudiar.
2. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según la reivindicación precedente caracterizado por que comprende además una etapa de visualización en tiempo real del indicador de posicionamiento en tiempo real (RT_IP).
3. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones 1 a 2 caracterizado por que la propiedad de propagación de la onda elástica medida es elegida entre la amplitud de la onda elástica, la fase de la onda elástica, la velocidad de fase de la onda elástica, la elasticidad del medio viscoelástico, el módulo de Young del medio viscoelástico y el módulo de cizalladura del medio viscoelástico.
4. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones precedentes caracterizado por que la etapa de aplicación de una vibración continua de baja frecuencia se activa sólo si la fuerza de contacto entre el vibrador y el medio viscoelástico es superior a un umbral inferior predeterminado.
5. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones precedentes caracterizado por que la etapa de aplicación de una vibración continua de baja frecuencia se activa sólo si la fuerza de contacto entre el vibrador y el medio viscoelástico es superior a un umbral inferior predeterminado e inferior a un umbral superior predeterminado.
6. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones precedentes caracterizado por que la serie de adquisiciones ultrasónicas está formada por una repetición de grupos que comprenden al menos dos adquisiciones ultrasónicas que tienen una tasa de repetición dentro del grupo (HPRF) comprendida entre 500 Hz y 10 kHz y una tasa de repetición (LPRF) comprendida entre 10 Hz y 10 kHz.
7. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones precedentes caracterizado por que la tasa de repetición (LPRF) es menor que la frecuencia de vibración continua (cSWF).
8. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según una de las reivindicaciones precedentes caracterizado por que el parámetro de atenuación ultrasónica es un parámetro instantáneo y se calcula un coeficiente de calidad asociado con la medición del parámetro de atenuación ultrasónica instantáneo a partir de una propiedad de la propagación de la onda elástica en el interior del medio viscoelástico y/o de las propiedades de la onda ultrasónica.
9. Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica según la reivindicación precedente caracterizado por que se calcula un parámetro de atenuación ultrasónica medio a partir de una pluralidad de parámetros de atenuación ultrasónica instantáneos y de los coeficientes de calidad asociados con cada parámetro de atenuación ultrasónica instantáneo.
10. Sonda de medición de un parámetro de atenuación ultrasónica guiada por elastografía armónica que comprende:
- un vibrador configurado para aplicar una vibración continua de baja frecuencia al medio viscoelástico, generando la vibración continua de baja frecuencia una onda elástica en el interior del medio viscoelástico;
- un transductor ultrasónico configurado para emitir una serie de adquisiciones ultrasónicas, comprendiendo la citada serie de adquisiciones ultrasónicas grupos de adquisiciones ultrasónicas, siendo generados los grupos de adquisiciones ultrasónicas con una tasa de repetición, comprendiendo cada grupo de adquisiciones ultrasónicas al menos una adquisición;
- medios de cálculo y de visualización en tiempo real de un indicador de posicionamiento de la sonda, siendo calculado el citado indicador de posicionamiento a partir de una propiedad de la propagación de la onda elástica, siendo determinada la citada propiedad de propagación de la onda elástica a partir de la serie de adquisiciones ultrasónicas; siendo medido el parámetro de atenuación ultrasónica a partir de las adquisiciones ultrasónicas realizadas durante la aplicación de la vibración continua de baja frecuencia.
11. Sonda de medición de un parámetro de atenuación ultrasónica guiada por elastografía armónica según la reivindicación precedente, estando configurada además la citada sonda para aplicar la vibración continua cuando la fuerza de contacto entre la sonda y el medio viscoelástico es superior a un valor predeterminado.
12. Sonda de medición de la atenuación ultrasónica guiada por elastografía armónica según las reivindicaciones 10 u 11, caracterizada por que el transductor es llevado por el vibrador.
13. Sonda de medición de atenuación ultrasónica guiada por elastografía armónica según una de las reivindicaciones 11 o 12, configurada para calcular un parámetro de atenuación ultrasónica medio, siendo calculado el citado parámetro de atenuación ultrasónica medio a partir de una multiplicidad de parámetros de atenuación ultrasónica instantáneos y de coeficientes de calidad asociados con los parámetros de atenuación instantáneos.
14. Dispositivo de medición de atenuación ultrasónica guiada por elastografía armónica que comprende:
- una sonda (PR) para la medición de atenuación ultrasónica guiada por elastografía armónica según las reivindicaciones 10 u 11;
- una unidad central (UC) conectada a la sonda (PR) y que comprende al menos medios de cálculo para el tratamiento de las señales ultrasónicas reflejadas, medios de visualización (SC) y medios de control y/o de entrada (ENT).
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1851822A FR3078484A1 (fr) | 2018-03-02 | 2018-03-02 | Procede de mesure d’un parametre d’attenuation ultrasonore guide par elastographie harmonique, sonde et dispositif pour la mise en œuvre du procede |
| PCT/EP2019/054656 WO2019166393A1 (fr) | 2018-03-02 | 2019-02-26 | Procede de mesure d'un parametre d'attenuation ultrasonore guide par elastographie harmonique, sonde et dispositif pour la mise en oeuvre du procede |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2986863T3 true ES2986863T3 (es) | 2024-11-12 |
Family
ID=63490519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES19706655T Active ES2986863T3 (es) | 2018-03-02 | 2019-02-26 | Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica; sonda y dispositivo para la puesta en práctica del procedimiento |
Country Status (10)
| Country | Link |
|---|---|
| US (1) | US11808898B2 (es) |
| EP (1) | EP3758610B1 (es) |
| JP (1) | JP7295126B2 (es) |
| KR (1) | KR20200128085A (es) |
| CN (1) | CN110831505B (es) |
| AU (1) | AU2019229314A1 (es) |
| ES (1) | ES2986863T3 (es) |
| FR (1) | FR3078484A1 (es) |
| IL (1) | IL277012A (es) |
| WO (1) | WO2019166393A1 (es) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114379178B (zh) | 2020-10-05 | 2023-05-12 | Skc株式会社 | 玻璃接合用薄膜、薄膜的制备方法以及透光层叠体 |
| KR102760427B1 (ko) * | 2022-06-23 | 2025-02-03 | 광주과학기술원 | 광신호를 이용한 탄성체의 성질 측정 장치 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2843290B1 (fr) | 2002-08-08 | 2005-06-24 | Echosens | Dispositif et procede pour la mesure de l'elasticite d'un organe humain ou animal |
| US8100831B2 (en) | 2006-11-22 | 2012-01-24 | General Electric Company | Direct strain estimator for measuring elastic properties of tissue |
| CN102469989B (zh) | 2009-07-07 | 2014-04-16 | 株式会社日立医疗器械 | 超声波诊断装置和超声波测量方法 |
| FR2949965B1 (fr) | 2009-09-17 | 2012-09-28 | Echosens | Procede pour la mesure d'au moins une propriete de tissu biologique |
| CN101699280B (zh) | 2009-10-15 | 2011-08-17 | 北京索瑞特医学技术有限公司 | 超声无损检测粘弹性介质弹性的方法及其装置 |
| JP5677757B2 (ja) | 2010-03-16 | 2015-02-25 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
| US20140147013A1 (en) | 2010-10-11 | 2014-05-29 | The Regents Of The University Of Colorado, A Body Corporate | Direct echo particle image velocimetry flow vector mapping on ultrasound dicom images |
| CN102151152A (zh) * | 2011-03-01 | 2011-08-17 | 深圳市一体医疗科技股份有限公司 | 一种用于测量粘弹性介质弹性的测量探头、系统及方法 |
| US8734350B2 (en) | 2011-03-04 | 2014-05-27 | Mayo Foundation For Medical Education And Research | System and method for correcting errors in shear wave measurements arising from ultrasound beam geometry |
| US8948534B2 (en) | 2011-04-22 | 2015-02-03 | Kabushiki Kaisha Toshiba | MRI Gibbs' ringing filtering with edge-protection |
| FR2978657B1 (fr) | 2011-08-03 | 2013-08-30 | Echosens | Procede pour determiner en temps-reel une probabilite de presence d'un tissu biologique cible en regard d'un transducteur ultrasonore |
| CN102283679B (zh) * | 2011-08-04 | 2014-05-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
| US11529121B2 (en) * | 2013-02-19 | 2022-12-20 | Echosens | Multipulse elastography method |
| US20150141821A1 (en) * | 2013-04-03 | 2015-05-21 | Hitachi Aloka Medical, Ltd. | Ultrasonic diagnostic apparatus and elastic evaluation method |
| JP2017533031A (ja) | 2014-10-29 | 2017-11-09 | メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ | 超音波トランスデューサの連続振動による超音波エラストグラフィのための方法 |
| CN105455851B (zh) * | 2015-12-24 | 2018-03-13 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
| FR3054123B1 (fr) * | 2016-07-25 | 2021-11-26 | Echosens | Procede de mesure d’un parametre viscoelastique d’un organe humain ou animal |
-
2018
- 2018-03-02 FR FR1851822A patent/FR3078484A1/fr not_active Ceased
-
2019
- 2019-02-26 ES ES19706655T patent/ES2986863T3/es active Active
- 2019-02-26 US US16/977,369 patent/US11808898B2/en active Active
- 2019-02-26 CN CN201980003389.6A patent/CN110831505B/zh active Active
- 2019-02-26 AU AU2019229314A patent/AU2019229314A1/en not_active Abandoned
- 2019-02-26 KR KR1020207027877A patent/KR20200128085A/ko not_active Withdrawn
- 2019-02-26 WO PCT/EP2019/054656 patent/WO2019166393A1/fr not_active Ceased
- 2019-02-26 JP JP2020545644A patent/JP7295126B2/ja active Active
- 2019-02-26 EP EP19706655.8A patent/EP3758610B1/fr active Active
-
2020
- 2020-08-30 IL IL277012A patent/IL277012A/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| EP3758610A1 (fr) | 2021-01-06 |
| FR3078484A1 (fr) | 2019-09-06 |
| AU2019229314A1 (en) | 2020-10-01 |
| RU2020132278A (ru) | 2022-04-05 |
| IL277012A (en) | 2020-10-29 |
| WO2019166393A1 (fr) | 2019-09-06 |
| US20210048518A1 (en) | 2021-02-18 |
| JP7295126B2 (ja) | 2023-06-20 |
| CN110831505A (zh) | 2020-02-21 |
| CN110831505B (zh) | 2024-12-31 |
| RU2020132278A3 (es) | 2022-05-04 |
| EP3758610B1 (fr) | 2024-04-03 |
| KR20200128085A (ko) | 2020-11-11 |
| JP2021514765A (ja) | 2021-06-17 |
| US11808898B2 (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| ES2926791T3 (es) | Procedimiento de elastografía híbrida, sonda y dispositivo para elastografía híbrida | |
| US11730448B2 (en) | Device and method for measuring the viscoelastic properties of a viscoelastic medium | |
| US12226260B2 (en) | Method and device for measuring a mean value of visco-elasticity of a region of interest | |
| ES2939246T3 (es) | Métodos y aparatos de formación de imágenes para realizar la formación de imágenes por elastografía por ondas de corte | |
| JP2005535420A (ja) | 膀胱体積の超音波瞬時測定 | |
| ES2986863T3 (es) | Procedimiento de medición de un parámetro de atenuación ultrasónica guiado por elastografía armónica; sonda y dispositivo para la puesta en práctica del procedimiento | |
| CN112075955B (zh) | 粘弹性介质超声波参数测量方法与装置 | |
| CN112294365B (zh) | 表征组织的系统及相关方法 | |
| WO2017009516A1 (es) | Dispositivo emisor de ondas ultrasónicas de torsión y transductor que lo comprende | |
| RU2794039C2 (ru) | Способ гибридной эластографии, зонд и устройство для гибридной эластографии | |
| RU2786332C2 (ru) | Способ для измерения параметра затухания ультразвука методом гармонической эластографии; зонд и устройство для осуществления способа | |
| JP6537706B2 (ja) | 超音波撮像装置および超音波受信信号の処理方法 | |
| US12484881B2 (en) | Device and method for measuring the viscoelastic properties of a viscoelastic medium | |
| CN112839589B (zh) | 用于自动选择用于计算粘弹性介质的性质的深度范围的方法 | |
| US9360460B2 (en) | Device for calibrating and testing echotomographic equipment | |
| RU2728681C1 (ru) | Способ стандартизации измерений в ультразвуковой двумерной сдвиговолновой эластографии | |
| HK40043279A (en) | Method for automatically selecting a depth range for calculating a property of a viscoelastic medium | |
| Demi̇rkiran | Bir Fantom Dokuda Ultrasonik Sinyal Karakterizasyonunun Sıcaklığa Bağlılığının Deneysel Olarak Incelenmesi | |
| Caenen | Numerical and Experimental Assessment of Supersonic Shear Wave Imaging as New Non-Invasive Ultrasound Technique for Arterial Stiffness Characterization |