Adding the GPR and / or geology / Añadiendo el georradar y/o la geología

Geología + Láser escáner terrestre (Geo+LST)
Geology + Laser scanner terrestrial (Geo + LST)

Otro tipo de datos que se puede añadir a nuestro GIS de carreteras en láser-escáner es el perfil del georradar.
Another type of data that can be added to our GIS road in laser-scanner is the GPR profile.

Georradar + Láser escáner terrerstre (GPR+LST)
GPR + LST

¿Qué es y cómo se obtiene el perfil del georradar?
What is it and how to obtain the GPR profile?

Una buena explicación la tiene la wikipedia:
A good explanation is held by the wikipedia:

http://es.wikipedia.org/wiki/Georradar

O estos artículos:
Or these items:

http://webs.uvigo.es/grupotf1/research/IWAGPR2005.pdf
http://webs.uvigo.es/grupotf1/research/rial_et_al_bruges.pdf

En esencia es la toma de reflexiones de ondas electromagnéticas que dependiendo de su longitud de onda profundizarán más o menos en el terreno. Cuanta más profundidad más vaga será la información o más nítida si el alcance resultara menor.

El aparato de emisión y lectura de señal suelen ser así:

In summary is taking reflections of electromagnetic waves depending on its wavelength deepen more or less on the ground. More would be more ambiguous information or if the range sharper lower result.

The emission device and signal reading usually as follows:

Georradares de 1GHz, 800MHz y 500MHz
GPR of 1GHz, 800MHz and 500MHz

Y las señales que reciben son:
And the signals they receive are:

Si le añadimos un remolque no metálico homologado de georradar de la empresa Enmacosa
If we add a trailer approved nonmetallic company GPR Enmacosa

Fuente Enmacosa: Remolque homologado no metálico para  GPR

Entonces tendremos un perfil muy largo del que tendremos que extraer el trozo que nos interese:
Then we will have a very long profile we have to extract the piece that interests us:

En la figura superior se aprecia un perfil muy largo correspondiente a la lectura de una carretera con el remolque anteriormente citado. En cada pk determinado tendremos que extraer un trozo del perfil correspondiente a nuestro pk y su entorno.
In the figure above shows very long profile corresponding to the reading of a road with the aforementioned trailer. In each particular KP have to remove a piece of the profile corresponding to our station and its surroundings.

¿Cómo se extrae el pk que estamos buscando?
How do I extract the KP we are looking for?

Hay dos formas:

1ª Desde los archivos fuente
2ª Desde el destino final


There are two ways:

From 1st source files
2nd Since the final destination 


1ª Desde los archivos fuente
From 1st source files

Los archivos fuente suelen ser una pareja de archivos. Un archivo de texto .rad:
The source files are usually a couple of files. A text file .rad:

SAMPLES:1024
FREQUENCY:24081.298464
FREQUENCY STEPS:18
SIGNAL POSITION:-0.031867
RAW SIGNAL POSITION:46580
DISTANCE FLAG:0
TIME FLAG:1
PROGRAM FLAG:0
EXTERNAL FLAG:0
TIME INTERVAL:0.100000
DISTANCE INTERVAL:0.000000
OPERATOR:
CUSTOMER:
SITE:
ANTENNAS:1000 MHz shielded
ANTENNA ORIENTATION:NOT VALID FIELD
ANTENNA SEPARATION:0.100000
COMMENT:
TIMEWINDOW:42.522624
STACKS:16
STACK EXPONENT:4
STACKING TIME:0.163840
LAST TRACE:2803
STOP POSITION:280.300000
SYSTEM CALIBRATION:0.0000023070
START POSITION:0.0000000000
SHORT FLAG:1
INTERMEDIATE FLAG:0
LONG FLAG:0
PREPROCESSING:0
HIGH:0
LOW:0
FIXED INCREMENT:0.0000000000
FIXED MOVES UP:0
FIXED MOVES DOWN:1
FIXED POSITION:0.0000000000
WHEEL CALIBRATION:2490.9000000000
POSITIVE DIRECTION:1

