Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira.

ISSN 0122-1701 y ISSN: 2344-7214 15

El ensayo no destructivo usando termografía

infrarroja en el mundo y en América Latina:
Una Revisión
Infrared thermal nondestructive testing in the world and Latin America: A Review

J. H. Erazo-Aux ; H. Loaiza-Correa ; A. D. Restrepo-Girón

DOI: [Link]
Artículo de revisión

Abstract—This paper presents a review concerning mayor y se evidenció una reducida cantidad de métricas para la
nondestructive testing by infrared thermography. The study evaluación objetiva de su rendimiento.
includes articles by well-known researchers from around the
world and Latin America. From a literature review, there is Palabras claves— ensayo no destructivo, materiales compuestos,
currently a growing trend in scientific production in this field. The revisión de literatura, termografía infrarroja.
above confirms its importance and the demand for constant
development. The works found were organized into four scope
lines: a. state of the art, b. thermal information processing I. INTRODUCCIÓN
techniques, c. performance metrics, and d. simulation, thermal
models, and applications. In addition, the paper presents the
traditional and emerging countries supporting the new
developments in the nondestructive testing area. Finally, there is a
L A Termografía Infrarroja (IRT, del inglés: infrared
thermography) es una técnica de ensayo no destructivo
que, a partir de su capacidad para capturar radiación infrarroja
trend towards processing thermal image sequences rather than
single images where are reported few developments. The number de una escena, es utilizada como método de inspección para
of methods is increasing, and there is a reduced number of metrics examinar una pieza, material o sistema sin afectar o alterar sus
for the quantitative performance evaluation. características, así como para evaluar su funcionamiento
presente y futuro. A pesar de existir una amplia variedad de
Index Terms—composite materials, infrared thermography, técnicas asociadas con Ensayos No Destructivos (NDT, del
nondestructive testing, review. inglés: nondestructive testing) que comparten algunas ventajas,
la técnica de ensayo no destructivo por termografía infrarroja
Resumen—En este artículo se sintetiza una revisión de los
principales trabajos relacionados con ensayo no destructivo
(TNDT, del inglés: infrared thermography nondestructive
usando termografía infrarroja incluyendo aportes realizados por testing) se caracteriza además por ser rápida, segura, no
investigadores de América Latina. Se encontró en este campo una invasiva y libre de contacto [1]. Características que la hacen aún
tendencia creciente de la producción académica en los últimos más atractiva en un amplio rango de aplicaciones especialmente
años, lo cual ratifica que es un área de interés y desarrollo en la industria aeroespacial, naval, automotriz y hace algunos
continuo. El resultado de la revisión se clasificó en cuatro líneas años en salud estructural entre otras [2], [3].
temáticas cubriendo principalmente: i) trabajos que recopilan el
estado del arte, ii) desarrollos en técnicas de procesamiento de
Este artículo tiene como objetivo presentar una revisión de
información térmica, iii) indicadores de desempeño y iv) literatura relacionada con ensayos no destructivos usando
simulación, modelos térmicos y campos de aplicación. termografía. Para el desarrollo de la revisión se definió el tema
Adicionalmente, se registraron los países tradicionales y central para la búsqueda, las estrategias de búsqueda, los
emergentes que están generando contribuciones en el sector del criterios de selección de la información y finalmente, se asocia
ensayo no destructivo usando termografía infrarroja. Finalmente, y presenta una síntesis de la información recopilada.
se identificó cómo el enfoque actual del procesamiento de
información térmica está orientado al procesamiento de
Este documento se encuentra organizado en cuatro secciones.
secuencias de imágenes, donde el número de técnicas es cada vez La primera hace un recorrido sobre los principales trabajos

Este manuscrito fue enviado el 29 de junio de 2021 y aceptado el 30 de H. Loaiza-Correa. Profesor Titular de la Escuela de Ingeniería Eléctrica y
marzo de 2022. Electrónica de la Universidad del Valle, Cali, VA 760032 Colombia (e-mail:
Este trabajo ha sido financiado por el Ministerio de Ciencia, Tecnología e [Link]@[Link]).
Innovación de Colombia (Minciencias); Convocatoria 727 de 2015. Además, A. D. Restrepo-Girón. Profesor Asociado de la Escuela de Ingeniería
por la Institución Universitaria Antonio José Camacho (UNIAJC). Proyectos Eléctrica y Electrónica de la Universidad del Valle, Cali, VA 760032 Colombia
de Investigación PD0916 y PD0121. (e-mail: [Link]@[Link])
J. H. Erazo-Aux. Profesor Asociado adscrito a la facultad de ingeniería de la
Institución Universitaria Antonio José Camacho, Cali, VA 760046 Colombia
(e-mail: jherazo@[Link]).
16 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

desarrollados que recopilan el estado del arte en el área. La la década de los 90 el concepto de Termografía Infrarroja
siguiente sección reúne información sobre los adelantos en Cuantitativa (QIRT, del inglés: quantitative infrared
diferentes tópicos relacionados al TNDT incluyendo los aportes thermography). Desde un enfoque experimental [15], se utiliza
más importantes hechos desde Latinoamérica. La sección de comúnmente una cámara IR que genera imágenes en el domino
resultados y análisis presenta información que permite espacial y temporal a partir de las cuales es posible aproximar
identificar los investigadores y ejecutores en el desarrollo del la profundidad o grosor de los defectos [4]. En la mayoría de
TNDT; y por último una sección de conclusiones que muestra las aplicaciones para TNDT no es necesario una medición
cómo el enfoque actual para procesamiento de información absoluta de temperatura, la cual generalmente se encuentra en
térmica se orienta al procesamiento de secuencias de imágenes, ejercicios de inspección con termografía pasiva. Por el
donde el número de técnicas está en crecimiento y se evidencia contrario, sólo se requiere realizar mediciones relativas usando
la necesidad de contar con métricas robustas para evaluar su el enfoque activo de la técnica [16]. En el marco de este
desempeño. procedimiento se resalta la importancia de conceptos como la
sensibilidad, la precisión, la resolución espacial y la frecuencia
II. CONCEPTOS Y EVOLUCIÓN DEL TNDT de adquisición. Estos términos se definen a continuación:
Las primeras patentes relacionadas con la identificación de  Sensibilidad Térmica (NETD, del inglés: noise
puntos u objetos calientes se documentan como evidencia de los equivalent temperature difference): ruido equivalente a
inicios del TNDT, al surgir de necesidades específicas a nivel la menor temperatura detectable. Cuanto más bajo es el
industrial, las cuales involucraban el análisis de propiedades NETD mejor se podrá detectar un contraste térmico
térmicas de materiales (principalmente metales). A mediados significativo observable - sobreponiéndose al ruido del
de los sesenta, se reporta la implementación de uno de los sistema.
primeros ejercicios de TNDT, que fue realizado a la carcasa del  Precisión: es la medida de la confiabilidad en la
motor de un proyectil y donde se identificaron problemas estimación de temperatura aparente en relación con la
relacionados con las propiedades de emisividad de los objetos temperatura verdadera. Normalmente depende de la
pertenecientes al sistema inspeccionado. Debido a la carrera tecnología del dispositivo de captura: en cámaras
bolométricas de tipo Arreglo de Plano Focal (FPA, del
espacial, el método se vuelve atractivo para investigadores
inglés: focal plane array) se encuentra entre 50−100 mK
relacionados con esa industria, y de forma paralela la Sociedad
y en detectores cuánticos refrigerados aprox. 20 mK.
Americana de Pruebas No Destructivas (ASNT, del inglés:
Casi todas las cámaras tienen una precisión de ± 2 °C ó
american society for nondestructive testing) establece el primer ± 2 % para rangos más amplios.
comité que daría origen a publicaciones periódicas en el área.  Resolución espacial: campo de visión instantáneo que
En su teoría, la evaluación térmica involucra la solución directa representa el área cubierta por un sensor remoto en un
e inversa de problemas de transferencia de calor, donde determinado momento. Se presenta en términos del
variaciones de las propiedades térmicas y la presencia de Campo de Visión Instantáneo (IFOV, del inglés:
defectos internos en el material inspeccionado generan instantaneous field of view) con unidades de
disturbios dinámicos sobre el patrón térmico del material de miliradianes y depende del tamaño del detector y de la
interés [4]. Hasta ese entonces el enfoque cualitativo del óptica del dispositivo de captura.
método limitaba su éxito frente a otras técnicas de inspección.  Frecuencia de la imagen: velocidad a la cual el
La inclusión de elementos de teoría de conducción de calor instrumento de captura puede adquirir imágenes. Puede
proporcionó un nuevo nivel al TNDT, lo cual fue explotado y ser de cientos de Hz, para tecnologías basadas en
desarrollado por autores como Balageas [5], Vavilov & Taylor detectores refrigerados y generalmente con la
[6], MacLuaghlin & Mirchandani [7], Popov & Karpelson [8] posibilidad de incrementarse si la cámara ofrece la
entre otros y quienes además en la década de los 80 introdujeron posibilidad de reducir el tamaño de las imágenes.
modelos de defectos 1D, 2D y 3D [9]; por último, desde la El material por inspeccionar se puede someter a diferentes
década de los 90 aparecen los aportes de Maldague [10]–[12]. procesos de estimulación de energía que son descritos
El TNDT incluye generalmente la estimulación térmica de la generalmente por funciones matemáticas en el dominio
muestra a inspeccionar registrando su respuesta térmica en temporal (impulso, pulso rectangular, sinusoides). Los
función del tiempo (enfoque activo). Este procedimiento generadores de estímulos térmicos de naturaleza óptica o
incrementa su complejidad al considerar más variables, aunque fototérmica (flash, lámparas halógenas o Xenón, láser, etc.) son
en primera instancia parezca sencillo. La obtención de los más utilizados. Las propiedades del material inspeccionado
temperaturas absolutas con una cámara Infrarroja (IR, del afectan en gran medida la cantidad de energía que puede
inglés: infrared) necesita relacionar variables como la absorber, por lo que se convierten en la principal condición que
temperatura del objeto de interés, su distancia con respecto a la limita la aplicación de cada alternativa de excitación [6, 45].
cámara, la humedad del ambiente, la temperatura del medio Las técnicas de TNDT se clasifican en función de cómo se
ambiente y de la atmósfera, la emisividad del objeto, el índice obtienen las imágenes y cómo se procesan. Las cuatro
de transmisión de la atmósfera, así como la cantidad de principales técnicas son: termografía pulsada, termografía
partículas de polvo presentes en el ambiente [13], [14]. modulada, vibrotermografía, encontrándose además la
Reuniendo los conceptos anteriores y complementándolos con termografía de fase pulsada que nace de la combinación de la
teorías fundamentales de transferencia de calor, tomó fuerza en
Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira. 17

