3ENGINE

Programación y otros cachivaches

Categoría: Tecnologia

Página 9/45

Tecnologia

Cómo crear un mapa con Leaflet y GeoJSON


Leaflet es una libreria de JavaScript de código abierto que permite construir aplicaciones de mapas web. Soporta HTML5, CSS3 y la mayoría de las plataformas móviles y de escritorio. Existen otras librerías, como OpenLayers o la API de Google Maps, pero Leaflet tiene la ventaja de ser muy compacta (carga rápido) y no requiere de conocimientos profundos de GIS.

leafleft

GeoJSON se encuentra dentro del grupo de los formatos de intercambio de datos geoespaciales mas extendidos. GeoJSON se basa en JSON. Permite la codificación de colecciones de estructuras de datos geográficos. Un objeto GeoJSON puede representar una geometría, una característica, o una colección de características. GeoJSON soporta los siguientes objetos geométricos: Point, MultiPoint, LineString, MultiLineString, Polygon, MultiPolygon, GeometryCollection, Feature, FeatureCollection.

geojson

Obtener los datos GeoJSON para la demo

Podemos crear nuestro propio geojson mediante la herramienta online geojson.io, pero para nuestra demo obtendremos los datos de ejemplo de algunos de los lugares que siguiendo la inicitiva Open Data, ofrecen datos públicos de forma libre para todo el mundo, sin restricciones de derechos de autor, de patentes o de otros mecanismos de control. Por ejemplo, descargamos los municipios de la isla de La Palma. Lo abrimos:

