From a7c20a0065844daf7528a6cb2d61507e520bf44a Mon Sep 17 00:00:00 2001
From: smitadolly <smit402@gmail.com>
Date: Thu, 6 Feb 2014 10:34:05 +0530
Subject: [PATCH 1/4] Create point_to_point_dijkstra_routing.sql

This script has been written based on the fact that reverse length of the unidirectional link is 100000.0 in the routable table. Based on your case in the table, you can edit the value. Also this script is only fot shortest_path and not for shooting start etc.
---
 .../point_to_point_dijkstra_routing.sql       | 624 ++++++++++++++++++
 1 file changed, 624 insertions(+)
 create mode 100644 functions-v2.0/point_to_point_dijkstra_routing.sql

diff --git a/functions-v2.0/point_to_point_dijkstra_routing.sql b/functions-v2.0/point_to_point_dijkstra_routing.sql
new file mode 100644
index 0000000..91376ca
--- /dev/null
+++ b/functions-v2.0/point_to_point_dijkstra_routing.sql
@@ -0,0 +1,624 @@
+-------------------------------------------------------------------
+--Author:-Smita Kumari Dolly
+--Contributor to the Algorithm: Sugandha
+--Date:-3/12/2013
+--Tested with Postgresql9.1, postgis1.5, pgrouting1.05 for 8.4
+-------------------------------------------------------------------
+
+-------------------------------------------------------------------
+--prerequisites for using the functions written in this file
+--tbl should have the following columns
+--start_id
+--end_id
+--the_geom
+--osm_name
+--speed_in_kmh
+--length 
+--reverse_length
+--row_flag-default value = false
+-------------------------------------------------------------------
+
+
+--CREATE TYPE path_result AS (vertex_id integer, edge_id integer, cost float8);
+
+-------------------------------------------------------------------
+-- This function delete all the virtual edge and the 
+-- tbl - table name return to oroginal format
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION cleanup_virtual_values_in_table(tbl varchar) 
+	RETURNS void AS 
+	$$
+    DECLARE
+         row record;
+         
+    BEGIN
+    	
+        EXECUTE
+	        'DELETE FROM '||tbl||' where row_flag = true';
+    END;
+$$ 
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+-------------------------------------------------------------------
+-- This function calculates shortest path between two point geometry along the edges specified in 
+-- tbl - table name 
+-------------------------------------------------------------------
+
+
+CREATE OR REPLACE FUNCTION point_to_point_shortest_path(source_point_geom varchar, target_point_geom varchar, tbl varchar, distance double precision, directed boolean, has_reverse_cost boolean) 
+	 RETURNS SETOF path_result AS
+	$$
+    DECLARE
+    	 sql text;
+         row record;
+         result path_result%rowtype;
+         srid int;
+         closest_linestring_id_from_source_geom int;
+         closest_linestring_id_from_target_geom int;
+         fraction_at_which_source_geom_on_closest_line_substring double precision;
+         fraction_at_which_target_geom_on_closest_line_substring double precision;
+         first_part_road_geom_where_source_geom_snapped text;
+         first_part_road_geom_where_source_geom_snapped_geom geometry;
+         second_part_road_geom_where_source_geom_snapped text;
+         second_part_road_geom_where_source_geom_snapped_geom geometry;
+         first_part_road_geom_where_target_geom_snapped text;
+         first_part_road_geom_where_target_geom_snapped_geom geometry;
+         second_part_road_geom_where_target_geom_snapped text;
+         second_part_road_geom_where_target_geom_snapped_geom geometry;
+         start_id int;
+         end_id int;
+         virtual_start_id int;
+         virtual_gid int;
+         virtual_end_id int;
+         is_source_edge_unidirectional int;
+         is_target_edge_unidirectional int;
+         vertex_id integer; 
+         edge_id integer; 
+         cost float8;
+         
+         
+    BEGIN
+    
+    	 EXECUTE 'select * from cleanup_virtual_values_in_table('''||tbl||''')';
+    
+    	 FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+	      END LOOP;
+	 srid:= row.srid;
+	
+	 FOR row IN EXECUTE 'select fetch_closest_edge_within_distance_xy('''||source_point_geom||''', '||distance||', '''||tbl||''') as gid' LOOP
+		      END LOOP;
+	 closest_linestring_id_from_source_geom:= row.gid;
+	
+	 FOR row IN EXECUTE 'select fetch_closest_edge_within_distance_xy('''||target_point_geom||''', '||distance||', '''||tbl||''') as gid' LOOP
+	 		      END LOOP;
+	 closest_linestring_id_from_target_geom:= row.gid;
+	 
+	 FOR row IN EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_source_geom||','''||source_point_geom||''', '''||tbl||''') as fraction' LOOP
+	 	 		      END LOOP;
+	 fraction_at_which_source_geom_on_closest_line_substring:= row.fraction;
+	 
+	 FOR row IN EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_target_geom||','''||target_point_geom||''', '''||tbl||''') as fraction' LOOP
+	 	 	 		      END LOOP;
+	 fraction_at_which_target_geom_on_closest_line_substring:= row.fraction;
+	 
+	 
+	 FOR row IN EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
+	 	 	 	 		      END LOOP;
+	 first_part_road_geom_where_source_geom_snapped:= row.the_geom;
+	 
+	 FOR row IN EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom'
+	 			    	 LOOP
+	 				 END LOOP;
+	 first_part_road_geom_where_source_geom_snapped_geom:= row.geom;
+	 
+	 FOR row IN EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
+	 	 	 	 	 		      END LOOP;
+	 second_part_road_geom_where_source_geom_snapped:= row.the_geom;
+	 
+	 FOR row IN EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom'
+	 	 			    	 LOOP
+	 	 				 END LOOP;
+	 second_part_road_geom_where_source_geom_snapped_geom:= row.geom;
+	 
+	 FOR row IN EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
+	 	 	 	 	 		      END LOOP;
+	 first_part_road_geom_where_target_geom_snapped:= row.the_geom;
+	 
+	 FOR row IN EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom'
+	 	 	 			    	 LOOP
+	 	 	 				 END LOOP;
+	 first_part_road_geom_where_target_geom_snapped_geom:= row.geom;
+	 	 
+	 FOR row IN EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
+	 	 	 	 	 	 		      END LOOP;
+	 second_part_road_geom_where_target_geom_snapped:= row.the_geom;
+	 
+	 FOR row IN EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom'
+	 	 	 	 			    	 LOOP
+	 	 	 	 				 END LOOP;
+	 second_part_road_geom_where_target_geom_snapped_geom:= row.geom;
+	 
+	
+	
+	
+	
+	 FOR row IN EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_source_geom||')  as val' LOOP
+	 	 		      END LOOP;
+	 is_source_edge_unidirectional:= row.val;
+	 
+	 FOR row IN EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_target_geom||')  as val' LOOP
+	 	 	 		      END LOOP;
+	 is_target_edge_unidirectional:= row.val;
+	 
+	
+	
+	 
+	 --virtual segments creation of the source point geom
+	 FOR row IN EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as start_id' LOOP
+	 		      END LOOP;
+	 start_id:= row.start_id;
+	 
+	 FOR row IN EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as end_id' LOOP
+	 	 		      END LOOP;
+	 end_id:= row.end_id;
+	 
+	 FOR row IN EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_start_id' LOOP
+	 	 	 		      END LOOP;
+	 virtual_start_id:= row.virtual_start_id;
+	 
+	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
+	 	 	 	 		      END LOOP;
+	 virtual_gid:= row.virtual_gid;
+	 
+	 --update first part of the road for source
+	 IF is_source_edge_unidirectional = 1  THEN
+	  FOR row IN EXECUTE
+	 	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||start_id||', '||virtual_start_id||','''||tbl||''','||is_source_edge_unidirectional||' )' LOOP
+	 	 	 	 		      END LOOP;
+	
+	 ELSE
+	 FOR row IN EXECUTE
+	 	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||start_id||','''||tbl||''', '||is_source_edge_unidirectional||')' LOOP
+	 	 	 	 		      END LOOP;
+
+	 END IF;
+	 
+	 --update second part of the road for source
+	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
+	 	 	 	 	 		      END LOOP;
+	 virtual_gid:= row.virtual_gid;
+	 
+	 FOR row IN EXECUTE
+	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||second_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||end_id||','''||tbl||''', '||is_source_edge_unidirectional||')' LOOP
+	 	 	 	 	 		      END LOOP;
+	
+	  --virtual segments creation of the target point geom
+	 
+	  FOR row IN EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as start_id' LOOP
+	 	 		      END LOOP;
+	 start_id:= row.start_id;
+
+	 FOR row IN EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as end_id' LOOP
+				      END LOOP;
+	 end_id:= row.end_id;
+
+	 FOR row IN EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_end_id' LOOP
+					      END LOOP;
+	 virtual_end_id:= row.virtual_end_id;
+
+	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
+						      END LOOP;
+	 virtual_gid:= row.virtual_gid;
+
+	 --update first part of the road for target
+	  FOR row IN EXECUTE
+	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||first_part_road_geom_where_target_geom_snapped||''', '||start_id||', '||virtual_end_id||','''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
+	  	 	 	 	 		      END LOOP;
+          
+          FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
+	  	 	 	 	 		      END LOOP;
+	 virtual_gid:= row.virtual_gid;
+          
+           IF is_target_edge_unidirectional = 1  THEN
+   	 --update second part of the road for target
+	  FOR row IN EXECUTE
+	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||virtual_end_id||', '||end_id||','''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
+	  	 	 	 	 		      END LOOP;
+
+           ELSE
+           FOR row IN EXECUTE
+	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||end_id||','||virtual_end_id||', '''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
+	  	 	 	 	 		      END LOOP;
+
+           END IF;
+	
+	
+	 
+	 
+	 
+	FOR result IN EXECUTE 'SELECT * FROM '||
+						  'shortest_path(''SELECT gid AS id, 
+								  start_id::int4 AS source, 
+								  end_id::int4 AS target, 
+								  length::float8 AS cost,
+								  reverse_length::float8 as reverse_cost
+								  FROM '||tbl||''',
+								 '||virtual_start_id||',
+								'||virtual_end_id||',
+								'||directed||',
+								'||has_reverse_cost||')'||''
+								LOOP
+									
+		 	 						  RETURN NEXT result;
+		 	 	 				      END LOOP;
+		 	 	
+					      RETURN;
+	 
+	 
+        
+    END;
+$$ 
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+-------------------------------------------------------------------
+-- This function will check whether a row have a reverse_length column in
+-- tbl - table name return to oroginal format
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION has_reverse_length(tbl varchar, gid int) 
+	RETURNS int AS 
+	$$
+    DECLARE
+         row record;
+         rownum int;
+         
+    BEGIN
+       rownum:= 0;
+       FOR row IN EXECUTE 'SELECT true as unidirectional FROM '||tbl||' WHERE gid = '||gid||' AND reverse_length = 1000000.0' LOOP
+       rownum:=rownum+1;
+       END LOOP;
+       
+       Return rownum;
+    END;
+$$ 
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+-------------------------------------------------------------------
+-- This function insert a virtual edge values in table
+-- tbl - table name 
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION create_update_virtual_edge(gid int, edge_id int, virtual_geom varchar, start_id int, end_id int, tbl varchar, is_unidirectional int) 
+	RETURNS void AS 
+	$$
+    DECLARE
+         row record;
+         osm_name varchar;
+         speed_in_kmh int;
+         srid int;
+         length double precision;
+         reverse_length double precision;
+         geom geometry;
+         
+    BEGIN
+    
+    	 FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+	      END LOOP;
+	srid:= row.srid;
+    
+    	 FOR row IN EXECUTE 'select osm_name as osm_name from '||tbl||' where gid = '||edge_id||'' 
+    	 LOOP
+	 END LOOP;
+	osm_name:= row.osm_name;
+	
+	 FOR row IN EXECUTE 'select speed_in_kmh as speed_in_kmh from '||tbl||' where gid = '||edge_id||'' 
+	    	 LOOP
+		 END LOOP;
+	speed_in_kmh:= row.speed_in_kmh;
+	
+	 FOR row IN EXECUTE 'select ST_Length(GeometryFromText('''||virtual_geom||''', '||srid||')) as length'
+		    	 LOOP
+			 END LOOP;
+	length:= row.length;
+	
+	FOR row IN EXECUTE 'select GeometryFromText('''||virtual_geom||''', '||srid||') as virtual_multigeom'
+			    	 LOOP
+				 END LOOP;
+	geom:= row.virtual_multigeom;
+	
+	IF is_unidirectional THEN
+	
+	reverse_length:= 100000;
+	
+	elSE
+	
+	reverse_length:= length;
+	
+	END IF;
+	
+	
+	
+	
+	
+        EXECUTE
+	        'INSERT INTO '||tbl||' (gid, the_geom, start_id, end_id, osm_name, speed_in_kmh, length, row_flag, reverse_length) VALUES ('||gid||', ST_GeomFromText('''||virtual_geom||''', '||srid||') ,'||start_id||','||end_id||', '''||osm_name||''', '||speed_in_kmh||', '||length||', true, '||reverse_length||')';
+        
+    END;
+$$ 
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+-------------------------------------------------------------------
+-- This function fetch the max of gid in
+-- tbl - table name 
+-- and increment it by 1
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION fetch_virtual_gid(tbl varchar)
+	RETURNS INT AS
+$$
+DECLARE
+    row record;
+    
+BEGIN
+	
+    -- get start id of specified edge
+
+    FOR row in EXECUTE 'select (max(gid)+1) as gid FROM '||tbl||''
+    LOOP
+    END LOOP;
+
+    IF row.gid IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.gid;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+-------------------------------------------------------------------
+-- This function fetch the max +1  of start_id and end_id column values in
+-- tbl - table name 
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION fetch_virtual_start_id_or_end_id(tbl varchar)
+	RETURNS INT AS
+$$
+DECLARE
+    row record;
+    
+BEGIN
+	
+    -- get start id of specified edge
+
+    FOR row in EXECUTE 'SELECT GREATEST(max(start_id), max(end_id))+1 as maxVal FROM '||tbl||''
+    LOOP
+    END LOOP;
+
+    IF row.maxVal IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.maxVal;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+-------------------------------------------------------------------
+-- This function fetch the end id of a specified edge
+-- edge_id - existing in 
+-- tbl - table name 
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION fetch_end_id_of_an_edge(edge_id int, tbl varchar)
+	RETURNS INT AS
+$$
+DECLARE
+    row record;
+    
+BEGIN
+	
+    -- get start id of specified edge
+
+    FOR row in EXECUTE 'SELECT end_id as end_id FROM '||tbl||' WHERE gid = '||edge_id||''
+    LOOP
+    END LOOP;
+
+    IF row.end_id IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.end_id;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+-------------------------------------------------------------------
+-- This function fetch the start id of a specified edge
+-- edge_id - existing in 
+-- tbl - table name 
+-------------------------------------------------------------------
+CREATE OR REPLACE FUNCTION fetch_start_id_of_an_edge(edge_id int, tbl varchar)
+	RETURNS INT AS
+$$
+DECLARE
+    row record;
+    
+BEGIN
+	
+    -- get start id of specified edge
+
+    FOR row in EXECUTE 'SELECT start_id as start_id FROM '||tbl||' WHERE gid = '||edge_id||''
+    LOOP
+    END LOOP;
+
+    IF row.start_id IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.start_id;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+-------------------------------------------------------------------
+-- This function gets the second half of the geom of an edge with
+-- edge_id - existing in 
+-- tbl - table name at which 
+-- fraction - divides the edge
+-------------------------------------------------------------------
+
+CREATE OR REPLACE FUNCTION fetch_second_half_of_the_geom_based_on_a_fraction(edge_id int, fraction double precision, tbl varchar)
+	RETURNS text AS
+$$
+DECLARE
+    row record;
+    srid integer;
+    
+BEGIN
+
+   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+      END LOOP;
+	srid:= row.srid;
+	
+    -- getting the fraction at which a point lie on an edge
+
+    FOR row in EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, '||fraction||', 1))) as geom_text
+						  FROM '||tbl||'
+						  WHERE gid = '||edge_id||''
+    LOOP
+    END LOOP;
+
+    IF row.geom_text IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.geom_text;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+-------------------------------------------------------------------
+-- This function gets the first half of the geom of an edge with
+-- edge_id - existing in 
+-- tbl - table name at which 
+-- fraction - divides the edge
+-------------------------------------------------------------------
+
+CREATE OR REPLACE FUNCTION fetch_first_half_of_the_geom_based_on_a_fraction(edge_id int, fraction double precision, tbl varchar)
+	RETURNS text AS
+$$
+DECLARE
+    row record;
+    srid integer;
+    
+BEGIN
+
+   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+      END LOOP;
+	srid:= row.srid;
+	
+    -- getting the fraction at which a point lie on an edge
+
+    FOR row in EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, 0, '||fraction||'))) as geom_text
+						  FROM '||tbl||'
+						  WHERE gid = '||edge_id||''
+    LOOP
+    END LOOP;
+
+    IF row.geom_text IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.geom_text;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+-------------------------------------------------------------------
+-- This function finds the fraction at which a node
+-- point - text representation of point lie on an edge with an 
+-- edge_id - existing in a 
+-- tbl - table name
+-------------------------------------------------------------------
+
+CREATE OR REPLACE FUNCTION fetch_snapped_loc_fraction_of_point_on_edge(edge_id int, point varchar, tbl varchar)
+	RETURNS double precision AS
+$$
+DECLARE
+    row record;
+    srid integer;
+    
+BEGIN
+
+   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+      END LOOP;
+	srid:= row.srid;
+	
+    -- getting the fraction at which a point lie on an edge
+
+    FOR row in EXECUTE 'SELECT ST_Line_Locate_Point(ST_LineMerge(the_geom),GeometryFromText('''||point||''', '||srid||')) as snapped_fraction
+						  FROM 
+						  '||tbl||'
+						  WHERE gid = '||edge_id||''
+    LOOP
+    END LOOP;
+
+    IF row.snapped_fraction IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.snapped_fraction;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+
+-------------------------------------------------------------------
+-- This function finds nearest edge to a given node
+-- point - text representation of point
+-- distance - function will search for an edge within this distance
+-- tbl - table name
+-------------------------------------------------------------------
+
+
+CREATE OR REPLACE FUNCTION fetch_closest_edge_within_distance_xy(point varchar,  distance double precision, tbl varchar)
+	RETURNS INT AS
+$$
+DECLARE
+    row record;
+    srid integer;
+    
+BEGIN
+
+   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
+      END LOOP;
+	srid:= row.srid;
+	
+    -- Searching for a link within the distance
+
+    FOR row in EXECUTE 'SELECT t1.gid  FROM  '||tbl||' as t1 WHERE
+						 ST_Expand(GeometryFromText('''||point||''', '||srid||'),'||distance||') && t1.the_geom 
+						 ORDER BY
+						 ST_Distance(GeometryFromText('''||point||''', '||srid||'),t1.the_geom)
+						 ASC LIMIT 1'
+    LOOP
+    END LOOP;
+
+    IF row.gid IS NULL THEN
+	    --RAISE EXCEPTION 'Data cannot be matched';
+	    RETURN NULL;
+    END IF;
+
+    RETURN row.gid;
+
+END;
+$$
+LANGUAGE 'plpgsql' VOLATILE STRICT;
+
+
+

From fe0b22772395bb6bdd5bc8f6fe2de0909ef88223 Mon Sep 17 00:00:00 2001
From: smitadolly <smit402@gmail.com>
Date: Thu, 6 Feb 2014 10:36:08 +0530
Subject: [PATCH 2/4] Update and rename null to README.md

---
 functions-v2.0/README.md | 63 ++++++++++++++++++++++++++++++++++++++++
 functions-v2.0/null      |  0
 2 files changed, 63 insertions(+)
 create mode 100644 functions-v2.0/README.md
 delete mode 100644 functions-v2.0/null

diff --git a/functions-v2.0/README.md b/functions-v2.0/README.md
new file mode 100644
index 0000000..556ac1f
--- /dev/null
+++ b/functions-v2.0/README.md
@@ -0,0 +1,63 @@
+point_to_point_routing
+======================
+
+This repository have library for point to point routing. 
+
+
+Prerequisites
+-------------
+
+1.	Install PostgreSQL - version 9.1
+
+2.	Install PostGIS - version 1.5
+
+3.	Install Pgrouting - version 1.05
+
+
+Installation
+-------------
+Start pgAdmin and create a new database based on your PostGIS template. Open a Query dialog, load and execute the 
+
+sql file available here. 
+
+
+Attributes Requirement in the routable table are as below:
+-----------------------------------------------------------
+1. start_id
+
+2. end_id
+
+3. the_geom
+
+4. osm_name
+
+5. speed_in_kmh
+
+6. length
+
+7. reverse_length
+
+8. row_flag (default value ‘false’)
+
+
+
+
+Function for point to point routing:
+------------------------------------
+point_to_point_shortest_path(source_point_geom varchar, target_point_geom varchar, tbl varchar, distance double precision, directed boolean, has_reverse_cost boolean) )
+
+SQL Query to get the point to point routing
+-------------------------------------------
+SELECT * from point_to_point_shortest_path('source_point_geometry_in_text_format','target_point_geometry_in_text_format','road_networks',#{0.003},false, false)"
+
+
+
+
+
+
+
+
+	
+
+
+
diff --git a/functions-v2.0/null b/functions-v2.0/null
deleted file mode 100644
index e69de29..0000000

From af91a18513ea690a99d73810eff3434088e5e83b Mon Sep 17 00:00:00 2001
From: smitadolly <smit402@gmail.com>
Date: Thu, 6 Feb 2014 13:16:39 +0530
Subject: [PATCH 3/4] Update point_to_point_dijkstra_routing.sql

Made changes based on the feedback got from Stephen
---
 .../point_to_point_dijkstra_routing.sql       | 341 ++++++------------
 1 file changed, 107 insertions(+), 234 deletions(-)

diff --git a/functions-v2.0/point_to_point_dijkstra_routing.sql b/functions-v2.0/point_to_point_dijkstra_routing.sql
index 91376ca..8439644 100644
--- a/functions-v2.0/point_to_point_dijkstra_routing.sql
+++ b/functions-v2.0/point_to_point_dijkstra_routing.sql
@@ -1,8 +1,5 @@
 -------------------------------------------------------------------
---Author:-Smita Kumari Dolly
---Contributor to the Algorithm: Sugandha
---Date:-3/12/2013
---Tested with Postgresql9.1, postgis1.5, pgrouting1.05 for 8.4
+--Tested with Postgresql8.4, postgis1.5, pgrouting1.05 for 9.1
 -------------------------------------------------------------------
 
 -------------------------------------------------------------------
@@ -78,163 +75,83 @@ CREATE OR REPLACE FUNCTION point_to_point_shortest_path(source_point_geom varcha
          
     BEGIN
     
-    	 EXECUTE 'select * from cleanup_virtual_values_in_table('''||tbl||''')';
+		EXECUTE 'select * from cleanup_virtual_values_in_table('''||tbl||''')';
     
-    	 FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-	      END LOOP;
-	 srid:= row.srid;
-	
-	 FOR row IN EXECUTE 'select fetch_closest_edge_within_distance_xy('''||source_point_geom||''', '||distance||', '''||tbl||''') as gid' LOOP
-		      END LOOP;
-	 closest_linestring_id_from_source_geom:= row.gid;
-	
-	 FOR row IN EXECUTE 'select fetch_closest_edge_within_distance_xy('''||target_point_geom||''', '||distance||', '''||tbl||''') as gid' LOOP
-	 		      END LOOP;
-	 closest_linestring_id_from_target_geom:= row.gid;
-	 
-	 FOR row IN EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_source_geom||','''||source_point_geom||''', '''||tbl||''') as fraction' LOOP
-	 	 		      END LOOP;
-	 fraction_at_which_source_geom_on_closest_line_substring:= row.fraction;
-	 
-	 FOR row IN EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_target_geom||','''||target_point_geom||''', '''||tbl||''') as fraction' LOOP
-	 	 	 		      END LOOP;
-	 fraction_at_which_target_geom_on_closest_line_substring:= row.fraction;
-	 
-	 
-	 FOR row IN EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
-	 	 	 	 		      END LOOP;
-	 first_part_road_geom_where_source_geom_snapped:= row.the_geom;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
+		
+		EXECUTE 'select fetch_closest_edge_within_distance_xy('''||source_point_geom||''', '||distance||', '''||tbl||''') as gid' INTO closest_linestring_id_from_source_geom;
+		     
+		EXECUTE 'select fetch_closest_edge_within_distance_xy('''||target_point_geom||''', '||distance||', '''||tbl||''') as gid' INTO closest_linestring_id_from_target_geom;
+		
+		EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_source_geom||','''||source_point_geom||''', '''||tbl||''') as fraction'  INTO  fraction_at_which_source_geom_on_closest_line_substring;
 	 
-	 FOR row IN EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom'
-	 			    	 LOOP
-	 				 END LOOP;
-	 first_part_road_geom_where_source_geom_snapped_geom:= row.geom;
+		EXECUTE 'select fetch_snapped_loc_fraction_of_point_on_edge('||closest_linestring_id_from_target_geom||','''||target_point_geom||''', '''||tbl||''') as fraction' INTO fraction_at_which_target_geom_on_closest_line_substring;
+		
+		EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' INTO first_part_road_geom_where_source_geom_snapped;
 	 
-	 FOR row IN EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
-	 	 	 	 	 		      END LOOP;
-	 second_part_road_geom_where_source_geom_snapped:= row.the_geom;
+		EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom' INTO first_part_road_geom_where_source_geom_snapped_geom;
 	 
-	 FOR row IN EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom'
-	 	 			    	 LOOP
-	 	 				 END LOOP;
-	 second_part_road_geom_where_source_geom_snapped_geom:= row.geom;
+		EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_source_geom||', '||fraction_at_which_source_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' INTO second_part_road_geom_where_source_geom_snapped;
 	 
-	 FOR row IN EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
-	 	 	 	 	 		      END LOOP;
-	 first_part_road_geom_where_target_geom_snapped:= row.the_geom;
+		EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_source_geom_snapped||''', '||srid||') as geom' INTO second_part_road_geom_where_source_geom_snapped_geom;
 	 
