<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Hi,</p>
<p>I am sending this small script as an inspiration. This is a
single-threaded the python program that<br>
</p>
<ol>
<li>Scans a directory for a list of .tif file and makes a list of
those tif files.<br>
</li>
<li>Creates an output directory</li>
<li>Build a gdal command line (export the tif to a xyz file)
inside the function called "gdal_translate2XYZ". <br>
</li>
<ol>
<li>shell = subprocess.Popen(r'c:\OSGeo4W64\OSGeo4W.bat',
stdin=subprocess.PIPE)</li>
<li>shell.communicate(gdal_translateCmd)</li>
<li>These two above line basically send the command to the
OSGeo4W shell. (Install Qgis if you don't have it).<br>
</li>
</ol>
<li>Runs this command in a loop (the loop is created by using
"map" or "p.map" function that operates on the list of tif
files)</li>
<li>If you comment out the (if __name__ == '__main__': ) section,
the commands will run in parallel depending on the size of the
pool variable. pool = 6 = 6 parallel commands. It is easier to
debug if you develop in single thread. <br>
</li>
</ol>
<p>The idea is to build your script using python and to run it
calling PDAL from python. The Fiona python library with the
Shapely library or the more complete geopandas library can help
manipulate the shapefile.</p>
<p>Hope this helps!</p>
<p>Nicolas <br>
</p>
<p><a href="https://fiona.readthedocs.io/en/stable/">https://fiona.readthedocs.io/en/stable/</a></p>
<p><a href="https://shapely.readthedocs.io/en/stable/manual.html">https://shapely.readthedocs.io/en/stable/manual.html</a></p>
<p><a href="http://geopandas.org/">http://geopandas.org/</a></p>
<p><br>
</p>
<div class="moz-cite-prefix">On 2019-12-08 8:09 p.m., Jason McVay
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAGHb011VeJ21_uyRU5aXrpZ0fwMs0sS13X6BD6eJtYqkb3NPJw@mail.gmail.com">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<div dir="ltr">I'm looking for some advice on the best way/how to
loop in thousands of bounding coordinates into a pdal pipeline.
<div><br>
</div>
<div>I have a csv (and a geojson) of several thousand min/max
x/y and a unique ID. The AOI's are not very big, so the
pipeline runs quickly, but there are a lot of AOIs to capture!
I'm querying an entwine dataset, the extent of which is
national, so I'm limiting the data with a bounding box of each
AOI.</div>
<div><br>
</div>
<div>My pipeline currently runs HAG and Ferry Z filter, then
uses the gdal.writer to make a GeoTiff at 1m resolution. It
works perfectly when I manually enter in a set of test
coordinates. How can I scale this to loop and update the
bounds automatically?</div>
<div><br>
</div>
<div>I'm running this locally on a MacBook Pro.</div>
<div><br>
</div>
<div>Thank you, any advice is appreciated!<br>
<div><br>
</div>
<div><br clear="all">
<div>
<div dir="ltr" class="gmail_signature"
data-smartmail="gmail_signature">
<div dir="ltr">
<div dir="ltr"><span>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Jason McVay<br>
</div>
</div>
<br>
</div>
MS Geography, Virginia Tech<br>
</div>
BA Environmental Studies, University of
Montana<br>
<a
href="http://www.linkedin.com/in/jasonmcvay86/"
title="View public profile"
name="SignatureSanitizer_SafeHtmlFilter_UNIQUE_ID_SafeHtmlFilter_14933b4cad52e0be_SafeHtmlFilter_webProfileURL"
target="_blank" moz-do-not-send="true">www.linkedin.com/in/jasonmcvay86/</a></div>
</div>
<div><br>
</div>
</span></div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
pdal mailing list
<a class="moz-txt-link-abbreviated" href="mailto:pdal@lists.osgeo.org">pdal@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/pdal">https://lists.osgeo.org/mailman/listinfo/pdal</a></pre>
</blockquote>
</body>
</html>