Desarrollo de técnicas específicas de procesado de imagen para su aplicación a la medida del grosor íntima-media de la arteria carótida sobre imágenes de ultrasonidos

Abstract

[ESP] Las enfermedades cardiovasculares son una de las principales causas de mortalidad del mundo. Tras la mayoría de muertes por enfermedad cardiovascular la principal responsable es la arteriosclerosis. La arteriosclerosis consiste en un engrosamiento progresivo del tejido vascular que reduce la elasticidad de los vasos sanguíneos afectados y puede llegar incluso a obstruirlos. Esta enfermedad se desarrolla en la infancia y adolescencia, pudiendo llegar a pasar desapercibida toda la vida o bien actuando como detonador de otras afecciones más serias, como infartos, derrames cerebrales o isquemias. Por todo ello, la detección precoz de la arteriosclerosis resulta de vital importancia. En la actualidad, se viene usando el grosor íntima-media o IMT de la arteria carótida común como indicador fiable y precoz de la arteriosclerosis. Este indicador mide el grosor entre las capas íntima y media de la carótida común en cortes longitudinales de la arteria en imágenes ecográficas. El uso de los ultrasonidos para extraer esta medida, además de barato, resulta no invasivo para el paciente. Sin embargo, esta modalidad de imagen no está exenta de desventajas, como el elevado nivel de ruido que presenta o la alta dependencia del operador. La medida del IMT se realiza de manera manual sobre imágenes de ultrasonidos. Para ello, un observador experto realiza de una a cinco mediciones manuales del grosor de la arteria, donde cada medición consiste en un par de puntos. La presente Tesis Doctoral pretende realizar una segmentación automática de las paredes de la arteria carótida común, de modo que, en lugar de un conjunto de puntos limitado, se pueda extraer el IMT en toda la longitud de la arteria presente en la ecografía. De forma adicional al IMT, se segmenta también la pared anterior de la arteria, pudiendo así proporcionar no sólo el grosor del vaso, sino también el diámetro del cauce de la arteria. El hecho de que sea una segmentación totalmente automática, evita la interacción con el usuario existente en otros métodos de ayuda a la medida del IMT, a la vez que elimina la subjetividad de la medida. Para llevar a cabo una delineación automática de las paredes arteriales, el método desarrollado se implementa en dos etapas básicas. La primera consiste en la detección automática del lumen o cauce por el que fluye la sangre. En una segunda etapa, tras la detección del lumen, se refina el resultado mediante el uso de contornos activos o snakes. En esta Tesis, se emplean contornos activos implementados en el dominio frecuencial. Esta implementación consigue una importante reducción en el coste computacional respecto a la formulación original. Como función de forma se emplean B-splines cúbicas, que presentan una excelente relación entre rendimiento y tiempo de ejecución. El uso de splines proporciona unos contornos finales suaves, lo que evita la típica rugosidad presente en las imágenes de ultrasonidos. Los resultados tras el algoritmo de contornos activos se validan automáticamente para evitar la inclusión de tramos erróneos en las medidas finales. Estos tramos erróneos están ocasionados, principalmente, por la presencia de huecos en la imagen que dificultan el proceso de segmentación, por no proporcionar ninguna información a la imagen de fuerzas externas. Esta validación de los resultados se basa en dos estrategias: estadística y de intensidad. Combinando ambas estrategias se descartan tramos en los que no haya información o en los que la medida resulte anatómicamente improbable. Se ha llevado a cabo una caracterización de los resultados extensa, empleando cuatro métricas distintas para la evaluación de las distancias entre contornos. Tomando el promedio de cuatro segmentaciones manuales como ground truth, se ha comparado con la segmentación automática calculada por el algoritmo de snakes. Tanto las medidas de IMT y del diámetro del lumen, como las distancias entre trazados manuales y automáticos se han evaluado para las cuatro métricas consideradas. En cada caso, se ha calculado el coeficiente de correlación de Pearson, la distribución del error, así como los diagramas de Bland-Altman. A diferencia de otros métodos, en los que la resolución espacial es la misma para todas las imágenes, la configuración del ecógrafo se ha dejado a criterio de los radiólogos. De esta manera, la resolución espacial varía de una imagen a otra, situándose en un rango de 0,029 mm/píxel a 0,082 mm/píxel. Para las 58 imágenes analizadas, se ha logrado medir el IMT en todas ellas. El error medio midiendo el IMT es comparable e incluso inferior al de otros métodos automáticos. El error de segmentación o la distancia directa entre trazados manuales y automáticos es muy bajo, siendo como máximo de unos 2 píxeles para cualquiera de los tres contornos delineados. Como ayuda a la transferencia de resultados, se ha implementado una interfaz de usuario instalable en cualquier PC que implementa el método desarrollado. Esto facilita su uso por parte de personal médico interesado en evaluar el riesgo cardiovascular, tanto de pacientes concretos como de un conjunto de personas para su posterior análisis. Así pues, el método aquí descrito es útil tanto en la práctica clínica diaria como para la realización de estudios del riesgo cardiovascular en la población general. [ENG] Cardiovascular diseases (CVDs) are one of the main causes of death worldwide. Amongst all deaths related to CVDs, atherosclerosis is responsible for the biggest amount of them. Atherosclerosis consists in a progressive thickening of vascular tissue, provoking a loss of elasticity and an increase of thickness of the blood vessels. This thickening can even occlude the affected vessels. This pathology, which is developed during childhood and adolescence, may be unnoticed for years before triggering other more serious conditions, such as infarction, stroke or ischemia. For these reasons, an early detection of atherosclerosis is of paramount interest. Nowadays, the intima-media thickness (IMT) of the common carotid artery is being used as a reliable and early atherosclerosis detector. IMT measures the thickness between intima and media layers of the common carotid artery in longitudinal cuts of the artery in ultrasound images. The use of ultrasound imaging is relatively cheap as well as being non invasive for the patient. However, this image modality presents some drawbacks as it is operator dependent and is quite affected by noise. IMT is manually measured on ultrasound images. With that purpose, an expert observer takes from one to five measurements of the artery thickness. Each manual measurement consists in placing a pair of markers on the artery wall. The present PhD thesis aims at automatically segmenting the layers of the common carotid artery. Thus, instead of a limited set of points, IMT can be measured along the artery length. In addition to IMT, near wall of the artery is also segmented, which provides not only the artery wall thickness, but also the diameter of the artery channel. The fact that the segmentation is completely automatic avoids the user interaction which is present in other IMT measurent methods and, thus, removes the subjetivity of the measurement. To automatically delineate the artery walls, the developed method is implemented in two main stages. The first of them comprises the automatic detection of the lumen, which is the channel where the blood flows. During the second stage, after the lumen detection, the result is refined by means of active contours or snakes. Automatic lumen detection is performed thanks to a correlation of the ultrasound image with a model of the intima-media complex. This step allows the detection of the far wall, which is placed at the bottom of the lumen. In this stage a median filtering is also used to reduce the characteristic speckle noise in ultrasound images. The median filtering gives an image with homogeneous regions, while maintaining the edges of the vessel. Therefore, the median filtering helps in the robustness of the method when there are blood turbulences or backscattering. From the edges of the lumen, active contours or snakes are initialized. The snakes refine the result of the previous stage and provide the different interfaces that define IMT and lumen diameter. Three curves are considered, one for the near wall (located above the lumen in the ultrasound image) and two for the far wall (located under the lumen in the ultrasound image), corresponding to the lumen-intima and media-adventitia interface that determine IMT. In the present work, a Fourier-domain implementation of active contours is used. This implementation achieves a considerable computational cost reduction with respect to the original formulation. Cubic B-splines have been chosen as shape function because of their excellent performance versus running time ratio. Moreover, B-splines provide soft final contours, dealing with the typical rugosity in ultrasound images. Besides the frequency implementation of snakes, the success of the refinement stage is based on the calculation of an adequate external forces image. To compute the external forces image, positive and negative transitions of the intensity gradient of the image are combined in a single image. This combination allows the simultaneous convergence of the near and far wall curves. Although the external force image contains plenty information of the edges to be detected, combined gradient image might not be clear enough. Hence, a morphological reconstruction is performed, being the mask image the combined gradient image. As marker, the cumulative result of morphological openings is used. The openings are performed over the combined gradient image, being the structuring elements lines oriented in the main directions of the gradient image. These main directions are extracted via a Hough transform. Finally, the reconstruction result gives a clearer edge image, in which small structures disappear, making the convergence of the snakes faster. The results after the active contour algorithm are automatically validated to avoid the inclusion of wrong segmented sections in the final measurements. These wrong sections are due to, mainly, the presence of gaps in the image which do not provide any information about the external forces image. This additional validation stage is based on two strategies, statistical and intensity-based. By combining both strategies, we avoid the inclusion of sections with big gaps or in which the measurements are unlikely in anatomical terms. An exhaustive result characterization has been accomplished. Four different metrics have been used to evaluate the distances between contours. Considering the average of four manual measurements as ground truth, automatic segmentation has been compared to manual segmentation. The four considered metrics have been used to evaluate the IMT, the lumen diameter and the distances between manual and automatic delineations. En each case, Pearson’s correlation coefficient, error distribution and Bland-Altman’s plots have been calculated. The ultrasound image dataset has been obtained with a single ultrasound scanner with two different probes, each of them working at different frequency range. Unlike other methods, in which the spatial resolution is not variable, the scanner configuration is flexible for the radiologist, at his/her own discretion. Thus, spatial resolution may vary from one image to another, ranging from 0.029 mm/pixel to 0.082 mm/pixel in our database. IMT has been automatically measured in all 58 images in our database. The images correspond to 35 patients, none of them previously diagnosed with atherosclerosis. The mean IMT error is comparable or even less than the error in other automatic methods. Besides, the distance between manual and automatic contours is very low, being under 2 pixels for any of the three detected interfaces. In order to disseminate and transfer the results, a user interface has been implemented. The program can be installed in any PC, which makes its use easier. Medical staff interested in the evaluation of cardiovascular risk cannot only use the interface for individual patients but also for population-based studies. Thus, the developed method is useful in daily clinical practice as well as in cardiovascular risk studies over population. The segmentation of the carotid artery layers allows the calculation of other measurements and statistics apart from the mean IMT. Maximum and minimum values, texture or ecogenicity in the intimamedia complex region can be extracted. In a similar way, it is possible to monitor the atherosclerosis evolution via registration methods. Therefore, the doctor could easily detect slight changes in the IMT and adjust the treatment and recommendations to the patient in a reliable way. Besides, since the developed technique is based on snakes, the method is easily extendable to 3D cases. This would lead to new applications of the segmentation, such as the analysis of the IMT during a whole systole and diastole cycle or the study of the elastogram of the carotid artery.