termografía pulsada y la termografía modulada [2], [3]. Pocas perspectivas en el campo del modelado se encuentra como meta
veces los defectos en el material inspeccionado son fácilmente el desarrollo de un paquete amigable y flexible con el usuario
detectables en las imágenes producto de los métodos que permita resolver problemas típicos 3D en TNDT
anteriormente citados, por lo que es indispensable aplicar involucrando intercambio de calor no adiabático (lineal y no
técnicas de procesamiento de imágenes para complementar el lineal), calentamiento arbitrario, transformación de fase en
estudio de ellas [2]. Tanto las técnicas de procesamiento de sustancias, entre otros. En cuanto al procesamiento de
información térmica clásicas y recientes se enfocan en reducir información los esfuerzos se concentran en el desarrollo de
el ruido antes y después de las etapas de análisis, buscando algoritmos inversos más precisos y robustos mejorando la
mejorar el contraste de las imágenes y entregar datos valiosos capacidad para caracterización de defectos en forma, grosor y
que faciliten la detección y/o caracterización de defectos. La profundidad.
sumatoria espacial de pixeles y sumatoria temporal de imágenes La siguiente producción académica [2], [22]–[24] recopila
se ha utilizado para el procesamiento de imágenes IR [16], así una cantidad considerable de fuentes de información sobre
como métodos para corrección de píxeles anómalos, efecto TNDT. Destacando de forma general aspectos como:
vignetting y atenuación de ruido [17]. Muchos procedimientos terminología y modelos físicos, formulación matemática y
y técnicas se han desarrollado y documentado en los últimos soluciones analíticas, características básicas de TNDT,
años [18], entre las que se destacan: la normalización del problemas inversos para caracterización de defectos con
contraste, el contraste absoluto diferencial, la reconstrucción de TNDT, algoritmos para procesamiento de información y
señales termográficas, termografía de componentes principales, eficiencia, hardware para TNDT y aplicaciones generales y
termografía de mínimos cuadrados parciales [2]. específicas. Lo anterior permite proyectar que los adelantos a
Muchas de las técnicas actuales de TNDT surgen de TABLE I
combinar métodos clásicos con información de propiedades NÚMERO DE PUBLICACIONES POR FUENTE DE INFORMACIÓN
físicas de los materiales, técnicas de procesamiento de
FUENTE DE INFORMACIÓN
imágenes y reconocimiento de patrones. En algunos casos se
Año
aplican las técnicas de forma independiente y luego se fusiona
Scopus Web of Science JEA-2021
su información para el análisis posterior; en otros casos las
técnicas son integradas y ejecutadas de forma simultánea. 2021 71 1 2
2020 307 7 9
Como ejemplos de este tipo de fusión se reporta la combinación 2019 339 12 2
de la IRT con la tecnología de Radar de Penetración Terrestre 2018 292 9 3
(GPR, del inglés: ground penetrating radar), los rayos X, la 2017 221 5 4
interferometría láser, las imágenes digitales acústicas, el 2016 216 6 17
2015 171 6 11
ultrasonido y las imágenes 3D entre otras [2]. Umar [19] en su 2014 148 7 18
trabajo realiza una descripción sobre la potencialidad de la 2013 146 5 7
técnica de Termografía Infrarroja Ultrasónica (UIT, del inglés: 2012 98 6 9
2011 90 4 2
ultrasonic infrared thermography), la cual se ha utilizado 2010 63 4 3
durante la última década, principalmente por investigadores del
área de termo-mecánica. un mediano plazo se encontrarán enfocados en las áreas de
En [20], [21] se realiza una descripción sobre los principales adquisición, estimulación, procesamiento y rangos de
avances en hardware y software desde el año 1990 alcanzados aplicaciones con nuevos materiales, que impondrán nuevos
por La Universidad Politécnica de Tomsk donde la principal desafíos en el área.
tarea se ha concentrado en modelado avanzado y procesamiento
de información. De su trabajo se destacan productos como III. LITERATURA Y TENDENCIAS ACTUALES
Multilayer, ThermoCalc: 2D-2DM-2D Build-6L-36L Source -
En el ámbito mundial existen algunas organizaciones que
Mine (para simulación de problemas TNDT) y ThermoFit Pro
integran a la comunidad en el área de TNDT. Su origen se
(para procesamiento de secuencia de imágenes IR sintéticas y
remonta a 1978 cuando en Estados Unidos nace el Thermosense
experimentales), paquetes que son continuamente actualizados.
respaldado por la Sociedad Internacional de Óptica y Fotónica
También presentan adelantos en la determinación de
(SPIE, del inglés: international society for optics and photonics)
propiedades térmicas de materiales anisótropos, evaluación de
y el ITC-InfraMation, evento que convocó a profesionales de
corrosión de superficies traseras, termografía infrarroja
diferentes campos y permitió conocer e intercambiar
ultrasónica y por inducción, y análisis de fractura en materiales
información valiosa en el área, que inicialmente se encontraba
entre otros. Por último, a nivel de hardware se presentan
dispersa entre diferentes conferencias y publicaciones sin
adelantos en el desarrollo de unidades portátiles de TNDT para
conexión directa. La existencia de una conferencia transversal
inspección de corrosión en cilindros metálicos y para
en Norte América ayudó a subsanar estos vacíos, ejercicio que
inspección de aeronaves, con la finalidad de aplicar TNDT
en 1992 se replicó en Europa, dando origen a la conferencia
dentro y fuera de las condiciones habituales de un laboratorio,
QIRT, la que posteriormente en 2014 también fue
principalmente en los materiales compuestos (CFRP y Plásticos
institucionalizada en Asia, facilitando así la contribución de
Reforzados con Fibra de Vidrio (GFRP, del inglés: glass fiber
expertos de estos continentes. Al año 2016 es notable el alto
reinforced plastic)) y estructuras en forma de panal. Como
18 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

