cnt := 0;
    dist_full := 0;
    FOR p1 IN select osm_id, public_transport,tags,way, ST_Distance_Spheroid(ST_Transform(way,4326), ST_GeomFromText('POINT('||$1||' '||$2||')',4326), 'SPHEROID["WGS 84",6378137,298.257223563]') as dist from opengeo.ru_psk_point where public_transport='stop_position' and tags -> 'bus' = 'yes' order by dist
    LOOP
       FOR p2 IN select osm_id, public_transport,tags,way,ST_Distance_Spheroid(ST_Transform(way,4326), ST_GeomFromText('POINT('||$3||' '||$4||')',4326), 'SPHEROID["WGS 84",6378137,298.257223563]') as dist from opengeo.ru_psk_point where public_transport='stop_position' and tags -> 'bus' = 'yes' order by dist
       LOOP
          EXECUTE 'select count(*) from opengeo.ru_psk_rels rels where rels.members::text similar to ''%'||p1.osm_id||'%'' and rels.tags::text similar to ''%route%'' and rels.tags::text similar to ''%ref%'' and rels.members::text similar to ''%'||p2.osm_id||'%''' INTO cnt;
          IF cnt > 0 AND (dist_full > (p1.dist + p2.dist) OR dist_full = 0) THEN
               dist_full := p1.dist + p2.dist;
               start_id := p1.osm_id;
               end_id := p2.osm_id;
               start_way := p1.way;
               end_way := p2.way;
          END IF;
       END LOOP;
    END LOOP;

-- Table: ru_psk_rels

-- DROP TABLE ru_psk_rels;

CREATE TABLE ru_psk_rels
(
  id bigint NOT NULL,
  way_off smallint,
  rel_off smallint,
  parts bigint[],
  members text[],
  tags text[],
  pending boolean NOT NULL,
  CONSTRAINT ru_psk_rels_pkey PRIMARY KEY (id)
)
WITH (
  OIDS=FALSE
);
ALTER TABLE ru_psk_rels
  OWNER TO opengeo;

-- Index: ru_psk_rels_idx

-- DROP INDEX ru_psk_rels_idx;

CREATE INDEX ru_psk_rels_idx
  ON ru_psk_rels
  USING btree
  (id)
  WHERE pending;

-- Index: ru_psk_rels_parts

-- DROP INDEX ru_psk_rels_parts;

CREATE INDEX ru_psk_rels_parts
  ON ru_psk_rels
  USING gin
  (parts);

select count(*) from opengeo.ru_psk_rels rels where '||p1.osm_id||'=ANY(rels.parts) and 'route'=ANY(rels.tags) and 'ref'=ANY(rels.tags) and '||p2.osm_id||'=ANY(rels.parts)

select count(*) from opengeo.ru_psk_rels rels where array['||p1.osm_id||'::bigint,'||p2.osm_id||'::bigint]<@rels.parts and array[''route'',''ref'']<@rels.tags limit 1

CREATE INDEX idx_ru_psk_rels_gin_parts_tags
   ON opengeo.ru_psk_rels USING gin (parts,tags);

"Function Scan on route_bus  (cost=0.25..10.25 rows=1000 width=108) (actual time=49380.514..49380.540 rows=4 loops=1)"
"Total runtime: 49380.726 ms"

"Function Scan on route_bus  (cost=0.25..10.25 rows=1000 width=108) (actual time=3570.939..3570.965 rows=4 loops=1)"
"Total runtime: 3571.130 ms"

CREATE OR REPLACE FUNCTION route_bus(character varying,character varying,character varying,character varying) -- lon_start_point,lat_start_point,lon_end_point,lat_end_point
  RETURNS
   table (
    way_tags text[],
    line_id integer,
    way_geom geometry,
    way_type text,
    way_length_meters double precision
   ) as
$BODY$
DECLARE 
 p1 record;
 p2 record;
 r1 opengeo.ru_psk_rels%rowtype;
 r2 opengeo.ru_psk_line%rowtype;
 cnt bigint;
 dist_full double precision;
 start_id bigint;
 end_id bigint;
 start_way geometry;
 end_way geometry;
 is_start boolean;
 is_end boolean;