-	 FOR row IN EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom'
-	 	 	 			    	 LOOP
-	 	 	 				 END LOOP;
-	 first_part_road_geom_where_target_geom_snapped_geom:= row.geom;
-	 	 
-	 FOR row IN EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' LOOP
-	 	 	 	 	 	 		      END LOOP;
-	 second_part_road_geom_where_target_geom_snapped:= row.the_geom;
-	 
-	 FOR row IN EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom'
-	 	 	 	 			    	 LOOP
-	 	 	 	 				 END LOOP;
-	 second_part_road_geom_where_target_geom_snapped_geom:= row.geom;
-	 
-	
-	
-	
-	
-	 FOR row IN EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_source_geom||')  as val' LOOP
-	 	 		      END LOOP;
-	 is_source_edge_unidirectional:= row.val;
+		EXECUTE 'select fetch_first_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' INTO first_part_road_geom_where_target_geom_snapped;
+		
+		EXECUTE 'select GeometryFromText('''||first_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom' INTO first_part_road_geom_where_target_geom_snapped_geom;
+	 	
+		EXECUTE 'select fetch_second_half_of_the_geom_based_on_a_fraction('||closest_linestring_id_from_target_geom||', '||fraction_at_which_target_geom_on_closest_line_substring||', '''||tbl||''') as the_geom' INTO second_part_road_geom_where_target_geom_snapped;
+		
+		EXECUTE 'select GeometryFromText('''||second_part_road_geom_where_target_geom_snapped||''', '||srid||') as geom' INTO second_part_road_geom_where_target_geom_snapped_geom;
+		
+		EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_source_geom||')  as val' INTO is_source_edge_unidirectional;
 	 
-	 FOR row IN EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_target_geom||')  as val' LOOP
-	 	 	 		      END LOOP;
-	 is_target_edge_unidirectional:= row.val;
-	 
-	
-	
-	 
-	 --virtual segments creation of the source point geom
-	 FOR row IN EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as start_id' LOOP
-	 		      END LOOP;
-	 start_id:= row.start_id;
-	 
-	 FOR row IN EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as end_id' LOOP
-	 	 		      END LOOP;
-	 end_id:= row.end_id;
-	 
-	 FOR row IN EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_start_id' LOOP
-	 	 	 		      END LOOP;
-	 virtual_start_id:= row.virtual_start_id;
-	 
-	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
-	 	 	 	 		      END LOOP;
-	 virtual_gid:= row.virtual_gid;
+		EXECUTE 'select has_reverse_length('''||tbl||''', '||closest_linestring_id_from_target_geom||')  as val' INTO is_target_edge_unidirectional;
+
+	--virtual segments creation of the source point geom
+		EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as start_id' INTO start_id;
+		
+		EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_source_geom||', '''||tbl||''') as end_id' INTO end_id;
+		
+		EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_start_id' INTO virtual_start_id;
+		
+		EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' INTO virtual_gid;
 	 
 	 --update first part of the road for source
-	 IF is_source_edge_unidirectional = 1  THEN
-	  FOR row IN EXECUTE
-	 	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||start_id||', '||virtual_start_id||','''||tbl||''','||is_source_edge_unidirectional||' )' LOOP
-	 	 	 	 		      END LOOP;
-	
-	 ELSE
-	 FOR row IN EXECUTE
-	 	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||start_id||','''||tbl||''', '||is_source_edge_unidirectional||')' LOOP
-	 	 	 	 		      END LOOP;
-
-	 END IF;
+	IF is_source_edge_unidirectional = 1  THEN
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||start_id||', '||virtual_start_id||','''||tbl||''','||is_source_edge_unidirectional||' )';
+	ELSE
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||first_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||start_id||','''||tbl||''', '||is_source_edge_unidirectional||')';
+	END IF;
 	 
 	 --update second part of the road for source
-	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
-	 	 	 	 	 		      END LOOP;
-	 virtual_gid:= row.virtual_gid;
-	 
-	 FOR row IN EXECUTE
-	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||second_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||end_id||','''||tbl||''', '||is_source_edge_unidirectional||')' LOOP
-	 	 	 	 	 		      END LOOP;
+		EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' INTO virtual_gid;
+		
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_source_geom||', '''||second_part_road_geom_where_source_geom_snapped||''', '||virtual_start_id||', '||end_id||','''||tbl||''', '||is_source_edge_unidirectional||')';
 	
-	  --virtual segments creation of the target point geom
+	
+	 --virtual segments creation of the target point geom
 	 
-	  FOR row IN EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as start_id' LOOP
-	 	 		      END LOOP;
-	 start_id:= row.start_id;
+		EXECUTE 'select fetch_start_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as start_id' INTO start_id;
 
-	 FOR row IN EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as end_id' LOOP
-				      END LOOP;
-	 end_id:= row.end_id;
+		EXECUTE 'select fetch_end_id_of_an_edge('||closest_linestring_id_from_target_geom||', '''||tbl||''') as end_id' INTO end_id;
 
-	 FOR row IN EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_end_id' LOOP
-					      END LOOP;
-	 virtual_end_id:= row.virtual_end_id;
+		EXECUTE 'select fetch_virtual_start_id_or_end_id('''||tbl||''') as virtual_end_id' INTO virtual_end_id;
 
-	 FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
-						      END LOOP;
-	 virtual_gid:= row.virtual_gid;
+		EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' INTO virtual_gid;
 
 	 --update first part of the road for target
-	  FOR row IN EXECUTE
-	 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||first_part_road_geom_where_target_geom_snapped||''', '||start_id||', '||virtual_end_id||','''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
-	  	 	 	 	 		      END LOOP;
-          
-          FOR row IN EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' LOOP
-	  	 	 	 	 		      END LOOP;
-	 virtual_gid:= row.virtual_gid;
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||first_part_road_geom_where_target_geom_snapped||''', '||start_id||', '||virtual_end_id||','''||tbl||''', '||is_target_edge_unidirectional||')';
           
