[pdal] assign projection info without creating new file

[Andrew Bell]

On Mon, Nov 5, 2018 at 5:32 PM Matt Beckley <beckley at unavco.org> wrote:

> Hello All,
>
> I have thousands of .laz files that have no projection information.  From
> my contact, I received the EPSG codes for the datasets.  I can apply them
> with a pipeline such as:
>
> {
>
>     "pipeline": [
>
>         {
>
>             "type" : "readers.las",
>
>                "filename": "input.laz"
>
>            },
>
>         {
>
>             "type" : "writers.las",
>
>                "filename": "output.laz",
>
>             "a_srs": "EPSG:2193+7839"
>
>            }
>
>     ]
>
> }
>
>
> I suppose I could just then do  "mv output.laz input.laz" in a bash script
> to end up with a single file with the proper projection in the header.
>
> My question: Is there a way I can do it without creating an intermediary
> file?  I tried the obvious of reading and writing to input.laz, but that
> didn't work.  With thousands of files, it seems like it could save some
> time.  Just curious.
>

Not with PDAL.  You may be able to do this with lastools, but you also may
not be able to do it at all, depending on how the file is formatted.  The
SRS information is stored in VLRs, and if the VLRs aren't large enough to
fit the new SRS, you're out of luck.

-- 
Andrew Bell
andrew.bell.ia at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20181105/155f26d6/attachment.html>