ES2685977T3 - Método y sistema para analizar datos de imágenes - Google Patents

Método y sistema para analizar datos de imágenes Download PDF

Info

Publication number
ES2685977T3
ES2685977T3 ES14184117.1T ES14184117T ES2685977T3 ES 2685977 T3 ES2685977 T3 ES 2685977T3 ES 14184117 T ES14184117 T ES 14184117T ES 2685977 T3 ES2685977 T3 ES 2685977T3
Authority
ES
Spain
Prior art keywords
image
segmentation
outliers
elements
processor
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
ES14184117.1T
Other languages
English (en)
Inventor
Saurabh Jain
Dirk Smeets
Diana Sima
Annemie Ribbens
Anke Maertens
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.)
Icometrix Nv
Original Assignee
Icometrix Nv
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 Icometrix Nv filed Critical Icometrix Nv
Application granted granted Critical
Publication of ES2685977T3 publication Critical patent/ES2685977T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2431Multiple classes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/143Segmentation; Edge detection involving probabilistic approaches, e.g. Markov random field [MRF] modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/174Segmentation; Edge detection involving the use of two or more images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10084Hybrid tomography; Concurrent acquisition with multiple different tomographic modalities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20128Atlas-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Data Mining & Analysis (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Geometry (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Un método de análisis de datos de imágenes, que comprende: - obtener (1) una primera imagen (10) de una primera parte de un objeto; - obtener (2) una segunda imagen (20) de una segunda parte del objeto, teniendo dicha segunda parte una superposición sustancial con dicha primera parte; - obtener un mapeo (110) entre dicha primera imagen (10) y dicha segunda imagen (20); - segmentar (200) dicha segunda imagen (20) para obtener una segmentación de tejidos; - detectar valores atípicos (300) en dicha primera imagen (10) identificando los valores de intensidad extrema de elementos dentro de una o más clases de elementos basándose en dicha segmentación; - reemplazar (500) los elementos de dicha segunda imagen (20) que corresponden, de acuerdo con dicho mapeo, a al menos algunos de dichos valores atípicos detectados de dicha primera imagen (10), con valores de reemplazo, para obtener una segunda imagen corregida, siendo los al menos algunos de dichos valores atípicos detectados de dicha primera imagen, una segmentación de la lesión; y - actualizar dicha segmentación realizando dicha segmentación (200) en dicha segunda imagen corregida; en el que dicha detección de los valores atípicos (300), dicho reemplazo (500) y dicha actualización (200) se realizan sucesivamente hasta que se cumpla un criterio de convergencia predeterminado (550), representando dicho criterio de convergencia predeterminado un punto en el cual no existen cambios significativos en el tejido y en la segmentación de la lesión.

Description

5
10
15
20
25
30
35
40
45
50
DESCRIPCION
Método y sistema para analizar datos de imágenes Campo de la invención
La presente invención pertenece al campo del análisis de la imagen digital, en particular al análisis de imágenes médicas.
Antecedentes
En muchas disciplinas médicas, el diagnóstico y el seguimiento de enfermedades se basa cada vez más en la interpretación de los datos obtenidos mediante técnicas de imágenes médicas. En el campo de la neurología, la imagen por resonancia magnética (MR) se ha vuelto particularmente importante, ya que tiene ventajas específicas sobre la tomografía computarizada (CT) basada en rayos X. en determinadas situaciones.
En particular, el volumen (estimado) de lesiones de materia blanca se usa como un bio marcador para algunas enfermedades neurológicas, en particular la esclerosis múltiple (MS). En los últimos años, se ha intentado automatizar la detección de tales lesiones basándose en imágenes MR.
En el artículo Automated Segmentation of Múltiple Sclerosis Lesions by Model Outlier Detection (IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL 20, N° 8, AGOSTO 2001), K. VAN LEEMPUT et al., se presenta un algoritmo para la segmentación de lesiones de esclerosis múltiple a partir de imágenes de resonancia magnética multi espectral. En el método descrito en ese artículo, las lesiones MS se detectan como valores atípicos con respecto a un modelo estadístico para intensidades del tejido cerebral normal en imágenes MR.
En el artículo An Automatic Segmentation of T2-FLAIR Multiple Sclerosis Lesions (The MIDAS Journal - MS Lesion Segmentation (MICCAI 2008 Workshop)), JC. SOUPLET et al., se presenta un método diseñado para detectar un área de señal hiper intensa en una secuencia T2-FLAIR. El algoritmo descrito utiliza tres secuencias MRI convencionales: T1, T2 y T2-FLAIR. En primer lugar, las imágenes están recortadas, no están polarizadas espacialmente y están desprovistas del cráneo. Se realiza una segmentación del cerebro en sus diferentes compartimentos en las secuencias T1 y T2. A partir de estas segmentaciones, se calcula automáticamente un umbral para la secuencia T2-FLAIR. Luego, las operaciones post procesamiento seleccionan las lesiones más plausibles en las hiper intensas señales obtenidas.
En el artículo Evaluating and Reducing the Impact of White Matter Lesions on Brain Volume Measurements (Human Brain Mapping, Volumen 33, Número 9, 2011), M. BATTAGLINI et al., se describe cómo la presencia de lesiones en la materia blanca afecta a ciertas mediciones del volumen cerebral basadas en la segmentación. El artículo indica que rellenar las lesiones con intensidades que coinciden con la materia blanca circundante de apariencia normal aseguraba mediciones precisas de la clase del tejido y por lo tanto representa un enfoque prometedor para la precisión de las mediciones de la clasificación del tejido y del volumen cerebral.
El artículo "White matter lesion extension to automatic brain tissue segmentation on MRI", DE BOER R ET AL., 2009, describe un método para segmentar lesiones cerebrales basándose en las imágenes T1 y FLAIR co- registradas. La imagen T1 primero se segmenta en materia blanca (WM), materia gris (GM) y líquido cefalorraquídeo (CSF). El resultado de esta segmentación se transfiere a la imagen FLAIR para detectar valores atípicos por detección del umbral.
Existe la necesidad de una forma automatizada para estimar con mayor precisión y fiabilidad los volúmenes respectivos de materia blanca, materia gris y/o líquido cefalorraquídeo (CSF), en presencia de lesiones, en particular lesiones de la materia blanca.
Resumen de la invención
De acuerdo con un aspecto de la presente invención, se proporciona un método para analizar datos de imágenes, de acuerdo con la reivindicación 1.
Una ventaja del método según la presente invención es que combina prudentemente información de diferentes tipos de imágenes para extraer información adicional que sería extremadamente difícil de discernir para un observador humano. La invención se basa, entre otras cosas, en la idea de los inventores de que diferentes tipos de imágenes pueden ser más adecuadas para la segmentación y la detección del valor atípico, respectivamente. La invención se basa además en una idea de los inventores de que diferentes tipos respectivos de informaciones obtenidas de ambos tipos de imágenes puede convenientemente combinarse si se dispone de un mapeo entre dos imágenes, lo que da lugar a una mejor segmentación y a una mejor detección de los valores atípicos. La invención se basa además en una idea de los inventores de que esta combinación de información puede realizarse sucesivamente para
obtener el equilibrio deseado entre el gasto de tiempo y la precisión de los resultados.
5
10
15
20
25
30
35
40
45
50
La primera imagen y la segunda imagen pueden ser bidimensionales o tridimensionales, por lo que los elementos de las imágenes son píxeles o vóxeles, respectivamente. El objeto de las imágenes puede comprender tejidos u órganos biológicos, en particular tejidos del sistema nervioso central, más en particular tejido cerebral. Las imágenes pueden haber sido adquiridas por medio de una técnica de imágenes médicas, tal como la resonancia magnética.
En una realización del método de acuerdo con la presente invención, los valores atípicos se detectan segmentando la primera imagen a partir de la segmentación de la segunda imagen con una clase de valor atípico adicional.
Una ventaja de esta realización es que los valores atípicos se detectan eficientemente por el mismo proceso general que proporciona la segmentación.
En una realización, el método de acuerdo con la presente invención comprende además iniciar la segmentación transfiriendo una segmentación de un ejemplo a la segunda imagen.
Una ventaja de esta realización es que la etapa de segmentación puede llevarse a cabo en un tiempo más corto usando una segmentación a modo de ejemplo, por ejemplo, un atlas.
En una realización en particular, la obtención del mapeo comprende realizar un registro de la primera imagen obtenida a la segunda imagen obtenida.
En los casos en que el mapeo no esté disponible de antemano, se puede obtener como parte de un método de acuerdo con la presente invención por medio de una etapa de registro.
En una realización del método de acuerdo con la presente invención, el objeto comprende al menos una parte de un cerebro, la primera imagen es una imagen FLAIR, la segunda imagen es una imagen T1; la segmentación comprende una clasificación de elementos de la imagen T1 como materia gris, materia blanca o líquido cefalorraquídeo; y los valores atípicos se detectan entre los elementos de la primera imagen que se clasifican como materia blanca.
Esta realización aplica los principios de la invención a imágenes MRI neurológicas, para obtener diagnóstica o terapéuticamente resultados de segmentación relevantes.
Una ventaja adicional de esta realización es que el tiempo requerido para el análisis de imagen con el método automatizado de acuerdo con la presente invención es mucho más corto que el tiempo que se requeriría para el análisis manual por expertos humanos. Por lo tanto, el método de la presente invención es particularmente adecuado para procesar la gran cantidad de imágenes implicadas en ensayos clínicos.
En una realización en particular, los valores de reemplazo se basan en los valores promedio de elementos clasificados como materia blanca en elementos no atípicos en la proximidad de dichos respectivos valores atípicos.
Los valores de reemplazo pueden ser los valores promedio reales de los elementos clasificados como materia blanca en elementos no atípicos en proximidad de dichos valores atípicos respectivos, o los valores derivados de estos promedios, por ejemplo mediante la adición de un término de ruido. Los inventores han encontrado que el uso de valores promedio, preferiblemente con ruido añadido, como valores de reemplazo, producen muy buenos resultados.
En una realización, el método de acuerdo con la presente invención comprende adicionalmente eliminar elementos que representan tejido no cerebral de la segunda imagen.
Una ventaja de esta realización es que las etapas posteriores del procesamiento de la imagen transcurren de una manera más eficiente si las imágenes no contienen ningún elemento que no sea relevante para el resultado final deseado.
En una realización del método de acuerdo con la presente invención, la segmentación se usa para calcular los respectivos volúmenes o áreas de materia gris, matera blanca y líquido cefalorraquídeo.
Estos cálculos del volumen o área son útiles para detectar anomalías o cambios en la distribución de estos tipos de materia en el cerebro.
En una realización del método de acuerdo con la presente invención, la detección de valores atípicos se usa para calcular volúmenes o áreas de valores atípicos.
Este cálculo del volumen o área es útil para detectar anomalías sospechosas tales como lesiones.
En una realización, el método de acuerdo con la presente invención comprende adicionalmente el cálculo de volúmenes o áreas en diferentes regiones anatómicas.
Una ventaja de esta realización es que pueden obtenerse estimaciones del volumen o área médicamente relevantes
para diferentes regiones anatómicas, cuyas anomalías pueden corresponder a diferentes condiciones médicas. Las diferentes regiones anatómicas se pueden definir de acuerdo con las divisiones que son de uso general en el campo de neurología, como por ejemplo la división del cerebro en los hemisferios izquierdo y derecho, y en un lóbulo frontal, un lóbulo parietal, un lóbulo temporal y un lóbulo occipital. Las diferentes regiones anatómicas se pueden 5 definir de acuerdo con una división que se usa preferentemente en el diagnóstico y seguimiento de un trastorno específico, tal como la esclerosis múltiple, para el que una división de ejemplo puede consistir en una región yuxtacortical, una región peri ventricular y una fosa posterior. El experto apreciará que se pueden utilizar otras divisiones en regiones anatómicas de acuerdo con el tipo de información deseada.
De acuerdo con un aspecto de la presente invención, se proporciona un producto de programa de ordenador que 10 comprende medios de código configurado para hacer que un procesador lleve a cabo el método descrito anteriormente.
De acuerdo con un aspecto de la presente invención, se proporciona un sistema de procesamiento de imágenes que comprende: un interfaz de entrada adaptado para recibir una primera imagen de una primera parte de un objeto y una segunda imagen de una segunda parte del objeto, estando la segunda parte 15 superpuesta sustancialmente con la primera parte; un procesador, conectado operativamente al interfaz de entrada; y un interfaz de salida, conectado operativamente al procesador y adaptado para presentar los resultados producidos por el procesador; y una memoria, operativamente conectada al procesador y dispuesta para interpretar el código a ejecutar por el procesador, comprendiendo dicha memoria el código configurado para hacer que el procesador lleve a cabo el método de cualquiera de las reivindicaciones adjuntas usando la primera imagen y 20 la segunda imagen como entradas.
Los efectos técnicos y las ventajas de las realizaciones del producto del programa de ordenador y del sistema de acuerdo con la presente invención corresponden mutatis mutandis a las de las realizaciones correspondientes del método de acuerdo con la presente invención.
Breve descripción de las figuras
25 Estos y otros aspectos técnicos y ventajas de las realizaciones de la presente invención se describirán ahora con más detalle con referencia a los dibujos adjuntos, en los que:
La figura 1 proporciona un diagrama de flujo de un método de acuerdo con una realización de la presente invención;
La figura 2 proporciona un diagrama de flujo de un método de acuerdo con otra realización de la presente invención;
La figura 3 proporciona un diagrama de flujo más detallado del método representado en la figura 2; y
30 La figura 4 ilustra esquemáticamente un sistema de acuerdo con una realización de la presente invención. Descripción detallada de las realizaciones
La figura 1 proporciona un diagrama de flujo de un método para analizar datos de imágenes de acuerdo con una realización general de la presente invención.
El método ilustrado comprende obtener 1°, una primera imagen 10 de una primera parte de un objeto y obtener 2°, 35 una segunda imagen 20 de una segunda parte del objeto, estando la segunda parte superpuesta sustancialmente con la primera parte. Las imágenes pueden obtenerse de la memoria, o directamente desde un dispositivo de imagen. El método comprende además obtener un mapeo 110 entre la primera imagen 10 y la segunda imagen 20. El mapeo puede ser disponible desde la memoria, o puede producirse sobre la marcha aplicando un algoritmo de registro.
40 La segunda imagen 20 está segmentada 200 para obtener una segmentación, es decir, diferentes elementos (píxeles o vóxeles) de la segunda imagen 20 se clasifican en diversas categorías predefinidas. La información de la segmentación se transfiere a la primera imagen 10 sobre la base del mapeo mencionado anteriormente. En la primera imagen 10, se detectaron valores atípicos 300 sobre la base de la segmentación; esta etapa consiste por lo tanto en identificar los valores de intensidad extrema de los elementos dentro de una o más clases de elementos. En 45 la segunda imagen 20, los elementos que corresponden, de acuerdo con el mapeo, a al menos algunos de los valores atípicos detectados de la primera imagen 10, se reemplazaron 500 con valores de reemplazo, a fin de obtener una
segunda imagen corregida. Los valores de reemplazo son preferiblemente valores promedio de elementos clasificados dentro de la misma clase en elementos no atípicos en la proximidad de los respectivos valores atípicos. 50 La segmentación se actualizó realizando la segmentación 200 en la segunda imagen corregida. La detección 300, el reemplazo 500 y la actualización 200 se realizan sucesivamente hasta que se cumple un criterio de convergencia predeterminado 550. El criterio de convergencia se define convenientemente para detener la sucesión en un punto en el que ya no existe un cambio significativo en la segmentación en el tejido y en la lesión con relación a la sucesión anterior. Una vez que se cumple el criterio, la información final de la segmentación se envía a la salida 55 deseada.
La presente solicitud describe más en particular un método que mide los volúmenes de materia blanca, materia gris y líquido cefalorraquídeo en presencia de lesiones de la materia blanca basado en dos exploraciones del cerebro MR (una imagen MR ponderada en T1 y una imagen MR FLAIR). Además, se miden los volúmenes de las lesiones presentes y del tejido cerebral parenquimal.
5 Las lesiones pueden deberse a enfermedades cerebrales tales como, por ejemplo, Esclerosis Múltiple (MS) Los volúmenes de materia blanca y gris son bio marcadores importantes para las enfermedades neurodegenerativas, el volumen de la lesión para las enfermedades neurológicas que involucran lesiones focales.
Una realización más en particular del método de acuerdo con la presente invención se describirá ahora con referencia a la figura 2 (descripción general) y a la figura 3 (diagrama detallado de flujo). La interdependencia de las 10 diversas etapas del método está indicada por medio de las flechas en los diagramas. El orden en el que las diferentes etapas del método descrito a continuación no pretende limitar la invención, a menos que una dependencia temporal en particular sea implícitamente o explícitamente dada por la relación causal entre las etapas. Ciertas etapas se describirán a continuación desde un punto de vista funcional. Los expertos podrán llevar a cabo estas etapas consultando la literatura publicada relevante, sin tener que recurrir a una experimentación indebida. La 15 literatura relevante incluye en particular las siguientes publicaciones:S.
- OURSELIN et al., Robust Registration of Multi-modal Images: Towards Real-Time Clinical Applications, Medical Image Computing and Computer-Assisted Intervention (MICCAI'02);
- M. MODAT et al., Fast free-form deformation using graphics processing units, Computer Methods and Programs in Biomedicine, Volumen 98, Número 3, páginas 278-284, junio de 2010;
20 - M. MODAT, Efficient dense non-rigid registration using the free-form deformation framework, Doctoral Thesis
University College Londres, 2012;
- M. CARDOSO et al., Adaptive neonate segmentation, 2011 Medical Image Computing and Computer- Assisted Intervention (MICCAI'11).
Los expertos apreciarán que la siguiente descripción no es limitativa y que no es esencial adherirse estrictamente a 25 las técnicas descritas en la literatura citada. Los elementos de la invención reivindicada se pueden realizar de una manera diferente sin apartarse del alcance de la invención tal como se define en las reivindicaciones.
La primera etapa 100 se denominará "pre procesamiento" y consiste en tres etapas:
- En la primera etapa 110, la imagen FLAIR 10 de entrada del paciente está rigurosamente co-registrada con la imagen de entrada ponderada-T1 20 (Ourselin et al., 2002).
30 - En la segunda etapa 120, la imagen de entrada ponderada-T1 20 está desprovista del cráneo clasificando cada
voxel como una región del cerebro o una región no cerebral basada en el registro afín de una máscara cerebral disponible de un atlas (MNI) usando NiftyReg (Modat et al., 2010). El registro se realiza utilizando una transformación afín multi resolución, basada en el esquema de Mínimos Cuadrados Recortados y un enfoque de emparejamiento de bloques (Ourselin et al., 2002), seguida de un registro multi resolución no rígido basado en la 35 deformación de Forma Libre (Modat et al., 2012).
- En la tercera etapa 130, los antecedentes anatómicos probabilísticos para la materia gris (GM), la materia blanca (WM) y el líquido cefalorraquídeo (CSF), que también están disponibles en el atlas cerebral M, se transfieren al espacio de imagen ponderada-T1 usando un registro afín (Ourselin et al., 2002) seguido de un registro no rígido (Modat et al., 2010).
40 En la segunda etapa 200, los tres antecedentes probabilísticos del tejido mencionado anteriormente, junto con la imagen T1 20 desprovista del cráneo, actúan como un conocimiento previo a un algoritmo de maximización de expectativas (EM) (Cardoso et al., 2011). El algoritmo modela las intensidades de cada clase de tejido como una distribución normal, supone un campo de polarización gaussiano distribuido para la corrección de la falta de uniformidad de la intensidad y contiene un modelo de consistencia espacial basado en el Campo Aleatorio Markov 45 (MRF). El algoritmo estima sucesivamente los parámetros de cada clase de tejido, así como los parámetros del campo de polarización, y mantiene la consistencia espacial hasta la convergencia. Después de la convergencia del algoritmo EM, la imagen T1 20 se corrige por polarización y se segmenta en las tres clases de tejido, es decir, GM, WM y CSF.
En la tercera etapa 301 (que corresponde a una primera parte de la etapa 300 de la figura 1), las tres 50 segmentaciones de clase del tejido de la etapa anterior 200, junto con la imagen FLAIR co-registrada del mismo paciente, se utilizan para estimar los valores atípicos en la imagen FLAIR 10. Esto se realiza utilizando el mismo algoritmo de maximización de expectativas (EM) (Cardoso et al., 2011), como se describe en la segunda etapa 200, pero ahora se incluye un mapa atípico y las segmentaciones de la imagen T1 20 se utilizan como información previa. En otras palabras, al modelar las intensidades de cada clase de tejido como una distribución normal, la desviación 55 de cada intensidad de estas distribuciones normales se estima y se modela también como un mapa de convicciones
5
10
15
20
25
30
35
40
45
50
55
atípicas. El algoritmo estima sucesivamente el mapa de convicciones atípicas y, después de la convergencia, se produce una imagen de convicción atípica. Esta imagen de convicción atípica es una inicialización para la segmentación de lesiones MS.
En la cuarta etapa 302 (que corresponde a una segunda parte de la etapa 300 de la figura 1), las lesiones se estiman a partir del mapa de valores atípicos ya que no todos los valores son una lesión (por ejemplo, el mapa atípico podría incluir efectos de volumen parcial, objetos, etc.). A fin de diferenciar las lesiones de tales valores atípicos que no son lesiones, se debe incorporar alguna información adicional a priori sobre la posición y la apariencia de las lesiones. En particular, el mapa atípico se enmascara. La máscara es un subconjunto de la máscara del cerebro donde solo se incluyen intensidades que sean híper intensas hacia GM. Además, los valores atípicos deben estar dentro de WM (se usa una máscara para excluir valores atípicos de GM y CSF, así como también regiones que comúnmente muestren objetos). Finalmente, cada valor atípico debe tener un número mínimo de vóxeles.
En una quinta etapa 500, esta segmentación de la lesión se usa luego para rellenar las lesiones en la imagen T1 corregida por polarización 20 con sus intensidades WM vecinas.
Posteriormente, se repiten las etapas 2 (200), 3 (301), 4 (302) y 5 (500) hasta que no haya un cambio significativo en el tejido y en la segmentación de la lesión. La idea de repetir la segunda y la tercera etapas es que las lesiones son principalmente WM, por lo tanto, el relleno de la lesión T1 dará como resultado una mejor segmentación de los tejidos del cerebro, lo que a su vez dará como resultado una mejor segmentación de las lesiones.
Opcionalmente, después de la última sucesión, las lesiones se recuperan del GM en caso de que la convicción atípica sea alta. Estas lesiones se agregan a las lesiones previamente encontradas y luego estas lesiones segmentadas finales de FLAIR 10 se dilatan para enfocar mejor el volumen de las lesiones. Posteriormente, la imagen T1 ponderada 20 se rellena una vez más y se segmenta, proporcionando nuevamente las segmentaciones finales WM, GM y CSF.
Opcionalmente, los volúmenes y el recuento de las lesiones se determinan dentro de diferentes regiones anatómicas 600, por ejemplo transformando las etiquetas de la región de un atlas anatómico. Esto conduce a los valores estimados finales de volumen WM, GM y CSF 710, y un recuento y volumen de las lesiones 720.
La presente invención también se refiere a un sistema que comprende un procesador de imágenes configurado para llevar a cabo los métodos descritos anteriormente. El procesador de imágenes puede realizarse en hardware exclusivo (por ejemplo, ASIC), hardware configurable (por ejemplo, FPGA), componentes programables (por ejemplo, un procesador DSP o un procesador de propósito general con el software apropiado), o cualquier combinación de los mismos. El(Los) mismo(s) componente(s) también puede(n) incluir otras funciones.
Un diseño de ejemplo de dicho sistema se ilustra esquemáticamente en la figura 4. El sistema 1000 comprende un interfaz de entrada 1030 adaptado para recibir las imágenes de entrada a utilizar a partir de un aparato de adquisición de imagen o de un elemento de almacenamiento. Sin pérdida de lo general, se ilustran dos elementos de almacenamiento 1010, 1020, que representan dispositivos de almacenamiento respectivos para la imagen FLAIR 10 y para la imagen T1 20. El sistema 1000 comprende además un procesador 1040, conectado operativamente al interfaz de entrada 1030 y a un interfaz de salida 1060. El procesador 1040 es capaz de ejecutar un programa almacenado en la memoria 1050, que está configurado para hacer que el procesador 1040 lleve a cabo el método de acuerdo con la presente invención.
El término "interfaz" se utiliza para designar la combinación de hardware y software o firmware necesarios para permitir un intercambio de datos entre el procesador 1040 y los componentes que proporcionan o reciben los datos correspondientes. El interfaz de entrada 1030 y el interfaz de salida 1060 pueden compartir hardware común. En particular, el interfaz puede ser un interfaz de red de área local (LAN), tal como un interfaz de acuerdo con la norma "Ethernet" IEEE 802.3, en la que se basan los protocolos de red y transporte apropiados,tal como una pila TCP/IP. Los interfaces pueden proporcionar acceso a una red de área de almacenamiento (SAN) o almacenamiento conectado a la red (NAS), destinadas a almacenar las imágenes de entrada que se utilizarán y/o el resultado del análisis. Los interfaces pueden proporcionar acceso a una red de área ancha (WAN), tal como Internet, que incluye otros ordenadores que proporcionan las imágenes de entrada a utilizar y/o recuperan el resultado del análisis.
La presente invención también se refiere a un producto de programa de ordenador que comprende un código configurado para hacer que un procesador configurado lleve a cabo los métodos descritos anteriormente. El producto de programa de ordenador puede comprender medios interpretables por ordenador, tales como cinta magnética, disco magnético, disco óptico, memoria de semiconductor, o similares, que tienen los medios de código almacenados en los mismos.
Aunque la invención se ha descrito anteriormente con referencia a realizaciones específicas, esto se ha hecho para aclarar y no limitar la invención. Los expertos apreciarán que son posibles diversas modificaciones y diferentes combinaciones de las características descritas sin apartarse del alcance de la invención tal como se define en las reivindicaciones

Claims (12)

  1. 5
    10
    15
    20
    25
    30
    35
    40
    45
    REIVINDICACIONES
    1. Un método de análisis de datos de imágenes, que comprende:
    - obtener (1) una primera imagen (10) de una primera parte de un objeto;
    - obtener (2) una segunda imagen (20) de una segunda parte del objeto, teniendo dicha segunda parte una superposición sustancial con dicha primera parte;
    - obtener un mapeo (110) entre dicha primera imagen (10) y dicha segunda imagen (20);
    - segmentar (200) dicha segunda imagen (20) para obtener una segmentación de tejidos;
    - detectar valores atípicos (300) en dicha primera imagen (10) identificando los valores de intensidad extrema de elementos dentro de una o más clases de elementos basándose en dicha segmentación;
    - reemplazar (500) los elementos de dicha segunda imagen (20) que corresponden, de acuerdo con dicho mapeo, a al menos algunos de dichos valores atípicos detectados de dicha primera imagen (10), con valores de reemplazo, para obtener una segunda imagen corregida, siendo los al menos algunos de dichos valores atípicos detectados de dicha primera imagen, una segmentación de la lesión; y
    - actualizar dicha segmentación realizando dicha segmentación (200) en dicha segunda imagen corregida;
    en el que dicha detección de los valores atípicos (300), dicho reemplazo (500) y dicha actualización (200) se realizan sucesivamente hasta que se cumpla un criterio de convergencia predeterminado (550), representando dicho criterio de convergencia predeterminado un punto en el cual no existen cambios significativos en el tejido y en la segmentación de la lesión.
  2. 2. El método de acuerdo con la reivindicación 1, en el que dichos valores atípicos se detectan (300) segmentando dicha primera imagen (10) empezando por dicha segmentación (200) de dicha segunda imagen (20) con una clase adicional de valores atípicos.
  3. 3. El método de acuerdo con cualquiera de las reivindicaciones precedentes, que comprende además iniciar dicha segmentación (200) transfiriendo (130) una segmentación de un ejemplo a dicha segunda imagen (20).
  4. 4. El método de acuerdo con cualquiera de las reivindicaciones precedentes, en el que dicha obtención de dicho mapeo (110) comprende realizar un registro de dicha primera imagen obtenida (10) a dicha segunda imagen obtenida (20).
  5. 5. El método de acuerdo con cualquiera de las reivindicaciones precedentes, en el que dicho objeto comprende al menos una parte de un cerebro, dicha primera imagen (10) es una imagen FLAIR, dicha segunda imagen (20) es una imagen T1;
    en el que dicha segmentación comprende una clasificación de elementos de dicha imagen T1 (20) como materia gris, materia blanca o líquido cefalorraquídeo; y
    en el que dichos valores atípicos se detectan (300) entre los elementos de dicha primera imagen (10) clasificados como materia blanca.
  6. 6. El método de acuerdo con la reivindicación 5, en el que dichos valores de reemplazo se basan en valores promedio de elementos clasificados como materia blanca en elementos no atípicos en la proximidad de dichos respectivos valores atípicos.
  7. 7. El método de acuerdo con cualquiera de las reivindicaciones 5 a 6, que comprende además eliminar (120) elementos que representan tejido no cerebral de dicha segunda imagen (20).
  8. 8. El método de acuerdo con cualquiera de las reivindicaciones 5 a 7, en el que dicha
    segmentación (200) se usa para calcular (600) volúmenes respectivos o áreas de materia gris, materia blanca y
    líquido cefalorraquídeo (710).
  9. 9. El método de acuerdo con cualquiera de las reivindicaciones 5 a 8, en el que dicha
    detección de valores atípicos (300) se usa para calcular (600) volúmenes o áreas de dichos valores atípicos (720).
  10. 10. El método de acuerdo con cualquiera de las reivindicaciones 8 a 9, que comprende adicionalmente calcular (600) los volúmenes o áreas en diferentes regiones anatómicas.
  11. 11. Un producto de programa de ordenador que comprende un medio de código configurado para hacer que un procesador lleve a cabo el método de acuerdo con cualquiera de las reivindicaciones precedentes.
  12. 12. Un sistema de procesamiento de imágenes (1000) que comprende:
    - un interfaz de entrada (1030) adaptado para recibir una primera imagen (10) de una primera parte de un objeto y una segunda imagen (20) de una segunda parte de dicho objeto, teniendo dicha segunda parte un solapamiento sustancial con dicha primera parte;
    - un procesador (1040), conectado operativamente a dicho interfaz de entrada (1030); y
    5 - un interfaz de salida (1060), conectado operativamente a dicho procesador (1040) y adaptado para presentar los
    resultados de salida producidos por dicho procesador (1040); y
    - una memoria (1050), operativamente conectada a dicho procesador (1040) y dispuesta para comprender el código a ejecutar por dicho procesador (1040), comprendiendo dicha memoria (1050) un código configurado para hacer que dicho procesador (1040) lleve a cabo el método de cualquiera de las reivindicaciones 1 a 10 usando dicha primera
    10 imagen (10) y dicha segunda imagen (20) como entradas.
ES14184117.1T 2014-09-09 2014-09-09 Método y sistema para analizar datos de imágenes Active ES2685977T3 (es)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
EP14184117.1A EP2996085B1 (en) 2014-09-09 2014-09-09 Method and system for analyzing image data

Publications (1)

Publication Number Publication Date
ES2685977T3 true ES2685977T3 (es) 2018-10-15

Family

ID=51542156

Family Applications (1)

Application Number Title Priority Date Filing Date
ES14184117.1T Active ES2685977T3 (es) 2014-09-09 2014-09-09 Método y sistema para analizar datos de imágenes

Country Status (4)

Country Link
US (1) US10198815B2 (es)
EP (1) EP2996085B1 (es)
ES (1) ES2685977T3 (es)
WO (1) WO2016038118A1 (es)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10635894B1 (en) * 2016-10-13 2020-04-28 T Stamp Inc. Systems and methods for passive-subject liveness verification in digital media
US11373449B1 (en) * 2016-10-13 2022-06-28 T Stamp Inc. Systems and methods for passive-subject liveness verification in digital media
CN107025457B (zh) * 2017-03-29 2022-03-08 腾讯科技(深圳)有限公司 一种图像处理方法和装置
US10474667B2 (en) * 2017-07-29 2019-11-12 Vmware, Inc Methods and systems to detect and correct outliers in a dataset stored in a data-storage device
US11885861B2 (en) * 2018-02-26 2024-01-30 The Board Of Regents Of The University Of Texas System 3-dimensional representations of post-contrast enhanced brain lesions
US11093771B1 (en) 2018-05-04 2021-08-17 T Stamp Inc. Systems and methods for liveness-verified, biometric-based encryption
US11496315B1 (en) 2018-05-08 2022-11-08 T Stamp Inc. Systems and methods for enhanced hash transforms
US11301586B1 (en) 2019-04-05 2022-04-12 T Stamp Inc. Systems and processes for lossy biometric representations
CN110047046B (zh) * 2019-04-15 2022-09-06 南京邮电大学 图像分割和偏移场矫正方法、可读存储介质和终端
EP4034241A4 (en) * 2019-09-24 2023-10-11 The Board Of Regents Of The University Of Texas System METHOD AND SYSTEMS FOR ANALYZING BRAIN LESIONS USING LONGITUDINAL 3D MRI DATA
US11967173B1 (en) 2020-05-19 2024-04-23 T Stamp Inc. Face cover-compatible biometrics and processes for generating and using same
CN112308853A (zh) * 2020-10-20 2021-02-02 平安科技(深圳)有限公司 电子设备、医学图像指标生成方法、装置及存储介质
US11610306B2 (en) 2020-12-16 2023-03-21 Industrial Technology Research Institute Medical image analysis method and device
US12079371B1 (en) 2021-04-13 2024-09-03 T Stamp Inc. Personal identifiable information encoder
CN113256654B (zh) * 2021-05-26 2023-09-12 西安交通大学 一种基于病变扩散算法的白质高信号提取方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5983251A (en) * 1993-09-08 1999-11-09 Idt, Inc. Method and apparatus for data analysis
US20160034638A1 (en) * 2013-03-14 2016-02-04 University Of Rochester System and Method for Detecting Population Variation from Nucleic Acid Sequencing Data

Also Published As

Publication number Publication date
EP2996085B1 (en) 2018-05-23
WO2016038118A1 (en) 2016-03-17
EP2996085A1 (en) 2016-03-16
US10198815B2 (en) 2019-02-05
US20170256056A1 (en) 2017-09-07

Similar Documents

Publication Publication Date Title
ES2685977T3 (es) Método y sistema para analizar datos de imágenes
US11699236B2 (en) Systems and methods for the segmentation of multi-modal image data
US10842445B2 (en) System and method for unsupervised deep learning for deformable image registration
Roy et al. Robust skull stripping using multiple MR image contrasts insensitive to pathology
WO2018001099A1 (zh) 一种血管提取方法与系统
Tourbier et al. Automated template-based brain localization and extraction for fetal brain MRI reconstruction
Jurek et al. CNN-based superresolution reconstruction of 3D MR images using thick-slice scans
Lee et al. Semi-automatic segmentation for 3D motion analysis of the tongue with dynamic MRI
Barman et al. MRI image segmentation using level set method and implement an medical diagnosis system
Kalavathi et al. Automatic segmentation of cerebral hemispheres in MR human head scans
Pszczolkowski et al. Automated segmentation of haematoma and perihaematomal oedema in MRI of acute spontaneous intracerebral haemorrhage
US20170046826A1 (en) High resolution intraoperative mri images
Yang et al. Geodesic active contours with adaptive configuration for cerebral vessel and aneurysm segmentation
Pei et al. Improved brain tumor growth prediction and segmentation in longitudinal brain MRI
KR101685821B1 (ko) 흉부 X-ray 영상에서 몸체 및 관심영역 분할 방법 및 시스템
Grandhe et al. Adaptive analysis & reconstruction of 3D DICOM images using enhancement based SBIR algorithm over MRI
CN113412502A (zh) 患者特定皮质表面密铺为偶极子贴片
Al-Shaikhli et al. Medical image segmentation using multi-level set partitioning with topological graph prior
Myint et al. Effective kidney segmentation using gradient based approach in abdominal CT images
Rivest-Hénault et al. Inverse-consistent rigid registration of CT and MR for MR-based planning and adaptive prostate radiation therapy
Kuijf et al. Assessment of methods to extract the mid-sagittal plane from brain MR images
Liu et al. Automatic detection of the anterior and posterior commissures on MRI scans using regression forests
Dore et al. Surface-Base approach using a multi-scale EM-ICP registration for statistical population analysis
Bergo et al. Fast, accurate and precise mid-sagittal plane location in 3D MR images of the brain
Korsager et al. The use of an active appearance model for automated prostate segmentation in magnetic resonance