Greetings all, I am working on a set of data where I have noded linestrings. All linestrings have just two points and the nodes were derived from the endpoints.
I've found one instance where st_intersects and st_equals return the incorrect value (AFAICT). ============================== select a.fid, b.fid, st_astext(a.g), st_astext(b.g), st_equals(st_startpoint(b.g),a.g ) start_equals, st_equals(st_endpoint(b.g),a.g )end_equals, st_astext(st_endpoint(b.g)) = st_astext(a.g) end_text_equals, st_astext(st_startpoint(b.g)) = st_astext(a.g) start_text_equals from gn.primaryline_nodes a join gn.primaryline b on /*st_intersects(a.g,b.g) --*/ st_dwithin(a.g,b.g,.05) where a.fid = 6097; --------------------------------------- fid fid st_astext st_astext start_equals end_equals end_text_equals start_text_equals 6097 16202 POINT(767513.95 656360.7) LINESTRING(767499.7 656360.55,767513.95 656360.7) False False True False 6097 16201 POINT(767513.95 656360.7) LINESTRING(767513.95 656360.7,767516.95 656361.6) False False False True =================== version info: POSTGIS="2.1.1 r12113" GEOS="3.4.2-CAPI-1.8.2 r3924" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.10.0, released 2013/04/24" LIBXML="2.7.8" LIBJSON="UNKNOWN" TOPOLOGY RASTER PostgreSQL 9.3.0, compiled by Visual C++ build 1600, 64-bit I did a select into a new schema and the same issue occurred. I did an `update table set g = st_astext(g)::geometry` and it works fine after that. I've posted a dump here: https://gist.githubusercontent.com/rhysallister/3d3c59c204db726f5068/raw/gistfile1.txt Maybe someone can reproduce it. The offending fid from the primarylines_node table is 6097 Regards, Rhys Peace & Love|Live Long & Prosper
_______________________________________________ postgis-users mailing list postgis-users@lists.osgeo.org http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users