grado de participación de diferentes países encontrando un para el diseño de la red se utiliza un modelo térmico sencillo,
incremento cercano al 130 % en los últimos 22 años (Europa encontrándose una limitación experimental entre la frecuencia
51,4 %, Continente Americano 10.8 %, Asia 24.3 % y otras de muestreo y la conductividad térmica del material
regiones 13.5 %) [25]. Lo anterior brinda un panorama de inspeccionado, evidenciando que con la metodología propuesta
referencia sobre la contribución por regiones al TNDT, lo que, en muestras de aluminio es necesario frecuencias de 140 Hz
complementado a un ejercicio de búsqueda sobre algunas de las mientras que para muestras plásticas parece ser suficiente sólo
principales fuentes de información reconocidas, permite de 14 Hz. Las ANNs también se han utilizado en la detección
recopilar las cifras consignadas en la Tabla 11 confirmando una de defectos y estimación de su profundidad en materiales
tendencia creciente de producción académica en el área (figura compuestos, fenómenos que están relacionados con problemas
inversos que son de naturaleza compleja y no lineal, y que por
1 columna 2: procesamiento de secuencias de imágenes). Con
lo tanto han representado una de las mayores limitantes para la
base en lo anterior, la información publicada en el área de
ejecución de TNDT. No obstante, las ANNs se utilizan para
TNDT y recopilada bajo la etiqueta JEA-2021 se clasificó en
proporcionar información cuantitativa sobre los defectos. El
los siguientes grupos: técnicas de procesamiento, métricas y trabajo [23] presenta un modelo térmico 3D para materiales no
desempeño, modelos térmicos y simulaciones y, por último, homogéneos de tipo CFRP, implementando dos redes
aplicaciones. neuronales para detección de defectos y estimación de su
A. Técnicas de procesamiento profundidad. Por último, se confrontan y analizan los resultados
teóricos obtenidos a partir del modelo térmico y los resultados
En muy pocas ocasiones los defectos o información de
obtenidos con ayuda de las redes neuronales.
interés presente en las imágenes térmicas de los materiales
Con el objetivo de extraer características y reducir la
inspeccionados a partir de TNDT es fácilmente evidenciable.
información redundante en secuencias de imágenes térmicas,
Además, las imágenes contienen implícitos efectos no deseados
como técnica para el procesamiento, es posible usar el Análisis
como vigneting, distorsión radial, reflexiones ambientales,
por Componentes Principales (PCA, del inglés: principal
variaciones de emisividad, calentamiento no uniforme, ruido,
component analysis), que proyecta los datos originales en un
etc. [22]. Para eliminar la información no deseada es necesario
sistema de componentes ortogonales, surgiendo la técnica PCT
el procesamiento de dichas imágenes. A través de los años se
[26].
han desarrollado múltiples algoritmos de procesamiento, cuya
En Colombia, en la Universidad del Valle se han realizado
selección depende del objetivo de la aplicación o campo de
avances a nivel software para el análisis de imágenes IR. En
investigación. En la figura 1 se presenta una reciente
[27] se presenta un sistema de visión artificial para imágenes
clasificación realizada por [113] donde se encuentran los
infrarrojas el cual permite el procesamiento y análisis de esas
métodos de procesamiento más utilizados en TND
imágenes a través de una interfaz gráfica desarrollada en JAVA,
En [12] se muestra la extracción de información cuantitativa
con la posibilidad de generar reportes de inspección y el
sobre resultados de la PPT, haciendo uso de la capacidad de las
almacenamiento de los resultados en una base de datos.
Redes Neuronales Artificiales (ANNs, del inglés: artificial
Posteriormente, la aplicación TermoIR [18] con un mayor
enfoque a TNDT permite la adquisición, procesamiento,
análisis y caracterización de información térmica. Por último,
en [28], se presenta un desarrollo software para fusión
multimodal de imágenes térmicas y de espectro visible.
En Estados Unidos, la excitación óptica a partir de flash de
potencia es el método más empleado en TNDT, ya que es rápido
en comparación con otros métodos, y facilita la estimación
cuantitativa de parámetros como grosor, profundidad o
difusividad térmica. Se ha demostrado que la técnica de
Reconstrucción de la Señal Termográfica (TSR, del inglés:
thermographic signal reconstruction) [29] es efectiva para la
reducción de ruido y, adicionalmente, su primera y segunda
derivada proveen información que es invariante a condiciones
ambientales, preparación de la superficie o energía de entrada;
ayudando a evidenciar la presencia de anomalías al interior de
las superficies inspeccionadas [30].
Fig. 1. Técnicas de procesamiento utilizadas en TNDT. Para mejorar las diferencias térmicas entre regiones sanas y
defectuosas de un material se han propuesto múltiples
definiciones de contraste térmico (absoluto, dinámico,
normalizado, estándar), las cuales necesitan utilizar como
neural network) para resolver problemas complejos y no referencia valores de temperatura de una región sana del
lineales a partir de información parcial y ruidosa. En este caso, material a inspeccionar, situación que en la práctica no es fácil

1
Ecuaciónes de búsqueda utilizadas: i) Scopus: ALL(“NDT” AND iii) para la información recopilada bajo la etiqueta JEA-2021 se utilizaron
“Thermography” AND “Composite materials”), ii) Web of Science: from All fuentes de información mixtas (Scopus, Web of Science, ScienceDirect,
Databases TS=(“NDT” AND “Thermography” AND “Composite materials”) y IEEEXplore, ResearchGate, QIRT Conference Open Archives, etc.).
Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira. 19

de identificar y que limita su aplicación. Este inconveniente fue efecto térmico en los materiales es de naturaleza tridimensional
resuelto a partir de la definición del Contraste Absoluto (3D). En las cámaras térmicas modernas el tamaño de las
Diferencial (DAC, del inglés: differential absolute contrast) imágenes se ha incrementado notablemente manteniendo una
[31], que recurre a la solución unidimensional de la ecuación de tendencia creciente, ante lo cual es necesario utilizar modelos
Fourier para materiales homogéneos y semi-infinitos unidimensionales para garantizar la implementación de algunos
estimulados por una función impulso; no obstante, el modelo métodos de procesamiento de la información térmica. A partir
no incluye el espesor de la muestra provocando que el método de lo anterior, con el objetivo de conciliar esta realidad y lograr
DAC pierda exactitud después de haber transcurrido un tiempo un procesamiento de secuencias de imágenes de una manera
considerable después del impulso de excitación. Con base en ágil y eficaz, se plantea un nuevo enfoque para la posible
este método Benítez et al. [32] mejoró esta propuesta al incluir detección y caracterización temprana. Enfoque que se ilustra en
en el análisis el grosor de la muestra a inspeccionar y el ancho varias aplicaciones y es comparado con métodos clásicos de
del pulso térmico, permitiendo extender el rango de validez del análisis de información térmica, demostrando que la
DAC a tiempos mayores después de aplicado el pulso de simplicidad de los modelos y cálculos es compatible con
excitación y generando así la versión modificada del DAC procesos de identificación exactos y eficientes.
(mDAC, del inglés: modified differential absolute contrast). Una técnica de procesamiento para mejorar el contraste
En [33] se emplean estadísticas de alto orden para procesar térmico de una secuencia de imágenes térmicas adquiridas a
imágenes IR comprimiendo la información disponible a partir partir de un ejercicio experimental de termografía activa se
de la secuencia típica en una sola imagen por cada inspección. presenta en [18]. El procedimiento puede compensar el efecto
Adicionalmente, se hace énfasis sobre la importancia del óptico de vigneting y el efecto de calentamiento no uniforme de
análisis automático de imágenes obtenidas en TNDT, la superficie producido por la fuente de excitación de energía
reduciendo la intervención humana y posibles fuentes de error inyectada a la muestra inspeccionada utilizando un filtro
propias de su subjetividad. mediana y sin depender de modelos de propagación de calor.
Restrepo y Loaiza [34] presentan un algoritmo para la Algunos de los resultados reportan una compensación efectiva
detección y corrección de píxeles anómalos en secuencias de por efectos de calentamiento no uniforme en comparación con
imágenes térmicas resultantes de aplicar un procedimiento de técnicas como DAC, alcanzando una reducción de no
termografía activa. Al contrario de técnicas clásicas de filtrado uniformidad de fondo de hasta un 74 % a partir de un parámetro
espacial, esta propuesta no afecta la apariencia global de cada de no uniformidad definido bajo este estudio. El método
imagen, ayudando a reducir errores de interpretación en etapas inicialmente propuesto se complementa posteriormente
de procesamiento posteriores. realizando el truncamiento de valores de contraste negativos,
Marinetti y Cesaratto [35] documentan un método para la alcanzando un desempeño por encima de 11 dB en Relación
estimación de emisividad en un material que se fundamenta en Señal Ruido (SNR, del inglés: signal-to-noise ratio) con
funciones de calibración individuales para un sensor IR respecto a las secuencias de imágenes analizadas [39].
específico. La temperatura reflejada se mantiene en un valor A partir del conocimiento de los múltiples avances en
constante durante la fase experimental sin necesidad de herramientas de pre y pos procesamiento desarrollados en las
estimarla o conocer su valor. Finalmente se realiza una dos últimas décadas para mejorar el desempeño de la
discusión sobre los resultados experimentales obtenidos en los termografía pulsada, en la detección y caracterización de
valores de temperatura absoluta y relativa, para diferentes tipos defectos, trabajos como [40], [41] reconocen a la TSR y la
de materiales usados comúnmente como referencia de detección temprana (del inglés: early detection) como dos de las
emisividad. más eficientes. Describen cómo los dos métodos son
En el trabajo de Ramirez, García y Castellanos [36] se complementarios y pueden ser aplicados también cuando se
implementa un algoritmo de segmentación en imágenes IR de utiliza la SHT. Teniendo en cuenta que la TSR se basa en el
un motor eléctrico, basado en el agrupamiento de información ajuste polinomial logarítmico de la evolución en el tiempo de la
a partir de maximizar los valores esperados (esperanza), y en temperatura, Roche, Leroy y Balageas [42], [43] resaltan el uso
donde el desempeño se realiza a partir de su comparación con de esta técnica para incrementar la eficiencia de la PT; aunque
otro método de segmentación basado en fronteras usando como no permita el análisis de secuencias completas de imágenes, lo
parámetro de evaluación el coeficiente de Dice. que es necesario en algunas aplicaciones como la inspección de
Ibarra, Tarpani y Maldague [37] realizan una diferenciación estructuras. Con base en la limitación expuesta anteriormente,
entre técnicas básicas de procesamiento de imágenes como se ha propuesto un nuevo enfoque de pos procesamiento [42]
promedio, substracción, división, filtrado, detección de bordes que sólo hace uso de los tres coeficientes polinomiales más
y técnicas avanzadas de procesamiento que son comúnmente relevantes entregados por el método TSR clásico, que además
aplicadas en etapas de pre y pos procesamiento en TNDT. Los mejora los tiempos de cómputo y produce una única imagen
métodos avanzados son clasificados en diferentes grupos de compuesta de defectos y sus profundidades, imagen
acuerdo con el manejo de la información, resaltando aquellos representada en el modelo (RGB, del inglés: red, green, blue).
basados en: contraste térmico, estadísticas, de factorización de Específicamente, en [43] se analizan secuencias de imágenes de
matrices, transformadas y por último los que utilizan materiales compuestos y de estructuras con daños reales,
inteligencia artificial (redes neuronales y computación producto del proceso de fabricación o de impactos, a través de
evolutiva). la generación de las imágenes RGB resultantes de la versión
Balageas [38] presenta que la técnica de termografía activa reducida de la TSR, y se comparan los resultados con las
permite caracterizar propiedades termo físicas y defectos en técnicas de Termografía de Fase Pulsada (PPT, del inglés:
estructuras complejas que presentan heterogeneidad, donde el pulsed pase thermography) [44], Termografía de Componentes
20 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

