pgRouting Manual (2.0.0)

pgr_alphaShape

«  pgr_drivingDistance   ::   Contents   ::   pgr_pointsAsPolygon  »


pgr_alphaShape

Name

pgr_alphashape — Core function for alpha shape computation.

Note

This function should not be used directly. Use pgr_drivingDistance instead.

Synopsis

Returns a table with (x, y) rows that describe the vertices of an alpha shape.

table() pgr_alphashape(text sql);

Description

sql:

text a SQL query, which should return a set of rows with the following columns:

SELECT id, x, y FROM vertex_table
id:int4 identifier of the vertex
x:float8 x-coordinate
y:float8 y-coordinate

Returns a vertex record for each row :

x:x-coordinate
y:y-coordinate

History

  • Renamed in version 2.0.0

Examples

In the alpha shape code we have no way to control the order of the points so the actual output you might get could be similar but different. The simple query is followed by a more complex one that constructs a polygon and computes the areas of it. This should be the same as the result on your system. We leave the details of the complex query to the reader as an exercise if they wish to decompose it into understandable pieces or to just copy and paste it into a SQL window to run.

SELECT * FROM pgr_alphashape('SELECT id, x, y FROM vertex_table');

 x | y
---+---
 2 | 0
 4 | 1
 4 | 2
 4 | 3
 2 | 4
 0 | 2
(6 rows)

SELECT round(st_area(ST_MakePolygon(ST_AddPoint(foo.openline, ST_StartPoint(foo.openline))))::numeric, 2) as st_area
from (select st_makeline(points order by id)  as openline from
(SELECT st_makepoint(x,y) as points ,row_number() over() AS id
FROM pgr_alphAShape('SELECT id, x, y FROM vertex_table')
) as a) as foo;

 st_area
---------
   10.00
(1 row)


SELECT * FROM pgr_alphAShape('SELECT id::integer, st_x(the_geom)::float as x, st_y(the_geom)::float as y  FROM edge_table_vertices_pgr');
  x  |  y
-----+-----
 0.5 | 3.5
   0 |   2
   2 |   0
   4 |   1
   4 |   2
   4 |   3
 3.5 |   4
   2 |   4
(8 rows)

SELECT round(st_area(ST_MakePolygon(ST_AddPoint(foo.openline, ST_StartPoint(foo.openline))))::numeric, 2) as st_area
from (select st_makeline(points order by id)  as openline from
(SELECT st_makepoint(x,y) as points ,row_number() over() AS id
FROM pgr_alphAShape('SELECT id::integer, st_x(the_geom)::float as x, st_y(the_geom)::float as y  FROM edge_table_vertices_pgr')
) as a) as foo;

 st_area
---------
   10.00
(1 row)

The queries use the Sample Data network.

See Also

«  pgr_drivingDistance   ::   Contents   ::   pgr_pointsAsPolygon  »