ES2803125B2 - Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos - Google Patents

Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos Download PDF

Info

Publication number
ES2803125B2
ES2803125B2 ES201930675A ES201930675A ES2803125B2 ES 2803125 B2 ES2803125 B2 ES 2803125B2 ES 201930675 A ES201930675 A ES 201930675A ES 201930675 A ES201930675 A ES 201930675A ES 2803125 B2 ES2803125 B2 ES 2803125B2
Authority
ES
Spain
Prior art keywords
soft solid
transducer
ultrasonic transducer
given
ultrasonic
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
Application number
ES201930675A
Other languages
English (en)
Other versions
ES2803125A1 (es
Inventor
González Noé Jiménez
Baviera José María Benlloch
Femenía Francisco Camarena
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Univ Polotecnica De Valencia
Consejo Superior de Investigaciones Cientificas CSIC
Original Assignee
Univ Polotecnica De Valencia
Consejo Superior de Investigaciones Cientificas CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Univ Polotecnica De Valencia, Consejo Superior de Investigaciones Cientificas CSIC filed Critical Univ Polotecnica De Valencia
Priority to ES201930675A priority Critical patent/ES2803125B2/es
Priority to PCT/ES2020/070457 priority patent/WO2021014040A1/es
Priority to US17/628,785 priority patent/US20220330919A1/en
Priority to EP20843425.8A priority patent/EP4005493A4/en
Publication of ES2803125A1 publication Critical patent/ES2803125A1/es
Application granted granted Critical
Publication of ES2803125B2 publication Critical patent/ES2803125B2/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • 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/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4494—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements
    • 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/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/043—Analysing solids in the interior, e.g. by shear waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22—Details, e.g. general constructional or apparatus details
    • G01N29/221—Arrangements for directing or focusing the acoustical waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22—Details, e.g. general constructional or apparatus details
    • G01N29/24—Probes
    • G01N29/2456—Focusing probes
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8925—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
    • 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
    • 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/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer

Description

DESCRIPCIÓN
PROCEDIMIENTO DE OBTENCIÓN DE PROPIEDADES ELÁSTICAS DE UN
SÓLIDO BLANDO QUE HACE USO DE VÓRTICES ACÚSTICOS
OBJETO DE LA INVENCIÓN
El objeto de la presente invención es un procedimiento de obtención propiedades elásticas de un sólido blando mediante ondas transversales quasi-omnidireccionales generadas por un haz ultrasónico focalizado de vórtice.
ANTECEDENTES DE LA INVENCIÓN
La imagen elastográfica es una modalidad de imagen médica que permite evaluar las propiedades elásticas de los tejidos blandos. Éstas permiten detectar cambios en la rigidez de los tejidos asociados a patologías subyacentes.
Los métodos elastográficos proponen la estimación de las propiedades elásticas de los tejidos midiendo las deformaciones que se producen cuando sobre ellos se aplica un determinado esfuerzo mecánico externo.
Por un lado, en los métodos quasi-estáticos el tejido es comprimido externamente de una manera análoga a la palpación, o aplicando una compresión oscilatoria externamente sobre el tejido.
Midiendo la amplitud de las deformaciones producidas se obtiene una imagen de la rigidez relativa del tejido como se propone en la técnica Sonoelasticity Imaging que figura en el documento Sono-elasticity imaging. De Lerner, R.M.; Parker, K. J. Kessler, L. W., ed. Acoustic imaging. New York: Plenum Co, 317-327. (1988). Dicho documento trata sobre el uso de técnicas de imagen Doppler para medir los desplazamientos del tejido que es sometido a un vibrador externo de baja frecuencia. Dichos desplazamientos aparecen en forma de ondas transversales y pueden ser empleados para obtener información sobre la elasticidad de tejidos biológicos.
Otros métodos, como la elastografía por compresión descrito en Elastography: a quantitative method for imaging the elasticity of biological tissues. de Ophir, J., Cespedes, I., Ponnekanti, H., Yazdi, Y., & Li, X.; Ultrasonic imaging, 13(2), 111-134 (1991) emplean una compresión mecánica para tomar dos imágenes de ultrasonidos en modo B-scan. La primera se toma como referencia y la segunda se obtiene tras aplicar una compresión mecánica mediante la misma sonda de ultrasonidos. Empleando técnicas de correlación cruzada se obtiene una imagen de las deformaciones producidas que proporcionan una imagen de la elasticidad del medio.
El inconveniente de los métodos anteriores es que la distribución de esfuerzos no es uniforme y depende de la geometría del medio, a priori desconocida, por lo que solo proporcionan una imagen cualitativa de la elasticidad. Por tanto, estos métodos quasiestáticos no permiten una evaluación cuantitativa de las propiedades elásticas de los tejidos.
Una segunda generación de métodos son los que emplean como mecanismo para generar el campo de esfuerzos la fuerza de radiación acústica producida por un haz de ultrasonidos focalizados, donde la transferencia de momento desde la onda al tejido es debida a la absorción y reflexión sobre las zonas no homogéneas del mismo.
Estos métodos pueden proporcionar una imagen cuantitativa de la elasticidad puesto que la deformación del tejido se realiza aplicando el esfuerzo en el interior del tejido y dicho esfuerzo es, en principio, conocido. En general, un haz primario de ultrasonidos se usa para producir una deformación en el tejido, mientras que un haz ultrasónico secundario en modo eco-impulso se emplea para adquirir un set de sucesivas imágenes en modo B-scan.
Empleando la correlación cruzada, con la que se comparan las distintas imágenes, se pueden detectar las deformaciones producidas en el interior del tejido. Comúnmente, dichos desplazamientos son de unos pocos micrómetros de amplitud. De esta manera, los parámetros elásticos del tejido pueden ser calculados midiendo las deformaciones producidas en el tejido cuando se conoce la fuerza de radicación aplicada.
Diferentes modalidades se han desarrollado empleando estos conceptos como la Acoustic Radiation Force Impulse imaging (ARFI), del documento On the feasibility of remóte palpation using acoustic radiation torce de Kathryn R. Rightingale, Mark L. Palmeri, Roger W. Nightingale, y Gregg E. Trahey, J.; Acoust. Soc. Am. 110 (1), de julio de 2001, en el que se propone una técnica capaz de obtener imágenes en las variaciones mecánicas de los tejidos. La técnica utiliza una fuerza de radiación en dirección axial (componentes lineales), y estudia los desplazamientos en la zona local del foco.
Otro ejemplo es amplitude-modulated Harmonic Motion Imaging (HMI), del documento Single-Element Focused Ultrasound Transducer Method for Harmonic Motion Imaging de Caroline Maleke, Mathieu Pernot y Elisa E. Konofagou; Ultrasonic Imaging 28, 144­ 158 (2006), en el que se propone el uso de un transductor focalizado excitado con una amplitud modulada a baja frecuencia para ejercer un empuje en el tejido localizado en su foco. Este documento parte de un haz que genera un empuje lineal, por tanto no helicoidal.
Otras modalidades emplean las ondas de cizalla (o transversales) que se generan tras la aplicación transitoria del haz ultrasónico primario. Puesto que las ondas de cizalla se propagan por los tejidos humanos con una velocidad lenta (alrededor de 1-10 m/s), las deformaciones producidas por las mismas pueden ser medidas por un haz de ultrasonidos secundario.
De esta manera, midiendo la velocidad de propagación local de las ondas transversales se obtiene un mapa del módulo elástico de cizalla del medio, ya que la velocidad de propagación de las ondas transversales viene determinada directamente por la rigidez de los tejidos.
Ejemplo de estas técnicas es Shear Wave Elastography Imaging (SWEI), del documento Shear wave elasticity imaging: a new ultrasonic technology of medical diagnostics. de Sarvazyan, A. P., Rudenko, O. V., Swanson, S. D., Fowlkes, J. B., & Emelianov, S. Y. Ultrasound in medicine & biology, 24(9), 1419-1435 (1998), que muestra una técnica para determinar las propiedades elásticas de un tejido basada en el uso de fuerza de radiación acústica para excitarlo. Al igual que en las técnicas previas, la fuerza de radiación acústica se ejerce en la dirección axial, pero en este caso se define un patrón de excitación espacial, es decir, se excita el tejido en diferentes puntos.
Además, Supersonic shear imaging (SSI), que figura en el documento Supersonic Shear Imaging: A New Technique for Soft Tissue Elasticity Mappin de Jéremy Bercoff, Mickáel Tanter, and Mathias Fink; ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 51, no. 4, abril de 2004. Este documento propone una nueva técnica para generar ondas transversales en el tejido biológico.
La técnica consiste en utilizar la fuerza de radiación ejercida por un haz ultrasónico que se focaliza en diferentes puntos del tejido a una velocidad superior a la de propagación de las ondas transversales que se generan en el tejido. Esto genera un cono de Match con ondas transversales de mayor amplitud, por la interferencia constructiva en los frentes de onda que se genera.
En todas las técnicas citadas las ondas transversales se generan en la dirección perpendicular al haz primario de ultrasonidos. De esta manera, el área sobre el eje del haz primario queda sin escanear puesto que las ondas transversales no se generan en esa dirección. Puesto que en estas técnicas se emplea un haz de ultrasonidos focalizado que solo transporta momento lineal, los esfuerzos producidos solo se producen en la dirección del haz axial.
Existen otros tipos de haces de ultrasonidos que, además de transportar momento lineal, transportan momento angular. Éstos son los haces de vórtice (o vórtices acústicos). Un ejemplo de dispositivo capaz de producir dichos haces de vórtice es el que se describe en el documento de Hefner, B. T., & Marston, P. L. (1999), An acoustical helicoidal wave transducer with applications for the alignment of ultrasonic and underwater systems; The Journal of the Acoustical Society of America, 106(6), 3313-3316. En este trabajo se presentan dos métodos, uno activo y otro pasivo, para generar vórtices acústicos. El método activo emplea cuatro transductores piezoeléctricos alimentados con diferente fase y el método pasivo emplea un único transductor con una superficie helicoidal.
Otros métodos incluyen rejillas en forma de espiral, como las descritas en los documentos Sharp acoustic vortex focusing by Fresnel-spiral zone plates. De Jiménez, N., Romero-García, V., García-Raffi, L. M., Camarena, F., & Staliunas, K. Applied Physics Letters, 112(20), 204101. (2018). y en el cocumento Formation of high-order acoustic Bessel beams by spiral diffraction gratings; de Jiménez, N., Picó, R., Sánchez Morcillo, V., Romero-García, V., García-Raffi, L. M., & Staliunas, K.; Physical Review E, 94(5), 053004. (2016).
No obstante, la fuerza de radiación aplicada en todas las técnicas elastográficas existentes en la actualidad tiene la dirección del haz de excitación generado, por lo que el patrón de radiación de las ondas transversales es limitado, no se puede definir la frecuencia de la excitación y la amplitud de las ondas transversales generadas es reducida.
DESCRIPCIÓN DE LA INVENCIÓN
La presente invención proporciona una mejora con respecto a los métodos anteriores, ya que amplía el patrón de radiación de las ondas transversales generadas, permite definir la frecuencia de excitación de las ondas transversales y aumenta la amplitud de las ondas transversales generadas.
La invención se basa en el uso de un haz ultrasónico de vórtice que produce un campo de esfuerzos de torsión, lo que genera un frente de ondas transversales en el interior de un sólido blando, y que sirve para determinar una serie de propiedades elásticas de dicho sólido blando, que es preferentemente un tejido, como el del hígado o la próstata, de un paciente al que se quiere diagnosticar.
Concretamente, con el presente procedimiento se genera un frente de ondas transversales quasi-omnidireccionales a partir de un haz ultrasónico focalizado, con un perfil de fase helicoidal, es decir, un vórtice acústico, propagándose el frente de ondas a través del sólido blando, con lo que se cubren todas las zonas de interés alrededor del haz ultrasónico.
Esto ofrece varias ventajas significativas respecto al resto de técnicas de imagen elastográficas del estado de la técnica. En primer lugar, las técnicas actuales sólo son capaces de generar el frente de ondas en la dirección perpendicular al eje del haz ultrasónico focalizado, de manera que no se propagan en la dirección axial al foco, por lo que resulta muy difícil extraer los parámetros elásticos de esa zona.
En segundo lugar, la amplitud de las ondas generadas es mayor para una misma intensidad acústica, lo que permite mejorar la relación señal-ruido de la imagen y reducir el nivel de amplitud del haz, reduciendo así efectos indeseados como el incremento de temperatura producida por el haz primario.
En tercer lugar, si se controla el sentido de giro angular del vórtice acústico y se varía en función del tiempo se puede controlar la frecuencia de las ondas transversales producidas y por tanto su longitud de onda en el interior del tejido. Finalmente, el control del sentido de giro permite controlar la polarización de las ondas transversales permitiendo evaluar la anisotropía de los tejidos, como por ejemplo ocurre en los tejidos fibrosos.
Preferentemente, este procedimiento se utiliza para obtener las propiedades elásticas de un tejido de un paciente. La información obtenida sobre las propiedades elásticas del tejido sirve para realizar un diagnóstico médico y detectar posibles anomalías en dicho tejido, que pueden ser consecuencia de un cáncer o de algún otro tipo de lesión, y que vienen acompañadas de cambios en las propiedades elásticas de los tejidos.
La primera etapa del procedimiento consiste en aplicar una señal pulsada o de amplitud modulada a un transductor ultrasónico que tiene una superficie destinada a contactar con el tejido que se va a estudiar. Esta señal está comprendida en el rango de los ultrasonidos (con una frecuencia portadora de entre 0.2 MHz y 20 MHz). Concretamente, se trataría de una señal sinusoidal de frecuencia de modulación igual a la del frente de ondas transversales que se desea generar.
Una vez aplicada la señal sobre el transductor ultrasónico, se genera un haz ultrasónico focalizado de vórtice, generando éste a su vez un frente de ondas transversales quasiomnidireccionales que se transmiten por el sólido blando.
El transductor ultrasónico utilizado puede ser de dos tipos distintos, y en función de cuál se utilice, la estrategia para generar el vórtice acústico será distinta.
En primer lugar, se puede utilizar un transductor de elemento simple, que comprende una lente holográfica. Dicha lente está destinada a posicionarse sobre la superficie (x0,y0) del transductor ultrasónico. El frente de ondas está caracterizado por una amplitud compleja A(x0,y0), que es modificada por la lente, de manera que se ajusta a la de un vórtice acústico focalizado. La fase viene dada por la ecuación:
Figure imgf000008_0001
(Ecuación 1)
dónde ^ (x 0,y0)es la amplitud compleja a lo largo de la superficie del transductor ultrasónico dada por x0,y0. El número de onda viene dado por k0 = 2 n f/c0, donde f es la frecuencia de la portadora y c0 es la velocidad del sonido en el sólido blando. F es la distancia focal de la lente y m la carga topológica del vórtice, que normalmente es un número entero. En función del signo de m, el vórtice gira en sentido horario o antihorario.
La lente, por tanto, ha de ser capaz de producir el perfil de fase A(x0,y0). Para ello, una estrategia es dividir la lente en píxeles y definir una altura para cada píxel como h(x0,y0) de manera que cumpla:
Figure imgf000008_0002
(Ecuación 2)
siendo Z = pLcL la impedancia de la lente, donde pL es la densidad de la lente y cL la velocidad de propagación de los ultrasonidos en la lente. kL = 2 n f/cL, y d es una distancia arbitraria que coincide con el plano de superficie del transductor ultrasónico. La obtención de las alturas h(x0,y0) en función d e ^(x0,y0) se realiza mediante inversión numérica de la ecuación 2.
En segundo lugar, se puede utilizar un transductor ultrasónico de elemento múltiple (o phased-array). En el caso de que el array sea plano, cada elemento del transductor se ajustará a una amplitud dada por |^|, y una fase dada por tan~1{¡m(A)/Re(A)), donde A viene dada por la Ecuación 2, Re(-) indica la parte real e Im(-) la parte imaginaria del valor complejo. Las coordenadas x0 e y0 vienen dadas por las posiciones espaciales en coordenadas cartesianas de cada elemento del transductor ultrasónico.
En el caso de que el array sea de focalización geométrica, como en el caso de un multielemento donde cada elemento se dispone sobre la superficie de una esfera de radio F, se aplica el mismo método que en el caso anterior, pero la expresión a emplear es:
A(x0,y0) = exp(-¿mtan 1(x0,y0)). (Ecuación 3)
Es decir, se trata de un perfil de fase que depende linealmente del ángulo polar que ocupa cada elemento del phased-array.
Por tanto, independientemente del tipo de transductor ultrasónico que se utilice, el haz ultrasónico emitido posee una intensidad acústica que rota con respecto de la coordenada angular, y transfiere al sólido blando tanto una cantidad de momento lineal en la dirección del haz ultrasónico, como un momento angular, en forma de toro alrededor del haz ultrasónico.
En particular se transmite al sólido blando un campo de fuerzas que se puede calcular como:
F(x,y,z) = i - 2 ^ ^ P i q L C q (pvP* - P*VP) (Ecuación 4)
Siendo F el campo vectorial de fuerzas, a (^ ) la absorción acústica del sólido blando,
m = 2n f la frecuencia angular, p0 la densidad, c0 la velocidad de las ondas longitudinales de ultrasonidos, P el campo de presión de ultrasonidos, y P* su conjugada compleja.
A la vista del campo de fuerzas anterior, se comprueba que cuando el campo de ultrasonidos es de tipo vórtice, en el sólido blando se produce una fuerza que es de tipo torque o torsional, con una pequeña componente axial. Puesto que el sólido blando absorbe gran parte de la energía del haz ultrasónico focalizado, la transferencia de momento angular en forma de torque al sólido blando provoca una deformación transitoria del mismo, enroscándolo.
De esta manera, se genera el frente de ondas transversales que se propaga no solo en la dirección perpendicular al haz ultrasónico, sino, además, en la misma dirección del haz ultrasónico, es decir, un frente de ondas quasi-omnidireccionales.
Otra de las ventajas que ofrece la presente invención es, que al poder controlar los parámetros que definen el haz ultrasónico, se puede controlar la polarización del frente de ondas que se genera en el interior del sólido blando. Concretamente, controlando el signo de la carga topológica del haz ultrasónico, se puede controlar el sentido de rotación del esfuerzo producido (horario/antihorario).
Como se ha indicado, la carga topológica será preferentemente igual a uno, aunque si se hace mayor a la unidad se generan campos de fuerza más anchos con un torque mayor.
Controlando el sentido de rotación, alternando entre uno y otro periódicamente, las ondas se excitan tanto en los ciclos positivos como en los negativos, lo que logra inducir en el sólido blando ciclos de deformación positivos y negativos, incrementando la amplitud del frente de ondas generadas y, por tanto, la robustez y sensibilidad de la técnica. De esta manera no hay necesidad de esperar a que el sólido blando se relaje para volver a empujarlo y generar ondas de manera continua.
En el caso de que el transductor ultrasónico sea de elemento simple existen diferentes estrategias que permiten controlar el signo de la carga topológica.
La primera consiste en utilizar una lente diseñada para trabajar con una carga topológica en una frecuencia de emisión y con otra carga topológica de signo contrario en otra frecuencia, para ir alternando entre ambas.
La segunda estrategia consiste en utilizar dos transductores ultrasónicos, posicionados como dos anillos concéntricos, teniendo cada uno de ellos una lente diferente, así como una carga topológica diferente, de signos contrarios, y alternando la emisión de uno y otro.
En el caso de que el transductor ultrasónico sea de elemento múltiple, el control del signo de la carga topológica es más sencillo, pues únicamente hay que invertir angularmente la fase de los elementos del array, es decir, invistiendo el signo del parámetro m en la Ecuación l ó e n l a Ecuación 3, respectivamente.
Por lo tanto, el control de la frecuencia y del sentido de rotación del haz ultrasónico permite el control del frente de ondas transversales generadas, lo que facilita la realización de estudios elastográficos a diferentes frecuencias.
La siguiente etapa del procedimiento, una vez generado el frente de ondas transversales, consiste en la adquisición de señales de radiofrecuencia que son reflejadas por el sólido blando en distintos instantes de tiempo, mientras se propaga dicho frente de ondas. Para ello se puede hacer uso de un segundo transductor de imagen médica por ultrasonidos, en modo eco-impulso. Este transductor secundario se emplea para obtener una serie de imágenes ultrasónicas en diferentes instantes de tiempo tomadas posteriormente o durante la activación del transductor primario.
Una vez obtenida la serie de imágenes, la siguiente etapa del procedimiento consiste en el cálculo de las deformaciones mediante métodos de correlación cruzada o técnicas Doppler entre las diferentes imágenes. Ello aporta una imagen de las deformaciones producidas en el tejido al paso de las ondas transversales.
A partir de las imágenes de las deformaciones, usando técnicas estándar de tracking, se calcula la velocidad de propagación del frente de ondas transversales.
Finalmente, a partir de las velocidades se obtiene el módulo de elasticidad transversal o de cizalla, que para un medio elástico se puede obtener a partir de la ecuación:
Figure imgf000011_0001
(Ecuación 5)
Siendo G el módulo de elasticidad transversal o de cizalla y p la densidad del sólido blando.
La densidad cambia muy poco respecto a la variación que sufre el módulo de elasticidad transversal del sólido blando, por lo que una diferencia en la velocidad es debida fundamentalmente a una variación en el módulo de elasticidad transversal de éste y, por lo tanto, a algún tipo de alteración en el sólido blando analizado.
Finalmente, a partir del módulo elástico transversal obtenido en distintos puntos del tejido, se obtienen las imágenes elastográficas que sirven para realizar un diagnóstico médico.
DESCRIPCIÓN DE LOS DIBUJOS
Para complementar la descripción que se está realizando y con objeto de ayudar a una mejor comprensión de las características de la invención, de acuerdo con un ejemplo preferente de realización práctica de la misma, se acompaña como parte integrante de dicha descripción, un juego de dibujos en donde con carácter ilustrativo y no limitativo, se ha representado lo siguiente:
Figura 1.- Muestra un diagrama de los transductores ultrasónicos primario y secundario de los que hace uso el procedimiento.
Figura 2.- Muestra un diagrama de bloques del proceso donde se muestra una posible secuencia a seguir.
Figura 3.- Muestra el campo acústico generado por el transductor ultrasónico.
Figura 4.- Muestra el campo de fuerza acústica de radiación generada sobre el sólido blando.
Figura 5.- Muestra el desplazamiento del sólido blando en la dirección z en diferentes instantes de tiempo.
REALIZACIÓN PREFERENTE DE LA INVENCIÓN
A la vista de las figuras descritas anteriormente, se puede observar un ejemplo de realización no limitativo del procedimiento de obtención de propiedades elásticas de un sólido blando, objeto de esta invención.
La primera etapa del procedimiento, del que se muestra un diagrama de bloques en la figura 2, consiste en aplicar una señal pulsada o de frecuencia modulada, de frecuencia portadora en torno a 1 MHz, comprendida en el rango de los ultrasonidos, y frecuencia moduladora en el rango de 1 Hz a 1000 Hz, a un transductor ultrasónico (1), como el de la figura 1, que comprende una superficie destinada a contactar con un sólido blando.
Una vez aplicada la señal sobre el transductor ultrasónico (1), se genera un haz ultrasónico (5) focalizado, con un perfil de fase helicoidal, es decir, un vórtice acústico, y que genera un frente de ondas (6) transversales quasi-omnidireccionales que se transmiten por el sólido blando. La frecuencia del frente de ondas (6) es igual a la frecuencia de modulación de la señal pulsada aplicada sobre el transductor ultrasónico
( 1).
El haz ultrasónico (5) focalizado de vórtice se genera mediante el transductor ultrasónico (1) que es de elemento múltiple (o phased-array), siendo el array de focalización geométrica. Para ello, se ajusta cada elemento del transductor ultrasónico (1) a una amplitud dada por:
\A(x0,y0)\ = exp(—¿ m • tan_1(x0,y0)), (Ecuación 6)
es decir, un perfil de fase que depende linealmente del ángulo polar que ocupa cada elemento del transductor ultrasónico (1).
Por tanto, el haz ultrasónico (5) emitido posee una intensidad acústica que rota con respecto de la coordenada angular, que transfiere al sólido blando tanto una cantidad de momento lineal en la dirección del haz ultrasónico, como un momento angular, en forma de toro alrededor del haz ultrasónico (5).
En la figura 3 se muestra el campo acústico generado por el transductor ultrasónico (1). En la imagen a) se representa la magnitud del campo en el plano sagital a la dirección de propagación y=0. En la b), la magnitud del campo en el plano transversal, sobre la distancia focal z=F. en la figura c) la fase del campo en el plano transversal, sobre la distancia focal z=F.
Se observa en la figura 3 cómo en el eje se produce una singularidad de fase que da lugar a un vórtice acústico. Así mismo, alrededor del foco la fase rota un número entero de veces.
La transferencia de momento lineal genera en el sólido blando un campo de fuerzas que se puede calcular como:
Figure imgf000014_0001
Siendo F el campo vectorial de fuerzas, a (^ ) la absorción del sólido blando, m la frecuencia angular, p0 la densidad, c0 la velocidad de las ondas transversales del frente, P el campo de presión producido, y P* su conjugada compleja.
Este campo de fuerzas se muestra en la figura 4. En la gráfica a) aparece una representación, en el plano transversal, de la componente de la fuerza en la dirección x, calculado en z=F. La b) es la representación en el plano transversal de la componente de la fuerza en la dirección y, calculado en z=F. La figura c) es una representación, en el plano sagital, de la componente de la fuerza en la dirección z, calculado en y=0. La d) es la representación en el plano transversal de la componente de torque de la fuerza, calculado en z=F. La subfigura e) es la representación del campo vectorial.
A partir del campo de fuerzas anterior, en el sólido blando se produce una fuerza que es de tipo torque, con una pequeña componente axial. Puesto que el sólido blando absorbe gran parte de la energía del haz ultrasónico (5), la transferencia de momento angular en forma de torque al sólido blando provoca una deformación transitoria del mismo, enroscándolo.
La siguiente etapa del procedimiento consiste en la adquisición de señales de radiofrecuencia que son reflejadas por el sólido blando en distintos instantes de tiempo, proceso que se repite mientras se propaga el frente de ondas (6) transversales. Para ello se usa un segundo transductor (2) de imagen médica por ultrasonidos, en modo eco-impulso.
Una vez obtenidas las señales, se calculan las deformaciones producidas en función del tiempo, a partir de la amplitud de los desplazamientos transversales que sufre el sólido blando, mediante correlación cruzada o técnicas Doppler. Estas deformaciones se pueden observar en las gráficas de la figura 5, en las que se refleja el desplazamiento del tejido en la dirección z, en diferentes instantes de tiempo, desde t=0.6 ms, hasta t=2.4 ms.
A partir de estas deformaciones, usando técnicas estándar de tracking, se calcula la velocidad de propagación del frente de ondas (6) transversales.
A partir de las velocidades se obtiene el módulo de elasticidad transversal o de cizalla, a partir de la ecuación:
v = J - (Ecuación 8)
Siendo G el módulo de elasticidad transversal o de cizalla y p la densidad del sólido blando.
Finalmente, a partir del módulo elástico transversal obtenido en distintos puntos del sólido blando, se pueden obtener unas imágenes elastográficas que sirvan para realizar undiagnósticomédico.

