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

Pierre Racine Pierre.Racine at sbf.ulaval.ca
Thu Aug 29 06:22:12 PDT 2013


Thanks Jean. I'll do this ASAP...

> -----Original Message-----
> From: postgis-users-bounces at lists.osgeo.org [mailto:postgis-users-
> bounces at lists.osgeo.org] On Behalf Of Jean Marchal
> Sent: Wednesday, August 28, 2013 6:52 PM
> To: PostGIS Users Discussion
> Subject: Re: [postgis-users] ST_Split fails on a quite simple geometry
> 
> 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


More information about the postgis-users mailing list