<div dir="ltr"><div>Hola Inti</div><div><br></div><div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 15 Sept 2022 at 22:18, inti luna <<a href="mailto:inti.luna.rsgis@gmail.com">inti.luna.rsgis@gmail.com</a>> wrote:<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"><p>Dear PROJ community,<br>
I am using a GNSS (Emlid REACH RS2 )to collect ground control points for drone mapping
and supporting a project of monitoring vegetation after fire. I have collected data in LAT,LON (CRS EPSG 4326) and ellipsoidal
height (m) and would like to transform it to another CRS. Target
horizontal is EPSG: 25831 and target height is based on EGM2008 (EPSG
3855).</p></div></blockquote><div>That RTK device should not produce vales in WGS84. The coordinates you provide are in Spain. I would assume that the reference system used by the base stations is ETRS89 (2D: EPSG:4258 , 3D EPSG:4937).</div><div>In addition to that, Spain has its own vertical CRS, Alicante height, EPSG:5782</div><div>So probably your source CRS would be EPSG:4937 and destination EPSG:25830+5782<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"><p>I just discovered PROJ and I am having issues with the vertical transformation. <br></p><p>Using Linux Ubuntu 11 and installed PROJ using anaconda with:</p></div></blockquote><div>Ubuntu 11? <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"><p><conda install -c conda-forge proj></p><p><conda install -c conda-forge proj-data></p><p>I tested with:</p><p><echo -2.69713764 43.37551925 | cs2cs +init=epsg:4326 +to +init=epsg:25830></p><p>which outputs: -2.69713764 43.37551925 which is fine.</p><p>However, I don't understand very well how to perform the vertical grid shift.</p><p>I tried:</p><p>1. <echo -2.69713764 43.37551925 102.635| cs2cs +init=epsg:4326 +vgridshift +to +init=epsg:25830 +grids=egm08_25.gtx></p><p>2.<echo -2.69713764 43.37551925 102.635 | proj +init=epsg:4326 +proj=vgridshift +grids=egm08_25.gtx><br></p></div></blockquote><div>This is working for me:</div><div>Using from WGS84 3D to EPSG:25830 + EGM2008<br></div><div><span style="font-family:monospace">echo 43.37551925 -2.69713764 102.635 | PROJ_NETWORK=ON cs2cs EPSG:4979 EPSG:25830+3855</span></div><div><span style="font-family:monospace">524534.95 4802561.48 53.07</span></div><div><br></div><div>Using from ETRS89 3D to EPSG:25830 + Alicante height<br></div><div><span style="font-family:monospace">echo 43.37551925 -2.69713764 102.635 | PROJ_NETWORK=ON cs2cs EPSG:4937 EPSG:25830+5782</span></div><div><span style="font-family:monospace">524534.95 4802561.48 53.72</span></div><div><br></div><div>PROJ_NETWORK=ON is not needed if proj-data is properly installed.<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"><p></p><p>But the first try ignores the height and the second got me an error message:</p><p>"can't initialize operations that produce angular output coordinates<br>program abnormally terminated"</p><p>Thank you for your help,</p><p>Inti<br></p><p><br></p><p><br></p><br>-- <br><div dir="ltr"><div dir="ltr"><div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div><font style="color:rgb(153,153,153);font-size:x-small" face="arial black, sans-serif">Inti Luna<br></font></div><div><font style="color:rgb(153,153,153);font-size:x-small" face="arial black, sans-serif">+
</font><font style="color:rgb(153,153,153);font-size:x-small" face="arial black, sans-serif">34 688756449<br></font></div><br><div><font size="1" face="arial black, sans-serif" color="#999999"><br></font></div>
</div>
</div>
</div>
</div></div></div></div>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br></blockquote><div><br></div><div><div dir="ltr"><br clear="all"><div><div dir="ltr" class="gmail_signature">.___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br><br></div></div></div> </div></div></div></div>