<p dir="ltr">No worries.  That may be why I am just slightly slower.  I will test next week.<br>
Thanks</p>
<div class="gmail_quote">On Aug 31, 2013 8:22 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">
Bruce,<br>
I may well have failed to document it. No excuses :(<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 8:16 AM, Bruce Rindahl wrote:<br>
<br>
> Paul<br>
> Is creating a spatial index like this a separate step? I missed that in your tutorial.<br>
> Thanks<br>
> Bruce<br>
> On Aug 31, 2013 8:00 AM, "Paul Ramsey" <<a href="mailto:pramsey@cleverelephant.ca">pramsey@cleverelephant.ca</a> (mailto:<a href="mailto:pramsey@cleverelephant.ca">pramsey@cleverelephant.ca</a>)> wrote:<br>

> > 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>) (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>
> > _______________________________________________<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>
</blockquote></div>