lunes, 23 de noviembre de 2020

Estilos QGIS para representar medidas de líneas (o polígonos).

Vamos hoy con otro poco de QGIS. En este caso queremos mostraros cómo representar las medidas de una línea o polígono mediante simbología de estilos para tener una salida gráfica atractiva y sencilla. Como una imagen vale más que mil palabras, dentro portada:
Lo primero que necesitamos es añadir los estilos que vamos a utilizar posteriormente en líneas y polígonos, y que son archivos XML con el código necesario para que automáticamente apliquen una simbología predefinida a los vectores. Para las líneas vamos a utilizar el estilo Measure Line disponible en el repositorio de estilos de QGIS.
Descargamos el archivo XML y lo añadimos a QGIS desde el Administrador de estilos:
Hacemos lo propio con el estilo para los polígonos. En nuestro caso hemos utilizado el estilo de Faunalia llamado Distance Measurement Styles, que podemos encontrar desde el repositorio del plugin QGIS Resource Sharing.
Como casi todo en QGIS hay diferentes opciones para un resultado similar, así que podéis utilizar el estilo que más os guste de los disponibles (u obviamente elaborarlo vosotros mismos).

El siguiente paso es crear un par de nuevas capas shapefile, una de líneas y otra de polígonos, a fin de poder probar nuestros estilos. Lo hacemos en un sistema de coordenadas proyectada en metros como unidad de medida. Usamos también como fondo la ortofoto PNOA para visualizar aquello que queremos medir.
Dibujamos después un par de líneas de ejemplo que nos servirán para la prueba (click derecho termina línea):
Seguimos el mismo proceso para crear la capa de polígonos y tras dibujar el polígono que usaremos de prueba ya tenemos el lienzo terminado con ambos elementos en sus capas (un par de líneas y un polígono).

El último paso es aplicar los estilos que hemos añadido a QGIS al comienzo del artículo a su capa correspondiente. Vemos el ejemplo con la capa de líneas:
Os enseñaría el paso con las medidas del polígono pero no ha habido manera, y después de muchos intentos con muchas proyecciones diferentes siempre me lo etiquetaba con medidas de 0 metros o similares... No entiendo porqué y me quedo con la sensación de que me ha funcionado cuando le ha dado la gana (como en la primera imagen de este artículo... ¿?¿?). Me rindo. Si alguno sabéis la clave para que funcione siempre...

Una última utilidad. Si queremos tener un campo con las medida de cada una de nuestras líneas en la tabla de atributos simplemente usaremos el comando $length en la calculadora de campos (en el ejemplo con dos decimales).
Para la superficie y perímetro de los polígonos podéis usar los algoritmos $area y $perimeter respectivamente.
Finalmente reseñar que los estilos son totalmente modificables en su simbología (colores, anchuras, distancias, etc...) y agradecer a sus creadores el ponerlos al alcance de los que nos queda muy grande programarlos. Saludos.


jueves, 8 de octubre de 2020

Nuevo producto IGN: cartografía antigua.

Disponible a descarga el nuevo producto Cartografía antigua con los mapas que se conservan en la Cartoteca del IGN desde el siglo XV hasta el siglo XIX de España y distintas partes del Mundo. Puede realizarse la descarga de los datos a través de la agrupación Documentación geográfica y cartografía antiguas del Centro de Descargas del CNIG.

Se trata de ficheros digitales, imágenes JPG y documentos en PDF, resultado del escaneado de cartografía y atlas antiguos (entre los siglos XV y XIX) pertenecientes a los fondos cartográficos del Instituto Geográfico Nacional. Son mapas, cartas náuticas, vistas topográficas y atlas de múltiples autores, tanto españoles como extranjeros, la mayoría originales y en algunos casos reproducciones facsímiles. Hay mapamundis, mapas de España y de distintos lugares del mundo.

