[postgis-users] The pcpatch datatype

2015-02-13 Thread Jonathan Moules
Hi List, I'm new to PostGIS and have inherited a database which is being used to store a pointcloud. The data seems to be in a "pcpatch" Data type, and this is being converted into points for extraction with: "select pc_explode(pa)::geometry ...". The problem is, I can find no reference to pcpa

Re: [postgis-users] The pcpatch datatype

2015-02-13 Thread Jonathan Moules
From: postgis-users-boun...@lists.osgeo.org<mailto:postgis-users-boun...@lists.osgeo.org> mailto:postgis-users-boun...@lists.osgeo.org>> on behalf of Jonathan Moules mailto:j.mou...@hrwallingford.com>> Sent: Friday, February 13, 2015 7:46 AM To: 'postgis-users@list

Re: [postgis-users] 2.1.6 Released

2015-03-26 Thread Jonathan Moules
Hi Paul, I suspect I know the answer to this (“no”), but to confirm, is the disk saving for points going to apply to points stored using the pointcloud extension in a postgis database? Thanks, Jonathan From: postgis-users-boun...@lists.osgeo.org [mailto:postgis-users-boun...@lists.osgeo.org] O

Re: [postgis-users] 2.1.6 Released

2015-03-26 Thread Jonathan Moules
patch and sensible xml schema, you should already have a 1:2 to 1:4 compression of the point size. (compared to binary ply in my case) Cheers, Rémi-C 2015-03-26 18:05 GMT+01:00 Jonathan Moules mailto:j.mou...@hrwallingford.com>>: Hi Paul, I suspect I know the answer to this (“no”), but to con

Re: [postgis-users] 2.1.6 Released

2015-03-27 Thread Jonathan Moules
pc_explode(pa) as pt ) SELECT block_id , pc_get(pt, 'X') as x, pc_get(pt, 'Y') as y, pc_get(pt, 'Z') as z FROM points Cheers, Rémi-C 2015-03-26 19:09 GMT+01:00 Jonathan Moules mailto:j.mou...@hrwallingford.com>>: Hi Rémi, Thanks for the reply and information.

[postgis-users] Distance between two furthest points of a group

2015-06-26 Thread Jonathan Moules
Hi List, I have sets of points (up to 250,000 in a set) and I want to get the furthest distance between any of them. In theory the simplest way is to use ST_MinimumBoundingCircle(ST_Collect(geography) and then get the diameter of that. The problem is - I don't seem to be able to get the diamet

Re: [postgis-users] Distance between two furthest points of a group

2015-07-13 Thread Jonathan Moules
tance(a.geom,b.geom) d from f a cross join f b order by d desc limit 1; ​should give you the id of your targets and the distance between them​ On Fri, Jun 26, 2015 at 2:54 PM, Jonathan Moules mailto:j.mou...@hrwallingford.com>> wrote: Hi List, I have sets of points (up to 250,0

[postgis-users] PointCloud to Delaunay Triangles

2015-08-05 Thread Jonathan Moules
Hi List, I have a pointcloud, and I want to turn it into a TIN (so I guess ST_DelaunayTriangles), and then measure the length of the sides of the triangles. Does anyone have any suggestions on how to do this, ideally in a performant fashion? I've tried a few things, but I'm currently struggli

[postgis-users] TOAST and PostGIS

2016-03-22 Thread Jonathan Moules
Hi List, Looking at - http://postgis.net/docs/performance_tips.html#small_tables_large_objects - item 7.1 describes a problem: "Current PostgreSQL versions (including 8.0) suffer from a query optimizer weakness regarding TOAST tables." I was wondering if the issue is resolved in the 9.x serie

[postgis-users] Spatially ordering data

2017-03-01 Thread Jonathan Moules
Hi List, I've got a database of about 60 million spatial features that I need to run through a process. The process can't do all of them at once, so I want to get the data from postgis in spatial chunks of up to say 500,000 features per query. The data itself isn't evenly spatially distributed,

Re: [postgis-users] Spatially ordering data

2017-03-02 Thread Jonathan Moules
ash-sorting-with-these-three-simple-tricks/ for simplest case, you can just order by your geom field. ср, 1 мар. 2017 г. в 14:44, Jonathan Moules <jonathan-li...@lightpear.com>: Hi List, I've got a database of about 60 million spatial features that I need to run through a process.

Re: [postgis-users] Spatially ordering data

2017-03-06 Thread Jonathan Moules
Hi List, Thanks for all the suggestions. I think I've solved this using a combination of methods from the responses and wanted to share for anyone in the future. Things I tried that didn't work for me (beyond what was in the previous post): * Convert to a point, then extract the Lat/Lon

[postgis-users] ST_Intersects

2017-03-17 Thread Jonathan Moules
Hi List, I'm doing a simple ST_Intersects: select geom from TABLENAME where ST_Intersects(geom, ST_GeomFromText('POLYGON((26 655000, 26 66, 27 66, 27 655000, 26 655000))', 27700)) TABLENAME does have a spatial index. But this query is still taking an

Re: [postgis-users] ST_Intersects

2017-03-20 Thread Jonathan Moules
athan ---- On Fri, 17 Mar 2017 14:19:55 + Giuseppe Broccolo<giuseppe.brocc...@2ndquadrant.it> wrote Hi Jonathan, 2017-03-17 13:18 GMT+01:00 Jonathan Moules <jonathan-li...@lightpear.com>: Hi List, I'm doing a simple ST_Intersects: select geom from TABLENAME