ES2523982T3 - Procedimiento y dispositivo de formación de imágenes por ultrasonido - Google Patents

Procedimiento y dispositivo de formación de imágenes por ultrasonido Download PDF

Info

Publication number
ES2523982T3
ES2523982T3 ES11799794.0T ES11799794T ES2523982T3 ES 2523982 T3 ES2523982 T3 ES 2523982T3 ES 11799794 T ES11799794 T ES 11799794T ES 2523982 T3 ES2523982 T3 ES 2523982T3
Authority
ES
Spain
Prior art keywords
differential
stage
along
measurement
designates
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
ES11799794.0T
Other languages
English (en)
Inventor
Olivier Couture
Mickael Tanter
Mathias Fink
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.)
Centre National de la Recherche Scientifique CNRS
Institut National de la Sante et de la Recherche Medicale INSERM
Universite Paris Diderot Paris 7
Original Assignee
Centre National de la Recherche Scientifique CNRS
Institut National de la Sante et de la Recherche Medicale INSERM
Universite Paris Diderot Paris 7
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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=44201042&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=ES2523982(T3) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Centre National de la Recherche Scientifique CNRS, Institut National de la Sante et de la Recherche Medicale INSERM, Universite Paris Diderot Paris 7 filed Critical Centre National de la Recherche Scientifique CNRS
Application granted granted Critical
Publication of ES2523982T3 publication Critical patent/ES2523982T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • G01S7/52041Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target detecting modification of a contrast enhancer, e.g. detecting the destruction of a contrast agent by an acoustic wave, e.g. loss of correlation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52046Techniques for image enhancement involving transmitter or receiver

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Nonlinear Science (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

Procedimiento de formación de imágenes por ultrasonido para obtener imágenes de un campo de observación (8) en un medio para formar imágenes (2) que contiene unos difusores (7), constando dicho procedimiento de varias etapas sucesivas de medición a lo largo de cada una de las cuales: - se hace que una red (3) de transductores (T1-Tn) emita una onda ultrasónica incidente en el campo de observación (8), con una longitud de onda central ; - a continuación se graban las señales en bruto Sj (i, t) captadas por cada transductor (T1-Tn) y representativas de una onda ultrasónica refleja reverberada por los difusores (7) del medio a partir de la onda incidente, siendo i un índice que designa cada transductor, siendo j un índice que designa cada etapa de medición y designando t el tiempo, extendiéndose la red (3) de transductores a lo largo de al menos una dimensión y propagándose principalmente las ondas incidentes en una dirección de propagación perpendicular a la red de transductores, caracterizado porque se genera un número máximo C de blancos diferenciales, diferentes de una etapa de medición a otra, siendo cada blanco diferencial un difusor (7) que está presente en el campo de observación durante una etapa de medición y ausente durante una etapa de medición inmediatamente contigua, siendo el número C como máximo igual a INT (A/(5λ)2)+1, donde A es un área del campo de observación, y porque dicho procedimiento consta, además, de las siguientes etapas: - una etapa de tratamiento diferencial a lo largo de la cual se comparan las señales en bruto Sj (i, t) correspondientes a las etapas de medición sucesivas de índice j para extraer de estas unas señales diferenciales Vj (i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición; - una etapa de ajuste a lo largo de la cual se determina al menos una función y >= Pj (x) correspondiente a cada señal diferencial Vj (i, t), donde x es una variable de espacio que designa una posición perpendicularmente a la dirección de propagación e y es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto t; - y una etapa de localización a lo largo de la cual se determina un vértice Aj (x0, y0) de dicha función Pj, que corresponde a la posición del blanco diferencial.

Description

15
25
35
45
55
65
E11799794
14-11-2014
DESCRIPCIÓN
Procedimiento y dispositivo de formación de imágenes por ultrasonido
Sector de la técnica
La presente invención se refiere a los procedimientos y dispositivos de formación de imágenes por ultrasonido.
De manera más particular, la invención se refiere a un procedimiento de formación de imágenes por ultrasonido para obtener imágenes de un campo de observación en un medio para formar imágenes que contiene unos difusores, comprendiendo dicho procedimiento varias etapas sucesivas de medición a lo largo de cada una de las cuales:
-se hace que una red de transductores emita una onda ultrasónica incidente en el campo de observación, con una
longitud de onda central λ; -a continuación se graban las señales en bruto Sj (i, t) captadas por cada transductor y representativas de una
onda ultrasónica refleja reverberada por los difusores del medio a partir de la onda incidente, siendo i un índice
que designa cada transductor, siendo j un índice que designa cada etapa de medición y designando t el tiempo,
extendiéndose la red de transductores a lo largo de al menos una dimensión y propagándose las ondas incidentes principalmente en una dirección de propagación perpendicular a la red de transductores.
Estado de la técnica
En los procedimientos conocidos de este tipo las señales en bruto de los diferentes transductores se combinan entre sí para cada disparo, por ejemplo mediante la formación de trazas, para obtener las imágenes correspondientes a cada disparo. La resolución de la imagen está tradicionalmente limitada a λ/2, donde λ es la longitud de onda de los ultrasonidos (para los ultrasonidos a 1,5 MHz, λ es del orden de 1 mm en los tejidos humanos blandos).
Objeto de la invención
La presente invención tiene en particular como objetivo obtener imágenes por ultrasonido con una resolución mejor que λ/2.
Para ello, de acuerdo con la invención, un procedimiento del tipo en cuestión se caracteriza por que genera un número máximo C de blancos diferenciales, diferentes de una etapa de medición a otra, siendo cada blanco diferencial un difusor que está presente en el campo de observación durante una etapa de medición y ausente durante una etapa de medición inmediatamente contigua, siendo el número C como máximo igual a INT (A/(5λ)2)+1, donde A es un área del campo de observación, y por que dicho dispositivo consta, además, de las siguientes etapas:
-una etapa de tratamiento diferencial a lo largo de la cual se comparan las señales en bruto Sj (i, t)
correspondientes a unas etapas de medición sucesivas de índice j para extraer de estas unas señales
diferenciales Vj (i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición; -una etapa de ajuste a lo largo de la cual se determina al menos una función y = Pj (x) correspondiente a cada
señal diferencial Vj(i, t), donde x es una variable de espacio que designa una posición perpendicularmente a la
dirección de propagación e y es una coordenada que designa la posición de un punto a lo largo de la dirección
de propagación correspondiente a un tiempo de trayecto t (y = c.t/2 si t es el tiempo de ida y vuelta entre la
emisión de la onda incidente y la detección de la onda refleja); -y una etapa de localización a lo largo de la cual se determina un vértice Aj(x0, y0) de dicha función Pj, que
corresponde a la posición del blanco diferencial.
Hay que señalar que los blancos diferenciales mencionados con anterioridad pueden corresponder por ejemplo a destrucciones sucesivas de microburbujas o elementos similares bajo el efecto de ondas ultrasónicas.
Por medio de estas disposiciones, los difusores que forman los blancos diferenciales mencionados con anterioridad, se activan individualmente y se identifican individualmente a partir de las señales en bruto captadas y no a partir de una imagen ya formada, de tal modo que es posible identificarlos de forma muy precisa en el espacio con una precisión muy inferior a λ/2, que puede llegar por ejemplo hasta λ/200. Las posiciones de los difusores así identificados pueden servir a continuación para construir una imagen del medio, con una resolución muy inferior a λ/2, por ejemplo del orden de entre λ/50 y λ/200. De este modo, la invención propuesta permite mejorar de forma muy significativa la resolución de la formación de imágenes por ultrasonido.
De manera más precisa, la invención permite activar sólo un número limitado de difusores a la vez con el fin de poder distinguir el frente de onda individual procedente de un blanco que se quiere localizar. Al ya no ser el frente de onda una superposición del eco de varios blancos, su forma está definida de forma precisa por la posición de un difusor diferente. Esta posición se puede establecer en prácticamente unos micrómetros. Con el objetivo de solo activar un número limitado de difusores, la invención puede utilizar diferentes efectos de umbral generando o
15
25
35
45
55
65
E11799794
14-11-2014
suprimiendo unos difusores puntuales de ultrasonidos: en particular, los agentes de contraste por ultrasonidos (microburbujas, microgotas o liposomas) se ven afectados por los ultrasonidos que pueden inducir su explosión o su vaporización. Esto genera un blanco que aparece y/o desaparece en un tiempo muy breve, del orden del milisegundo o de algunos milisegundos. El frente de onda vinculado a este difusor puntual se puede obtener en particular restando las señales retrodifundidas cuando el blanco está presente, las señales retrodifundidas cuando el blanco está ausente. En diferentes formas de realización del procedimiento de acuerdo con la invención eventualmente se puede recurrir a, además, una y/u otra de las siguientes disposiciones:
-el número C es como máximo igual a 2, de preferencia igual a 1;
-a lo largo de la etapa de ajuste, se determina la función y = Pj(x) ajustando dicha función para minimizar las diferencias con los puntos Dj (xi, yi), donde xi es una variable de espacio que designa una posición de cada transductor i perpendicularmente a la dirección de propagación e yi es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto ti característico de la señal Vj(i, t) (puede tratarse, por ejemplo, del máximo de la señal Vj(i, t) o del frente de onda de la señal Vj(i, t), u otra);
-dicha función P es parabólica; -la etapa de tratamiento diferencial consta de una subetapa de cálculo de señales diferenciales en bruto a lo largo de la cual se determinan unas señales diferenciales en bruto Vbj (i, t) = Sj (i, t) -Sj-1 (i, t);
-la etapa de tratamiento diferencial consta de una subetapa de cálculo de señales diferenciales en bruto a lo largo de la cual se determinan unas señales diferenciales en bruto Vbj (i, t) al menos mediante un filtrado de paso alto de las señales en bruto Sj (i, t) en j;
-la etapa de tratamiento diferencial consta, además, de una subetapa de determinación de envolvente a lo largo de la cual se determinan las señales diferenciales Vj (i, t) calculando una envolvente temporal de cada señal diferencial en bruto Vbj (i, t).
-a lo largo de la etapa de tratamiento diferencial, las señales diferenciales en bruto experimentan un filtrado temporal de paso alto en t antes de la subetapa de determinación de la envolvente;
-la subetapa de determinación de la envolvente consta de un cálculo de una envolvente temporal Vej (i, t) y a continuación de un filtrado de paso bajo de las envolventes temporales Vej (i, t) en i, para obtener las señales diferenciales Vj (i, t);
-el campo de observación consta de unas microburbujas que forman dichos difusores y se detectan las microburbujas que han desaparecido del campo de observación de una etapa de medición a otra, formando estas microburbujas desaparecidas dichos blancos diferenciales;
-la onda incidente emitida presenta una amplitud adaptada para destruir el número máximo C de microburbujas en cada etapa de medición;
-la onda incidente emitida en cada etapa de medición presenta una amplitud adaptada para no destruir microburbujas y el procedimiento consta, además, de forma alterna con las etapas de medición, unas etapas de destrucción a lo largo de las cuales se emite una onda ultrasónica destructiva con una amplitud adaptada para destruir el número máximo C de microburbujas en cada etapa de destrucción;
-las posiciones Aj(x0, y0) de los blancos diferenciales sucesivos se añaden en una imagen del campo de observación (obtenida mediante un método de formación de imágenes seleccionado entre la ecografía, la radiología, la IRM, etc.);
-dicha imagen del campo de observación se obtiene mediante ecografía a partir de dicha red de transductores.
Por otra parte, la invención también tiene por objeto un dispositivo para la aplicación de un procedimiento como el que se ha definido con anterioridad, constando este dispositivo de una red de transductores controlada por un dispositivo de control y tratamiento adaptado para obtener imágenes de un campo de observación en un medio para formar imágenes que contiene unos difusores, estando el dispositivo de control y tratamiento adaptado para, a lo largo de varias etapas sucesivas de medición:
-hacer que la red de transductores emita una onda ultrasónica incidente en el campo de observación en cada etapa de medición;
-a continuación grabar las señales en bruto Sj (i, t) captadas por cada transductor y representativas de una onda ultrasónica refleja reverberada por los difusores del medio a partir de la onda incidente, siendo i un índice que designa cada sensor, siendo j un índice que designa cada etapa de medición y designando t el tiempo,
extendiéndose la red de transductores a lo largo de al menos una dimensión y propagándose las ondas incidentes principalmente en una dirección de propagación perpendicular a la red de transductores, caracterizado por que el dispositivo de control y tratamiento está adaptado para generar un número máximo C de blancos diferenciales, diferentes de una etapa de medición a otra, siendo cada blanco diferencial un difusor que está presente en el campo de observación durante una etapa de medición y ausente durante una etapa de medición inmediatamente contigua, y por que el dispositivo de control y tratamiento está, además, adaptado para realizar las siguientes etapas:
-una etapa de tratamiento diferencial a lo largo de la cual se comparan las señales en bruto Sj (i, t) correspondientes a unas etapas de medición sucesivas de índice j para extraer de estas unas señales diferenciales Vj(i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición;
15
25
35
45
55
65
E11799794
14-11-2014
-una etapa de ajuste a lo largo de la cual se determina al menos una función y = Pj(x) correspondiente a cada
señal diferencial Vj(i, t), donde x es una variable de espacio que designa una posición perpendicularmente a la
dirección de propagación e y es una coordenada que designa la posición de un punto a lo largo de la dirección
de propagación correspondiente a un tiempo de trayecto t; -y una etapa de localización a lo largo de la cual se determina un vértice Aj(x0, y0) de dicha función Pj, que
corresponde a la posición del blanco diferencial.
Descripción de las figuras
Se mostrarán otras características y ventajas de la invención a lo largo de la siguiente descripción de una de sus formas de realización, que se da a título de ejemplo no limitativo, en relación a los dibujos adjuntos.
En los dibujos:
-la figura 1 es una vista esquemática de un dispositivo de formación de imágenes por ultrasonido de acuerdo con
una forma de realización de la invención; -la figura 2 es un esquema que muestra un detalle del medio para formar imágenes; -la figura 3 es una gráfica que representa unas señales diferenciales Vj obtenidas mediante el procedimiento de
formación de imágenes de acuerdo con una forma de realización de la invención; -y la figura 4 es una gráfica similar a la figura 3, que muestra una curva parabólica Pj ajustada en las señales
diferenciales Vj y la posición del vértice Aj de esta curva, que corresponde a un blanco diferencial que está
situado en el medio para formar imágenes con una precisión muy inferior a la longitud de onda de las ondas
ultrasónicas utilizadas.
Descripción detallada de la invención
La figura 1 muestra un ejemplo de dispositivo 1 adaptado para generar una imagen de un medio 2 (por ejemplo una parte del cuerpo de un paciente) mediante la emisión y recepción de ondas ultrasónicas de compresión, en particular en un campo de observación 8. El medio 2 es difusor para las ondas ultrasónicas, es decir que es heterogéneo y/o contiene unos difusores distribuidos de manera aleatoria y que pueden reflejar las ondas ultrasónicas emitidas en el medio 2.
El dispositivo de formación de imágenes que se representa en la figura 1 consta de una red 3 de transductores Ti-Tn, por ejemplo de una red lineal que se extiende a lo largo de un eje X (o a lo largo de dos ejes en el caso de que la red 3 fuera de dos dimensiones), o llegado el caso de una red de dos dimensiones. Esta red 3 consta de un número n de transductores, pudiendo estar n, por ejemplo, comprendido entre 50 y 500, e incluso superior a 5.000 para una red de dos dimensiones. Se puede utilizar, por ejemplo, una barra lineal de entre 60 y 100 transductores cada uno con una anchura inferior, por ejemplo, a 1 mm a lo largo del eje X.
La red 3 de transductores está controlada por un dispositivo de control y tratamiento 4 (UC) por ejemplo una unidad informática central o varias unidades centrales. El dispositivo de control y tratamiento 4 está adaptado para controlar individualmente los transductores de la red 3 y para grabar y tratar las señales retrodifundidas por el medio 2 de tal modo que se forme una imagen del medio.
En el ejemplo considerado, el dispositivo de formación de imágenes 1 consta, además, de una pantalla 5 o de cualquier otra interfaz de usuario que permita visualizar una imagen del campo de observación tal como se representa en la figura 2, que muestra unos vasos capilares 6 del cuerpo humano 2, vasos por los cuales circulan unas microburbujas 7 o elementos similares, previamente inyectadas en el medio 2 (por ejemplo en el sistema sanguíneo) y que forman unos difusores capaces de retrodifundir las ondas ultrasónicas emitidas por la red 3 de transductores.
Estas microburbujas pueden, por ejemplo ser del tipo que describen Dayton y otros [Molecular ultrasound imaging using microbubble contrast agent, en Frontiers in Bioscience 12, págs. 5.124-5.142 -septiembre de 2007].
El dispositivo de control y tratamiento 4 está adaptado (programado) para aplicar el procedimiento de formación de imágenes que se describe a continuación, que consta de:
(a)
una multitud de etapas de medición, a lo largo de las cuales el dispositivo de control y tratamiento 4 graba las señales en bruto; (a’) estando estas etapas de medición eventualmente separadas por unas etapas de destrucción de microburbujas;
(b)
al menos una etapa de tratamiento diferencial de las señales en bruto, para obtener las señales diferenciales representativas de blancos diferenciales formados en este caso mediante la destrucción de microburbujas;
(c)
al menos una etapa de ajuste a lo largo de la cual se ajusta una curva parabólica en cada señal diferencial;
(d)
y al menos una etapa de localización a lo largo de la cual se determina un vértice de cada curva parabólica,
15
25
35
45
55
65
E11799794
14-11-2014
que corresponde a la posición del blanco diferencial.
(a) Etapa de medición:
El procedimiento consta de varias etapas sucesivas de medición, un total de N, por ejemplo aproximadamente 400 etapas de medición repetidas aproximadamente cada ms o menos.
A lo largo de cada una de estas etapas de medición, se hace que la red 3 de transductores emita una onda ultrasónica incidente en el campo de observación 8, por ejemplo una onda plana o eventualmente una onda focalizada, a continuación se graban las señales en bruto Sj (i, t) captadas por cada transductor T1-Tn y representativas de una onda ultrasónica refleja reverberada por los difusores del medio a partir de la onda incidente, siendo i un índice que designa cada transductor, siendo j un índice que designa cada etapa de medición y designando t el tiempo.
Hay que señalar que las señales en bruto no son unas imágenes ecográficas, la realización de una imagen ecográfica precisa un tratamiento de las señales en bruto, por ejemplo mediante la formación de trazas.
Sin embargo, al menos una de las etapas de medición puede servir para establecer una imagen ecográfica del campo de observación 8 mediante el dispositivo de control y tratamiento 4, con una resolución milimétrica clásica, superponiéndose a continuación esta imagen a las posiciones de los blancos diferenciales obtenidos en la etapa (e) que se describirá a continuación. Eventualmente, se pueden realizar unas imágenes ecográficas mediante el dispositivo de control y tratamiento 4 en varias etapas de medición o en cada etapa de medición, para tener en cuenta eventuales deformaciones del medio 2 a lo largo del desarrollo del procedimiento.
De acuerdo con la invención, de una etapa de medición a otra, se destruye un reducido número C de microburbujas. C es un número máximo de microburbujas destruidas en el campo de observación 8 de una etapa de medición a otra. Cada microburbuja destruida se llamará aquí blanco diferencial. De manera más general, los blancos diferenciales pueden ser los difusores que desaparecen o aparecen de una etapa de medición a otra, es decir los difusores presentes en el campo de observación durante una etapa de medición y ausentes durante una etapa de medición inmediatamente contigua.
El número C es por lo general como máximo igual a 2, de preferencia como máximo igual a 1, De manera más general, el número C puede ser como máximo igual a INT(A/(5λ)2)+1, donde A es un área del campo de observación e INT designa la función parte entera.
Destrucción de microburbujas:
Para obtener esta destrucción de microburbujas, el dispositivo de control y tratamiento 4 hace que la red 3 de transductores emita unas ondas ultrasónicas (por ejemplo unas ondas planas o eventualmente unas ondas focalizadas) con una amplitud justa suficiente para destruir el número máximo C antes mencionado de microburbujas. Esta amplitud se puede determinar empíricamente durante una etapa inicial de ajuste (en este caso se determina mediante ecografía la cantidad de microburbujas destruidas en cada disparo de onda ultrasónica), o bien se puede determinar (empíricamente o mediante su cálculo) por adelantado en función de las condiciones de funcionamiento.
Estas ondas ultrasónicas que sirven para destruir una o varias microburbuja(s) en cada disparo pueden ser:
-o bien las ondas incidentes emitidas en cada etapa de medición, en cuyo caso el procedimiento no consta de
etapa específica de destrucción de microburbujas; -o bien unas ondas ultrasónicas destructivas con una amplitud superior a las ondas incidentes de las etapas de
medición, emitiéndose cada una de dichas ondas destructivas a lo largo de una etapa de destrucción de
microburbujas (a’), intercalada entre dos etapas de medición (a).
(b) Etapa de tratamiento diferencial:
Tras la grabación de las señales en bruto Sj (i, t), se procede a una etapa de tratamiento diferencial, que se puede iniciar tras la grabación completa de todas las señales en bruto o simplemente tras la grabación de las primeras señales en bruto, desarrollándose entonces parcialmente mientras se graban las otras señales en bruto.
A lo largo de esta etapa de tratamiento diferencial, el dispositivo de control y tratamiento 4 compara las señales en bruto Sj (i, t) correspondientes a unas etapas de medición sucesivas de índice j para extraer de estas unas señales diferenciales Vj (i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición. Estas señales, que corresponden a eventos puntuales aislados, son por lo general de tipo impulsional en cada canal de recepción, es decir en cada transductor.
Esta etapa de tratamiento diferencial puede a su vez constar de varias subetapas, que se encadenan, por ejemplo,
15
25
35
45
55
65
E11799794
14-11-2014
de la siguiente forma:
(b1) Subetapa de cálculo de señales diferenciales en bruto:
A lo largo de esta subetapa, el dispositivo de control y tratamiento 4 calcula las señales diferenciales en bruto Vbj (i, t).
Estas señales diferenciales en bruto Vbj (i, t) se pueden obtener, por ejemplo, de dos maneras a partir de las señales en bruto Sj (i, t):
-bien mediante el cálculo de diferencias: Vbj (i, t) = Sj (i, t) -Sj-1 (i, t), para j = 2 en N;
-o bien mediante un filtrado de paso alto de las señales en bruto Sj (i, t), realizándose el filtrado en j para conservar únicamente los cambios rápidos entre las señales en bruto obtenidas en las diferentes etapas de medición j.
(b2) Filtrado temporal de paso bajo en t:
A continuación el dispositivo de control y tratamiento 4 puede eventualmente proceder a una subetapa de filtrado temporal de paso bajo en t de las señales diferenciales en bruto Vbj (i, t).
(b3) Determinación de la envolvente:
A continuación el dispositivo de control y tratamiento 4 puede proceder a una subetapa de determinación de la envolvente a lo largo de la cual se determinan las señales diferenciales Vj (i, t):
-calculando una envolvente temporal Vej (i, t) de las señales diferenciales en bruto Vbj (i, t) filtradas resultantes de la subetapa (b2); -a continuación procediendo a un filtrado de paso bajo en i de las envolventes temporales Vej (i, t) para obtener las señales diferenciales Vj (i, t).
(c) Etapa de ajuste:
A continuación el dispositivo de control y tratamiento 4 procede a una etapa de ajuste a lo largo de la cual se determina al menos una función y = Pj(x) correspondiente a cada señal diferencial Vj (i, t).
Para ello, el dispositivo de control y tratamiento 4 determina en primer lugar unos valores máximos Vj (i, t) de Vj (i, t) para cada transductor i, y para cada etapa de medición j > 1.
A continuación el dispositivo de control y tratamiento 4 deduce de estos los puntos correspondientes Dj (xi, yi) del campo de observación 8, donde:
-xi es una variable de espacio que designa una posición de cada transductor i perpendicularmente a la dirección de propagación (es decir una variable unidimensional que consiste simplemente en una coordenada x a lo largo del eje X en el caso considerado aquí, o llegado el caso una variable bidimensional que comprende dos coordenadas en dos ejes en el caso de una red 3 de dos dimensiones);
-e yi es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto ti para el cual Vj (i, t) es máximo (yi = c.ti/2, donde c es la celeridad de la onda ultrasónica).
Hay que señalar que, de manera más general, estos puntos Dj (xi, yi) se podrían determinar con unos valores yi que designan la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto ti característico de la señal Vj (i, t): este tiempo puede, por ejemplo, corresponder al máximo de la señal Vj (i, t) como se ha descrito con anterioridad, o corresponder al frente de onda de la señal Vj (i, t), u otro.
A continuación, el dispositivo de control y tratamiento 4 determina para etapa de medición j > 1, una función y = Pj (x) que se ajusta para minimizar las diferencias con los puntos Dj (xi, yi) (por ejemplo, se minimiza la desviación cuadrática Σ (Pj (xi)-yi)2). La función y = Pj (x) es una función continua en x, de preferencia una función polinomial de segundo grado (parabólica). Esta función corresponde a una parábola para una red 3 lineal tal como se describe aquí, o un paraboloide para una red 3 bidimensional.
(d) Etapa de localización:
Por último, el dispositivo de control y tratamiento 4 procede a una etapa de localización del blanco diferencial de cada etapa de medición j > 1, es decir para simplificar, una localización de la microburbuja que ha desaparecido entre las etapas de medición j-1 y j. Para ello, el dispositivo de control y tratamiento 4 determina el vértice Aj (x0, y0) de la función Pj mencionada con anterioridad, y este vértice corresponde a la posición del blanco diferencial,
E11799794
14-11-2014
determinada con una resolución claramente inferior a la longitud de onda de las ondas ultrasónicas, por ejemplo del orden de entre λ/50 y λ/200.
A continuación las posiciones Aj (x0, y0) de los blancos diferenciales sucesivos pueden añadirse de manera ventajosa
5 en la imagen del campo de observación 8, obtenida por ejemplo mediante ecografía como se ha indicado con anterioridad (o en cualquier otra imagen del campo de observación, obtenida por ejemplo mediante radiografía, IRM u otro método), de tal modo que se obtenga una imagen de muy alta resolución del campo de observación 8.
Las etapas (c) y (d) mencionadas con anterioridad se pueden realizar una vez que se han determinado todas las 10 señales diferenciales, o a medida que se calculan estas señales diferenciales Vj (i, t).
Por otra parte, hay que señalar que los blancos diferenciales mencionados con anterioridad pueden corresponder a la creación de microburbujas o de otros difusores, por ejemplo burbujas de cavitación obtenidas entre las etapas de medición mediante la focalización de ultrasonidos o por efecto fotoacústico, utilizando en particular la red 3 de
15 transductores.

Claims (14)

  1. 5
    15
    25
    35
    45
    55
    65
    E11799794
    14-11-2014
    REIVINDICACIONES
    1. Procedimiento de formación de imágenes por ultrasonido para obtener imágenes de un campo de observación (8) en un medio para formar imágenes (2) que contiene unos difusores (7), constando dicho procedimiento de varias etapas sucesivas de medición a lo largo de cada una de las cuales:
    -se hace que una red (3) de transductores (T1-Tn) emita una onda ultrasónica incidente en el campo de observación (8), con una longitud de onda central λ; -a continuación se graban las señales en bruto Sj (i, t) captadas por cada transductor (T1-Tn) y representativas de una onda ultrasónica refleja reverberada por los difusores (7) del medio a partir de la onda incidente, siendo i un índice que designa cada transductor, siendo j un índice que designa cada etapa de medición y designando t el tiempo, extendiéndose la red (3) de transductores a lo largo de al menos una dimensión y propagándose principalmente las ondas incidentes en una dirección de propagación perpendicular a la red de transductores,
    caracterizado porque se genera un número máximo C de blancos diferenciales, diferentes de una etapa de medición a otra, siendo cada blanco diferencial un difusor (7) que está presente en el campo de observación durante una etapa de medición y ausente durante una etapa de medición inmediatamente contigua, siendo el número C como máximo igual a INT (A/(5λ)2)+1, donde A es un área del campo de observación, y porque dicho procedimiento consta, además, de las siguientes etapas:
    -una etapa de tratamiento diferencial a lo largo de la cual se comparan las señales en bruto Sj (i, t) correspondientes a las etapas de medición sucesivas de índice j para extraer de estas unas señales diferenciales Vj (i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición; -una etapa de ajuste a lo largo de la cual se determina al menos una función y = Pj (x) correspondiente a cada señal diferencial Vj (i, t), donde x es una variable de espacio que designa una posición perpendicularmente a la dirección de propagación e y es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto t; -y una etapa de localización a lo largo de la cual se determina un vértice Aj (x0, y0) de dicha función Pj, que corresponde a la posición del blanco diferencial.
  2. 2.
    Procedimiento de acuerdo con la reivindicación 1, en el que el número C es como máximo igual a 2, de preferencia igual a 1.
  3. 3.
    Procedimiento de acuerdo con la reivindicación 1 o la reivindicación 2, en el que, a lo largo de la etapa de ajuste, se determina la función y = Pj (x) ajustando dicha función para minimizar las diferencias con los puntos Dj (xi, yi), donde xi es una variable de espacio que designa una posición de cada transductor i perpendicularmente a la dirección de propagación e yi es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto ti característico de la señal Vj (i, t).
  4. 4.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones anteriores, en el que la función P es parabólica.
  5. 5.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones anteriores, en el que la etapa de tratamiento diferencial consta de una subetapa de cálculo de señales diferenciales en bruto a lo largo de la cual se determinan las señales diferenciales en bruto Vbj (i, t) = Sj (i, t) -Sj-1 (i, t).
  6. 6.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones 1 a 4, en el que la etapa de tratamiento diferencial consta de una subetapa de cálculo de señales diferenciales en bruto a lo largo de la cual se determinan las señales diferenciales en bruto Vbj (i, t) al menos mediante un filtrado de paso alto de las señales en bruto Sj (i, t) en j.
  7. 7.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones 5 y 6, en el que la etapa de tratamiento diferencial consta, además, de una subetapa de determinación de la envolvente a lo largo de la cual se determinan las señales diferenciales Vj (i, t) calculando una envolvente temporal de cada señal diferencial en bruto Vbj (i, t).
  8. 8.
    Procedimiento de acuerdo con la reivindicación 7, en el que la subetapa de determinación de la envolvente consta de un cálculo de una envolvente temporal Vej (i, t) y a continuación de un filtrado de paso bajo de las envolventes temporales Vej (i, t) en i, para obtener las señales diferenciales Vj (i, t).
  9. 9.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones anteriores, en el que el campo de observación (8) consta de unas microburbujas (7) y se detectan las microburbujas que han desaparecido del campo de observación de una etapa de medición a otra, formando estas microburbujas desaparecidas dichos blancos diferenciales.
  10. 10.
    Procedimiento de acuerdo con la reivindicación 9, en el que la onda incidente presenta una amplitud adaptada para destruir el número máximo C de microburbujas (7) en cada etapa de medición.
    8 5
    10
    15
    20
    25
    30
    35
    40
    E11799794
    14-11-2014
  11. 11.
    Procedimiento de acuerdo con la reivindicación 9, en el que la onda incidente emitida en cada etapa de medición presenta una amplitud adaptada para no destruir microburbujas (7) y el procedimiento consta, además, de forma alterna con las etapas de medición, de unas etapas de destrucción a lo largo de las cuales se emite una onda ultrasónica destructiva que tiene una amplitud adaptada para destruir el número máximo C de microburbujas (7) en cada etapa de destrucción.
  12. 12.
    Procedimiento de acuerdo con una cualquiera de las reivindicaciones anteriores, en el que las posiciones Aj (x0, y0) de los blancos diferenciales sucesivos se añaden en una imagen del campo de observación (8).
  13. 13.
    Procedimiento de acuerdo con la reivindicación 12, en el que dicha imagen del campo de observación (8) se obtiene mediante ecografía a partir de dicha red (3) de transductores.
  14. 14.
    Dispositivo para la aplicación de un procedimiento de acuerdo con una cualquiera de las reivindicaciones anteriores, constando este dispositivo de una red (3) de transductores (T1-Tn) controlada por un dispositivo de control y tratamiento (4) adaptado para obtener imágenes de un campo de observación (8) en un medio para formar imágenes (2) que contiene unos difusores (7), estando el dispositivo de control y tratamiento (4) adaptado para, a lo largo de varias etapas sucesivas de medición:
    -hacer que la red de transductores (3) emita una onda ultrasónica incidente en el campo de observación (8) en cada etapa de medición; -a continuación grabar las señales en bruto Sj (i, t) captadas por cada transductor (T1-Tn) y representativas de una onda ultrasónica refleja reverberada por los difusores (7) del medio a partir de la onda incidente, siendo i un índice que designa cada sensor, siendo j un índice que designa cada etapa de medición y designando t el tiempo, extendiéndose la red (3) de transductores a lo largo de al menos una dimensión y propagándose las ondas incidentes principalmente en una dirección de propagación perpendicular a la red de transductores,
    caracterizado porque el dispositivo de control y tratamiento (4) está adaptado para generar un número máximo C de blancos diferenciales, diferentes de una etapa de medición a otra, siendo cada blanco diferencial un difusor (7) que está presente en el campo de observación durante una etapa de medición y ausente durante una etapa de medición inmediatamente contigua, siendo el número C como máximo igual a INT (A/(5λ)2)+1, donde A es un área del campo de observación, y porque el dispositivo de control y tratamiento (4) está, además, adaptado para realizar las siguientes etapas:
    -una etapa de tratamiento diferencial a lo largo de la cual se comparan las señales en bruto Sj (i, t) correspondientes a las etapas de medición sucesivas de índice j para extraer de estas unas señales diferenciales Vj (i, t) representativas de variaciones entre señales en bruto de las etapas sucesivas de medición; -una etapa de ajuste a lo largo de la cual se determina al menos una función y = Pj (x) correspondiente a cada señal diferencial Vj (i, t), donde x es una variable de espacio que designa una posición perpendicularmente a la dirección de propagación e y es una coordenada que designa la posición de un punto a lo largo de la dirección de propagación correspondiente a un tiempo de trayecto t; -y una etapa de localización a lo largo de la cual se determina un vértice Aj (x0, y0) de dicha función Pj, que corresponde a la posición del blanco diferencial.
    9
ES11799794.0T 2010-12-16 2011-11-29 Procedimiento y dispositivo de formación de imágenes por ultrasonido Active ES2523982T3 (es)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1060633 2010-12-16
FR1060633A FR2969350B1 (fr) 2010-12-16 2010-12-16 Procede et dispositif d'imagerie sonore.
PCT/FR2011/052810 WO2012080614A1 (fr) 2010-12-16 2011-11-29 Procede et dispositif d'imagerie ultrasonore

Publications (1)

Publication Number Publication Date
ES2523982T3 true ES2523982T3 (es) 2014-12-03

Family

ID=44201042

Family Applications (1)

Application Number Title Priority Date Filing Date
ES11799794.0T Active ES2523982T3 (es) 2010-12-16 2011-11-29 Procedimiento y dispositivo de formación de imágenes por ultrasonido

Country Status (7)

Country Link
US (1) US9329260B2 (es)
EP (1) EP2652517B1 (es)
JP (1) JP5947312B2 (es)
ES (1) ES2523982T3 (es)
FR (1) FR2969350B1 (es)
IL (1) IL226916A (es)
WO (1) WO2012080614A1 (es)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10751028B2 (en) 2016-04-01 2020-08-25 The Board Of Trustees Of The Leland Stanford Junior University Coherence-based beamforming for improved microbubble detection in contrast enhanced ultrasound
WO2018050817A1 (en) 2016-09-16 2018-03-22 Institut National De La Sante Et De La Recherche Medicale (Inserm) Method for imaging a sample with blood and associated devices
EP3600065A1 (en) 2017-03-22 2020-02-05 Institut National de la Sante et de la Recherche Medicale (INSERM) Method for imaging an area of a medium with ultrasound contrast agents and associated device
JP7606662B2 (ja) 2019-01-03 2024-12-26 コーニンクレッカ フィリップス エヌ ヴェ コントラスト強調画像のためのシステム及び方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7841982B2 (en) * 1995-06-22 2010-11-30 Techniscan, Inc. Apparatus and method for imaging objects with wavefields
JP3862838B2 (ja) * 1997-11-26 2006-12-27 株式会社東芝 超音波診断装置
US6398735B1 (en) * 2000-03-07 2002-06-04 Koninklijke Philips Electronics N.V. Detecting a relative level of an ultrasound imaging contrast agent
ATE466596T1 (de) * 2004-01-20 2010-05-15 Sunnybrook & Womens College Hochfrequenz-ultraschall-darstellung mit kontrastmitteln
JP4521204B2 (ja) * 2004-02-27 2010-08-11 株式会社東芝 超音波診断装置、画像処理装置及び超音波画像撮影方法
US7850611B2 (en) * 2004-09-20 2010-12-14 Innervision Medical Technologies Inc. System and methods for improved ultrasound imaging
US7862508B2 (en) * 2004-09-20 2011-01-04 Innervision Medical Technologies Inc. Systems and methods for ultrasound imaging
US9439571B2 (en) * 2006-01-20 2016-09-13 Washington University Photoacoustic and thermoacoustic tomography for breast imaging
US20070287912A1 (en) * 2006-05-31 2007-12-13 Khuri-Yakub Butrus T Functional imaging using capacitive micromachined ultrasonic transducers
US20080312581A1 (en) * 2007-06-06 2008-12-18 Biovaluation & Analysis, Inc. Peptosomes for Use in Acoustically Mediated Intracellular Drug Delivery in vivo
US8306293B2 (en) * 2008-05-15 2012-11-06 University Of Virginia Patent Foundation Reduction of echo decorrelation facilitating motion estimation

Also Published As

Publication number Publication date
JP5947312B2 (ja) 2016-07-06
FR2969350A1 (fr) 2012-06-22
IL226916A (en) 2017-01-31
FR2969350B1 (fr) 2013-01-11
US9329260B2 (en) 2016-05-03
JP2013545568A (ja) 2013-12-26
EP2652517A1 (fr) 2013-10-23
US20130301382A1 (en) 2013-11-14
WO2012080614A1 (fr) 2012-06-21
EP2652517B1 (fr) 2014-10-08

Similar Documents

Publication Publication Date Title
ES2631914T3 (es) Mapeado y caracterización de actividad de cavitación
US9743909B1 (en) Imaging bubbles in a medium
US6334846B1 (en) Ultrasound therapeutic apparatus
US6454713B1 (en) Ultrasound therapeutic apparatus
EP1451593B1 (en) System and method for sensing and locating disturbances in an energy path of a focused ultrasound system
US7587291B1 (en) Focusing of broadband acoustic signals using time-reversed acoustics
ES2325919T3 (es) Procedimiento y dispositivo de formacion de imagenes que usa ondas de cizalla.
KR101337250B1 (ko) 초음파 진단 및 치료 장치
JP5905080B2 (ja) オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成
CN107205720B (zh) 一种超声适应性波束形成方法及其对经颅成像的应用
US20150320384A1 (en) Methods and Systems for Detecting an Object in a Subject with Ultrasound
ES2523982T3 (es) Procedimiento y dispositivo de formación de imágenes por ultrasonido
JP2009500082A (ja) 凹状フェーズドアレイ・撮像カテーテル
Ding et al. Ultrasound line-by-line scanning method of spatial–temporal active cavitation mapping for high-intensity focused ultrasound
Bawiec et al. Inertial cavitation behaviors induced by nonlinear focused ultrasound pulses
CN112023283B (zh) 基于高阶孔径自相关的环形多阵列超声被动成像方法及系统
US20210113192A1 (en) Shear wave amplitude reconstruction for tissue elasticity monitoring and display
US20110208095A1 (en) Method for modifying glomerular permeability and function with focused ultrasound
Yeats et al. Soft tissue aberration correction for histotripsy using acoustic emissions from cavitation cloud nucleation and collapse
Song et al. Quantitative assessment of boiling histotripsy progression based on color Doppler measurements
US20220175406A1 (en) SYSTEM AND METHOD FOR DETECTING AND ALIGNING ACOUSTIC BEAM IN SITU TO A TARGET USING WIDE-BEAM, LOW FREQUENCY (<1 MHz) ULTRASOUND
Zhou et al. Acoustic beam mapping for guiding HIFU therapy in vivo using sub-therapeutic sound pulse and passive beamforming
WO2018037130A1 (en) Method and system for localizing a region of interest in a medium in which cavitation occurs
Jeong et al. A novel approach for the detection of every significant collapsing bubble in passive cavitation imaging
ES2982024T3 (es) Procedimiento y aparato para formación de imágenes en tiempo real de la propagación de una onda mecánica en un material de propagación acústica