Principales (PCT, del inglés: principal component no uniforme y difusión lateral del calor. Lo anterior se logra
thermography) [45] y con Estadísticos de Orden Superior realizando el proceso de división de la secuencia de imágenes
(HOS, del inglés: high order statistics). experimentales por una secuencia de imágenes sintéticas
Vavilov y Chulkov en su trabajo [46] confirman la tendencia obtenida a partir de un modelo de transferencia de calor 3D.
creciente en el campo del procesamiento de información para Zhang y Maldague [52] describen la implementación de una
TNDT en los dominios temporal y espacial, gracias al cual es Transformada de Contorno sin Submuestreo (NSCT, del inglés:
posible atenuar considerablemente muchos de los factores que nonsubsampled contourlet transform) que se soporta bajo un
producen ruido. En ese orden de ideas, se presenta el desarrollo enfoque multimodal de fusión de imágenes (IR y espectro
de dos técnicas cuantitativas para TNDT: i) mapas de inercia visible). En sus resultados experimentales y a partir de la
térmica y ii) difusividad térmica, métodos que buscan percepción humana, el proceso de fusión propuesto evidencia
establecer relaciones cuantitativas entre la severidad de su efectividad y ventajas en relación con aspectos como escalas
defectos encontrados en una inspección y parámetros para toma múltiples, resolución múltiples e invariancia a la traslación,
de decisiones. frente a métodos de análisis como la transformada Wavelet.
Restrepo y Loaiza [47] introducen dos enfoques de Adicionalmente, su eficiencia y robustez se analizan a partir de
discretización por diferencias finitas para el modelo de diferentes indicadores como desviación estándar, entropía de
propagación de calor 3D de Fourier, con lo cual se propone una Shannon, error cuadrático medio, índice de información mutua
nueva técnica para el mejoramiento del contraste térmico de e índice de similitud basada en bordes.
secuencias de imágenes térmicas adquiridas usando PT en
muestras de CFRP. Los dos enfoques se apoyan en modelos
B. Métricas de desempeño
adaptables a una estructura de filtro espacial que puede ser
aplicable a cada imagen térmica con el fin de obtener un mejor Para evaluar las técnicas de procesamiento de información
contraste térmico entre las posibles fallas y regiones sanas del térmica es importante disponer de indicadores de desempeño
material inspeccionado, generando así una mejor probabilidad que permitan cuantificar sus resultados con criterios objetivos.
para detección de defectos. El desempeño de la técnica A continuación se describen algunos de los indicadores de
propuesta es evaluado usando secuencias térmicas sintéticas desempeño que se encontraron en la literatura científica
generadas con ThermoCalc6L y con secuencias de imágenes consultada:
reales de láminas de CFRP. Se obtiene un mejor contraste entre Teniendo en cuenta la existencia de Ruido No Uniforme (NUN,
defectos y el fondo de la imagen comparado con mDAC, y una del inglés: nonuniformity noise) en los sensores de captura de
ejecución más rápida que las técnicas que se apoyan en una las cámaras IR, San Martin, Torres y Pezoa [53] proponen el
operación de reconstrucción de la distribución térmica del índice de Aspereza Efectivo (ERo, del inglés: effective-
material como el método de filtrado térmico 3D [48]. Sobre esta roughness) que no necesita una imagen de referencia para
línea Restrepo [49] propone un método robusto para evaluar el desempeño de métodos de corrección de no
caracterizar la profundidad de defectos contrastados a partir del uniformidad propuestos en la literatura.
mejoramiento de contraste en secuencias de imágenes térmicas. Vavilov [54] reporta que en TNDT uno de los principales
La técnica propuesta entrega buenos resultados para defectos problemas inherentes a la técnica es la presencia de ruido en la
superficiales y con mejor contraste. Sin embargo, disminuye su señal IR, donde se introducen dos tipos de ruido: (1) de tipo
confiabilidad para aquellos que se encuentran a una mayor aditivo, definido por las reflexiones del entorno y el detector IR;
profundidad en el material y con menor contraste. Los y (2) de tipo multiplicativo, que se atribuye a características
resultados evidencian una reducción del error relativo promedio propias del material en términos de variaciones en su
de aproximadamente la mitad. emisividad. La realización de varios experimentos a diferentes
Akhloufi y Verney [50] presentan avances en tecnologías 3D materiales ha mostrado que en ninguno de estos se encuentra un
que facilitan el acceso a información de alta precisión a costos ruido aditivo puro o un ruido multiplicativo puro. Sin embargo,
razonables. Las técnicas de visión 3D permiten la inspección con calentamiento óptico, algunos materiales compuestos con
avanzada, el análisis y la metrología en procesos de superficies de alta emisividad (generalmente negras)
manufactura. Estas características combinadas con los evidenciaron ruido multiplicativo con un índice de contraste de
desarrollos comprobados del TNDT para detección de defectos ruido entre 2 y 5 %. Este índice se define como la relación entre
debajo de una superficie, genera un amplio potencial en la la desviación estándar del ruido y el valor medio de la
adquisición y fusión multimodal de información. El trabajo temperatura en una zona no defectuosa del material
documenta resultados prometedores, mostrando cómo estas inspeccionado.
nuevas técnicas pueden ser aceptadas e implementadas de Por su parte, en [55] se afirma que normalmente la
forma confiable y eficiente en procesos de manufactura a nivel evaluación de los métodos de procesamiento dependen del
aeroespacial y automotriz. concepto subjetivo emitido por el personal experto, quien se
Pawar y Vavilov [51] expresan que las falsas detecciones encarga de interpretar las imágenes y determina con qué técnica
relacionadas a las propiedades ópticas no uniformes sobre un se evidencian de forma visible los defectos. La SNR es un
material inspeccionado con termografía pulsada son un indicador cuantitativo relacionado con la visibilidad de los
problema típico que se puede atenuar haciendo uso de técnicas defectos, pero es un parámetro local y variable en función del
de normalización. Para este caso se presenta una versión 3D del tiempo. En ese sentido el autor propone un método para la
método conocido como normalización 1D (división de toda la Evaluación de la Calidad de Imagen (IQA, del inglés: image
secuencia de imágenes térmicas por una imagen), presentando quality assessment) soportado en métricas que permiten
como resultados la reducción de los efectos por calentamiento
Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira. 21