Claims (8)

REIVINDICACIONES
1.- Procedimiento de obtención de propiedades elásticas de un sólido blando sobre el que se ejerce una fuerza de radiación acústica que provoca unas deformaciones en dicho sólido blando, caracterizado por que comprende las etapas de:
- aplicación de una señal pulsada o de amplitud modulada a un transductor ultrasónico (1),
- generación, en el transductor ultrasónico (1), de un haz ultrasónico (5) focalizado de vórtice, que genera un frente de ondas (6) transversales quasiomnidireccionales, caracterizado por una velocidad, que se transmite por el sólido blando,
- adquisición de imágenes del sólido blando mientras se propaga el frente de ondas (6), haciendo uso de un segundo transductor ultrasónico (2) en contacto con el sólido blando,
- cálculo de las deformaciones producidas en el sólido blando por correlación cruzada o técnicas Doppler a partir de las imágenes,
- cálculo de la velocidad de propagación del frente de ondas (6) a partir de las deformaciones, usando técnicas estándar de tracking,
- cálculo del módulo de elasticidad transversal (G) del sólido blando a partir de la ecuación:
G
v =
M P
en la que v es la velocidad de propagación del frente de ondas (6) y p la densidad del sólido blando, y
- obtención de imágenes elastográficas a partir del módulo de elasticidad transversal en distintos puntos del sólido blando.
2.- El procedimiento de la reivindicación 1, en el que la fuerza de radiación acústica está dada por un campo de fuerzas que viene definido por:
a(w)
F(x,y,z) = i- (PVP* - P*VP)
2 up 0c0
siendo F el campo vectorial de fuerzas, a (^ ) la absorción del sólido blando, m la frecuencia angular, p0 la densidad del sólido blando, c0 la velocidad del frente de ondas (6), P el campo de presión producido, y P* su conjugada compleja.
3. - El procedimiento de la reivindicación 1, en el que el transductor ultrasónico (1) es un transductor de elemento simple, que comprende una lente (4) holográfica, posicionada sobre la superficie del transductor ultrasónico, modificando la lente holográfica (4) la fase del frente de ondas de manera que se ajuste a la de un vórtice acústico focalizado, dado por:
A(x0, y 0) = exp(-¿A:0Vxo 7o + f 2 ) exp (-ím ta n_1(y0,x 0)) siendo A(x0,y0) la fase a lo largo de la superficie del transductor ultrasónico (1) dada por x0, y0,k 0 = 2 n f/c0 el número de onda, donde f es la frecuencia y c0 es la velocidad del sonido en el sólido blando, F es la distancia focal de la lente (4) y m la carga topológica del vórtice.
4. - El procedimiento de la reivindicación 1, en el que el transductor ultrasónico (1) es un transductor de elemento múltiple plano, que comprende unos elementos que se ajustan a una amplitud dada por \A\ y una fase dada por tan~1(lm(A)/Re(A)) , siendo:
Figure imgf000017_0001
siendo A(x0,y0) la fase a lo largo de la superficie del transductor ultrasónico (1) dada por x0,y0 que representan la posición espacial en coordenadas cartesianas de cada elemento del transductor primario (1), k0 = 2 n f/c0 el número de onda, donde fe s la frecuencia y c0 es la velocidad del sonido en el sólido blando, F es la distancia focal de la lente (4) y m la carga topológica del vórtice.
5. - El procedimiento de la reivindicación 1, en el que el transductor ultrasónico (1) es un transductor de elemento múltiple con focalización geométrica, en que cada elemento del transductor ultrasónico (1) se ajusta a una amplitud dada por \A\ y una fase dada por tan~1(¡m(A)/Re(A)), siendo:
A{x0,y0) = exp(-ím tan“ 1(x0,y0))
donde A(x0,y0) es la fase a lo largo de la superficie del transductor ultrasónico (1) dada por x0,y0 que representan la posición espacial en coordenadas cartesianas de cada elemento del transductor primario (1) ,ym la carga topológica del vórtice.
6. - El procedimiento de la reivindicación 3, en el que se varía el signo de la carga topológica m modificando la fase del frente de ondas (6) a la salida del transductor primario (1) utilizando la lente (4), que trabaja con una carga topológica m positiva a una primera frecuencia y con una carga topológica m negativa a una segunda frecuencia.
7. - El procedimiento de la reivindicación 3, en el que se varía el signo de la carga topológica m utilizando dos transductores ultrasónicos (1) posicionados en forma de anillos concéntricos en el que cada uno comprende una lente (4) diferente con una carga topológica m distinta, siendo una positiva y otra negativa, y alternando la emisión entre un transductor ultrasónico y otro (1).
8. - El procedimiento de las reivindicaciones 4 o 5, en el que se varía el signo de la carga topológica m invirtiendo angularmente la fase de los elementos del transductor ultrasónico (1).
ES201930675A 2019-07-22 2019-07-22 Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos Active ES2803125B2 (es)

