Comparación de DTM usando Global Mapper 17.0.1

Hagamos hoy algo sencillo, comparar, primero cualitativamente (visualmente) y después cuantitativamente dos DTM. Por un lado elegimos una fuente muy usual, SRTM de 3 arc sec (aproximadamente 90m) con un DTM derivado de Fotogrametría Stereo.

  • Comparación CUALITATIVA (i.e visual)
  • Comparación CUANTITATIVA (i.e RMSE)

Abrimos por un lado un DTM cuya fuente sea SRTM, en este caso me he conectado via WMS (Web Mapping Service) a través del data online disponible dentro de la misma aplicación Global Mapper (File/Download Online Imagery/data). La resolución es de aproximadamente 90m (3 arc sec).

DTM-COMPARISON-20160212

Por otro lado he encontrado este DTM cuya fuente conozco (Stereo Photogrammetry). La resolución es de 5m.

DTM-COMPARISON-20160212-02

A través de la herramienta ‘digitizer tool’ (Tools/Digitizer) seleccionamos una línea dibujada al azar sobre los dos. Botón derecho del ratón-> analysis/measurement/path profile. Exporto ambas imágenes (es importante en path setup definir un mismo mínimo y máximo para poder compararlas adecuadamente).

Con Photoshop superpongo (Layer display/ multiply) ambas imágenes y veo cuán diferente son.

DTM-COMPARISON-20160212-03

Esto nos da una primera idea de la comparación, pero vayamos un poco más allá: ¿Cuál es el RMSE (Error medio cuadrático, Root Mean Square Error) entre ambas bases de datos?.

DTM-COMPARISON-20160212-04

Esta es una medida de desviación que nos va a definir mucho más exactamente que una simple visualización. Podéis ver algo más desarrollado este punto en este link de esta misma página:

https://geovisualization.net/2010/06/30/using-excel-to-calculate-the-rmse-for-lidar-vertical-ground-control-points/

DTM-COMPARISON-20160212-05

Ahora tan solo hemos de verificar que esta cifra sea la correcta teniendo en cuenta los valores de precisión prometidos en la entrega.

Espero que os haya resultado interesante, si así es, no olvidéis comentar, compartir o simplemente decir Hola. Cualquiera de estas opciones es apreciada.

Un saludo cordial,
Alberto C.
MSc GIS and Remote Sensing

Descargas del CNIG. Open Source bien hecho!

Hola amigos del GIS,
Por motivos de trabajo que no vienen al caso, he tenido que bucear de manera sistemática la web de descargas del CNIG. http://centrodedescargas.cnig.es/CentroDescargas/inicio.do
Una maravilla.

cnig-20160208-01

Por motivos que tampoco viene al caso, he de hacer esto mismo de vez en cuando en todos los Institutos cartográficos del mundo y el del CNIG es sin duda en el que me resulta más fácil, en el que el modelo de datos en más lógico y en el que los links son más fiables de todo el mundo. La única obligación es la atribución obligatoria de los datos. ¿No es mucho pedir, no? Desde el día 27 de diciembre, los datos del IGN son libres CC By 4.0.
https://creativecommons.org/licenses/by/4.0/

Por tanto es obligatorio que mencione la procedencia a pie de imagen, créditos, etc.., sobre todo en publicaciones, usos comerciales, artículos, etc… (Por ejemplo puede poner “<tal dato> CC by instituto Geográfico Nacional” o más bien “derivado de <tal dato” CC by ign.es” o similares…).

cnig-20160208-02

Ya sea porque necesitemos las imágenes del PNOA (Plan Nacional de Ortofotografía Aérea), un modelo digital del terreno de alta resolución o imágenes históricas de nuestro pueblo… tan solo hay que bucear un poco en el catálogo de geodatos del Instituto Geográfico Nacional (Centro Nacional de Información Geográfica) y los conseguiremos.