cuantificar el mejoramiento del contraste obtenido después de transferencia de calor en 3D sobre un medio heterogéneo
aplicar técnicas de procesamiento de imágenes en TNDT. cuando un pulso térmico es aplicado a la superficie del material.
García, Benitez, Castellano y Ramirez [56] proponen un A partir de la estimulación por ultrasonido y óptica, se han
algoritmo de evaluación de calidad que no requiere definir generado estudios que documentan ejemplos de caracterización
referencias (regiones o imágenes) para la evaluación de de daños por impacto, presentando la solución analítica de un
desempeño de diferentes técnicas de procesamiento en TNDT. problema termo mecánico de fricción interna en fisuras de
Con el estudio se encuentra una fuerte correlación entre el grosor casi nulo [63]. Posteriormente en [64] se muestran
método de evaluación propuesto y la evaluación subjetiva soluciones numéricas y analíticas de modelos donde los
realizada por un humano. Adicionalmente, se reporta una materiales defectuosos han sido objeto de estimulación térmica
relación con la densidad espectral de potencia espacial y e incluyen estudios de comparación entre modelos de
algunos coeficientes Wavelet. transferencia de calor 1D, 2D y 3D para definir el nivel de
Lopez, Maldague e Ibarra [57] examinan diferentes métodos exactitud de las soluciones.
de procesamiento: DAC, TSR, PPT y PLSR. El desempeño de En [65], se utiliza el criterio Tanimoto y la información de
los métodos se analiza en términos de la relación señal ruido y simulaciones por elementos finitos con el software COMSOL,
el contraste máximo. Este trabajo se complementa en [57] y para analizar imágenes térmicas obtenidas usando diferentes
usando los indicadores de desempeño del trabajo previo, se medios de excitación (pulsada, modulada, vibrotermografía e
evalúan las técnicas TSR y PLSR aplicadas a una muestra de iluminación con tecnología de Diodos Emisores de Luz (LED,
CFRP con inserciones de Teflón para generar los defectos. Los del inglés: light emmiting diode)). Lo anterior con la finalidad
autores discuten sobre el rendimiento de los métodos, sus de detectar inserciones en Teflón sobre láminas de CFRP. Por
limitaciones y sus ventajas. También se mencionan los último, en [66] se desarrolla un enfoque inverso para el
principales aspectos que afectan el ejercicio de la PT problema de transferencia de calor con la finalidad de
(reflexiones ópticas, variaciones de emisividad, calentamiento caracterizar el tamaño y profundidad de defectos
no uniforme y difusión de calor no uniforme). subsuperficiales en materiales compuestos. El procedimiento se
En [58] se realiza una revisión de los alcances obtenidos con basa en la reconstrucción del perfil de fase usando una
la TSR para la detección de defectos en estructuras, simulación por elementos finitos de un experimento de LT. La
comparando su desempeño con métodos de procesamiento técnica aplica un método híbrido entre un algoritmo de
reconocidos ampliamente en la literatura como: PPT, PCT y simulación y otro de búsqueda para determinar el tamaño y
HOS. Adicionalmente, ante la necesidad de realizar una profundidad de los defectos en láminas de CFRP. Lo anterior,
comparación objetiva sobre la efectividad de estos métodos de encontrando una solución óptima para una función objetivo
procesamiento usados en TNDT, Balageas, Roche y Leroy [59] construida con información del perfil de fase del experimento
utilizan como métricas de desempeño la relación contaste ruido, LT y una solución numérica. Con esta estrategia se reportan
la forma de los bordes y la extensión lateral de los defectos, errores menores al 5 % y 4 % para tamaño y profundidad,
teniendo en cuenta que algunos trabajos previos sólo respectivamente.
mencionan variaciones de la relación señal ruido presente en las
D. Aplicaciones
imágenes térmicas. Las definiciones anteriores se aplican
seleccionando manualmente filas o columnas de píxeles que se Los rangos de aplicación relacionados con TNDT son
encuentren en los centros de los defectos identificados. notablemente amplios y en expansión. Comprenden áreas como
Por último, Liu, Zhang, Fernandes y Maldague [60] la medicina y la biología, la electrónica, los semiconductores,
presentan la aplicación de diferentes métodos de procesamiento la inspección de patrimonio histórico, las aplicaciones
de imágenes: sustracción de imagen fría, PCT, TSR y la industriales, los procesos de fabricación y control de calidad,
transformada de Fourier a imágenes térmicas obtenidas a partir así como tareas de mantenimiento, entre otros. A continuación,
de vibrotermografía y PT. La discusión de los resultados se presentan algunos ejemplos de trabajos asociados con
obtenidos se realiza en función de los métodos de excitación algunas áreas de aplicación del TNDT.
utilizados. En obras civiles es viable la detección de anormalidades a
partir de una revisión del comportamiento y estado térmico de
C. Simulación y modelos térmicos las estructuras [67], relacionando el estado transitorio de
Los modelos analíticos y numéricos son comúnmente tensión y estrés del material inspeccionado ocasionado por la
utilizados en simulaciones relacionadas con TNDT, por lo que relación existente entre fenómenos mecánicos y efectos
a continuación se relacionan algunos de los trabajos que se térmicos. Bajo otro enfoque, más relacionado con el TNDT
consultaron. clásico, los daños se detectan, localizan y caracterizan por las
En [61] se compara el desempeño de cuatro paquetes perturbaciones que se introducen en las propiedades térmicas
software para simulación numérica de problemas en TNDT locales de la estructura (conducción, difusión, capacidad de
(ThermoCalc-2D, ThermoCalc-6L, ThermoModel y Comsol calor, efusividad, resistencia térmica). En los últimos años han
Multiphysics), en términos de cálculos de temperaturas aparecido nuevas opciones de materiales para tuberías,
absolutas y algunos parámetros térmicos (señales de principalmente basados en polímeros compuestos, sustituyendo
temperatura diferenciales y contraste dinámico). así los materiales tradicionales, ante lo cual es posible realizar
En [62] se realiza la descripción de una metodología para la inspección de tuberías con TNDT [68]. Por otra parte,
modelar y simular un ejercicio de termografía pulsada. Se Milovanovic y Pecur [69] realizan un recorrido a través de los
utiliza un modelo térmico basado en la ecuación de principales avances en el uso de TNDT para detección y
caracterización de defectos bajo el enfoque activo sobre
22 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

concreto reforzado. También, presentan una revisión acerca de disminución, el número de publicaciones del año 2020 (307
los principios físicos, equipos de inspección y métodos de artículos, 14.2 %) sigue siendo mayor que el alcanzado en el
procesamiento involucrados. año 2018 (292 artículos, 13.5 %), fecha anterior al pico
En el área biomédica algunas aplicaciones derivan sus máximo.
resultados del análisis de imágenes térmicas; por ejemplo, en la
identificación biométrica utilizando imágenes IR de la red
vascular de la cara dorsal de la mano [70]. El uso de IRT
representa una herramienta valiosa para detectar síndrome de
Túnel Carpiano [71], al aplicar técnicas para registro y
extracción de características de las imágenes térmicas que
permite realizar una comparación entre sujetos sanos y aquellos
que presentan el síndrome. De manera similar en [72] se
presenta la aplicación del método TSR para observar el
comportamiento térmico en regiones vasculares de un ser
humano, encontrando diversas evidencias sobre los patrones de
calentamiento entre venas y válvulas venosas, permitiendo así
su identificación y funcionalidad. También se han logrado
avances como apoyo en la detección de anormalidades en las
imágenes IR de zonas mamarias [73], [74]. Y se han
desarrollado sistemas de visión por computador
multiespectrales en tiempo real para el análisis morfo-térmico
Fig. 2. Publicaciones en TNDT por año a nivel mundial entre 2010 y 2021.
de la huella plantar [75].
En aplicaciones industriales, Roche y Balageas [76]
describen una implementación práctica en una línea de Las figuras 3 y 4 muestran la producción académica
producción para representar mapas de profundidad de defectos organizada por país en la misma ventana de tiempo (2010-
al aplicar TNDT para la inspección de materiales compuestos. 2021). Para este ejercicio se utilizó como fuente de información
El resultado es similar a los producidos por un escáner los registros etiquetados como JEA-2021 (Tabla 1 columna 4).
ultrasónico, permitiendo la detección de daños en el material En estas gráficas se identifica la contribución que diferentes
originados por defectos de fabricación o por eventos fortuitos. naciones a nivel mundial están realizando en el área de TNDT.
También se encuentran trabajos relacionados con el uso de Con barras se representa la contribución en publicaciones por
termografía para la medición de temperatura de bombillos LED Autor Principal (AP). La línea continua complementa el
de alta definición, comúnmente utilizados en los faros indicador anterior incluyendo las coautorías (Co). A partir del
delanteros de los automóviles. Las lecturas térmicas resultantes primer criterio (AP), se identifican países como Canadá,
son comparadas con datos de temperatura generados a partir de Francia y Rusia con 15, 13 y 11 publicaciones, respectivamente
simulaciones [77]. Por otra parte, se ha propuesto utilizar (Figura 3). Usando el mismo criterio, únicamente se identifica
termografía activa conjuntamente con métodos de en América Latina la contribución de Colombia, Brasil, Chile y
procesamiento de información térmica para la inspección de México con 17, 5, 2 y 1 publicaciones, respectivamente (Figura
tuberías de acero y estimación de la cantidad de agua remanente 4). Con el segundo criterio, incluyendo las coautorías (AP +
en estas estructuras [78]. Por último, en [79] se expone un Co), Canadá incrementa sus registros a 59, Francia a 28 y Rusia
procedimiento para la reconstrucción de impedancia térmica a a 21 (Figura 3). Por su parte, en América Latina, Colombia
partir de mediciones infrarrojas de la superficie de un
dispositivo electrónico.

IV. RESULTADOS Y ANÁLISIS

Con la revisión de la literatura y la organización de la
información relacionada en las secciones anteriores, es posible
documentar algunas tendencias a nivel de publicaciones en el
área de TNDT, así como la producción académica por país y
por autor en la última década, obteniendo así los resultados
registrados en las figuras 2, 3, 4, 5 y 6.
La figura 2 presenta, en una ventana de tiempo entre el año
2010 y el 2021, el número de publicaciones a nivel mundial en
el área de TNDT usando como fuente de información Scopus y
la ecuación de búsqueda que acompaña la Tabla 1 (columna 2).
La tendencia creciente en el número de registros coincide con
la reportada en [25]. El número máximo de artículos publicados
se presenta en el año 2019 con 339 publicaciones lo cual
equivale al 15.7 % del total de registros dentro de la ventana de Fig. 3. Número de publicaciones en TNDT por país entre 2010 y 2021 en el
observación; y aunque para el siguiente año se aprecia una leve contexto mundial.
Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira. 23

Florez continúa con 2, lo que significa que sólo ha contribuido

como AP (Figura 6). El nivel de participación de la mayoría de
los investigadores como AP (gráfico de barras) se incrementa
cuando la información se complementa con las coautorías
(información que se representa en el gráfico como una línea
continua).

Fig. 4. Número de publicaciones en TNDT por país entre 2010 y 2021 en el

contexto de América Latina.

suma 50 registros en publicaciones académicas, seguido de

Brasil con 10, Chile con 2 y por último México con 2 (Figura
4). Con base en la información anterior, es notoria una mayor Fig. 5. Número de publicaciones en TNDT por autor entre 2010 y 2021 en el
participación de autores con afiliación canadiense, francesa o contexto mundial.
rusa con respecto al resto de países. Estos investigadores han
permitido que otros científicos de países emergentes en el
campo de TNDT, como es el caso de América Latina,
contribuyan y propongan desarrollos en los últimos años.
También es evidente el pobre desarrollo y contribuciones
encontrados en el resto de los países que conforman América
Latina. Lo anterior, posiblemente se debe a que la ecuación de
búsqueda utilizada en esta investigación haya excluido
evidencias de la producción académica, que efectivamente no
existan desarrollos publicados en el área de TNDT, que no se
están realizando ejercicios bajo este enfoque o que las prácticas
en las otras regiones del continente americano en términos de
NDT aún no involucren o promuevan la termografía infrarroja
(activa).

Las figuras 5 y 6 presentan la producción académica

organizada por autor para la misma ventana de tiempo (2010- Fig. 6. Número de publicaciones en TNDT por autor entre 2010 y 2021 en el
contexto de América Latina.
2021) y utilizando la misma fuente de información que el
ejercicio anterior, los registros etiquetados como JEA-2021
(Tabla 1 columna 4). De igual forma, con barras se representa
la contribución en publicaciones por (AP) y la línea continua
complementa el indicador anterior incluyendo además las
coautorías. Como autor principal aparecen investigadores V. CONCLUSIONES
ampliamente conocidos en la literatura académica relacionada Actualmente se presenta una tendencia creciente en producción
con TNDT. A nivel mundial se destacan nombres como académica relacionada con ensayos no destructivos por
Vavilov, Balageas, Maldague, Roche, Marinetti e Ibarra con 10, termografía infrarroja. Esta tendencia confirma que es un área
9, 6, 4, 4 y 2 publicaciones (Figura 5). Por su parte, a nivel de de conocimiento vigente y con un amplio potencial de
América Latina aparecen investigadores como Benítez, desarrollo en el futuro inmediato. Además, con alto interés para
Restrepo, Lopez, Erazo, Flórez y Fernandez con 6, 6, 4, 3, 2 y diferentes industrias, principalmente en sectores como el
1 artículos (Figura 6). Incluyendo las coautorías (AP + Co), se militar, el aeroespacial, el naval, el automotriz, en biomédica y
observa que investigadores como Maldague incrementan el
más recientemente en salud estructural.
número de su producción a 29, Ibarra y Vavilov a 15, Balageas
a 13, Roche a 9 y Marinetti a 5 (Figura 5). A nivel de América
Existe una división bien definida para clasificar los métodos
Latina, aparece Loaiza con 14 artículos en coautoría, Benítez y
utilizados en el procesamiento de información térmica. La
Restrepo con 11, López con 6, Erazo con 4, Fernándes con 2 y
24 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

división consta de dos grupos: el primero que reúne todos las [5] D. L. Balageas, J. C. Krapez, and P. Cielo, “Pulsed photothermal
modeling of layered materials,” J. Appl. Phys., vol. 59, no. 2, pp. 348–
técnicas orientadas al procesamiento de imágenes individuales
357, Jan. 1986, doi: 10.1063/1.336690.
y el segundo grupo que está conformado por los métodos [6] V. P. Vavilov and R. Taylor, “Theoretical and practical aspects of the
orientados al procesamiento de secuencias de imágenes thermal NDT of bonded structures,” Res. Tech. NDT, vol. 5, no. N.A.,
térmicas. Actualmente, este segundo conjunto es el que pp. 239–280, 1982.
[7] P. V. MacLaughlin and H. G. Mirchandani, “Aerostructure NDT
presenta un desarrollo mayor incluyendo trabajos relacionados
evaluation by thermal field detection (Phase II). Final Rep.,”
al mejoramiento de contraste térmico, estadísticas de orden AIRTASK, Nav. Air Syst. Command AIR- 310G, vol. N.A., no. N.A.,
superior, factorización de matrices, transformaciones p. N.A, 1984.
matemáticas y el uso de técnicas avanzadas de inteligencia [8] Y. A. Popov, K. A. E., V. A. Strokov, and E. Al., “Thermal NDT of
multi-layer structures,” Defectoscopiya (Russ. J. Nondestr. Test.),
artificial [23][46]. Entre las técnicas más reconocidas y citadas
vol. 3, no. N. A., pp. 76–81, 1976.
por los autores consultados se encuentran la reconstrucción de [9] V. Vavilov, “Thermal Non-Destructive Testing: Short History, State-
señales termográficas, el contraste absoluto diferencial, la of-the-Art and Trends,” Image Process., 1995.
termografía de fase pulsada, la termografía de componentes [10] X. Maldague, J. C. Krapez, and D. Poussart, “Thermographic
Nondestructive Evaluation (NDE): An Algorithm for Automatic
principales, la reconstrucción por mínimos cuadrados parciales
Defect Extraction in Infrared Images,” IEEE Trans. Syst. Man
y la detección temprana. Por su parte, autores de origen Cybern., vol. 20, no. 3, pp. 722–725, 1990, doi: 10.1109/21.57287.
latinoamericano han propuesto técnicas como el contraste [11] X. Maldague and S. Marinetti, “Pulse phase infrared thermography,”
absoluto diferencial modificado, el contraste térmico por J. Appl. Phys., vol. 79, no. 5, p. 2694, 1996, doi:
[Link]
diferencias finitas, la compensación térmica de fondo por
[12] X. Maldague, Y. Largouët, and J.-P. Couturier, “A study of defect
filtrado y por reconstrucción. A la par de lo anterior, se depth using neural networks in pulsed phase thermography:
identifica que existe una tendencia para que las nuevas técnicas modelling, noise, experiments,” Rev. Générale Therm., vol. 37, pp.
de procesamiento se dirijan a procesos de análisis automáticos 704–717, 1998, doi: 10.1016/S0035-3159(98)80048-2.
[13] D. Pan, Z. Jiang, W. Gui, K. Jiang, and X. Maldague, “Compensation
y cuantitativos sobre secuencias de imágenes. Method for the Influence of Dust in Optical Path on Infrared
Temperature Measurement,” IEEE Trans. Instrum. Meas., vol. 70,
La evaluación de desempeño para las técnicas de procesamiento pp. 1–11, 2021, doi: 10.1109/TIM.2020.3012015.
es un tema que ha presentado algunos avances en los últimos [14] D. Pan, Z. Jiang, X. Maldague, and W. Gui, “Research on the
Influence of Multiple Interference Factors on Infrared Temperature
años [32][80][56][36][17][53] utilizando enfoques como la Measurement,” IEEE Sens. J., pp. 1–1, Jan. 2021, doi:
relación señal-ruido, el criterio Tanimoto y la probabilidad de 10.1109/jsen.2021.3055757.
detección. Sin embargo, la creciente aparición de nuevas [15] J. Erazo-Aux, H. Loaiza-Correa, A. D. Restrepo-Giron, C. Ibarra-
técnicas de procesamiento exige desarrollos en indicadores de Castanedo, and X. Maldague, “Thermal imaging dataset from
composite material academic samples inspected by pulsed
desempeño que permitan inferir de manera objetiva, thermography,” Data Br., vol. 32, no. 106313, pp. 1–6, 2020, doi:
preferiblemente cuantitativa acerca de la eficiencia de los 10.1016/[Link].2020.106313.
métodos evaluados, sus ventajas y limitaciones. [16] D. L. Balageas, “Termografía Infrarroja : una técnica multifacética
para la Evaluación No Destructiva ( END ),” IV Conf. Panam. END,
no. January 2007, p. 14, 2007.
Países como Canadá, Francia, Rusia y Estados Unidos tienen [17] H. D. Benitez, “Contribución a la Caracterización de Defectos en
una trayectoria amplia y reconocida en prácticas relacionados Termografía Infrarroja mediante Máquinas de Aprendizaje,”
con termografía infrarroja con lo cual generan avances y Universidad del Valle, 2008.
publicaciones de forma continua. Los investigadores adscritos [18] A. D. Restrepo-Girón, “Contribucion a la Solucion del Problema de
Detección y Caracterizacion de Defectos en Termografia Activa
a universidades, grupos y/o centros de investigación con sede mediante Procesamiento Espacio-Temporal de Secuencias de
en estos países han sido los responsables de fomentar y realizar Imagenes Infrarrojas,” 2013.
transferencia de conocimiento a investigadores de otras [19] M. Z. Umar, V. Vladimir, H. Abdullah, and A. K. Ariffin, “Ultrasonic
naciones emergentes como Brasil, Colombia, México y Chile. infrared thermography in non- destructive testing : A review,” Russ.
J. Nondestruct. Test., vol. 52, no. April, pp. 212–219, 2016, doi:
10.1134/S1061830916040082.
AGRADECIMIENTOS [20] A. O. Chulkov and V. P. Vavilov, “Hardware and Software for
Thermal Nondestructive Testing of Metallic and Composite
Los autores agradecen al Departamento Administrativo de Materials,” J. Phys. Conf. Ser., vol. 671, no. January, p. 012011,
Ciencia, Tecnología e Innovación de Colombia (Colciencias); 2016, doi: 10.1088/1742-6596/671/1/012011.
Convocatoria 727 de 2015. Igualmente, a la Institución [21] V. P. Vavilov, A. O. Chulkov, D. A. Derusova, and Y. Pan, “Thermal
Universitaria Antonio José Camacho (UNIAJC). Proyectos de NDT research at Tomsk Polytechnic University,” Quant. Infrared
Thermogr. J., no. January, pp. 1–16, 2016, doi:
Investigación PD0916 y PD0121. 10.1080/17686733.2015.1131855.
[22] F. Khodayar, S. Sojasi, and X. Maldague, “Infrared thermography
REFERENCIAS and NDT: 2050 horizon,” Quant. Infrared Thermogr. J., vol. 6733,
no. July, pp. 1–22, 2016, doi: 10.1080/17686733.2016.1200265.
[1] X. Maldague, Theory and Practice of Infrared Technology for [23] V. P. Vavilov and D. D. Burleigh, “Review of pulsed thermal NDT:
Nondestructive Testing. New York: John Wiley-Interscience, 2001. Physical principles, theory and data processing,” NDT E Int., vol. 73,
[2] D. Balageas et al., “Thermal (IR) and Other NDT Techniques for pp. 28–52, 2015, doi: 10.1016/[Link].2015.03.003.
Improved Material Inspection,” J. Nondestruct. Eval., vol. 35, no. 1, [24] V. Vavilov, “Thermal NDT: historical milestones, state-of-the-art
pp. 1–17, 2016, doi: 10.1007/s10921-015-0331-7. and trends,” Quant. Infrared Thermogr. J., vol. 11, no. 1, pp. 66–83,
[3] C. Ibarra, “Quantitative subsurface defect evaluation by pulsed phase 2014, doi: 10.1080/17686733.2014.897016.
thermography: depth retrieval with the phase,” Laval University, [25] D. L. Balageas, “A brief history of QIRT concept, structures and
2005. community,” Quant. Infrared Thermogr. J., vol. 13, no. 1, pp. 109–
[4] V. P. Vavilov, “Thermal non destructive testing: short history and 125, 2016, doi: 10.1080/17686733.2015.1131852.
state-of-art,” in QIRT 92 - Eurotherm Series 27, 1992, pp. 179–194, [26] S. Marinetti et al., “Statistical analysis of IR thermographic
doi: [Link] sequences by PCA,” Infrared Phys. Technol., vol. 46, no. 1-2 SPEC.
Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira. 25

ISS., pp. 85–91, 2004, doi: 10.1016/[Link].2004.03.012. the Aero Space Industry : A Quantitative Approach,” in 11th
[27] H. Benitez, H. Loaiza, B. Bacca, and J. Acevedo, “Sistema de Vision European Conference on Non-Destructive Testing (ECNDT 2014),
Artificial para el Procesamiento y Analisis de Imagenes Infrarrojas,” 2014, no. Ecndt, pp. 1–10.
XI Symp. Image, Signal Process. Artif. Vision, STSIVA [Disponible [47] A. D. Restrepo-Giron and H. Loaiza-Correa, “New 3D Finite
en CD]. Difference Method for Thermal Contrast Enhancement in Slabs
[28] R. E. Ospina, S. D. Cardona, and B. Bacca-Cortes, “Software tool for Pulsed Thermography Inspection,” J. Nondestruct. Eval., pp. 62–73,
thermographic inspection using multimodal fusing of thermal and 2013, doi: 10.1007/s10921-013-0203-y.
visible images,” Ing. y Compet., vol. 19, no. 1, pp. 53–68, 2017, [48] E. Grinzato, Bison P.G., Marinetti S., and Vavilov V., “Thermal NDE
[Online]. Available: [Link] enhanced by 3D numerical modeling applied to works of art,” 15th
[Link]. World Conference on Non-destructive Testing, Rome, Italy, 2000.
[29] S. M. Shepard, “Advances in pulsed thermography,” in Thermosense [Link]
XXIII, 2001, pp. 511–516, doi: 10.1117/12.421032. (accessed Mar. 07, 2017).
[30] S. M. Shepard, “Flash Thermography of Aerospace Composites 2 . [49] A. D. Restrepo Girón, “Enhanced method for flaws depth estimation
Thermographic Signal Reconstruction,” in IV Conferencia in CFRP slabs from FDTC thermal contrast sequences,” Ing. e
Panamericana de END Buenos, 2007, vol. IV, no. Octubre, pp. 2–7, Investig., vol. 35, no. 3, pp. 61–68, 2015, doi:
[Online]. Available: 10.15446/[Link].v35n3.50552.
[Link] [50] M. A. Akhloufi and B. Verney, “Multimodal registration and fusion
[31] M. Pilla, M. Klein, X. Maldague, and A. Salerno, “New absolute for 3D thermal imaging,” Math. Probl. Eng., 2015, doi:
contrast for pulsed Thermography,” in 6th Conference on 10.1155/2015/450101.
Quantitative Infrared Thermography, 2002, pp. 53–58. [51] S. S. Pawar and V. P. Vavilov, “Applying the heat conduction-based
[32] H. Benítez, X. Maldague, C. Ibarra-Castanedo, H. Loaiza, A. 3D normalization and thermal tomography to pulsed infrared
Bendada, and E. Caicedo, “Modified differential absolute contrast thermography for defect characterization in composite materials,”
using thermal quadrupoles for the nondestructive testing of finite Int. J. Heat Mass Transf., vol. 94, pp. 56–65, 2016, doi:
thickness specimens by infrared thermography,” Can. Conf. Electr. 10.1016/[Link].2015.11.018.
Comput. Eng., no. May, pp. 1039–1042, 2007, doi: [52] Q. Zhang and X. Maldague, “An adaptive fusion approach for
10.1109/CCECE.2006.277741. infrared and visible images based on NSCT and compressed sensing,”
[33] F. J. Madruga, C. Ibarra-Castanedo, O. M. Conde, J. M. López- Infrared Phys. Technol., vol. 74, pp. 11–20, 2016, doi:
Higuera, and X. Maldague, “Infrared thermography processing based 10.1016/[Link].2015.11.003.
on higher-order statistics,” NDT E Int., vol. 43, no. 8, pp. 661–666, [53] C. San Martin, S. N. Torres, and J. E. Pezoa, “An effective reference-
2010, doi: 10.1016/[Link].2010.07.002. free performance metric for non-uniformity correction algorithms in
[34] A. D. Restrepo Girón and H. Loaiza Correa, “A new algorithm for infrared imaging system,” Conf. Proc. - Lasers Electro-Optics Soc.
detecting and correcting bad pixels in infrared images,” Ing. E Annu. Meet., pp. 576–577, 2007, doi: 10.1109/LEOS.2007.4382537.
Investig., vol. 30, no. 2, pp. 197–207, 2010. [54] V. Vavilov, “Noise-limited thermal/infrared nondestructive testing,”
[35] S. Marinetti and P. G. Cesaratto, “Emissivity estimation for accurate NDT E Int., vol. 61, pp. 16–23, 2014, doi:
quantitative thermography,” NDT E Int., vol. 51, pp. 127–134, 2012, 10.1016/[Link].2013.09.002.
doi: 10.1016/[Link].2012.06.001. [55] H. D. Benitez Restrepo, “Objective Image Quality Assessment in
[36] T. J. Ramirez-Rozo, J. C. Garcia-Alvarez, and C. G. Castellanos- infrared non-destructive testing,” in 2012 XVII Symposium of Image,
Dominguez, “Infrared thermal image segmentation using Signal Processing, and Artificial Vision (STSIVA), Sep. 2012, pp.
expectation-maximization-based clustering,” STSIVA 2012 - 17th 122–127, doi: 10.1109/STSIVA.2012.6340568.
Symp. Image, Signal Process. Artif. Vis., pp. 223–226, 2012, doi: [56] J. C. Garcia Alvarez, H. D. Benitez Restrepo, G. Castellano
10.1109/STSIVA.2012.6340586. Dominguez, and T. Ramirez, “Non referenced Quality Assessment of
[37] C. Ibarra-Castanedo, J. R. Tarpani, and X. P. V Maldague, Image Processing Methods in Infrared Non-destructive Testing
“Nondestructive testing with thermography,” Eur. J. Phys., vol. 34, ICIAP 2013,” Lect. Notes Comput. Sci. (including Subser. Lect. Notes
no. November, pp. S91–S109, 2013, doi: 10.1088/0143- Artif. Intell. Lect. Notes Bioinformatics), vol. 8156, no. PART 1, pp.
0807/34/6/S91. 722–732, 2013, doi: 10.1007/978-3-642-41181-6.
[38] D. Balageas, “In search of early time ... - Article 2013,” Adv. Opt. [57] F. Lopez, X. Maldague, and C. Ibarra-Castanedo, “Enhanced image
Technol., vol. 2013, no. January 2013, pp. 1–13, 2013, doi: processing for infrared non-destructive testing,” Opto-Electronics
[Link] Rev., vol. 22, no. 4, Jan. 2014, doi: 10.2478/s11772-014-0202-2.
[39] A. D. Restrepo-Girón and H. Loaiza-Correa, “Background Thermal [58] D. Balageas, J. Roche, and F. Leroy, “Quantitative Assessment of
Compensation by Filtering for Contrast Enhancement in Active Defect Characterisation Using the TSR Coefficient Images,” in
Thermography,” J. Nondestruct. Eval., vol. 35, no. 1, pp. 1–11, 2016, QIRT-Asia, 2015, no. October, pp. 1–11.
doi: 10.1007/s10921-016-0336-x. [59] D. L. Balageas, J.-M. Roche, and F.-H. Leroy, “Comparison and
[40] D. Balageas and J. Roche, “Common tools for quantitative time ranking procedure for an objective assessment of thermographic NDE
resolved pulse and step heating thermography part 1 theoretical methods,” Proc. 2016 Int. Conf. Quant. InfraRed Thermogr., no. July,
basics,” 2014, no. July, doi: 10.13140/2.1.2645.4402. 2016, doi: 10.21611/qirt.2016.003.
[41] J. M. Roche and D. Balageas, “2014-Common tools for quantitative [60] B. Liu, H. Zhang, H. Fernandes, and X. Maldague, “Quantitative
time resolved pulse and step heating thermography part II: Evaluation of Pulsed Thermography, Lock-in Thermography and
experimental validation,” QIRT 2014 Conf., no. July, 2014, doi: Vibrothermography on Foreign Object Defect (FOD) in CFRP,”
10.13140/2.1.1334.7207. Sensors, vol. 16, no. 5, p. 743, 2016, doi: 10.3390/s16050743.
[42] J. Roche, F. Leroy, and D. L. Balageas, “Information condensation in [61] V. Vavilov, S. Marinetti, and D. Nesteruk, “Accuracy issues in
defect detection using TSR coefficients images,” QIRT2014 modeling thermal NDT problems,” in Proceedings Volume 6939,
Conférence, no. 1, 2014, doi: 10.13140/2.1.1596.8648. Thermosense XXX;, 2008, no. March, p. 693913, doi:
[43] J.-M. Roche and D. L. Balageas, “Detection and characterization of 10.1117/12.775684.
composite real-life damage by the TSR-polynomial coefficients [62] F. Lopez, V. De Paulo Nicolau, C. Ibarra-Castanedo, and X.
RGB-projection technique,” in Proceedings of the 2014 International Maldague, “Thermal-numerical model and computational simulation
Conference on Quantitative InfraRed Thermography, 2014, vol. 0, of pulsed thermography inspection of carbon fiber-reinforced
no. July, pp. 1–10, doi: 10.21611/qirt.2014.010. composites,” Int. J. Therm. Sci., vol. 86, pp. 325–340, 2014, doi:
[44] X. Maldague, F. Galmiche, and A. Ziadi, “Advances in pulsed phase 10.1016/[Link].2014.07.015.
thermography,” Infrared Phys. Technol., vol. 43, no. 3–5, pp. 175– [63] V. Vavilov, W. Ś, and D. Derusova, “Ultrasonic and optical
181, 2002, doi: 10.1016/S1350-4495(02)00138-X. stimulation in IR thermographic NDT of impact damage in carbon
[45] N. Rajic, “Principal component thermography for flaw contrast composites,” Quant. infrared Thermogr., pp. 2–6, 2014, doi:
enhancement and flaw depth characterisation in composite 10.1080/17686733.2015.1046678.
structures,” Compos. Struct., vol. 58, no. 4, pp. 521–528, Dec. 2002, [64] V. P. Vavilov, “Modeling thermal NDT problems,” Int. J. Heat Mass
doi: 10.1016/S0263-8223(02)00161-7. Transf., vol. 72, pp. 75–86, 2014, doi:
[46] V. P. Vavilov and A. O. Chulkov, “Thermal NDT of Composites in 10.1016/[Link].2013.12.084.
26 Scientia et Technica Año XXVII, Vol. 27, No. 01, enero-marzo de 2022. Universidad Tecnológica de Pereira

