Cartogramando - Hazlo fácil. Usa un mapa

El uso de un sistema de gestión de base de datos con soporte geospacial es fundamental para cualquier entorno GIS en todos los niveles. Este post tiene como propósito recopilar las funciones que bajo mi criterio son más utilizadas (o al menos son la que más utilizo) y por lo tanto, es altamente recomendable que cualquier gestor de información geográfica deba saber.

Si por contra no sabes que es esto de Postgis entra aquí y échale un vistazo antes de seguir leyendo.

También dejo el enlace al manual completo de Postgis 2.2 en PDF y un enlace a un [curso] (http://workshops.boundlessgeo.com/postgis-intro/) de inciación a Postgis.

####St_AsText y AsEWKT Ambas funciones obtienen la geometría en formato WKT. La diferencia es que la segunda incluye el SRID de la misma.

text ST_AsText ( geometry g1)
text ST_AsEWKT ( geometry g1)

select St_AsText(the_geom) from  myTable;
select St_AsTEWKT(the_geom) from  myTable;

####St_GeomFromText y GeomFromEWKT En este caso estas funciones generan una geometría en función a las cadenas en formato WKT o EWKT respectivamente.

geometry ST_GeomFromText(text WKT);
geometry ST_GeomFromText(text WKT, integer srid);
geometry ST_GeomFromEWKT(text EWKT);

select ST_GeomFromText('POINT (0 0)') from  myTable;
select ST_GeomFromText('POINT (0 0)', 4326) from  myTable;
select ST_GeomFromEWKT('SRID=4326;POINT(0 0)') from  myTable;

####St_Buffer Una de las funciones más utilizadas en los sistemas de información geográfica. Esta función genera un polígono alrededor de la geometría pasada por parámetro. Dicho polígono puede definirse según los parámetros adicionales de la función.

geometry ST_Buffer(geometry g1, float radio buffer);
geometry ST_Buffer(geometry g1, float radio buffer, integer num_seg_quarter);
geometry ST_Buffer(geometry g1, float radio buffer, text params);
geography ST_Buffer(geography g1, float radio buffer);

SELECT ST_Buffer( ST_GeomFromText('POINT(0 0)'), 100, 'quad_segs=8');
SELECT ST_Buffer( ST_GeomFromText('LINESTRING(50 50,150 150,150 50)' ), 10, 'endcap=round join=round');
SELECT ST_Buffer( ST_GeomFromText('LINESTRING(50 50,150 150,150 50)' ), 10, 'join=bevel');

####St_Distance Calcula la distancia existente entre las dos geometrías pasadas por parámetro

float ST_Distance(geometry g1, geometry g2);
float ST_Distance(geography gg1, geography gg2, boolean use_spheroid);

SELECT ST_Distance(ST_GeomFromText('POINT(0 -1)',4326),ST_GeomFromText('LINESTRING(0 -3, -2 0)', 4326));

####St_Intersects Calcula si dos geometrías intersectan espacialmente. Devuelve ‘TRUE’ si intersecta y ‘FALSE’ en caso contrario.

boolean ST_Intersects( geometry geomA , geometry geomB );

SELECT ST_Intersects('POINT(0 0)', 'LINESTRING ( 2 0, 0 2 )');

####St_Union Si estás familiarizado con un GIS de escritorio lo estarás con la función ‘merge’. St_Union realiza esta operación de unión de objetos espaciales.

geometry ST_Union(geometry set g1field);
geometry ST_Union(geometry g1, geometry g2);
geometry ST_Union(geometry[] g1_array)

####St_ClosestPoint Una de mis favoritas. Existen multitud de situaciones en las que se necesita ‘asociar’ un punto sobre una línea o algo similar. Esta función obtiene el punto más cercano sobre la geometría g1 desde a la geometría g2.

geometry ST_ClosestPoint(geometry g1, geometry g2);

SELECT ST_AsText(ST_ClosestPoint(ST_GeomFromText('LINESTRING(0 0, 10 0, 10 10, 20 10, 20 20)'),ST_GeomFromText('POINT(15 15)')));

####St_Transform Función que transforma geometrías entre diferentes sistemas de referencia.

geometry ST_ClosestPoint(geometry g1, geometry g2);

SELECT St_AsText(St_Trasnform(St_GeometryFromText(('POINT (0 0)'), 4326,25830))

####St_Interpolate Obtiene la posición del punto sobre la geometría indicada en el primer parámetro, que corresponde con el segundo argumento de la función. Este argumento indica el porcentaje sobre la geometría sobre la que se quiere obtener la posición. Esta función es muy usada para referenciación lineal.

geometry ST_Line_Interpolate_Point(geometry line, float fraction);

SELECT ST_AsEWKT(ST_Line_Interpolate_Point(the_line, 0.50)) FROM (SELECT ST_GeomFromEWKT('LINESTRING(0 0, 10 10, 20 20)') as the_line)

En el ejemplo se quiere obtener la posición sobre el LINESTRING que se encuentre en el punto medio de la misma (0.5).