A GDAL tile index is an OGR-readable data source of boundary information. PDAL provides a similar concept for PDAL-readable point cloud data. You can use the tindex command application to generate tile index files in any format that OGR supports writing. Once you have the tile index, you can then use the readers.tindex driver to automatically merge and query the data described by the tiles.
Given a tile index that was generated with the following scenario:
pdal tindex index.sqlite \
"/Users/hobu/dev/git/pdal/test/data/las/interesting.las" \
-f "SQLite" \
--lyr_name "pdal" \
--t_srs "EPSG:4326"
Use the following Pipeline example to read and automatically merge the data.
{
"pipeline":[
{
"type":"readers.tindex",
"sql":"SELECT * from pdal",
"filter_srs":"+proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=399999.9999999999 +y_0=0 +ellps=GRS80 +units=ft +no_defs",
"merge":"true",
"filename":"index.sqlite",
"where":"location LIKE \'%nteresting.las%\'",
"boundary":"([635629.85, 638982.55], [848999.70 , 853535.43])",
"polygon":"POLYGON ((635629.85000000 848999.70000000, 635629.85000000 853535.43000000, 638982.55000000 853535.43000000, 638982.55000000 848999.70000000, 635629.85000000 848999.70000000))"
},
{
"type":"writers.las",
"filename":"outputfile.las"
}
]
}
wkt
option.wkt
or boundary
option to this
coordinate system before filtering or reading data.lyr_name