I have a table that contains the shortest path between destinations for survey respondents on an osm2po created network. When calculating the shortest path result, I store the geometry in a column. Now, using the same osm2po network table and each shortest path geometry, I am trying to determine whether the shortest path crosses or intersects a highway. I created another table that just contains the osm2po network with clazz 11-16 and run an st_intersects(shortestpathgeom,highwaygeom)
and get all false results. The shortest path geometry is multilinestring and the highway geometry is linestring, both are SRID 4326 and have been created from the same table. Why am I getting no intersections? Here is a picture of the shortest path in purple, highway in green, osm2po in red.
here is my code where po_routing.route is the shortest path, po_highway.the_geom is highway only network, and obstacle is the intersection test.
update analysis.po_routing
set obstacle = st_intersects(route,b.the_geom)
from analysis.po_highway as b;
I don't understand why this fundamental function is not working as expected.