El análisis, llevado a cabo en R, emplea diferentes paquetes: solaR para el cálculo de geometría solar, irradiación global en el plano horizontal y efectiva en el plano inclinado; raster y rasterVis para la lectura, manipulación y visualización de imágenes raster; gstat y sp para los métodos de interpolación geoestadística. El primer paso es realizar un mapa del promedio anual de la irradiación global diaria en el plano horizontal de los datos de CMSAF. Será este mapa el que se compare con las estaciones meteorológicas (sin interpolar) y con el mapa obtenido de la interpolación de las estaciones terrestres, distribuidas de manera no uniforme. Se aplican técnicas no estadísticas como el Inverse Distance Weighted (IDW) y el ajuste por superficie que ofrecen aproximaciones groseras. Se ensayan también técnicas de análisis estadístico como el Ordinary Kriging, solución poco fiable en aquellos emplazamientos con poca densidad de estaciones. Finalmente, se opta por el método de kriging with external drift, empleando la irradiación anual según CMSAF como variable externa.
Este procedimiento se aplica tanto para la estimación de irradiación global en el plano horizontal, como para la irradiación efectiva en un sistema fotovoltaico fijo, en un sistema con un eje de seguimiento y en un sistema con doble eje de seguimiento.