Cómo reescribir funciones aleatorias multivariadas como univariadas para hacer cokrigeage. Un algoritmo flexible

Resumen

La mayoría de los software para geoestadísticas poseen un limitado número de modelos de krigeage y no son lo suficientemente flexibles para definir nuevos o complicados modelos. Con el objetivo de obtener un algoritmo flexible de krigeage se empleó la notación de funciones aleatorias regionalizadas (RF) multivariadas escritas como una RF univariada Z(x, i), con un drift m(x, i, j), donde i y j son identificadores de variables. Se programó en C++ una librería geoestadística llamada OpenKriging, con una función que se puede llamar desde código escrito en Python, la cual es capaz de construir el sistema de krigeage de la RF mencionada, sin restricciones en el número de dimensiones espaciales, variables y monomios en el drift. Se implementaron dos tipos de drift: los geográficos y los drift externos, estos pueden ser independientes o no. La función responsable de construir el sistema de ecuaciones krigeage no se programó para un modelo predefinido, el modelo es definido por el usuario pasando desde el código en Python los datos y argumentos apropiados y modificando el sistema antes de resolverlo. La librería resultante puede asimilar muchos de los modelos conocidos de krigeage y también modelos nuevos o poco comunes; esta es fácil de modificar para incrementar su funcionalidad, la cual se mostró a través de dos ejemplos de modelos de krigeage: uno con definición complicada y el otro poco común.