Con la información concerniete al tipo de antena y a la forma de guardar datos. Y un archivo binario .rd3. para leer este primero hay que leer el anterior. El primer dato SAMPLES nos dice que cada lectura (columna) tiene 1024 valores de tipo entero (2 bytes). Si leyéramos estos valores como columna y los pusiéramos uno tras otro en fila y a cada valor le hiciéramos corresponder un pixel en gris (dividiendo el valor entero por 256) tendríamos la típica gráfica del georradar.

With the information concerning the type of antenna and how to save data. And a binary file. Rd3. to read this first you have to read the last. The first data, SAMPLES, tells us that each reading (column) has 1024 integer values ​​(2 bytes). If we read these values ​​as column and we put one after another in a row and each value did you map a gray pixel (integer value dividing by 256) would the typical plot of GPR.

La forma de seleccionar un trozo sería muy sencilla si conocemos el pk de inicio de la lectura y el pk final. una simple regla de 3 nos colocaría en el punto exacto del pk buscado. Por ejemplo, supongamos que el tramo dibujado fuera el tramo entre el pk 0+000 al pk 0+100 y tiene 1.000 columnas de datos. Si queremos tener desde el pk 0+030 al 0+040 tendríamos que leer el archivo desde la posición 1024x2x300 a la 1024x2x400:

How to select a piece would be very simple if you know the home KP and KP reading end. a simple rule of three would put us in the exact point of pk sought. For example, suppose that the section drawn out the stretch between KP 0 KP 0 +000 to +100 and has 1,000 columns of data. If we want from KP 0 +030 to 0 +040 would have to read the file from position 1024x2x400 1024x2x300 to:

Un caso particular consiste en tener además el archivo de coordenada de la captación de datos desde GPS. El archivo .cor tiene este formato:
A particular case is also have coordinate file uptake from GPS. The .cor has this format:

1,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

2,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

3,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

4,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

5,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

6,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

7,2013-5-31,13:36:23,78:01:00.50 NN,33:58:52.89 WW,2097.699219M,1.459554

8,2013-5-31,13:36:24,78:01:00.50 NN,33:58:52.88 WW,2097.699219M,1.459554

[...]

A partir de él se puede colocar cada columna de datos en un plano/mapa georreferenciado.
From it you can place each column of data in a plane / georeferenced map.

2ª Desde el destino final
2nd Since the final destination

Puede que nos interese aplicar filtros y obtener un resultado deducido de los programas de tratamiento de archivos de georradar. En ese caso hay que trabajar con el resultado que suele ser la imagen tratada. el sistema es el mismo sólo que en lugar de conseguir el extracto de un archivo hay que jugar con píxeles. La técnica es la misma pero no el medio. Es aconsejable usar ImageImagick o su mezcla son PHP.
It can be required to apply filters and get a result deduced from processing programs for GPR files. In that case you have to work with the result that the image is usually treated. the system is the same except that instead of getting the extract of a file you have to play with pixels. The technique is the same but not the medium. It is advisable to use ImageImagick or their mixture are PHP.

¿Qué es y cómo se obtiene el perfil geológico?

What is it and how to obtain the geological profile?

Sencillamente, con una máquina de sondeos, sus laborantes, un geólogo, un delineante (en algunos casos coinciden las personas) y sobre todo conocimiento, trabajo y paciencia.

A partir de:

Simply a polling machine, their laboring, a geologist, a draftsman (in some cases match people) and all knowledge, work and patience.

Starting from:

Se crea: (No olvidar utilizar experiencia y esfuerzo. Esta parte no la puede hacer una máquina.) It is created: (Do not forget to use experience and effort. This part can not be done by a machine.)

Resultados de incorporar el georradar:

GPR results incorporate:

Fuente: http://sitegi.enmacosa.com. Visualizando el georradar
From: http://sitegi.enmacosa.com. Visualizando el georradar

Resultados de incorporar la geología:

Results of incorporating the geology:

Fuente: http://sitegi.enmacosa.com, Añadiendo la geología
From: http://sitegi.enmacosa.com, "Añadiendo la geología"