ES2685977T3 - Método y sistema para analizar datos de imágenes - Google Patents
Método y sistema para analizar datos de imágenes Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000011218 segmentation Effects 0.000 claims abstract description 53
- 230000003902 lesion Effects 0.000 claims abstract description 34
- 238000013507 mapping Methods 0.000 claims abstract description 13
- 238000013450 outlier detection Methods 0.000 claims abstract description 3
- 210000004885 white matter Anatomy 0.000 claims description 24
- 210000004884 grey matter Anatomy 0.000 claims description 17
- 210000001175 cerebrospinal fluid Anatomy 0.000 claims description 15
- 210000001519 tissue Anatomy 0.000 claims description 15
- 210000004556 brain Anatomy 0.000 claims description 14
- JXSJBGJIGXNWCI-UHFFFAOYSA-N diethyl 2-[(dimethoxyphosphorothioyl)thio]succinate Chemical compound CCOC(=O)CC(SP(=S)(OC)OC)C(=O)OCC JXSJBGJIGXNWCI-UHFFFAOYSA-N 0.000 claims description 10
- 210000003484 anatomy Anatomy 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 7
- 210000005013 brain tissue Anatomy 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000000977 initiatory effect Effects 0.000 claims description 2
- 238000007405 data analysis Methods 0.000 abstract 1
- 201000006417 multiple sclerosis Diseases 0.000 description 12
- 238000002595 magnetic resonance imaging Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 9
- 206010072731 White matter lesion Diseases 0.000 description 5
- 208000027418 Wounds and injury Diseases 0.000 description 5
- 230000006378 damage Effects 0.000 description 5
- 208000014674 injury Diseases 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 230000010287 polarization Effects 0.000 description 3
- 210000003625 skull Anatomy 0.000 description 3
- 208000012902 Nervous system disease Diseases 0.000 description 2
- 208000025966 Neurological disease Diseases 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 239000000090 biomarker Substances 0.000 description 2
- 230000002490 cerebral effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 208000014644 Brain disease Diseases 0.000 description 1
- 206010051290 Central nervous system lesion Diseases 0.000 description 1
- 208000034189 Sclerosis Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 210000003169 central nervous system Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011496 digital image analysis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 210000001652 frontal lobe Anatomy 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000008570 general process Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004770 neurodegeneration Effects 0.000 description 1
- 208000015122 neurodegenerative disease Diseases 0.000 description 1
- 230000000926 neurological effect Effects 0.000 description 1
- 210000000869 occipital lobe Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 210000001152 parietal lobe Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 210000003478 temporal lobe Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2431—Multiple classes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/143—Segmentation; Edge detection involving probabilistic approaches, e.g. Markov random field [MRF] modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/174—Segmentation; Edge detection involving the use of two or more images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10084—Hybrid tomography; Concurrent acquisition with multiple different tomographic modalities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20128—Atlas-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; 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)
- 51015202530354045REIVINDICACIONES1. 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. 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. 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. 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. 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; yen el que dichos valores atípicos se detectan (300) entre los elementos de dicha primera imagen (10) clasificados como materia blanca.
- 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. 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. El método de acuerdo con cualquiera de las reivindicaciones 5 a 7, en el que dichasegmentación (200) se usa para calcular (600) volúmenes respectivos o áreas de materia gris, materia blanca ylíquido cefalorraquídeo (710).
- 9. El método de acuerdo con cualquiera de las reivindicaciones 5 a 8, en el que dichadetección de valores atípicos (300) se usa para calcular (600) volúmenes o áreas de dichos valores atípicos (720).
- 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. 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. 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); y5 - un interfaz de salida (1060), conectado operativamente a dicho procesador (1040) y adaptado para presentar losresultados 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 primera10 imagen (10) y dicha segunda imagen (20) como entradas.
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)
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)
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 |
-
2014
- 2014-09-09 EP EP14184117.1A patent/EP2996085B1/en active Active
- 2014-09-09 ES ES14184117.1T patent/ES2685977T3/es active Active
-
2015
- 2015-09-09 WO PCT/EP2015/070663 patent/WO2016038118A1/en active Application Filing
- 2015-09-09 US US15/510,151 patent/US10198815B2/en active Active
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 |