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,
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
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,
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
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
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
, 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
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
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
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
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
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
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
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
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:
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
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
...@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
/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
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
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
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
...@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
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
: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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
41 matches
Mail list logo