[pdal] tindex creating empty geometry column

[Howard Butler]

Adam,

We're just using OGR to write the SQLite file, which I think defaults to using spatialite, but I'm not exactly sure.

You might try writing to a shapefile or something to see if it is merely the OGR driver falling over:

> $ pdal tindex act_8ppm_index4326.shp  "LAS/AHD/*.las" -f "ESRI Shapefile" --lyr_name "pdal" --t_srs "EPSG:4326" —fast_boundary

Hope this helps,

Howard

> On Feb 23, 2017, at 7:17 PM, Adam Steer <Adam.Steer at anu.edu.au> wrote:
> 
> I’m making tile indexes for a set of LAS tiles, and getting no geometry column:
> 
> # pdal tindex act_8ppm_index4326.sqlite "LAS/AHD/*.las" -f "SQLite" --lyr_name "pdal" --t_srs "EPSG:4326" —fast_boundary
> 
> # sqlite3 act_8ppm_index4326.sqlite
> 
> sqlite> PRAGMA table_info(pdal);
> 0|OGC_FID|INTEGER|0||1
> 1|GEOMETRY|BLOB|0||0
> 2|location|VARCHAR(254)|0||0
> 3|srs|VARCHAR(254)|0||0
> 4|modified|TIMESTAMP|0||0
> 5|created|TIMESTAMP|0||0
> 
> sqlite> select * from pdal;
> 34||/mnt/pointcloud/act_all/8ppm/LAS/AHD/ACT2015_8ppm-C3-AHD_6886100_55_0002_0002.las|EPSG:4326|2016-11-25T01:08:31|2017-02-07T00:02:54
> 35||/mnt/pointcloud/act_all/8ppm/LAS/AHD/ACT2015_8ppm-C3-AHD_6866076_55_0002_0002.las|EPSG:4326|2016-11-25T00:53:45|2017-02-07T00:02:54
> …
> 
> …here, the second column GEOMETRY is empty - it should contain a blob of geometry data.
> 
> Where have I gone astray?
> 
> On this topic - do hooks exist in the tindex system to extract these data (tile boundaries) from a postGIS store?
> 
> Thanks
> 
> Adam
> 
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal