Control de las termografías

Problemática inicial:

Cuando empezamos a hacer termografías, tuvimos los siguientes problemas:

• Algunas cámaras venían con programas que no permitían exportar termografías modificadas con falso color.

Ejemplo de programa de tratamiento termográfico sin opciones de exportación

• La termografías realizadas tenían un pobre contraste inicial

• Si realizaban el paso a falso color no se podían controlar los límites de temperatura.
• Por supuesto, no podría existir correlación entre los datos Lidar del láser escáner terrestre con las termografías.

Mejora del contraste:

Tras mucho buscar encontramos la forma más sencilla que consiste, o bien en utilizar la biblioteca CV de código abierto o las utilidades Image Imagick.

La orden es muy simple: convertir original -normalize destino. El resultado es el siguiente:

Aplicando patrones termográficos:

Ya al menos se puede ver. Ahora toco el hacerla interactuar con un patrón. Por ejemplo: azul, verde, rojo.

Ejemplo con azul-verde-rojo: dará tres grupos de temperaturas.

Ejemplo con violeta-azul-cyan-verde-amarillo-naranja-rojo: 7 grupos de temperaturas

Dado que estamos en Galicia podemos ajustar el rango de lectura de 0 a 35ºC. de esta manera que ayudaría a ver tres temperaturas diferentes

• .... 5ºC: violeta
• 5-10ºC: azul
• 15-20ºC: cyan
• 20-25ºC: verde
• 25-30ºC: amarillo
• 30-35ºC: naranja
• 35 ... ªC: rojo

De este modo, la imagen tratada no sólo será más vistosa:

Sino que además sabemos la temperatura de cada elemento por el color.

Esto además nos deja abiertos a cualquier otro tipo de patrón. Imaginemos un patrón un poco más amplio: negro, marrón, rojo, naranja, amarillo, verde, azul claro, azul oscuro, violeta, gris, blanco:

Ecd Vl resultado sería:

Pudiendo llegar a crear efectos que rayan la psicodelia:

Resumiendo, podemos tener una amplia tipología de espectros relacionados directamente con la temperatura real obtenida de una termografía:

Correlacionado datos lidar y termografías:

