[en] This is a multidisciplinary thesis in Geographical Information Science. Methodologically, it draws on Geostatistics, Remote Sensing, Computer Science and Geographical Information Systems. The thesis studies certain algorithms for spatial analysis and image processing, and proposes new ways of automatize them. This automatization is meant to enable and improve the efficiency of massive processing of big geospatial data. The automatization of these processes is based on: 1) The use of metadata in automated decision making. This metadata can be generated either by the supplier of the original data or by the intermediate results derived from the processes involved in the algorithms. Allowing them to drive the processes of spatial analysis, gives metadata a geographic dimension that is missing in more common uses, such as catalogue search. 2) The study and specific adaptation of the algorithm to automatization. Each algorithm has been studied and redesigned to achieve an automatic solution while preserving the highest possible quality: kriging, variogram analysis, optimal interpolation, multivariant or logistic regression models, and radiometric corrections. In addition, automatization allows us to distribute computing amongst different processing units, significantly reducing execution time (speed-ups close to 60 have been achieved in some instances, while in others behaviour was nearly linear). These computational improvements result from parallel programming in C language and from sending data between processors through Message Passing Interface (MPI) in a High Performance Computing (HPC) environment. Automatization has been applied to various fields and various types of data, first and foremost to temporal series of Remote Sensing images, our main domain of research and application. For example, this thesis proposes a method that improves the selection of high quality images (Terra MODIS reflectance) with regard to the selection based on the quality masks accompanying these products. The new method adds an automatic variogram analysis capable of detecting possible radiometric and geometric anomalies with respect to a spatial pattern type, defined with geostatistical criteria. Another aim of this thesis is to automatize the radiometric correction of Landsat images from various sensors, on the basis of radiometric reference values from pseudo-invariant polygons that are also automatically generated. This contribution makes it possible to process long temporal series of high spatial resolution images covering large geographic areas, and to do so with high radiometric accuracy. This could be applied to geographical analysis at a level of detail not feasible until now. Finally, in the field of lossy compression for Remote Sensing images, the thesis proposes a view that is more spatial (and therefore geographical) than more traditional views in the domain of image processing. Image compression is widely accepted today as a methodology to efficiently manage big data. Here again, geostatistical analysis has allowed us to compress images with a minimum loss of their spatial properties.

[es] Se presenta una tesis multidisciplinar en el contexto de la Ciencia de la Información Geográfica dónde convergen sinérgicamente metodologías de la Geoestadística, la Teledetección, la Ciencia Computacional y los Sistemas de Información Geográfica con el objetivo de estudiar y ofrecer nuevas soluciones que permitan la automatización de determinados algoritmos de análisis espacial y de procesamiento de imágenes. La automatización de estos algoritmos está enfocada a posibilitar y mejorar la eficiencia del procesamiento masivo de grandes volúmenes de datos geoespaciales. La automatización de estos procesos está fundamentada en dos pilares principales: • El aprovechamiento adecuado de los metadatos en la toma de decisiones automáticas. Estos metadatos pueden ser generados tanto por el proveedor de los datos originales como por los resultados intermedios derivados de los procesos implicados en los algoritmos estudiados. Este importante rol de conductor de los procesos de análisis espacial añade a los metadatos una componente más geográfica que otros usos más habituales como son las buscas sobre catálogos. • El estudio y adaptación específico del algoritmo a la automatización. Se ha analizado y rediseñado cada algoritmo hasta conseguir una solución automática que mantenga la máxima calidad posible: kriging, análisis del variograma, interpoladores óptimos, modelos de regresión multivariante o logística adecuados y correcciones radiométricas. Adicionalmente, la automatización ha permitido aplicar las metodologías de distribución de cómputo entre diferentes unidades de procesamiento y reducir de forma significativa los tiempos de ejecución (se ha conseguido speed-ups próximos a 60 en algún caso, y próximos a un comportamiento lineal en otras). Estas mejoras computacionales se han llevado a cabo mediante la paralelización de los códigos en lenguaje C y el envío de datos entre procesadores mediante Message Passing Inteface (MPI) en un entorno de elevadas prestaciones computacionales (HPC). Las automatizaciones se han aplicado a varios campos y a varios tipos de datos, siendo las series temporales de imágenes de Teledetección el principal ámbito de investigación y aplicación. Por ejemplo, en esta tesis se propone un método que mejora la selección de imágenes (reflectividades de Terra MODIS) de elevada calidad respecto a la que sólo se efectúa en función de las máscaras de calidad que acompañan estos productos. El nuevo método añade un análisis automático del variograma que detecta las posibles anomalías radiométricas y geométricas respecto a un patrón espacial tipo, definido también con criterios geoestadísticos. Conseguir corregir radiométricamente de forma automática imágenes Landsat de varios sensores en base a valores radiométricos de referencia desde polígonos pseudo-invariantes, también generados automáticamente, es otro de los logros de esta tesis. Esta aportación permite procesar largas series temporales de imágenes de elevada resolución espacial para ámbitos geográficos extensos con un elevado rigor radiométrico y, por lo tanto, abre la puerta a realizar análisis del territorio con un niveles de detalle hasta ahora no viables. Finalmente, en el campo de la compresión con pérdida de imágenes de Teledetección se ha realizado aportaciones que añaden una visión más espacial (y por lo tanto territorial) sobre visiones más clásicas en el ámbito general del tratamiento de imágenes. La compresión de imágenes es una de las metodologías ampliamente aceptadas hoy para gestionar de forma eficiente grandes volúmenes de datos y, nuevamente, el análisis geoestadístico ha permitido conducir las compresiones hacia resultados que conserven mejor las propiedades espaciales de las imágenes comprimidas.