Por ejemplo, la semana pasada tuve que encontrar datos sobre algunas ciudades españolas para hacer varios escenarios 3D para un cliente y aquí encontré por un lado un DSM 5m elaborado con fuentes LIDAR, por otro lado me bajé de Cartociudad los datos relativos a vectores lineales, manzanas y luego desde la web de CATASTRO (https://www.sedecatastro.gob.es/OVCFrames.aspx?TIPO=TIT&a=masiv) me bajé las geometrías de todos los edificios de la ciudad (que planeo geoprocesar para eliminar las formas no deseadas y para adjudicar alturas precisas gracias al LIDAR bajado con anterioridad).

Por qué no añadir geometrías de Open Street Maps (https://www.openstreetmap.org/export) o de la propia Base Topográfica Nacional BTN25 para completar dicho escenario?

barcelona-bldg-osm-capture-20160112
MADRID-GISDATA

La verdad amigos es que desde que empezó a funcionar el Open Data, los Geógrafos y derivados tenemos mucho con lo que ‘jugar’ para hacer nuestros análisis.
http://idee.es/

Espero que os resulte interesante.

Un saludo cordial,

Alberto
Geógrafo/ Máster SIG UAH/ Diseñador Multimedia

Risk exposure. Geoprocessing using Open Source Data!! Next steps!!

Now that we have completed a first example, let’s continue with a real-world one. Its important working on a Data Model to define what we understand as a Risk and how important this is. Meaning. High voltage power lines are an actual risk but the closer we are, i guess the bigger the risk is, meaning i.e 3 if we are within 50m and 1 if we are 150m away… It’s only a guess.

risk3

Same thing related to antennas, Petrol stations, etc.

This is my Data Model defined over the city of Madrid, Spain.

1 LINES- Roads speed >50 km/h within 100m risk=3
2 LINES- Power lines within 100m risk=3

3 POINTS-

Antenna,
High voltage towers,
Petrol stations:

risk if within 50m=3; risk if within 100m=2; risk if within 150m=1; 

4 AREAS-

Cement factories,
Electric Sub-stations,
Waste storage facilities:

risk if within 50m=3; risk if within 100m=2; risk if within 150m=1; 

(NOTE: You can choose your own risk thresholds and importance. Also note these information downloaded from Open Source data (Cartociudad, CNIG) has not been double checked and it has been used as is).

risk2risk1

How is this risk, or these combination of risks impacting in the population of Madrid?

risk4

Can we extrapolate these patterns to other cities in the world?
We will definitely continue  this analysis shortly.

You can also visuallize this analysis using CartoDB, the field regarding “risk exposure level” is called ALL2, and ranges from 2 to 12:

Software: ArcGIS 10.3, Global Mapper 17, CartoDB

Please share if you enjoyed it… or just to say hello!

Alberto C
MSc GIS and remote sensing UAH

Jugando con CartoDB

Hace ya mucho tiempo que he oído hablar de CartoDB y que vengo practicando en su página web a visualizar bases de datos sencillas.

  1. Crea una cuenta
  2. Incorpora tus datos o tómalos de la galería
  3. Selecciona en modo datos la columna que quieres simbolizar/visualizar
  4. Conviértela en NUMBER si estuviera en STRING
  5. Selecciona en modo mapa en WIZARD
  6. COROPLETAS, columna _población
  7. Visualiza el resultado

cartodb03cartodb01cartodb02

Change detection – Detección de cambios en polígonos

change-detection-bogota-telemediciones-20151023-02
THE IDEA: DEMONSTRATING HOW DYNAMIC A CITY IS, THUS HOW IMPORTANT IS HAVING AN UPDATED DATASET
bogota-change-detection-20151105-02

THE FACTS: THE CITY OF BOGOTÁ IN COLOMBIA 2012-2014

Overall growth rate: -0.12% ONLY HAVING INTO ACCOUNT THE DIFFERENCE OF BUILDINGS CAPTURED BETWEEN 2012 AND 2014 (We can do this because we have used the same data capture model in both years)

(De acuerdo al censo catastral, para 2015 la ciudad incorporó 51.531 predios nuevos urbanos. En total, hay 2’402.581 predios en la ciudad, de esos, 266,9 millones de metros cuadrados son de área totalmente edificada. Source: http://www.eltiempo.com/bogota/crecimiento-bogota-/15394797)

bogota-change-detection-20151105

THE PROCEDURE: Centroids of buildings; Spatial join showing presence-absence, considering a 10m accuracy threshold, meaning if the centroid has not moved more than 10m, its the same building. If the centroid in 2012 is not in 2014, its considered as demolished. If a new centroid appears its considered new building.

DENSITY MAPS+3D buildings
Help to quickly focus on the highlights
bogota-change-detection-news-20151021

La geográfica cabecera de ‘Up in the air’ con George Clooney

Para alguien que maneja mapas a diario es una sorpresa interesante cruzarse con esta cabecera de la película de Jason Reirman protagonizada por George Clooney ‘Up in the air’. Nubes, campos de cultivo, ciudades en 2D, 3D. Qué bien he elegido mi profesión:-)

Para alguien que viaja mucho, salvando las distancias, este vídeo, también del mismo fil ‘Up in the air’ es un simpático recuerdo de lo que es un viaje pasando por uno y mil auropuertos, arcos voltaicos, empaquetando, desempaquetando…

Espero que os guste!

Alberto

(Fuente: Jose Ignacio Sánchez de Nosolosig)

Tipologías usuarios Madrid Río. Estadísticas y tendencias interesantes

Después de 6 semanas en Madrid, no ha habido un solo día en que no tuviera que meterme en Madrid Rio o bien para ir a trabajar, para llevar a la niña a la guardería, volver a casa o simplemente para pasear… Madrid Río se ha convertido en la espina dorsal de mis comunicaciones por la ciudad. Una obra con la que originalmente estaba en profundo desacuerdo (por su planificación y ejecución) se ha convertido en, de alguna manera en el eje que articula mis movimientos.

2015-08-21 08.38.20

Para la multitud de personas a las que como a mí, les ha sorprendido esta infrastructura en frente de sus casas hay posibilidad de encontrarse en el mismo metro cuadrado a alguien corriendo, alguien andando, una pareja con un carrito de niños, unos amigos en patines, alguna persona en bici… unos más rápido y unos más lento, todos hemos de convivir en un trazado de unos 7,5km de largo y no más unos metros de ancho, con lo que unas mínimas normas de circulación se imponen.

Partiendo de la máxima de que el peatón tiene prioridad, no se debería pensar que este pueda hacer de su capa un sayo y moverse a su antojo por el recorrido. Otra máxima debe ser que las bicicletas respeten un límite máximo de velocidad (así como la gente en patines, patinetes, segways, etc).

Mi punto de partida ha sido medir desde un mismo punto la pasada de los usuarios y tipificarlos de acuerdo a su sexo, edad aproximada, tipo de deporte que practican y si estaban ubicados de manera correcta en el recorrido de manera que pudieran interactuar de manera normal con los otros usuarios, minimizando al posibilidad de encontronazos, golpes, caídas, etc. Entendiendo como ‘correcto’ si los usuarios circulan por su derecha.

51243949b

He tipificado a 100 usuarios en dos momentos diferentes del día y en el mismo lugar, para poder estableceer comparaciones. Ahí van los datos y posteriormente los resultados del análisis y algunas preguntas abiertas para cuando haya más tiempo o más interés.

  • #1 Avenida de Manzanares 204/ Madrid RIO 20 de Agosto 2015 entre las 16:28 y las 17:45. 35º centígrados

estadisticas-01
estadisticas-01B

sexo: 0= varón, 1=mujer
tipo: 0=andando, 1=corriendo, 2=bicicleta, 3=otro (segway, patín, patinete, etc.)
resultado: 0=correcto, 1=incorrecto

edad mediana= 28 años
moda sexo=hombre
moda tipo deporte=bicicleta
densidad usuarios=78 usuarios/hora

porcentaje posición correcta: 72%
CORRECTO

correlación sexo-corrección?= 0.28, débil
correlación edad-corrección?= -0.29, hay correlación negativa (débil)
correlación tipo deporte-corrección?= -0.33, hay correlación negativa (débil)
(ver post sobre correlación de variables)

(…)
r=1, correlation is PERFECT
0.75<r<1, correlation is STRONG
0.5<r<0.75, correlation is MODERATE
0.25<r<0.5, correlation is WEAK
<0.25, almost NO correlation, both variables are hardy related
(…)

Resumiendo, a esta hora de la tarde, las 4 y pico del mes de agosto con unos 35 grados celsius, la densidad es de 78 personas a la hora, de las cuales el 72% circula de manera correcta.

El perfil tipo de usuario a esta hora es el de VARON, CICLISTA, 28 AÑOS, POSICIÓN EN LA VÍA CORRECTA

Hay una correlación débil entre sexo y posición correcta, lo que quiere decir que las mujeres y hombres se ubican de manera incorrecta sin seguir ningún patrón definido o lo que es lo mismo, entre los mal colocados, casi el mismo número eran mujeres que hombres.

tampoco hay una relación clara de correlación de acuerdo a la edad o el tipo con respecto a la corrección o no de la ubicación.

  • #2 Avenida de Manzanares 204/ Madrid RIO 21 de Agosto 2015 entre las 09:20 y las 9:42. 28º centígrados

estadisticas-02
estadisticas-02b

sexo: 0= varón, 1=mujer
tipo: 0=andando, 1=corriendo, 2=bicicleta, 3=otro (segway, patín, patinete, etc.)
resultado: 0=correcto, 1=incorrecto

edad mediana= 35 años
moda sexo=hombre
moda tipo deporte=bicicleta
densidad usuarios=273 usuarios/hora

porcentaje posición correcta: 90% 

CORRECTO

correlación sexo-corrección?= 0.08, débil
correlación edad-corrección?= -0.07, hay correlación negativa (débil)
correlación tipo deporte-corrección?= -0.22, hay correlación negativa (débil)
(ver post sobre correlación de variables)

(…)
r=1, correlation is PERFECT
0.75<r<1, correlation is STRONG
0.5<r<0.75, correlation is MODERATE
0.25<r<0.5, correlation is WEAK
<0.25, almost NO correlation, both variables are hardy related
(…)

Resumiendo, a esta hora de la mañana, las 9 y pico del mes de agosto con unos 28 grados celsius, la densidad es de 273 personas a la hora, de las cuales el 90% circula de manera correcta.

El perfil tipo de usuario a esta hora es el de VARON, CICLISTA, 35 AÑOS, POSICIÓN EN LA VÍA CORRECTA

Hay una correlación casi inexistente entre sexo y posición correcta, entre los mal colocados, hay casi el mismo número eran mujeres que hombres, tampoco hay una correlación de acuerdo a la edad o el tipo con respecto a la corrección o no de la ubicación.

Ahora dejo algunas pregutas en el aire, es siempre el perfil tipo el de varón en bici de mediana edad o por el contrario hay picos horarios o ubicaciones donde este perfil es diferente. Podríamos encontrar alguna correlación mayor entre la posición correcta en el recorrido y alguno de los tipos estudiados?. Hay algún otro tipo (por ejemplo nivel de estudios o algún rango específico de edad) en el que veamos una relación clara con la correcta/incorrecta ubicación?.

El estudio específico de estas correlaciones podría permitir informar adecuadamente a los usuarios a através de paneles informativos y de esta manera reducir los potenciales golpes entre las personas que disfrutan de Madrid Río pero también ayudaría a integrar correctamente a los diferentes grupos de usuarios para que todos disfrutemos más adecuadamente de estas instalaciones.

Espero que te haya parecido interesante, si necesitas o quieres más información o aclaración, no dudes en ponerte en contacto conmigo por email.
Un saludo cordial!.

Alberto C.
MSc GIS

Freelancing is good?. Well, kind of!

Hi guys… fed up with this never ending heat we have these days in Madrid (almost 4 weeks already with more than 40º celsius max!). Just wanted to share with you one of the reasons why its good freelancing… you can wear hawaian t-shirts without any problem!

freelancing-20150729

The bad side is working on your own most of the time can get you completely mad. Thank god I have a nice bar in front of my co-working office so i can forget about GIS for a short while!.

Cheers!