<p dir="ltr">Paul<br>
Is creating a spatial index like this a separate step?  I missed that in your tutorial.<br>
Thanks<br>
Bruce</p>
<div class="gmail_quote">On Aug 31, 2013 8:00 AM, "Paul Ramsey" <<a href="mailto:pramsey@cleverelephant.ca">pramsey@cleverelephant.ca</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Roger,<br>
<br>
I'm going to guess you don't have an "appropriate" spatial index on the table, so it has to scan the whole thing twice: once to figure out how big the buffer is (PDAL inefficiency) and once to pull the data.<br>

<br>
CREATE INDEX pa_gix ON all_returns USING GIST ( geometry(pa) );<br>
<br>
P.<br>
<br>
--<br>
Paul Ramsey<br>
<a href="http://cleverelephant.ca/" target="_blank">http://cleverelephant.ca/</a><br>
<a href="http://postgis.net/" target="_blank">http://postgis.net/</a><br>
<br>
<br>
On Saturday, 31 August, 2013 at 6:55 AM, Roger Bedell wrote:<br>
<br>
> Hello again,<br>
><br>
> Ok, finally got a reasonable number of points to import into our database. Took a while because of a hardware failure, but now it is done.<br>
><br>
> I imported about 2500 LAS files, each one about 8mb into a single table. The resulting table has 55 million records (counted not estimated). Each record has 400 points. Database size is about 375gb.<br>
><br>
> I am attempting to do the following from pcpipeline:<br>
><br>
> <?xml version="1.0" encoding="utf-8"?><br>
> <Pipeline version="1.0"><br>
> <Writer type="drivers.las.writer"><br>
> <Option name="filename">test_small.las</Option><br>
> <Option name="spatialreference">EPSG:4326</Option><br>
> <Reader type="drivers.pgpointcloud.reader"><br>
> <Option name="connection">host='ogisrv102' port='8080' dbname='NRCS_LIDAR' user='postgres'password='5Tb3CwBu353w' </Option><br>
> <Option name="table">all_returns</Option><br>
> <Option name="column">pa</Option><br>
> <Option name="srid">4326</Option><br>
> <Option name="where">PC_Intersects(pa, ST_MakeEnvelope(-97.27021, 34.46225, -97.26987, 34.46251, 4326))</Option><br>
> </Reader><br>
> </Writer><br>
> </Pipeline><br>
><br>
><br>
> Start time 8:24, end time 8:30 for the program to get here:<br>
><br>
> F:\PDAL64>pdal64\pcpipeline -i export_pg_to_las_test.xml<br>
> Requested to read 700 points<br>
> Requested to write 700 points<br>
> Buffer capacity is 700<br>
> 0<br>
><br>
><br>
> Finally completes at 08:35:<br>
><br>
> F:\PDAL64>pdal64\pcpipeline -i export_pg_to_las_test.xml<br>
> Requested to read 700 points<br>
> Requested to write 700 points<br>
> Buffer capacity is 700<br>
> 0.100<br>
><br>
><br>
><br>
><br>
> The new database server has reasonable horsepower. 2 Xeon e5-2690 CPUs, 32gb RAM. Postgresql is 9.2.4, 64 bit on a Windows 2008 Server OS. The 64 bit pgpointcloud build for Windows was provided by Regina. I built the 64 bit PDAL.<br>

><br>
> 11 minutes seems excessive for extracting 700 points. Am I missing something, or is this to be expected?<br>
><br>
> Thanks for any insight,Roger<br>
><br>
> --<br>
> Roger Bedell<br>
> Coordinate Solutions Inc.<br>
> _______________________________________________<br>
> pgpointcloud mailing list<br>
> <a href="mailto:pgpointcloud@lists.osgeo.org">pgpointcloud@lists.osgeo.org</a> (mailto:<a href="mailto:pgpointcloud@lists.osgeo.org">pgpointcloud@lists.osgeo.org</a>)<br>
> <a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/pgpointcloud" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/pgpointcloud</a><br>
<br>
<br>
<br>
_______________________________________________<br>
pgpointcloud mailing list<br>
<a href="mailto:pgpointcloud@lists.osgeo.org">pgpointcloud@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/pgpointcloud" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/pgpointcloud</a><br>
</blockquote></div>