[65] S. Sojasi et al., “Infrared Testing of CFRP Components : Jorge H. Erazo Aux. Ingeniero
Comparisons of Approaches using the Tanimoto Criterion,” in NDT
Electrónico (2004), Magister en Ingeniería
in Canada 2015 Conference, 2015, pp. 1–8.
[66] G. Jinlong, L. Junyan, W. Fei, and W. Yang, “Inverse heat transfer con énfasis en Ingeniería Eléctrica y
approach for nondestructive estimation the size and depth of Electrónica (2010) y Doctor en Ingeniería
subsurface defects of CFRP composite using lock-in thermography,” con énfasis en Ingeniería eléctrica y
Infrared Phys. Technol., vol. 71, pp. 439–447, 2015, doi:
Electrónica (2021) de la Universidad del
10.1016/[Link].2015.06.005.
[67] D. L. Balageas, “Thermal Imaging Methods,” Encycl. Struct. Heal. Valle, Cali – Colombia. Profesor
Monit., no. August, 2008, doi: 10.1002/9780470061626.shm018. Asociado de la Facultad de Ingeniería de
[68] A. Pereira F., J. Diniz M., C. Paskocimas, S. Nóbrega, F. de Araujo la Institución Universitaria Antonio
M., and A. Ricalde R., “Addition of magnetic markers for non-
Camacho, Cali – Colombia. Termógrafo Nivel I (2007) y Nivel
destructive evaluation of polymer composites,” Mater. Res., vol. 14,
no. 4, pp. 508–513, Dec. 2011, doi: 10.1590/S1516- II (2011) del Infrared Training Center – ITC. Áreas de interés:
14392011005000071. Termografía, visión artificial y procesamiento digital de
[69] B. Milovanović and I. B. Pečur, “Review of Active IR Thermography señales.
for Detection and Characterization of Defects in Reinforced
ORCID: [Link]
Concrete,” J. Imaging, vol. 2, no. 2, p. 11, 2016, doi:
10.3390/jimaging2020011.
[70] O. Motato and H. Loaiza, “Identificación biométrica utilizando Humberto Loaiza Correa. Ingeniero
imágenes infrarrojas de la red vascular de la cara dorsal de la mano,” Electricista (1990) y Magister en
Ing. e Investig., vol. 29, no. 1, pp. 90–100, 2009, Accessed: Mar. 10,
Automática (1995) de la Universidad del
2017. [Online]. Available:
[Link] Valle. Doctor en Robótica y Visión
[71] K. Estupinan Roldan, M. A. Ortega Piedrahita, and H. D. Benitez, Artificial (1999) de L’Université d’Evry,
“Design and implementation of a protocol for acquisition and Francia. Profesor Titular de la Escuela de
processing of infrared images obtained from hands,” Symp. Signals,
Ingeniería Eléctrica y Electrónica de la
Images Artif. Vis. - 2013, STSIVA 2013, pp. 1–7, 2013, doi:
10.1109/STSIVA.2013.6644944. Universidad del Valle. Codirector del
[72] W. M. Liu, K. Chang, S. Yoon, and A. M. Gorbach, “Reconstruction Grupo de Investigación en Percepción y
of Thermal Signals in Infrared Images Reveals Temperature Sistemas Inteligentes (PSI). Áreas de interés: robótica móvil,
Perturbations during Full Forearm Occlusion,” QIRT 2014 Conf.,
visión artificial y el procesamiento digital de señales.
2014, doi: 10.21611/qirt.2014.135.
[73] R. Gonzalez-Leal, M. Kurban, F. J. Gonzalez, and O. Cruz, ORCID: [Link]
“Quantitative human interpretation for breast thermography,” in 15th
Quantitative InfraRed Thermography Conference, 2020, pp. 1–5, doi:
10.21611/qirt.2020.078.
Andrés David Restrepo Girón. Ingeniero
[74] S. Guerrero, H. Loaiza, and A. D. Restrepo, “Segmentación
automática de imágenes térmicas de la mama como apoyo a la Electrónico (1999), Magister en
detección de cáncer,” Sist. Telemática, vol. 12, no. 30, pp. 25–34, Automática (2005) y Doctor en Ingeniería
2014, Accessed: Mar. 10, 2017. [Online]. Available: (2014) de la Universidad del Valle.
[Link]
Profesor Asociado de la Escuela de
[75] C. D. Ferrin and H. Loaiza, “A Real-Time Multispectral Computer
Vision System for Morpho-Thermal Analysis of Footprint Plantar,” Ingeniería Eléctrica y Electrónica de la
IEEE Lat. Am. Trans., vol. 13, no. 8, pp. 2680–2686, 2015, Accessed: Universidad del Valle. Áreas de interés:
Mar. 10, 2017. [Online]. Available: procesamiento digital de señales e
[Link]
imágenes, ultrasonido e instrumentación
015/13TLA8_32Ferrin.pdf.
[76] J. Roche and D. Balageas, “Advanced thermal imaging techniques for inteligente.
quantitative damage mapping of composite materials,” in 6th ORCID: [Link]
International Symposium on NDT in Aerospace, 2014, no. November,
pp. 12–17.
[77] C. Rongier, R. Gilblas, Y. Le Maoult, L. Redjem-Saad, and F.
Schmidt, “Infrared thermography applied to the validation of thermal
simulation of highluminance LED used in automotive front lighting,”
in 15th Quantitative InfraRed Thermography Conference, 2020, pp.
1–2, doi: 10.21611/qirt.2020.002.
[78] Y. Kawai, M. Ishikawa, H. Nishino, H. Ishigaki, and K. Ogawa,
“Estimation of remaining water in steel pipes using active
thermography,” in 15th Quantitative InfraRed Thermography
Conference, 2020, pp. 1–6, doi: 10.21611/qirt.2020.094.
[79] M. Kopeć and B. Więcek, “AC temperature estimation of power
electronic devices using 1D thermal modeling and IR thermography
measurements,” in 15th Quantitative InfraRed Thermography
Conference, 2020, pp. 1–7, doi: 10.21611/qirt.2020.161.
[80] R. Usamentiaga, C. Ibarra-Castanedo, and X. Maldague, “More than
Fifty Shades of Grey: Quantitative Characterization of Defects and
Interpretation Using SNR and CNR,” J. Nondestruct. Eval., vol. 37,
no. 2, pp. 1–17, 2018, doi: 10.1007/s10921-018-0479-z.