[pdal] Clipping from multiple LAS tiles

Adam Steer Adam.Steer at anu.edu.au
Mon Oct 24 21:56:41 PDT 2016


Hi PDAL Ninjas

I’m wondering how we can replicate LAStool’s multi-tile clipping fiunction, eg: lasclip -i folder/*.las -poly polygon.shp -o newfile.las

I’ve got a pretty clunky pipeline that works, based on the github pipeline examples:

{
    "pipeline": [
        { "filename": "Merimbula2013-C3-AHD_7585912_55_0002_0002.las",
          "tag": "A"
        },
	{ "filename": "Merimbula2013-C3-AHD_7585914_55_0002_0002.las",
          "tag": "B"
        },
	{ "filename": "Merimbula2013-C3-AHD_7605914_55_0002_0002.las",
          "tag": "C"
        },
	{
            "inputs": ["A", "B", "C"],
            "type": "filters.crop",
            "polygon": "POLYGON((759094.480855234 5913008.2715937095,758464.6999094139 5912716.199270982,757743.6463627518 5912898.7444$
            "outside": false
        },
	"./merimbulatown.las"
    ]
}


…but would prefer to not have to list all my tiles. The obvious:

    "pipeline": [
         "../merimbula2013/Tiles_2k_2k/\*.las",
        {
            "type": "filters.crop”,

…didn’t work.

I’ve read the examples from the workshop/tutorial samples - and can’t put a finger on how to extend those either.

I see that tindex might be a way, but I’d like to do the job without constructing another file. Hmm. I feel like I’m missing something really obvious.


Any suggestions will be appreciated!

Thanks

Adam


More information about the pdal mailing list