{
"type": "FeatureCollection",
"crs": {
"type": "name",
"properties": {
"name": "urn:ogc:def:crs:OGC:1.3:CRS84"
}
},
"features": [{
"type": "Feature",
"properties": {
"OBJECTID": 1,
"SUPERFICIE": 135741394.267612,
"PERIMETRO": 66918.002370,
"ID": 1,
"MUNICIPIO": "EL PASO",
"CODIGO": 38027,
[...]
},
"geometry": {
"type": "Polygon",
"coordinates": [[
[-17.884219, 28.616156], 
[-17.884263, 28.616166], 
[-17.884306, 28.616177],
[...]

«crs» indica el sistema de coordenadas de referencia utilizado. Usa WGS 84 como sistema de coordenadas, que es ampliamente utilizado por Google Maps, OpenStreetMaps, etc. «features» es una colección de características. Cada característica se compone de una colección de propiedades y una geometría, en este caso, un polígono.

Opcionalmente, para validar el GeoJson podemos utilizar por ejemplo geojsonlint. Esta herramienta online carga los datos en un mapa y valida que sean correctos.

Cómo crear un mapa con Leafleft y GeoJSON

Visualizar el contenido en un mapa con Leaflet

La manera mas sencilla de que Leaflet reconozca nuestros datos es convertirlos a código javascript. Aprovechando que JSON es compatible con javascript, esto es tan sencillo como editar el fichero, hacer que los datos GeoJSON estén asociados a una variable municipios y guardar el fichero con extensión .js:

var municipios = {
"type": "FeatureCollection",
[...]

El código es el siguiente, mas abajo comento el código:

<!DOCTYPE html>
<html>
<head>
<title>Municipios de La Palma</title>
<meta charset="utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.css"/>
</head>
<body>
<div id="map" style="width: 600px; height: 400px"></div>
<script src="demo.js" type="text/javascript"></script>
<script src="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.js"></script>
<script>
var map = L.map('map').setView([28.68, -17.85], 10);

L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', {
   maxZoom: 18,
   attribution: 'Map data © <a href="http://openstreetmap.org">OpenStreetMap</a> contributors, '+
   '<a href="http://creativecommons.org/licenses/by-sa/2.0/">CC-BY-SA</a>, '+
   'Imagery © <a href="http://cloudmade.com">CloudMade</a>',
   id: 'mapbox.light'
}).addTo(map);

function onEachFeature(feature, layer) {
  if (feature.properties && feature.properties.MUNICIPIO) {
     var popupContent = "<p>Municipio: " + feature.properties.MUNICIPIO + "</p>";
     layer.bindPopup(popupContent);
  };
};

function style(feature) {
   return {
     weight: 2,
     opacity: 1,
     color: 'white',
     dashArray: '3',
     fillOpacity: 0.7,
     fillColor: '#FEB24C'
   };
};

L.geoJson(municipios, {
   style: style,
   onEachFeature: onEachFeature
}).addTo(map);

</script>
</body>
</html>
  • Línea 11: aquí referenciamos al geojson con los municipios.
  • Línea 14: map es la clase principal de Leafleft. Indicamos las coordenadas iniciales y el nivel de zoom por defecto.
  • Línea 16: añadimos una capa (layer) con el mapa. En este caso usamos las imágenes (tiles) de OpenStreetMap (OSM). Además especificamos el nivel máximo de zoom. Existen otros servidores de mapas de OSM, hay una lista de ellos en la Wiki de OSM o utiliza un comparador de mapas como map compare service de bbbike.org
  • Línea 24: función que muestra un popup con el nombre del municipio.
  • Línea 31: función con el estilo de cada uno de los features de la capa GeoJSON.
  • Línea 42: aquí creamos una capa de tipo GeoJSON, donde le pasamos la variable con los datos, le indicamos el estilo y le pasamos la función a llamar cuando el ratón haga click en una feature.

Y así quedaría el resultado:

Cómo crear un mapa con Leaflet y GeoJSON

Leer directamente GeoJSON con jQuery

En el anterior ejemplo hemos convertido el GeoJSON en código javascript. Pero dependiendo de nuestras necesidades, a veces esto no es posible y necesitaremos leer directamente el fichero en formato GeoJSON. Para hacer esto haremos uso de jQuery. El código es el siguiente, mas abajo comento el código:

<!DOCTYPE html>
<html>
<head>
<title>Municipios de La Palma</title>
<meta charset="utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.css"/>
</head>
<body>
<div id="map" style="width: 600px; height: 400px"></div>
<script src="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.js"></script>
<script src="https://code.jquery.com/jquery-2.1.4.min.js"></script>
<script>
var map = L.map('map').setView([28.68, -17.85], 10);

L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', {
   maxZoom: 18,
   attribution: 'Map data © <a href="http://openstreetmap.org">OpenStreetMap</a> contributors, '+
   '<a href="http://creativecommons.org/licenses/by-sa/2.0/">CC-BY-SA</a>, '+
   'Imagery © <a href="http://cloudmade.com">CloudMade</a>',
   id: 'mapbox.light'
}).addTo(map);

function onEachFeature(feature, layer) {
  if (feature.properties && feature.properties.MUNICIPIO) {
     var popupContent = "<p>Municipio: " + feature.properties.MUNICIPIO + "</p>";
     layer.bindPopup(popupContent);
  };
};

function style(feature) {
   return {
     weight: 2,
     opacity: 1,
     color: 'white',
     dashArray: '3',
     fillOpacity: 0.7,
     fillColor: '#FEB24C'
   };
};

$.getJSON("demo.json", function(data){
   L.geoJson(data, {
      style: style,
      onEachFeature: onEachFeature
   }).addTo(map);
});

</script>
</body>
</html>
  • Línea 12: aquí referenciamos a la libreria jQuery.
  • Línea 42: esta función lee el GeoJson, y despues, como el ejemplo anterior, creamos una capa de tipo GeoJSON, donde le pasamos la variable con los datos, le indicamos el estilo y le pasamos la función a llamar cuando el ratón haga click en una feature.



Tecnologia

Cómo guardar coordenadas geográficas en SQL Server


Introducción

SQL Server soporta dos tipos de datos espaciales:

  • geography: Superficies elipsoidales. Tiene en cuenta la superficie curva de la Tierra para sus cálculos. Se trabaja en longitud y latitud.
  • geometry: Superficies planas. La Tierra se trata como una proyección plana, no tiene en cuenta la forma elipsoidal de la Tierra. Adecuado para distancias cortas (edificios, calles,..)

Los tipos geography y geometry se crean a partir de objetos vectoriales, especificados en formato Well-Known Text (WKT) o Well-Known Binary (WKB).

Los objetos vectoriales admitidos son:

  • Point: una ubicación.
  • MultiPoint: una serie de puntos.
  • LineString: una serie de cero o más puntos conectados por líneas.
  • MultiLineString: un conjunto de linestrings.
  • Polygon: una región contigua descrita por un conjunto de linestrings.
  • MultiPolygon: un conjunto de polígonos.
  • GeometryCollection: una recopilación de tipos geométricos.

El modo mas sencillo de crear estos objetos es mediante el uso de los métodos de geografía
estáticos de Open Geospatial Consortium (OGC) incluidos en SQL Server que devuelven una instancia de geography o geometry a partir de una representación WKT:

  • STGeomFromText
  • STPointFromText
  • STMPointFromText
  • STLineFromText
  • STMLineFromText
  • STPolyFromText
  • STMPolyFromText
  • STGeomCollFromText

El método mas polivalente es STGeomFromText. Un ejemplo:

geography::STGeomFromText('LINESTRING(-122.360 47.656, -122.343 47.656)', 4326);

Donde el primer parámetro es un WKT con el objeto vectorial y el segundo parámetro es el identificador de referencia espacial (SRID).

El SRID corresponde a un sistema de referencia espacial. Sirve para identificar unívocamente el sistema de coordenadas utilizado para definir columnas de información espacial u objetos espaciales individuales en una columna espacial. Para averiguar la lista de SRID soportadas por SQL Server tenemos la vista sys.spatial_reference_systems:

select * from sys.spatial_reference_systems

Normalmente esto se compone de un sistema de coordenadas geodésicas, el más utilizado es el WGS84, con un patrón matemático de tres dimensiones que representa la tierra por medio de un elipsoide, mas un sistema de coordenadas cartesiano que pasa el modelo 3D a uno en 2D llamado proyección.

proy

Los dos mas comunes son el EPSG:3857 y el EPSG:4326, el primero es el usado en Google Maps, OpenStreetMap y otros mientras que el segundo es usado en GPS. La proyección de Mercator es las más utilizada y es en la que se basan EPSG:3857 y EPSG:4326. Representa la superficie esférica terrestre sobre una superficie cilíndrica, tangente al ecuador, que al desplegarse genera un mapa terrestre plano.

413px-Usgs_map_mercator.svg

Caso práctico

Despues de esta pequeña introducción 😉 un caso práctico.

1. Para guardar las coordenadas gráficas, he creado una tabla con un campo de tipo geography:

CREATE TABLE ejemplo (
[id] [int] NOT NULL PRIMARY KEY,
[geom] [geography] NOT NULL 
)

2. Para obtener coordenadas de prueba. En bing maps he añadido a Mis Lugares un área y dos puntos.

Cómo guardar coordenadas geográficas en SQL Server

Uno de los puntos se encuentra dentro del área y otro fuera. Por si alguien siente curiosidad el área delimita una parte del monasterio budista de Plana Novella. En concreto la zona donde están los mantras.

Ahora exporto las coordenadas a KML.

Cómo guardar coordenadas geográficas en SQL Server

Del KML me quedo con las coordenadas del área y monto el siguiente script para guardar las coordenadas geográficas en SQL Server:

INSERT INTO ejemplo (id, geom)
VALUES (1, 'POLYGON (( 1.8534866159243046 41.291925369663, 1.8540418331903874 41.291659349930804, 1.8537119214815556 41.291254272345995, 1.8531727974695622 41.29151021718217, 1.8534866159243046 41.291925369663))')
GO

3. Ahora quiero saber si hay una relación espacial entre el área y los dos puntos. En concreto quiero saber si alguno de los dos puntos se encuentra dentro del área. Hay un buen artículo que habla de este tema y que explica los métodos OGC disponibles.
SELECT id 
FROM ejemplo 
WHERE geom.STIntersects(geography::STGeomFromText('POINT (1.8535724466127812 41.29165531932045)', 4326)) = 1

SELECT id
FROM ejemplo
WHERE geom.STIntersects(geography::STGeomFromText('POINT (1.85385139635033 41.291979782656405)', 4326)) = 1



Tecnologia

MySQL: Error Code: 2006 – Mysql Server Has Gone Away


Al ejecutar un script SQL muy grande MySQL puede devolver un Error Code: 2006 – Mysql Server Has Gone Away. o bien un Error Code: 1153 Got a packet bigger than ‘max_allowed_packet’ bytes

Esto pasa especialmente al intentar restaurar copias de seguridad generadas con mysqldump.
Por defecto esta herramienta genera copias de seguridad que permite restaurar la base de datos. Genera un script que incluye las sentencias de creación de las tablas y las sentencias inserts para cada una de los registros que forman parte de las mismas.

Al tener que procesar un script SQL muy grande es conveniente restaurar desde terminal:

mysql --user=usuario --password=password < archivo.sql

O bien ya dentro del terminal MySQL:
mysql -u=root [-p=password]
mysql> create database mi_base_de_datos;
mysql> use mi_base_de_datos;
mysql> source archivo.sql;

Para solucionar el error debemos modificar el fichero de configuración de MySQL y aumentar el tamaño del parámetro max_allowed_packet (el tamaño máximo es de 1000MB):
max_allowed_packet=16M

En sistemas Linux el fichero de configuración se encuentra en /etc/my.cnf. Sino está ahí utiliza el comando locate my.cnf. En Windows el fichero se llama my.ini. Para averiguar donde se encuentra abrir un terminal y escribir services.msc:
> services.msc

Ya dentro del terminal buscar el servicio de MySQL (MySQL56) click al botón derecho del ratón y seleccionar propiedades. En el campo Ruta de acceso al ejecutable tenemos la ruta al fichero de configuración:

Error Code: 2006 - Mysql Server Has Gone Away

Si utilizamos XAMPP. Click en botón Config para acceder al fichero de configuración:

Error Code: 2006 - Mysql Server Has Gone Away




Tecnologia

MySQL: Error 1064 – You have an error in your SQL syntax … near ‘TYPE=MyISAM’


Al ejecutar el siguiente script para crear una tabla en un MySQL 5.6:

CREATE TABLE my_table(
id INT(11) NOT NULL AUTO_INCREMENT ,
name VARCHAR(60) NOT NULL DEFAULT '',
PRIMARY KEY (id),
KEY id(id)
) TYPE=MYISAM;

MySQL retorna un error 1064 – You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near ‘TYPE=MyISAM’.

El problema está en TYPE=MYISAM; porque TYPE quedó desfasada a partir de MySQL 4.0 y se eliminó definitivamente en MySQL 5.5. Para que funcione hay que sustituir TYPE por ENGINE.

CREATE TABLE my_table(
id INT(11) NOT NULL AUTO_INCREMENT ,
name VARCHAR(60) NOT NULL DEFAULT '',
PRIMARY KEY (id),
KEY id(id)
) ENGINE=MYISAM;