<div dir="ltr"><div>Well, after reading and asking and your response to Howard and Jim and that of Even from PROJ, it seems that I have achieved it, thank you.<br>It would still be necessary to take into account Jim Klassen's considerations about what he says in point 6), but at least the change in heights is now concrete:</div><div><br></div><div>[<br>    "./input_WGS84_UTM14N.las",<br>    {<br>        "type":"filters.projpipeline",<br>        "coord_op":"+proj=pipeline +step +inv +proj=utm +zone=11 +ellps=WGS84 +step +proj=vgridshift +grids=/nubes/geoides/GGM10_GEO.gtx +step +proj=utm +zone=11 +ellps=WGS84"<br>    }, <br>    {<br>        "type":"writers.las",<br>        "scale_x":"0.001",<br>        "scale_y":"0.001",<br>        "scale_z":"0.001",<br>        "offset_x":"570000",<br>        "offset_y":"3590000",<br>        "offset_z":"0",<br>        "global_encoding":1,<br>        "a_srs":"EPSG:32611+5703",<br>        "filename":"./output_WGS84_UTM14N_GGM2010.las"<br>    }<br>]</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 14 mar 2024 a las 21:22, Ulises Ibarra (<<a href="mailto:ulisesmartinibarra@gmail.com">ulisesmartinibarra@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Thank you very much for the answers:<br><br>In previous versions of pdal, the result of the vertical datum change was an output.las file with a new height system, still with the legend "Warning 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order."<br><br>But currently I have not been able to achieve it, I am using pdal version 2.6.3. And the output.las file does not reflect the vertical datum change.<br><br>I did several pipeline tests, for example:<br><br>[<br>     "./input_WGS84_UTM14N.las",<br>     {<br>         "type":"filters.reprojection",<br>         "in_srs":"EPSG:32614",<br>         "out_srs":"+proj=utm +zone=14 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +geoidgrids=C:/clouds/geoids/GGM10_GEO.gtx"<br>     },<br>     {<br>         "type":"writers.las",<br>          "a_srs":"EPSG:32614+5703",<br>         "filename":"./output_WGS84_UTM14N_GGM2010.las"<br>     }<br>]<br><br>and there is no vertical datum change in the output.<br><br>Later I tried to use the filters.projpipeline filter, but I don't know how to tell it "change vertical datum using the geogrid.gtx achvio"<br><br>I used this pipeline:<br>[<br>     "./input_WGS84_UTM14N.las",<br>     {<br>         "type":"filters.projpipeline",<br>         "coord_op":"+proj=pipeline +step +proj=vgridshift +grids=C:/clouds/geoids/GGM10_GEO.gtx",<br>         "out_srs":"EPSG:32614+5703"<br>     },<br>     {<br>         "type":"writers.las",<br>         "a_srs":"EPSG:32614+5703",<br>         "filename":"./output_WGS84_UTM14N_GGM2010.las"<br>     }<br>]<br><br>and I got this:<br><br>(pdal pipeline filters.projpipeline Error) GDAL failure (1) PROJ: vgridshift: Invalid latitude<br><br><br>Could you guys give me a clue what I'm doing wrong?<br><br>Cordially<br><br>Ulises Ibarra</div>
</blockquote></div>