Son ficheros sin georreferenciar con una resolución de 250 ppp para los JPG y cuya unidad de descarga es el documento individual. A fecha de este artículo la colección cuenta con un total de nada más y nada menos que 1752 documentos (1675 en JPG y 77 PDF's). Pensamos que se trata de trasladar los fondos de la Cartoteca al Centro de Descargas, quizá para unificar aún más el material disponible dependiente del mismo organismo.

Ante semejante aluvión de documentos podemos al menos establecer algún pequeño filtro por siglo en la búsqueda. También por tipo de fichero (JPG o PDF) e incluso el habitual buscador genérico que puede delimitarnos la búsqueda inicial (por ejemplo: "África"). El localizador espacial con el área que ocupa el mapa también puede sernos especialmente útil para este producto.


PDF con Plano del Real Sitio de Aranjuez. 1810
En el Centro de Descargas simplemente tendremos la imagen o documento, sin ninguna información accesoria más, por lo que para todos los detalles de cada documento recomendamos seguir accediendo a su referencia en la Cartoteca (ejemplo anterior).
JPG con Cartas Naúticas de Cádiz. 1584
Desde luego para los amantes de la cartografía antigua se configura un acceso sencillo, potente y rápido a una buena colección de documentos que harán la delicia de cualquier aficionado. Somos de los que compartimos esa opinión de que los mapas antiguos tienen un alma diferente, algo que atrae por enigmático, remoto y trazos particularmente bellos. Agradecemos al IGN esta nueva iniciativa en la que seguramente "perderemos" unas cuantas horas en el futuro.


miércoles, 30 de septiembre de 2020

Asignar coordenadas y alturas a puntos en QGIS.

Hoy vamos a describir un par de operaciones básicas que podemos necesitar realizar sobre una capa de puntos vectorial. Esta capa es posible que se encuentre desnuda de atributos y podamos enriquecer sus datos con unos campos tan importantes como las coordenadas geográficas y la altitud de esos puntos.
Partimos de un ejemplo práctico en el que tenemos cargado en QGIS un archivo de relieve sobre cuya extensión generaremos un determinado número de puntos de forma aleatoria. Usamos para ello la herramienta Random points in extent (podemos buscarla en la cajita inferior izquierda del programa):
Rellenamos el Input con la extensión de nuestra capa MDT.
Establecemos en 30 el número de puntos a crear a lo largo de la extensión de la capa y ejecutamos el algoritmo.
Tenemos nuestra nueva capa temporal de puntos llamada Random points con el único atributo de su valor numérico correlativo ID.

El siguiente paso será darle los valores de altitud a cada punto, proveniente de nuestra capa de relieve MDT, utilizando el complemento Point Sampling Tool.
Con nuestra capa de puntos seleccionada en el complemento marcamos el MDT como fuente de datos para los valores y le damos al archivo una ubicación y un nombre (en el ejemplo puntos Z.shp). En la pestaña Fields podemos personalizar el nombre del campo atributo para el dato de altitud (en el ejemplo Z_MDT25).

Ejecutamos el complemento y la capa puntos Z.shp se añade a nuestro proyecto. El siguiente paso será darle unas coordenadas a cada punto para lo que utilizaremos el algoritmo Saga Add coordinates to points

Una nueva capa temporal Points with coordinates se añade a nuestro proyecto. En este caso ya conteniendo todos los atributos que buscábamos: la altitud y las coordenadas de cada punto.

Solamente nos restaría guardar nuestra capa en el formato que deseemos antes de cerrar el proyecto.

PD: en este caso hemos usado un random points para generar una capa de puntos de ejemplo sobre la que trabajar, pero obviamente podemos partir de una capa de puntos ya dada.
PD2: los algoritmos disponibles en QGIS son muchos y a veces dependen de los proveedores instalados, por lo que hay otras formas de insertar datos en los puntos. Ésta es solamente una.

miércoles, 9 de septiembre de 2020

Mosaicos nacionales históricos Landsat.

Disponible nuevo producto en el Centro de Descargas del Centro Nacional de Información Geográfica (CNIG). Bajo la categoría Ortofotos e Imágenes de Satélite se ofrece la serie de mosaicos nacionales históricos de ortoimágenes por satélite Landsat 5 y 8, a 30 y 15 metros de resolución de pixel respectivamente. La unidad de descarga es un mosaico para la Península, Islas Baleares, Ceuta y Melilla, y otro mosaico para las Islas Canarias. Desde 1986 se generan cada década dos modalidades: color natural y falso color infrarrojo.
El Sistema Geográfico de Referencia (SGR) es ETRS89 en la Península, Islas Baleares, Ceuta y Melilla, y WGS84 en las Islas Canarias. Proyección UTM en su huso correspondiente.

Las imágenes empleadas en los mosaicos son las procesadas por el United States Geological Survey (USGS) a nivel L1TP. Por lo tanto, tienen corrección geométrica pero no corrección atmosférica. Cada mosaico Landsat 5 se presenta en falso color infrarrojo (bandas 4-3-2), y una composición de color natural y pseudocolor (bandas 4, 3, 2 y 1). Cada mosaico Landsat 8 se presenta en falso color infrarrojo (bandas 8, 6, 5 y 4), y una composición de color natural y pseudocolor (bandas 8, 4, 3 y 2). Los datos están corregidos de nubosidad y la periodicidad de la serie es de un producto por década.
Listado completo de ficheros disponibles a fecha de este artículo. Quizá se echa de menos especificar si se trata de imagen Landsat 5 u 8, cosa que no sabemos hasta comenzar la descarga del archivo y ver su nomenclatura, aunque por el tamaño de los archivo tampoco es difícil deducir que las imágenes Landsat 8 son las de un peso mucho mayor al tener más resolución (2014).
Los archivos se ofrecen según el IGN en formato COG. Y si os ha pasado lo mismo que a mi que al leerlo me quedé igual de extrañado por no reconocerlo, se trata de una especie de acrónimo de Cloud Optimized GeoTIFF, un formato de imágenes para el procesamiento geoespacial nativo de la nube. Se trata de un archivo GeoTIFF normal, destinado a ser alojado en un servidor de archivos HTTP, con una organización interna que permite flujos de trabajo más eficientes en la nube. Presentan una nomenclatura como la de este ejemplo:
pnt_landsat8_2014_mosaico_islas_canarias_b432_hu28_8bits_COG_90.tif

La información auxiliar del producto nos proporciona la distribución de las teselas Landsat, tanto en formato PDF como SHP.
En el caso de la versión de las teselas en SHP las mismas ofrecen una buena cantidad de atributos que podrían sernos útiles a la hora de trabajar con la teselas territoriales de las huellas de imágenes Landsat.
Un vistazo por último a una de las imágenes de las Islas Canarias, en este caso Landsat 8 de 2014 a 15 metros de resolución y color natural. Debajo detalle a zoom 100% o pixel real sobre la zona de Santa Cruz de Tenerife.
Esto es todo por hoy. Saludos y hasta la próxima entrada.

lunes, 3 de agosto de 2020

Modelo Digital de Pendientes MDP05

Disponible en el Centro de Descargas CNIG el Modelo Digital de Pendientes con paso de malla de 5 m, generador a partir del MDT-LIDAR 1ª cobertura. Ficheros cortados por hojas del Mapa Topográfico Nacional a escala 1: 50.000. En las zonas de solape entre husos se ha calculado cada hoja en los dos husos. Cada celda de la malla almacena el valor de la pendiente en grados y fracción de grados sexagesimales. Realice la descarga a través de la agrupación Modelos Digitales de Elevaciones.

Según la información auxiliar del producto se han seguido los siguientes pasos para obtener el modelo digital de pendientes: 
  1. Generación de un catálogo con las hojas del modelo digital de elevaciones MDE05, con 5 metros de paso de malla, cortado por hojas 50 mil y obtenido a partir de la clasificación de los datos lídar 1ª cobertura. 
  2. Cálculo de la pendiente para cada celda. Para ello se calcula la pendiente entre cada celda y cada una de las celdas colindantes, tomándose de todas ellas la mayor en valor absoluto.
  3. El mosaico formado por todas las hojas del mapa de pendientes se recorta por hojas 50 mil para cada huso, manteniendo un solape entre husos de 1 hoja.
Representación del MDP05 con un degradado de color

Formato de archivo ASCII matriz ESRI (asc). Sistema geodésico de referencia ETRS89 (en Canarias REGCAN95, compatible con ETRS89) y proyección UTM en el huso correspondiente a cada hoja del MTN50. En Canarias el huso UTM es el 28. 
El MDP es un producto derivado del MDT el cual fue obtenido a partir de la clasificación de la nube de puntos LIDAR de la 1ª cobertura del proyecto PNOA-LIDAR. En el proceso de cálculo del MDP a cada celda del fichero asc se le asigna el valor de la mayor de las pendientes calculadas entre la celda homóloga del MDT, con iguales coordenadas, y cada una de las celdas colindantes. No se dispone de MDP de Ceuta, Melilla ni Gibraltar ya que no existen datos LIDAR en estas zonas.
Obtención del dato de pendiente en una coordenada concreta
Uso de la calculadora ráster para extraer las zonas con pendientes entre 15 y 30 grados


miércoles, 24 de junio de 2020

WMS de los Modelos Digitales de Andalucía 2016.

El artículo de hoy quiere desglosar el servicio web interoperable de la Junta de Andalucía que, a través del Instituto de Estadística y Cartografía de Andalucía, ofrece el Modelo Digital del Terreno (MDT) y los Modelos Derivados generados a partir de datos de dos vuelos fotogramétricos realizados en el año 2016. Uno de ellos de la mitad norte de Andalucía con resolución de 0.5 metros/pixel y otro para la zona Sur de Andalucía con resolución 0.25 metros/pixel.

Decir antes de nada que este MDT 2016 de Andalucía solamente parece estar disponible a través del servicio web, ya que en el portal cartográfico andaluz Line@ solamente se ofrece para descarga en formato TIFF el MDT 2008-2009 (10 metros/pixel). Imaginamos que está en fase de producción y que su resolución estimada dada la de sus vuelos origen será de 5 metros/pixel.
Se trata pues de un servicio WMS cuya URL es https://www.ideandalucia.es/wms/mdt_2016
Cargamos el WMS en QGIS y según observamos en la imagen superior dicho servicio consta de cuatro capas diferenciadas. La primera de ellas, denominada modelo_digital_terreno_2016_color, muestra el Modelo Digital del Terreno del año 2016 representado mediante un gradiente de tintas hipsométricas. Esta clásica representación de colores por alturas ofrece una paleta muy atractiva.

La segunda de las capas se denomina modelo_digital_terreno_2016_escala_grises y muestra el Modelo Digital del Terreno del año 2016 representado mediante un gradiente de escala de grises.

La tercera capa es la llamada modelo_digital_pendientes_2016 y muestra el Modelo Digital de Pendientes obtenido a partir del MDT.

Y la última de las capas, denominada sombreado_orografico_2016, nos muestra el sombreado orográfico de Andalucía obtenido a partir del modelo digital del terreno del año 2016.

Es ésta la única de las capas totalmente opaca, ya que las tres primeras vienen predefinidas con cierto grado de transparencia, ideal para ser superpuestas sobre la del sombreado y conseguir un efecto 3D mucho más realista en cada caso. En ninguna de las capas hemos detectado función GetFeatureInfo por la que podamos obtener la combinación RGB, el dato de altitud o alguna otra información.

En el visor geográfico de la IDE Andalucía está presente, como capa de fondo, una combinación de las capas 1 y 4 de este servicio (MDT color+sombreado), sin duda la combinación más clásica y aparente para representaciones cartográficas, por lo que es también la que nosotros hemos escogido para agregar a nuestro lote de mapas para SASPlanet. Ya es cuestión de los gustos de cada usuario utilizarla como capa, como mapa o superponer ésta o aquella información para conseguir productos más atractivos o útiles.

viernes, 22 de mayo de 2020

De QGIS al mundo Google.

Hoy traemos una rápida entrada con dos complementos muy sencillos que conectan nuestro lienzo de QGIS con el  mundo cartográfico Google, lo que siempre puede resultar útil e interesante.

El primero de ellos es Send2GE (Send to Google Earth), simple herramienta que graba la posición del mapa en la que pinchemos y la envía a Google Earth, creando una marca de posición. Si tenemos Google Earth cerrado el complemento lo arranca automáticamente. Y si marcamos varias posiciones éstas se van acumulando en los lugares temporales de Earth.

El segundo complemento es Street View, que abre la legendaria vista de pie de calle de Google en nuestro navegador, justo en el punto donde hayamos pinchado en QGIS. Si al pinchar estiramos con el ratón, una línea roja nos permite escoger la dirección de la vista.

Ambos complementos generan al instalarse tanto una entrada en el menú como un icono en la barra de botones.

Por cierto que probando en lugares donde no había calle física (grupo de edificios en este ejemplo) me he llevado alguna sorpresa y Street View me ha abierto vistas como esta....
Pero obviamente si pinchamos en lugares alejados de toda cobertura de Street View el navegador nos devolverá una imagen en negro. Esto se echa un poco en falta, algún indicador que nos asegure o no si Street View cubre la zona en la que vamos a pinchar.

En resumen dos complementos muy sencillos y sin ninguna pretensión ni complicación que unen nuestro escritorio de QGIS con el mundillo Google a través de Earth o Street View, y que muchas veces pueden abrirnos las puertas a trabajar con diferentes puntos de vista desde una misma posición. Saludos.

miércoles, 13 de mayo de 2020

Mapa LIDAR

El Centro de Descargas del CNIG contiene en su categoría Modelos Digitales de Elevaciones los llamados mapas LIDAR. Se trata de unos modelos digitales de superficies sombreados a los que se les han superpuesto dos capas procedentes de la rasterización de las clases edificación y vegetación de las nubes de puntos en formato .las y una tercera de hidrografía
Tenemos los derivados de la primera cobertura LIDAR que cubren todo el territorio nacional y ya llevan un tiempo disponibles para descarga. Se ofrecen en el Sistema Geográfico de Referencia WGS84, proyección Web Mercator (EPSG:3857). Y el motivo de este artículo es que el IGN ha comenzado la publicación del mapa LIDAR derivado de la segunda cobertura, para el que en esta ocasión ha optado por un Sistema Geográfico de Referencia ETRS89 en la Península, Islas Baleares, Ceuta y Melilla, y REGCAN95 en las Islas Canarias (ambos sistemas compatibles con WGS84) y proyección UTM en el huso correspondiente.
Datos comunes para ambas coberturas son el formato de archivo de descarga ECW y la unidad de dicha descarga (hojas del MTN50). También el proceso de elaboración que el mismo IGN nos especifica:
Este producto se genera por composición de cuatro capas, debidamente superpuestas. A partir de las nubes de puntos LiDAR correspondientes a la segunda cobertura del proyecto PNOA-LiDAR, debidamente clasificadas, se genera un modelo digital de superficies al que se le aplica un sombreado. Sobre este modelo se superponen tres capas: dos modelos digitales de superficies normalizados, uno de la clase vegetación y otro de la clase edificación, y la rasterización de la clase "agua" de las nubes de puntos LiDAR. Se aplica una paleta de color a cada capa, estableciendo un degradado en la información correspondiente a la vegetación y a la edificación en función de su altura sobre el terreno.
A fecha de este artículo las hojas disponibles de la segunda cobertura del mapa LIDAR son 669, con esta extensión geográfica (similar en cierta lógica a la del MDT02 en misma fecha):

La información auxiliar proporciona una cuadrícula en formato shape con las hojas disponibles, su nomenclatura y las fechas de vuelo.

En lo que respecta al mapa propiamente dicho se trata de una composición muy atractiva. La combinación del relieve sombreado con el coloreado de los elementos más importantes del paisaje (edificaciones, vegetación y aguas) resultan en un lienzo realista que no deja indiferente y con el que es más sencillo realizar escenas 3D visualmente buenas.

El mapa LIDAR cuenta así mismo con sendos servicios web disponibles que podemos usar en nuestros softwares favoritos. Además se puede consultar la altura de los edificios y vegetación a partir de cierta escala (1:20.000):
WMS........... http://wms-mapa-lidar.idee.es/lidar?request=GetCapabilities&service=WMS
WMTS......... http://wmts-mapa-lidar.idee.es/lidar?request=GetCapabilities&service=WMTS

Gracias a estos servicios web y SASPlanet pudimos elaborar por ejemplo los fondos JNX para Garmin que tenemos subidos en nuestra sección Garmin JNX del blog y que son una buena alternativa para complementar con nuestros mapas vectoriales en BaseCamp.
Ejemplo con TopoValencia para Garmin y su fondo LIDAR en BaseCamp.
Por último decir sobre los archivos físicos a descarga, los ECW's, que están elaborados a una resolución de 2,5 metros y con una compresión 1:6. La nomenclatura de estos archivos es un libro abierto, por ejemplo: MapaLidar-COB2-H30-0509.ecw Es decir nombre de producto, cobertura 2, Huso 30, hoja del MTN50 y formato. Una imagen de su aspecto a zoom 100%:

En definitiva un mapa LIDAR fantástico como fondo o mapa base para aplicaciones web y para el trabajo en entornos SIG. Un mapa además que cuenta con todas las versiones deseables en un producto TOP: servicios web y archivos físicos, para que cada usuario elija la forma de trabajo que mejor le convenga. Enhorabuena y gracias una vez más a nuestro IGN.