[pgrouting-users] Re: Pgrouting-users Digest, Vol 40, Issue 2
woodbri at swoodbridge.com
Mon Jan 9 20:55:49 EST 2012
On 1/9/2012 8:26 PM, Steve Horn wrote:
> > Hello all!
> > I am creating drive time radii (or isochrones) by allowing the
> user to
> > click on a map surface and entering a number of minutes.
> > Using the driving_distance function in pgrouting I am able to get
> > close to a good solution, but it is somewhat limited in that you can
> > only route from node to node (vertex to vertex). In order to
> improve the
> > accuracy, I thought I would try to insert a new node/vertex as
> close as
> > possible to the coordinates that my user clicked. That way my
> start node
> > would always be as accurate as it can be.
> > Has anyone done this? I'm not sure I understand the assign_vertex_id
> > function and how the graph is constructed well enough to come up
> with a
> > solution myself.
> > Any help or direction is appreciated!
> What you would want to do in this case is:
> 1. take you position and find the nearest edge to that location (EID,
> ESOURCE, ETARGET, ECOST, EREVERSE_COST)
> 2. compute the ratio of distance along that edge (RATIO)
> 3. then add a create a temporary node at that location (NEWNODE)
Right as Daniel noted, so the NEWNODE is just a number that would be
MAX_NODE_ID+1 If you have the vertices_tmp then you can get this number with
select max(id)+1 as newnode from vertices_tmp;
> 4. then add edges from the temp node to each of the end points to that
> edge, paying attention to oneway-ness of that edge and assigning costs
> as a ratio of the original cost.
> So select edges for your query might look like and you would pass
> NEWNODE as the starting node.
> select gid, source, target, cost, reverse_cost from roads where the_geom && <bbox>
> union all
> select maxgid+1, ESOURCE, NEWNODE, ECOST*RATIO, EREVERSE_COST*(1-RATIO)
> union all
> select maxgid+2, NEWNODE, ETARGET, ECOST*(1-RATIO), EREVERSE_COST*RATIO
This query is providing all the edges in the "roads" table and adding
two additional edges for the split nearest edge with NEWNODE at you
> I get lost on your 3rd step: "Create a temporary node"
> How do you insert just a node? The "roads" table is made up of edges, so
> wouldn't I need to insert an edge?
> Pgrouting-users mailing list
> Pgrouting-users at lists.osgeo.org
More information about the Pgrouting-users