Priority Applications (4)

Application Number Priority Date Filing Date Title
ES201930675A ES2803125B2 (es) 2019-07-22 2019-07-22 Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos
PCT/ES2020/070457 WO2021014040A1 (es) 2019-07-22 2020-07-14 Procedimiento de obtención de propiedades elásticas de un sólido blando que hace uso de vórtices acústicos
US17/628,785 US20220330919A1 (en) 2019-07-22 2020-07-14 Method for obtaining elastic properties of a soft solid, which uses acoustic vortices
EP20843425.8A EP4005493A4 (en) 2019-07-22 2020-07-14 PROCESS FOR OBTAINING ELASTIC PROPERTIES OF A SOFT SOLID USING ACOUSTIC VORTEXES

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
ES201930675A ES2803125B2 (es) 2019-07-22 2019-07-22 Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos

Publications (2)

Publication Number Publication Date
ES2803125A1 ES2803125A1 (es) 2021-01-22
ES2803125B2 true ES2803125B2 (es) 2021-08-13

Family

ID=74186287

Family Applications (1)

Application Number Title Priority Date Filing Date
ES201930675A Active ES2803125B2 (es) 2019-07-22 2019-07-22 Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos

Country Status (4)

Country Link
US (1) US20220330919A1 (es)
EP (1) EP4005493A4 (es)
ES (1) ES2803125B2 (es)
WO (1) WO2021014040A1 (es)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI803021B (zh) * 2021-10-20 2023-05-21 國立成功大學 用於生物組織彈性量測的方法、超音波裝置與分析裝置
CN114815017A (zh) * 2022-04-24 2022-07-29 武汉大学 基于声与非均匀流体作用的可调液体透镜及制备方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105392428B (zh) * 2013-06-26 2019-06-04 皇家飞利浦有限公司 用于映射超声剪切波弹性成像测量的系统和方法
BR112015032724B1 (pt) * 2013-06-26 2022-03-08 Koninklijke Philips N.V. Sistema de elastografia ultrassônica e método de elastografia ultrassônica para inspecionar um sítio anatômico
BR112017008908B1 (pt) * 2014-10-28 2022-10-18 Supersonic Imagine Métodos e aparelhos de formação de imagem para realizar formação de imagem por elastografia com onda de cisalhamento
EP3236857B1 (en) * 2014-12-24 2020-05-06 SuperSonic Imagine Shear wave elastrography method and apparatus for imaging an anisotropic medium
RU2702090C2 (ru) * 2015-03-31 2019-10-04 Конинклейке Филипс Н.В. Калибровка ультразвукового, основанного на эластичности, отображения границы очага поражения
EP3475992A4 (en) * 2016-06-22 2020-02-19 Duke University ULTRASONIC TRANSDUCERS FOR CONSTRUCTIVE SHEAR WAVE INTERFERENCE AND RELATED METHODS AND SYSTEMS
CN109261472B (zh) * 2018-08-30 2019-11-08 西安交通大学 一种空间聚焦涡旋声场的产生装置及方法

