<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 28, 2017 at 12:11 AM, Howard Butler <span dir="ltr"><<a href="mailto:howard@hobu.co" target="_blank">howard@hobu.co</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-"><br>
> On Feb 27, 2017, at 10:31 PM, Adam Steer <<a href="mailto:Adam.Steer@anu.edu.au">Adam.Steer@anu.edu.au</a>> wrote:<br>
><br>
> Thanks Andrew, Howard.<br>
><br>
> I can create geometries which can be loaded into QGIS - so clearly the lack of some binary blob symbol in sqlite’s 'select * from…’ is not a real issue.<br>
><br>
> It appears the OGR drivers are working, and I’ve learned a lot in the past two days - nearly answering my own question about PostGIS. Short summary: OGR has a method for injecting data straight into a postgres store - but PDAL seems to want a filename to write to. I’ll revisit that later.<br>
<br>
</span>I think you mean the ExecuteSQL method of OGR datasource, and no we don't expose that to the readers.tindex at the moment.<br>
<span class="gmail-"><br>
<br>
><br>
> So why is merging inside a clipping polygon from a tindex failing then? Following this example: <a href="https://www.pdal.io/stages/readers.tindex.html" rel="noreferrer" target="_blank">https://www.pdal.io/stages/<wbr>readers.tindex.html</a><br>
><br>
> I get:<br>
> PDAL: readers.tindex: Unexpected argument 'merge’.<br>
<br>
</span>hmm, something must be broken here. I would expect this to work..</blockquote><div><br></div><div>The error is correct.  The example in the documentation is wrong.  You should probably use the kernel to extract data based on the index.  It builds a pipeline and invokes the stage "readers.tindex" automatically.  I'll try to update the documentation.</div><div><br></div><div>From your previous tindex creation I think it should be something like:</div><div><br></div><div>pdal tindex --merge --tindex "/data/act_8ppm.sqlite" --lyr_name pdal --filename "/data/outputfile.laz" --polygon "POLYGON((149.0896224975586 -35.25812600834376,149.08790588378906 -35.26485397359959,149.08721923828125 -35.27326314463213,149.08653259277344 -35.27942931528233,149.0899658203125 -35.284193761969895,149.09820556640625 -35.28363325338022,149.10472869873047 -35.28111091671122,149.10438537597656 -35.27046018460613,149.10850524902344 -35.26933897343947,149.1126251220703 -35.26289170810056,149.10747528076172 -35.25980805202397,149.09889221191406 -35.260368725492896,149.0896224975586 -35.25812600834376))"</div><div><br></div><div>If you want to filter the layer, you'll have to add --readers.tindex.where "location LIKE \"%ACT2015_8ppm-C3-AHD%\"</div><div><br></div></div>-- <br><div class="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</div></div>