Re: [postgis-users] [postgis-devel] Problem in st_StriaghtSkeleton - Hangs!

2018-01-17 Thread Stephen Mather
For sure you can get a much less expensive result than the SFCGAL approach by using voronoi on a densified version of your geometry. It’s a few orders of magnitude faster that SFCGAL. Stephen V. Mather +1 (216) 339-6347 (Signal, Telegram, Cell, and WhatsApp) > On Jan 18, 2018, at 02:12,

Re: [postgis-users] Preserving Repeated Points

2015-03-17 Thread Stephen Mather
No problems here... On Mon, Mar 16, 2015 at 11:58 AM, Paul Ramsey pram...@cleverelephant.ca wrote: Users, Per this ticket: http://trac.osgeo.org/postgis/ticket/3061 We currently strip out repeated points from GML, KML and GeoJSON inputs. We allow them on WKT and WKB input. This seems

Re: [postgis-users] raster2pgsql: -t TILE_SIZE

2015-03-16 Thread Stephen Mather
Hi Ramon, I am not an expert on this by any measure. Better to check with the postgis user group (cc'd). I know the better size choices have changed a lot since Bborie optimized / rewrote in C. Best, Steve On Mar 16, 2015 4:24 AM, Ramon de Leon r.a.deleo...@gmail.com wrote: Hi Steve,

Re: [postgis-users] PostgreSQL 9.3.4/PostGIS 2.1.1 function

2014-12-05 Thread Stephen Mather
Hi Freddie, select st_astext(st_setsrid(st_scale('POINT(1 1)'::geometry,180.0/pi(),180.0/pi()),4326)) FROM table_containing_polygon; Best, Steve On Fri, Dec 5, 2014 at 3:33 PM, Burgess, Freddie fburg...@radiantblue.com wrote: How would I do this on a polygon? Thanks

Re: [postgis-users] PostgreSQL 9.3.4/PostGIS 2.1.1 function

2014-12-05 Thread Stephen Mather
Sorry typos and such, try this again: select st_astext(st_setsrid(st_scale(geom,180.0/pi(),180.0/pi()),4326)) FROM table_containing_polygons; Best, Steve On Fri, Dec 5, 2014 at 4:30 PM, Stephen Mather step...@smathermather.com wrote: Hi Freddie, select st_astext(st_setsrid(st_scale('POINT(1

Re: [postgis-users] ST_Within equivalent for TIN 3D

2014-08-19 Thread Stephen Mather
Hi Tsolmon, If memory serves, there's no support for TINs with ST_Within yet. I suspect (without testing the suspicion) that you can cast or convert to 3D polygons and do the test, but you'd need to use a 3D equivalent, probably ST_3DDFullyWithin, but the Oslandia folks might be able to answer

Re: [postgis-users] A question on PostGIS Cookbook

2014-07-16 Thread Stephen Mather
, 2014 at 4:33 AM, Cyril Poitevin cyril.poite...@univ-lr.fr wrote: Oops Hi Steve Sorry. Message original Objet: Re: A question on PostGIS Cookbook Date: 15.07.2014 10:08 De: Cyril Poitevin cyril.poite...@univ-lr.fr À: Stephen Mather step...@smathermather.com Hi Stephen

Re: [postgis-users] A question on PostGIS Cookbook

2014-07-16 Thread Stephen Mather
Hmm, Cyril, (going back to Paolo's earlier question) did you enable plpython3u or plpython2u on your database? Best, Steve ___ postgis-users mailing list postgis-users@lists.osgeo.org http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

Re: [postgis-users] 'Z' Difference between DEM file PG DEM layer

2014-05-26 Thread Stephen Mather
Remember the phrase lies, d*mn lies, and statistics? It applies here... . Most statistics for images are sampled, meaning the value of every pixel will not be read. In other words, you can't ordinarily trust the stats to tell you whether the data have changed on import. I believe gdalinfo will

Re: [postgis-users] PostGIS 2.1.2 and 2.0.5 released

2014-03-31 Thread Stephen Mather
I thought you were joking at first, and then read your affiliation... . On Mon, Mar 31, 2014 at 3:45 PM, Gerry Creager - NOAA Affiliate gerry.crea...@noaa.gov wrote: I'm not sure I'd call it God's country (on either your part or mine). I chose to live in Oklahoma because... that's where the

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-16 Thread Stephen Mather
to PostGIS makes no difference to my use of PostGIS. Best, Steve On Mon, Dec 16, 2013 at 12:52 PM, Mike Toews mwto...@gmail.com wrote: On 15 Dec 2013 21:50, Stephen Mather step...@smathermather.com wrote: Hmm, I know so little about Euler, strike, dip and rake. What are the advantages

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-15 Thread Stephen Mather
be interested in using Quaternion math for your rotations. You can even do rotations on an arbitrary axis. Take a look at this: http://3dgep.com/?p=1815 -Bob On Dec 14, 2013 7:57 PM, Stephen Mather step...@smathermather.com wrote: Ahem-- bugs fixed: -- Function: st_rotatex(geometry, double

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-15 Thread Stephen Mather
as part of make check, you add MYTEST to regress/Makefile.in -bborie On Sat, Dec 14, 2013 at 10:29 PM, Stephen Mather step...@smathermather.com wrote: I have ticketed. It's a little messy of a ticket (actually two), so apologies for that. http://trac.osgeo.org/postgis/ticket/2570 http

Re: [postgis-users] Old question resurfacing

2013-12-15 Thread Stephen Mather
That's nice-- it's the first one I've seen to focus on visualization of point clouds specifically. Best, Steve On Sun, Dec 15, 2013 at 3:10 PM, Frank Henze he...@tu-cottbus.de wrote: I dont know, at which point the downsampling should take place. An interesting solution for webbased

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-15 Thread Stephen Mather
Hmm, I know so little about Euler, strike, dip and rake. What are the advantages? What I do know is that it's easy now to construct a 3 axis rotation function (which also might be better handled with ST_Affine): http://trac.osgeo.org/postgis/attachment/ticket/2575/ -- Function:

[postgis-users] ST_RotateX with pointOrigin!

2013-12-14 Thread Stephen Mather
Hi All, I think I avoided doing linear algebra, which is good since I never studied it... . This is my cludgy patch for making a version of geometry ST_RotateX(geometry geomA, float rotRadians, geometry pointOrigin) It's not pretty enough to be a real patch ('cause my brain couldn't do that

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-14 Thread Stephen Mather
home rotTrans AS ( SELECT ST_Translate(the_geom, ST_X(pointOrigin), ST_Y(pointOrigin), ST_Z(pointOrigin)) AS the_geom FROM rotated ) - profit SELECT the_geom from rotTrans ; $BODY$ LANGUAGE sql VOLATILE COST 100; On Sat, Dec 14, 2013 at 8:37 PM, Stephen Mather step

Re: [postgis-users] ST_RotateX with pointOrigin!

2013-12-14 Thread Stephen Mather
...@gmail.com wrote: Sweet. Can you create a ticket and attach the patch? Some unit tests would be nice so that expected behavior is maintained... and bugs are caught. Thanks! -bborie On Sat, Dec 14, 2013 at 7:57 PM, Stephen Mather step...@smathermather.com wrote: Ahem-- bugs fixed

Re: [postgis-users] Intersections on polyhedral surfaces....

2013-12-12 Thread Stephen Mather
/ST_3DIntersects.html Mateusz Loskot, http://mateusz. loskot.net (Sent from mobile, apology for top-posting or broken quotes) On 12 Dec 2013 03:31, Stephen Mather step...@smathermather.com wrote: Hi All, I assume that ST_Intersection is not yet able to do intersections on polyhedralsurfaces

Re: [postgis-users] Intersections on polyhedral surfaces....

2013-12-12 Thread Stephen Mather
Brilliant. I'll give it a whirl then. On Thu, Dec 12, 2013 at 3:20 AM, Oliver Courtin olivier.cour...@oslandia.com wrote: Le 12 déc. 2013 à 08:42, Rémi Cura a écrit : waht would you mean by intersection in 3D space? volumetric intersection, or 2D intersection of colinear geometries? Both

Re: [postgis-users] Intersections on polyhedral surfaces....

2013-12-12 Thread Stephen Mather
The docs don't have ST_3DIntersection, just ST_3DIntersects: http://postgis.net/docs/manual-2.1/reference.html#Geometry_Constructors http://postgis.net/docs/manual-dev/reference.html#Geometry_Constructors And, my 2.1.2 Ubuntu build (with SFCGAL) doesn't have the function. Would I have it if I

Re: [postgis-users] ST_Tesselate

2013-12-09 Thread Stephen Mather
Ah, very cool. So hard breaklines plus a point cloud could be fed in creating a DTM... . On Mon, Dec 9, 2013 at 3:26 AM, Oliver Courtin olivier.cour...@oslandia.com wrote: Le 8 déc. 2013 à 00:23, Stephen Mather a écrit : What does ST_Tesselate do? Triangulation, a constrained Delaunay

Re: [postgis-users] ST_Tesselate

2013-12-09 Thread Stephen Mather
...@shaw.ca wrote: Why would such a function be named Tesselate? On 09/12/2013 12:26 AM, Oliver Courtin wrote: Le 8 déc. 2013 à 00:23, Stephen Mather a écrit : What does ST_Tesselate do? Triangulation, a constrained Delaunay one. http://doc.cgal.org/latest/Triangulation_2/index.html

Re: [postgis-users] postgis.backend -- how to set?

2013-12-09 Thread Stephen Mather
Thanks Remi, I'll proceed cautiously. Best, Steve On Mon, Dec 9, 2013 at 9:59 AM, Rémi Cura remi.c...@gmail.com wrote: I had a bug with this so be a little carefull ;-) Cheers, Rémi-C 2013/12/8 Stephen Mather step...@smathermather.com Awesome. Thank you. On Sun, Dec 8, 2013

Re: [postgis-users] ST_MinkowskiSum

2013-12-09 Thread Stephen Mather
:24 AM, Oliver Courtin olivier.cour...@oslandia.com wrote: Le 8 déc. 2013 à 05:41, Stephen Mather a écrit : Hi Stephen, but I can't say I quite understand the product it is creating... . Did you have a look on: http://en.wikipedia.org/wiki/Minkowski_addition http://doc.cgal.org/latest

[postgis-users] ST_MinkowskiSum

2013-12-07 Thread Stephen Mather
Hi All, Digging SFCGAL... . Quick question-- is ST_MinkowskiSum a 2D calculation of the Minkowski Sum when non-polyhedral surfaces are fed in? For example, I tried this: WITH pointt AS ( SELECT ST_MakePoint(0,0,0) AS the_geom ), buffered AS ( SELECT ST_Buffer(the_geom, 30, 20) AS the_geom

Re: [postgis-users] ST_MinkowskiSum

2013-12-07 Thread Stephen Mather
On Sat, Dec 7, 2013 at 5:51 PM, Stephen Mather step...@smathermather.comwrote: Hi All, Digging SFCGAL... . Quick question-- is ST_MinkowskiSum a 2D calculation of the Minkowski Sum when non-polyhedral surfaces are fed in? For example, I tried this: WITH pointt AS ( SELECT

[postgis-users] ST_Tesselate

2013-12-07 Thread Stephen Mather
What does ST_Tesselate do? Thanks, Best, Steve ___ postgis-users mailing list postgis-users@lists.osgeo.org http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

Re: [postgis-users] ST_MinkowskiSum

2013-12-07 Thread Stephen Mather
Hmm, partially answering my own question, but it's obvious the output is not 2D, but I can't say I quite understand the product it is creating... . Best, Steve On Sat, Dec 7, 2013 at 5:55 PM, Stephen Mather step...@smathermather.comwrote: Hit send too soon... Hi All, Digging SFCGAL

[postgis-users] -

2013-12-03 Thread Stephen Mather
If I create an nd aware index, say: CREATE INDEX lidar_the_geom_3dx ON lidar USING gist(the_geom gist_geometry_ops_nd); Can - and # be made to use that? If so, how? Thanks, Best, Steve ___ postgis-users mailing list postgis-users@lists.osgeo.org

Re: [postgis-users] -

2013-12-03 Thread Stephen Mather
3, 2013 at 12:05 PM, Stephen Mather step...@smathermather.com wrote: If I create an nd aware index, say: CREATE INDEX lidar_the_geom_3dx ON lidar USING gist(the_geom gist_geometry_ops_nd); Can - and # be made to use that? If so, how? Thanks, Best, Steve

[postgis-users] ST_DWithin with nd index

2013-12-03 Thread Stephen Mather
Hi All, Will ST_DWithin work against an nd index? I'm getting back the same number of records, whether I use a standard 2D gist or nd, only running it with the nd index takes 5 times as long. Thanks, Best, Steve ___ postgis-users mailing list

Re: [postgis-users] ST_DWithin with nd index

2013-12-03 Thread Stephen Mather
Ahem, dumb question:tired brain. ST_3DDWithin http://postgis.net/docs/ST_3DDWithin.html On Tue, Dec 3, 2013 at 5:06 PM, Stephen Mather step...@smathermather.comwrote: Hi All, Will ST_DWithin work against an nd index? I'm getting back the same number of records, whether I use a standard 2D

[postgis-users] Readable reference for - the distance operators

2013-11-20 Thread Stephen Mather
Hi All, What's the best (normal) human readable reference for the PostGIS distance operators-- but deep enough to talk about the use of index structure in the use of the operators? Thanks, Best, Steve ___ postgis-users mailing list

Re: [postgis-users] Readable reference for - the distance operators

2013-11-20 Thread Stephen Mather
On Wed, Nov 20, 2013 at 10:49 AM, Stephen Mather step...@smathermather.com wrote: Hi All, What's the best (normal) human readable reference for the PostGIS distance operators-- but deep enough to talk about the use of index structure in the use of the operators? Thanks, Best

Re: [postgis-users] The first release of the PostGIS Add-ons is out!

2013-11-19 Thread Stephen Mather
How could we avoid a license discussion... :) . I think there are two relevant points from the discussion: 1) It should be licensed such that it can go into core if such an option for some or all functions reveals itself in the future 2) It should be flexible enough for wide reuse. #1

Re: [postgis-users] The first release of the PostGIS Add-ons is out!

2013-11-19 Thread Stephen Mather
Oh, I prefer GPL, and in this case, since it's uncompiled code it probably doesn't matter. I just know lots of people belly ache about the GPL, so I'm trying to insert that perspective... . :) On Tue, Nov 19, 2013 at 2:42 PM, Pierre Racine pierre.rac...@sbf.ulaval.cawrote: IMHO, It's fine

Re: [postgis-users] Σχετ: postgis

2013-05-17 Thread Stephen Mather
That said, there is value to implementing dependency-free routing algorithms in the core of PostGIS, at least simpler stuff anyway. I have dijkstra on my todo for PostGIS, but haven't gotten started on it yet. Best, Steve On May 17, 2013, at 11:00 AM, Stephen Woodbridge

Re: [postgis-users] Need a method for noding a street network

2013-05-09 Thread Stephen Mather
A couple quick thoughts on approach 2: 1) The union need only apply to intersecting geometries, rather than the whole dataset. This helps considerably with memory footprint, but in some edge cases could still be a real problem. 2) We need Mr. Davis to hurry up and demonstrate streaming geometry

Re: [postgis-users] Count Points in Polygon with Postgis

2013-03-16 Thread Stephen Mather
Try unioning with the opposite condition. Sent from my iPod On Mar 16, 2013, at 11:53 AM, Carsten Hogertz carsten.hoge...@gmail.com wrote: I've got a simple problem: I want to count the number of points within a set of polygons. I have a SQL already but it only gives back the gid of the

[postgis-users] ST_DelaunayTriangles

2013-01-21 Thread Stephen Mather
Hi All, Saw this post: http://www.bostongis.com/blog/index.php?/archives/195-Waiting-for-PostGIS-2.1-Delaunay-Triangulation-with-3D.html on Regina and Leo's blog on the addition of Delauney triangulation to PostGIS. Quick question about this-- for 3D this is predicated on GEOS 3.4. Is