Also Published As

Publication number Publication date
EP4005493A1 (en) 2022-06-01
ES2803125A1 (es) 2021-01-22
EP4005493A4 (en) 2023-01-25
WO2021014040A1 (es) 2021-01-28
US20220330919A1 (en) 2022-10-20

Similar Documents

Publication Publication Date Title
US11413007B2 (en) Non-contact acoustic radiation force based (ARF-based) generation of broad bandwidth mechanical waves using air-coupled ultrasound
Ambroziński et al. Acoustic micro-tapping for non-contact 4D imaging of tissue elasticity
US9084559B2 (en) Imaging method, displacement measurement method and apparatus
Montaldo et al. Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography
Ambroziński et al. Air-coupled acoustic radiation force for non-contact generation of broadband mechanical waves in soft media
Chen et al. Comparison of stress field forming methods for vibro-acoustography
Iwasaki et al. Monitoring of high-intensity focused ultrasound treatment by shear wave elastography induced by two-dimensional-array therapeutic transducer
ES2803125B2 (es) Procedimiento de obtencion de propiedades elasticas de un solido blando que hace uso de vortices acusticos
Callé et al. Temporal analysis of tissue displacement induced by a transient ultrasound radiation force
CN104323762A (zh) 一种基于光声显微成像的鲜红斑痣血管定量化检测装置
JP2022048984A (ja) 媒体の超音波特性評価用の方法とシステム
JP2017003587A (ja) ハイブリッド光学音響断層撮影及び超音波検査に関するデバイスおよび方法
Nabavizadeh et al. Multi-source and multi-directional shear wave generation with intersecting steered ultrasound push beams
Karzova et al. Shock formation and nonlinear saturation effects in the ultrasound field of a diagnostic curvilinear probe
Jiménez et al. A new elastographic technique using acoustic vortices
CN111449629B (zh) 一种光学相干弹性成像方法及装置
Burgess et al. Fast qualitative two-dimensional mapping of ultrasound fields with acoustic cavitation-enhanced ultrasound imaging
CN103313649A (zh) 测量装置
US20140121490A1 (en) Acousto-optic imaging device
US10080550B2 (en) Ultrasonic apparatus and control method for the same
Graham et al. Theoretical application of short-lag spatial coherence to photoacoustic imaging
JP2017047185A (ja) 超音波装置
González-Mateo et al. Quasi-omnidirectional shear wave generation using acoustic vortices for elastography
JP6816898B2 (ja) イメージング装置
Maleke et al. Single-element focused transducer method for harmonic motion imaging

Legal Events

Date Code Title Description
BA2A Patent application published

Ref document number: 2803125

Country of ref document: ES

Kind code of ref document: A1

Effective date: 20210122

FG2A Definitive protection

Ref document number: 2803125

Country of ref document: ES

Kind code of ref document: B2

Effective date: 20210813