BEGIN
    cnt := 0;
    dist_full := 0;
    FOR p1 IN select osm_id, public_transport,tags,way, ST_Distance_Spheroid(ST_Transform(way,4326), ST_GeomFromText('POINT('||$1||' '||$2||')',4326), 'SPHEROID["WGS 84",6378137,298.257223563]') as dist from opengeo.ru_psk_point where public_transport='stop_position' and tags -> 'bus' = 'yes' order by dist 
    LOOP
       FOR p2 IN select osm_id, public_transport,tags,way,ST_Distance_Spheroid(ST_Transform(way,4326), ST_GeomFromText('POINT('||$3||' '||$4||')',4326), 'SPHEROID["WGS 84",6378137,298.257223563]') as dist from opengeo.ru_psk_point where public_transport='stop_position' and tags -> 'bus' = 'yes' order by dist
       LOOP
          select count(*) INTO cnt from opengeo.ru_psk_rels rels where array[p1.osm_id::bigint,p2.osm_id::bigint]<@rels.parts and array['route','ref']<@rels.tags limit 1;
          IF cnt > 0 AND (dist_full > (p1.dist + p2.dist) OR dist_full = 0) THEN
               dist_full := p1.dist + p2.dist;
               start_id := p1.osm_id;
               end_id := p2.osm_id;
               start_way := p1.way;
               end_way := p2.way;
          END IF;
       END LOOP;
    END LOOP;
    IF dist_full > 0 THEN
         FOR r1 IN select * from opengeo.ru_psk_rels rels where array[start_id::bigint,end_id::bigint]<@rels.parts and array['route','ref']<@rels.tags
         LOOP
            FOR r2 IN select * from opengeo.ru_psk_line where osm_id=('-'||r1.id)::bigint
            LOOP
               is_start := ST_contains(r2.way,start_way);
               is_end := ST_contains(r2.way,end_way);
               IF is_start IS TRUE AND is_end IS TRUE THEN
                    way_tags := r1.tags;
                    line_id := r1.id;
                    IF ST_Line_Locate_Point(ST_MakeLine(r2.way), start_way) > ST_Line_Locate_Point(ST_MakeLine(r2.way), end_way) THEN
                         way_geom := ST_Line_Substring(ST_MakeLine(r2.way), ST_Line_Locate_Point(ST_MakeLine(r2.way), end_way), ST_Line_Locate_Point(ST_MakeLine(r2.way), start_way));
                    ELSE
                         way_geom := ST_Line_Substring(ST_MakeLine(r2.way), ST_Line_Locate_Point(ST_MakeLine(r2.way), start_way), ST_Line_Locate_Point(ST_MakeLine(r2.way), end_way));
                    END IF;
                    way_type := GeometryType(way_geom);
                    way_length_meters := ST_Length_Spheroid(ST_Transform(way_geom,4326), 'SPHEROID["WGS 84",6378137,298.257223563]');
                    RETURN NEXT;
               END IF;
            END LOOP;
         END LOOP;
    END IF;
END;
$BODY$
  LANGUAGE plpgsql stable;

"Limit  (cost=24.02..24.03 rows=1 width=0) (actual time=0.310..0.319 rows=1 loops=1)"
"  ->  Aggregate  (cost=24.02..24.03 rows=1 width=0) (actual time=0.287..0.287 rows=1 loops=1)"
"        ->  Bitmap Heap Scan on ru_psk_rels rels  (cost=20.00..24.02 rows=1 width=0) (actual time=0.087..0.186 rows=7 loops=1)"
"              Recheck Cond: ('{2314760553,2314760557}'::bigint[] <@ parts)"
"              Filter: ('{route,ref}'::text[] <@ tags)"
"              ->  Bitmap Index Scan on idx_ru_psk_rels_gin_parts  (cost=0.00..20.00 rows=1 width=0) (actual time=0.048..0.048 rows=7 loops=1)"
"                    Index Cond: ('{2314760553,2314760557}'::bigint[] <@ parts)"
"Total runtime: 0.550 ms"