[postgis-users] ST_Split fails on a quite simple geometry

[Jean Marchal]

Hi Pierre,

seems to work with the following version of PostGIS

SELECT postgis_full_version()

"POSTGIS="2.1.0 r11822" GEOS="3.4.1-CAPI-1.8.0 r0" PROJ="Rel. 4.8.0, 6
March 2012" GDAL="GDAL 1.10.0, released 2013/04/24" LIBXML="2.7.8"
LIBJSON="UNKNOWN" RASTER"

Ouput of* ST_AsText (*SELECT ST_Split(ST_GeomFromEWKT('
POLYGON((252049.803755834 5230978.87794769,
                                          252049.803755834 5230977.67794767,
                                          252050.003755837 5230977.67794767,
                                          252050.003755837 5230978.27794768,
                                          252050.20375584  5230978.27794768,
                                          252050.20375584  5230978.47794768,
                                          252050.003755837 5230978.47794768,
                                          252050.003755837 5230978.87794769,
                                          252049.803755834
5230978.87794769))'),
                           ST_MakeLine(ST_MakePoint(252050.003755837,
 5230978.87794769),

 ST_MakePoint(252049.803755834,  5230977.67794767))
                          )*) *

"GEOMETRYCOLLECTION(POLYGON((252049.803755834
5230977.67794767,252049.803755834 5230978.87794769,252050.003755837
5230978.87794769,252049.803755834
5230977.67794767)),POLYGON((252050.003755837
5230978.87794769,252050.003755837 5230978.47794768,252050.20375584
5230978.47794768,252050.20375584 5230978.27794768,252050.003755837
5230978.27794768,252050.003755837 5230977.67794767,252049.803755834
5230977.67794767,252050.003755837 5230978.87794769)))"

2 polygons are returned here. Try an upgrade ?

HTH,

Jean


2013/8/28 Pierre Racine <Pierre.Racine at sbf.ulaval.ca>

> Hi,
>
> This query should return the geometry splitted in two but it returns only
> one of the two sides.
>
> SELECT ST_Split(ST_GeomFromEWKT('POLYGON((252049.803755834
> 5230978.87794769,
>                                           252049.803755834
> 5230977.67794767,
>                                           252050.003755837
> 5230977.67794767,
>                                           252050.003755837
> 5230978.27794768,
>                                           252050.20375584
>  5230978.27794768,
>                                           252050.20375584
>  5230978.47794768,
>                                           252050.003755837
> 5230978.47794768,
>                                           252050.003755837
> 5230978.87794769,
>                                           252049.803755834
> 5230978.87794769))'),
>                            ST_MakeLine(ST_MakePoint(252050.003755837,
>  5230978.87794769),
>
>  ST_MakePoint(252049.803755834,  5230977.67794767))
>                           ) geom
>
> SELECT postgis_full_version()
>
> "POSTGIS="2.1.0SVN r11202" GEOS="3.4.0dev-CAPI-1.8.0 r0" PROJ="Rel. 4.8.0,
> 6 March 2012" GDAL="GDAL 1.9.2, released 2012/10/08" LIBXML="2.7.8"
> LIBJSON="UNKNOWN" RASTER"
>
> Bug? Or it works on earlier version?
>
> Thanks,
>
> Pierre
> _______________________________________________
> postgis-users mailing list
> postgis-users at lists.osgeo.org
> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users
>



-- 
Jean Marchal

Etudiant au doctorat / PHD student
Université Laval - Pavillon Abitibi-Price
Faculté de foresterie, de géographie et de géomatique
2405, rue de la Terrasse
Québec (Québec) G1V 0A6
Tél: 418-656-2131 poste 2620
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20130828/7b29e688/attachment.html>