-           IF is_target_edge_unidirectional = 1  THEN
-   	 --update second part of the road for target
-	  FOR row IN EXECUTE
-	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||virtual_end_id||', '||end_id||','''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
-	  	 	 	 	 		      END LOOP;
-
-           ELSE
-           FOR row IN EXECUTE
-	  'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||end_id||','||virtual_end_id||', '''||tbl||''', '||is_target_edge_unidirectional||')' LOOP
-	  	 	 	 	 		      END LOOP;
-
-           END IF;
-	
+		EXECUTE 'select fetch_virtual_gid('''||tbl||''') as virtual_gid' INTO virtual_gid;
+         
+	--update second part of the road for target	 
+    IF is_target_edge_unidirectional = 1  THEN
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||virtual_end_id||', '||end_id||','''||tbl||''', '||is_target_edge_unidirectional||')';
+	ELSE
+		EXECUTE 'select create_update_virtual_edge('||virtual_gid||','||closest_linestring_id_from_target_geom||', '''||second_part_road_geom_where_target_geom_snapped||''', '||end_id||','||virtual_end_id||', '''||tbl||''', '||is_target_edge_unidirectional||')';
+	END IF;
 	
-	 
-	 
-	 
+ 
 	FOR result IN EXECUTE 'SELECT * FROM '||
 						  'shortest_path(''SELECT gid AS id, 
 								  start_id::int4 AS source, 
@@ -298,46 +215,23 @@ CREATE OR REPLACE FUNCTION create_update_virtual_edge(gid int, edge_id int, virt
          
     BEGIN
     
-    	 FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-	      END LOOP;
-	srid:= row.srid;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
     
-    	 FOR row IN EXECUTE 'select osm_name as osm_name from '||tbl||' where gid = '||edge_id||'' 
-    	 LOOP
-	 END LOOP;
-	osm_name:= row.osm_name;
-	
-	 FOR row IN EXECUTE 'select speed_in_kmh as speed_in_kmh from '||tbl||' where gid = '||edge_id||'' 
-	    	 LOOP
-		 END LOOP;
-	speed_in_kmh:= row.speed_in_kmh;
+    	EXECUTE 'select osm_name as osm_name from '||tbl||' where gid = '||edge_id||'' INTO osm_name;
 	
-	 FOR row IN EXECUTE 'select ST_Length(GeometryFromText('''||virtual_geom||''', '||srid||')) as length'
-		    	 LOOP
-			 END LOOP;
-	length:= row.length;
+		EXECUTE 'select speed_in_kmh as speed_in_kmh from '||tbl||' where gid = '||edge_id||'' INTO speed_in_kmh;
 	
-	FOR row IN EXECUTE 'select GeometryFromText('''||virtual_geom||''', '||srid||') as virtual_multigeom'
-			    	 LOOP
-				 END LOOP;
-	geom:= row.virtual_multigeom;
+		EXECUTE 'select ST_Length(GeometryFromText('''||virtual_geom||''', '||srid||')) as length' INTO length;
+		
+		EXECUTE 'select GeometryFromText('''||virtual_geom||''', '||srid||') as virtual_multigeom' INTO geom;
 	
 	IF is_unidirectional THEN
-	
-	reverse_length:= 100000;
-	
+		reverse_length:= 100000;
 	elSE
-	
-	reverse_length:= length;
-	
+		reverse_length:= length;
 	END IF;
 	
-	
-	
-	
-	
-        EXECUTE
-	        'INSERT INTO '||tbl||' (gid, the_geom, start_id, end_id, osm_name, speed_in_kmh, length, row_flag, reverse_length) VALUES ('||gid||', ST_GeomFromText('''||virtual_geom||''', '||srid||') ,'||start_id||','||end_id||', '''||osm_name||''', '||speed_in_kmh||', '||length||', true, '||reverse_length||')';
+        EXECUTE 'INSERT INTO '||tbl||' (gid, the_geom, start_id, end_id, osm_name, speed_in_kmh, length, row_flag, reverse_length) VALUES ('||gid||', ST_GeomFromText('''||virtual_geom||''', '||srid||') ,'||start_id||','||end_id||', '''||osm_name||''', '||speed_in_kmh||', '||length||', true, '||reverse_length||')';
         
     END;
 $$ 
@@ -352,21 +246,20 @@ CREATE OR REPLACE FUNCTION fetch_virtual_gid(tbl varchar)
 $$
 DECLARE
     row record;
+	gid int;
     
 BEGIN
 	
     -- get start id of specified edge
 
-    FOR row in EXECUTE 'select (max(gid)+1) as gid FROM '||tbl||''
-    LOOP
-    END LOOP;
+		EXECUTE 'select (max(gid)+1) as gid FROM '||tbl||'' INTO gid;
 
-    IF row.gid IS NULL THEN
+    IF gid IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.gid;
+    RETURN gid;
 
 END;
 $$
@@ -380,21 +273,19 @@ CREATE OR REPLACE FUNCTION fetch_virtual_start_id_or_end_id(tbl varchar)
 $$
 DECLARE
     row record;
-    
+    maxVal int;
 BEGIN
 	
     -- get start id of specified edge
 
-    FOR row in EXECUTE 'SELECT GREATEST(max(start_id), max(end_id))+1 as maxVal FROM '||tbl||''
-    LOOP
-    END LOOP;
+		EXECUTE 'SELECT GREATEST(max(start_id), max(end_id))+1 as maxVal FROM '||tbl||'' INTO maxVal;
 
-    IF row.maxVal IS NULL THEN
+    IF maxVal IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.maxVal;
+    RETURN maxVal;
 
 END;
 $$
@@ -410,21 +301,19 @@ CREATE OR REPLACE FUNCTION fetch_end_id_of_an_edge(edge_id int, tbl varchar)
 $$
 DECLARE
     row record;
-    
+    end_id int;
 BEGIN
 	
     -- get start id of specified edge
 
-    FOR row in EXECUTE 'SELECT end_id as end_id FROM '||tbl||' WHERE gid = '||edge_id||''
-    LOOP
-    END LOOP;
+		EXECUTE 'SELECT end_id as end_id FROM '||tbl||' WHERE gid = '||edge_id||'' INTO end_id;
 
-    IF row.end_id IS NULL THEN
+    IF end_id IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.end_id;
+    RETURN end_id;
 
 END;
 $$
@@ -439,21 +328,19 @@ CREATE OR REPLACE FUNCTION fetch_start_id_of_an_edge(edge_id int, tbl varchar)
 $$
 DECLARE
     row record;
-    
+    start_id int;
 BEGIN
 	
     -- get start id of specified edge
 
-    FOR row in EXECUTE 'SELECT start_id as start_id FROM '||tbl||' WHERE gid = '||edge_id||''
-    LOOP
-    END LOOP;
+		EXECUTE 'SELECT start_id as start_id FROM '||tbl||' WHERE gid = '||edge_id||'' INTO start_id;
 
-    IF row.start_id IS NULL THEN
+    IF start_id IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.start_id;
+    RETURN start_id;
 
 END;
 $$
@@ -472,27 +359,24 @@ $$
 DECLARE
     row record;
     srid integer;
+	geom_text text;
     
 BEGIN
 
-   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-      END LOOP;
-	srid:= row.srid;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
 	
     -- getting the fraction at which a point lie on an edge
 
-    FOR row in EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, '||fraction||', 1))) as geom_text
+		EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, '||fraction||', 1))) as geom_text
 						  FROM '||tbl||'
-						  WHERE gid = '||edge_id||''
-    LOOP
-    END LOOP;
+						  WHERE gid = '||edge_id||'' INTO geom_text;
 
-    IF row.geom_text IS NULL THEN
+    IF geom_text IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.geom_text;
+    RETURN geom_text;
 
 END;
 $$
@@ -511,27 +395,24 @@ $$
 DECLARE
     row record;
     srid integer;
+	geom_text text;
     
 BEGIN
 
-   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-      END LOOP;
-	srid:= row.srid;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
 	
     -- getting the fraction at which a point lie on an edge
 
-    FOR row in EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, 0, '||fraction||'))) as geom_text
+		EXECUTE 'SELECT ST_AsText(ST_GeomFromEWKT(ST_Line_Substring(the_geom, 0, '||fraction||'))) as geom_text
 						  FROM '||tbl||'
-						  WHERE gid = '||edge_id||''
-    LOOP
-    END LOOP;
+						  WHERE gid = '||edge_id||'' INTO geom_text;
 
-    IF row.geom_text IS NULL THEN
+    IF geom_text IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.geom_text;
+    RETURN geom_text;
 
 END;
 $$
@@ -550,28 +431,24 @@ $$
 DECLARE
     row record;
     srid integer;
-    
+    snapped_fraction double precision;
 BEGIN
 
-   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-      END LOOP;
-	srid:= row.srid;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
 	
     -- getting the fraction at which a point lie on an edge
 
-    FOR row in EXECUTE 'SELECT ST_Line_Locate_Point(ST_LineMerge(the_geom),GeometryFromText('''||point||''', '||srid||')) as snapped_fraction
+		EXECUTE 'SELECT ST_Line_Locate_Point(ST_LineMerge(the_geom),GeometryFromText('''||point||''', '||srid||')) as snapped_fraction
 						  FROM 
 						  '||tbl||'
-						  WHERE gid = '||edge_id||''
-    LOOP
-    END LOOP;
+						  WHERE gid = '||edge_id||'' INTO snapped_fraction;
 
-    IF row.snapped_fraction IS NULL THEN
+    IF snapped_fraction IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.snapped_fraction;
+    RETURN snapped_fraction;
 
 END;
 $$
@@ -592,29 +469,25 @@ $$
 DECLARE
     row record;
     srid integer;
-    
+    gid int;
 BEGIN
 
-   FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' LOOP
-      END LOOP;
-	srid:= row.srid;
+		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
 	
     -- Searching for a link within the distance
 
-    FOR row in EXECUTE 'SELECT t1.gid  FROM  '||tbl||' as t1 WHERE
+		EXECUTE 'SELECT t1.gid  FROM  '||tbl||' as t1 WHERE
 						 ST_Expand(GeometryFromText('''||point||''', '||srid||'),'||distance||') && t1.the_geom 
 						 ORDER BY
 						 ST_Distance(GeometryFromText('''||point||''', '||srid||'),t1.the_geom)
-						 ASC LIMIT 1'
-    LOOP
-    END LOOP;
+						 ASC LIMIT 1' INTO gid;
 
-    IF row.gid IS NULL THEN
+    IF gid IS NULL THEN
 	    --RAISE EXCEPTION 'Data cannot be matched';
 	    RETURN NULL;
     END IF;
 
-    RETURN row.gid;
+    RETURN gid;
 
 END;
 $$

From 8b8544bcf634e20ef907ac09d7ff4a9ca482a0b8 Mon Sep 17 00:00:00 2001
From: smitadolly <smit402@gmail.com>
Date: Tue, 11 Feb 2014 12:39:51 +0530
Subject: [PATCH 4/4] Update point_to_point_dijkstra_routing.sql

replaced the single quote by double quote if osm_name have any single quote
---
 functions-v2.0/point_to_point_dijkstra_routing.sql | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/functions-v2.0/point_to_point_dijkstra_routing.sql b/functions-v2.0/point_to_point_dijkstra_routing.sql
index 8439644..5a39114 100644
--- a/functions-v2.0/point_to_point_dijkstra_routing.sql
+++ b/functions-v2.0/point_to_point_dijkstra_routing.sql
@@ -16,7 +16,7 @@
 -------------------------------------------------------------------
 
 
---CREATE TYPE path_result AS (vertex_id integer, edge_id integer, cost float8);
+
 
 -------------------------------------------------------------------
 -- This function delete all the virtual edge and the 
@@ -217,7 +217,7 @@ CREATE OR REPLACE FUNCTION create_update_virtual_edge(gid int, edge_id int, virt
     
 		EXECUTE 'select getsrid(the_geom) as srid from '||tbl||' where gid = (select min(gid) from '||tbl||')' INTO srid;
     
-    	EXECUTE 'select osm_name as osm_name from '||tbl||' where gid = '||edge_id||'' INTO osm_name;
+    		EXECUTE 'select replace(osm_name,'''''''','''''''''''') as osm_name from '||tbl||' where gid = '||edge_id||'' INTO osm_name;
 	
 		EXECUTE 'select speed_in_kmh as speed_in_kmh from '||tbl||' where gid = '||edge_id||'' INTO speed_in_kmh;
 	
@@ -231,6 +231,9 @@ CREATE OR REPLACE FUNCTION create_update_virtual_edge(gid int, edge_id int, virt
 		reverse_length:= length;
 	END IF;
 	
+	
+	
+	
         EXECUTE 'INSERT INTO '||tbl||' (gid, the_geom, start_id, end_id, osm_name, speed_in_kmh, length, row_flag, reverse_length) VALUES ('||gid||', ST_GeomFromText('''||virtual_geom||''', '||srid||') ,'||start_id||','||end_id||', '''||osm_name||''', '||speed_in_kmh||', '||length||', true, '||reverse_length||')';
         
     END;