Partimos de datos lidar y sus proyecciones (se pueden obtener desde la aplicación web: http://sitegi.enmacosa.com)

Tenemos una nube de puntos proyectada ortogonalmente:

Y/o la imagen en color (aunque éste se note poco)

También tenemos las imágenes termográficas:

Si cruzamos los datos de falso color y posición real de los puntos Lidar:
Aquí mostramos tres opciones posibles:

Nos decantamos por la tercera opción y ésta fue la aplicabilidad que le unimos a nuestra aplicación web:

A continuación comentaremos como realizar estas nubes de puntos termográficas desde la aplicación de escritorio desarrollada por Extraco, Misturas, Logica y Enmacosa.

Utilizando nuestro programa MEEL

Las particularidades que debe tener son:

• Creación de cualquier espectro:

Si deseáramos cambiar por ejemplo a <rosa, violeta, rojo, naranja, amarillo verde cian, azul gris negro> (pink blueviolet red orange yellow green1 cyan blue gray black) el espectro sería:

• Lectura de un grupo de archivos: Fácil, con un obketo "file".
• Geolocalización de termografías:

Para la correlación de datos tenemos los siguientes puntos de partida:

1. Podemos tener un GPS y con sus datos guardar la posición y el tiempo en que se realiza la lectura.
2. Al guardar las fotografías ya tenemos el tiempos de su creación
3. Como ya utilizamos con el programa los datos del eje

Conclusión de cada punto del eje (pk) podemos tener las coordenadas de las coordenadas el momento más cercano de los datos tomados del GPS y del momento del GPS buscaremos la foto más cercana creada en el tiempo. Una "simple" triple relación.

Y por último:

• Previsualización de la foto normalizada y la relativa al espectro.

Veamos más de cerca estos ejemplos:

Claro está, podemos así variar más fácilmente el espectro:

Con el espectro convencional:

Como conclusión final, podremos tomar cualquier termografía, georreferencia o mezclar sus datos con una nube de puntos con la simple ayuda de un GPS.

Photogrammetry and Geotechnics with VisualSFM / Fotogrametría y Geotecnia con VisualSFM

Manual fotogrametría práctica con VisualSFM para su uso en geología:

Descargable desde http://ccwu.me/vsfm. Se instala y se abre el programa: (antes hemos hecho unas cuantas fotos de un talud para tener un ejemplo)

Manual practice VisualSFM photogrammetry for use in geology: 

Downloadable from http://ccwu.me/vsfm. It installs and the program opens (before we made a few pictures of a slope for an example)

Se seleccionan las fotos y se le deja trabajar:
We select the photos and let we work at the computer:

A continuación se pica en : después en    y se le deja trabajar:
Then, we clik in:                          then wue click in        and let we work at the computer:

Ahora hay que poner los puntos en sus coordenadas correctas. IMPORTANTÍSISMO, EL PROGRAMA ES MUY BUENO PERO NO HACE MAGIA NO ES ADIVINO. Para ello se pica en

Para seleccionar un punto se busca ese punto y se pica con control+ratón(derecho):

Now we have to put points in their correct coordinates. VERY IMPORTANT, THE PROGRAM IS VERY GOOD BUT IT DOES NOT MAGIC, IT IS NOT SEER. For it, click in
To select a point and that point is looking chop with control + mouse (right):

En este caso les pondremos valores sencillitos a tres puntos. Posteriormente se pica de nuevo en  para decirle al programa que reacalcule de nuevo todos los puntos. Se pica con el botón de mayúsculas pulsado. (si no se pica con Mayúsculas pulsado no hace nada y hay que repetir el proceso.

Ahora lo que queda es guardar lo cambiado. No lo anterior.

Si se pulsa: SFM > SaveNViewMatch:

In this case we will put simple values ​​to three points. Then you click again to tell the program that reacalcule again all points. you click the button with shift pressed. (if you click with Shift pressed does nothing and have to repeat the process. 

Now need to save what changed. Not before. 

Pressing: SFM> SaveNViewMatch:

Lo guardamos como cosa.ply (EXTENSIÓN PLY)
Save it with the extension ply.

Otra forma de guardar en con CMVS: 
Another way is to save it with CMVS:

Y guardará todo como lo estamos dejando:
 Guardará el proyecto como archivo cantos.nvm y directorio cantos.nvm.cmvs:
All will be saved as file .nvm and .nvm.cmvs

Si se desea se puede abrir con MeshLab:
You can open it with MeshLab:

Queda chulo:
That's OK:

Ya en nuestra aplicación se le cambia de .ply a .las:
With our application, we can change format .ply to .las.

Y ya se puede leer:
Then we can read it:

Y el resto es ya trabajar con nuestro programa.
Now, we can work with our program.

Geotechnics, fractals and laser scanner. Can these three concepts have something in common? / Geotécnia, fractales y láser escáner. ¿Pueden estos tres conceptos tener algo en común?

Hace unos días, cuando empezamos a investigar si sería posible la complementación de conceptos como el láser escáner y la geotecnia, descubrimos que sí sería posible tal concomitancia. Al profundizar en cual sería el sistema de cálculo más adecuado, más cercano al sistema utilizado en campo, fuimos elaborando diversos procedimientos de cálculo con la ayuda de nuestro programa. 

A few days ago, when we began to investigate whether supplementation of concepts such as laser scanners and geotechnics be possible, such conjunction discovered that yes it would be possible. By delving into the system which would be more appropriate computer closer to the system used in the field, we were developing various computational procedures with the help of our program.

Al cartearnos con algunos de nuestros habituales seguidores de nuestro blog, en particular con R.R., surgieron varias dudas en torno a la forma de realizar estas mediciones. La duda era debida a que el sistema habitual era el de trepar físicamente por el talud y elegir planos característicos para la definición de los parámetros de buzamiento y rumbo. Estos planos característicos tienen estructura fractal y esta estructura es la que decidimos intentar reproducir con nuestro programa.

Some of our regular followers of our blog, in particular RR, several questions about how to perform these measurements arose. The doubt was due to the usual system was to physically climb the slope and choose characteristic to define the parameters of dip and direction planes. These planes are characteristic fractal structure and this structure is that we decided to try to play with our program.

El proceso definido anteriormente es bastante difícil de entender sin ejemplos claros y una detallada explicación. A continuación expondremos los avances que fuimos alcanzando en orden cronológico e intentaremos explicarlos con la mayor claridad posible.

The process defined above is quite difficult to understand without clear examples and a detailed explanation. Then we will present the progress that we were reaching in chronological order and we will try to explain as clearly as possible.

El ejemplo:
Hemos escaneado el siguiente talud:
We've scanned the following slope:



Y, por tanto, obtenido la nube de puntos:
And, therefore, obtained the point cloud:

Triangulamos en malla con nuestro programa:
Triangulated mesh with our program:

Una vez triangulado lo podemos ver en 3D:
Once triangulated, we can see it in 3D:

El cálculo de planos sería:
The calculation of planes would be this:

Y con la opción de sombras:
And with the option of shadows:

Si en lugar de dar estos pasos acometemos la dimensión fractal:
If instead of taking these steps we tackle the fractal dimension:

La selección fractal es un subconjuto (mayor cuanto mayor es el número de puntos) de aquellos planos claros que hubiera elegido un geólogo para la valoración de los diversos planos de definición.
The fractal is a subset selection (higher the greater the number of points) of those clear shots that had chosen a geologist for assessing various levels of definition.

Y de este subconjunto se hace la estadística de las familias de planos. Por polos:
And with this subset, the statistics of the families of planes are made. For poles:

Y por sombras:
And for shadows:

Si hubíeramos ido a tomar los datos de campo habríamos tomados los siguiente rumbos y buzamientos:
If we had gone to take the field data we had taken the following courses and dips:

155 48
122 35
126 48
282 66
325 88

Con los que hubiéramos obtenidos los siguientes diagramas:
With that we had obtained the following diagrams:

Se puede concluir que es más fiable el método fractal que la triangulación completa. Parece lógico. no todos los planos que aparecen en la triangulación tienen que ser planos de fractura característicos de las rocas. Un geólogo también preselecciona que pano tiene que medir antes de hacerlo.

La otra opción es seguir el procedimiento de la anterior entrada y decidir que puntos nos valen desde el 3D.


It can be concluded that it is more reliable to complete the triangulation fractal method. It seems logical. Not all levels shown in the triangulation must be flat characteristic fractured rock. A geologist also preselected to have to measure cloth before doing so.

The other option is to follow the procedure of the previous post and decide that we are worth points from the 3D.